數(shù)學建模-微分方程的穩(wěn)定性

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

17、QbVrVdQcbcrccVQdtdc10)(bVQdcrdtdb)(2設設V=1.6104 m3 和和 3104 m3, 用數(shù)值方法解微分方程組用數(shù)值方法解微分方程組:有害物質(zhì)濃度將有有害物質(zhì)濃度將有約約1300小時超過小時超過2c*, 最高達到最高達到5c*單池模型的動態(tài)過程單池模型的動態(tài)過程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)有害物質(zhì)濃度將有有害物質(zhì)濃度將有約約

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

19、穩(wěn)態(tài)狀況11110121111)(,:crVccQbrVdVQcPdcrQVcc02101當當 時時P1穩(wěn)定穩(wěn)定c0=c01, V18103 m3池池方程方程的平衡點的平衡點)(4)(21,)(:21222111221112222122122VQdcrVQdbrcrVQdbrcrrccrVccQbP*11*2*1*12*2)()(cbrdcrccccQVcc要求穩(wěn)態(tài)下要求穩(wěn)態(tài)下雙池模型的穩(wěn)態(tài)狀況雙池模型的穩(wěn)態(tài)狀況*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相配合的方案相配