常微分方程的非线性系统

时间:2016-11-26 02:11:22

标签: numerical-methods differential-equations

是否有可能在给定最终状态的情况下求解系统的初始状态,以及描述系统各方面变化率的非线性微分方程系统?

例如:

V_x' = n * cos(b)
V_y' = (n * sin(b)) - (g_0 * (r_0 / r)^2)
b' = ((g_0 * (r_0 / r)^2) * cos(b)) / ||V||

b_final = 0
V_y_final = 0
V_x_final = √((G×M)/r^2)
  • V_x是水平速度
  • V_y是垂直速度
  • b是地平线以上的角度
  • n是推力/重量比
  • r_0是行星的半径

1 个答案:

答案 0 :(得分:1)

是的,它就像向前整合一样容易。事实上,对于y'=f(t,y)t in [t0,tf] y(tf)=yf的{​​{1}},您可以

  1. )通过考虑z(s)=y(tf-s)来正式改变时间方向,然后z'(s)=-f(tf-s,z(s))满足微分方程T, Z = ode45(odefunc, [tf, t0], yf),或

  2. )采用负步长,在大多数标准积分器中,可以通过按照<div id="icons">的降序传递采样时间列表来实现。我在文档中发现的唯一限制是时间列表必须是一个方向,列表中没有反转。