方程如下:
dXp/dt=μm*S*Xp/(Ks+S)-k1*Xp (1)
dXs/dt= k1*Xp-kd*Xs (2)
X=Xp+Xs (3)
dX/dt=μm*S*Xp/(Ks+S)-kd*Xs (4)
dEt/dt=k2*Xs-k3*Et (5)
-dS/dt=(k4*(1-x)^k5+k6)*0.9062*S*Ef/(1+0.9062*Ef) (6)
其中式(4)是由前面(1)-(3)推导得到
方程组中X, Et, Ef, S, x为已知数据的变量(见下表),μm, kd, Ks, k1, k2, k3, k4, k5, k6为待拟合求解的参数。请高人们给帮帮忙吧,冰天雪地里跪谢啦!
t(h) 0 6 12 18 24 30 37 42 48 54 61 66 72 78 84 90 96 102 108 114 120 126 132
X 2.76 3.93 6.01 8.58 11.25 12.29 13.42 18.16 22.55 22.53 22.52 20.81 17.95 17.94 17.94 17.29 16.63 16.35 15.98 14.21 12.78 12.43 11.96
S 50 48.66 46.2 42.78 38.3 32.25 26.4 26.13 25.72 23.28 21.69 20.55 22.4 19.83 19.07 19.26 20.38 19.13 19.06 18.18 17.24 16.93 16.56
Ef 0 0.025 0.076 0.129 0.198 0.95 1.89 2.62 3.46 4.58 5.72 6.38 6.91 7.27 7.58 9.05 10.63 10.66 10.48 10.59 10.69 11.72 12.25
Et 0 0.046 0.133 0.215 0.31 1.237 2.21 2.973 3.835 4.94 6.069 6.717 7.281 7.601 7.9 9.38 10.985 10.993 10.811 10.906 10.99 12.017 12.542
x 0 0.0268 0.076 0.1444 0.234 0.355 0.472 0.4774 0.4856 0.5344 0.5662 0.589 0.552 0.6034 0.6186 0.6148 0.5924 0.6174 0.6188 0.6364 0.6552 0.6614 0.6688