/Modules/Sensors/ALTIMET01A/DOC/SRC/mereni/plot.gp
1,26 → 1,41
#set term svg size 800,640 #enhanced mouse
set terminal png size 800,640
#set samples 1000
set xlabel "Time [s]"
set ylabel "T [deg C]"
set y2label "P [pres Pa]"
 
set ylabel "press DPI145 [Pa]"
set xlabel "press cidla [Pa]"
set output "KorekceTlaku.png"
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
fit f1(x) "mereni_tlaku_16122013.log" using ($10):($17*100000) via k1, q1
fit f2(x) "mereni_tlaku_16122013.log" using ($11):($17*100000) via k2, q2
fit f3(x) "mereni_tlaku_16122013.log" using ($12):($17*100000) via k3, q3
fit f4(x) "mereni_tlaku_16122013.log" using ($14):($17*100000) via k4, q4
fit f5(x) "mereni_tlaku_16122013.log" using ($15):($17*100000) via k5, q5
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)
 
 
set xlabel "time [s]"
 
#set y2label "P [pres Pa]"
set key under
#set grid on
set xdata time
set timefmt "%s"
set format x "%M:%S"
set logscale y2
set y2tics
#set logscale y2
#set y2tics
 
set y2range [3000:95000]
#set y2range [3000:95000]
#set xrange [500:1500]
set output "Temp.png"
g(x)=k*x+q*x**2+a
fit g(x) "mereni_tlaku_02122013.log" using 1:10 via k, q, a
set output "Press.png"
 
 
plot "mereni_tlaku_02122013.log" using 1:3 with linespoints title "T1", "" using 1:4 with linespoints title "T2", "" using 1:5 with linespoints title "T3", "" using 1:6 with linespoints title "T4", "" using 1:7 with linespoints title "T5", "" using 1:8 with linespoints title "T6", "" using 1:9 with linespoints title "TRef", "" using 1:10 axes x1y2 with linespoints title "P1", "" using 1:11 axes x1y2 with linespoints title "P2", "" using 1:12 axes x1y2 with linespoints title "P3", "" using 1:13 axes x1y2 with linespoints title "P4", "" using 1:14 axes x1y2 with linespoints title "P5", "" using 1:15 axes x1y2 with linespoints title "P6", "" using 1:16 axes x1y2 with linespoints title "PRef", "" using 1:($17*100000) axes x1y2 with linespoints title "PDPI145", g(x)
 
#plot "mereni_tlaku_16122013.log" using 1:3 with linespoints title "T1", "" using 1:4 with linespoints title "T2", "" using 1:5 with linespoints title "T3", "" using 1:6 with linespoints title "T4", "" using 1:7 with linespoints title "T5", "" using 1:8 with linespoints title "T6", "" using 1:9 with linespoints title "TRef", "" using 1:10 axes x1y2 with linespoints title "P1", "" using 1:11 axes x1y2 with linespoints title "P2", "" using 1:12 axes x1y2 with linespoints title "P3", "" using 1:13 axes x1y2 with linespoints title "P4", "" using 1:14 axes x1y2 with linespoints title "P5", "" using 1:15 axes x1y2 with linespoints title "P6", "" using 1:16 axes x1y2 with linespoints title "PRef", "" using 1:($17*100000) axes x1y2 with linespoints title "PDPI145"
 
plot "mereni_tlaku_16122013.log" using 1:10 axes x1y1 with linespoints title "P1", "" using 1:11 axes x1y1 with linespoints title "P2", "" using 1:12 axes x1y1 with linespoints title "P3", "" using 1:14 axes x1y1 with linespoints title "P5", "" using 1:15 axes x1y1 with linespoints title "P6", "" using 1:($17*100000) axes x1y1 with linespoints title "PDPI145"
 
set output "Pdif.png"
plot "mereni_tlaku_02122013.log" using 1:($10-$17*100000) with linespoints title "P1", "" using 1:($11-$17*100000) with linespoints title "P2", "" using 1:($12-$17*100000) with linespoints title "P3", "" using 1:($13-$17*100000) with linespoints title "P4", "" using 1:($14-$17*100000) with linespoints title "P5", "" using 1:($15-$17*100000) with linespoints title "P6"
plot "mereni_tlaku_16122013.log" using 1:($10-$17*100000) with linespoints title "P1", "" using 1:($11-$17*100000) with linespoints title "P2", "" using 1:($12-$17*100000) with linespoints title "P3", "" using 1:($14-$17*100000) with linespoints title "P5", "" using 1:($15-$17*100000) with linespoints title "P6"
 
 
/Modules/Sensors/ALTIMET01A/DOC/SRC/mereni/vyska
0,0 → 1,25
set terminal png size 800,640
 
set ylabel "alt [m]"
set xlabel "time [min]"
#set xdata time
#set timefmt "%s"
#set format x "%H:%M:%S"
set output "KorekceTlaku.png"
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
fit f1(x) "mereni_tlaku_16122013.log" using ($10):($17*100000) via k1, q1
fit f2(x) "mereni_tlaku_16122013.log" using ($11):($17*100000) via k2, q2
fit f3(x) "mereni_tlaku_16122013.log" using ($12):($17*100000) via k3, q3
fit f4(x) "mereni_tlaku_16122013.log" using ($14):($17*100000) via k4, q4
fit f5(x) "mereni_tlaku_16122013.log" using ($15):($17*100000) via k5, q5
T0=288.15; tau=-6.5e-3; R=29.271; p0=101325
h1(x)=(T0/tau)*((p0/f1(x))**(R*tau)-1)
h2(x)=(T0/tau)*((p0/f2(x))**(R*tau)-1)
h3(x)=(T0/tau)*((p0/f3(x))**(R*tau)-1)
h4(x)=(T0/tau)*((p0/f4(x))**(R*tau)-1)
h5(x)=(T0/tau)*((p0/f5(x))**(R*tau)-1)
 
#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)
 
set output "Stoupani.png"
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"