matlab数值计算pdf_Gnuplot科学绘图（九）——栅格以及方程数值解估算

最新推荐文章于 2021-05-21 14:40:47 发布

weixin_39674414

最新推荐文章于 2021-05-21 14:40:47 发布

阅读量257

点赞数

文章标签： matlab数值计算pdf

Gnuplot科学绘图系列内容

Gnuplot科学绘图(一)——从安装到简单函数绘图(文末有彩蛋)

Gnuplot科学绘图(二)——坐标取值范围及刻度(文末有彩蛋)

Gnuplot科学绘图(三)——点线风格

Gnuplot科学绘图(四)——多组数据绘图

Gnuplot科学绘图(五)——输出eps 图片

Gnuplot科学绘图(六)——输出pdf 和png 图片

Gnuplot科学绘图(七)——Enhanced 文本模式

Gnuplot科学绘图(八)——插入LATEX 公式

【作者前言】

对于科技工作者来说，gnuplot 是一个非常好用的绘图软件。因为感慨于gnuplot 中文资料和文档的缺乏，我于数月之前在科学网博客开始撰写“谈谈gnuplot”系列博文，至今已写了四十五篇。

虽然还有不少细节问题没有涉及，但是我觉得这些博文已经可以构成一个比较完整的gnuplot 中文教程。对大多数用户来说，这一系列博文已经涵盖了日常应用的需要。为了方便大家阅读，我现在将这一系列博文集结编辑成为此教程，以CC BY-NC-SA 知识共享协议发布。

我们现在来画一个0 阶贝塞尔函数J0(x)：

gnuplot> set term wxt enhanced
gnuplot> set xlabel “X”
gnuplot> set ylabel “Y”
gnuplot> set xrange [ 0 : 1 0 ]
gnuplot> set xtics 0 ,1 ,10
gnuplot> unset key
gnuplot> set title “0阶贝塞尔函数J_0( x ) ”
gnuplot> plot besj0 ( x )

这里的besj0(x) 就是gnuplot 里面预定义的0 阶贝塞尔函数。如果现在请您从这个图上估计出[0, 10] 内J0(x) 的零点数值，也就是方程J0(x) = 0 的解，恐怕您很难说的准确。但是如果为这个图加上栅格(grid)，就容易多了：

gnuplot> set grid
gnuplot> replot

这时我们很容易估计出三个零点的数值：2.4, 5.5, 8.6。通过查表我们可以知道，这三个零点比较精确的数值分别为2.4048, 5.5201, 8.6537。这和我们的估计值差不太多。如果我们想更精确的估计数值，可以尝试改一下xrange：

gnuplot> set xrange [ 8 : 9 ]
gnuplot> set xtics 8 , 0 . 1 , 9
gnuplot> replot

这相当于把图像在零点附近放大了。把鼠标放在画图区域，画图框左下角就会显示出鼠标所在位置的标。现在我们把鼠标放在函数图线和X轴的交叉点上，左下角显示的横坐标为8.65243，这和我们查表所得的数值更接近了。

如果想进一步让结果精确一些，我们可以利用gnuplot 的计算功能。我们可以通过尝试计算的方法获得方程的数值解：

gnuplot> print besj0 ( 8 . 6 5 )
0.00101216621937318
gnuplot> print besj0 ( 8 . 6 6 )
-0.0017019446057587
gnuplot> print besj0 ( 8 . 6 5 3 7 )
7.5770361108123 e-06
gnuplot> print besj0 ( 8 . 6 5 3 6 )
3.47225104115535 e-05
gnuplot> print besj0 ( 8 . 6 5 3 8 )
-1.95681245811775 e-05

所以在8.6 附近，J0(x) = 0 精确到小数点后4 位的数值解为8.6537，这和我们查表的结果一模一样。由于我们已经通过图像知道了数值解的大概位置，再加上合理利用线性插值，我们可以很快得到精确的结果。

本文整理自马欢老师科学网博客，特此感谢原作者的分享。

Materials Studio系列内容中科大 Materials Studio 培训快速入门教程(包你学会！)精选众多论坛讨论贴，Materials-Studio 论坛问答全集怎样把Materials studio能带图做得漂亮，用Matlab绘制呀！Materials Studio手把手教程丨界面操作、应用与NaCl晶体构建实例怎么利用Materials Studio建立晶体模型？就是这么简单！手把手教你用Materials Studio计算碳纳米管的能带结构基于Materials Studio的异质结构建模技巧如何分析第一原理计算结果：电荷密度图、能带结构、态密度的分析70张图，搞定能带理论基础！

【长篇图文，包教包会】Materials Studio软件——CASTEP使用教程详解