sympy求微分方程

孙士东 2019-08-21 11:04:39

求大神帮忙：

dφ/dθ=-k/ω/V*(ω*φ*dV/dθ-sqrt(2*κ*R*Ts/(κ-1))*N*Av*av*h/H*φ**(1/κ)*sqrt(1-φ**((κ-1)/κ)))

dV/dθ=Vs*(sinθ/2+λ/4*sin2θ/sqrt(1-λ**2*(sinθ)**2)

以上是微分方程组，其他数都是变量，自行输入。

想求φ=f(θ)在0~180度范围的变化曲线。

...全文

126 1 打赏收藏转发到动态举报

写回复

1 条回复

切换为时间正序

请发表友善的回复…

发表回复

祥瑞Coding 2020-03-30

打赏
举报

回复

我也比较想知道，如何用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

diff(f,x)diff(f, x)diff(f,x)求导数可引入求微分方程 sympy 求微分方程.(点击可跳转) 1.一阶导数基本格式 print(diff(f, x)) # f为所求导函数，x为对x进行求导例：求该函数的导数 f(x)=cos(x)f(x)=cos(x)f(x)=cos(x) 程序，如下 from sympy import * x = symbols('x') print(diff(cos(x), x)) 结果：-sin(x) ~ 2.二阶导数例：求该函数的二阶导数 f(x)=cos2(x)f(x)=cos^{2}(x)f(x)=cos2(x) 程序，如

差异性 微分方程求解器和Sympy Grapher

sympy是python的一个用于科学计算的库，可解方程和方程组，计算极限，积分，微分，微积分，矩阵，级数等

最优控制理论广泛应用于对交通工程、经济管理等领域的复杂动态系统优化。由于最优控制模型包括控制决策、状态、协态乘子等诸多变量，通常在求解上十分复杂。在对最优控制模型基本概念介绍基础上，重点针对固定终端时间下的三类问题，通过构建哈密顿函数，得到相应的最优解必要条件，推导出状态与协态所满足的微分方程组，进一步利用Python的符号工具库(Sympy)，对微分方程组以及系数进行了求解，最后绘制了相应的最优控制策略以及状态轨迹图，案例求解为控制类专业相关课程的教学提供了素材。

37,720

社区成员

34,239

社区内容

发帖

与我相关

我的任务

社区管理员

加入社区

近7日
近30日
至今

加载中

查看更多榜单

试试用AI创作助手写篇文章吧

+ 用AI写文章