百度360必应搜狗淘宝本站头条
当前位置:网站首页 > 技术资源 > 正文

Python数学编程 Python第四章 第四节 Sympy包绘图

off999 2024-10-26 11:57 35 浏览 0 评论

在第二章中我们学习了用已知的数据绘制图形。例如,为了绘制两物体间引力与物体之间距离的关系图,你必须计算每一个距离值对应的引力,并将包含引力和距离的列表传递给matplotlib。与此不同的是,Sympy可以直接绘制你创建的方程的图形。下面看看如何绘制方程y = 2x + 3的图形:


 >>> from sympy.plotting import plot
 >>> from sympy import Symbol
 >>> x = Symbol('x')
 >>> plot(2*x+3)

我们只需要导入plot和Symbol,创建符号对象x并且将函数表达式传入plot()函数,Sympy就会绘制出这和函数的图像,如图所示:

图形显示的x值是默认自动选择的范围:-10.0到10.0。这里的窗口和第二章第三章的窗口相似,这是因为Sympy在幕后自动调用matplotlib绘制图形。注意,我们没有调用show()函数来显示图形,Sympy会自动完成。

假如你希望设置x的区间范围在-5到5的范围内,方法如下:


 >>> plot((2*x + 3), (x, -5, 5))

在这里,包含符号、范围的下限和上限的元组(x, -5, 5)作为plot()函数的第二个参数。现在图形就只显示x的范围为-5到5的图像。

也可以在plot()函数中使用其他关键字,如title输入标题,或者分别以xlabel和ylabel标记x轴和y轴。以下的plot()函数设定了这三个关键字参数:


 >>> plot((2*x + 3), (x, -5, 5), title='A Line', xlabel='x',ylabel='2x+3')


上图显示了标题、x轴以及y轴上的标签。还可以通过plot()函数的其他关键字指定其绘制图形的细节,show参数可以控制是否显示图形,调用plot()函数时设置show=False将不显示图形:


 >>> p = plot((2*x + 3), (x, -5, 5), title='A Line', xlabel='x',ylabel='2x+3', show=False)

我们使用标签p来指代创建的图形对象,执行此条语句,我们会发现图形并没有显示。我们可以通过调用p.show()函数来显示图形,还可以使用save()函数将图形保存为一个图像文件:


 >>> p.show()
 >>> p.save('C://line.png')

这样就可以将图形对象保存到指定路径。

4.4.1 绘制用户输入的表达式

传递给用户的函数表达式只包含x,例如前面绘制的y = 2x + 3,输入的绘制函数就只能是2x + 3,如果表达式不是这种形式就要进行重写,当然在程序执行之前我们可以手动完成。但是假如要编写一个程序,允许用户绘制任何表达式的图形,那该怎么办?假设用户输入类似2x + 3y - 6的表达式,我们需要先转换它,solve()函数能帮助我们实现转换。让我们来看一个例子:


 >>> expr = input('Enter an expression : ')
 Enter an Expression : 2*x + 3*y - 6
 >>> expr = sympify(expr) # 将输入的字符串转为表达式
 >>> y = Symbol('y')
 >>> solve(expr, y)
 [2 - 2*x/3]

我们调用sympify()函数将输入的表达式转换为一个Sympy对象,然后创建一个Symbol对象y,这样就能告诉Sympy我们要为哪个变量解方程,回顾之前解方程使用的solve()函数,传入的表达式等于0,然后解出指定的符号对应的表达式。可以看到我们解得 y = 2 - 2*x/3,我们只需要提取保存即可,因为它是以列表形式返回的。


 >>> solutions = solve(expr, y)
 >>> expr_y = solutions[0]
 >>> expr_y
 2 - 2*x/3

因为solve()函数返回的是只包含一项的y对应的表达式,因此我们只需要提取第一个元素即可。现在我们就可以调用plot()函数来绘制图形。下述代码显示了一个完整的图形绘制程序:


 '''
 Plot the graph of an input expression
 '''
 
 from sympy import Symbol, sympify, solve, SympifyError
 from sympy.plotting import plot
 
 def plot_expression(expr):
     y = Symbol('y')
     solutions = solve(expr, y)
     expr_y = solutions[0]
     plot(expr_y, xlabel='x', ylabel=expr_y, title='B Line')
 if __name__ == '__main__':
     expr = input('Enter your expression : ')
     try:
         expr = sympify(expr)
     except SympifyError:
         print('Invalid input')
     else:
         plot_expression(expr)

搭配try...except模块,检查用户的无效输入,与之前在sympify()函数完成的类似。运行这个程序,它会让你输入一个表达式,然后将相应的图形绘制出来。

4.4.2 多函数图像绘制

你可以输入多个表达式,并调用Sympy的plot()函数在同一张图上绘制它们。例如,下面的代码可以实现同时绘制两条线,如图所示:


 >>> from sympy.plotting import plot
 >>> from sympy import Symbol
 >>> x = Symbol('x')
 >>> plot(2*x+3, 3*x+1)

我们可以为这个图形添加图例来区别这两条线


 >>> from sympy.plotting import plot
 >>> from sympy import Symbol
 >>> x = Symbol('x')
 >>> p = plot(2*x+3, 3*x+1, legend=True, show=False)
 >>> p[0].line_color = 'b'
 >>> p[1].line_color = 'r' 
 >>> p.show()

可以看到,我们为图像添加的图例与所绘制的表达式相匹配,并且,我们还为每条线设置了线条颜色,'b' 为blue,'r'为red,'g'为green,'c'为cyan(青色),'m'为magenta(紫红色),'y'为yellow,'w'为white(使用每种颜色的第一个字母表示,如果颜色首字母相同,可以输入颜色整个词以进行区分,如black和blue)。

相关推荐

戴尔官网保修查询入口(戴尔售后保质期查询)

可以按照以下步骤查询戴尔笔记本电脑的保修期:1.打开戴尔官网:https://www.戴尔.com/zh-cn/售后服务/保修政策.html2.点击页面上方的“服务与支持”按钮,进入戴尔的服务支持...

手机号邮箱登录入口(手机号邮箱官网)

手机163邮箱登录入口如下:163邮箱官网入口:https://smart.mail.163.com/login.htm点击进入登录或者注册邮箱即可。手机浏览器访问进入官网http://www.123...

sd卡(sd卡无法读取怎么修复)

  SD卡是大卡,相机用的;普通的手机内存卡,是小卡,正规的名称是macrosd卡,也就是微型SD卡。可以通过卡套转为普通的SD卡的大小。  其实就是大小不同。但手机上的内存卡,人们经常也俗称为SD...

路由器连接图(网络路由器连接图)
  • 路由器连接图(网络路由器连接图)
  • 路由器连接图(网络路由器连接图)
  • 路由器连接图(网络路由器连接图)
  • 路由器连接图(网络路由器连接图)
windows7蓝牙功能在哪里打开

点击搜索框在windows7系统主界面点击开始菜单,点击打开搜索框。输入命令输入services.msc后回车,在列表中找到并右击BluetoothSupportS...点击属性选择进入属性菜单,...

2010激活密钥(microsoft2010激活密钥)
2010激活密钥(microsoft2010激活密钥)

步骤/方式1officeprofessionalplus2010:(office专业版)6QFdx-pYH2G-ppYFd-C7RJM-BBKQ8Bdd3G-xM7FB-Bd2HM-YK63V-VQFdKVYBBJ-TRJpB-QFQ...

2025-11-19 04:03 off999

联想官方刷新bios工具(联想电脑刷新bios)

刷新BIOS需要使用联想的官方网站或授权维修中心来进行操作。以下是一些基本步骤:1.访问联想的官方网站,找到BIOS更新程序并下载。在下载过程中,请确保选择与您计算机型号匹配的版本。2.将下载的B...

苹果ios14系统下载(苹果ios14.1下载)
苹果ios14系统下载(苹果ios14.1下载)

1方法一步骤/方式一打开Appstore。步骤/方式二在搜索栏点击搜索框。步骤/方式三搜索并点击需要下载的软件。步骤/方式四点击获取。步骤/方式五最后验证ID密码即可。1.在应用商店搜索你要下载的应用名称。2.点击下载按钮,如果要求登...

2025-11-19 03:03 off999

office2010怎么免费永久激活密钥

用这个试试,一个KMS激活工具可以激活2010到2019的Office自家的目前用的就是这个microsoft6477.moe/1716.html直接使用这个Microsoftoffice2010...

类似爱加速的国内ip(类似爱加速的app)
类似爱加速的国内ip(类似爱加速的app)

推荐“V8盒子”。这一款免费无广告的模拟器,不同于其它软件盒子,而是类似于X8沙箱,满足游戏多开,画中画,悬浮球操作,熄屏后台运行等多功能的沙箱盒子.支持一键root,一键安装xposed框架,能在安卓/苹果手机上运行多个安卓/ios虚拟系...

2025-11-19 02:03 off999

阿里旺旺手机客户端(阿里旺旺手机app)

手机淘宝的旺旺在打开商品后,会看到左下角有个旺旺的图标,点击就可以联系了。  阿里旺旺是将原先的淘宝旺旺与阿里巴巴贸易通整合在一起的一个新品牌。它是淘宝和阿里巴巴为商人量身定做的免费网上商务沟通软件,...

最纯净的pe装机工具(pe工具哪个纯净)

U盘装系统步骤:1.制作U盘启动盘。这里推荐大白菜U盘启动盘制作工具,在网上一搜便是。2.U盘启动盘做好了,我们还需要一个GHOST文件,可以从网上下载一个ghost版的XP/WIN7/WIN8系统,...

装一个erp系统多少钱(wms仓库管理软件)

现在主流有客户端ERP和云端ERP两种客户端通常一次买断,价格在万元左右,但是还有隐性费用,你需要支付服务器、数据管理员,此外如果系统需要更新维护,你还需要支付另外一笔不菲的费用。云端ERP:优势...

cad2014序列号和密钥永久(autocad2014序列号和密钥)

1在cad2014中修改标注样式后,需要将其保存2单击“样式管理器”按钮,在弹出的窗口中选择修改后的标注样式,然后单击“设置为当前”按钮,再单击“保存当前样式”按钮,将其保存为新的样式名称3为了...

qq修改密保手机号(qq修改密保手机号是什么意思)

QQ更改绑定的手机号码操作步骤如下:1、打开手机主界面,找到“QQ”软件点击打开。2、输入正确的QQ账户和密码登录到qq主界面。3、点击左上角的头像“图片”,进入到个人中心界面。4、进入到个人中心界面...

取消回复欢迎 发表评论: