iccsa-21-wind

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

acf.gnuplot (4113B)


      1 # ACF X min. RMSE = 0.003198 (1614920400), mean speed = 3.230443
      2 # ACF Y min. RMSE = 0.004612 (1615388400), mean speed = 1.435405
      3 # ACF Z min. RMSE = 0.001059 (1614135600), mean speed = 2.748603
      4 # ACF X max. RMSE = 0.178318 (1616058000), mean speed = 3.347379
      5 # ACF Y max. RMSE = 0.222642 (1614430800), mean speed = 3.308343
      6 # ACF Z max. RMSE = 0.058600 (1614452400), mean speed = 2.643957
      7 
      8 set terminal svg size 450,450/3*2 font 'Free Serif, 10' enhanced round dynamic
      9 set xrange [-40:40]
     10 set xtics 0,20,40 nomirror out offset 0,0.5 add ('20' -20, '40' -40)
     11 set ytics nomirror out offset 0.5,0
     12 set border 1+2 back
     13 #set key top center outside Left reverse
     14 unset key
     15 
     16 set output 'build/gnuplot/acf.svg'
     17 
     18 set xlabel offset 0,1.0
     19 set xlabel 'Lag'
     20 set title offset 0,-1
     21 set lmargin 5
     22 set rmargin 0
     23 set tmargin 2
     24 set multiplot layout 2,3
     25 
     26 round(x) = floor(x+0.5)
     27 filename(axis,timestamp) = sprintf('build/acf/%s/%d', axis, timestamp)
     28 speed(x) = sqrt(x)
     29 positive(x,y) = x>=0 ? y : NaN
     30 positive(x,y) = x<=0 ? y : NaN
     31 
     32 load "build/gnuplot/acf_rmse.gnuplot"
     33 
     34 ###################################
     35 # first row
     36 ###################################
     37 set yrange [0:2]
     38 
     39 set title sprintf("%s–%s\nNRMSE_x = %.1f%%",\
     40     strftime('%F %H:%M', acf_x_rmse_max_timestamp),\
     41     strftime('%H:%M', acf_x_rmse_max_timestamp+2*60*60),\
     42     (acf_x_rmse_max*100))
     43 plot \
     44 filename('x',acf_x_rmse_max_timestamp) index 0 using (-$1):2 with points pt 6 lc '#404040',\
     45 '' index 0 using (-$1):3 with lines lw 2 lc '#4040c0',\
     46 '' index 1 using ($1):2 with points pt 6 lc '#404040',\
     47 '' index 1 using ($1):3 with lines lw 2 lc '#c04040'
     48 
     49 set yrange [0:0.25]
     50 set title sprintf("%s–%s\nNRMSE_y = %.1f%%",\
     51         strftime('%F %H:%M', acf_y_rmse_max_timestamp),\
     52         strftime('%H:%M', acf_y_rmse_max_timestamp+2*60*60),\
     53         (acf_y_rmse_max*100))
     54 plot \
     55 filename('y',acf_y_rmse_max_timestamp) index 0 using (-$1):2 with points pt 6 lc '#404040',\
     56 '' index 0 using (-$1):3 with lines lw 2 lc '#4040c0',\
     57 '' index 1 using ($1):2 with points pt 6 lc '#404040',\
     58 '' index 1 using ($1):3 with lines lw 2 lc '#c04040'
     59 
     60 set yrange [0:2]
     61 set title sprintf("%s–%s\nNRMSE_z = %.1f%%",\
     62         strftime('%F %H:%M', acf_z_rmse_max_timestamp),\
     63         strftime('%H:%M', acf_z_rmse_max_timestamp+2*60*60),\
     64         (acf_z_rmse_max*100))
     65 plot \
     66 filename('z',acf_z_rmse_max_timestamp) index 0 using (-$1):2 with points pt 6 lc '#404040',\
     67 '' index 0 using (-$1):3 with lines lw 2 lc '#4040c0',\
     68 '' index 1 using ($1):2 with points pt 6 lc '#404040',\
     69 '' index 1 using ($1):3 with lines lw 2 lc '#c04040'
     70 
     71 ###################################
     72 # second row
     73 ###################################
     74 set bmargin 3
     75 set yrange [0:0.25]
     76 
     77 set title sprintf("%s–%s\nNRMSE_x = %.1f%%",\
     78         strftime('%F %H:%M', acf_x_rmse_min_timestamp),\
     79         strftime('%H:%M', acf_x_rmse_min_timestamp+2*60*60),\
     80         (acf_x_rmse_min*100))
     81 plot \
     82 filename('x',acf_x_rmse_min_timestamp) index 0 using (-$1):2 with points pt 6 lc '#404040',\
     83 '' index 0 using (-$1):3 with lines lw 2 lc '#4040c0',\
     84 '' index 1 using ($1):2 with points pt 6 lc '#404040',\
     85 '' index 1 using ($1):3 with lines lw 2 lc '#c04040'
     86 
     87 set yrange [0:0.25]
     88 set title sprintf("%s–%s\nNRMSE_y = %.1f%%",\
     89         strftime('%F %H:%M', acf_y_rmse_min_timestamp),\
     90         strftime('%H:%M', acf_y_rmse_min_timestamp+2*60*60),\
     91         (acf_y_rmse_min*100))
     92 plot \
     93 filename('y',acf_y_rmse_min_timestamp) index 0 using (-$1):2 with points pt 6 lc '#404040',\
     94 '' index 0 using (-$1):3 with lines lw 2 lc '#4040c0',\
     95 '' index 1 using ($1):2 with points pt 6 lc '#404040',\
     96 '' index 1 using ($1):3 with lines lw 2 lc '#c04040'
     97 
     98 set yrange [0:0.1]
     99 set title sprintf("%s–%s\nNRMSE_z = %.1f%%",\
    100         strftime('%F %H:%M', acf_z_rmse_min_timestamp),\
    101         strftime('%H:%M', acf_z_rmse_min_timestamp+2*60*60),\
    102         (acf_z_rmse_min*100))
    103 plot \
    104 filename('z',acf_z_rmse_min_timestamp) index 0 using (-$1):2 with points pt 6 lc '#404040',\
    105 '' index 0 using (-$1):3 with lines lw 2 lc '#4040c0',\
    106 '' index 1 using ($1):2 with points pt 6 lc '#404040',\
    107 '' index 1 using ($1):3 with lines lw 2 lc '#c04040'