第2講房室模型

1、1第二講第二講 房室模型房室模型一、藥物在體內(nèi)的分布與排除（二室模型）一、藥物在體內(nèi)的分布與排除（二室模型）二、三種群二、三種群Volterra 模型（三室模型）模型（三室模型）三、三、SARS模型模型2房室模型房室模型 藥物動力學(xué)通常用房室模擬人體，只要體內(nèi)某些部位接受藥物動力學(xué)通常用房室模擬人體，只要體內(nèi)某些部位接受或消除藥物的速率相似，即可歸入一個房室。房室模型僅是進或消除藥物的速率相似，即可歸入一個房室。房室模型僅是進行藥動學(xué)分析的一種抽象概念，并不一定代表某一特定解剖部行藥動學(xué)分析的一種抽象概念，并不一定代表某一特定解剖部位。位。 把機體劃分為一個或多個獨立單元，可對藥物在體內(nèi)吸收、

2、把機體劃分為一個或多個獨立單元，可對藥物在體內(nèi)吸收、分布、消除的特性作出模式圖，以建立數(shù)學(xué)模型，揭示其動態(tài)分布、消除的特性作出模式圖，以建立數(shù)學(xué)模型，揭示其動態(tài)變化規(guī)律。變化規(guī)律。 1. 假設(shè)機體給藥后，藥物立即在全身各部位達到動態(tài)平衡，假設(shè)機體給藥后，藥物立即在全身各部位達到動態(tài)平衡，這時把整個機體視為一個房室，稱為這時把整個機體視為一個房室，稱為一室模型一室模型。 2. 假設(shè)藥物進入機體后，瞬時就可在血液供應(yīng)豐富的組織假設(shè)藥物進入機體后，瞬時就可在血液供應(yīng)豐富的組織（如血液、肝、腎等）分布達到動態(tài)平衡，（如血液、肝、腎等）分布達到動態(tài)平衡， 然后再在血液供應(yīng)然后再在血液供應(yīng)較少或血流較慢的

3、組織（如脂肪、皮膚、骨骼等）分布達到動較少或血流較慢的組織（如脂肪、皮膚、骨骼等）分布達到動態(tài)平衡，態(tài)平衡， 此時可把這些組織分別稱為此時可把這些組織分別稱為中央室中央室和和周邊室周邊室，即，即二室二室模型模型。3infectiverecoveredsusceptiblekl (1) (2)( )( )( )1 (3)0, (0)0odiksilidtdrlidts ti tr tni( )i r 將人群劃分為三類（見右圖）：易感染者、已感染將人群劃分為三類（見右圖）：易感染者、已感染者和已恢復(fù)者（者和已恢復(fù)者（recovered）。分別記）。分別記 t 時刻的三類人數(shù)時刻的三類人數(shù)為為s(t

4、)、i(t)和和r(t)，則可建立下面的，則可建立下面的三房室模型三房室模型： 三房室模型三房室模型SIRSIR模型：模型：單房室模型單房室模型MalthusMalthus模型與模型與LogisticLogistic模型模型二房室模型二房室模型捕食（捕食（P-PP-P）模型）模型4一、藥物在體內(nèi)的分布與排除（二室模型） 藥物進入機體形成藥物進入機體形成血藥濃度血藥濃度(單位體積血液的藥物量單位體積血液的藥物量) 血藥濃度需保持在一定范圍內(nèi)血藥濃度需保持在一定范圍內(nèi)給藥方案設(shè)計給藥方案設(shè)計 藥物在體內(nèi)吸收、分布和排除過程藥物在體內(nèi)吸收、分布和排除過程 藥物動力學(xué)藥物動力學(xué) 建立建立房室模型房室模

5、型藥物動力學(xué)的基本步驟藥物動力學(xué)的基本步驟 房室房室機體的一部分，藥物在一個房室內(nèi)均勻分布機體的一部分，藥物在一個房室內(nèi)均勻分布(血藥血藥濃度為常數(shù)濃度為常數(shù))，在房室間按一定規(guī)律轉(zhuǎn)移，在房室間按一定規(guī)律轉(zhuǎn)移 本節(jié)討論本節(jié)討論二室模型二室模型中心室中心室(心、肺、腎等心、肺、腎等)和和周邊室周邊室(四四肢、肌肉等肢、肌肉等)5 中心室中心室周邊室周邊室給藥給藥排除排除)(0tf111)(),(Vtxtc222)(),(Vtxtc12k21k13k)()(02211131121tfxkxkxktx模型假設(shè)模型假設(shè) 中心室中心室(1)和周邊室和周邊室(2), ,容積不變?nèi)莘e不變 藥物在房室間藥物在

6、房室間轉(zhuǎn)移速率轉(zhuǎn)移速率及向體外及向體外排除速率排除速率，與該室血藥濃度成正比與該室血藥濃度成正比 藥物從體外進入中心室，在二室間藥物從體外進入中心室，在二室間相互轉(zhuǎn)移相互轉(zhuǎn)移, ,從中心室排出體外從中心室排出體外模型建立模型建立2 , 1)()(iVtctxiii容積濃度藥量給藥速率0f2211122)(xkxktx6tttteBeAtceBeAtc222111)()(1321132112kkkkk2211122121022112113121)()()()(ckckVVtcVtfckVVckktc2 , 1),()(itcVtxiii線性常系數(shù)線性常系數(shù)非齊次方程非齊次方程對應(yīng)齊次對應(yīng)齊次方程

7、通解方程通解模型建立模型建立7)()()()()()()(212022121101tttteeVkDtcekekVDtc0)0(,)0(,0)(21010cVDctf幾種常見的給藥方式幾種常見的給藥方式1. .快速靜脈注射快速靜脈注射t=0 瞬時瞬時注射劑量注射劑量D0的藥物進入中心室的藥物進入中心室, ,血血藥濃度立即為藥濃度立即為D0/V12211122121022112113121)()()()(ckckVVtcVtfckVVckktc1321132112kkkkk給藥速率給藥速率 f0(t) 和初始條件和初始條件初始條件：初始條件：812211312121221131212213210

8、122221130111)(,)(0,)(0,)(BVkkkVBAVkkkVATtVkkkkeBeAtcTtVkkeBeAtctttt0)0(, 0)0(,)(2100ccktf2. .恒速靜脈滴注恒速靜脈滴注2211122121022112113121)()()()(ckckVVtcVtfckVVckktct T, c1(t)和和 c2(t)按指數(shù)規(guī)律趨于零按指數(shù)規(guī)律趨于零藥物以速率k0進入中心室0Tt 90010 xkf )(0tx吸收室中心室000010)0()(Dxxktx tkttEeBeAetc01)(1tkeDtx0100)(tkekDtxktf010100010)()(3. .

9、口服或肌肉注射口服或肌肉注射相當(dāng)于藥物相當(dāng)于藥物( 劑量劑量D0)先進入吸收室，吸收后進入中心室先進入吸收室，吸收后進入中心室吸收室藥量吸收室藥量x0(t)2211122121022112113121)()()()(ckckVVtcVtfckVVckktcEBAcc,0)0(, 0)0(2110ttBeAetctc)()(11先估計參數(shù)：先估計參數(shù)：A, B, , 由各種給藥方式下的由各種給藥方式下的 c1(t), c2(t) 確定參數(shù)確定參數(shù)k12, k21, k13, V1,V2t=0快速靜脈注射快速靜脈注射D0 , ,在在ti (i=1,2,n)測得測得c1(ti)()()()(2121

10、101ttekekVDtc充分大設(shè)t ,由較大的由較大的 用最小二乘法定用最小二乘法定A A, , )(,1iitct由較小的由較小的 用最小二乘法定用最小二乘法定B, , )(,1iitctttAeeVkDtc)()()(12101反問題：反問題：11211312kkkBAVDc101)0(011130)(dttcVkD0,21cct1321132112kkkkkBAVkD1130ABBAk)(131321kk再估計參數(shù)再估計參數(shù): 進入中心室的藥物全部排除進入中心室的藥物全部排除132112,kkk12注：注： 建立房室模型的目的是研究體內(nèi)血藥濃度的建立房室模型的目的是研究體內(nèi)血藥濃度的變

11、化過程，確定諸如轉(zhuǎn)移和排除速率系數(shù)等參數(shù)，變化過程，確定諸如轉(zhuǎn)移和排除速率系數(shù)等參數(shù)，為制定給藥方案和計量大小提供依據(jù)。建模過程是為制定給藥方案和計量大小提供依據(jù)。建模過程是將機理分析和測試分析相結(jié)合，先由機理分析確定將機理分析和測試分析相結(jié)合，先由機理分析確定方程的形式，再由測試數(shù)據(jù)估計參數(shù)。方程的形式，再由測試數(shù)據(jù)估計參數(shù)。 可根據(jù)需要選用一室模型、二室模型或多室可根據(jù)需要選用一室模型、二室模型或多室模型，甚至非線性房室模型。模型，甚至非線性房室模型。 13二、三種群Volterra 模型 自然環(huán)境中的某一種生物的群體自然環(huán)境中的某一種生物的群體, 生態(tài)學(xué)上稱為種群。生態(tài)學(xué)上稱為種群。如果

12、一個自然環(huán)境中有兩個或兩個以上種群生存如果一個自然環(huán)境中有兩個或兩個以上種群生存, 那么它們那么它們之間就要存在著或是相互競爭之間就要存在著或是相互競爭, 或是相互依存或是相互依存, 或是弱肉強或是弱肉強食食( 食餌與捕食者食餌與捕食者) 的關(guān)系的關(guān)系, 自然界中不同種群之間還存在自然界中不同種群之間還存在著一種非常有趣的既有依存、又有制約的生存方式著一種非常有趣的既有依存、又有制約的生存方式: 種群甲種群甲靠豐富的自然資源生長靠豐富的自然資源生長, 而而種群乙種群乙靠捕食種群甲為生靠捕食種群甲為生, 種群種群丙丙又靠捕食種群乙為生又靠捕食種群乙為生, 類似的現(xiàn)象還存在很多。類似的現(xiàn)象還存在很

13、多。 假設(shè)假設(shè): 一個島嶼上棲居著食肉爬行動物和哺乳動物一個島嶼上棲居著食肉爬行動物和哺乳動物,又又長著茂盛的植物長著茂盛的植物, 爬行動物以哺乳動物為食爬行動物以哺乳動物為食, 哺乳動物又依哺乳動物又依賴植物生存賴植物生存, 由此建立描述三種群數(shù)量變化規(guī)律的微分方程由此建立描述三種群數(shù)量變化規(guī)律的微分方程模型模型三房室模型三房室模型。141. 模型的建立模型的建立 當(dāng)植物、哺乳動物、爬行動物在一個自然環(huán)境中生存時當(dāng)植物、哺乳動物、爬行動物在一個自然環(huán)境中生存時, 把把植物、哺乳動物、爬行動物植物、哺乳動物、爬行動物的數(shù)量分別記作的數(shù)量分別記作x1(t), x2(t), x3(t) 。若不考慮

14、自然資源對植物的限制。若不考慮自然資源對植物的限制, 植物獨立生存時以指植物獨立生存時以指數(shù)規(guī)律增長數(shù)規(guī)律增長, 相對增長率為相對增長率為r1, 即即x(t)=r1x1, 而哺乳動物的存在而哺乳動物的存在使植物的增長率減小使植物的增長率減小, 設(shè)減小的程度與捕食者數(shù)量成正比設(shè)減小的程度與捕食者數(shù)量成正比, 于于是是植物的模型植物的模型為為:)(21111xrxdtdx比例系數(shù)比例系數(shù)1 反映哺乳動物掠取植物的能力。反映哺乳動物掠取植物的能力。(1)15 哺乳動物離開植物無法生存哺乳動物離開植物無法生存, 設(shè)它獨自存在時死亡率設(shè)它獨自存在時死亡率為為r2, 即即 x2(t) = -r2 x2,

15、而植物的存在又為哺乳動物提供了食而植物的存在又為哺乳動物提供了食物物, 植物的存在相當(dāng)于使哺乳動物的死亡率降低植物的存在相當(dāng)于使哺乳動物的死亡率降低, 且促使哺且促使哺乳動物增長乳動物增長, 設(shè)這種作用與植物的數(shù)量成正比設(shè)這種作用與植物的數(shù)量成正比, 則有則有:)(12222xrxdtdx比例系數(shù)比例系數(shù)2 反映植物對哺乳動物的供養(yǎng)能力。反映植物對哺乳動物的供養(yǎng)能力。(2)16 哺乳動物又為爬行動物提供了食物哺乳動物又為爬行動物提供了食物, , 爬行動物的存爬行動物的存在使哺乳動物的增長率減小在使哺乳動物的增長率減小, , 設(shè)減小的程度與爬行動物設(shè)減小的程度與爬行動物的數(shù)量成正比的數(shù)量成正比,

16、 , 于是于是(2) (2) 式右端應(yīng)減去爬行動物對哺乳式右端應(yīng)減去爬行動物對哺乳動物增長的阻滯作用動物增長的阻滯作用, , 于是于是哺乳動物的模型哺乳動物的模型應(yīng)為應(yīng)為: :)(312222xxrxdtdx(3)比例系數(shù)比例系數(shù)反映爬行動物掠取哺乳動物的能力。反映爬行動物掠取哺乳動物的能力。17 爬行動物離開哺乳動物無法生存爬行動物離開哺乳動物無法生存, 設(shè)它獨自存在時死設(shè)它獨自存在時死亡率為亡率為r3, 即即x3(t)= - r3 x3, 而哺乳動物的存在又為爬行動物而哺乳動物的存在又為爬行動物提供了食物提供了食物, 相當(dāng)于使爬行動物的死亡率降低相當(dāng)于使爬行動物的死亡率降低, 且促使爬且促

17、使爬行動物的增長行動物的增長, 于是于是爬行動物的模型爬行動物的模型為為:比例系數(shù)比例系數(shù)3 反映哺乳動物對爬行動物的供養(yǎng)能力。反映哺乳動物對爬行動物的供養(yǎng)能力。)(23333xrxdtdx(4)18 方程方程(1) 、(3) 、(4) 構(gòu)成植物、哺乳動物、爬行動物構(gòu)成植物、哺乳動物、爬行動物三者依存、制約現(xiàn)象的數(shù)學(xué)模型三者依存、制約現(xiàn)象的數(shù)學(xué)模型, 即即)(21111xrxdtdx)(312222xxrxdtdx(5)(23333xrxdtdx記植物、哺乳動物、爬行動物的初始數(shù)量分別為記植物、哺乳動物、爬行動物的初始數(shù)量分別為:x1(0)=x10 , x2(0)=x20 , x3(0)=x3

18、0 (6)19 微分方程組微分方程組(5) 沒有解析解沒有解析解, 可利用可利用MatLab 求微分求微分方程組方程組(5) 的數(shù)值解的數(shù)值解, 通過對數(shù)值結(jié)果和圖形的觀察通過對數(shù)值結(jié)果和圖形的觀察, 猜測猜測它的解析解的構(gòu)造它的解析解的構(gòu)造; 為求微分方程組為求微分方程組(5) 及初始條件及初始條件( 6) 的數(shù)值解的數(shù)值解x1(t), x2(t), x3(t) ( 并作圖并作圖) , 設(shè)設(shè) r1=1, r2=0.5, r3=0.6, 1=0.1,2 = 0.02, 3= 0.06, = 0.1, x10= 100, x20= 40, x30= 6, 用用MatLab 軟件編制程序如下軟件編

19、制程序如下:2. 模型的求解模型的求解20function f=fun1(t,x); r1 =1; r2 =0.5; r3 =0.6; lambda1 =0.1; lambda2 =0.02; lambda3=0.06; mu=0.1; f= x(1)*(r1 - lambda1*x(2); x(2)*(- r2 + lambda2*x(1) - mu*x(3); x(3)*(- r3 + lambda3*x(2) ;t, x = ode45 (fun1, 0, 20, 100, 40, 6);subplot(1,2,1)plot(t,x(:,1),- ,t,x(:,2),- .,t,x(:,

20、3),:)legend(x1(t),x2(t),x3(t)gridsubplot(1,2,2)plot3(x(:,1),x(:,2),x(:,3)grid2105101520050100150200250 x1(t)x2(t)x3(t)0100200300010203040500510152025從圖中可以猜測從圖中可以猜測x1(t), x2(t), x3(t)是周期函數(shù)是周期函數(shù), 從數(shù)值解近似從數(shù)值解近似定出周期為定出周期為6.25, 用數(shù)值積分可以算出用數(shù)值積分可以算出x1(t), x2(t), x3(t)在一個在一個周期的平均值周期的平均值x1=71, x2=13, x3=11.223

21、. 模型的改進：模型的改進： (略略)4. 模型的評價模型的評價 若不考慮植物、哺乳動物、爬行動物對自身的若不考慮植物、哺乳動物、爬行動物對自身的阻滯增長阻滯增長作用作用, 則從基本模型中可以得出則從基本模型中可以得出, 植物、哺植物、哺乳動物、爬行動物的數(shù)量都是呈周期變化的。若考乳動物、爬行動物的數(shù)量都是呈周期變化的。若考慮植物、哺乳動物、爬行動物對自身的阻滯增長作慮植物、哺乳動物、爬行動物對自身的阻滯增長作用用, 則從改進模型可望得出則從改進模型可望得出, 植物、哺乳動物、爬行植物、哺乳動物、爬行動物的數(shù)量當(dāng)達到一定程度時動物的數(shù)量當(dāng)達到一定程度時, 它們的數(shù)量都穩(wěn)定它們的數(shù)量都穩(wěn)定在一定

22、的范圍內(nèi)。說明植物、哺乳動物、爬行動物在一定的范圍內(nèi)。說明植物、哺乳動物、爬行動物在滿足一定的條件下在滿足一定的條件下, 它們相互依存的數(shù)量變化最它們相互依存的數(shù)量變化最終都將趨于穩(wěn)定。達到現(xiàn)實生活中的生態(tài)平衡。終都將趨于穩(wěn)定。達到現(xiàn)實生活中的生態(tài)平衡。23題目：題目： SARS（嚴(yán)重急性呼吸道綜合癥（嚴(yán)重急性呼吸道綜合癥, 俗稱：非典型肺炎）俗稱：非典型肺炎）是是21世紀(jì)第一個在世界范圍內(nèi)傳播的傳染病。世紀(jì)第一個在世界范圍內(nèi)傳播的傳染病。SARS的爆的爆發(fā)和蔓延給我國的經(jīng)濟發(fā)展和人民生活帶來了很大影響，發(fā)和蔓延給我國的經(jīng)濟發(fā)展和人民生活帶來了很大影響，我們從中得到了許多重要的經(jīng)驗和教訓(xùn)，認(rèn)識

23、到定量地我們從中得到了許多重要的經(jīng)驗和教訓(xùn)，認(rèn)識到定量地研究傳染病的傳播規(guī)律、為預(yù)測和控制傳染病蔓延創(chuàng)造研究傳染病的傳播規(guī)律、為預(yù)測和控制傳染病蔓延創(chuàng)造條件的重要性。請你們對條件的重要性。請你們對SARS 的傳播建立數(shù)學(xué)模型，具的傳播建立數(shù)學(xué)模型，具體要求如下：體要求如下： （1）對附件）對附件1所提供的一個早期的模型，評價其合理所提供的一個早期的模型，評價其合理性和實用性。性和實用性。三、三、SARS模型模型24（2）建立你們自己的模型，說明為什么優(yōu)于）建立你們自己的模型，說明為什么優(yōu)于附件附件1中的模型；中的模型；特別要說明怎樣才能建立一個真正能夠預(yù)測以及能為預(yù)防特別要說明怎樣才能建立一個

24、真正能夠預(yù)測以及能為預(yù)防和控制提供可靠、足夠的信息的模型，這樣做的困難在哪和控制提供可靠、足夠的信息的模型，這樣做的困難在哪里？對于衛(wèi)生部門所采取的措施做出評論，如：提前或延里？對于衛(wèi)生部門所采取的措施做出評論，如：提前或延后后5天采取嚴(yán)格的隔離措施，對疫情傳播所造成的影響做出天采取嚴(yán)格的隔離措施，對疫情傳播所造成的影響做出估計。估計。附件附件2提供的數(shù)據(jù)供參考。提供的數(shù)據(jù)供參考。（3）收集）收集SARS對經(jīng)濟某個方面影響的數(shù)據(jù)，建立相應(yīng)的數(shù)對經(jīng)濟某個方面影響的數(shù)據(jù)，建立相應(yīng)的數(shù)學(xué)模型并進行預(yù)測。學(xué)模型并進行預(yù)測。附件附件3提供的數(shù)據(jù)供參考。提供的數(shù)據(jù)供參考。（4）給當(dāng)?shù)貓罂瘜懸黄ㄋ锥涛模f

25、明建立傳染病數(shù)學(xué)模型）給當(dāng)?shù)貓罂瘜懸黄ㄋ锥涛模f明建立傳染病數(shù)學(xué)模型的重要性。的重要性。 題目：題目：25基本假設(shè)1)假設(shè)所考查人群的總數(shù)恒定，且無病源的輸入和輸出。假設(shè)所考查人群的總數(shù)恒定，且無病源的輸入和輸出。2)將所考查人群分為將所考查人群分為現(xiàn)有病人現(xiàn)有病人、治愈者治愈者、死亡者死亡者、正常人正常人四類。四類。3)假設(shè)已治愈的患者二度感染的概率為假設(shè)已治愈的患者二度感染的概率為0，即患者具有免疫能，即患者具有免疫能 力，不考慮其再感染。力，不考慮其再感染。4)假設(shè)所有患者均為假設(shè)所有患者均為“他人輸入型他人輸入型”患者，即不考慮人群個體患者，即不考慮人群個體自身發(fā)病。自身發(fā)病。5)假

26、設(shè)各類人群在人群總體中分布均勻。假設(shè)各類人群在人群總體中分布均勻。6)假設(shè)已被隔離的人群之間不會發(fā)生交叉感染。假設(shè)已被隔離的人群之間不會發(fā)生交叉感染。7) 不考慮隱性不考慮隱性SARS患者，即只要感染上患者，即只要感染上SARS病毒的患者最終病毒的患者最終都會表現(xiàn)出癥狀都會表現(xiàn)出癥狀.26符號說明)(tX)(tY)(tR)(tDT1L2Lp)(tr符號符號說明現(xiàn)有病人數(shù)累計病人數(shù)累計治愈人數(shù)累計死亡人數(shù)采取強制措施的時間病人的死亡率病人的治愈率采取控制措施后的隔離強度未被隔離的病人平均每人每天感染的人數(shù)27問題分析與準(zhǔn)備 )(tH)(tX)(tR)(tD 把人群分為把人群分為四類四類：正常人群

27、正常人群、患病人群患病人群、治愈人群治愈人群和和死亡死亡人人群，分別用群，分別用 、和表示。 SARS的傳播規(guī)律可分為的傳播規(guī)律可分為“控前控前”和和“控控后后”兩個階段，兩個階段，如圖所示。如圖所示。 近乎自然的傳播模式控 制 前控 制 后政府控制后的傳播模式28 控前模型為近似于自然傳播時的控前模型為近似于自然傳播時的SIR模型，控后模型為介模型，控后模型為介入隔離強度后的微分方程模型，兩個模型中各類人的轉(zhuǎn)化關(guān)入隔離強度后的微分方程模型，兩個模型中各類人的轉(zhuǎn)化關(guān)系如圖所示。系如圖所示。29建立建立SIR和微分方程模型，先作一些數(shù)據(jù)上的準(zhǔn)備。和微分方程模型，先作一些數(shù)據(jù)上的準(zhǔn)備。 SARS的

28、死亡率和治愈率兩個參數(shù)，一般只能通過醫(yī)學(xué)界的死亡率和治愈率兩個參數(shù)，一般只能通過醫(yī)學(xué)界對治病機理的進一步研究加以控制，在短期內(nèi)不會發(fā)生變化。對治病機理的進一步研究加以控制，在短期內(nèi)不會發(fā)生變化。根據(jù)附錄根據(jù)附錄2的所給的累計病人數(shù)、累計死亡人數(shù)、累計治愈人的所給的累計病人數(shù)、累計死亡人數(shù)、累計治愈人數(shù)，我們可以對數(shù)，我們可以對 和和 作作最小平方誤差估計最小平方誤差估計。 1L2L累計病人數(shù)累計死亡人數(shù)1L累計病人數(shù)累計治愈人數(shù)2L死亡率死亡率治愈率治愈率作作線性回歸線性回歸，得到，得到 0530. 01L0695. 02L301.1.控前模型控前模型現(xiàn)有病人數(shù)現(xiàn)有病人數(shù) 假設(shè)某地區(qū)產(chǎn)生第一例

29、假設(shè)某地區(qū)產(chǎn)生第一例SARSSARS病人的時間為病人的時間為T T0 0，在，在（T T0 0,T,T）時段，是近乎于自由傳播的時段，隔離強度為時段，是近乎于自由傳播的時段，隔離強度為0 0，每個病，每個病人每天感染人數(shù)為一常數(shù)。人每天感染人數(shù)為一常數(shù)。 考察考察(t, (t, t)t)時段內(nèi)現(xiàn)有病人數(shù)的變化，應(yīng)該等于時段內(nèi)現(xiàn)有病人數(shù)的變化，應(yīng)該等于t t時時間段新增的病人數(shù)減去死亡和治愈的人數(shù)。間段新增的病人數(shù)減去死亡和治愈的人數(shù)。新增病人 現(xiàn)有病人死亡和治愈病人31 現(xiàn)有病人數(shù)的變化新增病人數(shù)現(xiàn)有病人數(shù)的變化新增病人數(shù)(死亡人數(shù)治死亡人數(shù)治愈人數(shù)愈人數(shù))。我們設(shè)。我們設(shè)r為每個未被隔離的病

30、人每天感為每個未被隔離的病人每天感染的人數(shù)，染的人數(shù)，L1和和L2分別為治愈率和死亡率。則有分別為治愈率和死亡率。則有 )()( )()( )()( t2211ttXLttXLtttXLttXLtttrXtrtX病人數(shù)治愈人數(shù)治愈率病人數(shù)死亡人數(shù)死亡率時間內(nèi)感染人數(shù)每人在新增病人數(shù)病人數(shù)32 于是有 當(dāng)t0時， 累計死亡人數(shù)累計死亡人數(shù) 死亡累計人數(shù)的變化新增死亡人數(shù) 當(dāng)t0時)()()()()()()()()()(2121tXLLtrXttXttXttXLLttrXtXttX)()()()(21tXLLtrXdttdX ttXLtDttD)()()(1)()(1tXLdttdD33 累計治愈

31、人數(shù)累計治愈人數(shù) 治愈累計人數(shù)的變化新增治愈人數(shù)。 累計病人數(shù)累計病人數(shù) 累計病人數(shù)現(xiàn)有病人數(shù)累計死亡人數(shù)累計治愈人數(shù) ttXLtRttR)()()(2)()(2tXLdttdR)()()()(tRtDtXtY34SARS傳播的傳播的控前模型控前模型 )()()()()()()()()()()()(2121tRtDtXtYtXLdttdRtXLdttdDtXLLtrXdttdX0) 0(0) 0(1) 0(1) 0(RDYXTeXtXtLLrt ,)0()()(211)0(0695. 0053. 055. 021XLLr352. 控后模型 控后隔離強度從控前的0變?yōu)閜。未被隔離的病人平均每人每

32、天感染的人數(shù)r隨時間逐漸變化，它從初始的最大值逐漸減小至客觀存在的最小值，從參考文獻中查到。設(shè)每個未被隔離的病人每天感染的人數(shù)設(shè)T為實施強力控制的時間（以天為單位）。 )()(Ttbeatr)(tr用來反映變化 快慢，可以用附件2中的數(shù)據(jù)估計出它的大小。 36SARSSARS傳播的傳播的控后模型控后模型： TttRtDtXtYtXLdttdRtXLdttdDtXLLtXtrpdttdX,)()()()()()()()()()()()()1 ()(2121其中， 6 . 0245. 0)(babeatrt37TteTXtXTtebPTtLLap,)()()1 ()1 ()()1()(21476

33、. 0245. 0Tbaba、47T為一客觀參數(shù)，可以從文獻中查到。由于3月5日第一例SARS進入北京，是記時的起點；4月20日即為38p和為待估計的參數(shù) 根據(jù)附件2中的數(shù)據(jù)，將各時刻累計病人數(shù)減去累計治愈人數(shù)再減去死亡人數(shù)，可得到現(xiàn)有病人數(shù)，估計p和的值。估計時我們按均方最小誤差原則，用SPSS軟件計算出其估計值分別為%65P02. 0 根據(jù)以上求出的解，作出了現(xiàn)有病人數(shù)、累計死亡人數(shù)、累計治愈人數(shù)、累計病人數(shù)的曲線圖，如圖所示。其中，打點的是實際公布數(shù)據(jù)。 3940414243模型檢驗與結(jié)模型檢驗與結(jié)果分析果分析(1) 靈敏度分析靈敏度分析根據(jù)我們所建的根據(jù)我們所建的模型，衛(wèi)生部門通模型，

34、衛(wèi)生部門通常可以采取兩種方常可以采取兩種方案對疫情進行有效案對疫情進行有效控制。一是改變控控制。一是改變控制時間點制時間點T；二是改；二是改變控制強度變控制強度P?，F(xiàn)在。現(xiàn)在我們分別考察他們我們分別考察他們對模型的影響。對模型的影響。 隔離強度隔離強度P對對的模型影響的模型影響圖圖5 隔離強度對的模型影響隔離強度對的模型影響44隔離強度隔離強度累計病人數(shù)累計病人數(shù)555569966996656528272827757513391339p 表表1 由圖由圖5和表和表1可以看出：可以看出： 隔離強度隔離強度75%與隔離強度與隔離強度65%相比，可使發(fā)病總?cè)藬?shù)減小相比，可使發(fā)病總?cè)藬?shù)減小1500人左

35、右。人左右。 隔離強度隔離強度65%與隔離強度與隔離強度55%相比，可使發(fā)病總?cè)藬?shù)減小相比，可使發(fā)病總?cè)藬?shù)減小4000人左右。人左右。 說明隔離強度，對疫情的傳播具有極大的敏感度和相關(guān)性。說明隔離強度，對疫情的傳播具有極大的敏感度和相關(guān)性。 45 控制時間控制時間T T對的模型影響對的模型影響 46圖圖6 控制時間對的模型影響控制時間對的模型影響47 表表2控制時間累計病人數(shù)累計病人數(shù)延后延后5 5天天53825382延后延后4 4天天47294729延后2天373337334月20日28792879提前2天2764提前提前4 4天天1576提前提前5 5天天1621T 由圖由圖6和表和表2可

36、以看出：控制時間的提前或延后，對累計病人可以看出：控制時間的提前或延后，對累計病人影響顯著。影響顯著。 說明控制時間說明控制時間T，對疫情的傳播具有極大的敏感度和相關(guān)性。，對疫情的傳播具有極大的敏感度和相關(guān)性。48 計算機模擬檢驗從以上曲線可以看出：計算機模擬結(jié)果與模型計算結(jié)果有著良從以上曲線可以看出：計算機模擬結(jié)果與模型計算結(jié)果有著良好的一致性。本模型是可以信賴的好的一致性。本模型是可以信賴的SARS傳播模型。傳播模型。 49模型的評價模型的評價模型的模型的優(yōu)點優(yōu)點： 本文中所建立的是一個連續(xù)的微分方程模型，它從機理上本文中所建立的是一個連續(xù)的微分方程模型，它從機理上準(zhǔn)確地描述了每一時刻的現(xiàn)

37、有病人、治愈者、死亡者的變化規(guī)準(zhǔn)確地描述了每一時刻的現(xiàn)有病人、治愈者、死亡者的變化規(guī)律，消除了離散模型在處理非整數(shù)天數(shù)時的困難，機理合理、律，消除了離散模型在處理非整數(shù)天數(shù)時的困難，機理合理、方法直觀、實用，結(jié)果與實際數(shù)據(jù)擬合的很好。方法直觀、實用，結(jié)果與實際數(shù)據(jù)擬合的很好。 該模型根據(jù)附錄給出的數(shù)據(jù)設(shè)置變量，各變量之間相互該模型根據(jù)附錄給出的數(shù)據(jù)設(shè)置變量，各變量之間相互影響，關(guān)系明確；同時設(shè)定的參數(shù)合情合理，意義明確，消除影響，關(guān)系明確；同時設(shè)定的參數(shù)合情合理，意義明確，消除了人為因素對模型結(jié)果的影響。了人為因素對模型結(jié)果的影響。 建立的微分方程穩(wěn)定性較好，給出了模型的收斂性條件，建立的微分

38、方程穩(wěn)定性較好，給出了模型的收斂性條件，即隔離強度達到多少才能控制疫情，對政府的決策有指導(dǎo)意義。即隔離強度達到多少才能控制疫情，對政府的決策有指導(dǎo)意義。 該模型針對不同隔離強度進行分段研究，能夠方便有效的該模型針對不同隔離強度進行分段研究，能夠方便有效的預(yù)測疫情趨勢。欲對某疫區(qū)進行預(yù)測，只需對參數(shù)進行估計，預(yù)測疫情趨勢。欲對某疫區(qū)進行預(yù)測，只需對參數(shù)進行估計，給出初值帶入方程即可。給出初值帶入方程即可。50模型的模型的缺點缺點： 為了簡化模型的復(fù)雜性，設(shè)定隔離強度，治愈率、死亡為了簡化模型的復(fù)雜性，設(shè)定隔離強度，治愈率、死亡率等參數(shù)在一定階段不發(fā)生變化，而實際情況下，隨著感染率等參數(shù)在一定階段

39、不發(fā)生變化，而實際情況下，隨著感染人數(shù)的減少，其會發(fā)生變化，還需要針對具體情況做具體分人數(shù)的減少，其會發(fā)生變化，還需要針對具體情況做具體分析。析。 模型給出的把人群的每一個個體、每一個地區(qū)視為相同模型給出的把人群的每一個個體、每一個地區(qū)視為相同的，忽略了性別、年齡結(jié)構(gòu)以及地區(qū)差異對隔離措施強度、的，忽略了性別、年齡結(jié)構(gòu)以及地區(qū)差異對隔離措施強度、控制時間等參數(shù)的影響等，而事實上，個體免疫力與個體年控制時間等參數(shù)的影響等，而事實上，個體免疫力與個體年齡因素有關(guān)的，同時不同地域?qū)σ咔榈内厔菀灿杏绊懀写g因素有關(guān)的，同時不同地域?qū)σ咔榈内厔菀灿杏绊懀写倪M。改進。 忽略了人口流動給該地區(qū)傳染病帶

40、來的影響，而實際上忽略了人口流動給該地區(qū)傳染病帶來的影響，而實際上SARS的傳染源多為輸入性病人。如果考慮人口流動，模型的傳染源多為輸入性病人。如果考慮人口流動，模型要加以改進。要加以改進。51問題的推廣與應(yīng)用：問題的推廣與應(yīng)用： 傳染病對人類的威脅與禍害由來已久，自從人類傳染病對人類的威脅與禍害由來已久，自從人類開始向文明社會邁進，病毒就已不斷的襲擊人類。當(dāng)開始向文明社會邁進，病毒就已不斷的襲擊人類。當(dāng)某種傳染傳染病病菌首次侵入缺乏患病經(jīng)驗的種群時，某種傳染傳染病病菌首次侵入缺乏患病經(jīng)驗的種群時，往往會爆發(fā)大規(guī)模的傳入病，造成嚴(yán)重后果。雖然隨往往會爆發(fā)大規(guī)模的傳入病，造成嚴(yán)重后果。雖然隨著人

41、類的醫(yī)學(xué)研究的發(fā)展與突破，已經(jīng)能夠有效的防著人類的醫(yī)學(xué)研究的發(fā)展與突破，已經(jīng)能夠有效的防治和控制許多傳染病，但是由于病毒的遺傳與變異，治和控制許多傳染病，但是由于病毒的遺傳與變異，可能會出現(xiàn)新的突發(fā)性傳染病。可能會出現(xiàn)新的突發(fā)性傳染病。52 大面積、大規(guī)模突發(fā)性傳染病具有蔓延迅速、來勢兇猛、大面積、大規(guī)模突發(fā)性傳染病具有蔓延迅速、來勢兇猛、難以預(yù)防與治療的特點。難以預(yù)防與治療的特點。 傳染病流行過程的研究與其它學(xué)科傳染病流行過程的研究與其它學(xué)科有所不同，不能通過在人群中進行科學(xué)試驗的方式獲得科學(xué)準(zhǔn)有所不同，不能通過在人群中進行科學(xué)試驗的方式獲得科學(xué)準(zhǔn)確的數(shù)據(jù)。在人群中作傳染病試驗，來取得傳染病

42、流行的數(shù)據(jù)確的數(shù)據(jù)。在人群中作傳染病試驗，來取得傳染病流行的數(shù)據(jù)的作法是極不人道也是不可行的。數(shù)學(xué)模型是研究傳染病的重的作法是極不人道也是不可行的。數(shù)學(xué)模型是研究傳染病的重要工具它有助于研究影響疾病傳播的社會和生物機理的相互作要工具它有助于研究影響疾病傳播的社會和生物機理的相互作用，能使我們判斷流行病傳播過程各種因素的相互作用；能夠用，能使我們判斷流行病傳播過程各種因素的相互作用；能夠幫助政府、醫(yī)學(xué)界和科學(xué)界提供治療和控制措施由于上述原因，幫助政府、醫(yī)學(xué)界和科學(xué)界提供治療和控制措施由于上述原因，我們通常主要依據(jù)機理的方法來建力數(shù)學(xué)模型。我們通常主要依據(jù)機理的方法來建力數(shù)學(xué)模型。 我們可以通過收集分析從已有