微分方程的應用

微分方程的應用第1頁，共38頁，2023年，2月20日，星期六2006ICM

題目：

如何在抵抗艾滋病的過程中均衡利用資源。

總任務：為艾滋病的幾種重要情況建立模型，并根據(jù)你的模型給出財政資源的分配方案。第2頁，共38頁，2023年，2月20日，星期六任務描述任務1：在每個大洲中選取一個艾滋病病情嚴重的國家，建立模型，未來50年中，預測在沒有任何干預措施的情況下你所選定的國家中艾滋病感染者的數(shù)量的變化率。并對模型及模型假設作詳細的解釋。第3頁，共38頁，2023年，2月20日，星期六任務描述任務2：估計所選定國家中，從2006年到2050年每年對抗艾滋病所需的國際支援的水平，利用這些財政資源和任務1中建立的模型來預測在所選定的國家中，從2006年到2050年，在下面三種情況下艾滋病感染者數(shù)量的變化率。（1）實施抗逆轉錄酶(ARV)藥物治療（2）注射抗艾滋病毒疫苗（3）實施抗逆轉錄酶藥物治療和注射疫苗第4頁，共38頁，2023年，2月20日，星期六任務描述任務3：重新闡述任務2中建立的3個模型，把抗藥性的產(chǎn)生考慮在內。任務4：寫一份白皮書給聯(lián)合國，在以下三個方面提出建議（1）在實施抗逆轉錄酶藥物治療和注射抗艾滋病毒疫苗兩種方案之間如何合理分配已有資源；（2）跟其他國際事務相比較，抵抗艾滋病應該有怎樣的優(yōu)先級別；（3）如何協(xié)調為抵抗艾滋病的國際捐助。第5頁，共38頁，2023年，2月20日，星期六答案簡介通過對給定的資料和收集到的資料的分析，以下國家被選出來作為分析對象：Zambia、theU.S.、Spain、Australia、Thailand、Brazil。其中對美國和贊比亞做了詳細的分析。第6頁，共38頁，2023年，2月20日，星期六重點講述建立模型，未來50年中，預測在沒有任何干預措施的情況下美國艾滋病感染者的數(shù)量變化。并對模型及模型假設作詳細的解釋。第7頁，共38頁，2023年，2月20日，星期六問題分析預測在數(shù)學上的幾類方法：

第一類：分形、灰度、馬爾科夫、時間序列，適用于短期預測，對歷史數(shù)據(jù)量的要求比較大。

第二類：根據(jù)已有大量數(shù)據(jù)用toolbox工具箱畫圖通過對數(shù)據(jù)的圖像走勢分析進行預測，適用于長期預測，對歷史數(shù)據(jù)量的要求比較大。第三類：微分方程法適用于短期預測和長期預測，且對歷史數(shù)據(jù)量要求不高，既可定性又可定量分析問題，但是對模型準確度要求高，稍有誤差將導致預測結果不準確。

第8頁，共38頁，2023年，2月20日，星期六模型假設(I)模型假設的基本要求：要求一：恰如其分（不能與常識公理現(xiàn)有成熟理論相違背）；

要求二：要自圓其說（每個模型都有不足與缺點需要大膽假設且要在合理假設下嚴密論證自己結論使自己結論合理化）；要求三：假設要明確，切不可模棱兩可讓評委有歧義。第9頁，共38頁，2023年，2月20日，星期六模型假設(II)本題的模型假設：1、艾滋病人不可能被治愈；2、忽略艾滋病的具體分類；3、艾滋病人喪失大部分常人所具有的能力，如：性能力，所以我們假設他們不能傳播艾滋??；4、只考慮三種傳播艾滋病的方式：性交傳播，毒品注射傳播，母嬰傳播。忽略其他傳播方式；第10頁，共38頁，2023年，2月20日，星期六模型假設(III)5、我們不考慮吸毒者戒毒的情況，以簡化艾滋病傳播過程的復雜度；7、不考慮該國的遷入和遷出人口，因為這些人口相對于國家的總人口是很小的，以便簡化模型。6、吸毒傳播的渠道只能是從有吸毒習慣的易感人群中直接轉變到有吸毒習慣的hiv攜帶者；第11頁，共38頁，2023年，2月20日，星期六符號說明(I)符號說明的基本要求：寫作者可參考往年outstanding文章中符號說明部分的格式。確保符號完整無遺漏；定義準確無誤；格式優(yōu)美。第12頁，共38頁，2023年，2月20日，星期六符號說明(II)符號定義S艾滋病易感人群的數(shù)量；S1艾滋病易感且不吸毒人群的數(shù)量；S2艾滋病易感且吸毒人群的數(shù)量；H艾滋病攜帶者數(shù)量;H1

攜帶艾滋病且不吸毒人群的數(shù)量；H2

攜帶艾滋病且吸毒人群的數(shù)量；N

該國人口數(shù)量；A艾滋病患者的數(shù)量；該國人口的出生率；Pinfection艾滋病攜帶者占總人口的比例；每人每次性交感染艾滋病的概率；從S1到S2通過毒品注射的轉化率；第13頁，共38頁，2023年，2月20日，星期六符號說明(III)符號定義dS

S中的人口的死亡率；dH

H中的人口的死亡率；dAA中人口的自然死亡率；C(N)

一名艾滋病攜帶者在單位時間內與S人群中的人發(fā)生性交的次數(shù)，一般與總人口N有關。S中人口占總人的比例；H到A的轉化率；該國家單位時間內的人口增量；單位時間內一個人通過吸毒感染艾滋病的概率；A中艾滋病造成的死亡率；H中人口占總人的比例；第14頁，共38頁，2023年，2月20日，星期六模型分析(I)第15頁，共38頁，2023年，2月20日，星期六模型分析(II)對于上圖中各種傳播方式的定量分析：1、艾滋病通過性交的傳播：假設每次性接觸傳染的概率為，我們把賦有傳染概率的接觸率稱為有效接觸率，即：。它表示一個hiv病毒攜帶者傳染他人hiv病毒的能力，反映了hiv病毒攜帶者的活動能力、環(huán)境條件以及hiv病毒的毒力等因素。假設單位時間內一個hiv攜帶者與其他成員的性接觸次數(shù)稱之為性接觸率，它通常依賴于環(huán)境中總成員數(shù)N，記作。如果被接觸者為易感者，就有可能被傳染。第16頁，共38頁，2023年，2月20日，星期六模型分析(III)但是，只有當hiv攜帶者與總成員中的易感者發(fā)生性接觸時才會使病毒發(fā)生傳染，與總成員中的其他成員接觸并不會發(fā)生病毒的傳染（比如其他的hiv攜帶者），易感者在總成員中所占的比例為S/N。因此，每一個患者對易感者的平均有效接觸率應為，它就是每一個hiv攜帶者平均對易感者的傳播率，簡稱為傳染率。假設hiv攜帶者的總數(shù)量為，從而在單位時間內被所有患者傳染的新成員數(shù)為：第17頁，共38頁，2023年，2月20日，星期六模型分析(IV)2、艾滋病通過毒品注射的傳播：用表示在單位時間內S2類成員與H2

類成員之間的轉化率，則單位時間內從S2類轉化到H2類的總數(shù)為：

無吸毒習慣的易感人群可以間接向有吸毒習慣的hiv攜帶者轉化，因為無吸毒習慣的易感人群可以先變?yōu)橛形玖晳T易感人群再變?yōu)橛形玖晳T的hiv攜帶者，用表示從S1類轉變到S2類的轉化率則單位時間內從S1類轉變到S2類的總數(shù)為:

第18頁，共38頁，2023年，2月20日，星期六模型分析(V)3、艾滋病通過母嬰的傳播：我們知道如果一位母親是hiv攜帶者那么她的孩子將有百分之30的概率在出生時已經(jīng)成為hiv攜帶者[ChinaHIV/AIDSinformationnetwork2006]，用Pinfection表示艾滋病攜帶者占總人口的比例，則每年hiv攜帶者自己繁殖出的攜帶者數(shù)量為：第19頁，共38頁，2023年，2月20日，星期六模型建立

第20頁，共38頁，2023年，2月20日，星期六模型求解(I)模型結果大多用圖來描述，一幅好圖勝過千言萬語。圖的要求：自明。舉例1：fig*

Proportionsofpopulationofeachagegroupofelephants.第21頁，共38頁，2023年，2月20日，星期六模型求解(II)舉例2：Figure*.Aexponentialfitoftheelephantpopulationofeachageset.

第22頁，共38頁，2023年，2月20日，星期六模型求解(III)求解微分方程通常的方法：

對于相對簡單的方程，我們通常使用解析法；對于相對復雜的方程，通常使用數(shù)值法。第23頁，共38頁，2023年，2月20日，星期六模型求解(IV)本題的方程組比較復雜，我們使用數(shù)值法求解。求解結果如下：Figure5PopulationtrendsineachgroupvstimeinU.S.第24頁，共38頁，2023年，2月20日，星期六合理性分析(I)

對自己的結果自圓其說重中之重！技巧：用已有的比較成熟的結論來佐證自己的結論。本問題的合理性分析如下：第25頁，共38頁，2023年，2月20日，星期六合理性分析(II)在傳染病生物學中，人們關心的已感人群和易感人群的變化，在艾滋病傳染模型中，我們關心的是HIV易感人群和HIV攜帶者數(shù)量的變化趨勢。通過研究我們發(fā)現(xiàn)，HIV攜帶者占總人口的比例隨著時間的變化最后趨于穩(wěn)定，如下圖：Figure9TheproportionofHIVcarrierstothetotalpopulationoftheU.S.第26頁，共38頁，2023年，2月20日，星期六合理性分析(III)對于HIV攜帶者來說，自然環(huán)境對它的承載能力由很多因素決定，比如：易感人群的數(shù)量，總人口的數(shù)量，出生率，死亡率等，其中易感人群的數(shù)量是主要的決定因素．由于HIV攜帶者是由易感者轉化來的，所以HIV攜帶者數(shù)量的增長會抑制易感人群的增長。從圖9可以看出，美國hiv攜帶者的數(shù)量占人群總數(shù)的比例隨時間的變化分為兩個階段：增長期和穩(wěn)定期．在增長期中，該比例隨時間以羅杰斯迪曲線的形式增加，最后達到一個峰值．這是因為在增長期的前段，環(huán)境承受能力大于HIV攜帶者的數(shù)量，所以ＨＩＶ攜帶者數(shù)量以指數(shù)增長，而在增長其后段HIV攜帶者的數(shù)量越來越接近環(huán)境的承載能力，所以其增長速率減慢，以至為０．在穩(wěn)定期，該比例隨時間以羅杰斯迪曲線的形式減少，最后趨近一個穩(wěn)定值．第27頁，共38頁，2023年，2月20日，星期六合理性分析(IV)這是因為種群數(shù)量超過了環(huán)境的承載能力．但是為什么當種群數(shù)量達到環(huán)境承載能力時，它不立即停止增長以至于超過了環(huán)境的承載能力呢？為了對這個問題作出解釋，我們引入人口慣性的概念．所謂人口慣性就是人口抵抗自身變化的趨勢。正是由于人口慣性的存在，當群數(shù)量達到環(huán)境承載能力時，種群數(shù)量增長不會立即停止，而是在一個穩(wěn)定值周圍以羅杰斯迪曲線的速率來回擺動，正如物理學中波的傳播．圖示如下：第28頁，共38頁，2023年，2月20日，星期六合理性分析(V)Figure10Thecurvesofoscillatingconvergenceinthestablephas。

第29頁，共38頁，2023年，2月20日，星期六合理性分析(VI)由于這種擺動的幅度越來越小，所以在穩(wěn)定期的后期，它可以被忽略．這也就是穩(wěn)定期存在的原因．在自然界中實事也確實如此，就以羅杰斯地曲線為例當種群數(shù)量一開始無論是增長還是減少的變化率在一開始時卻都是很慢的，隨后才慢慢增大的，這就好比用一個外力在拉一個靜止的物體，物體的速度也是由慢到快增長的，并且物體的運動方向是與外力同方向的。那么在種群變化中是誰扮演了外力的角色，我們隊伍經(jīng)過討論認為應當是環(huán)境能夠承載此物種數(shù)量的上限扮演了這種外力，當環(huán)境能夠承載此物種數(shù)量的上限大于物種現(xiàn)有的數(shù)量時物種的數(shù)量會增大，環(huán)境能夠承載此物種數(shù)量的上限小于物種現(xiàn)有的數(shù)量時物種的數(shù)量會減小。但是按照達爾文的進化理論物種應當都是希望能夠壯大自己的數(shù)量好保存自己在自然界中競爭的實力，第30頁，共38頁，202