| 851 |
kaklik |
1 |
set table 'tohleprepishned.txt'
|
|
|
2 |
splot 'tohleprepishned.avs' binary format='%4uint8' record=180x180 skip=8 using 2
|
|
|
3 |
unset table
|
|
|
4 |
|
|
|
5 |
set xlabel 'pixels'
|
|
|
6 |
set ylabel 'pixels'
|
|
|
7 |
set zlabel 'intensity' rotate by 90;
|
|
|
8 |
set hidden3d
|
|
|
9 |
unset key
|
|
|
10 |
set terminal png
|
|
|
11 |
set output 'tohleprepishned_3d.png'
|
|
|
12 |
splot 'tohleprepishned.txt' with lines
|
|
|
13 |
unset output
|
|
|
14 |
|
|
|
15 |
unset key
|
|
|
16 |
set samples 5000
|
|
|
17 |
set xlabel 'pixels'
|
|
|
18 |
set ylabel 'intensity'
|
|
|
19 |
f(x)=a*exp(-(x-x0)**2/(2*w**2))
|
|
|
20 |
a=100
|
|
|
21 |
x0=90
|
|
|
22 |
w=30
|
|
|
23 |
fit f(x) 'tohleprepishned.txt' every ::90::100 using 2:3 via a,x0,w
|
|
|
24 |
set output 'tohleprepishned_fit1.png'
|
|
|
25 |
plot f(x) w lines linewidth 3, 'tohleprepishned.txt' every ::90::100 using 2:3
|
|
|
26 |
|
|
|
27 |
|
|
|
28 |
unset output
|
|
|
29 |
unset key
|
|
|
30 |
f(x)=a*exp(-(x-x0)**2/(2*w**2))
|
|
|
31 |
a=200
|
|
|
32 |
x0=80
|
|
|
33 |
w=20
|
|
|
34 |
fit f(x) 'tohleprepishned.txt' every :::90::100 using 1:3 via a,x0,w
|
|
|
35 |
set output 'tohleprepishned_fit2.png'
|
|
|
36 |
plot f(x) w lines linewidth 3, 'tohleprepishned.txt' every :::90::100 using 1:3
|
|
|
37 |
unset output
|
|
|
38 |
|