Subversion Repositories svnkaklik

Rev

Rev 780 | Rev 969 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log

Rev 780 Rev 781
Line 3... Line 3...
3
set key off
3
set key off
4
set grid
4
set grid
5
set angles degrees
5
set angles degrees
6
 
6
 
7
f(x)=I0 * (sin(x))**4 + c
7
f(x)=I0 * (sin(x))**4 + c
-
 
8
g(x)=I0 * (sin(x) * cos(x))**2 + a
-
 
9
 
-
 
10
 
8
#g(x)=p*x+t
11
set xlabel "fi[°]"
-
 
12
set ylabel "I[-]"
9
 
13
 
10
fit f(x) "./data/polarizace_vertikalne_stupne.txt" using 1:2 via I0, c
14
fit f(x) "./data/polarizace_vertikalne_stupne.txt" using 1:2 via I0, c
11
set output "polarizace.png"
15
set output "polarizace.png"
12
set xlabel "Číslo kanálu"
-
 
13
set ylabel "Energie [keV]"
-
 
14
plot "./data/polarizace_vertikalne_stupne.txt" using 1:2 with points, f(x)
16
plot "./data/polarizace_vertikalne_stupne.txt" using 1:2 with points, f(x)
15
 
17
 
-
 
18
fit g(x) "./data/polarizace_horizontalne_stupne.txt" using 1:2 via I0, a
-
 
19
set output "polarizace_horizontalne.png"
-
 
20
plot "./data/polarizace_horizontalne_stupne.txt" using 1:2 with points, g(x)
-
 
21
 
-
 
22
set output "podelny_rez.png"
-
 
23
set xlabel "x[-]"
-
 
24
set ylabel "I[-]"
-
 
25
plot "podelny_rez.txt" using 1:2 with points
-
 
26
 
-
 
27
set samples 400
-
 
28
set angles radians
-
 
29
a=4
-
 
30
d=100
-
 
31
c=6
-
 
32
I0=3
-
 
33
 
-
 
34
h(x)=I0 * sin(c*x+d) + a
-
 
35
 
-
 
36
fit h(x) "./data/stojata_vlna_vzdalenost.txt" using 1:2 via a,d,c
-
 
37
set output "stojata_vlna.png"
-
 
38
plot "./data/stojata_vlna_vzdalenost.txt" using 1:2 with points, h(x)
-
 
39
 
-
 
40
 
-
 
41
I0=2
-
 
42
 
-
 
43
fit h(x) "./data/stojata_vlna_vzdalenost.txt" using 1:3 via a,d,c
-
 
44
set output "stojata_vlna_deska.png"
-
 
45
plot "./data/stojata_vlna_vzdalenost.txt" using 1:3 with points, h(x)
-
 
46
 
-
 
47
set output "hrana.png"
-
 
48
plot "./data/hrana.txt" using 1:(abs($2)) with points
-
 
49
 
-
 
50
set output "pasek.png"
-
 
51
plot "./data/pasek.txt" using 1:(abs($2)) with points
16
 
52
 
-
 
53
set key under
-
 
54
set xlabel "fi[°]"
-
 
55
set ylabel "I[-]"
-
 
56
set output "sterbina.png"
-
 
57
plot "sterbina.txt" using 1:(abs($2)) with points title "40mm", "" using 1:(abs($3)) with points title "60mm"
17
 
58