loop inside ode45 solver

Question

0 个投票

Hello, I am using loop inside ode45 solver with shooting method ,and want to store data from each iteration. My code only save the last iteration. Please help.

    function F=solver(x)
    m=-2:1:-1;
    u1=zeros(length(m));
    for k=1:length(m);
    [t{k},u{k}]=ode45(@(t,u)equation(t,u,m(k)),[0 10],init1,options);
    s=length(t{k});
    F=[u{k}(s,2)-1,u{k}(s,4)];
    u1(k)=u{k}(1,3);
    end
    end
     for k=1:length(m);
       Final=[m(k),u1(k)]
     end

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Fayeza Fayeza 2017-10-27

I am still wondering why this loop saved only the last iteration? or I should call ode45 solver many times? but this is useless if there are many values for a parameter m. Please help!

请先登录，再进行评论。

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

Follow Question

Answer 1

KSSV 2017-10-24

编辑：KSSV 2017-10-24

在 MATLAB Online 中打开

0 个投票

function F=solver(x)
m=-2:1:-1;
u1=zeros(length(m));
S = cell(length(m),1) ;
for k=1:length(m);
    [t,u]=ode45(@(t,u)equation(t,u,m(k)),[0 10],init1,options);
    S{k} = u(:,3) ;
    s=length(t{k});
    F=[u{k}(s,2)-1,u{k}(s,4)];
    u1(k)=u{k}(1,3);
end
end
for k=1:length(m);
    Final=[m(k),u1(k)]
end

3 个评论
显示 1更早的评论隐藏 1更早的评论

KSSV 2017-10-24

Edited..is that you are looking for?

Fayeza Fayeza 2017-10-24

yes. After calling ode45 solver I got u with 5 columns u(n,5) with n=1 to 10. I want to edit only u(1,3) for different parameters m. Because of that I am using k loop.

请先登录，再进行评论。

loop inside ode45 solver

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

回答（1 个）

3 个评论
显示 1更早的评论隐藏 1更早的评论

类别

标签

Community Treasure Hunt

loop inside ode45 solver

1 个评论 显示 -1更早的评论 隐藏 -1更早的评论

回答（1 个）

3 个评论 显示 1更早的评论 隐藏 1更早的评论

类别

标签

另请参阅

Community Treasure Hunt

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

3 个评论
显示 1更早的评论隐藏 1更早的评论