Using ODE 45 coupled ODEs

Question

Andreas Skovhøj 2019-10-23

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/486981-using-ode-45-coupled-odes

回答： Star Strider 2019-10-23

采纳的回答： Star Strider

在 MATLAB Online 中打开

Hi!

I am supposed to solve these two linear ODEs

k = 0.0440;
C_A0 = 0.3045;
epsilon = 2;
alpha = 0.001;
F_A0 = 2.5;
rho = 0.001;
ode1 = diff(y) == -(alpha/2*y)*(1+epsilon*X)*rho;
ode2 = diff(X) == k*C_A0*(1-X)/(F_A0*1+epsilon*X)*y*rho;

I cant figure out how to use ode45 but i have tried something like this by using some other answers which i dont understand..

function dz = myode2(v,z)
syms v z
alpha = 0.001;
C0 = 0.3;
esp = 2;
k = 0.044;
f0 = 2.5;
dz = zeros(2,1);
dz(1) = k*C0/f0*(1-z(1)).*z(2)./(1-esp*z(1));
dz(2) = -alpha*(1+esp*z(1))./(2*z(2));
ode45(@myode2,[0 500],[0 1])
end

This does not work tho..

Anyone knows how to solve this by ode45?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Star Strider 2019-10-23

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/486981-using-ode-45-coupled-odes#answer_397865

在 MATLAB Online 中打开

Start with what you already have:

syms X(t) y(t) Y 
k = 0.0440;
C_A0 = 0.3045;
epsilon = 2;
alpha = 0.001;
F_A0 = 2.5;
rho = 0.001;
ode1 = diff(y) == -(alpha/2*y)*(1+epsilon*X)*rho;
ode2 = diff(X) == k*C_A0*(1-X)/(F_A0*1+epsilon*X)*y*rho;