Subversion Repositories svnkaklik

Rev

Details | Last modification | View Log

Rev Author Line No. Line
851 kaklik 1 
set table 'gauss2+.txt'
2 
splot 'gauss2+.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 'gauss2+_3d.png'
12 
splot 'gauss2+.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) 'gauss2+.txt' every ::90::100 using 2:3 via a,x0,w
24 
set output 'gauss2+_fit1.png'
25 
plot f(x) w lines linewidth 3, 'gauss2+.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) 'gauss2+.txt' every :::90::100 using 1:3 via a,x0,w
35 
set output 'gauss2+_fit2.png'
36 
plot f(x) w lines linewidth 3, 'gauss2+.txt' every :::90::100 using 1:3
37 
unset output
38