Subversion Repositories svnkaklik

Rev

Rev 969 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log

Rev 969 Rev 971
Line 1... Line 1...
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 
 
Line 32... Line 33...
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 
 
Line 59... Line 61...
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