Subversion Repositories svnkaklik

Rev

Go to most recent revision | Details | Last modification | View Log

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