數(shù)學(xué)建模微分方程模型 數(shù)學(xué)建模微分方程實(shí)驗(yàn)

1、數(shù)學(xué)建模微分方程模型 數(shù)學(xué)建模微分方程實(shí)驗(yàn) 2 微分方程實(shí)驗(yàn) 1、微分方程穩(wěn)定性分析繪出下列自治系統(tǒng)相應(yīng)的軌線，并標(biāo)出隨t增加的運(yùn)動(dòng)方向，確定平衡點(diǎn)，并按穩(wěn)定的、漸近穩(wěn)定的、或不穩(wěn)定的進(jìn)行分類：dxdxdxdx=x,=-x,=y,=-x+1,dtdtdtdt(1)(2)(3)(4) dydydydy=y;=2y;=-2x;=-2y.dtdtdtdt 解：（1）由f（x）=x=0，f（y）=y=0；可得平衡點(diǎn)為（0,0），1系數(shù)矩陣A=00，求得特征值11=1，2=1；p=-(1+2)=-2<0，q=12=1>0；對照穩(wěn)定性的情況表，可知平衡點(diǎn)（0，0）是不穩(wěn)定的。圖形如下： （2）

2、如上題可求得平衡點(diǎn)為（0,0），特征值1=-1，2=2；p=-(1+2)=-1<0，q=12=-2<0；對照穩(wěn)定性的情況表，可知平衡點(diǎn)（0，0）是不穩(wěn)定的。 其圖形如下： （3）如上題可求得平衡點(diǎn)為（0,0），特征值1=0 + 1.4142i，2=0 - 1.4142i； p=-(1+2)= 0，q=12=1.4142>0；對照穩(wěn)定性的情況表，可知平衡點(diǎn)（0，0）是不穩(wěn)定的。其圖形如下： （4）如上題可求得平衡點(diǎn)為（1,0），特征值1=-1，2=-2；p=-(1+2)= 3>0，q=12=2>0；對照穩(wěn)定性的情況表，可知平衡點(diǎn)（1，0）是穩(wěn)定的。 其圖形如下： 2

3、、種群增長模型一個(gè)片子上的一群病菌趨向于繁殖成一個(gè)圓菌落.設(shè)病菌的數(shù)目為N，單位成員的增長率為r1，則由Malthus生長律有dNdt=r1N，但是，處于周界表面的那些病菌由于寒冷而受到損傷，它們死亡的數(shù)量與N1/2成比例，其比例系數(shù)為r2，求N滿足的微分方程.不用求解，圖示其解族.方程是否有平衡解，如果有，是否為穩(wěn)定的?解：由題意很容易列出N滿足的微分方程：dNdt1=r1N-r2N2=f(N)令f(N)=0，可求得方程的兩個(gè)平衡點(diǎn)N1=0,N2=r22/r12 進(jìn)而求得dNdt22dNdt22=(r1-12r2N-121)(r1N-r2N2)令=0可求得N=r2/4r122則N=N1，N=

4、N2，N=r22/4r12可以把第一象限劃為三部分，且從下到上三部分中分別有dNdt0,dNdt220,；dNdt0,dNdt220；dNdt0。 則可以畫出N（t）的圖形，即微分方程的解族，如下圖所示： 由圖形也可以看出，對于方程的兩個(gè)平衡點(diǎn)，其中N1=0是不穩(wěn)定的；N2=r2/r122是穩(wěn)定的。 3、有限資源競爭模型1926年Volterra提出了兩個(gè)物種為共同的、有限的食物來源而競爭的模型 dx1=b1-l1(h1x1+h2x2)x1dtdx2=b-l(hx+hx)x2211222dt假設(shè)b1l1b2l2，稱bili為物種i對食物不足的敏感度，（1）證明當(dāng)x1(t0)>0時(shí)，物種2

5、最終要滅亡; （2）用圖形分析方法來說明物種2最終要滅亡. 解：（1）由上述方程組f（x1）=b1-l1(h1x1+h2x2)x1=0，f（x2）=b2-l2(h1x1+h2x2)x2=0，可得方程的平衡點(diǎn)為P0（0,0），P1（b2b1l1h1,0），P2（0,l2h2）.對平衡點(diǎn)P0（0,0），b1-h1l1x1-h2l1x2系數(shù)矩陣A=-h1l2x2b1=b2-h1l2x1-h1l2x20-h2l2x10 b2則p=-（b1+b2）<0，所以該平衡點(diǎn)不穩(wěn)定。對平衡點(diǎn)P1（b1l1h1,0），-b1-h2l2x1=b2-h1l2x1-h1l2x20b1l2)系數(shù)矩陣A=b1-h1l1

6、x1-h2l1x2-h1l2x2h2b1h1 blb2-12l1-則p= b1-b2+b1b2b1l2l1，q= -b1(b2-l1，由題意l1l2，x1(t0)>0，可以得出p>0,q>0,因此該平衡點(diǎn)是穩(wěn)定的。b1即t時(shí)，(x1(t),x2(t)(l1h1,0)，說明物種2最終要滅亡。對平衡點(diǎn)P2（0,b1b2b2l2h2），同理可以得到q<0，在該平衡點(diǎn)不穩(wěn)定。因此，在l1l2，x1(t0)>0的條件下，物種2最終要滅亡。b1-l1(h1x1+h2x2)=0（2）對于線性方程組b2-l2(h1x1+h2x2)=0在平面上匹配兩條直線l1和l2，由題意b1l1

7、b2l2，x1(t0)>0，可將第一象限分為三個(gè)區(qū)域。在最左邊區(qū)域，x1,x2都大于0；在中間區(qū)域，x1,x2都小于0，在最右邊區(qū)域，x1,x2分別是大于0和小于0.，由各區(qū)域中x1,x2的取值可得到如下圖形： 由圖也可以看出，隨著時(shí)間的增加，物種1最終能達(dá)到穩(wěn)定值終要滅亡。4、蝴蝶效應(yīng)與混沌解 b1l1h1，物種2最考慮Lorenz模型x1(t)=-bx1(t)+x2(t)x3(t)x2(t)=-sx2(t)+sx3(t) x(t)=-x(t)x(t)+rx(t)-x(t)12233其中=10，=28，=8/3，且初值為,x1（0）=x2（0）=0，x3（0）=，為一個(gè)小常數(shù)，假設(shè)=1

8、0-10，且0t100。（1）用函數(shù)ode45求解，并畫出x2x1,x2x3,x3x1的平面圖；（2）適當(dāng)?shù)卣{(diào)整參數(shù)，值，和初始值x1（0），x2（0）=0，x3（0），重復(fù)一的工作，看有什么現(xiàn)象發(fā)生。解：（1）編寫Lorenz函數(shù)，function xdot=lorenz1(t,x,b,a,c)xdot=-b*x(1)+x(2)*x(3);-a*x(2)+a*x(3);-x(1)*x(2)+c*x(2)-x(3);對各參數(shù)賦值并用ode45函數(shù)求解，可得數(shù)值解：Columns 1 through 90 0.1250 0.2500 0.3750 0.5000 0.5352 0.5705 0.6

9、057 0.64090 0.0000 0.0000 -0.0000 -0.0000 0.0000 0.0000 0.0000 0.00000 0.0000 0.0000 -0.0000 -0.0000 0.0000 0.0000 0.0000 -0.00000.0000 -0.0000 -0.0000 0.0000 0.0000 0.0000 -0.0000 0.0000 0.0000 Columns 10 through 180.6761 0.7114 0.7466 0.7818 0.8308 0.8797 0.9286 0.9776 1.01050.0000 0.0000 0.0000 0

10、.0000 0.0000 0.0000 0.0000 0.0000 0.00000.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.00000.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 Columns 19 through 271.0434 1.0763 1.1092 1.1421 1.1750 1.2079 1.24091.2797 1.31860.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.00

11、00 0.00000.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0001 0.00010.0000 0.0000 0.0000 0.0000 0.0000 0.0001 0.0001 0.0002 0.0002Columns 28 through 361.3575 1.3964 1.4246 1.4528 1.4810 1.5092 1.5374 1.5656 1.59380.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.00000.0002 0.0003 0.0004 0

12、.0005 0.0008 0.0011 0.0015 0.0021 0.00290.0004 0.0006 0.0008 0.0012 0.0016 0.0023 0.0045 0.0063.Columns 5590 through 559899.5751 99.5921 99.6090 99.6260 99.6462 99.6664 99.7069 99.733816.9457 16.5261 16.2010 15.9854 15.8978 16.0348 17.1933 18.7894-3.3551 -3.7119 -4.1098 -4.5568 -5.1636 -5.8601 -7.54

13、38 -8.8677-5.3476 -5.9274 -6.5941 -7.3519 -8.3766 -9.5370 -12.1944 -14.0725Columns 5599 through 560799.7608 99.7878 99.8148 99.8306 99.8465 99.8623 99.8940 99.909921.2186 24.4709 28.2629 30.5354 32.6533 34.4645 36.7149 37.0528-10.3044 -11.7459 -12.9932 -13.5361 -13.8800 -13.9854 -13.4302 -12.7919-15

14、.7415 -16.8309 -16.9274 -16.3648 -15.3302 -13.8707 -10.1013 -8.1005Columns 5608 through 561399.9257 99.9416 99.9562 99.9708 99.9854 100.0000 36.8957 36.3239 35.5248 34.5435 33.4481 32.2971 -11.9606 -10.9899 -10.0237 -9.0332 -8.0563 -7.1223 -6.2187 -4.5596 -3.2818 -2.2590 -1.4835 -0.9275x2x1,x2x3,x3x

15、1的平面圖分別如下：0.0032 99.6867 16.4476 -6.6527 -10.8209 99.8782 35.8466 -13.8351 -12.0805 （2）令參數(shù)，值各減1，初始值x1（0），x2（0）不變，x3（0）=10-8 分別得到得到x2x1,x2x3,x3x1的平面圖如下： 可以看出，參數(shù)，值各減1，初始值x1（0），x2（0）不變，x3（0）數(shù)值變?yōu)閤3（0）=10-8，參數(shù)和初始值很小的改變，就會(huì)導(dǎo)致最后圖形發(fā)生很大的變化。 5、用微分方程考察共振現(xiàn)象 設(shè)物體沿x軸運(yùn)動(dòng)(如圖所示)其平衡位置取為原點(diǎn)0，物體的質(zhì)量為1，在時(shí)間t物體的位置為x(t)其所受的恢復(fù)為(

16、如彈性力等)與物體所在位置的坐標(biāo)成正比，即k2x,其中常數(shù)k稱為恢復(fù)系數(shù)，運(yùn)動(dòng)過程所受的阻力(由于介質(zhì)及摩擦等)設(shè)與速度成正比，即2hdxdt，h>0，稱為阻尼系數(shù)。（1）根據(jù)Newton第二定律，建立相應(yīng)的微分方程.不妨設(shè)初始位置為1，初始速度為0,取k=2, h=0(當(dāng)h = 0稱為簡諧振動(dòng)的方程)和h=0.1，用Matlab軟件得到相應(yīng)的數(shù)值解，并在t-x平面上畫出x（t）的圖形。（2）如果物體還受到附加外力的干擾，且外力是一個(gè)依據(jù)時(shí)間t的函數(shù)f(t)(設(shè)f (t)=B sinwt)，建立相應(yīng)的微分方程(該方程稱為強(qiáng)迫振動(dòng)方程).在上述參數(shù)不變的情況下，取振幅B=1，分別取w=1,

17、 1.2, 1.4, 1.6, 1.8, 2.0, 2.2,2.4, 2.6, 2.8, 3.0，用Matlab軟件得到相應(yīng)的數(shù)值解，并在t-x平面上畫出x(t)的圖形。（3）分別對上述圖形講行分析，并解釋為什么會(huì)出現(xiàn)這些現(xiàn)象。 解：（1）根據(jù)Newton第二定律，F(xiàn)=ma，可得微分方程：m由題意知邊界值：x(0)=1,x(0)=0, 令y1=x，y2=dy1dt=y2dy可將二階方程轉(zhuǎn)化為：2=mg-k2y1-2hy2dty1(0)=1,y2(0)=0dxdtdxdt22=mg-kx-2h2dxdt ， 其中，m=1；g=9.8；k=2 當(dāng)h=0時(shí)，由matlab解得數(shù)值解：Columns

18、1 through 9 0 0.0000 0.0000 0.0000 0.0000 0.0001 0.0001 0.0002 0.00021.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000Columns 10 through 18 0.0004 0.0006 0.0009 0.0011 0.0022 0.0032 0.0043 0.0054 0.01081.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0001 1.0001 1.0003Columns 19 through 27 0

19、.0162 0.0216 0.0271 0.0541 0.0812 0.1083 0.1353 0.2155 0.29561.0008 1.0014 1.0021 1.0085 1.0191 1.0339 1.05281.1326 1.2461.Columns 100 through 108 8.2321 8.3593 8.4615 8.5636 8.6658 8.7679 8.87668.9852 9.09393.5019 3.2161 2.9505 2.6641 2.3688 2.0769 1.78341.5214 1.3034Columns 109 through 117 9.2025

20、9.3142 9.4260 9.5377 9.6494 9.7371 9.82479.9124 10.00001.1393 1.0345 1.0001 1.0383 1.1467 1.2773 1.44381.6412 1.8634 h=0時(shí)，x（t）圖形如下所示： 當(dāng)h=0.1時(shí)，由matlab解得數(shù)值解：Columns 1 through 90 0.0000 0.0000 0.0000 0.0000 0.0001 0.0001 0.0002 0.00021.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000Columns 10 through 18

21、0.0004 0.0006 0.0009 0.0011 0.0022 0.0054 0.01081.0000 1.0000 1.0000 1.0000 1.0000 1.0001 1.0003Columns 19 through 270.0162 0.0216 0.0271 0.0541 0.0812 0.2156 0.29581.0008 1.0014 1.0021 1.0085 1.01901.1308 1.2417.Columns 109 through 1178.5631 8.6690 8.7749 8.8809 8.98689.3934 9.52892.5857 2.4563 2.3

22、294 2.2106 2.10491.8880 1.8962Columns 118 through 1219.6467 9.7645 9.8822 10.00001.9353 2.0016 2.0909 2.1979 h=0.1時(shí)，x（t）圖形如下所示：1.0000 1.0000 0.0032 0.0043 1.0000 1.0001 0.1083 0.1353 1.0336 1.0523 9.1223 9.2579 1.9952 1.9213 （2）如果還受到外力干擾，則微分方程變?yōu)椋?mdx2dxdt2=mg-kx-2hdt+Bsinwt將其化為一階方程組：dy1dt=y2dy2=mg-k

23、2y1-2hy2+Bsinwt dty1(0)=1,y2(0)=0其中，m=1；g=9.8；k=2；B=1；h=0.1； 用Matlab解得數(shù)值解：w=1時(shí)Columns 1 through 90 0.0000 0.0000 0.0000 0.0002 0.00021.0000 1.0000 1.0000 1.0000 1.0000 1.0000.Columns 118 through 1219.8207 9.8804 9.9402 10.00001.8973 1.9210 1.9503 1.9846w=1.2時(shí)Columns 1 through 90 0.0000 0.0000 0.0000

24、 0.0002 0.00021.0000 1.0000 1.0000 1.0000 1.0000 1.0000.Columns 118 through 1219.8281 9.8854 9.9427 10.00001.6999 1.7568 1.8199 1.8886w=1.4時(shí)Columns 1 through 90 0.0000 0.0000 0.0000 0.0000 0.0001 1.0000 1.0000 0.0000 0.0001 1.0000 1.0000 0.0000 0.0001 0.0001 1.0000 0.0001 1.0000 0.00010.0002 0.00021

25、.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000.Columns 118 through 1259.6236 9.7309 9.8382 9.9454 9.9591 9.9727 9.9864 10.00002.2401 2.3171 2.4084 2.5103 2.5238 2.5373 2.5510 2.5647w=1.4時(shí)Columns 1 through 90 0.0000 0.0000 0.0000 0.0002 0.00021.0000 1.0000 1.0000 1.0000 1.0000 1.0000.C

26、olumns 118 through 1219.7823 9.8549 9.9274 10.0000 2.1150 2.0676 2.0285 1.9979w=1.6時(shí)Columns 1 through 90 0.0000 0.0000 0.0000 0.0002 0.00021.0000 1.0000 1.0000 1.0000 1.0000 1.0000.Columns 118 through 1219.6933 9.7955 9.8978 10.00000.8062 0.7503 0.7562 0.8231w=1.8時(shí)Columns 1 through 90 0.0000 0.0000

27、0.0000 0.0002 0.00021.0000 1.0000 1.0000 1.0000 1.0000 1.0000.Columns 118 through 1219.7752 9.8501 9.9251 10.0000 0.8629 1.0850 1.3383 1.61740.0000 0.0001 1.0000 1.0000 0.0000 0.0001 1.0000 1.0000 0.0000 0.0001 1.0000 1.0000 0.0001 1.0000 0.0001 1.0000 0.0001 1.0000w=2.0時(shí)Columns 1 through 90 0.0000

28、0.0000 0.0000 0.0000 0.0001 0.0001 0.0002 0.00021.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000.Columns 118 through 1219.7538 9.8359 9.9179 10.00002.3278 2.6023. w=3.0時(shí)Columns 1 through 90 0.00000.0002 0.00021.0000 1.00001.0000 1.0000.Columns 118 through 1219.7495 9.83302.2414 2.3303

29、2.8705 3.1243 0.0000 0.0000 0.0000 0.0001 0.0001 1.0000 1.0000 1.0000 1.0000 1.0000 9.9165 10.0000 2.4145 2.4915 如上圖所示，w分別為1.0；1.2；.；3.0時(shí)x（t）的圖形。 當(dāng)h=0時(shí)，即無阻尼振動(dòng)的情況下，得到t（x）圖形如下所示： （3）由上述圖形可以看出，不考慮外加作用力時(shí)，當(dāng)沒有摩擦等阻力的時(shí)候，簡諧運(yùn)動(dòng)一直進(jìn)行下去，當(dāng)有摩擦等阻力的時(shí)候，其振幅會(huì)逐漸的減小，最后趨近于零，但是周期不會(huì)發(fā)生變化。在外加一個(gè)作用力之后，整個(gè)運(yùn)動(dòng)情況就發(fā)生了很大的變化。首先剛開始的振幅和周期

30、都變化很大，后來隨著時(shí)間的推移，他們也逐漸開始趨于穩(wěn)定。但是，穩(wěn)定之后的振幅和周期都各不相同，與外力的頻率有直接的關(guān)系，這就是受迫振動(dòng)，不僅跟系統(tǒng)本身的性質(zhì)有關(guān)，還和外界干擾的情況有很大關(guān)系?；謴?fù)系數(shù)k與w的值越接近時(shí)，隨著時(shí)間增加，物體的振動(dòng)振幅會(huì)一直增大，當(dāng)w=k=2時(shí)，物體的振幅增大的速度最快，可以預(yù)見，在t趨于無窮大時(shí)，振幅也趨于無窮大，這就是共振現(xiàn)象。也就是當(dāng)外力與系統(tǒng)處于共振時(shí)，會(huì)引起振幅無限增大的振動(dòng)，這在機(jī)械和建筑中一般是必須嚴(yán)格避免的。 6、人口模型與曲線擬合問題表2.1列出的是美國1790一1980年的人口統(tǒng)計(jì)表.（1）試用Malthus人口模型，按三段時(shí)間(1790-18

31、50,1850-1910,1910-1970)分別確定其增長率r。并將數(shù)據(jù)和不同r值的三段曲線畫在同一張圖上.（2）利用Logistic模型，重新確定固有增長率r和最大容量Nm，作圖，再利用該模型得到的結(jié)果預(yù)測1990年的美國人口數(shù)。解：（1）分段研究，我們先求出增長率，編寫命令如下：t = 10 : 10 : 70;x = 3.9 5.3 7.2 9.6 12.9 17.1 23.2;p = polyfit(t, log(x), 1); p得到結(jié)果p = 0.0296即17901850年時(shí)間段的增長率是0.0296同理可以得到另外兩段的增長率，1850-1910年 為0.0226191019

32、70年為0.0129每次畫出圖形，使用以下命令plot(t, x, bx);hold on在以上每一次求時(shí)順便畫出圖形，得到最后的總圖如下 (2)利用Logistics模型，將以上所有數(shù)據(jù)進(jìn)行擬合，編寫程序如下t = 0 : 10 : 190;x = 3.9 5.3 7.2 9.6 12.9 17.1 23.2 31.4 38.6 50.2 62.9 76.0 92.0 106.5 123.2 131.7 150.7 179.3 204.0 226.5;f = inline(p(1)./(1+p(2)*exp(-p(3).*t),p,t);p = lsqcurvefit(f, 300, 50,

33、 0.02, t, x); 得到結(jié)果人口最大容量Nm是360.3560（百萬），增長率r是0.0234 如下圖所示： 輸入t_pre = 200;x_pre = p(1) ./ (1 + p(2) * exp(-p(3)*t_pre) );可得最后的結(jié)果，即1990年的人口總數(shù)x_pre =241.7704萬人。如下圖： 7、加分實(shí)驗(yàn)（餐廳廢物的堆肥優(yōu)化問題） 一家環(huán)保餐廳用微生物將剩余的食物變成肥料。餐廳每天將剩余的食物制成槳狀物并與蔬菜下腳及少量紙片混合成原料，加入真菌菌種后放入容器內(nèi)。真菌消化這此混合原料，變成肥料，由于原料充足，肥料需求旺盛，餐廳希望增加肥料產(chǎn)量。由于無力購置新設(shè)備，餐

34、廳希望用增加真菌活力的辦法來加速肥料生產(chǎn).試通過分析以前肥料生產(chǎn)的記錄(如表2.2所示)，建立反映肥料生成機(jī)理的數(shù)學(xué)模型，提出改善肥料生產(chǎn)的建議。 解答：1.問題條件假設(shè)1）每次堆肥的質(zhì)量相等2）所給的幾次堆肥混合物的比例僅由當(dāng)天的實(shí)際情況決定3）所有分離肥倉工作條件相同4）每磅蔬菜可提供的氧氣量相等5）細(xì)菌消耗的溶解氧氣完全由蔬菜葉提供6）每天提供的廢物混合物中的化學(xué)成分大致相同7）廢物混合物在喂給細(xì)菌前混合均勻并保持良好的通氣環(huán)境。2.問題分析堆肥是利用微生物的分解作用將有機(jī)廢物轉(zhuǎn)化成無害穩(wěn)定形式的生物化學(xué)過程，要提高堆肥產(chǎn)量方法之一就是通過增強(qiáng)細(xì)菌的生長繁殖能力以提高分解率。 細(xì)菌群體的

35、增長一般要經(jīng)歷延滯期，加速生長期、對數(shù)生長期、減數(shù)生長期、穩(wěn)定期、加速死亡期和對數(shù)死亡期。另外，對數(shù)生長期培養(yǎng)基中所有養(yǎng)分都過剩，細(xì)菌可以充分繁殖，其倍增速率恒定，取決于底物濃度、溫度、水活度、供氧量。對于當(dāng)前該餐廳來說底物濃度由每天的剩余食物、蔬菜葉和碎紙屑決定。碎紙屑是吸收水分的調(diào)理劑，微菌呼吸所需要的溶解氧由蔬菜也提供，水活度可以通過測定相對濕度來決定，其關(guān)系式是：相對濕度：B=P/P0*100%水活度： aw=P/P0 其中P為該溶液蒸汽壓，P0為純水蒸汽壓， 在這個(gè)堆肥系統(tǒng)中，可供微菌消耗的營養(yǎng)底物和溶解氧都是有限的，他們的消耗會(huì)對微菌生長率產(chǎn)生重要影響，一種微菌消耗營養(yǎng)底物的速率和

36、底物濃度之間的關(guān)系式為：ds/dt=KmSX/（Ks+S），式中ds/dt表示為底物的有效消耗率，x表示為微菌濃度，Km表示為最大有效系數(shù)，在高濃度營養(yǎng)底物中最大的物料消耗率（物質(zhì)質(zhì)量/微菌一天的質(zhì)量）；Ks表示為半速系數(shù)（質(zhì)量/體積），S表示為有限底物的濃度（質(zhì)量/體積）微菌生長過程是一個(gè)生物化學(xué)過程，其到經(jīng)驗(yàn)最佳比例，由于假設(shè)各次堆肥后肥料質(zhì)量相同，因此堆肥時(shí)間較短就對應(yīng)了較好的廢料配比，將12組數(shù)據(jù)按其堆肥日期及完全堆肥時(shí)間分成三組，每組模型二：營養(yǎng)底物和氧氣是細(xì)菌生長的兩種底物，物耗公式為：dSi/dt=KmiXSi/(Ksi+Si) （1）i=1時(shí)代表營養(yǎng)底物中可降解得有機(jī)物。i=2時(shí)代表氧氣。1）在高濃度底物中，物料的轉(zhuǎn)化過程很迅速，進(jìn)一步增加底物濃度就不再引起底物轉(zhuǎn)化率的提高，可假設(shè)S1>>Ks1,S2>>Ks2，