數(shù)學(xué)建模微分方程模型市公開課一等獎(jiǎng)省名師優(yōu)質(zhì)課賽課一等獎(jiǎng)?wù)n件

微分方程模型2怎樣預(yù)報(bào)人口增加3

怎樣施救藥品中毒4

人口預(yù)測(cè)和控制模型1目標(biāo)跟蹤問題動(dòng)態(tài)模型

描述對(duì)象特征隨時(shí)間(空間)演變過程.

分析對(duì)象特征改變規(guī)律.

預(yù)報(bào)對(duì)象特征未來性態(tài).

研究控制對(duì)象特征伎倆.

依據(jù)函數(shù)及其改變率之間關(guān)系確定函數(shù).微分方程建模根據(jù)建模目和問題分析作出簡(jiǎn)化假設(shè).

按照內(nèi)在規(guī)律或用類比法建立微分方程.1目標(biāo)跟蹤問題

設(shè)位于坐標(biāo)原點(diǎn)甲艦向位于x軸上點(diǎn)A(1,0)處乙艦發(fā)射導(dǎo)彈，導(dǎo)彈頭一直對(duì)準(zhǔn)乙艦．假如乙艦以最大速度v0(常數(shù))沿平行于y軸直線行駛，導(dǎo)彈速度是5v0，求導(dǎo)彈運(yùn)行曲線方程．乙艦行駛多遠(yuǎn)時(shí)，導(dǎo)彈將它擊中？由(1),(2)消去t,整理得模型:解法二(數(shù)值解法)1．建立M文件eq1．m

functiondy=eq1(x,y)dy=zeros(2,1);dy(1)=y(2);dy(2)=1/5*sqrt(1+y(1)^2)/(1-x);2．取x0=0，xf=0．9999，建立主程序以下：

x0=0，xf=0．9999[x,y]=ode15s('eq1',[x0xf],[00]);plot(x,y(:,1),’b．')holdony=0:0．01:2;plot(1,y,’b*')

結(jié)論:導(dǎo)彈大致在（1，0．2）處擊中乙艦.令y1=y,y2=y1’，將方程（3）化為一階微分方程組．結(jié)果見圖導(dǎo)彈大致在（1，0．2）處擊中乙艦，與前面結(jié)論一致．返回

結(jié)論：時(shí)刻t=0．21時(shí)，導(dǎo)彈在（1，0．21）處擊中乙艦．背景

年份1625183019301960197419871999人口(億)5102030405060世界人口增加概況中國(guó)人口增加概況

年份1908193319531964198219901995人口(億)3.04.76.07.210.311.312.013.0研究人口改變規(guī)律控制人口過快增加2怎樣預(yù)報(bào)人口增加做出較準(zhǔn)確預(yù)報(bào)建立人口數(shù)學(xué)模型指數(shù)增加模型——馬爾塞斯1798年提出慣用計(jì)算公式x(t)~時(shí)刻t人口基本假設(shè)

:人口(相對(duì))增加率r

是常數(shù)今年人口x0,年增加率rk年后人口伴隨時(shí)間增加，人口按指數(shù)規(guī)律無(wú)限增加.與慣用公式一致?指數(shù)增加模型應(yīng)用及不足

與19世紀(jì)以前歐洲一些地域人口統(tǒng)計(jì)數(shù)據(jù)吻合.

適合用于19世紀(jì)后遷往加拿大歐洲移民后代.

可用于短期人口增加預(yù)測(cè).

不符合19世紀(jì)后多數(shù)地域人口增加規(guī)律.

不能預(yù)測(cè)較長(zhǎng)久人口增加過程.19世紀(jì)后人口數(shù)據(jù)人口增加率r不是常數(shù)(逐步下降)阻滯增加模型——Logistic模型人口增加到一定數(shù)量后，增加率下降原因：資源、環(huán)境等原因?qū)θ丝谠黾幼铚饔?且阻滯作用隨人口數(shù)量增加而變大假設(shè)r~固有增加率(x很小時(shí))xm~人口容量（資源、環(huán)境能容納最大數(shù)量）r是x減函數(shù)dx/dtxOxmxm/2txOx增加先快后慢xmx0xm/2阻滯增加模型(Logistic模型)指數(shù)增加模型Logistic模型應(yīng)用

經(jīng)濟(jì)領(lǐng)域中增加規(guī)律(耐用消費(fèi)品售量).

種群數(shù)量模型(魚塘中魚群,森林中樹木).S形曲線參數(shù)預(yù)計(jì)用指數(shù)增加模型或阻滯增加模型作人口預(yù)報(bào)，必須先預(yù)計(jì)模型參數(shù)r或r,xm.模型參數(shù)預(yù)計(jì)、檢驗(yàn)和預(yù)報(bào)

指數(shù)增加模型阻滯增加模型由統(tǒng)計(jì)數(shù)據(jù)用線性最小二乘法作參數(shù)預(yù)計(jì)例：美國(guó)人口數(shù)據(jù)(百萬(wàn))t186018701880…1960197019801990x31.438.650.2…179.3204.0226.5251.4281.4r=0.2022/，x0=6.0450模型參數(shù)預(yù)計(jì)、檢驗(yàn)和預(yù)報(bào)

指數(shù)增加模型阻滯增加模型r=0.2557/，xm=392.0886年實(shí)際人口計(jì)算人口(指數(shù)增加模型)計(jì)算人口(阻滯增加模型)17903.96.03.918005.37.45.0…………1960179.3188.0171.31970204.0230.1196.21980226.5281.7221.21990251.4344.8245.3422.1指數(shù)增加模型阻滯增加模型用模型計(jì)算美國(guó)人口誤差約2.5%與實(shí)際數(shù)據(jù)比較(281.4)=274.5模型參數(shù)預(yù)計(jì)、檢驗(yàn)和預(yù)報(bào)

為作模型檢驗(yàn)在參數(shù)預(yù)計(jì)時(shí)未用實(shí)際數(shù)據(jù)加入數(shù)據(jù)重估模型參數(shù)r=0.2490，xm=434.0x()=306.0預(yù)報(bào)美國(guó)人口美國(guó)人口普查局12月21日公布：截止到4月1日美國(guó)總?cè)丝跒?.087億.預(yù)報(bào)誤差不到1%！場(chǎng)景3怎樣施救藥品中毒兩位家長(zhǎng)帶著孩子急急忙來到醫(yī)院急診室.訴說兩小時(shí)前孩子一次誤吞下11片治療哮喘病、劑量100mg/片氨茶堿片，已出現(xiàn)嘔吐、頭暈等不良癥狀.按照藥品使用說明書，氨茶堿每次用量成人是100~200mg，兒童是3~5mg/kg.過量服用可使血藥濃度(單位血液容積中藥量)過高，100μg/ml濃度會(huì)出現(xiàn)嚴(yán)重中毒,200μg/ml濃度可致命.醫(yī)生需要判斷：孩子血藥濃度會(huì)不會(huì)到達(dá)100~200μg/ml；假如會(huì)到達(dá)，應(yīng)采取怎樣緊急施救方案.調(diào)查與分析轉(zhuǎn)移率正比于x排除率正比于y胃腸道血液系統(tǒng)口服藥品體外認(rèn)為血液系統(tǒng)內(nèi)藥品分布，即血藥濃度是均勻，能夠?qū)⒀合到y(tǒng)看作一個(gè)房室，建立“一室模型”

.藥量x(t)藥量y(t)血液系統(tǒng)對(duì)藥品吸收率(胃腸道到血液系統(tǒng)轉(zhuǎn)移率)和排除率能夠由半衰期確定.半衰期能夠從藥品說明書上查到.通常，血液總量約為人體體重7%

~8%，體重50~60kg成年人有4000ml左右血液.目測(cè)這個(gè)孩子體重約為成年人二分之一，可認(rèn)為其血液總量約為ml.調(diào)查與分析血藥濃度=藥量/血液總量

口服活性炭來吸附藥品，可使藥品排除率增加到原來（人體本身）2倍.臨床施救方法：

體外血液透析，藥品排除率可增加到原來6倍，不過安全性不能得到充分確保.模型假設(shè)

1.胃腸道中藥品向血液轉(zhuǎn)移率與x(t)成正比，百分比系數(shù)λ(>0)，總劑量1100mg藥品在t=0瞬間進(jìn)入胃腸道.2.血液系統(tǒng)中藥品排除率與y(t)成正比，百分比系數(shù)μ(>0)，t=0時(shí)血液中無(wú)藥品.3.氨茶堿被吸收半衰期為5h，排除半衰期為6h.4.孩子血液總量為ml.胃腸道中藥量x(t),血液系統(tǒng)中藥量y(t)，時(shí)間t以孩子誤服藥時(shí)刻為起點(diǎn)（t=0）.模型建立x(t)下降速度與x(t)成正比(百分比系數(shù)λ),總劑量1100mg藥品在t=0瞬間進(jìn)入胃腸道.轉(zhuǎn)移率正比于x排除率正比于y胃腸道血液系統(tǒng)口服藥品體外藥量x(t)藥量y(t)y(t)由吸收而增加速度是λx，由排除而降低速度與y(t)成正比(百分比系數(shù)μ),t=0時(shí)血液中無(wú)藥品.模型求解

藥品吸收半衰期為5h藥品排除半衰期為6h只考慮血液對(duì)藥品排除血液總量ml血藥濃度200μg/ml結(jié)果及分析胃腸道藥量血液系統(tǒng)藥量血藥濃度100μg/mly(t)=200mg嚴(yán)重中毒y(t)=400mg致命t=1.62t=4.87t=7.89y=442孩子抵達(dá)醫(yī)院前已嚴(yán)重中毒，如不及時(shí)施救，約3h后將致命！y(2)=236.5施救方案

口服活性炭使藥品排除率μ增至原來2倍.

孩子抵達(dá)醫(yī)院(t=2)就開始施救，血液中藥量記作z(t)λ=0.1386(不變)，μ=0.1155×2=0.2310施救方案

t=5.26z=318

施救后血液中藥量z(t)顯著低于y(t).z(t)最大值低于致命水平.

要使z(t)在施救后馬上下降，可算出μ最少應(yīng)為0.4885.若采取體外血液透析，μ可增至0.1155×6=0.693，血液中藥量下降更加快；臨床上是否需要采取這種方法，當(dāng)由醫(yī)生綜合考慮并征求病人家眷意見后確定.偏微分方程與數(shù)學(xué)模型/10/10濟(jì)南大學(xué)數(shù)學(xué)科學(xué)學(xué)院24偏微分方程偏微分方程（PartialDifferentialEquations）指在物理學(xué)、力學(xué)、工程技術(shù)以及其它自然科學(xué)、技術(shù)科學(xué)、管理科學(xué)、甚至社會(huì)科學(xué)等研究中歸納出來一些含有未知函數(shù)及其偏導(dǎo)數(shù)方程/10/10濟(jì)南大學(xué)數(shù)學(xué)科學(xué)學(xué)院25什么是偏微分方程？

/10/10濟(jì)南大學(xué)數(shù)學(xué)科學(xué)學(xué)院26物理量(如位移、溫度等)----時(shí)間、空間位置

---------------物理量改變規(guī)律(偏微分方程)例子

/10/10濟(jì)南大學(xué)數(shù)學(xué)科學(xué)學(xué)院27研究?jī)?nèi)容

/10/10濟(jì)南大學(xué)數(shù)學(xué)科學(xué)學(xué)院28普通規(guī)律+定解條件(初始條件、邊界條件)定解問題定解問題適定性：存在性（Existence）唯一性（Uniqueness）穩(wěn)定性（Stability）+附加條件方程4人口預(yù)測(cè)和控制年紀(jì)分布對(duì)于人口預(yù)測(cè)主要性.只考慮自然出生與死亡，不計(jì)遷移.人口發(fā)展方程人口發(fā)展方程一階偏微分方程人口發(fā)展方程p0(r)~已知函數(shù)(人口調(diào)查)f(t)~生育率(控制人口伎倆)0tr生育率分解~總和生育率h~生育模式0人口發(fā)展方程和生育率~總和生育率——控制生育多少~生育模式——控制生育早晚和疏密正反饋系統(tǒng)滯后作用很大人口指數(shù)1）人口總數(shù)2）平均年紀(jì)3）平均壽命t時(shí)刻出生人，死亡率按(r,t)計(jì)算平均存活時(shí)間4）老齡化指數(shù)控制生育率控制N(t)不過大控制(t)不過高5

煙霧擴(kuò)散與消失現(xiàn)象和問題炮彈在空中爆炸，煙霧向四面擴(kuò)散，形成圓形不透光區(qū)域.不透光區(qū)域不停擴(kuò)大，然后區(qū)域邊界逐步明亮，區(qū)域縮小，最終煙霧消失.建立模型描述煙霧擴(kuò)散和消失過程，分析消失時(shí)間與各原因關(guān)系.問題分析無(wú)窮空間由瞬時(shí)點(diǎn)源造成擴(kuò)散過程，用二階偏微分方程描述煙霧濃度改變.觀察到煙霧消失與煙霧對(duì)光線吸收、以及儀器對(duì)明暗靈敏程度相關(guān).模型假設(shè)1）煙霧在無(wú)窮空間擴(kuò)散，不受大地和風(fēng)影響；擴(kuò)散服從熱傳導(dǎo)定律.2）光線穿過煙霧時(shí)光強(qiáng)相對(duì)降低與煙霧濃度成正比；無(wú)煙霧大氣不影響光強(qiáng).3）穿過煙霧進(jìn)入儀器光線只有明暗之分，明暗界限由儀器靈敏度決定.模型建立1）煙霧濃度改變規(guī)律熱傳導(dǎo)定律：?jiǎn)挝粫r(shí)間經(jīng)過單位法向面積流量與濃度梯度成正比.