Is it possible to extract the time and two output variables from ODE45?
8 views (last 30 days)
Show older comments
Hello,
I have written a function which holds state space equations for ODE45 to solve, being in the nature of:
function zdot = odess(t,z)
In my mainscript which calls upon this function and has the initial conditions for the ode, I am recalling the outputs by:
[t,z] = ode45(@odess, t_step, z0, opts);
However, is it possible to get a second variable from the 2nd code listed above? For example:
[t,z,new_output] = ode45(@odess, t_step, z0, opts);
Thanks for your time.
2 Comments
Answers (2)
Torsten
on 17 Jan 2017
The easiest way is to reculculate new_output from t and z after the call
[t,z] = ode45(@odess, t_step, z0, opts);
Best wishes
Torsten.
8 Comments
Torsten
on 20 Jan 2017
If you set up your code as above, the variable z(9) is z(1), integrated over time, i.e.
z_9(t)=integral_{0}^{t} z_1(t) dt.
If you want to output z(1), just use your previous code without any changes and use the command
plot (t,z(:,1));
after the call to ode45, as I already suggested.
Your assertion
z(1) isn't included in the output matrix "z" from [t,z] = ode45(@odess,t_step,z0, opts), i think what is included are only zdots.
is simply wrong.
Best wishes
Torsten.
Jan
on 20 Jan 2017
@Tyler: An integrator as ODE45 gets the start position (or "state vector") as input, than calls the function to be integrated, which replies the derivatives and accumulates them. For each successful time step, the new position is replied. Therefore the output of the integrator contains the positions (state varaibales) and the first element is its first component.
You can simply check this manually: Computer your odess(t,z) for the inital time and start vector and check the first row of the output of the integrator.
Walter Roberson
on 17 Jan 2017
Edited: Walter Roberson
on 17 Jan 2017
No, that cannot be done. See though
0 Comments
See Also
Products
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!