Subversion Repositories svnkaklik

Rev

Rev 969 | Only display areas with differences | Ignore whitespace | Details | Blame | Last modification | View Log

Rev 969 Rev 971
1 
set terminal png
 1 
set terminal png
2 
 
 2 
 
3 
set key off
 3 
set key off
4 
set grid
 4 
set grid
5 
set angles degrees
 5 
set angles degrees
- 
 
 6 
set pointsize 3
6 
 
 7 
 
7 
f(x)=I0 * (sin(x))**4 + c
 8 
f(x)=I0 * (sin(x))**4 + c
8 
g(x)=I0 * (sin(x) * cos(x))**2 + a
 9 
g(x)=I0 * (sin(x) * cos(x))**2 + a
9 
 
 10 
 
10 
 
 11 
 
11 
set xlabel "fi[°]"
 12 
set xlabel "fi[°]"
12 
set ylabel "I[-]"
 13 
set ylabel "I[-]"
13 
 
 14 
 
14 
fit f(x) "./data/Malusuv_zakon/polarizace_vertikalne_stupne.txt" using 1:2 via I0, c
 15 
fit f(x) "./data/Malusuv_zakon/polarizace_vertikalne_stupne.txt" using 1:2 via I0, c
15 
set output "polarizace.png"
 16 
set output "polarizace.png"
16 
plot "./data/Malusuv_zakon/polarizace_vertikalne_stupne.txt" using 1:2 with points, f(x)
 17 
plot "./data/Malusuv_zakon/polarizace_vertikalne_stupne.txt" using 1:2 with points, f(x)
17 
 
 18 
 
18 
fit g(x) "./data/Malusuv_zakon/polarizace_horizontalne_stupne.txt" using 1:2 via I0, a
 19 
fit g(x) "./data/Malusuv_zakon/polarizace_horizontalne_stupne.txt" using 1:2 via I0, a
19 
set output "polarizace_horizontalne.png"
 20 
set output "polarizace_horizontalne.png"
20 
plot "./data/Malusuv_zakon/polarizace_horizontalne_stupne.txt" using 1:2 with points, g(x)
 21 
plot "./data/Malusuv_zakon/polarizace_horizontalne_stupne.txt" using 1:2 with points, g(x)
21 
 
 22 
 
22 
set output "podelny_rez.png"
 23 
set output "podelny_rez.png"
23 
set xlabel "x[-]"
 24 
set xlabel "x[-]"
24 
set ylabel "I[-]"
 25 
set ylabel "I[-]"
25 
plot "podelny_rez.txt" using 1:2 with points
 26 
plot "podelny_rez.txt" using 1:2 with points
26 
 
 27 
 
27 
set output "pricny_rez.png"
 28 
set output "pricny_rez.png"
28 
set xlabel "x[-]"
 29 
set xlabel "x[-]"
29 
set ylabel "I[-]"
 30 
set ylabel "I[-]"
30 
plot "./data/rozlozeni_pole/pricny_rez.txt" using 1:2 with points
 31 
plot "./data/rozlozeni_pole/pricny_rez.txt" using 1:2 with points
31 
 
 32 
 
32 
set samples 400
 33 
set samples 400
33 
set angles radians
 34 
set angles radians
34 
a=4
 35 
a=4
35 
d=100
 36 
d=100
36 
c=6
 37 
c=6
- 
 
 38 
 
37 
I0=3
 39 
I0=3
38 
 
 40 
 
39 
h(x)=I0 * sin(c*x+d) + a
 41 
h(x)=I0 * sin(c*x+d) + a
40 
 
 42 
 
41 
fit h(x) "./data/stojata_vlna/stojata_vlna_vzdalenost.txt" using 1:2 via a,d,c
 43 
fit h(x) "./data/stojata_vlna/stojata_vlna_vzdalenost.txt" using 1:2 via I0,a,d,c
42 
set output "stojata_vlna.png"
 44 
set output "stojata_vlna.png"
43 
plot "./data/stojata_vlna/stojata_vlna_vzdalenost.txt" using 1:2 with points, h(x)
 45 
plot "./data/stojata_vlna/stojata_vlna_vzdalenost.txt" using 1:2 with points, h(x)
44 
 
 46 
 
45 
 
 47 
 
46 
I0=2
 48 
I0=2
47 
 
 49 
 
48 
fit h(x) "./data/stojata_vlna/stojata_vlna_vzdalenost.txt" using 1:3 via a,d,c
 50 
fit h(x) "./data/stojata_vlna/stojata_vlna_vzdalenost.txt" using 1:3 via I0,a,d,c
49 
set output "stojata_vlna_deska.png"
 51 
set output "stojata_vlna_deska.png"
50 
plot "./data/stojata_vlna/stojata_vlna_vzdalenost.txt" using 1:3 with points, h(x)
 52 
plot "./data/stojata_vlna/stojata_vlna_vzdalenost.txt" using 1:3 with points, h(x)
51 
 
 53 
 
52 
 
 54 
 
53 
 
 55 
 
54 
set output "hrana.png"
 56 
set output "hrana.png"
55 
plot "./data/difrakce/hrana.txt" using 1:(abs($2)) with points
 57 
plot "./data/difrakce/hrana.txt" using 1:(abs($2)) with points
56 
 
 58 
 
57 
set output "pasek.png"
 59 
set output "pasek.png"
58 
plot "./data/difrakce/pasek.txt" using 1:(abs($2)) with points
 60 
plot "./data/difrakce/pasek.txt" using 1:(abs($2)) with points
59 
 
 61 
 
60 
set key under
 62 
set key under
61 
set xlabel "fi[°]"
 63 
set xlabel "fi[°]"
62 
set ylabel "I[-]"
 64 
set ylabel "I[-]"
63 
set output "sterbina.png"
 65 
set output "sterbina.png"
- 
 
 66 
 
- 
 
 67 
 
- 
 
 68 
g(x)= I1*(sin(d1*pi/l1*sin(x))/(d1*pi/l1*sin(x)))**2
- 
 
 69 
 
- 
 
 70 
sinc(x) = sin(pi*x)/(pi*x)
- 
 
 71 
i(x)=I2*(sinc(d2/l2*sin(x)))**2
- 
 
 72 
 
- 
 
 73 
set angles degrees
- 
 
 74 
 
- 
 
 75 
I1=10
- 
 
 76 
l1=0.03
- 
 
 77 
d1=0.04
- 
 
 78 
 
- 
 
 79 
I2=10
- 
 
 80 
l2=0.03
- 
 
 81 
d2=0.04
- 
 
 82 
 
- 
 
 83 
 
- 
 
 84 
fit [5:60] g(x) "sterbina.txt" using 1:2 via I1,d1,l1
- 
 
 85 
fit [5:60] i(x) "sterbina.txt" using 1:3 via I2,d2,l2
64 
plot "sterbina.txt" using 1:(abs($2)) with points title "40mm", "" using 1:(abs($3)) with points title "60mm"
 86 
plot "sterbina.txt" using 1:2 with points title "40mm", "" using 1:(abs($3)) with points title "60mm", i(x), g(x)
- 
 
 87 
 
65 
 
 88