Need to do ode45 for various values iteratively

2 次查看（过去 30 天）
Muhammad Usman
Muhammad Usman 2014-7-21
编辑： Mischa Kim 2014-7-21
Hi i am writing a code for ode45 for which i have to variate the value of "aplha",after each time the ode45 evaluates and the use the next value. Say i have to use alpha=10:15
function xp=check(t,x)
xp=zeros(2,1);
xp(1)=x(2);
xp(2)=alpha*x(2)-x(1);
end
running in another m-file
[T,Y]=ode45(@check,[0 0.28],[0.087 0]);
plot(T,Y(:,1),'-o')
grid on;
i need the plot of all in the iteration in a same plot,please help me Thanks

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

采纳的回答

Mischa Kim
Mischa Kim 2014-7-21
编辑：Mischa Kim 2014-7-21
Muahmmad, you are almost there. Use a loop and something like
function myode()
alpha = 10:15;
figure
hold all
for ii = 1:numel(alpha)
[T,Y] = ode45(@check,[0 0.28],[0.087 0],[],alpha(ii));
plot(T,Y(:,1),'-o')
end
hold off
grid on;
end
function xp = check(t,x,alpha)
xp = zeros(2,1);
xp(1) = x(2);
xp(2) = alpha*x(2) - x(1);
end

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by