3,6 → 3,7 |
set key off |
set grid |
set angles degrees |
set pointsize 3 |
|
f(x)=I0 * (sin(x))**4 + c |
g(x)=I0 * (sin(x) * cos(x))**2 + a |
34,11 → 35,12 |
a=4 |
d=100 |
c=6 |
|
I0=3 |
|
h(x)=I0 * sin(c*x+d) + a |
|
fit h(x) "./data/stojata_vlna/stojata_vlna_vzdalenost.txt" using 1:2 via a,d,c |
fit h(x) "./data/stojata_vlna/stojata_vlna_vzdalenost.txt" using 1:2 via I0,a,d,c |
set output "stojata_vlna.png" |
plot "./data/stojata_vlna/stojata_vlna_vzdalenost.txt" using 1:2 with points, h(x) |
|
45,7 → 47,7 |
|
I0=2 |
|
fit h(x) "./data/stojata_vlna/stojata_vlna_vzdalenost.txt" using 1:3 via a,d,c |
fit h(x) "./data/stojata_vlna/stojata_vlna_vzdalenost.txt" using 1:3 via I0,a,d,c |
set output "stojata_vlna_deska.png" |
plot "./data/stojata_vlna/stojata_vlna_vzdalenost.txt" using 1:3 with points, h(x) |
|
61,5 → 63,26 |
set xlabel "fi[°]" |
set ylabel "I[-]" |
set output "sterbina.png" |
plot "sterbina.txt" using 1:(abs($2)) with points title "40mm", "" using 1:(abs($3)) with points title "60mm" |
|
|
g(x)= I1*(sin(d1*pi/l1*sin(x))/(d1*pi/l1*sin(x)))**2 |
|
sinc(x) = sin(pi*x)/(pi*x) |
i(x)=I2*(sinc(d2/l2*sin(x)))**2 |
|
set angles degrees |
|
I1=10 |
l1=0.03 |
d1=0.04 |
|
I2=10 |
l2=0.03 |
d2=0.04 |
|
|
fit [5:60] g(x) "sterbina.txt" using 1:2 via I1,d1,l1 |
fit [5:60] i(x) "sterbina.txt" using 1:3 via I2,d2,l2 |
plot "sterbina.txt" using 1:2 with points title "40mm", "" using 1:(abs($3)) with points title "60mm", i(x), g(x) |
|
|