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

Réponses (1)

Jon
Jon le 12 Nov 2020

1 vote

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 commentaires

Zachary Braida
Zachary Braida le 12 Nov 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 le 12 Nov 2020
Modifié(e) : James Tursa le 12 Nov 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 le 12 Nov 2020
Thank you! That second one works. If you add it as an answer outside the comment, I will accept it as the working!

Connectez-vous pour commenter.

Catégories

En savoir plus sur Programming dans Centre d'aide et File Exchange

Produits

Version

R2020a

Community Treasure Hunt

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

Start Hunting!

Translated by