gnuplot画图初步探究和命令讲解(处理变密度 的reactive chemistry 得到的算例结果)

原理:命令含义

  1. plot ‘POST_Xs_D.dat’
    根据数据文件绘图

  2. #
    注释

  3. u 1:2
    using 的简写,表示使用第一列和第二列的数据

  4. w lp
    with linespoints的简写,表示使用点线图

  5. set title ‘mass frac of 4 species’
    设置图片的标题

  6. ,
    如果在同一个图画出多条曲线,要用逗号分割,除了plot其他的都要重写一遍

  7. \
    换行 表示下行实际上和上一行是同一行

  8. t ‘XD DEM’
    title 的简写,这里的title其实是图例legend的意思

  9. set key bottom
    设置图例放在右下角

  10. set multiplot layout 2,2
    设置子图

  11. unset multiplot
    取消设置子图,如果不写有可能会发生错误

  12. set xlabel ‘time’
    设置X轴标签

  13. set xrange [0:2.1]
    设置X轴区间

设置双Y轴的方法

set y2range [0:6e-5] #第二个Y轴的坐标范围
set y2tics#必须写,否则第二个Y轴无刻度显示
set ytics nomirror#必须写,否则第二个Y轴刻度会和第一个Y轴的刻度重叠显示
set y2label 'relative err'#必须写,第二个Y轴的图例

在plot的选项中加入

axis x1y2

表示这列数据是绘制在第二个Y轴上的

操作方法

  1. 终端输入gnuplot
  2. load 'gnuplot’画出图片并保存
  3. load 'gnuplot_pmass’画出图片并保存

脚本代码

脚本1:名字gnuplot

set multiplot layout 2,2
#set key bottom
set xrange [0:2.1]
set yrange [0:1]
set y2range [0:6e-5]
set y2tics
set ytics nomirror
set y2label 'relative err'

#DBI
set title 'mass frac of 4 species'
set xlabel 'time'
set ylabel 'mass frac'
plot 'POST_Xs_D.dat'  u 1:3 w lp lw 2 t 'XD DEM', \
'POST_Xs_B.dat'  u 1:3 w lp lw 2 dt 2 t 'XB DEM', \
'POST_Xs_I.dat'  u 1:3 w lp  t 'XI DEM'

#D
set title 'mass frac of D'
set xlabel 'time'
set ylabel 'Xs_D'
plot 'POST_Xs_D.dat'  u 1:2 w lp lw 2 t 'Xs_D', \
     'POST_Xs_D.dat'  u 1:3 w lp lw 2 dt 2 t 'Xs_D DEM', \
     'POST_Xs_D.dat'  u 1:4 w p pt 7 lc  'red' ps 1 axis x1y2 t 'abs err'

#B
set title 'mass frac of B'
set xlabel 'time'
set ylabel 'Xs_B'
plot 'POST_Xs_B.dat'  u 1:2 w lp lw 2 t 'Xs_B', \
     'POST_Xs_B.dat'  u 1:3 w lp lw 2 dt 2 t 'Xs_B DEM', \
     'POST_Xs_B.dat'  u 1:4 w p pt 7 lc  'red' ps 1 axis x1y2 t 'abs err'

#I
set title 'mass frac of I'
set xlabel 'time'
set ylabel 'Xs_I'
plot 'POST_Xs_I.dat'  u 1:2 w lp lw 2 t 'Xs_I', \
     'POST_Xs_I.dat'  u 1:3 w lp lw 2 dt 2 t 'Xs_I DEM', \
     'POST_Xs_I.dat'  u 1:4 w p pt 7 lc  'red' ps 1 axis x1y2 t 'abs err'

unset multiplot
reset

脚本2:名字gnuplot _pmass

reset
set xlabel 'time'
set ylabel 'pmass'
set y2label 'relative err'
set y2tics 
set ytics nomirror
set y2range [0 : 5e-6]
set yrange [0.03 : 0.036]
#set key bottom
plot 'POST_MASS.dat' skip 9 u 1:2 w lp lw 2 t 'pmass', \
'POST_MASS.dat' skip 9 u 1:3 w lp lw 2 dt 2 t 'pmass DEM', \
'POST_MASS.dat' skip 9 u 1:4 w p pt 7 lc  'red' ps 2 axis x1y2 t 'abs err'

结果

1
在这里插入图片描述
2
在这里插入图片描述

数据

POST_MASS.dat

   Initial Conditions:
               Analytic       MFIX-DEM
      Mass     0.03534292     0.03534290
      MW_Bs:  15.00000000
      MW_Ds:  25.00000000
      Reaction Rate:   0.00039270


       S_TIME          Mp           Mp-DEM        ABS ERR
     0.00000836     0.03534290     0.03534290     0.00000000
     0.10000706     0.03514655     0.03514636     0.00000020
     0.20000706     0.03495020     0.03495001     0.00000020
     0.30000706     0.03475385     0.03475366     0.00000020
     0.40000706     0.03455750     0.03455731     0.00000020
     0.50000706     0.03436115     0.03436096     0.00000020
     0.60000706     0.03416480     0.03416461     0.00000020
     0.70000706     0.03396845     0.03396826     0.00000020
     0.80000706     0.03377210     0.03377191     0.00000020
     0.90000706     0.03357575     0.03357556     0.00000020
     1.00000706     0.03337940     0.03337921     0.00000020
     1.10000706     0.03318305     0.03318286     0.00000020
     1.20000706     0.03298670     0.03298651     0.00000020
     1.30000706     0.03279035     0.03279016     0.00000020
     1.40000706     0.03259400     0.03259381     0.00000020
     1.50000706     0.03239767     0.03239778     0.00000010
     1.60000706     0.03239767     0.03239778     0.00000010
     1.70000706     0.03239767     0.03239778     0.00000010
     1.80000706     0.03239767     0.03239778     0.00000010
     1.90000706     0.03239767     0.03239778     0.00000010
     2.00000706     0.03239767     0.03239778     0.00000010

POST_Xs_B.dat

       S_TIME          XB           XB-DEM        ABS ERR
     0.00000836     0.49999745     0.49999745     0.00000000
     0.10000706     0.46927149     0.46924059     0.00003090
     0.20000706     0.43819990     0.43816865     0.00003125
     0.30000706     0.40677721     0.40674561     0.00003160
     0.40000706     0.37499744     0.37496548     0.00003196
     0.50000706     0.34285448     0.34282216     0.00003233
     0.60000706     0.31034206     0.31030936     0.00003270
     0.70000706     0.27745377     0.27742069     0.00003308
     0.80000706     0.24418306     0.24414959     0.00003346
     0.90000706     0.21052321     0.21048936     0.00003386
     1.00000706     0.17646737     0.17643311     0.00003426
     1.10000706     0.14200850     0.14197383     0.00003466
     1.20000706     0.10713940     0.10710432     0.00003508
     1.30000706     0.07185270     0.07181720     0.00003550
     1.40000706     0.03614087     0.03610494     0.00003593
     1.50000706     0.00000000     0.00000000     0.00000000
     1.60000706     0.00000000     0.00000000     0.00000000
     1.70000706     0.00000000     0.00000000     0.00000000
     1.80000706     0.00000000     0.00000000     0.00000000
     1.90000706     0.00000000     0.00000000     0.00000000
     2.00000706     0.00000000     0.00000000     0.00000000

POST_Xs_D.dat

   S_TIME          XD           XD-DEM        ABS ERR
 0.00000836     0.00000232     0.00000232     0.00000000
 0.10000706     0.02793501     0.02796310     0.00002809
 0.20000706     0.05618191     0.05621032     0.00002841
 0.30000706     0.08474799     0.08477672     0.00002873
 0.40000706     0.11363869     0.11366774     0.00002906
 0.50000706     0.14285956     0.14288895     0.00002939
 0.60000706     0.17241631     0.17244604     0.00002973
 0.70000706     0.20231475     0.20234483     0.00003007
 0.80000706     0.23256086     0.23259128     0.00003042
 0.90000706     0.26316072     0.26319149     0.00003078
 1.00000706     0.29412057     0.29415172     0.00003114
 1.10000706     0.32544682     0.32547833     0.00003151
 1.20000706     0.35714600     0.35717789     0.00003189
 1.30000706     0.38922482     0.38925709     0.00003227
 1.40000706     0.42169012     0.42172278     0.00003266
 1.50000706     0.45454545     0.45454715     0.00000170
 1.60000706     0.45454545     0.45454715     0.00000170
 1.70000706     0.45454545     0.45454715     0.00000170
 1.80000706     0.45454545     0.45454715     0.00000170
 1.90000706     0.45454545     0.45454715     0.00000170
 2.00000706     0.45454545     0.45454715     0.00000170

POST_Xs_I.dat

   S_TIME          XI           XI-DEM        ABS ERR
 0.00000836     0.50000023     0.50000023     0.00000000
 0.10000706     0.50279350     0.50279631     0.00000281
 0.20000706     0.50561819     0.50562103     0.00000284
 0.30000706     0.50847480     0.50847767     0.00000287
 0.40000706     0.51136387     0.51136677     0.00000291
 0.50000706     0.51428596     0.51428889     0.00000294
 0.60000706     0.51724163     0.51724460     0.00000297
 0.70000706     0.52023148     0.52023448     0.00000301
 0.80000706     0.52325609     0.52325913     0.00000304
 0.90000706     0.52631607     0.52631915     0.00000308
 1.00000706     0.52941206     0.52941517     0.00000311
 1.10000706     0.53254468     0.53254783     0.00000315
 1.20000706     0.53571460     0.53571779     0.00000319
 1.30000706     0.53892248     0.53892571     0.00000323
 1.40000706     0.54216901     0.54217228     0.00000327
 1.50000706     0.54545455     0.54545285     0.00000170
 1.60000706     0.54545455     0.54545285     0.00000170
 1.70000706     0.54545455     0.54545285     0.00000170
 1.80000706     0.54545455     0.54545285     0.00000170
 1.90000706     0.54545455     0.54545285     0.00000170
 2.00000706     0.54545455     0.54545285     0.00000170
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值