如何求解这组涉及积分的方程组?

计算科学 一体化
2021-11-26 11:55:25

我有以下一组方程:

x(t)=x0ψ,y(t)=κlnψx0ψ+1,z(t)=κlnψ,

tt0=ψ0ψdηη(1+κlnη+x0η).
这里t代表时间。我想在 Python 中为时间范围绘制这些方程t[0,1000]. 当然,绘图x,y要么z应该没问题。

但是我正在努力生成这个情节。我不希望有人真正为我做这件事,但我能得到一个关于如何解决这个问题的指针吗?特别是,我应该如何生成所需的ψ对于给定的这个问题ψ0,x0κ?

1个回答

求解微分方程总是比求解积分方程更容易。

您可以轻松区分最后一个等式 wrt 时间变量t, 并设置初始条件ψ(0)=ψ0对于以下微分方程:

dψdt=ψ(1+κln(ψ)+x0ψ).

一旦有了上述解决方案,请将其替换为您的轨迹x,yz