scipy solve_ivp使用无法表示为导数的微分方程

rslzwgfq  于 2022-11-10  发布在  其他
关注(0)|答案(1)|浏览(158)

我希望解决系统
y' = y + u'
u' = y' + u
看起来solve_ivp要求我求解导数并用原始函数来表达它,但我不知道在这种情况下如何做。我不想找到显式结果,我想知道如何让solve_ivp合作,如果可能的话,这样我就可以用它来解决同类的更难的问题,这些问题是不容易处理的。

arknldoa

arknldoa1#

Your equations imply y + u = 0, so you can do a bit algebra to get y' = y/2. You don't need a numerical solver for that; the solution is

y(t) = y0*exp(t/2)

Then

u(t) = -y0*exp(t/2)

相关问题