最新下载
热门教程
- 1
- 2
- 3
- 4
- 5
- 6
- 7
- 8
- 9
- 10
python中sympy库求常微分方程的用法代码解析
时间:2020-04-28 编辑:袖梨 来源:一聚教程网
本篇文章小编给大家简单分享一下python中sympy库求常微分方程的用法代码解析,文章代码介绍的很详细,小编觉得挺不错的,现在分享给大家供大家参考,有需要的小伙伴们可以来看看。
问题1:
程序,如下
from sympy import * f = symbols('f', cls=Function) x = symbols('x') eq = Eq(f(x).diff(x, x) - 2*f(x).diff(x) + f(x), sin(x)) print(dsolve(eq, f(x)))
结果
Eq(f(x), (C1 + C2*x)*exp(x) + cos(x)/2)
附:布置考试中两题
1.利用python的Sympy库求解微分方程的解y=f(x),并尝试利用matplotlib绘制函数图像
程序,如下
from sympy import * f = symbols('f', cls=Function) x = symbols('x') eq = Eq(f(x).diff(x,1)+f(x)+f(x)**2, 0) print(dsolve(eq, f(x))) C1 = symbols('C1') eqr = -C1/(C1 - exp(x)) eqr1 = eqr.subs(x, 0) print(solveset(eqr1 - 1, C1)) eqr2 = eqr.subs(C1, 1/2) # 画图 import matplotlib.pyplot as plt import numpy as np x_1 = np.arange(-5, 5, 0.1) y_1 = [-0.5/(0.5 - exp(x)) for x in x_1] plt.plot(x_1, y_1) plt.axis([-6,6,-10,10]) plt.grid() plt.show()
结果
Eq(f(x), -C1/(C1 - exp(x)))
FiniteSet(1/2)
2.利用python的Sympy库求解微分方程的解y=y(x),并尝试利用matplotlib绘制函数图像
程序,如下
from sympy import * y = symbols('y', cls=Function) x = symbols('x') eq = Eq(y(x).diff(x,1), y(x)) print(dsolve(eq, y(x))) C1 = symbols('C1') eqr = C1*exp(x) eqr1 = eqr.subs(x, 0) print(solveset(eqr1 - 1, C1)) eqr2 = eqr.subs(C1, 1) # 画图 import matplotlib.pyplot as plt import numpy as np x_1 = np.arange(-5, 5, 0.01) y_1 = [exp(x) for x in x_1] plt.plot(x_1, y_1, color='orange') plt.grid() plt.show()
结果
Eq(y(x), C1*exp(x))
FiniteSet(1)
相关文章
- Golang ProtoBuf的基本语法详解 10-20
- Python识别MySQL中的冗余索引解析 10-20
- Python+Pygame绘制小球代码展示 10-18
- Python中的数据精度问题介绍 10-18
- Python随机值生成的常用方法介绍 10-18
- python3解压缩.gz文件分析 09-27