iccsa-21-wind

Wind Simulation Using High-Frequency Velocity Component Measurements
git clone https://git.igankevich.com/iccsa-21-wind.git
Log | Files | Refs

hold-peak.gnuplot (1763B)


      1 set terminal svg size 450,450/2 font 'Free Serif, 10' enhanced round dynamic
      2 set xtics nomirror out offset 0,0.5
      3 set ytics nomirror out offset 0.5,0
      4 set border 1+2 back
      5 unset key
      6 
      7 set output 'build/gnuplot/hold-peak.svg'
      8 
      9 set xlabel offset 0,1.0
     10 set xlabel 'Time, s'
     11 set ylabel 'Wind speed, m/s' offset 1.5,0
     12 set title offset 0,-1.5
     13 set lmargin 5
     14 set rmargin 2
     15 set tmargin 4.0
     16 set multiplot layout 1,2
     17 
     18 round(x) = floor(x+0.5)
     19 filename(timestamp) = sprintf('build/hold-peak/%d', timestamp)
     20 
     21 load "build/gnuplot/hold_peak_rmse.gnuplot"
     22 
     23 set yrange [0:5]
     24 set ytics 0,1,5
     25 set xtics 0,30
     26 
     27 set title sprintf("%s–%s\nNRMSE = %.1f%%",\
     28         strftime('%F %H:%M', speed_rmse_max_timestamp),\
     29         strftime('%H:%M', speed_rmse_max_timestamp+2*60*60),\
     30         (speed_rmse_max*100))
     31 plot \
     32 filename(speed_rmse_max_timestamp) index 0 using 1:2 with points pt 6 lc '#404040',\
     33 '' index 0 using 1:3 with lines lw 2 lc '#c04040',\
     34 '' index 1 using 1:2 with lines lw 2 lc '#c0c040'
     35 
     36 set key Left reverse at screen 0.80, screen 1.02
     37 
     38 set title sprintf("%s–%s\nNRMSE = %.1f%%",\
     39         strftime('%F %H:%M', speed_rmse_min_timestamp),\
     40         strftime('%H:%M', speed_rmse_min_timestamp+2*60*60),\
     41         (speed_rmse_min*100))
     42 plot \
     43 filename(speed_rmse_min_timestamp) index 0 using 1:2 with points pt 6 lc '#404040' title 'Load cell data',\
     44 '' index 0 using 1:3 with lines lw 2 lc '#c04040' title 'HP-866A data',\
     45 '' index 1 using 1:2 with lines lw 2 lc '#c0c040' title 'Smoothed load cell data'
     46 
     47 # set title sprintf('NRMSE_x = %.1f%%', (speed_rmse_longest*100))
     48 # plot \
     49 # filename(speed_rmse_longest_timestamp) index 0 using 1:2 with points pt 6 lc '#404040',\
     50 # '' index 0 using 1:3 with lines lw 2 lc '#c04040',\
     51 # '' index 1 using 1:2 with lines lw 2 lc '#40c040'