3479 |
pomiceva |
1 |
set terminal png size 800,640 |
|
|
2 |
|
|
|
3 |
set ylabel "alt [m]" |
|
|
4 |
set xlabel "time [min]" |
|
|
5 |
#set xdata time |
|
|
6 |
#set timefmt "%s" |
|
|
7 |
#set format x "%H:%M:%S" |
|
|
8 |
set output "KorekceTlaku.png" |
|
|
9 |
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 |
|
|
10 |
fit f1(x) "mereni_tlaku_16122013.log" using ($10):($17*100000) via k1, q1 |
|
|
11 |
fit f2(x) "mereni_tlaku_16122013.log" using ($11):($17*100000) via k2, q2 |
|
|
12 |
fit f3(x) "mereni_tlaku_16122013.log" using ($12):($17*100000) via k3, q3 |
|
|
13 |
fit f4(x) "mereni_tlaku_16122013.log" using ($14):($17*100000) via k4, q4 |
|
|
14 |
fit f5(x) "mereni_tlaku_16122013.log" using ($15):($17*100000) via k5, q5 |
|
|
15 |
T0=288.15; tau=-6.5e-3; R=29.271; p0=101325 |
|
|
16 |
h1(x)=(T0/tau)*((p0/f1(x))**(R*tau)-1) |
|
|
17 |
h2(x)=(T0/tau)*((p0/f2(x))**(R*tau)-1) |
|
|
18 |
h3(x)=(T0/tau)*((p0/f3(x))**(R*tau)-1) |
|
|
19 |
h4(x)=(T0/tau)*((p0/f4(x))**(R*tau)-1) |
|
|
20 |
h5(x)=(T0/tau)*((p0/f5(x))**(R*tau)-1) |
|
|
21 |
|
|
|
22 |
#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) |
|
|
23 |
|
|
|
24 |
set output "Stoupani.png" |
|
|
25 |
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" |