| 1 |
set terminal png size 800,640 |
1 |
set terminal png size 800,640 |
| 2 |
|
2 |
|
| 3 |
set ylabel "alt [m]" |
3 |
set ylabel "alt [m]" |
| 4 |
set xlabel "time [min]" |
4 |
set xlabel "time [min]" |
| 5 |
#set output "KorekceTlaku.png" |
5 |
#set output "KorekceTlaku.png" |
| 6 |
#f1(x)=k1*x+q1; f2(x)=k2*x+q2; f3(x)=k3*x+q3; f4(x)=k4*x+q4; f5(x)=k5*x+q5 |
6 |
f1(x)=k1*x+q1; f2(x)=k2*x+q2; f3(x)=k3*x+q3; f4(x)=k4*x+q4; f5(x)=k5*x+q5; f6(x)=k6*x*q6 |
| 7 |
#fit f1(x) "mereni_tlaku_16122013.log" using ($10):($17*100000) via k1, q1 |
7 |
fit f1(x) "mereni_tlaku_s_DPI145_16122013.log" using ($10):(($17+$18)*50000) via k1, q1 |
| 8 |
#fit f2(x) "mereni_tlaku_16122013.log" using ($11):($17*100000) via k2, q2 |
8 |
fit f2(x) "mereni_tlaku_s_DPI145_16122013.log" using ($11):(($17+$18)*50000) via k2, q2 |
| 9 |
#fit f3(x) "mereni_tlaku_16122013.log" using ($12):($17*100000) via k3, q3 |
9 |
fit f3(x) "mereni_tlaku_s_DPI145_16122013.log" using ($12):(($17+$18)*50000) via k3, q3 |
| 10 |
#fit f4(x) "mereni_tlaku_16122013.log" using ($14):($17*100000) via k4, q4 |
10 |
fit f4(x) "mereni_tlaku_s_DPI145_16122013.log" using ($14):(($17+$18)*50000) via k4, q4 |
| 11 |
#fit f5(x) "mereni_tlaku_16122013.log" using ($15):($17*100000) via k5, q5 |
11 |
fit f5(x) "mereni_tlaku_s_DPI145_16122013.log" using ($15):(($17+$18)*50000) via k5, q5 |
| - |
|
12 |
#fit f6(x) "mereni_tlaku_s_DPI145_16122013.log" using ($15):(($17+$18)*50000) via k6, q6 |
| 12 |
T0=288.15; tau=-6.5e-3; R=29.271; p0=101325 |
13 |
T0=288.15; tau=-6.5e-3; R=29.271; p0=101325 |
| 13 |
h1(x)=(T0/tau)*((p0/f1(x))**(R*tau)-1) |
14 |
h1(x)=(T0/tau)*((p0/f1(x))**(R*tau)-1) |
| 14 |
h2(x)=(T0/tau)*((p0/f2(x))**(R*tau)-1) |
15 |
h2(x)=(T0/tau)*((p0/f2(x))**(R*tau)-1) |
| 15 |
h3(x)=(T0/tau)*((p0/f3(x))**(R*tau)-1) |
16 |
h3(x)=(T0/tau)*((p0/f3(x))**(R*tau)-1) |
| 16 |
h4(x)=(T0/tau)*((p0/f4(x))**(R*tau)-1) |
17 |
h4(x)=(T0/tau)*((p0/f4(x))**(R*tau)-1) |
| 17 |
h5(x)=(T0/tau)*((p0/f5(x))**(R*tau)-1) |
18 |
h5(x)=(T0/tau)*((p0/f5(x))**(R*tau)-1) |
| - |
|
19 |
#h6(x)=(T0/tau)*((p0/f6(x))**(R*tau)-1) |
| 18 |
|
20 |
|
| 19 |
#plot "mereni_tlaku_16122013.log" using ($10):($17*100000) with linespoints title "P1", f1(x), "" using ($11):($17*100000) axes x1y1 with linespoints title "P2", f2(x), "" using ($12):($17*100000) axes x1y1 with linespoints title "P3", f3(x), "" using ($14):($17*100000) axes x1y1 with linespoints title "P4", f4(x), "" using ($15):($17*100000) axes x1y1 with linespoints title "P5", f5(x) |
21 |
#plot "mereni_tlaku_16122013.log" using ($10):(($17+$18)*50000) with linespoints title "P1", f1(x), "" using ($11):(($17+$18)*50000) axes x1y1 with linespoints title "P2", f2(x), "" using ($12):($($17+$18)*50000) axes x1y1 with linespoints title "P3", f3(x), "" using ($14):(($17+$18)*50000) axes x1y1 with linespoints title "P4", f4(x), "" using ($15):(($17+$18)*50000) axes x1y1 with linespoints title "P5", f5(x) |
| 20 |
|
22 |
|
| 21 |
a=1 |
- |
|
| 22 |
b=1 |
- |
|
| 23 |
c=1 |
- |
|
| 24 |
|
- |
|
| 25 |
f1(x)=a*x+b*x**2+c |
- |
|
| 26 |
fit f1(x) "mereni_tlaku_16122013.log" using ($10/1000):(h1($10)) via a, b, c |
- |
|
| 27 |
#fit f2(x) "mereni_tlaku_16122013.log" using ($11):($17*100000) via k2, q2 |
- |
|
| 28 |
#fit f3(x) "mereni_tlaku_16122013.log" using ($12):($17*100000) via k3, q3 |
- |
|
| 29 |
#fit f4(x) "mereni_tlaku_16122013.log" using ($14):($17*100000) via k4, q4 |
- |
|
| 30 |
#fit f5(x) "mereni_tlaku_16122013.log" using ($15):($17*100000) via k5, q5 |
- |
|
| 31 |
#tlak-vyska |
23 |
###tlak-vyska |
| 32 |
set output "Tlak_vyska.png" |
24 |
set output "Tlak_vyska.png" |
| 33 |
set ylabel "alt [m]" |
25 |
set ylabel "alt [m]" |
| 34 |
set xlabel "press [kPa]" |
26 |
set xlabel "press [kPa]" |
| 35 |
plot "mereni_tlaku_16122013.log" using ($10/1000):(h1($10)) with linespoints title "s1", "" using ($11/1000):(h2($11)) with linespoints title "s2", "" using ($12/1000):(h3($12)) with linespoints title "s3", "" using ($14/1000):(h4($14)) with linespoints title "s4", "" using ($15/1000):(h5($15)) with linespoints title "s5", f1(x) |
27 |
#plot "mereni_tlaku_s_DPI145_16122013.log" using ($10/1000):(h1($10)) with linespoints title "s1", "" using ($11/1000):(h2($11)) with linespoints title "s2", "" using ($12/1000):(h3($12)) with linespoints title "s3", "" using ($13/1000):(h4($13)) with linespoints title "s4", "" using ($14/1000):(h5($14)) with linespoints title "s5", "" using ($15/1000):(h6($15)) with linespoints title "s6" |
| 36 |
|
28 |
|
| 37 |
|
29 |
|
| 38 |
|
30 |
|
| 39 |
|
31 |
|
| 40 |
#cas-vyska |
32 |
###stoupani: cas-vyska |
| 41 |
set output "Stoupani.png" |
33 |
set output "Stoupani.png" |
| 42 |
set xdata time |
34 |
set xdata time |
| 43 |
set timefmt "%s" |
35 |
set timefmt "%s" |
| 44 |
set format x "%H:%M:%S" |
36 |
set format x "%H:%M:%S" |
| 45 |
set key under |
37 |
set key under |
| - |
|
38 |
set xlabel "time" |
| 46 |
plot "mereni_tlaku_16122013.log" using 1:(h1($10)) with linespoints title "h1", "" using 1:(h2($11)) with linespoints title "h2", "" using 1:(h3($12)) with linespoints title "h3", "" using 1:(h4($14)) with linespoints title "h4", "" using 1:(h5($15)) with linespoints title "h5" |
39 |
plot "mereni_tlaku_16122013.log" using 1:(h1($10)) with linespoints title "s1", "" using 1:(h2($11)) with linespoints title "s2", "" using 1:(h3($12)) with linespoints title "s3", "" using 1:(h4($14)) with linespoints title "s4", "" using 1:(h5($15)) with linespoints title "s5" |
| 47 |
|
40 |
|
| 48 |
|
41 |
|