Simultaneously numerically integrate multiple ODE's in MATLAB

4 次查看（过去 30 天）
Rong Song
Rong Song 2019-10-5
回答： Pavel Osipov 2019-10-6
I was trying to solve a system of muptiple ODEODE at the same time in Matlab but wasn't sure how should I solve for that.
So if I have
dy1/dx=sum of yi , i from 1 to j
dy2/dx=y1^2-y2*sum of yi, i from 1 to j
dy3/dx=y1*y2+y2^2- y3*sum of yi, i from 1 to j
and given initial condition that y1(0)=1 other yi=0 (i>1)
Could anyone have me with this please? Thank you!
  2 个评论
Rong Song
Rong Song 2019-10-5
K=10;
% Initial Conditio
for k=1:K
if k==1
yi(k)=1;
else
yi(k)=0;
end
end
b=0;
d=20;
% y include y1 y2 y3 y4...to yK
y=zeros(1,100);
for k=1:K
if k==1
y=e(k);
dydd{k}=@(y,t) -y(k)*sum(y);
else
y=e(k);
for j=1:k-1
b=b+y(j)*y(k-j);
end
dydd{k}=@(y,t) 1/(2*y(k))*b-y(k)*sum(y);
end
end
interval=linspace(0,d,100);
[t,y]=ode45(dydd,interval,yi);

请先登录，再进行评论。

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

回答（2 个）

Pavel Osipov
Pavel Osipov 2019-10-5
"I was trying to solve a system of multiple orders at the same time in Matlab but wasn't sure how should I solve for that". "is that right?"
Then imagine the left part as a system of equations of the 1st order and solve in Matlab in vector form. Good luck!
Pavel Osipov
Pavel Osipov 2019-10-6
ode45.PNG
odefun.PNG

类别

在 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