How do I use ode45 for descriptive form?

1 Ansicht (letzte 30 Tage)
Tawsif Khan
Tawsif Khan am 11 Aug. 2014
Kommentiert: Tawsif Khan am 12 Aug. 2014
What would be the most efficient way to solve,
A1x'(t) = A2x(t), where both A1 and A2 are nxn matrices.
Both are sparse matrices and hence I want to avoid inversion.

Akzeptierte Antwort

Yu Jiang
Yu Jiang am 11 Aug. 2014
Bearbeitet: Yu Jiang am 11 Aug. 2014
You may solve it as a DAE (differential algebraic equation), instead of an ODE.
This documentation (Link) mentions the DAE in the form of
M(t,y)y′ = f(t,y)
Your example is a special case of this form. It can be solved by ode15s and ode23t solvers. Using those solvers, you can directly specify the M matrix without the need to invert it.
Basically, you define a function as
function dx = mySys(t,x)
dx = A2*x;
end
Then, before you solve it using ode solvers ode15s and ode23t, specify
options = odeset('Mass',A1);
Next, apply the option when using the solver as follows
[t,y] = ode15s(@mySys,tspan,y0,options);
More examples can be found in the documentation .
  3 Kommentare
Yu Jiang
Yu Jiang am 12 Aug. 2014
Bearbeitet: Yu Jiang am 12 Aug. 2014
I think it is because they use different numerical methods. See the following link for details:
-Yu
Tawsif Khan
Tawsif Khan am 12 Aug. 2014
Thanks, this really helped.

Melden Sie sich an, um zu kommentieren.

Weitere Antworten (0)

Community Treasure Hunt

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

Start Hunting!

Translated by