Subversion Repositories svnkaklik

Compare Revisions

Ignore whitespace Rev 970 → Rev 971

/dokumenty/skolni/PRA2/Mikrovlny/plot.gp
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)