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