How do I use ode45 for descriptive form?

2 views (last 30 days)
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.

Accepted Answer

Yu Jiang
Yu Jiang on 11 Aug 2014
Edited: Yu Jiang on 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 Comments
Yu Jiang
Yu Jiang on 12 Aug 2014
Edited: Yu Jiang on 12 Aug 2014
I think it is because they use different numerical methods. See the following link for details:
-Yu

Sign in to comment.

More Answers (0)

Community Treasure Hunt

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

Start Hunting!