Loop to iterate ode45 function

8 Ansichten (letzte 30 Tage)
Hami the Penguin
Hami the Penguin am 7 Mär. 2021
Kommentiert: Hami the Penguin am 7 Mär. 2021
Hello, I'm having some issue using the ode45 function. In the function F, I assigned e = 0, but I want to make a loop that assigns e = 0:0.2:1 and use the ode45 function to plot the equation for each e in one plot. The code below is what I currently have, but I keep on getting errors when I try looping the function. Any ideas?
%e = 0:0.2:1;
[t,x]=ode45(@F,[0,20],[0,1]);
plot(t,x(:,1))
function xp = F(t,x)
e = 0;
xp = zeros(2,1);
xp(1) = x(2);
xp(2) = -x(1) - e*x(1)^3;
end

Akzeptierte Antwort

Jan
Jan am 7 Mär. 2021
Bearbeitet: Jan am 7 Mär. 2021
figure;
axes('NextPlot', 'add');
for e = 0:0.2:1
[t, x] = ode45(@(t, y) F(t, y, e), [0, 20], [0, 1]);
plot(t, x(:,1))
end
function xp = F(t, x, e)
xp = zeros(2,1);
xp(1) = x(2);
xp(2) = -x(1) - e*x(1)^3;
end
Here you create an "anonymous function", which appends the 3rd input argument e.

Weitere Antworten (0)

Kategorien

Mehr zu Programming finden Sie in Help Center und File Exchange

Tags

Produkte


Version

R2019b

Community Treasure Hunt

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

Start Hunting!

Translated by