中國(guó)石油大學(xué)畢業(yè)論文答辯ppt

1、,,,,,,,,,,,,,,,,,,,中國(guó)石油大學(xué)理學(xué)院,,,姓名：韋吉黛 學(xué)號(hào)：12092207,一類具有耐藥性的HIV模型,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,目 錄,CONTENTS,,,,,,,,,,,,,,,,,,,,,,1.HIV: HIV被稱為人類免疫缺陷病毒，英文為Human Immunodeficiency Virus ，HIV是一種專門感染人類免疫系統(tǒng)細(xì)胞的慢病毒，是逆轉(zhuǎn)

2、錄病毒。HIV大致呈球狀，其直徑為1200納米左右，由核殼蛋白和脂質(zhì)膜包裹。HIV病毒主要攻擊人體的T淋巴細(xì)胞系統(tǒng)，病毒侵入機(jī)體細(xì)胞后將會(huì)和細(xì)胞壁整合在一起，并且終生難以消除。人體被HIV病毒感染后，被感染者的血液、腦脊液、乳汁、陰道分泌物、精液和有神經(jīng)癥狀的腦組織液中。所以HIV病毒主要通過(guò)性接觸、血液、和母嬰等方式進(jìn)行感染。2.艾滋病：又名獲得性免疫缺陷綜合癥，英文為Acquired Immunodeficiency Syndro

3、me。目前國(guó)際上公認(rèn)的艾滋病特征為HIV病毒靶向性侵入CD4+T淋巴細(xì)胞（為方便書寫，以下簡(jiǎn)稱T細(xì)胞）引起該細(xì)胞數(shù)量減少和功能破壞以及全身免疫異常激活，從而引起各種機(jī)會(huì)性感染和腫瘤發(fā)生。3.HIV病毒顆粒感染T細(xì)胞的過(guò)程： 病毒進(jìn)入人體后，首先借助于反轉(zhuǎn)錄酶的作用，反轉(zhuǎn)錄成DNA前病毒，然后整合到宿主的DNA中，其轉(zhuǎn)化的過(guò)程包括吸附和穿入、反轉(zhuǎn)錄和整合、轉(zhuǎn)錄、翻譯、病毒核心顆粒的裝配和釋放。4.治療艾滋病的藥物：目前治療艾滋病的藥物

4、主要有四種。 （1）逆轉(zhuǎn)錄酶抑制劑(RTI)、（2）蛋白酶抑制劑(PI)、（3） HIV整合酶抑制劑（4）抑制HIV病毒的進(jìn)入抑制劑（ELS） 當(dāng)今世界上，已經(jīng)研究出許多療法來(lái)治療艾滋病，如高效抗逆轉(zhuǎn)錄療法、免疫調(diào)節(jié)療法、瘧疾療法和中醫(yī)藥療法等等。其中，在中國(guó)，應(yīng)用最廣泛最有效的是高效抗逆轉(zhuǎn)錄病毒療法（簡(jiǎn)稱HAART，也叫做雞尾酒療法），是將抗逆轉(zhuǎn)錄酶抑制劑和蛋白酶抑制劑的結(jié)合的治療方法。,,,,,,,,,,,,,,,,

5、,,,,,,近幾年來(lái)，很多學(xué)者對(duì)HIV病毒感染的動(dòng)力學(xué)模型做了大量的研究，患者感染HIV病毒以后，體內(nèi)的T細(xì)胞濃度及艾滋病病毒顆粒濃度都會(huì)變化。把體內(nèi)的T細(xì)胞的濃度和HIV病毒的濃度作為研究對(duì)象，建立數(shù)學(xué)模型來(lái)分析艾滋病在人體內(nèi)的進(jìn)展，基本病毒動(dòng)力學(xué)模型可以用方程(1-1)描述[6][7] 其中 代表了t時(shí)刻健康T細(xì)胞的濃度、 為t時(shí)刻受感染的T細(xì)胞的濃度、 為血液里游離的HIV病毒顆粒濃度。s表示健康T細(xì)胞的細(xì)胞生

6、成率、 表示健康細(xì)胞的死亡率、k是HIV病毒的感染率、 為被感染細(xì)胞的死亡率、 帶表了每個(gè)被感染T細(xì)胞死亡時(shí)釋放的HIV病毒顆粒個(gè)數(shù)、 表示病毒的清除率。,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,在前面的內(nèi)容中，我們已經(jīng)介紹了HIV感染動(dòng)力學(xué)的基本模型（1-1）。在此模型的基礎(chǔ)上，我們將T細(xì)胞的自身增殖和不同的感染源考慮進(jìn)模型中，通過(guò)對(duì)HIV感染動(dòng)力

7、學(xué)基本模型的適當(dāng)修改，得到具有耐藥性的HIV模型（2-1） （2-1）

8、 其中第一個(gè)方程表示的是健康T細(xì)胞的來(lái)源是由胸腺等免疫器官產(chǎn)生及T細(xì)胞的自身增殖，損失是T細(xì)胞的死亡和被感染。第二個(gè)方程描述了受敏感株感染的T細(xì)胞的來(lái)源是被敏感株感染，損失是受敏感株感染T細(xì)胞的自然死亡以及突變。第四個(gè)方程描述了受耐藥株感染的T細(xì)胞的來(lái)源是被耐藥株感染和由受敏感株感染的T細(xì)胞突變而來(lái)，損失是受耐藥株感染T細(xì)胞的自然死亡。第三個(gè)和第五個(gè)方程描述了病毒顆粒的來(lái)源是由受感染T細(xì)胞中的病毒自我復(fù)制產(chǎn)生，損

9、失是病毒顆粒被清除,,,,,,,,,,,,,,,,,,,,,解該方程組可得模型（2-1）的三個(gè)平衡點(diǎn) 、,,,,,,,,,,,,,,,,,,,,,,,2.3.1 當(dāng) 時(shí)平衡點(diǎn)的局部穩(wěn)定性分析 定理2.1 當(dāng)基本再生數(shù) 且 時(shí)，模型（2-1）的未感染平衡點(diǎn) 是局部漸近穩(wěn)定的；當(dāng) 或 時(shí)，未感染平衡點(diǎn) 是不穩(wěn)定的。 定

10、理2.2 當(dāng) ， 且 時(shí)， 是局部漸近穩(wěn)定的；當(dāng) 或 或 時(shí)， 是不穩(wěn)定的。 定理2.3 當(dāng) 且 且 時(shí)，平衡點(diǎn) 是局部漸近穩(wěn)定的；當(dāng) 或 或 時(shí)，平衡點(diǎn) 是不穩(wěn)定的。2.3.2 當(dāng) 時(shí)平衡點(diǎn)的局部穩(wěn)定性分析 前面的三個(gè)定理中， 的取值范圍為 ，接下來(lái)我們討論在 （也就是不考慮敏感株到耐藥株的轉(zhuǎn)化率）的情況下，系統(tǒng)（2-

11、1）的平衡點(diǎn)的穩(wěn)定性。 定理2.4 當(dāng) 時(shí)，系統(tǒng)（2-1）的平衡點(diǎn)的穩(wěn)定性有以下三種情況：（1）對(duì)于未感染平衡點(diǎn) ，當(dāng) 且 成立時(shí)， 是局部漸近穩(wěn)定的；當(dāng) 或 時(shí)， 是不穩(wěn)定的。（2）對(duì)于僅存在耐藥株平衡點(diǎn) ，當(dāng) 且 且 時(shí)， 是局部漸近穩(wěn)定的；當(dāng) 或 或 時(shí)， 是不穩(wěn)定的。（3）對(duì)于平衡點(diǎn) ， 使得由平衡點(diǎn) 可以得出僅有敏感株存在的平衡點(diǎn) 。并且

12、當(dāng) ， 和 同時(shí)成立時(shí)， 是局部漸近穩(wěn)定的；當(dāng) 或 或 時(shí)， 是不穩(wěn)定的。,,,,,,,,,,,,,,,,,,,,,在前面的內(nèi)容中，我們已經(jīng)求得模型的平衡點(diǎn)以及平衡點(diǎn)具有漸近穩(wěn)定性的充分條件。接下來(lái)我們將利用Matlab進(jìn)行數(shù)值模擬，驗(yàn)證前面四個(gè)定理中關(guān)于模型平衡點(diǎn)穩(wěn)定性結(jié)論的合理性。首先模擬未感染平衡點(diǎn) 𝐸 0 的穩(wěn)定性：其參數(shù)值為表3-1中的數(shù)據(jù)Data1，基本再生數(shù) 𝑅 &

13、#119903; =0.0686<1 ， 𝑅 𝑠 =0.0571<1，滿足定理2.1的充分條件，所以未感染平衡點(diǎn) 𝐸 0 是局部漸近穩(wěn)定的。選取各變量的初值為𝑇 0 =1000, 𝑇 𝑠 ? = 10 ?3 , 𝑉 𝑠 = 10 ?3 , 𝑇 𝑟 ? = 10 ?3

14、 , 𝑉 𝑟 = 10 ?3 。運(yùn)用Matlab進(jìn)行編程，得出結(jié)果如圖3-1所示。從圖3-1可以看出，未感染T細(xì)胞的濃度隨著時(shí)間的增長(zhǎng)不斷增大，當(dāng)時(shí)間達(dá)到200days時(shí)，未感染T細(xì)胞的濃度基本穩(wěn)定在1300ul-1，此時(shí)其它變量已經(jīng)達(dá)到濃度為0 ul-1的穩(wěn)定狀態(tài)。結(jié)合上面計(jì)算出的理論結(jié)果及數(shù)值模擬結(jié)果可知，未感染平衡點(diǎn) 𝐸 0 在 𝑅 𝑟 =0.068

15、6<1， 𝑅 𝑠 =0.0571<1的條件下是局部漸近穩(wěn)定的。,,,,,,,,,,,,,,,,,,,,,,其次，模擬僅存在耐藥株的平衡點(diǎn) 𝐸 𝑟 的穩(wěn)定性。其參數(shù)值為表3-1中的Data2，選取初值為𝑇 0 =1000, 𝑇 𝑠 ? = 10 ?3 , 𝑉 𝑠 = 10 ?3 ,

16、𝑇 𝑟 ? = 10 ?3 , 𝑉 𝑟 = 10 ?3 。運(yùn)用Matlab軟件進(jìn)行編程，得出結(jié)果如圖3-2所示。從圖3-2中可以看出，在400days以后，所有變量都基本達(dá)到穩(wěn)定狀態(tài)。利用給定的參數(shù)算出理論結(jié)果 𝑅 𝑟 =1.7139>1,𝜎=0.04<1。滿足定理2.2的充分條件，所以平衡點(diǎn) ⻒

17、4; 𝑟 是局部漸近穩(wěn)定的。由此可知理論結(jié)果和數(shù)值模擬結(jié)果是符合的。,,,,,,,,,,,,,,,,,,,,,,最后，模擬同時(shí)存在敏感株和耐藥株的平衡點(diǎn) 𝐸 𝑒 的穩(wěn)定性。其參數(shù)值為表3-1中的Data3，選取初值為𝑇 0 =1000, 𝑇 𝑠 ? = 10 ?3 , 𝑉 𝑠 = 10 ?3 , w

18、879; 𝑟 ? = 10 ?3 , 𝑉 𝑟 = 10 ?3 。運(yùn)用Matlab軟件進(jìn)行編程，得出結(jié)果如圖3-3所示。從圖3-3可以看出，在500days之后，所有的變量都基本處于穩(wěn)定狀態(tài)。利用所給參數(shù)計(jì)算出理論結(jié)果 𝑅 𝑠 =2.1381>1,𝜎=1.0667>1,𝑝? 2𝑝 𝑇

19、; 𝑒 𝑇 𝑚𝑎𝑥 ?𝑑=0.0056<0，滿足定理2.3的充分條件，所以平衡點(diǎn) 𝐸 𝑒 是局部漸近穩(wěn)定的。由此可知理論結(jié)果和數(shù)值模擬結(jié)果是符合的。,,,,,,,,,,,,,,,,,,,,,,特別地，當(dāng)𝜇=0時(shí)，我們模擬平衡點(diǎn) 𝐸 𝑠 的穩(wěn)定性，其參數(shù)

20、值為表3-1中的Data4，選取初值為𝑇 0 =1000, 𝑇 𝑠 ? = 10 ?3 , 𝑉 𝑠 = 10 ?3 , 𝑇 𝑟 ? = 10 ?3 , 𝑉 𝑟 = 10 ?3 。運(yùn)用Matlab軟件進(jìn)行編程，得出結(jié)果如圖3?4所示。從圖3-4中可以看出，在500days之后，所有變量都基本

21、處于穩(wěn)定狀態(tài)。利用所給參數(shù)算出理論結(jié)果 𝑅 𝑠 =2.14>1,𝜎=36>1,𝑝? 2𝑝 𝑇 𝑠 𝑇 𝑚𝑎𝑥 ?𝑑=?0.005<0，滿足定理2.4的充分條件，所以平衡點(diǎn) 𝐸 𝑠 是局部漸近穩(wěn)

22、定的。由此可知理論結(jié)果和數(shù)值模擬結(jié)果是符合的。,,,,,,,,,,,,,,,,,,,,,,當(dāng)不考慮T細(xì)胞的自我增值時(shí)，Logistic增長(zhǎng)將不考慮進(jìn)模型中，則可得另一類具有耐藥性的HIV模型，具體模型如下：接下來(lái)，我們運(yùn)用Matlab進(jìn)行數(shù)值模擬來(lái)比較模型（2-1）和模型（3-1）的區(qū)別。其參數(shù)值為𝜆=6，𝑝=0.01, 𝑇 𝑚𝑎𝑥 =1

23、500,𝑑=0.01, 𝑘 𝑠 =2.4× 10 ?5 , 𝑘 𝑟 =2× 10 ?5 ,𝜇=3× 10 ?8 ,𝛿=1, 𝑁 𝑠 =1500, 𝑁 𝑟 =1800,𝑐=23，運(yùn)用Matlab軟件進(jìn)行編程，得出結(jié)果如圖

24、3-5所示。由圖3-5可以看出，Logistic增長(zhǎng)模型的各變量的取值都明顯比非Logistic增長(zhǎng)模型的取值要大。這說(shuō)明，不考慮T細(xì)胞自身增殖的模型會(huì)低估健康T細(xì)胞濃度，被感染細(xì)胞濃度，敏感株的濃度和耐藥株的濃度，從而不能真實(shí)反映HIV感染者的實(shí)際情況。,,,,,,,,,,,,,,,,,,,,,,本文在基本的HIV病毒感染動(dòng)力學(xué)模型的基礎(chǔ)上建立了一類具有耐藥性的HIV感染模型。第二章中求出了模型（2-1)的三類平衡點(diǎn)，通過(guò)微分方程穩(wěn)

25、定性分析方法分別證明了各平衡點(diǎn)穩(wěn)定的充分條件。在第三章中結(jié)合醫(yī)學(xué)上的實(shí)際數(shù)據(jù)對(duì)模型（2-1）進(jìn)行數(shù)值模擬，驗(yàn)證了平衡點(diǎn) 的穩(wěn)定性并與非Logistic增長(zhǎng)模型進(jìn)行比較。通過(guò)第三章中給出的比較結(jié)果，可以看出考慮Logistic增長(zhǎng)的具有耐藥性的HIV模型比非Logistic增長(zhǎng)模型更具有現(xiàn)實(shí)意義，更接近真實(shí)的HIV感染過(guò)程。,,,,,,,,,,,,,,,,,,,,,,,四年的大學(xué)生活即將結(jié)束，在此，我要感謝所有教過(guò)我的老師，謝謝他們授予我

26、知識(shí)；感謝所有關(guān)心我、幫助我的老師和同學(xué)，他們的關(guān)懷溫暖了我。 首先，我要特別感謝我的導(dǎo)師王艷老師。大二時(shí)王艷老師是《常微分方程》這門課的任課老師，我很喜歡她的上課方式，她以自己的親身經(jīng)歷以及耐心的講解和傾聽(tīng)教給我們很多知識(shí)。課上她讓我收獲了專業(yè)知識(shí)，課下和她的交流讓我對(duì)自己未來(lái)的發(fā)展有了一定的規(guī)劃。因此在畢業(yè)設(shè)計(jì)選題時(shí)我選擇了王老師的課題。在論文的撰寫過(guò)程中，王艷老師給我講解了很多新知識(shí)，并介紹了很多和課題有關(guān)的思路，讓

27、我對(duì)畢設(shè)的總體有了大概的把握。在畢設(shè)的過(guò)程中，我有很多不懂的地方，但老師每次都很耐心地講解，直到我能領(lǐng)悟問(wèn)題的思想。最后，我想對(duì)王艷老師再次說(shuō)聲謝謝，謝謝你的的教導(dǎo)、幫助以及你給予我的理解和支持！ 在這即將畢業(yè)的季節(jié)，我對(duì)中國(guó)石油大學(xué)這一個(gè)培養(yǎng)了我四年的校園充滿不舍和感動(dòng)，我感激這所大學(xué)， 感激大學(xué)里的很多人、事、物。最后，我要感激遠(yuǎn)方的家人，感謝他們一直默默地為我遮風(fēng)擋雨！,,,,,,,,,,,,,,,,,,,,,,,

28、,,,,,,,,,,第一章,引言,1.HIV: HIV被稱為人類免疫缺陷病毒，英文為Human Immunodeficiency Virus ，HIV是一種專門感染人類免疫系統(tǒng)細(xì)胞的慢病毒，是逆轉(zhuǎn)錄病毒。HIV大致呈球狀，其直徑為1200納米左右，由核殼蛋白和脂質(zhì)膜包裹。HIV病毒主要攻擊人體的T淋巴細(xì)胞系統(tǒng)，病毒侵入機(jī)體細(xì)胞后將會(huì)和細(xì)胞壁整合在一起，并且終生難以消除。人體被HIV病毒感染后，被感染者的血液、腦脊液、乳汁、陰道分泌物、精

29、液和有神經(jīng)癥狀的腦組織液中。所以HIV病毒主要通過(guò)性接觸、血液、和母嬰等方式進(jìn)行感染。2.艾滋?。河置@得性免疫缺陷綜合癥，英文為Acquired Immunodeficiency Syndrome。目前國(guó)際上公認(rèn)的艾滋病特征為HIV病毒靶向性侵入CD4+T淋巴細(xì)胞（為方便書寫，以下簡(jiǎn)稱T細(xì)胞）引起該細(xì)胞數(shù)量減少和功能破壞以及全身免疫異常激活，從而引起各種機(jī)會(huì)性感染和腫瘤發(fā)生。,HIV和艾滋病,,,,,,,,,,,,,,,,,,,,,

30、,,,,,,,,,,,,第一章,引言,3.HIV病毒顆粒感染T細(xì)胞的過(guò)程：病毒進(jìn)入人體后，首先借助于反轉(zhuǎn)錄酶的作用，反轉(zhuǎn)錄成DNA前病毒，然后整合到宿主的DNA中，其轉(zhuǎn)化的過(guò)程包括吸附和穿入、反轉(zhuǎn)錄和整合、轉(zhuǎn)錄、翻譯、病毒核心顆粒的裝配和釋放。4.治療艾滋病的藥物：目前治療艾滋病的藥物主要有四種。（1）逆轉(zhuǎn)錄酶抑制劑(RTI)（2）蛋白酶抑制劑(PI) （3） HIV整合酶抑制劑（4）抑制HIV病毒的進(jìn)入抑制劑（ELS）

31、 當(dāng)今世界上，已經(jīng)研究出許多療法來(lái)治療艾滋病，如高效抗逆轉(zhuǎn)錄療法、免疫調(diào)節(jié)療法、瘧疾療法和中醫(yī)藥療法等等。其中，在中國(guó)，應(yīng)用最廣泛最有效的是高效抗逆轉(zhuǎn)錄病毒療法（簡(jiǎn)稱HAART，也叫做雞尾酒療法），是將抗逆轉(zhuǎn)錄酶抑制劑和蛋白酶抑制劑的結(jié)合的治療方法。,HIV和艾滋病,,,,,,,,,,,,,,,,,,,,,,,,,,,引言,第一章,預(yù)備知識(shí),,,,,,,,,,,,,,,,,,,,,預(yù)備知識(shí),,,,,,,,,,,,,,,,,,,,,

32、,,,,,,,,,,,,第一章,此處插入文本框,,,,,引理1.1Rout-Hurwitz判據(jù),,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,2,此處插入文本框,,,,,引理1.1Rout-Hurwitz判據(jù),,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,2,此處插入文本框,,此處插入文本框,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,3,,,,,此處插入文本框此處插入文本框此

33、處插入文本框此處插入文本框此處插入文本框,此處插入文本框此處插入文本框此處插入文本框此處插入文本框,此處插入文本框此處插入文本框此處插入文本框此處插入文本框此處插入文本框,此處插入文本框此處插入文本框此處插入文本框,此處插入文本框,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,3,此處插入文本框此處插入文本框此處插入文本框,此處插入文本框此處插入文本框此處插入文本框,此處插入文本框此處插入文本框此處插入文本框,