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

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

off999 2024-10-26 11:57 47 浏览 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)。

相关推荐

安全教育登录入口平台(安全教育登录入口平台官网)

122交通安全教育怎么登录:122交通网的注册方法是首先登录网址http://www.122.cn/,接着打开网页后,点击右上角的“个人登录”;其次进入邮箱注册,然后进入到注册页面,输入相关信息即可完...

大鱼吃小鱼经典版(大鱼吃小鱼经典版(经典版)官方版)

大鱼吃小鱼小鱼吃虾是于谦跟郭麒麟的《我的棒儿呢?》郭德纲说于思洋郭麒麟作诗的相声,最后郭麒麟做了一首,师傅躺在师母身上大鱼吃小鱼小鱼吃虾虾吃水水落石出师傅压师娘师娘压床床压地地动山摇。...

谷歌地球下载高清卫星地图(谷歌地球地图下载器)
  • 谷歌地球下载高清卫星地图(谷歌地球地图下载器)
  • 谷歌地球下载高清卫星地图(谷歌地球地图下载器)
  • 谷歌地球下载高清卫星地图(谷歌地球地图下载器)
  • 谷歌地球下载高清卫星地图(谷歌地球地图下载器)
哪个软件可以免费pdf转ppt(免费的pdf转ppt软件哪个好)
哪个软件可以免费pdf转ppt(免费的pdf转ppt软件哪个好)

要想将ppt免费转换为pdf的话,我们建议大家可以下一个那个wps,如果你是会员的话,可以注册为会员,这样的话,在wps里面的话,就可以免费将ppt呢转换为pdfpdf之后呢,我们就可以直接使用,不需要去直接不需要去另外保存,为什么格式转...

2026-02-04 09:03 off999

电信宽带测速官网入口(电信宽带测速官网入口app)

这个网站看看http://www.swok.cn/pcindex.jsp1.登录中国电信网上营业厅,宽带光纤,贴心服务,宽带测速2.下载第三方软件,如360等。进行在线测速进行宽带测速时,尽...

植物大战僵尸95版手机下载(植物大战僵尸95 版下载)

1可以在应用商店或者游戏平台上下载植物大战僵尸95版手机游戏。2下载教程:打开应用商店或者游戏平台,搜索“植物大战僵尸95版”,找到游戏后点击下载按钮,等待下载完成即可安装并开始游戏。3注意:确...

免费下载ppt成品的网站(ppt成品免费下载的网站有哪些)

1、Chuangkit(chuangkit.com)直达地址:chuangkit.com2、Woodo幻灯片(woodo.cn)直达链接:woodo.cn3、OfficePlus(officeplu...

2025世界杯赛程表(2025世界杯在哪个国家)

2022年卡塔尔世界杯赛程公布,全部比赛在卡塔尔境内8座球场举行,2022年,决赛阶段球队全部确定。揭幕战于当地时间11月20日19时进行,由东道主卡塔尔对阵厄瓜多尔,决赛于当地时间12月18日...

下载搜狐视频电视剧(搜狐电视剧下载安装)

搜狐视频APP下载好的视频想要导出到手机相册里方法如下1、打开手机搜狐视频软件,进入搜狐视频后我们点击右上角的“查找”,找到自已喜欢的视频。2、在“浏览器页面搜索”窗口中,输入要下载的视频的名称,然后...

pubg免费下载入口(pubg下载入口官方正版)
  • pubg免费下载入口(pubg下载入口官方正版)
  • pubg免费下载入口(pubg下载入口官方正版)
  • pubg免费下载入口(pubg下载入口官方正版)
  • pubg免费下载入口(pubg下载入口官方正版)
永久免费听歌网站(丫丫音乐网)

可以到《我爱音乐网》《好听音乐网》《一听音乐网》《YYMP3音乐网》还可以到《九天音乐网》永久免费听歌软件有酷狗音乐和天猫精灵,以前要跳舞经常要下载舞曲,我从QQ上找不到舞曲下载就从酷狗音乐上找,大多...

音乐格式转换mp3软件(音乐格式转换器免费版)

有两种方法:方法一在手机上操作:1、进入手机中的文件管理。2、在其中选择“音乐”,将显示出手机中的全部音乐。3、点击“全选”,选中所有音乐文件。4、点击屏幕右下方的省略号图标,在弹出菜单中选择“...

电子书txt下载(免费的最全的小说阅读器)

1.Z-library里面收录了近千万本电子书籍,需求量大。2.苦瓜书盘没有广告,不需要账号注册,使用起来非常简单,直接搜索预览下载即可。3.鸠摩搜书整体风格简洁清晰,书籍资源丰富。4.亚马逊图书书籍...

最好免费观看高清电影(播放免费的最好看的电影)

在目前的网上选择中,IMDb(互联网电影数据库)被认为是最全的电影网站之一。这个网站提供了各种类型的电影和电视节目的海量信息,包括剧情介绍、演员表、评价、评论等。其还提供了有关电影制作背后的详细信息,...

孤单枪手2简体中文版(孤单枪手2简体中文版官方下载)

要将《孤胆枪手2》游戏的征兵秘籍切换为中文,您可以按照以下步骤进行操作:首先,打开游戏设置选项,通常可以在游戏主菜单或游戏内部找到。然后,寻找语言选项或界面选项,点击进入。在语言选项中,选择中文作为游...

取消回复欢迎 发表评论: