Python+SymPy实现秒解微积分详解
之前我们分享过很多有用有趣的Python库，今天继续介绍一个：
SymPy 是一个Python库，专注于符号数学，它的目标是成为一个全功能的计算机代数系统，同时保持代码简洁、易于理解和扩展。
举一个简单的例子，比如说展开二次方程：
from sympy import *
x = Symbol('x')
y = Symbol('y')
d = ((x+y)**2).expand()
print(d)
# 结果：x**2 + 2*x*y + y**2你可以随便输入表达式，即便是十次方，它都能轻易的展开，非常方便：
from sympy import *
x = Symbol('x')
y = Symbol('y')
d = ((x+y)**10).expand()
print(d)
# 结果：x**10 + 10*x**9*y + 45*x**8*y**2 + 120*x**7*y**3 + 210*x**6*y**4 + 252*x**5*y**5 + 210*x**4*y**6 + 120*x**3*y**7 + 45*x**2*y**8 + 10*x*y**9 + y**10下面就来讲讲这个模块的具体使用方法和例子。
1.准备
请选择以下任一种方式输入命令安装依赖：
1. Windows 环境 打开 Cmd (开始-运行-CMD)。
2. MacOS 环境 打开 Terminal (command+空格输入Terminal)。
3. 如果你用的是 VSCode编辑器 或 Pycharm，可以直接使用界面下方的Terminal.
pip install Sympy
2.基本使用
简化表达式(化简)
sympy支持三种化简方式，分别是普通化简、三角化简、指数化简。
普通化简 simplify( )：
from sympy import *
x = Symbol('x')
d = simplify((x**3 + x**2 - x - 1)/(x**2 + 2*x + 1))
print(d)
# 结果：x - 1三角化简 trigsimp( )：
from sympy import *
x = Symbol('x')
d = trigsimp(sin(x)/cos(x))
print(d)
# 结果：tan(x)指数化简 powsimp( )：
from sympy import *
x = Symbol('x')
a = Symbol('a')
b = Symbol('b')
d = powsimp(x**a*x**b)
print(d)
# 结果：x**(a + b)解方程 solve()
第一个参数为要解的方程，要求右端等于0，第二个参数为要解的未知数。
如一元一次方程：
from sympy import *
x = Symbol('x')
d = solve(x * 3 - 6, x)
print(d)
# 结果：[2]二元一次方程：
from sympy import *
x = Symbol('x')
y = Symbol('y')
d = solve([2 * x - y - 3, 3 * x + y - 7],[x, y])
print(d)
# 结果：{x: 2, y: 1}求极限 limit()
dir=’+’表示求解右极限，dir=’-‘表示求解左极限：
from sympy import *
x = Symbol('x')
d = limit(1/x,x,oo,dir='+')
print(d)
# 结果：0
d = limit(1/x,x,oo,dir='-')
print(d)
# 结果：0求积分 integrate( )
先试试求解不定积分：
from sympy import *
x = Symbol('x')
d = integrate(sin(x),x)
print(d)
# 结果：-cos(x)再试试定积分：
from sympy import *
x = Symbol('x')
d = integrate(sin(x),(x,0,pi/2))
print(d)
# 结果：1求导 diff()
使用 diff 函数可以对方程进行求导：
from sympy import *
x = Symbol('x')
d = diff(x**3,x)
print(d)
# 结果：3*x**2
d = diff(x**3,x,2)
print(d)
# 结果：6*x解微分方程 dsolve( )
以 y′=2xy 为例:
from sympy import *
x = Symbol('x')
f = Function('f')
d = dsolve(diff(f(x),x) - 2*f(x)*x,f(x))
print(d)
# 结果：Eq(f(x), C1*exp(x**2))3.实战一下
今天群里有同学问了这个问题，“大佬们，我想问问，如果这个积分用Python应该怎么写呢，谢谢大家”：
# Python 实用宝典
from sympy import *
x = Symbol('x')
y = Symbol('y')
d = integrate(x-y, (y, 0, 1))
print(d)
# 结果：x - 1/2为了计算这个结果，integrate的第一个参数是公式，第二个参数是积分变量及积分范围下标和上标。
运行后得到的结果便是 x - 1/2 与预期一致。
如果大家也有求解微积分、复杂方程的需要，可以试试sympy，它几乎是完美的存在。
到此这篇关于Python+SymPy实现秒解微积分详解的文章就介绍到这了,更多相关Python SymPy微积分内容请搜索脚本之家以前的文章或继续浏览下面的相关文章希望大家以后多多支持脚本之家！
相关文章
今天小编就为大家分享一篇使用python+whoosh实现全文检索，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧2019-12-12
pytest解读fixtures之Teardown处理yield和addfinalizer方案
这篇文章主要为大家介绍了pytest解读fixtures之Teardown处理yield和addfinalizer的方案实例，有需要的朋友可以借鉴参考下，希望能够有所帮助，祝大家多多进步，早日升职加薪2022-06-06
这篇文章主要为大家详细介绍了pygame实现飞机大战，文中示例代码介绍的非常详细，具有一定的参考价值，感兴趣的小伙伴们可以参考一下2020-03-03
这篇文章主要介绍了python plt如何保存为emf图像问题,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教2023-09-09
这篇文章主要介绍了Python列表与元组的异同详解,“列表(list)与元组(tuple)两种数据类型有哪些区别”这个问题在初级程序员面试中经常碰到，超出面试官预期的答案往往能加不少印象分，也会给后续面试顺利进行提供一定帮助,需要的朋友可以参考下2019-07-07
requests.gPython 用requests.get获取网页内容为空 ’ ’问题
这篇文章主要介绍了requests.gPython 用requests.get获取网页内容为空 ’ ’，温行首先举例说明，具有一定得参考价值，需要的小伙伴可以参考一下2022-01-01
这篇文章主要介绍了python多线程并发实例及其优化,threading是扩展模块，在thread的基础上进行了封装及改进。所以只需要使用threading这个模块就能完成并发的测试,需要的朋友可以参考下2019-06-06
这篇文章主要介绍了python工具dtreeviz决策树可视化和模型可解释性，决策树是梯度提升机和随机森林的基本构建块，在学习这些模型的工作原理和模型可解释性时，可视化决策树是一个非常有帮助，下文相关资料，需要的小伙伴可任意参考一下2022-03-03
基于本地知识的问答机器人langchain-ChatGLM 大语言模型实现方法详解
这篇文章主要介绍了基于本地知识的问答机器人langchain-ChatGLM 大语言模型实现方法,结合具体实例形式详细分析了langchain-ChatGLM的功能、原理、部署方法与操作注意事项,需要的朋友可以参考下2023-07-07
这篇文章主要介绍了django_orm查询性能优化方法，小编觉得挺不错的，现在分享给大家，也给大家做个参考。一起跟随小编过来看看吧2018-08-08
最新评论