1、穩定性模型穩定性模型1 微分方程穩定性的理論知微分方程穩定性的理論知識識2 捕魚業的持續收獲捕魚業的持續收獲3 軍備競賽軍備競賽4 廢水的生物處理廢水的生物處理 穩定性模型穩定性模型 對象仍是動態過程，而建模目的是研究時對象仍是動態過程，而建模目的是研究時間充分長以后過程的變化趨勢間充分長以后過程的變化趨勢 平衡狀平衡狀態是否穩定態是否穩定. 不求解微分方程，而是用微分方程穩定性不求解微分方程，而是用微分方程穩定性理論研究平衡狀態的穩定性理論研究平衡狀態的穩定性.1 微分方程穩定性的理論知識微分方程穩定性的理論知識1.1 一階微分方程的平衡點及其穩定性一階微分方程的平衡點及其穩定性1.2 線性

2、常系數微分方程組的平衡點及其穩定性線性常系數微分方程組的平衡點及其穩定性1.1 一階微分方程的平衡點及其穩定性一階微分方程的平衡點及其穩定性) 1 ()(xFx 一階非線性自治一階非線性自治(右端不含右端不含t)方程方程F(x)=0的根的根x0 微分方程的微分方程的平衡點平衡點000 xxxxx設設x(t)是方程的解，若從是方程的解，若從x0 某鄰域的任一初值出發，某鄰域的任一初值出發，都有都有,)(lim0 xtxt稱稱x0是方程是方程(1)的的穩定平衡點穩定平衡點.不求不求x(t), 判斷判斷x0穩定性的方法穩定性的方法直接法直接法)2()(00 xxxFx(1)的近似線性方程的近似線性方

3、程)1 (),2(0)(00對穩定xxF)1 (),2(0)(00對不穩定xxF1.2 線性常系線性常系數微分方程組數微分方程組dycxtybyaxtx)()(的平衡點及其穩定性的平衡點及其穩定性平衡點平衡點P0(x0,y0)=(0,0) 代數方程代數方程00dycxbyax的根的根若從若從P0某鄰域的任一初值出發，都有某鄰域的任一初值出發，都有,)(lim0 xtxt稱稱P0是微分方程的是微分方程的穩定平衡點穩定平衡點,)(lim0ytyt記系數矩陣記系數矩陣dcbaA特征方程特征方程0)det( IAAqdapqpdet)(02特征根特征根2/ )4(22, 1qpp線性常系數線性常系數微

4、分方程組微分方程組dycxtybyaxtx)()(的平衡點及其穩定性的平衡點及其穩定性特征根特征根2/ )4(22, 1qpp平衡點平衡點 P0(0,0)微分方程一般解形式微分方程一般解形式ttecec2121平衡點平衡點 P0(0,0)穩定穩定平衡點平衡點 P0(0,0)不穩定不穩定 1,2為負數或有負實部為負數或有負實部p 0 且且 q 0p 0 或或 q 0cErEpNEESETER)1 ()()()(R(E)=0時的捕撈強度時的捕撈強度(臨界強度臨界強度) Es=2ER)1 (rENxsspc臨界強度下的漁場魚量臨界強度下的漁場魚量 cp ,ER)1(2pNcrERE*令令=0)1 (

5、pNcrEsssxE,xs由成本由成本價格比決定價格比決定捕撈過度捕撈過度捕撈捕撈過度過度T(E)0rES(E)Es2Es1S(E)pNEE*pNE/2)()()(ESETERpNcpN2/)1 ()(rEpNEET收入收入cEES)(支出支出利潤利潤臨界強度臨界強度Es=0)/2/(NcpNc 經濟學捕撈過度經濟學捕撈過度*1EEEss2/pNc )/2(Ncp*2EEEss 生態學捕撈過度生態學捕撈過度捕魚業的捕魚業的持續收獲持續收獲在自然增長和捕撈情況的合理假設下建模在自然增長和捕撈情況的合理假設下建模.用平衡點穩定性分析確定漁場魚量穩定條件用平衡點穩定性分析確定漁場魚量穩定條件,討論產

6、量、效益和捕撈過度討論產量、效益和捕撈過度3個模型個模型.3 軍備競賽軍備競賽 描述雙方描述雙方(國家或國家集團國家或國家集團)軍備競賽過程軍備競賽過程. 解釋解釋(預測預測)雙方軍備競賽的結局雙方軍備競賽的結局.假設假設 1）由于相互不信任，一方軍備越大，另一）由于相互不信任，一方軍備越大，另一方軍備增加越快；方軍備增加越快； 2）由于經濟實力限制，一方軍備越大，對）由于經濟實力限制，一方軍備越大，對自己軍備增長的制約越大；自己軍備增長的制約越大； 3）由于相互敵視或領土爭端，每一方都存）由于相互敵視或領土爭端，每一方都存在增加軍備的潛力在增加軍備的潛力.進一步進一步假設假設 1）2）的作用

7、為線性；）的作用為線性；3）的作用為常數）的作用為常數.目的目的gkyxtx)( 建模建模軍備競賽的結局軍備競賽的結局微分方程的平衡點及其穩定性微分方程的平衡點及其穩定性x(t)甲方軍備數量，甲方軍備數量， y(t)乙方軍備數量乙方軍備數量hylxty)( , 本方經濟實力的制約；本方經濟實力的制約； k, l 對方對方軍備數量的刺激；軍備數量的刺激；g, h 本方本方軍備競賽的潛力軍備競賽的潛力.t 時的時的x(t)，y(t)klAqpdet0)(klhglyklgkhx00,平衡點平衡點穩定性判斷穩定性判斷lkA系數系數矩陣矩陣平衡點平衡點(x0, y0)穩定的條件穩定的條件0, 0qpk

8、lhylxtygkyxtx)()(模型模型軍備競賽軍備競賽模型的定性解釋模型的定性解釋kl雙方軍備穩定雙方軍備穩定(時間充分時間充分長后趨向有限值長后趨向有限值)的條件的條件1) 雙方經濟制約大于雙方軍備刺激時，軍備競賽雙方經濟制約大于雙方軍備刺激時，軍備競賽 才會穩定，否則軍備將無限擴張才會穩定，否則軍備將無限擴張.平衡點平衡點klhglyklgkhx00,2) 若若g=h=0, 則則 x0=y0=0, 在在 kl 下下 x(t), y(t)0, 即友好鄰國通過裁軍可達到永久和平即友好鄰國通過裁軍可達到永久和平.hylxtygkyxtx)()(模型模型 , 本方經濟實力的制約；本方經濟實力的

9、制約； k, l 對方對方軍備數量的刺激；軍備數量的刺激；g, h 本方本方軍備競賽的潛力軍備競賽的潛力.3）若）若 g,h 不為零，即便雙方一時和解，使某時不為零，即便雙方一時和解，使某時x(t), y(t)很小，但因很小，但因 ，也會重整軍備，也會重整軍備.0, 0yx4）即使某時一方）即使某時一方(由于戰敗或協議由于戰敗或協議)軍備大減軍備大減, 如如 x(t)=0， 也會因也會因 使該方重整軍備，使該方重整軍備，gkyx 即存在互不信任即存在互不信任( ) 或固有爭端或固有爭端( ) 的單方面的單方面裁軍不會持久裁軍不會持久.0k0g模型的定性解釋模型的定性解釋 , 本方經濟實力的制約

10、；本方經濟實力的制約； k, l 對方對方軍備數量的刺激；軍備數量的刺激；g, h 本方本方軍備競賽的潛力軍備競賽的潛力.hylxtygkyxtx)()(模型模型4 廢水的生物處理廢水的生物處理 廢水處理廢水處理 (去掉有害的有機物去掉有害的有機物) 通常通常有生物化學與物理化學兩種方法有生物化學與物理化學兩種方法.背景與問題背景與問題 生物處理生物處理 利用微生物利用微生物(主要是細菌主要是細菌)的生命活動過的生命活動過程程, 把廢水中的有機物轉化為簡單的無機物把廢水中的有機物轉化為簡單的無機物. 已知廢水中有害物質濃度為已知廢水中有害物質濃度為10-310-2g/m3, 要將濃度降要將濃度

11、降至至510-4g/m3以下以下, 需建立廢水與微生物混合的處理池需建立廢水與微生物混合的處理池. 設廢水將以設廢水將以10m3/h的流量進入處理池的流量進入處理池, 確定處理池的確定處理池的容積容積, 使排出廢水中使排出廢水中有害物質的濃度達到規定的標準有害物質的濃度達到規定的標準.模型假設模型假設 生物化學提供了有機物分解、轉化和生物化學提供了有機物分解、轉化和微生物增殖、衰亡的規律及相關參數微生物增殖、衰亡的規律及相關參數2. 微生物依于有害物質分解、轉化的能量而增殖的速微生物依于有害物質分解、轉化的能量而增殖的速率與有害物質濃度成正比，比例系數率與有害物質濃度成正比，比例系數r2=1.

12、26m3/g . h4. 處理池內有害物質和微生物任何時候都均勻混合處理池內有害物質和微生物任何時候都均勻混合,排出廢水中有害物質和微生物的濃度與池內相同排出廢水中有害物質和微生物的濃度與池內相同.3. 微生物的自然死亡率為常數微生物的自然死亡率為常數 d=10-5/h1. 有害物質被微生物分解、轉化而消失的速率與微生有害物質被微生物分解、轉化而消失的速率與微生物濃度成正比，比例系數物濃度成正比，比例系數r1=0.1m3/g . hc(t) 時刻時刻 t 有害物質的濃度有害物質的濃度 b(t) 時刻時刻 t 微生物的濃度微生物的濃度模型假設模型假設 生物化學提供了有機物分解、轉化和生物化學提供

13、了有機物分解、轉化和微生物增殖、衰亡的規律及相關參數微生物增殖、衰亡的規律及相關參數6. 進入處理池的廢水中有害物質濃度為進入處理池的廢水中有害物質濃度為c0, c01 c0 c02, c01= 10-3g/m3, c02= 10-2g/m3, c0可以改變可以改變, 最壞情況是最壞情況是c0由由c01突然增加到突然增加到c027. 環境保護法規定的廢水中有害物質濃度為環境保護法規定的廢水中有害物質濃度為c*=510-4 g/m3, 它是長期穩定排放的標準它是長期穩定排放的標準, 如果是短期排放并超如果是短期排放并超標不大標不大, 可以用處罰等方法解決可以用處罰等方法解決.5. 忽略蒸發等因素

14、忽略蒸發等因素, 廢水進入處理池和排出處理池的廢水進入處理池和排出處理池的流量均為常數流量均為常數Q=10m3/h; 廢水滿池廢水滿池, 池的容積為池的容積為V單池模型單池模型建立一個處理池建立一個處理池 (t, t+ t) 內內池內有害物質的平衡池內有害物質的平衡改變量改變量 = 進入量進入量 排出量排出量 分解轉化量分解轉化量)()(tcttcVttQctQc)(0ttctbVr)()(1c(t) 有害物質濃度有害物質濃度b(t) 微生物的濃度微生物的濃度V池的容積池的容積bcrccVQdtdc10)(t, t+ t) 內內池內微生物的平衡池內微生物的平衡)()(tbttbVttQbttV

15、dbttctbVr)()()()(2bVQdcrdtdb)(2Q流量流量非線性方程組無解析解非線性方程組無解析解單池模型的穩態狀況單池模型的穩態狀況平衡點平衡點用微分方程穩定性理論可以驗證用微分方程穩定性理論可以驗證:bcrccVQdtdc10)(bVQdcrdtdb)(2cVrccQbVrVdQcP1021)(,:0,:02bccPVQdcr/02微生物的增殖率大于死亡和排除率微生物的增殖率大于死亡和排除率dcrQVcc020當當 時時P1穩定穩定, P2不穩定不穩定單池模型的穩態狀況單池模型的穩態狀況c01 c0 c02Q=10m3/hr2=1.26m3/g . hd=10-5/hc0=c

16、01= 10-3g/m3V8103 m3V1.6104 m3c0=c*= 510-4 g/m3Vc0dcrQV02平衡點平衡點P1穩定條件穩定條件為使穩定狀況下有害物質濃度達到規定標準為使穩定狀況下有害物質濃度達到規定標準c*, 處理池的容積至少需要達到處理池的容積至少需要達到1.6104 m3.一個長寬各一個長寬各100 m, 深深1.6m的的池子池子!考察最壞情況考察最壞情況取穩定平衡點取穩定平衡點P1為初值為初值, 即即單池模型的動態過程單池模型的動態過程當當c0=c01時池內濃度已處于穩態時池內濃度已處于穩態, c0突然增加到突然增加到c02)0()0()0(,)0(102cVrccQ

17、bVrVdQcbcrccVQdtdc10)(bVQdcrdtdb)(2設設V=1.6104 m3 和和 3104 m3, 用數值方法解微分方程組用數值方法解微分方程組:有害物質濃度將有有害物質濃度將有約約1300小時超過小時超過2c*, 最高達到最高達到5c*單池模型的動態過程單池模型的動態過程01000200030004000500000.511.522.53x 10-3c*V=1.6104m3c(10-3g/m3)t(h)01000200030004000500000.511.522.53x 10-3c*V=3104m3c( 10-3g/m3)t(h)有害物質濃度將有有害物質濃度將有約約9

18、00小時超過小時超過2c*, 最高達到最高達到3c*13002c*要達到規定的標準需要太大的池子要達到規定的標準需要太大的池子!長寬各長寬各100 m, 深深3m的的池子池子兩個串接的池子兩個串接的池子雙池模型雙池模型V1, c1, b1Q, c0V2, c2, b2Q, c1, b1池池池池1111011)(cbrccVQdtdc1222222)(bVQbVQdcrdtdb2212122)(cbrccVQdtdc11121)(bVQdcrdtdb與單池模型相同與單池模型相同(只是加上下標只是加上下標1)增加從池增加從池的流入量的流入量池池方程方程的平衡點的平衡點雙池模型的穩態狀況雙池模型的穩

19、態狀況11110121111)(,:crVccQbrVdVQcPdcrQVcc02101當當 時時P1穩定穩定c0=c01, V18103 m3池池方程方程的平衡點的平衡點)(4)(21,)(:21222111221112222122122VQdcrVQdbrcrVQdbrcrrccrVccQbP*11*2*1*12*2)()(cbrdcrccccQVcc要求穩態下要求穩態下雙池模型的穩態狀況雙池模型的穩態狀況*11*2*1*12)()(cbrdcrccccQV在在c0=c01, V18103 m3下取值計算下取值計算c1, b1, V2V1(103m3)c1(10-3g/m3)b1(10-3g/m3)V2(103m3)81.000.0216.03100.802.509.63120.674.145.39140.575.312.37160.506.180.10,)(,1111012111crVccQbrVdVQc應選擇較大的應選擇較大的V1和較小的和較小的V2相配合的方案相配