这篇文章我们来了解“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库来解微积分”就介绍到这了，如果大家觉得不错可以参考了解看看，如果想要了解更多，欢迎关注群英网络，小编每天都会为大家更新不同的知识。

免责声明：本站发布的内容（图片、视频和文字）以原创、转载和分享为主，文章观点不代表本网站立场，如果涉及侵权请联系站长邮箱：mmqy2019@163.com进行举报，并提供相关证据，查实之后，将立刻删除涉嫌侵权内容。

猜你喜欢

• Python中Addict模块的安装和用法是什么，有哪些要点

  Addit是一个Python模块，除了提供标准的字典语法外，Addit 生成的字典的值既可以使用属性来获取，也可以使用属性进行设置。本文将详细讲讲它的使用方法，需要的可以参考一下

• Python进行Socket网络编程如何实现主机进行通话

  这篇文章主要介绍了python实现客户端和服务端之间进行通信，文章通过python利用socket展开详情介绍，具有一定的参考价值，需要的小伙伴可以参考一下

• numpy中的dot函数用法是什么？

  这篇文章主要给大家分享numpy中的dot函数用法的内容，下文有具体的实例，具有一定的借鉴和学习价值，感兴趣的朋友可以参考一下，希望大家阅读完这篇文章能有所收获，下面我们一起来学习一下吧。

• python列表的增删操作怎样实现？

  python列表的增删操作怎样实现？对列表做增删操作是比较常见的需求，也是很基础的内容，对此这篇文章就给大家分享一些python列表的增删操作的实例，感兴趣的朋友可以参考。

• 利用Python求解微分方程怎么实现，用什么方法

  这篇文章主要介绍了python 解决微分方程的操作(数值解法)，具有很好的参考价值，希望对大家有所帮助。如有错误或未考虑完全的地方，望不吝赐教