python数学建模导论2.2 常微分方程符号解与数值解的计算（未放入章节为微分方程介绍章节，自行百度）

最新推荐文章于 2024-09-03 11:31:52 发布

weixin_42787080

最新推荐文章于 2024-09-03 11:31:52 发布

阅读量398

点赞数

分类专栏： tool for python 文章标签： python numpy 开发语言动态规划算法

本文链接：https://blog.csdn.net/weixin_42787080/article/details/128951735

版权

tool for python 专栏收录该内容

10 篇文章 0 订阅

订阅专栏

文章展示了如何利用Python的sympy库来解决二阶线性微分方程。首先定义了符号变量x和y，然后构建了一个方程，并用dsolve函数求解。解出的方程包含了常数C1和C2，以及指数函数和三角函数的组合。接着计算了解的一阶导数，用于进一步分析。

摘要由CSDN通过智能技术生成

在这里插入图片描述

from sympy import *
y = symbols('y', cls=Function)
x = symbols('x')
eq = Eq(y(x).diff(x,2)+2*y(x).diff(x,1)+y(x),x*x)
print(dsolve(eq, y(x)))

结果
在这里插入图片描述

from sympy import *
y = symbols('y', cls=Function)
x = symbols('x')
eq = Eq(y(x).diff(x,2)+4*y(x).diff(x,1)+29*y(x), 0)
print(dsolve(eq, y(x)))

C1=symbols('C1')
C2=symbols('C2')
f=(C1*sin(5*x) + C2*cos(5*x))*exp(-2*x)    #第五行eq的解
print(f.diff(x,1)) #eq的一阶导的形式

'''
(0*C1+1*C2)*1=0
-2*(C1*0+C2*1)*1+(5*C1*1-5*C2*0)*1=15
'''