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

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

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

相关推荐

面试官:来,讲一下枚举类型在开发时中实际应用场景!

一.基本介绍枚举是JDK1.5新增的数据类型,使用枚举我们可以很好的描述一些特定的业务场景,比如一年中的春、夏、秋、冬,还有每周的周一到周天,还有各种颜色,以及可以用它来描述一些状态信息,比如错...

一日一技:11个基本Python技巧和窍门

1.两个数字的交换.x,y=10,20print(x,y)x,y=y,xprint(x,y)输出:102020102.Python字符串取反a="Ge...

Python Enum 技巧,让代码更简洁、更安全、更易维护

如果你是一名Python开发人员,你很可能使用过enum.Enum来创建可读性和可维护性代码。今天发现一个强大的技巧,可以让Enum的境界更进一层,这个技巧不仅能提高可读性,还能以最小的代价增...

Python元组编程指导教程(python元组的概念)

1.元组基础概念1.1什么是元组元组(Tuple)是Python中一种不可变的序列类型,用于存储多个有序的元素。元组与列表(list)类似,但元组一旦创建就不能修改(不可变),这使得元组在某些场景...

你可能不知道的实用 Python 功能(python有哪些用)

1.超越文件处理的内容管理器大多数开发人员都熟悉使用with语句进行文件操作:withopen('file.txt','r')asfile:co...

Python 2至3.13新特性总结(python 3.10新特性)

以下是Python2到Python3.13的主要新特性总结,按版本分类整理:Python2到Python3的重大变化Python3是一个不向后兼容的版本,主要改进包括:pri...

Python中for循环访问索引值的方法

技术背景在Python编程中,我们经常需要在循环中访问元素的索引值。例如,在处理列表、元组等可迭代对象时,除了要获取元素本身,还需要知道元素的位置。Python提供了多种方式来实现这一需求,下面将详细...

Python enumerate核心应用解析:索引遍历的高效实践方案

喜欢的条友记得关注、点赞、转发、收藏,你们的支持就是我最大的动力源泉。根据GitHub代码分析统计,使用enumerate替代range(len())写法可减少38%的索引错误概率。本文通过12个生产...

Python入门到脱坑经典案例—列表去重

列表去重是Python编程中常见的操作,下面我将介绍多种实现列表去重的方法,从基础到进阶,帮助初学者全面掌握这一技能。方法一:使用集合(set)去重(最简单)pythondefremove_dupl...

Python枚举类工程实践:常量管理的标准化解决方案

本文通过7个生产案例,系统解析枚举类在工程实践中的应用,覆盖状态管理、配置选项、错误代码等场景,适用于Web服务开发、自动化测试及系统集成领域。一、基础概念与语法演进1.1传统常量与枚举类对比#传...

让Python枚举更强大!教你玩转Enum扩展

为什么你需要关注Enum?在日常开发中,你是否经常遇到这样的代码?ifstatus==1:print("开始处理")elifstatus==2:pri...

Python枚举(Enum)技巧,你值得了解

枚举(Enum)提供了更清晰、结构化的方式来定义常量。通过为枚举添加行为、自动分配值和存储额外数据,可以提升代码的可读性、可维护性,并与数据库结合使用时,使用字符串代替数字能简化调试和查询。Pytho...

78行Python代码帮你复现微信撤回消息!

来源:悟空智能科技本文约700字,建议阅读5分钟。本文基于python的微信开源库itchat,教你如何收集私聊撤回的信息。[导读]Python曾经对我说:"时日不多,赶紧用Python"。于是看...

登录人人都是产品经理即可获得以下权益

文章介绍如何利用Cursor自动开发Playwright网页自动化脚本,实现从选题、写文、生图的全流程自动化,并将其打包成API供工作流调用,提高工作效率。虽然我前面文章介绍了很多AI工作流,但它们...

Python常用小知识-第二弹(python常用方法总结)

一、Python中使用JsonPath提取字典中的值JsonPath是解析Json字符串用的,如果有一个多层嵌套的复杂字典,想要根据key和下标来批量提取value,这是比较困难的,使用jsonpat...

取消回复欢迎 发表评论: