Loop to iterate ode45 function

Question

Hami the Penguin 2021-3-7

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/765331-loop-to-iterate-ode45-function

评论： Hami the Penguin 2021-3-7

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 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Jan 2021-3-7

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/765331-loop-to-iterate-ode45-function#answer_641181

编辑：Jan 2021-3-7

在 MATLAB Online 中打开

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