MATLAB爱好者论坛-LabFans.com

MATLAB爱好者论坛-LabFans.com (https://www.labfans.com/bbs/index.php)
-   MATLAB论坛 (https://www.labfans.com/bbs/forumdisplay.php?f=6)
-   -   【求助】请高人帮我解一下下面的微分方程组,很急切,非常感谢! (https://www.labfans.com/bbs/showthread.php?t=1616)

cinderliu 2008-03-09 21:06

【求助】请高人帮我解一下下面的微分方程组,很急切,非常感谢!
 
:cry: :cry:
dx/dt=u*(1-x/xm)^n1*(1-p/pm)^n2*(1-s/sm)^n3*x;
dp/dt=a*dx/dt+b;
ds/dt=(-1/Y1)*dp/dt-(1/Y2)*dx/dt-m*s
非常感谢,这是描述微生物生长的方程组,分别代表细胞的生长,产物的生成以及底物的消耗

cinderliu 2008-03-12 11:02

大家指点一下,应该用什么求,谢谢,非常感谢!

slgu 2008-03-12 15:34

[y z w]=Dsolve('Dy=u*(1-x/xm)^n1*(1-p/pm)^n2*(1-s/sm)^n3*x','Dz=a*Dy+b','Dw=(-1/Y1)*Dz-(1/Y2)*Dy-m*s')
y =
-(t*m*s*Y1*Y2+t*a*u*(1-x/xm)^n1*(1-p/pm)^n2*(1-s/sm)^n3*x*Y2+t*b*Y2+u*(1-x/xm)^n1*(1-p/pm)^n2*(1-s/sm)^n3*x*t*Y1-C1*Y1*Y2)/Y1/Y2
z =
u*(1-x/xm)^n1*(1-p/pm)^n2*(1-s/sm)^n3*x*t+C3
w =
t*a*u*(1-x/xm)^n1*(1-p/pm)^n2*(1-s/sm)^n3*x+t*b+C2


所有时间均为北京时间。现在的时间是 05:00

Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.