How does ode45 handle highly coupled ODEs (where y' depends on x' and x' depends on y', etc.)?

Ali
2013 5 19
1 个回答
回答已采纳
1 次查看（30 天）

0 个投票

If I have two sets of highly coupled ODEs such as:
dy/dt = A*dx/dt + B*y + C*x
dx/dt = D*dy/dt + E*y + F*x
or if I have a system like:
U = function(dy/dt)
dy/dt = function(U,...)
How would ode45 normally handle this? If I were to do the following:
DY(1) = A*DY(2) + B*y(1) + C*y(2)
DY(2) = D*DY(1) + E*y(1) + F*y(2)
or similarly for the second problem:
U = A*DY(1) + ...
DY(1) = B*U + ...
and run ode45, MATLAB produces results without complaint but I'm slightly hesitant as I'm not sure how it's handling such a coupled system of equations.
MATLAB shouldn't be able to solve this sort of problem with ode45 as it's generally an implicit problem (example number 2 provided above). It should give me an error when I try to run it, however, it runs smoothly and that's what's confusing me. Realistically I should only be able to use ode15i for this sort of problem...
And suggestions?

请先登录，再回答此问题。

 采纳的回答

Ali
Ali 2013-5-20

0 个投票

After looking through many lines of code, I realize that the reason MATLAB is accepting an implicit equation without complaint is that when you're running a system of equations, you have to initialize the vector with DY = zeros(N,1) where N is the number of equations. It's using this initial zero value for the integration and effectively giving incorrect results.

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息

产品

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by