Differential equations solver for a chemostat

7 次查看（过去 30 天）

Kurt 2019-4-22

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/457823-differential-equations-solver-for-a-chemostat

mu=0.72;  %Constant
Ks=0.4;   %Constant
Yxs=0.44; %Constant
 
D(0)=0.6;                     %D @ t=0 (range 0.6 to 0.25)
D(t)=((mu*S(t))/(Ks+S(t)));               %D @ t=t
X(0)=0.13                       %X @ t=0 (range 0.13 to 1.67)
X(t)=(Yxs*(S(0)-((Ks*D(t))/(mu-D(t)))));  %X @ t=t
S(0)=3.7;                        %S @ t=0  (range 3.7 to .2)
S(t)=((D(t)*Ks)/(mu-D(t)));               %S @ t=t
 
ode1==dX/dt = (((mu*S)/(Ks+S))-D)*X);
ode2==dS/dt = (D*S(1)-(X/Yxs)*((mu*S)/(Ks+S))-D*S);
tspan = [0 10]; %time range, cond1
y0 = [0.2 ; 3.7]; %range of S, cond2
[T, Y] = ode45(odes,tspan,y0);
plot(t,X,t,S);
xSol(t) = dsolce(ode1,cond1,ode2,cond2);