如果我有两组高度耦合的 ODE,例如:
dy/dt = A*dx/dt + B*y + C*x
dx/dt = D*dy/dt + E*y + F*x
或者如果我有这样的系统:
U = function(dy/dt)
dy/dt = function(U,...)
ode45 通常会如何处理这个问题?如果我要执行以下操作:
DY(1) = A*DY(2) + B*y(1) + C*y(2)
DY(2) = D*DY(1) + E*y(1) + F*y(2)
或类似的第二个问题:
U = A*DY(1) + ...
DY(1) = B*U + ...
并运行 ode45,MATLAB 会毫无怨言地产生结果,但我有点犹豫,因为我不确定它是如何处理这样一个耦合方程组的。和建议?