Python数学编程 Python第四章 第四节 Sympy包绘图
off999 2024-10-26 11:57 28 浏览 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)。
相关推荐
- Alist 玩家请进:一键部署全新分支 Openlist,看看香不香!
-
Openlist(其前身是鼎鼎大名的Alist)是一款功能强大的开源文件列表程序。它能像“万能钥匙”一样,解锁并聚合你散落在各处的云盘资源——无论是阿里云盘、百度网盘、GoogleDrive还是...
- 白嫖SSL证书还自动续签?这个开源工具让我告别手动部署
-
你还在手动部署SSL证书?你是不是也遇到过这些问题:每3个月续一次Let'sEncrypt证书,忘了就翻车;手动配置Nginx,重启服务,搞一次SSL得花一下午;付费证书太贵,...
- Docker Compose:让多容器应用一键起飞
-
CDockerCompose:让多容器应用一键起飞"曾经我也是一个手动启动容器的少年,直到我的膝盖中了一箭。"——某位忘记--link参数的运维工程师引言:容器化的烦恼与...
- 申请免费的SSL证书,到期一键续签
-
大家好,我是小悟。最近帮朋友配置网站HTTPS时发现,还有人对宝塔面板的SSL证书功能还不太熟悉。其实宝塔早就内置了免费的Let'sEncrypt证书申请和一键续签功能,操作简单到连新手都能...
- 飞牛NAS部署TVGate Docker项目,实现内网一键转发、代理、jx
-
前面分享了两期TVGate:Q大的转发代理工具TVGate升级了,操作更便捷,增加了新的功能跨平台内网转发神器TVGate部署与使用初体验现在项目已经开源,并支持Docker部署,本文介绍如何通...
- Docker Compose 编排实战:一键部署多容器应用!
-
当项目变得越来越复杂,一个服务已经无法满足需求时,你可能需要同时部署数据库、后端服务、前端网页、缓存组件……这时,如果还一个一个手动dockerrun,简直是灾难这就是DockerCompo...
- 深度测评:Vue、React 一键部署的神器 PinMe
-
不知道大家有没有这种崩溃瞬间:领导突然要看项目Demo,客户临时要体验新功能,自己写的小案例想发朋友圈;找运维?排期?还要走工单;自己买服务器?域名、SSL、Nginx、防火墙;本地起服务?断电、关...
- 超简单!一键启动多容器,解锁 Docker Compose 极速编排秘籍
-
想要用最简单的方式在本地复刻一套完整的微服务环境?只需一个docker-compose.yml文件,你就能一键拉起N个容器,自动组网、挂载存储、环境隔离,全程无痛!下面这份终极指南,教你如何用...
- 日志文件转运工具Filebeat笔记_日志转发工具
-
一、概述与简介Filebeat是一个日志文件转运工具,在服务器上以轻量级代理的形式安装客户端后,Filebeat会监控日志目录或者指定的日志文件,追踪读取这些文件(追踪文件的变化,不停的读),并将来自...
- K8s 日志高效查看神器,提升运维效率10倍!
-
通常情况下,在部署了K8S服务之后,为了更好地监控服务的运行情况,都会接入对应的日志系统来进行检测和分析,比如常见的Filebeat+ElasticSearch+Kibana这一套组合...
- 如何给网站添加 https_如何给网站添加证书
-
一、简介相信大家都知道https是更加安全的,特别是一些网站,有https的网站更能够让用户信任访问接下来以我的个人网站五岁小孩为例子,带大家一起从0到1配置网站https本次配置的...
- 10个Linux文件内容查看命令的实用示例
-
Linux文件内容查看命令30个实用示例详细介绍了10个Linux文件内容查看命令的30个实用示例,涵盖了从基本文本查看、分页浏览到二进制文件分析的各个方面。掌握这些命令帮助您:高效查看各种文本文件内...
- 第13章 工程化实践_第13章 工程化实践课
-
13.1ESLint+Prettier代码规范统一代码风格配置//.eslintrc.jsmodule.exports={root:true,env:{node...
- 龙建股份:工程项目中标_龙建股份有限公司招聘网
-
404NotFoundnginx/1.6.1【公告简述】2016年9月8日公告,公司于2016年9月6日收到苏丹共和国(简称“北苏丹”)喀土穆州基础设施与运输部公路、桥梁和排水公司出具的中标通知书...
- 福田汽车:获得政府补助_福田 补贴
-
404NotFoundnginx/1.6.1【公告简述】2016年9月1日公告,自2016年8月17日至今,公司共收到产业发展补助、支持资金等与收益相关的政府补助4笔,共计5429.08万元(不含...
你 发表评论:
欢迎- 一周热门
- 最近发表
- 标签列表
-
- python计时 (73)
- python安装路径 (56)
- python类型转换 (93)
- python进度条 (67)
- python吧 (67)
- python的for循环 (65)
- python格式化字符串 (61)
- python静态方法 (57)
- python列表切片 (59)
- python面向对象编程 (60)
- python 代码加密 (65)
- python串口编程 (77)
- python封装 (57)
- python写入txt (66)
- python读取文件夹下所有文件 (59)
- python操作mysql数据库 (66)
- python获取列表的长度 (64)
- python接口 (63)
- python调用函数 (57)
- python多态 (60)
- python匿名函数 (59)
- python打印九九乘法表 (65)
- python赋值 (62)
- python异常 (69)
- python元祖 (57)