常微分方程李淑娟

常微分方程李淑娟第1頁(yè)，共41頁(yè)，2023年，2月20日，星期一

函數(shù)，是反映客觀現(xiàn)實(shí)世界運(yùn)動(dòng)過程中量與量之間的一種關(guān)系。

現(xiàn)實(shí)中的量與量之間的關(guān)系稍微復(fù)雜一些，不能直接寫出它們之間的關(guān)系，但容易地建立這些變量和它們的導(dǎo)數(shù)（或微分）直接的關(guān)系式。數(shù)學(xué)建模工作室歡迎您第2頁(yè)，共41頁(yè)，2023年，2月20日，星期一表示未知函數(shù)、未知函數(shù)的導(dǎo)數(shù)(微分)以及自變量之間的關(guān)系的方程。數(shù)學(xué)建模工作室歡迎您未知函數(shù)是多元函數(shù)的微分方程未知函數(shù)是一元函數(shù)的微分方程。微分方程常微分方程偏微分方程第3頁(yè)，共41頁(yè)，2023年，2月20日，星期一引例：曲線方程

已知曲線上任意一點(diǎn)處切線的斜率等于該點(diǎn)橫坐標(biāo)3倍，且過點(diǎn)(-1,3)，求此曲線方程。數(shù)學(xué)建模工作室歡迎您第4頁(yè)，共41頁(yè)，2023年，2月20日，星期一（1）根據(jù)規(guī)律列方程利用數(shù)學(xué)、力學(xué)、物理、化學(xué)等學(xué)科中的定理或經(jīng)過實(shí)驗(yàn)檢驗(yàn)的規(guī)律等來建立微分方程模型。（2）微元分析法利用已知的定理與規(guī)律尋找微元之間的關(guān)系式，與第一種方法不同的是對(duì)微元而不是直接對(duì)函數(shù)及其導(dǎo)數(shù)應(yīng)用規(guī)律。建立微分方程模型的方法數(shù)學(xué)建模工作室歡迎您第5頁(yè)，共41頁(yè)，2023年，2月20日，星期一在生物、經(jīng)濟(jì)等學(xué)科的實(shí)際問題中，許多現(xiàn)象的規(guī)律性不很清楚，即使有所了解也是極其復(fù)雜的，建模時(shí)在不同的假設(shè)下去模擬實(shí)際的現(xiàn)象，建立能近似反映問題的微分方程，然后從數(shù)學(xué)上求解或分析所建方程及其解的性質(zhì)，再去同實(shí)際情況對(duì)比，檢驗(yàn)此模型能否刻畫、模擬某些實(shí)際現(xiàn)象。建立微分方程模型的方法（3）模擬近似法數(shù)學(xué)建模工作室歡迎您第6頁(yè)，共41頁(yè)，2023年，2月20日，星期一

微分方程是數(shù)學(xué)建模中用到的最廣泛的模型之一。它可以應(yīng)用于預(yù)測(cè)人口的走勢(shì)，種群的增長(zhǎng)，污染物的擴(kuò)散等實(shí)際問題的解決中。數(shù)學(xué)建模工作室歡迎您第7頁(yè)，共41頁(yè)，2023年，2月20日，星期一數(shù)學(xué)建模中常用的微分方程模型人口增長(zhǎng)模型1傳染病傳播模型2數(shù)學(xué)建模工作室歡迎您第8頁(yè)，共41頁(yè)，2023年，2月20日，星期一數(shù)學(xué)建模工作室歡迎您人口增長(zhǎng)模型種群的數(shù)量本應(yīng)取離散值，但由于種群數(shù)量一般較大，短時(shí)間內(nèi)改變的是少數(shù)個(gè)體，與整體數(shù)量相比，這種變化是很微小的。我們通常假定大規(guī)模種群的個(gè)體數(shù)量是時(shí)間的連續(xù)可微函數(shù)，由此引起的誤差將是十分微小的。第9頁(yè)，共41頁(yè)，2023年，2月20日，星期一馬爾薩斯(Malthus)模型人口增長(zhǎng)模型阻滯增長(zhǎng)(Logistic)模型數(shù)學(xué)建模工作室歡迎您第10頁(yè)，共41頁(yè)，2023年，2月20日，星期一數(shù)學(xué)建模工作室歡迎您模型1：馬爾薩斯（Malthus）模型第11頁(yè)，共41頁(yè)，2023年，2月20日，星期一數(shù)學(xué)建模工作室歡迎您等式兩邊同時(shí)除以，有由初始條件，即為初始時(shí)刻的人口數(shù)再運(yùn)用極限的思想，令有

故解方程得第12頁(yè)，共41頁(yè)，2023年，2月20日，星期一數(shù)學(xué)建模工作室歡迎您當(dāng)r>0,人口將以指數(shù)規(guī)律增長(zhǎng)。當(dāng)r<0,人口將以指數(shù)規(guī)律減少。當(dāng)r=0,人口將保持常數(shù)。第13頁(yè)，共41頁(yè)，2023年，2月20日，星期一數(shù)學(xué)建模工作室歡迎您

馬爾薩斯模型的一個(gè)顯著特點(diǎn)：

種群數(shù)量翻一番所需的時(shí)間是固定的第14頁(yè)，共41頁(yè)，2023年，2月20日，星期一數(shù)學(xué)建模工作室歡迎您幾何級(jí)數(shù)增長(zhǎng)馬爾薩斯模型人口預(yù)測(cè)圖第15頁(yè)，共41頁(yè)，2023年，2月20日，星期一數(shù)學(xué)建模工作室歡迎您模型檢驗(yàn)人口增長(zhǎng)的實(shí)際情況與馬爾薩斯模型的預(yù)報(bào)結(jié)果基本相符。第16頁(yè)，共41頁(yè)，2023年，2月20日，星期一數(shù)學(xué)建模工作室歡迎您第17頁(yè)，共41頁(yè)，2023年，2月20日，星期一數(shù)學(xué)建模工作室歡迎您Malthus模型實(shí)際上只有在群體總數(shù)不太大時(shí)才合理，到總數(shù)增大時(shí)，生物群體的各成員之間由于有限的生存空間，有限的自然資源及食物等原因，就可能發(fā)生生存競(jìng)爭(zhēng)等現(xiàn)象。所以Malthus模型假設(shè)的人口凈增長(zhǎng)率不可能始終保持常數(shù)，它應(yīng)當(dāng)與人口數(shù)量有關(guān)。得出結(jié)論：Malthus模型只適用于短時(shí)期預(yù)測(cè)第18頁(yè)，共41頁(yè)，2023年，2月20日，星期一數(shù)學(xué)建模工作室歡迎您模型2：阻滯增長(zhǎng)（Logistic）模型

人口凈增長(zhǎng)率應(yīng)與人口數(shù)量有關(guān)，即反應(yīng)了自然因素對(duì)人口增長(zhǎng)的影響，令r=r(N)

第19頁(yè)，共41頁(yè)，2023年，2月20日，星期一數(shù)學(xué)建模工作室歡迎您1r(N)是未知函數(shù)，但根據(jù)實(shí)際背景，它無(wú)法用擬合方法來求.2為了得到實(shí)際意義的模型，我們不妨采用一下工程師原則.工程師們?cè)诮?shí)際問題的數(shù)學(xué)模型時(shí)，總是采用盡可能簡(jiǎn)單的方法.3r(N)最簡(jiǎn)單的形式是常數(shù)，此時(shí)得到Malthus模型。對(duì)Malthus模型的最簡(jiǎn)單的改進(jìn)就是引入一次項(xiàng)（競(jìng)爭(zhēng)項(xiàng)）.分析：第20頁(yè)，共41頁(yè)，2023年，2月20日，星期一數(shù)學(xué)建模工作室歡迎您注：設(shè)環(huán)境能供養(yǎng)的種群數(shù)量的上界為K（近似地將K看成常數(shù)），N表示當(dāng)前的種群數(shù)量。故第21頁(yè)，共41頁(yè)，2023年，2月20日，星期一數(shù)學(xué)建模工作室歡迎您故滿足初始條件N(0)=N0的解為：易見：

N(0)=N0

，第22頁(yè)，共41頁(yè)，2023年，2月20日，星期一數(shù)學(xué)建模工作室歡迎您不同初始條件下的N(t)的圖形第23頁(yè)，共41頁(yè)，2023年，2月20日，星期一數(shù)學(xué)建模工作室歡迎您大量實(shí)驗(yàn)資料表明用Logistic模型來描述種群的增長(zhǎng)規(guī)律，效果還是相當(dāng)不錯(cuò)的。例如，1945年克朗皮克（Crombic）做了一個(gè)人工飼養(yǎng)小谷蟲的實(shí)驗(yàn)；數(shù)學(xué)生物學(xué)家高斯（E·F·Gauss）也做了一個(gè)原生物草履蟲實(shí)驗(yàn)，實(shí)驗(yàn)結(jié)果都和Logistic曲線十分吻合。

模型檢驗(yàn)第24頁(yè)，共41頁(yè)，2023年，2月20日，星期一數(shù)學(xué)建模工作室歡迎您高斯把5只草履蟲放進(jìn)一個(gè)盛有0.5cm3營(yíng)養(yǎng)液的小試管，他發(fā)現(xiàn)，開始時(shí)草履蟲以每天230.9%的速率增長(zhǎng)，此后增長(zhǎng)速度不斷減慢，到第五天達(dá)到最大量375個(gè)，實(shí)驗(yàn)數(shù)據(jù)與r=2.309，N(0)=5的Logistic曲線：

幾乎完全吻合，見圖3.6。

圖3-6第25頁(yè)，共41頁(yè)，2023年，2月20日，星期一

Malthus模型和Logistic模型均為對(duì)微分方程所作的模擬近似方程。前一模型假設(shè)了種群增長(zhǎng)率r為一常數(shù)。后一模型則假設(shè)環(huán)境只能供養(yǎng)一定數(shù)量的種群，從而引入了一個(gè)競(jìng)爭(zhēng)項(xiàng)。

Malthus模型呈現(xiàn)的是J型增長(zhǎng)，只適應(yīng)于短期內(nèi)，并無(wú)外界因素影響。而Logistic模型呈現(xiàn)S型，適應(yīng)于中長(zhǎng)期且有外界因素影響。Malthus模型和Logistic模型的總結(jié)數(shù)學(xué)建模工作室歡迎您用模擬近似法建立微分方程來研究實(shí)際問題時(shí)必須對(duì)求得的解進(jìn)行檢驗(yàn)，看其是否與實(shí)際情況相符或基本相符。相符性越好則模擬得越好。否則，就得找出不相符的主要原因，對(duì)模型進(jìn)行修改。第26頁(yè)，共41頁(yè)，2023年，2月20日，星期一數(shù)學(xué)建模工作室歡迎您

Malthus模型與Logistic模型雖然都是為了研究種群數(shù)量的增長(zhǎng)情況而建立的，但它們也可用來研究其他實(shí)際問題，只要這些實(shí)際問題的數(shù)學(xué)模型有相同的微分方程即可。Malthus模型和Logistic模型的推廣如：?jiǎn)挝蝗藛T管理問題漁業(yè)管理問題單位資金管理問題森林管理問題第27頁(yè)，共41頁(yè)，2023年，2月20日，星期一數(shù)學(xué)建模工作室歡迎您單位人員管理問題合理安排進(jìn)人速度和出人速度，使得單位人員的利用率達(dá)到最高。如：?jiǎn)挝蝗藛T管理問題漁業(yè)管理問題單位資金管理問題森林管理問題第28頁(yè)，共41頁(yè)，2023年，2月20日，星期一數(shù)學(xué)建模工作室歡迎您

單位資金管理問題當(dāng)收入資金速率一定時(shí)，合理安排支出，使得在某段時(shí)間內(nèi)資金積累達(dá)到所需要求。如：?jiǎn)挝蝗藛T管理問題漁業(yè)管理問題單位資金管理問題森林管理問題第29頁(yè)，共41頁(yè)，2023年，2月20日，星期一數(shù)學(xué)建模工作室歡迎您森林管理問題主要協(xié)調(diào)植樹和用材的關(guān)系，使得森林發(fā)揮其應(yīng)有的作用。如：?jiǎn)挝蝗藛T管理問題漁業(yè)管理問題單位資金管理問題森林管理問題第30頁(yè)，共41頁(yè)，2023年，2月20日，星期一數(shù)學(xué)建模工作室歡迎您

漁業(yè)管理問題每年捕撈的速率控制在多少時(shí)，既能保持持續(xù)發(fā)展，還能有較大的收獲量。如：?jiǎn)挝蝗藛T管理問題漁業(yè)管理問題單位資金管理問題森林管理問題第31頁(yè)，共41頁(yè)，2023年，2月20日，星期一數(shù)學(xué)建模工作室歡迎您

漁業(yè)管理問題每年捕撈的速率控制在多少時(shí)，既能保持持續(xù)發(fā)展，還能有較大的收獲量。如：?jiǎn)挝蝗藛T管理問題漁業(yè)管理問題單位資金管理問題森林管理問題傳染病傳播第32頁(yè)，共41頁(yè)，2023年，2月20日，星期一數(shù)學(xué)建模工作室歡迎您傳染病經(jīng)常在世界各地流行，如霍亂，天花，艾滋病，SARS,H5N1病毒等，建立傳染病的數(shù)學(xué)模型,分析其變化規(guī)律，防止其蔓延是一項(xiàng)艱巨的任務(wù).一個(gè)地區(qū)有m個(gè)人，一名群眾不慎患傳染病，t小時(shí)后有n人發(fā)病，由于此地區(qū)不能及時(shí)隔離，問經(jīng)過t1小時(shí)、t2小時(shí)，患此傳染病的人數(shù)有多少？第33頁(yè)，共41頁(yè)，2023年，2月20日，星期一數(shù)學(xué)建模工作室歡迎您微分方程模型可以解決問題描述傳染病的傳播過程分析受感染人數(shù)的變化規(guī)律預(yù)報(bào)傳染病高潮到來的時(shí)刻第34頁(yè)，共41頁(yè)，2023年，2月20日，星期一數(shù)學(xué)建模工作室歡迎您傳染病模型每個(gè)病人每天有效接觸(足以使人致病)人數(shù)為模型1假設(shè)建模t時(shí)刻已感染人數(shù)(病人)

為i(t)第35頁(yè)，共41頁(yè)，2023年，2月20日，星期一數(shù)學(xué)建模工作室歡迎您若有效接觸的是病人，則不能使病人數(shù)增加必須區(qū)分已感染者(病人)和未感染者(健康人)?等式兩邊同時(shí)除以，有由初始條件，即為初始時(shí)刻的已感染疾病的人口數(shù).再運(yùn)用極限的思想，令有

故解方程得第36頁(yè)，共41頁(yè)，2023年，2月20日，星期一數(shù)學(xué)建模工作室歡迎您模型2區(qū)分已感染者(病人)和未感染者(健康人)假設(shè)1）總?cè)藬?shù)N不變，t時(shí)刻病人和健康人的比例分別為2）每個(gè)病人每天有效接觸(使接觸的健康人致病)人數(shù)為建模稱為SI模型為日接觸率即為L(zhǎng)ogistic模型第37頁(yè)，共41頁(yè)，2023年，2月20日，星期一數(shù)學(xué)建模工作室歡迎您(日接觸率)tm病人可以治愈！?t=tm,di/dt最大求解得

tm表示傳染病高潮到來時(shí)刻

1/2tmii010t第38頁(yè)，共41頁(yè)，2023年，2月20日，星期一數(shù)學(xué)建模工作室歡迎您微分方程可以解決的相關(guān)問題經(jīng)濟(jì)增長(zhǎng)的預(yù)測(cè)正規(guī)戰(zhàn)與游擊戰(zhàn)地震破壞程度的估計(jì)種群的相互競(jìng)爭(zhēng)藥物在體內(nèi)的分布與排除香煙過濾嘴煙霧的擴(kuò)散與消失萬(wàn)有引力定律的發(fā)現(xiàn)火箭的發(fā)射范.梅格倫(VanMeegren)偽造名畫案第39頁(yè)，共41頁(yè)，2023年，2月20日，星期一數(shù)學(xué)建模工作室期待您的加入謝謝?。?！第40頁(yè)，共41頁(yè)，2023年，2