Rev Author Line No. Line
3479 pomiceva 1 set terminal png size 800,640
2  
3 set ylabel "alt [m]"
4 set xlabel "time [min]"
3482 pomiceva 5 #set output "KorekceTlaku.png"
3487 pomiceva 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_s_DPI145_16122013.log" using ($10):(($17+$18)*50000) via k1, q1
8 fit f2(x) "mereni_tlaku_s_DPI145_16122013.log" using ($11):(($17+$18)*50000) via k2, q2
9 fit f3(x) "mereni_tlaku_s_DPI145_16122013.log" using ($12):(($17+$18)*50000) via k3, q3
10 fit f4(x) "mereni_tlaku_s_DPI145_16122013.log" using ($14):(($17+$18)*50000) via k4, q4
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
3479 pomiceva 13 T0=288.15; tau=-6.5e-3; R=29.271; p0=101325
14 h1(x)=(T0/tau)*((p0/f1(x))**(R*tau)-1)
15 h2(x)=(T0/tau)*((p0/f2(x))**(R*tau)-1)
16 h3(x)=(T0/tau)*((p0/f3(x))**(R*tau)-1)
17 h4(x)=(T0/tau)*((p0/f4(x))**(R*tau)-1)
18 h5(x)=(T0/tau)*((p0/f5(x))**(R*tau)-1)
3487 pomiceva 19 #h6(x)=(T0/tau)*((p0/f6(x))**(R*tau)-1)
3479 pomiceva 20  
3487 pomiceva 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)
3479 pomiceva 22  
3487 pomiceva 23 ###tlak-vyska
3482 pomiceva 24 set output "Tlak_vyska.png"
25 set ylabel "alt [m]"
3487 pomiceva 26 set xlabel "press [kPa]"
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"
3482 pomiceva 28  
29  
30  
31  
3487 pomiceva 32 ###stoupani: cas-vyska
3479 pomiceva 33 set output "Stoupani.png"
3481 kaklik 34 set xdata time
35 set timefmt "%s"
36 set format x "%H:%M:%S"
3482 pomiceva 37 set key under
3487 pomiceva 38 set xlabel "time"
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"
3482 pomiceva 40  
41