Subversion Repositories svnkaklik

Rev

Rev 969 | Details | Compare with Previous | Last modification | View Log

Rev Author Line No. Line
779 kaklik 1 
set terminal png
2 
 
3 
set key off
4 
set grid
5 
set angles degrees
971 kaklik 6 
set pointsize 3
779 kaklik 7 
 
780 kaklik 8 
f(x)=I0 * (sin(x))**4 + c
781 kaklik 9 
g(x)=I0 * (sin(x) * cos(x))**2 + a
779 kaklik 10 
 
781 kaklik 11 
 
12 
set xlabel "fi[°]"
13 
set ylabel "I[-]"
14 
 
969 kaklik 15 
fit f(x) "./data/Malusuv_zakon/polarizace_vertikalne_stupne.txt" using 1:2 via I0, c
779 kaklik 16 
set output "polarizace.png"
969 kaklik 17 
plot "./data/Malusuv_zakon/polarizace_vertikalne_stupne.txt" using 1:2 with points, f(x)
779 kaklik 18 
 
969 kaklik 19 
fit g(x) "./data/Malusuv_zakon/polarizace_horizontalne_stupne.txt" using 1:2 via I0, a
781 kaklik 20 
set output "polarizace_horizontalne.png"
969 kaklik 21 
plot "./data/Malusuv_zakon/polarizace_horizontalne_stupne.txt" using 1:2 with points, g(x)
779 kaklik 22 
 
781 kaklik 23 
set output "podelny_rez.png"
24 
set xlabel "x[-]"
25 
set ylabel "I[-]"
26 
plot "podelny_rez.txt" using 1:2 with points
779 kaklik 27 
 
969 kaklik 28 
set output "pricny_rez.png"
29 
set xlabel "x[-]"
30 
set ylabel "I[-]"
31 
plot "./data/rozlozeni_pole/pricny_rez.txt" using 1:2 with points
32 
 
781 kaklik 33 
set samples 400
34 
set angles radians
35 
a=4
36 
d=100
37 
c=6
971 kaklik 38 
 
781 kaklik 39 
I0=3
40 
 
41 
h(x)=I0 * sin(c*x+d) + a
42 
 
971 kaklik 43 
fit h(x) "./data/stojata_vlna/stojata_vlna_vzdalenost.txt" using 1:2 via I0,a,d,c
781 kaklik 44 
set output "stojata_vlna.png"
969 kaklik 45 
plot "./data/stojata_vlna/stojata_vlna_vzdalenost.txt" using 1:2 with points, h(x)
781 kaklik 46 
 
47 
 
48 
I0=2
49 
 
971 kaklik 50 
fit h(x) "./data/stojata_vlna/stojata_vlna_vzdalenost.txt" using 1:3 via I0,a,d,c
781 kaklik 51 
set output "stojata_vlna_deska.png"
969 kaklik 52 
plot "./data/stojata_vlna/stojata_vlna_vzdalenost.txt" using 1:3 with points, h(x)
781 kaklik 53 
 
969 kaklik 54 
 
55 
 
781 kaklik 56 
set output "hrana.png"
969 kaklik 57 
plot "./data/difrakce/hrana.txt" using 1:(abs($2)) with points
781 kaklik 58 
 
59 
set output "pasek.png"
969 kaklik 60 
plot "./data/difrakce/pasek.txt" using 1:(abs($2)) with points
781 kaklik 61 
 
62 
set key under
63 
set xlabel "fi[°]"
64 
set ylabel "I[-]"
65 
set output "sterbina.png"
66 
 
971 kaklik 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
86 
plot "sterbina.txt" using 1:2 with points title "40mm", "" using 1:(abs($3)) with points title "60mm", i(x), g(x)
87 
 
88