test/plot/plotpolar.gnuplot

   1
   2 set nokey
   3 set title "Glide polar at M=1.7"
   4 set xlabel "V (m/s)"
   5 set ylabel "w (m/s)"
   6 plot 'output/results/res-polar.txt' using 2:($3) with lines, \
   7   'output/results/res-polar-17-best.txt' using 2:($3) with linespoints, \
   8   'output/results/res-polar-m.txt' using 4:(-$5) with points, \
   9   'output/results/res-polar-m.txt' using 6:(-$7) with points
  10 pause -1
  11
  12 set ylabel "LD"
  13 plot 'output/results/res-polar.txt' using 2:4 with lines, \
  14   'output/results/res-polar-17-best.txt' using 2:4 with linespoints
  15 pause -1
  16
  17 set title "Optimal solution variation with mc"
  18 set xlabel "M (m/s)"
  19 set ylabel "VbestLD (m/s)"
  20 plot 'output/results/res-polar-m.txt' using 1:2 with lines notitle
  21 pause -1
  22
  23 set xlabel "M (m/s)"
  24 set ylabel "LD"
  25 plot 'output/results/res-polar-m.txt' using 1:3 with lines notitle
  26 pause -1
  27
  28 set key
  29 set title "Glide solution variation with aircraft alt"
  30 set xlabel "h (m)"
  31 set ylabel "alt diff (m)"
  32 plot 'output/results/res-polar-h-00.txt' using 1:2 with lines title "W=0.0", \
  33      'output/results/res-polar-h-50.txt' using 1:2 with lines title "W=5.0"
  34 pause -1
  35
  36 set ylabel "time elapsed (s)"
  37 plot 'output/results/res-polar-h-00.txt' using 1:3 with lines title "W=0.0", \
  38      'output/results/res-polar-h-50.txt' using 1:3 with lines title "W=5.0"
  39 pause -1
  40
  41 set ylabel "V opt (m/s)"
  42 plot 'output/results/res-polar-h-00.txt' using 1:4 with lines title "W=0.0", \
  43      'output/results/res-polar-h-50.txt' using 1:4 with lines title "W=5.0"
  44 pause -1
  45
  46
  47 set nokey
  48 set title "Glide solution variation with wind speed"
  49
  50 set xlabel "W (m/s)"
  51 set ylabel "time elapsed (s)"
  52 plot 'output/results/res-polar-w.txt' using 5:3 with lines
  53 pause -1
  54
  55 set ylabel "V opt (m/s)"
  56 plot 'output/results/res-polar-w.txt' using 5:4 with lines
  57 pause -1
  58
  59 set nokey
  60 set title "Glide solution variation with wind angle, W=10.0 m/s"
  61
  62 set xlabel "Wangle (deg)"
  63 set ylabel "time elapsed (s)"
  64 plot 'output/results/res-polar-a.txt' using 6:3 with lines
  65 pause -1
  66
  67 set ylabel "V opt (m/s)"
  68 plot 'output/results/res-polar-a.txt' using 6:4 with lines
  69 pause -1
  70
  71 set key
  72 set title "Speed to fly"
  73 set xlabel "h (m)"
  74 plot 'output/results/res-polar-s0.txt' using 1:4 with lines title "STF", \
  75      'output/results/res-polar-s0.txt' using 1:3 with lines title "VOpt"
  76 pause -1
  77
  78 set xlabel "S gross (m/s)"
  79 plot 'output/results/res-polar-s1.txt' using 7:4 with lines title "STF", \
  80      'output/results/res-polar-s1.txt' using 7:3 with lines title "VOpt"
  81 pause -1
  82
  83 set xlabel "S gross (m/s)"
  84 plot 'output/results/res-polar-s1.txt' using 7:4 with lines title "STF (below) no wind", \
  85      'output/results/res-polar-s2.txt' using 7:4 with lines title "STF (above) no wind", \
  86      'output/results/res-polar-s3.txt' using 7:4 with lines title "STF (below) wind 10 m/s", \
  87      'output/results/res-polar-s4.txt' using 7:4 with lines title "STF (above) wind 10 m/s"
  88 pause -1
  89
  90 set xlabel "Wangle (deg)"
  91 set ylabel "bearing (deg)"
  92 plot 'output/results/res-polar-cb.txt' using 2:3 with lines title "target bearing", \
  93      'output/results/res-polar-cb.txt' using 2:4 with lines title "cruise track bearing"
  94 pause -1
  95
  96