如何解决如何用数值方法求解该系统的任意数量的微分方程?
我如何解决一个由k个微分方程组成的系统,并且每个方程中都包含导数?我正在尝试使用Scipy的resolve_ivp。
所有等式均具有以下形式:
如何使用任何求解器对方程组进行数值求解?使用solve_ivp,看来您应该能够独立于其他方程式编写每个方程式,在这种情况下,当我们有两个以上的方程式时,这似乎是不可能的。
解决方法
如果设置C[i]=B[i,i]
,则可以将方程式变换为线性系统B*z'=A
。可以解决
zdot = numpy.linalg.solve(B,A)
使得导数是一个恒定线性系统的常数解,而z
的结果解是线性z(t)=z(0)+zdot*t
。
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。