Loop to iterate ode45 function
3 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Hami the Penguin
le 7 Mar 2021
Commenté : Hami the Penguin
le 7 Mar 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
0 commentaires
Réponse acceptée
Jan
le 7 Mar 2021
Modifié(e) : Jan
le 7 Mar 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.
Plus de réponses (0)
Voir également
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!