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

Antworten (1)

Jon
Jon am 12 Nov. 2020

1 Stimme

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 Kommentare

Zachary Braida
Zachary Braida am 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 am 12 Nov. 2020
Bearbeitet: James Tursa am 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 am 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!

Melden Sie sich an, um zu kommentieren.

Kategorien

Mehr zu Programming finden Sie in Hilfe-Center und File Exchange

Produkte

Version

R2020a

Gefragt:

am 12 Nov. 2020

Kommentiert:

am 12 Nov. 2020

Community Treasure Hunt

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

Start Hunting!

Translated by