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

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 )

4585658b9beecce61c88ce6dd5a33b92.png

这里的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使用教程详解

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

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

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

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

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

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

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

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

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/542616.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

dda算法_计算机图形学中的DDA(数字差分分析仪)算法

dda算法DDA(数字差分分析仪)算法 (DDA (Digital Differential Analyzer) Algorithm) In computer graphics, the DDA algorithm is the simplest algorithm among all other line generation algorithms. Here, the DDA is an abbreviation that stands for "Digital Diff…

购物商城框架java_基于jsp的购物商城-JavaEE实现购物商城 - java项目源码

基于jspservletpojomysql实现一个javaee/javaweb的购物商城, 该项目可用各类java课程设计大作业中, 购物商城的系统架构分为前后台两部分, 最终实现在线上进行购物商城各项功能,实现了诸如用户管理, 登录注册, 权限管理等功能, 并实现对各类购物商城相关的实体进行管理。该购物…

c语言++数组名【数字】_C ++程序在数组中打印所有非重复数字

c语言数组名【数字】Problem statement: Write a C program to print all the non-repeated numbers in an array in minimum time complexity. 问题陈述:编写一个C 程序, 以最小的时间复杂度将所有未重复的数字打印在数组中 。 Input Example: 输入示例…

java最接近对点及距离_最接近点对问题_分治法

一、问题描述给定平面上的n个点,找其中的一对点,使得在n个点组成的所有点对中该点对间的距离最小。二、解题思路及所选算法策略的可行性分析思路:利用分治法来解决问题。递归子结构求最接近点对总体可分为几个步骤:1、当问题规模小…

python return用法_初学Python要了解什么 装饰器知识汇总有哪些

初学Python要了解什么?装饰器知识汇总有哪些?在Python学习过程中,有多种方法对函数和类进行加工,相对于其它方式,装饰器语法简单,代码可读性高。因此,装饰器在Python项目中有广泛的应用&#xf…

android emulator虚拟设备分析第三篇之pipe上的qemud service

一、概述 本篇和第二篇是强相关的&#xff0c;需要结合第二篇一起看。 以boot-properties为例&#xff0c;注意不需要看ANDROID-QEMUD.TXT&#xff0c;这个是和guest os中的qemud进行相关的&#xff0c;已废弃。 启动emulator时&#xff0c;有一个参数-prop <key><val…

c#异常处理_C#异常处理能力问题和解答 套装4

c#异常处理1) Which is not a valid keyword used in the context of exception handling? trycatchfinalfinally Answer & Explanation Correct answer: 3final The final keyword is not used to handle exceptions in C#.NET. 1)在异常处理的上下文中使用哪个无效关键字…

Castor xsd生成java_java – Castor可以处理从基础XSD导入的多个XSD生成类吗?

注意&#xff1a;我是EclipseLink JAXB (MOXy)领导者,也是JAXB 2 (JSR-222)专家组的成员.Can Castor do this? If so, what would be the Ant task syntax for it.If not, would perhaps JAXB be a better alternative?下面是如何使用JAXB完成此操作的示例&#xff1a;产品xm…

串口通信 校验码_一文读懂S7-200 SMART自由口通信!

学习S7-200 SMART时了解到&#xff0c;基于RS485接口可实现一下几种通信&#xff1a;1&#xff09;modbus RTU通信2&#xff09;PPI协议通信3&#xff09;USS协议通信4&#xff09;自由口通信何为自由口通信呢&#xff1f;前三种通信必须要PLC和与其通信的设备支持相同的通信协…

hbase 学习(十三)集群间备份原理

集群建备份&#xff0c;它是master/slaves结构式的备份&#xff0c;由master推送&#xff0c;这样更容易跟踪现在备份到哪里了&#xff0c;况且region server是都有自己的WAL 和HLog日志&#xff0c;它就像mysql的主从备份结构一样&#xff0c;只有一个日志来跟踪。一个master集…

python expect模块_Python基础教程:用Python怎么telnet到网络设备

Python基础教程&#xff1a;用Python怎么telnet到网络设备0.前言Telnet协议属于TCP/IP协议族里的一种&#xff0c;对于我们这些网络攻城狮来说&#xff0c;再熟悉不过了&#xff0c;常用于远程登陆到网络设备进行操作&#xff0c;但是&#xff0c;它的缺陷太明显了&#xff0c;…

Java实现动态加载页面_[Java教程]动态加载页面数据的小工具 javascript + jQuery (持续更新)...

[Java教程]动态加载页面数据的小工具 javascript jQuery (持续更新)0 2014-05-07 18:00:06使用该控件&#xff0c;可以根据url&#xff0c;参数&#xff0c;加载html记录模板(包含json参数对应&#xff0c;以及具体记录位置Index根据参数描述加载对应的属性&#xff0c;并可以…

马哥linux第六周作业

1、复制/etc/rc.d/rc.sysinit文件至/tmp目录&#xff0c;将/tmp/rc.sysinit文件中的以至少一个空白字符开头的行的行首加#&#xff1b;[rootmageedu tmp]# cp /etc/rc.d/rc.sysinit . [rootmageedu tmp]# vim rc.sysinit :% s/^[[:space:]]/#&/ #按Esc进入vi…

Java ObjectInputStream enableResolveObject()方法与示例

ObjectInputStream类enableResolveObject()方法 (ObjectInputStream Class enableResolveObject() method) enableResolveObject() method is available in java.io package. enableResolveObject()方法在java.io包中可用。 enableResolveObject() method is used to enable th…

pygame render怎么显示中文_PyGame开发游戏(2D)02.基础图元

这节将介绍PyGame的基础架构。并学习如何在PyGame里绘制各种几何图形和显示加载图片。01.应用框架上一节的示例程序里&#xff0c;我们用到一个PyGame的应用程序框架。这是一个基础框架&#xff0c;利用它我们可以很轻松的添加各类图型绘制&#xff0c;键盘鼠标输入处理和各类逻…

word+增加水印+java_为Word2019文档添加水印的两种方法

水印的类型包括文字水印和图片水印两种。在Word文档中添加文字水印时&#xff0c;可以使用程序中预设的水印效果&#xff0c;而图片水印则需要自定义添加。一、使用程序预设的文字水印Word 2019中预设了机密、紧急、免责声明三种类型的文字水印&#xff0c;用户可根据文件的类型…

如何设置CentOS 7获取动态及静态IP地址

自动获取动态IP地址1.输入“ip addr”并按回车键确定&#xff0c;发现无法获取IP(CentOS 7默认没有ifconfig命令)&#xff0c;记录下网卡名称&#xff08;本例中为ens33&#xff09;。2.输入“cd /etc/sysconfig/network-scripts/”按回车键确定&#xff0c;继续输入“ls”按回…

请求列出指定服务器上的可用功能失败_滥用 ESI 详解(上)

在进行安全性评估时&#xff0c;我们注意到了标记语言 Edge Side Includes (ESI)中的一个意外行为&#xff0c;这种语言用于许多流行的 HTTP 代理(反向代理、负载平衡器、缓存服务器、代理服务器)。我们发现成功的 ESI 攻击可以导致服务器端请求伪造(SSRF)、各种绕过 HTTPOnly …

Java ClassLoader setPackageAssertionStatus()方法与示例

ClassLoader类setPackageAssertionStatus()方法 (ClassLoader Class setPackageAssertionStatus() method) setPackageAssertionStatus() method is available in java.lang package. setPackageAssertionStatus()方法在java.lang包中可用。 setPackageAssertionStatus() metho…

java上传kafka的方法_哪种方法是将所有数据从Kafka主题复制到接收器(文件或Hive表)的最佳方法?...

我正在使用Kafka Consumer API将所有数据从Kafka主题复制到Hive表 . 为此&#xff0c;我使用HDFS作为中间步骤 . 我使用唯一的组ID并将偏移重置为“最早”&#xff0c;以便从头开始获取所有数据&#xff0c;并在执行后忽略提交 . 然后我遍历Kafka主题中的记录&#xff0c;并将每…