velocity-dist-weibull2.gnuplot (4483B)
1 load 'gnuplot/paired.pal' 2 set terminal svg size 500,500/3*2 font 'Free Serif, 10' enhanced round dynamic 3 set xtics nomirror out offset 0,0.25 4 set ytics nomirror out offset 0.5,0 5 set border 1+2 back 6 #set key top center outside Left reverse 7 unset key 8 9 set output 'build/gnuplot/velocity-dist-weibull2.svg' 10 11 # 258 1615914000 0.08700399 0.06517805 0.06370091 1.660843 12 # timestamp x_rmse y_rmse z_rmse speed 13 # 295 1616187600 0.03527403 0.08648245 0.03217789 1.637611 14 # timestamp x_rmse y_rmse z_rmse speed 15 # 399 1616958000 0.02882994 0.03539276 0.1744659 0.8539361 16 # timestamp x_rmse y_rmse z_rmse speed 17 # 292 1616166000 0.06300805 0.02371221 0.094087 2.052675 18 # [1] "MIN RMSE" 19 # timestamp x_rmse y_rmse z_rmse speed 20 # 322 1616382000 0.008555768 0.009697918 0.01404791 3.981209 21 # timestamp x_rmse y_rmse z_rmse speed 22 # 115 1614877200 0.01305041 0.009301931 0.02022811 4.356801 23 # timestamp x_rmse y_rmse z_rmse speed 24 # 39 1614330000 0.01867408 0.01648803 0.01133067 4.552804 25 26 27 set xlabel 'υ_x' offset 0,1.0 28 set yrange [0:1.2] 29 set rmargin 0 30 set tmargin 1 31 set multiplot layout 2,3 32 33 positive(x) = x<0 ? (1/0) : x 34 negative(x) = x>0 ? (1/0) : x 35 36 set lmargin 5 37 set title offset 0,-1 38 set title sprintf('NRMSE_x = %f', 0.08700399) 39 set xrange [-4:4] 40 plot \ 41 'build/velocity/x/weibull2/1615914000' index 0 using 1:2 with points pt 6 lc '#404040' title 'Observed PDF of υ_x',\ 42 '' index 0 using (positive($1)):3 with lines lw 2 lc '#c04040' title 'Estimated PDF of positive υ_x',\ 43 '' index 0 using (negative($1)):3 with lines lw 2 lc '#4040c0' title 'Estimated PDF of negative υ_x' 44 #'' index 1 using 1:2 with lines dashtype 1 lw 2 lc '#808080' notitle, \ 45 46 #set key top center outside Left reverse 47 set lmargin 5 48 #set xrange [-3:3] 49 set xlabel 'υ_y' 50 set title sprintf('NRMSE_y = %f', 0.08648245) 51 set xrange [-3:5] 52 plot \ 53 'build/velocity/y/weibull2/1616187600' index 0 using 1:2 with points pt 6 lc '#404040' title 'Observed PDF of υ_y',\ 54 '' index 0 using (positive($1)):3 with lines lw 2 lc '#c04040' title 'Estimated PDF of positive υ_y',\ 55 '' index 0 using (negative($1)):3 with lines lw 2 lc '#4040c0' title 'Estimated PDF of negative υ_y' 56 #'' index 1 using 1:2 with lines dashtype 1 lw 2 lc '#808080' notitle, \ 57 #unset key 58 59 set rmargin 2 60 set xrange [-4:4] 61 #set xtics -2,1,2 62 set xlabel 'υ_z' 63 set title sprintf('NRMSE_z = %f', 0.1744659) 64 plot \ 65 'build/velocity/z/weibull2/1616958000' index 0 using 1:2 with points pt 6 lc '#404040' title 'Observed PDF of υ_z',\ 66 '' index 0 using (positive($1)):3 with lines lw 2 lc '#c04040' title 'Estimated PDF of positive υ_z',\ 67 '' index 0 using (negative($1)):3 with lines lw 2 lc '#4040c0' title 'Estimated PDF of negative υ_z' 68 #'' index 1 using 1:2 with lines dashtype 1 lw 2 lc '#808080' notitle, \ 69 70 ################################### 71 # second row 72 ################################### 73 set bmargin 3 74 set xrange [-10:10] 75 set yrange [0:0.4] 76 set ytics 0,0.1 77 #set xtics -10,1 78 79 set lmargin 5 80 set title sprintf('NRMSE_x = %f', 0.008555768) 81 plot \ 82 'build/velocity/x/weibull2/1616382000' index 0 using 1:2 with points pt 6 lc '#404040' title 'Observed PDF of υ_x',\ 83 '' index 0 using (positive($1)):3 with lines lw 2 lc '#c04040' title 'Estimated PDF of positive υ_x',\ 84 '' index 0 using (negative($1)):3 with lines lw 2 lc '#4040c0' title 'Estimated PDF of negative υ_x' 85 #'' index 1 using 1:2 with lines dashtype 1 lw 2 lc '#808080' notitle, \ 86 87 set lmargin 5 88 #set xrange [-3:3] 89 set xlabel 'υ_y' 90 set title sprintf('NRMSE_y = %f', 0.009301931) 91 #set xrange [-3:5] 92 plot \ 93 'build/velocity/y/weibull2/1614877200' index 0 using 1:2 with points pt 6 lc '#404040' title 'Observed PDF of υ_y',\ 94 '' index 0 using (positive($1)):3 with lines lw 2 lc '#c04040' title 'Estimated PDF of positive υ_y',\ 95 '' index 0 using (negative($1)):3 with lines lw 2 lc '#4040c0' title 'Estimated PDF of negative υ_y' 96 #'' index 1 using 1:2 with lines dashtype 1 lw 2 lc '#808080' notitle, \ 97 98 set rmargin 2 99 #set xtics -2,1,2 100 set xlabel 'υ_z' 101 set title sprintf('NRMSE_z = %f', 0.01133067) 102 plot \ 103 'build/velocity/z/weibull2/1614330000' index 0 using 1:2 with points pt 6 lc '#404040' title 'Observed PDF of υ_z',\ 104 '' index 0 using (positive($1)):3 with lines lw 2 lc '#c04040' title 'Estimated PDF of positive υ_z',\ 105 '' index 0 using (negative($1)):3 with lines lw 2 lc '#4040c0' title 'Estimated PDF of negative υ_z' 106 #'' index 1 using 1:2 with lines dashtype 1 lw 2 lc '#808080' notitle, \ 107