Subversion Repositories svnkaklik

Rev

Rev 1055 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log

Rev 1055 Rev 1056
Line 1... Line 1...
1
set terminal png
1
set terminal png
2
set output "chadba_utlum.png"
2
set output "chodba_utlum.png"
3
 
3
 
4
set xlabel "vzdalenost [m]"
4
set xlabel "vzdalenost [m]"
5
set ylabel "intenzita signalu [dBm]"
5
set ylabel "intenzita signalu [dBm]"
6
 
6
 
7
#f(x) = Lo + 10*n*log(x/d0) + a
7
#f(x) = Lo + 10*n*log(x/d0) + a
8
g(x) = k*x + a
8
g(x) = k*x + a
9
q(x) = j*x + b
9
q(x) = j*x + b
10
h(x) = l*x + c
10
h(x) = l*x + c
11
 
11
 
-
 
12
 
12
#fit q(x) "chodba2" using (-30.0/126.0 * $1+30.0):($3) via j,b
13
fit g(x) "chodba1" using (-30.0/151.0 * $1+31.0):($3) via k,a
13
#fit h(x) "chodba3" using (-30.0/111.0 * $1+30.0):($3) via l,c
14
#fit h(x) "chodba3" using (-30.0/111.0 * $1+30.0):($3) via l,c
14
 
15
 
-
 
16
set xrange [1:50]
15
 
17
 
-
 
18
l=k
16
fit g(x) "chodba1" using (-30.0/151.0 * $1+31.0):($3), "chodba2" using (-30.0/126.0 * $1+30.0):($3) via k,a
19
c=a - 3*0.65
17
 
20
 
18
 
21
 
19
plot "chodba1" using (-30.0/151.0 * $1+31.0):($3) with points title "1.pokus", "chodba2" using (-30.0/126.0 * $1+30.0):($3) with points title "2.pokus", "chodba3" using (-30.0/111.0 * $1+30.0):($3) with points title "3.pokus", g(x), q(x), h(x)
22
plot "chodba1" using (-30.0/151.0 * $1+31.0):($3) with points title "1.pokus", "chodba2" using (-30.0/126.0 * $1+30.0):($3) with points title "2.pokus", "chodba3" using (-30.0/111.0 * $1+30.0):($3) with points title "3.pokus", g(x) title "mean", h(x) title "90% probability treshold" 
20
 
23
 
21
 
24
 
22
set output "chodba_braun.png"
25
set autoscale x
23
plot "Chodba_s_lidmi.dat" using (-30.0/151.0 * $1+30.0):($3) with points notitle 
-
 
24
 
26
 
-
 
27
fit q(x) "mistnost1" using (-12.0/58.0 * $1+17.0):($3) via j,b
25
 
28
 
26
set output "mistnost.png"
29
set output "mistnost.png"
27
plot "mistnost1" using (-12.0/58.0 * $1+17.0):($3) with points notitle 
30
plot "mistnost1" using (-12.0/58.0 * $1+17.0):($3) with points notitle, q(x) 
-
 
31
 
-
 
32
 
-
 
33
set xlabel "cas [s]"
-
 
34
set ylabel "intenzita signalu [dBm]"
-
 
35
 
-
 
36
bin(x,size) = floor(1.0*x/size)*size
-
 
37
 
-
 
38
set output "chodba_braun.png"
-
 
39
plot "Chodba_s_lidmi.dat" using (0.2 * $1):($3) with points notitle 
-
 
40
set output "chodba_braun_rozdeleni.png"
-
 
41
plot "Chodba_s_lidmi.dat" using (bin($3, 1.0)):($3) with boxes notitle 
-
 
42
 
-
 
43
 
28
 
44