#********************************************* #f(x)=a*(0.5+0.5*tanh((x-b)/c)) #******************************************** set terminal png size 450,600 set output "fig_CHtest.png" set datafile separator "," set multiplot # # f(x)=a*(0.5+0.5*tanh((x-b)/c)) set fit logfile "fit_CHtest.txt" fit f(x) "inp_CHtest.txt" using 1:2 via a,b,c # set size 1.0,0.65 set origin 0,0 set xrange [-160:20] set yrange [0:300] set xtics 20 set ytics 50 set xlabel "Testing temperature (deg.C)" set ylabel "Charpy absorbed energy (J)" set label 1 "vTrE=%.1f",b at graph 0.02,0.95 left set label 2 "USE=%.1f",a at graph 0.02,0.90 left set grid xtics ytics mxtics mytics plot f(x) linetype -1 notitle, \ "inp_CHtest.txt" using 1:2 with points pointsize 1.5 pointtype 4 linetype -1 notitle # # #********************************************* #g(x)=100.0-100.0*(0.5+0.5*tanh((x-d)/e)) #******************************************** g(x)=100.0-100.0*(0.5+0.5*tanh((x-d)/e)) set fit logfile "fit_CHtest.txt" d=-70 fit g(x) "inp_CHtest.txt" using 1:3 via d,e # set size 1.0,0.35 set origin 0,0.65 set xrange [-160:20] set yrange [0:100] set xtics 20 set ytics 25 set noxlabel set ylabel "Percent brittle fracture (%)" set label 3 "vTrs=%.1f",d at graph 0.98,0.9 right unset label 1 unset label 2 set grid xtics ytics mxtics mytics #グリッド表示 plot g(x) linetype -1 notitle, \ "inp_CHtest.txt" using 1:3 with points pointsize 1.5 pointtype 4 linetype -1 notitle # unset multiplot reset #