Python数据处理篇之Sympy系列(五)---解方程


Posted in Python onOctober 12, 2019

前言

sympy不仅在符号运算方面强大,在解方程方面也是很强大。

本章节学习对应官网的:Solvers

官方教程

https://docs.sympy.org/latest/tutorial/solvers.html

(一)求解多元一次方程-solve()

1.说明:

解多元一次方程可以使用solve(),在sympy里,等式是用Eq()来表示,

例如:2x=42x=4 表示为:Eq(x*2, 4)

2.源代码:

"""
 解下列二元一次方程
 2x-y=3
 3x+y=7
"""
# 导入模块
from sympy import *
# 将变量符号化
x = Symbol('x')
y = Symbol('y')
z = Symbol('z')
# 解一元一次方程
expr1 = x*2-4
r1 = solve(expr1, x)
r1_eq = solve(Eq(x*2, 4), x)
print("r1:", r1)
print("r1_eq:", r1_eq)
# 解二元一次方程
expr2 = [2*x-y-3, 3*x+y-7]
r2 = solve(expr2, [x, y])
print("r1:", r2)
# 解三元一次方程
f1 = x+y+z-2
f2 = 2*x-y+z+1
f3 = x+2*y+2*z-3
r3 = solve([f1, f2, f3], [x, y, z])
print("r3:", r3)

3.输出:

Python数据处理篇之Sympy系列(五)---解方程

(二)解线性方程组-linsolve()

1.说明:

在sympy中,解线性方程组有三种形式:

默认等式为0的形式:linsolve(eq, [x, y, z])

矩阵形式:linsolve(eq, [x, y, z])

增广矩阵形式:linsolve(A,b, x, y, z)

2.源代码:

"""
  x+y+z-2=0
  2x-y+z+1=0
  x+2y+2z-3=0
"""
from sympy import *
x, y, z = symbols("x y z")
# 默认等式为0的形式
print("======默认等式为0的形式 =======")
eq = [x+y+z-2, 2*x-y+z+1, x+2*y+2*z-3]
result = linsolve(eq, [x, y, z])
print(result)
print(latex(result))
# 矩阵形式
print("======矩阵形式 =======")
eq = Matrix(([1, 1, 1, 2], [2, -1, 1, -1], [1, 2, 2, 3]))
result = linsolve(eq, [x, y, z])
print(result)
print(latex(result))
# 增广矩阵形式
print("======增广矩阵形式 =======")
A = Matrix([[1, 1, 1], [2, -1, 1], [1, 2, 2]])
b = Matrix([[2], [-1], [3]])
system = A, b
result = linsolve(system, x, y, z)
print(result)
print(latex(result))

3.输出:

Python数据处理篇之Sympy系列(五)---解方程

(三)解非线性方程组-nonlinsolve()

1.说明:

nonlinsolve()用于求解非线性方程组,例如二次方,三角函数,,,等方程

2.源代码:

"""
  x**2+y**2-2=0
  x**3+y**3=0
"""
import sympy as sy
x, y = sy.symbols("x y")
eq = [x**2+y**3-2, x**3+y**3]
result = sy.nonlinsolve(eq, [x, y])
print(result)
print(sy.latex(result))

3.输出:

Python数据处理篇之Sympy系列(五)---解方程

Python数据处理篇之Sympy系列(五)---解方程

(四)求解微分方程-dsolve()

1.说明:

求解微分方程使用dsolve(),注意:

f = symbols('f', cls=Function)的作用是声明f()是一个函数。

2.源代码:

from sympy import *
# 初始化
x = symbols('x')
f = symbols('f', cls=Function)
# 表达式
expr1 = Eq(f(x).diff(x, x) - 2*f(x).diff(x) + f(x), sin(x))
# 求解微分方程
r1 = dsolve(expr1, f(x))
print(r1)
print("原式:", latex(expr1))
print("求解后:", latex(r1))

3.输出:

原式:

f(x)−2ddxf(x)+d2dx2f(x)=sin(x)
f(x)−2ddxf(x)+d2dx2f(x)=sin⁡(x)

解微分后:

f(x)=(C1+C2x)ex+cos(x)2
f(x)=(C1+C2x)ex+cos⁡(x)2

Python数据处理篇之Sympy系列(五)---解方程

总结

以上所述是小编给大家介绍的Python数据处理篇之Sympy系列(五)---解方程,希望对大家有所帮助,如果大家有任何疑问请给我留言,小编会及时回复大家的。在此也非常感谢大家对三水点靠木网站的支持!
如果你觉得本文对你有帮助,欢迎转载,烦请注明出处,谢谢!

Python 相关文章推荐
python网络爬虫采集联想词示例
Feb 11 Python
Python类属性与实例属性用法分析
May 09 Python
Django rest framework实现分页的示例
May 24 Python
Python找出微信上删除你好友的人脚本写法
Nov 01 Python
使用Django2快速开发Web项目的详细步骤
Jan 06 Python
Python反爬虫技术之防止IP地址被封杀的讲解
Jan 09 Python
python实现QQ邮箱/163邮箱的邮件发送
Jan 22 Python
python仿evething的文件搜索器实例代码
May 13 Python
Python使用scrapy爬取阳光热线问政平台过程解析
Aug 14 Python
tensorflow入门:TFRecordDataset变长数据的batch读取详解
Jan 20 Python
python属于跨平台语言码
Jun 09 Python
Python创建简单的神经网络实例讲解
Jan 04 Python
详解Python绘图Turtle库
Oct 12 #Python
Python中的list与tuple集合区别解析
Oct 12 #Python
Python 生成器,迭代,yield关键字,send()传参给yield语句操作示例
Oct 12 #Python
Python 类,property属性(简化属性的操作),@property,property()用法示例
Oct 12 #Python
python获取全国城市pm2.5、臭氧等空气质量过程解析
Oct 12 #Python
Python 迭代,for...in遍历,迭代原理与应用示例
Oct 12 #Python
Python 网络编程之TCP客户端/服务端功能示例【基于socket套接字】
Oct 12 #Python
You might like
Windows下XDebug 手工配置与使用说明
2010/07/11 PHP
php 面试碰到过的问题 在此做下记录
2011/06/09 PHP
php表单习惯用的正则表达式
2017/10/11 PHP
Gambit vs ForZe BO3 第一场 2.13
2021/03/10 DOTA
js parentElement和offsetParent之间的区别
2010/03/23 Javascript
JS获得选取checkbox整行数据的方法
2015/01/28 Javascript
javascript实现倒计时并弹窗提示特效
2015/06/05 Javascript
JS右下角广告窗口代码(可收缩、展开及关闭)
2015/09/04 Javascript
初步使用bootstrap快速创建页面
2016/03/03 Javascript
封装的dialog插件 基于bootstrap模态对话框的简单扩展
2016/08/10 Javascript
jQuery 常见小例汇总
2016/12/14 Javascript
JS实现搜索框文字可删除功能
2016/12/28 Javascript
微信小程序的生命周期的详解
2017/10/19 Javascript
JavaScript基于对象方法实现数组去重及排序操作示例
2018/07/10 Javascript
javascript中innerHTML 获取或替换html内容的实现代码
2020/03/17 Javascript
[33:17]OG vs VGJ.T 2018国际邀请赛小组赛BO2 第二场 8.18
2018/08/19 DOTA
学习python (1)
2006/10/31 Python
Python实现数据库编程方法详解
2015/06/09 Python
使用Python读写及压缩和解压缩文件的示例
2016/07/08 Python
python dict 字典 以及 赋值 引用的一些实例(详解)
2017/01/20 Python
python机器学习理论与实战(四)逻辑回归
2018/01/19 Python
python编程使用selenium模拟登陆淘宝实例代码
2018/01/25 Python
python增加矩阵维度的实例讲解
2018/04/04 Python
详解python中__name__的意义以及作用
2019/08/07 Python
Python2.6版本pip安装步骤解析
2020/08/17 Python
HTML5如何实现元素拖拽
2016/03/11 HTML / CSS
The Hut德国站点:时装、家居用品、美容等
2016/09/23 全球购物
俄罗斯金苹果网上化妆品和香水商店:Goldapple
2019/12/01 全球购物
英文自我鉴定
2013/12/10 职场文书
银行门卫岗位职责
2013/12/29 职场文书
女方回门宴答谢词
2014/01/14 职场文书
运动会入场解说词
2014/02/07 职场文书
中学生社会实践教育活动总结
2015/05/06 职场文书
Golang Gob编码(gob包的使用详解)
2021/05/07 Golang
python 利用PyAutoGUI快速构建自动化操作脚本
2021/05/31 Python
游戏《我的世界》澄清Xbox版暂无计划加入光追
2022/04/03 其他游戏