微分方程模型()第一頁，共五十五頁，2022年，8月28日微分方程是研究變化規(guī)律的有力工具，在科技、工程、經(jīng)濟(jì)管理、生態(tài)、環(huán)境、人口和交通各個(gè)領(lǐng)域中有廣泛的應(yīng)用。不少實(shí)際問題當(dāng)我們采用微觀眼光觀察時(shí)都遵循著下面的模式：凈變化率＝輸入率－輸出率（守恒原理）一、微分方程模型簡介第二頁，共五十五頁，2022年，8月28日引例一在凌晨1時(shí)警察發(fā)現(xiàn)一具尸體，測得尸體溫度是29oC，當(dāng)時(shí)環(huán)境的溫度是21oC。1h后尸體溫度下降到27oC，若人體的正常溫度是37oC，估計(jì)死者的死亡時(shí)間。解：設(shè)T(t)為死者在被殺害后t時(shí)刻尸體的溫度；k為比例系數(shù)。由牛頓冷卻定律，得則通解為第三頁，共五十五頁，2022年，8月28日由已知，由因此死者大約是在前一天的夜晚10:35被害的。可得微分方程的特解：，代入解得圖1尸體的溫度下降曲線第四頁，共五十五頁，2022年，8月28日建立微分方程的常用方法1、按變化規(guī)律直接列方程，如：利用人們熟悉的力學(xué)、數(shù)學(xué)、物理、化學(xué)等學(xué)科中的規(guī)律，如牛頓第二定律，放射性物質(zhì)的放射規(guī)律等。對(duì)某些實(shí)際問題直接列出微分方程．2、利用微元分析方法建模根據(jù)已知的規(guī)律或定理，通過尋求微元之間的關(guān)系式得出微分方程。3、模擬近似法，如：在生物、經(jīng)濟(jì)等學(xué)科中，許多現(xiàn)象所滿足的規(guī)律并不很清楚，而且現(xiàn)象也相當(dāng)復(fù)雜，因而需根據(jù)實(shí)際資料或大量的實(shí)驗(yàn)數(shù)據(jù)，提出各種假設(shè)，在一定的假設(shè)下，給出實(shí)際現(xiàn)象所滿足的規(guī)律，然后利用適當(dāng)?shù)臄?shù)學(xué)方法得出微分方程。第五頁，共五十五頁，2022年，8月28日微分方程的建模步驟1、翻譯或轉(zhuǎn)化：在實(shí)際問題中許多表示導(dǎo)數(shù)的常用詞，如“速率”、‘增長”(在生物學(xué)以及人口問題研究中)，“衰變”(在放射性問題中)，以及“邊際的”(在經(jīng)濟(jì)學(xué)中)等．2、建立瞬時(shí)表達(dá)式：根據(jù)自變量有微小改變△t時(shí)，因變量的增量△W，建立起在時(shí)段△t上的增量表達(dá)式，令△t→0，即得到的表達(dá)式．二、微分方程模型第六頁，共五十五頁，2022年，8月28日3、配備物理單位：在建模中應(yīng)注意每一項(xiàng)采用同樣的物理單位．4、確定條件：這些條件是關(guān)于系統(tǒng)在某一特定時(shí)刻或邊界上的信息，它們獨(dú)立于微分方程而成立，用以確定有關(guān)的常數(shù)。為了完整充分地給出問題的數(shù)學(xué)陳述，應(yīng)將這些給定的條件和微分方程一起列出。第七頁，共五十五頁，2022年，8月28日案例1：以為女士每天攝入2500cal食物，1200cal用于基本新陳代謝(即自動(dòng)消耗)，并以每千克體重消耗16cal用于日常鍛煉，其他的熱量轉(zhuǎn)變?yōu)樯眢w的脂肪（設(shè)10000cal可轉(zhuǎn)換成1kg脂肪）。星期天晚上，該女士的體重是57.1526kg，星期四那天她飽餐了一頓，共攝入了3500cal的食物，要求建立一個(gè)通過時(shí)間預(yù)測體重函數(shù)W(t)的數(shù)學(xué)模型，并用它估計(jì)：（1）星期六該女士的體重？（2）為了不增重，每天她最多的攝入量是多少？（3）若不進(jìn)食，N周后她的體重是多少？二、微分方程案例分析第八頁，共五十五頁，2022年，8月28日解1、翻譯或轉(zhuǎn)化：2、配備物理單位：3、建立表達(dá)式：4、確定條件：第九頁，共五十五頁，2022年，8月28日1、“每天”：體重的變化＝輸入一輸出其中輸入指扣除了基本新陳代謝之后的凈重量吸收；輸出是進(jìn)行健身訓(xùn)練時(shí)的消耗．2、上述陳述更好的表示結(jié)構(gòu)式：取天為計(jì)時(shí)單位，記W(t)為t天時(shí)體重(kg)，則：每天的凈吸收量＝2500–1200＝1300（cal)每天的凈輸出量＝16(cal)×W＝16W(cal)轉(zhuǎn)換成脂肪量＝1300–16W（cal）3、體重的變化／天＝(千克／天)第十頁，共五十五頁，2022年，8月28日1、翻譯或轉(zhuǎn)化：2、配備物理單位：3、建立表達(dá)式：4、確定條件：第十一頁，共五十五頁，2022年，8月28日有些量是用能量(cal)的形式給出的，而另外一些量是用重量的形式(cal)給出，考慮單位的匹配，利用單位匹配第十二頁，共五十五頁，2022年，8月28日1、翻譯或轉(zhuǎn)化：2、配備物理單位：3、建立表達(dá)式：4、確定條件：第十三頁，共五十五頁，2022年，8月28日建立表達(dá)式積分后可求得其通解為：（1）當(dāng)時(shí)，每天體重的變化：初始條件為：，代入解出則第十四頁，共五十五頁，2022年，8月28日積分后可求得其通解為：（2）當(dāng)時(shí)，每天體重的變化：初始條件為：，代入解出則第十五頁，共五十五頁，2022年，8月28日積分后可求得其通解為：（2）當(dāng)時(shí)，食物的攝入量恢復(fù)正常初始條件為：，代入解出則第十六頁，共五十五頁，2022年，8月28日最后得到不同階段的微分方程是：第十七頁，共五十五頁，2022年，8月28日（1）代入對(duì)應(yīng)方程，求得現(xiàn)回答上述問題（2）要滿足體重不增，即所以因此每天總卡路里攝取量是1200+914＝2114cal(cal)（3）由于每天不攝取能量，所以解得因此，n周后的體重為第十八頁，共五十五頁，2022年，8月28日案例2在一個(gè)巴基斯坦洞穴里，發(fā)現(xiàn)了具有古代尼安德特人特征的人骨碎片，科學(xué)家們把它們帶到實(shí)驗(yàn)室，作碳14年代測定。分析表明C14與C12的比例僅僅是活組織內(nèi)的6.24％，此人生活在多少年前？（碳14年代測定：活體中的碳有一小部分是放射性同位素C14。這種放射性碳是由于宇宙射線在高層大氣中的撞擊引起的，經(jīng)過一系列交換過程進(jìn)入活組織中，直到在生物體中達(dá)到平衡濃度。這意味著在活體中，C14的數(shù)量與穩(wěn)定的C12的數(shù)量成定比。生物體死亡后，交換過程就停止了，放射性碳便以每年八千分之一的速度減少）第十九頁，共五十五頁，2022年，8月28日（1）問題分析與模型的建立1、放射性衰變的這種性質(zhì)還可描述為“放射性物質(zhì)在任意時(shí)刻的衰變速度都與該物質(zhì)現(xiàn)存的數(shù)量成比例”。而C14的比例數(shù)為每年八千分之一。2、碳14年代測定可計(jì)算出生物體的死亡時(shí)間；所以，我們問題實(shí)際上就是：“這人死去多久了？”若設(shè)t為死后年數(shù)，y(t)為比例數(shù)，則y(t)=C14/C12(mgC14/mgC12)，則上文中最后一句話就給出了我們的微分方程，單位為mgC14/mgC12/yr(與關(guān)鍵詞“速率”相當(dāng))第二十頁，共五十五頁，2022年，8月28日（2）解微分方程的通解為：由初始條件，故有由問題，當(dāng)，代入原方程第二十一頁，共五十五頁，2022年，8月28日案例3、追線問題我緝私艦雷達(dá)發(fā)現(xiàn)，在其正西方距c海里處有一艘走私船正以勻速度a沿直線向北行駛，緝私艦立即以最大的速度b追趕，若用雷達(dá)進(jìn)行跟蹤，保持船的瞬時(shí)速度方向始終指向走私船，試求緝私艦追逐路線和追上的時(shí)間。第二十二頁，共五十五頁，2022年，8月28日?qǐng)D2走私船與緝私艦的位置關(guān)系(c,0)xD(x,y)走私船R(0,at)緝私艇O第二十三頁，共五十五頁，2022年，8月28日幾何關(guān)系第二十四頁，共五十五頁，2022年，8月28日如何消去時(shí)間t?1、速度與路程的關(guān)系：2、分解得：3、求導(dǎo)：4、將第2、3步代入第1步，可得模型第二十五頁，共五十五頁，2022年，8月28日追線模型：模型的解：第二十六頁，共五十五頁，2022年，8月28日解的進(jìn)一步討論（1）若a<b，從而k<1，由積分式得當(dāng)x=0時(shí)，即走私船被緝私艦捕捉前所花的時(shí)間為所跑過的距離為（2）若a=b，即k=1，由積分式得顯然x不能取零值，即緝私艦不可能追上走私船。（3）若a>b，即k>1，顯然緝私艦也不可能追上走私船。第二十七頁，共五十五頁，2022年，8月28日如圖所示一個(gè)容量為2000m3的小湖的示意圖，通過小河A水以0.1m3/s的速度流入，以相同的流量湖水通過B流出。在上午11:05時(shí)，因交通事故一個(gè)盛有毒性化學(xué)物質(zhì)的容器傾翻，圖中X點(diǎn)處注入湖中。在采取緊急措施后，于11:35事故得到控制，但數(shù)量不詳案例4湖泊污染問題的化學(xué)物質(zhì)Z已瀉入湖中，初步估計(jì)Z的量在5~20m3之間。建立一個(gè)模型，通過它來估計(jì)湖水污染程度隨時(shí)間的變化并估計(jì)：（1）湖水何時(shí)到達(dá)污染高峰；（2）何時(shí)污染程度可降至安全水平(<0.05%)圖3小湖示意圖第二十八頁，共五十五頁，2022年，8月28日湖泊污染問題分析設(shè)湖水在t時(shí)的污染程度為C(t)，即每立方米受污染的水中含有Cm3的化學(xué)物質(zhì)和(1-C)m3的清潔水。用分鐘作為時(shí)間t的單位。在0<t<30的時(shí)間內(nèi)，污染物流入湖中小湖示意圖的速率是Z/30(m3/min)，而排出湖外的污染物的速率是60×0.1C(m3/min)，因?yàn)槊苛⒎搅髯叩乃泻蠧m3的污染物，而湖水始終保持2000m3的容積不變，所以可列方程：第二十九頁，共五十五頁，2022年，8月28日由初始條件：，可得微分方程的特解為顯然，t=30時(shí)，污染達(dá)到高峰，所以因污染源被截?cái)啵饰⒎址匠套優(yōu)樗奶亟鉃椋汉泻廴疚锏乃矔r(shí)變化率=污染物流入量－污染物排出量第三十頁，共五十五頁，2022年，8月28日當(dāng)達(dá)到安全水平，即C(t)=0.0005時(shí)，可求出此時(shí)的t=T，即解得Z取不同值時(shí)的濃度C(30)和時(shí)間T51015200.002390.004780.007170.009565527389181014第三十一頁，共五十五頁，2022年，8月28日何為房室系統(tǒng)？在用微分方程研究實(shí)際問題時(shí)，人們常常采用一種叫“房室系統(tǒng)”的觀點(diǎn)來考察問題。根據(jù)研究對(duì)象的特征或研究的不同精度要求，我們把研究對(duì)象看成一個(gè)整體（單房室系統(tǒng)）或?qū)⑵淦史殖扇舾蓚€(gè)相互存在著某種聯(lián)系的部分（多房室系統(tǒng)）。房室具有以下特征：它由考察對(duì)象均勻分布而成，（注：考察對(duì)象一般并非均勻分布，這里采用了一種簡化方法一集中參數(shù)法）；房室中考察對(duì)象的數(shù)量或濃度（密度）的變化率與外部環(huán)境有關(guān)，這種關(guān)系被稱為“交換”且交換滿足著總量守衡。在本節(jié)中，我們將用房室系統(tǒng)的方法來研究藥物在體內(nèi)的分布。在下一節(jié)中，我們將用多房室系統(tǒng)的方法來研究另一問題。兩者都很簡單，意圖在于介紹建模方法。交換環(huán)境內(nèi)部單房室系統(tǒng)均勻分布案例5藥物在體內(nèi)的變化（房室模型）第三十二頁，共五十五頁，2022年，8月28日藥物的分解與排泄（輸出）速率通常被認(rèn)為是與藥物當(dāng)前的濃度成正比的，即：藥物分布的單房室模型單房室模型是最簡單的模型，它假設(shè)：體內(nèi)藥物在任一時(shí)刻都是均勻分布的，設(shè)t時(shí)刻體內(nèi)藥物的總量為x(t)；系統(tǒng)處于一種動(dòng)態(tài)平衡中，即成立著關(guān)系式：藥物的輸入規(guī)律與給藥的方式有關(guān)。下面，我們來研究一下在幾種常見的給藥方式下體內(nèi)藥體的變化規(guī)律。機(jī)體環(huán)境藥物總量圖3-8

假設(shè)藥物均勻分布第三十三頁，共五十五頁，2022年，8月28日情況1快速靜脈注射機(jī)體環(huán)境只輸出不輸入房室其解為：藥物的濃度：

與放射性物質(zhì)類似，醫(yī)學(xué)上將血漿藥物濃度衰減一半所需的時(shí)間稱為藥物的血漿半衰期：負(fù)增長率的Malthus模型在快速靜脈注射時(shí)，總量為D的藥物在瞬間被注入體內(nèi)。設(shè)機(jī)體的體積為V，則我們可以近似地將系統(tǒng)看成初始總量為D，濃度為D/V，只輸出不輸入的房室，即系統(tǒng)可看成近似地滿足微分方程：（3.12）

第三十四頁，共五十五頁，2022年，8月28日情況2恒速靜脈點(diǎn)滴機(jī)體環(huán)境恒定速率輸入房室藥物似恒速點(diǎn)滴方式進(jìn)入體內(nèi)，即:則體內(nèi)藥物總量滿足：（x(0)=0）（3.13）這是一個(gè)一階常系數(shù)線性方程，其解為：或易見：稱為穩(wěn)態(tài)血藥濃度對(duì)于多次點(diǎn)滴，設(shè)點(diǎn)滴時(shí)間為T1，兩次點(diǎn)滴之間的間隔時(shí)間設(shè)為T2，則在第一次點(diǎn)滴結(jié)束時(shí)病人體內(nèi)的藥物濃度可由上式得出。其后T2時(shí)間內(nèi)為情況1。故：（第一次）0≤t≤T1

T1≤t≤T1

+T2

類似可討論以后各次點(diǎn)滴時(shí)的情況，區(qū)別只在初值上的不同。第二次點(diǎn)滴起，患者體內(nèi)的初始藥物濃度不為零。第三十五頁，共五十五頁，2022年，8月28日情況3口服藥或肌注y(t)x(t)k1ykx環(huán)境機(jī)體外部藥物口服藥或肌肉注射時(shí)，藥物的吸收方式與點(diǎn)滴時(shí)不同，藥物雖然瞬間進(jìn)入了體內(nèi)，但它一般都集中與身體的某一部位，靠其表面與肌體接觸而逐步被吸收。設(shè)藥物被吸收的速率與存量藥物的數(shù)量成正比，記比例系數(shù)為K1，即若記t時(shí)刻殘留藥物量為y(t)，則y滿足：D為口服或肌注藥物總量因而：所以：解得：從而藥物濃度：在通常情況下，總有k1>k（藥物未吸收完前，輸入速率通?？偞笥诜纸馀c排泄速率），但也有例外的可能（與藥物性質(zhì)及機(jī)體對(duì)該藥物的吸收、分解能力有關(guān)）。當(dāng)k1>k時(shí)，體內(nèi)藥物量均很小，這種情況在醫(yī)學(xué)上被稱為觸發(fā)翻轉(zhuǎn)（flip-flop）。當(dāng)k1=k時(shí)，對(duì)固定的t，令k→k1取極限（應(yīng)用羅比達(dá)法則），可得出在這種情況下的血藥濃度為：第三十六頁，共五十五頁，2022年，8月28日如下圖給出了上述三種情況下體內(nèi)血藥濃度的變化曲線。容易看出，快速靜脈注射能使血藥濃度立即達(dá)到峰值，常用于急救等緊急情況；口服、肌注與點(diǎn)滴也有一定的差異，主要表現(xiàn)在血藥濃度的峰值出現(xiàn)在不同的時(shí)刻，血藥的有效濃度保持時(shí)間也不盡相同，（注：為達(dá)到治療目的，血藥濃度應(yīng)達(dá)到某一有效濃度，并使之維持一特定的時(shí)間長度）。房室系統(tǒng)我們已求得三種常見給藥方式下的血藥濃度C(t)，當(dāng)然也容易求得血藥濃度的峰值及出現(xiàn)峰值的時(shí)間，因而，也不難根據(jù)不同疾病的治療要求找出最佳治療方案。第三十七頁，共五十五頁，2022年，8月28日國家質(zhì)量監(jiān)督檢驗(yàn)檢疫局2004年5月31日發(fā)布了新的《車輛駕駛?cè)藛T血液、呼氣酒精含量閾值與檢驗(yàn)》國家標(biāo)準(zhǔn)，新標(biāo)準(zhǔn)規(guī)定，車輛駕駛?cè)藛T血液中的酒精含量大于或等于20毫克／百毫升，小于80毫克／百毫升為飲酒駕車（原標(biāo)準(zhǔn)是小于100毫克／百毫升），血液中的酒精含量大于或等于80毫克／百毫升為醉酒駕車（原標(biāo)準(zhǔn)是大于或等于100毫克／百毫升）。五、微分方程綜合案例分析第三十八頁，共五十五頁，2022年，8月28日大李在中午12點(diǎn)喝了一瓶啤酒，下午6點(diǎn)檢查時(shí)符合新的駕車標(biāo)準(zhǔn)，緊接著他在吃晚飯時(shí)又喝了一瓶啤酒，為了保險(xiǎn)起見他呆到凌晨2點(diǎn)才駕車回家，又一次遭遇檢查時(shí)卻被定為飲酒駕車，這讓他既懊惱又困惑，為什么喝同樣多的酒，兩次檢查結(jié)果會(huì)不一樣呢？請(qǐng)你參考下面給出的數(shù)據(jù)（或自己收集資料）建立飲酒后血液中酒精含量的數(shù)學(xué)模型，對(duì)大李碰到的情況做出解釋.第三十九頁，共五十五頁，2022年，8月28日參考數(shù)據(jù)1.人的體液占人的體重的65%至70%，其中血液只占體重的7%左右；而藥物（包括酒精）在血液中的含量與在體液中的含量大體是一樣的。2.體重約70kg的某人在短時(shí)間內(nèi)喝下2瓶啤酒后，隔一定時(shí)間測量他的血液中酒精含量（毫克／百毫升），得到數(shù)據(jù)如下：時(shí)間(小時(shí))0.250.50.7511.522.533.544.55酒精含量306875828277686858515041時(shí)間(小時(shí))678910111213141516酒精含量3835282518151210774第四十頁，共五十五頁，2022年，8月28日問題分析一個(gè)人的血中酒精含量取決于他喝了多少酒、他體內(nèi)原有的酒精含量以及喝酒方式等。由科普知識(shí)知道，酒精是經(jīng)胃腸（主要是肝臟）的吸收與分解進(jìn)體液的。因此本文把酒精的從胃腸（含肝臟）向體液轉(zhuǎn)移情況用如下簡圖直觀地表示：k11為酒精從胃腸滲透到（除體液外）其它地方的速率系數(shù)；k12為酒精從胃腸進(jìn)入體液的速率系數(shù)；k21為酒精在體液中消耗（向外排除或分解或吸收）的速率系數(shù)；f(t)為酒精進(jìn)入胃腸的速率。第四十一頁，共五十五頁，2022年，8月28日由題意，參照房室模型，可建立如下微分方程組：（1）大李在中午12點(diǎn)喝一瓶啤酒時(shí)，即在t=0時(shí)，胃腸中的酒精量x1(0)為一瓶酒中的酒精a與飲酒瓶數(shù)N的乘積Na，而此時(shí)體液中的酒精量y1(0)為0。因此初始條件為體液（或血液）中的酒精的濃度為第四十二頁，共五十五頁，2022年，8月28日（2）大李第二次喝酒時(shí)胃腸和體液中已經(jīng)有酒精，所以在第二次喝酒即t=0時(shí)胃腸中的酒精量x2(0)為N瓶酒中的酒精質(zhì)量Na與第一次喝酒后殘留在胃腸中的酒精質(zhì)量x1(T1)之和，而此時(shí)體液中的酒精量y1(0)為第一次喝酒后殘留在胃腸中的酒精質(zhì)量y1(T1)，因此大李第二次喝酒的模型如下：第四十三頁，共五十五頁，2022年，8月28日解以上微分方程組，得令，，解可轉(zhuǎn)化為第四十四頁，共五十五頁，2022年，8月28日N＝2，運(yùn)用最小二乘擬合法，求解得作圖如下：第四十五頁，共五十五頁，2022年，8月28日將以上數(shù)據(jù)代入問題一的模型中，可求得大李在中午12點(diǎn)飲一瓶啤酒，即N=1時(shí)，到下午6點(diǎn)第一次檢查時(shí)體液中的酒精含量（即血液中的酒精含量）所以大李通過了第一次檢查。第四十六頁，共五十五頁，2022年，8月28日大李第二次喝酒模型的方程解為：考慮到大李在下午6點(diǎn)接受檢查，之后由于停車等待等原因耽誤了大約半個(gè)小時(shí)，假設(shè)大李從第一次檢驗(yàn)到第二次喝酒之間間隔0.5小時(shí)，代入數(shù)據(jù)計(jì)算可得第二次檢驗(yàn)時(shí)，大李血液中酒精含量為：20.2448(毫克/百毫升)。這就解釋了大李在第一次喝酒通過檢查，第二次喝同樣的酒且經(jīng)過更長的時(shí)間檢查卻被定為飲酒駕車的情況，因?yàn)榈诙魏染茣r(shí)有第一次喝酒的殘留量。第四十七頁，共五十五頁，2022年，8月28日求微分方程（組）解析解的命令:dsolve(‘方程1’,‘方程2’,…,‘方程n’,‘初始條件’,‘自變量’)ToMATLAB（ff1）結(jié)果：u=tan(t-c)五、微分方程的MATLAB求解第四十八頁，共五十五頁，2022年，8月28日解輸入命令:y=dsolve('D2y+4*Dy+29*y=0','y(0)=0,Dy(0)=15','x')結(jié)果為:y=3e-2xsin（5x）ToMATLAB（ff2）第四十九頁，共五十五頁，2022年，8月28日解輸入命令：

[x,y,z]=dsolve('Dx=2*x-3*y+3*z','Dy=4*x-5*y+3*z','Dz=4*x-4*y+2*z','t')；x=simple(x)%將x化簡y=simple(y)z=simple(z)結(jié)果為：x=(c1-c2+c3+c2e-3t-c3e-3t)e2t

y=-c1e-4t+c2e-4t+c2e-3t-c3e-3t+c1-c2+c3)e2t

z=(-c1e-4t+c2e-4t+c1-c2+c3)e2t

ToMATLAB（ff3）返回第五十頁，共五十五頁，2022年，8月28日微分方程的數(shù)值解（一）常微分方程數(shù)值解的定義在生產(chǎn)和科研中所處理的微分方程往往很復(fù)雜,且大多得不出一般解．而實(shí)際中的對(duì)初值問題，一般是要求得到解在若干個(gè)點(diǎn)上滿足規(guī)定精確度的近似值，或者得到一個(gè)滿足精確度要求的便于計(jì)算的表達(dá)式．因此，研究常微分方程的數(shù)值解法是十分必要的．返回第五十一頁，共五十五頁，2022年，8月28日（二）用MATLAB軟件求常微分方程的數(shù)值解[t，x]=solver(’f’,ts,x0,options)ode45ode23ode113ode15sode23s由待解方程寫成的M文件名ts=[t0,tf]，t0、tf為自變量的初值和終值函數(shù)的初值ode23：組合的2/3階龍