For loop data extraction

I am having trouble finding a way to extract data at then end of a for loop with function.
for i = 1:length(height)
% solve coupled drag ODE
options = odeset('RelTol',1e-5,'AbsTol',1e-7); % ODE options
f = @(t,x) horizEqns(Isp,Wo,Wb,Wdot,msl,1,height(i),t,x); % the equations that will be solved
[t,soln] = ode45(f,[0 maxtime],[0. v0],options);
end

Respuestas (1)

Jon
Jon el 12 de Nov. de 2020

1 voto

You overwrite your t and soln values with each loop iteration. If you want to save those to use them later you need to put them into arrays. Assuming t and soln have the same length each time
[t(:,i),soln(:,i)] = ode45(f,[0 maxtime],[0. v0],options);

3 comentarios

Zachary Braida
Zachary Braida el 12 de Nov. de 2020
I initially tried that but get this as the error:
Unable to perform assignment because the indices on the left side are not compatible
with the size of the right side.
James Tursa
James Tursa el 12 de Nov. de 2020
Editada: James Tursa el 12 de Nov. de 2020
Probably because soln has more than one column.
Try this
[t,soln(:,:,i)] = ode45(f,[0 maxtime],[0. v0],options);
or this
[t,soln{i}] = ode45(f,[0 maxtime],[0. v0],options);
Zachary Braida
Zachary Braida el 12 de Nov. de 2020
Thank you! That second one works. If you add it as an answer outside the comment, I will accept it as the working!

Iniciar sesión para comentar.

Categorías

Más información sobre Programming en Centro de ayuda y File Exchange.

Productos

Versión

R2020a

Preguntada:

el 12 de Nov. de 2020

Comentada:

el 12 de Nov. de 2020

Community Treasure Hunt

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

Start Hunting!

Translated by