《微分方程模型》PPT课件

1、精选ppt1微分方程模型 马忠明 精选ppt2动态动态模型模型 描述对象特征随时间描述对象特征随时间(空间空间)的演变过程的演变过程 分析对象特征的变化规律分析对象特征的变化规律 预报对象特征的未来性态预报对象特征的未来性态 研究控制对象特征的手段研究控制对象特征的手段 根据函数及其变化率之间的关系确定函数根据函数及其变化率之间的关系确定函数微分微分方程方程建模建模 根据建模目的和问题分析作出简化假设根据建模目的和问题分析作出简化假设 按照内在规律或用类比法建立微分方程按照内在规律或用类比法建立微分方程精选ppt3背景背景 年年 1625 1830 1930 1960 1974 1987 19

2、99人口人口(亿亿) 5 10 20 30 40 50 60世界人口增长概况世界人口增长概况中国人口增长概况中国人口增长概况 年年 1908 1933 1953 1964 1982 1990 1995 2000人口人口(亿亿) 3.0 4.7 6.0 7.2 10.3 11.3 12.0 13.0研究人口变化规律研究人口变化规律控制人口过快增长控制人口过快增长一、一、 人口增长模型人口增长模型精选ppt4指数增长模型指数增长模型马尔萨斯提出马尔萨斯提出 ( (1798) )常用的计算公式常用的计算公式kkrxx)1 (0 x(t) 时刻时刻t的人口的人口基本假设基本假设 : 人口人口(相对相对

3、)增长率增长率 r 是常数是常数trtxtxttx)()()(今年人口今年人口 x0, 年增长率年增长率 rk年后人口年后人口0)0(,xxrxdtdxrtextx0)(trextx)()(0trx)1 (0随着时间增加，人口按指数规律无限增长随着时间增加，人口按指数规律无限增长精选ppt5指数增长模型的应用及局限性指数增长模型的应用及局限性 与与19世纪以前欧洲一些地区人口统计数据吻合世纪以前欧洲一些地区人口统计数据吻合 适用于适用于19世纪后迁往加拿大的欧洲移民后代世纪后迁往加拿大的欧洲移民后代 可用于短期人口增长预测可用于短期人口增长预测 不符合不符合19世纪后多数地区人口增长规律世纪后

4、多数地区人口增长规律 不能预测较长期的人口增长过程不能预测较长期的人口增长过程1919世纪后人口数据世纪后人口数据人口增长率人口增长率r r不是常数不是常数( (逐渐下降逐渐下降) )精选ppt6精选ppt7阻滞增长模型阻滞增长模型( (Logistic模型模型) ) 人口增长到一定数量后，增长率下降的原因：人口增长到一定数量后，增长率下降的原因：资源、环境等因素对人口增长的阻滞作用资源、环境等因素对人口增长的阻滞作用 且阻滞作用随人口数量增加而变大且阻滞作用随人口数量增加而变大假设假设) 0,()(srsxrxrr固有增长率固有增长率(x很小时很小时)xm人口容量（资源、环境能容纳的最大数量

5、）人口容量（资源、环境能容纳的最大数量）)1 ()(mxxrxrr是是x的减函数的减函数mxrs 0)(mxr精选ppt8rxdtdx)1 ()(mxxrxxxrdtdxdx/dtx0 xmxm/2xmx txxxemmrt( )()110tx0 x(t)S形曲线形曲线, x增加先快后慢增加先快后慢x0 xm/2阻滞增长模型阻滞增长模型( (Logistic模型模型) )精选ppt9参数估计参数估计用指数增长模型或阻滞增长模型作人口用指数增长模型或阻滞增长模型作人口预报，必须先估计模型参数预报，必须先估计模型参数 r 或或 r, xm 利用统计数据用最小二乘法作拟合利用统计数据用最小二乘法作拟

6、合例：美国人口数据（单位例：美国人口数据（单位百万）百万） 1860 1870 1880 1960 1970 1980 1990 31.4 38.6 50.2 179.3 204.0 226.5 251.4专家估计专家估计阻滞增长模型阻滞增长模型( (Logistic模型模型) )r=0.2557, xm=392.1精选ppt10模型检验模型检验用模型计算用模型计算2000年美国人口，与实际数据比较年美国人口，与实际数据比较/ )1990(1)1990()1990()1990()2000(mxxrxxxxx实际为实际为281.4 (百万百万)5 .274)2000(x模型应用模型应用预报美国预

7、报美国2010年的人口年的人口加入加入2000年人口数据后重新估计模型参数年人口数据后重新估计模型参数Logistic 模型在经济领域中的应用模型在经济领域中的应用( (如耐用消费品的售量如耐用消费品的售量) )阻滞增长模型阻滞增长模型( (Logistic模型模型) )r=0.2490, xm=434.0 x(2010)=306.0精选ppt11Leslie模型与人口发展模型模型与人口发展模型l中国人口增长预测中国人口增长预测论文论文1l中国人口增长预测中国人口增长预测论文论文2精选ppt12人口预测和控制人口预测和控制)(),(,0),0(tNtrFtFmrFtrp),( 年龄分布对于人口

8、预测的重要性年龄分布对于人口预测的重要性 只考虑自然出生与死亡，不计迁移只考虑自然出生与死亡，不计迁移人口人口发展发展方程方程的人口）年龄人口分布函数rtrF(),(人口密度函数),(trp人口总数)(tN最高年龄)(mr精选ppt13),(),(trptrtprp11,),(),(),(),(),(),(drdtdttrptrtrpdttrpdttrpdttdrrp人口发展方程人口发展方程死亡率),(trdrtrp),(人数年龄,drrrt死亡人数内),(dttt人数年龄,11drdrrdrrdtt1drdt 一阶偏微分方程一阶偏微分方程drdttrptr),(),(drdttdrrp),(

9、1精选ppt140),(),0(0),()0 ,(),(),(0ttftprrprptrptrtprp人口发展方程人口发展方程已知函数（人口调查）已知函数（人口调查）生育率（控制人口手段）生育率（控制人口手段）0tr)(0rprt )(tfrt rt )(),(rtrrtertfrtetrptrprrtrdssdss,)(0,)(),(0)()(0rdstsptrF0),(),(mrdstsptN0),()(精选ppt1521),(),(),()(rrdrtrptrktrbtf),()(),(trhttrb211),(rrdrtrh21),()(rrdrtrbt生育率的分解生育率的分解性别比函

10、数女性 )(),(trk生育数女性 )(),(trb育龄区间,21rr21),(),(),()()(rrdrtrptrktrhttf 总和生育率总和生育率h生育模式生育模式)(),(rhtrh01r2rr精选ppt16rtertfrtetrptrprrtrdssdss,)(0,)(),(0)()(021),(),(),()()(rrdrtrptrktrhttf人口发展方程和生育率人口发展方程和生育率)(t总和生育率总和生育率控制生育的多少控制生育的多少),(trh生育模式生育模式控制生育的早晚和疏密控制生育的早晚和疏密),(),(trptrtprp)(tf)(0rp),(trp)(t 正反馈系

11、统正反馈系统 滞后作用很大滞后作用很大精选ppt17mrdrtrrptNtR0),()(1)(tdrtrdetSt0),()()(/ )()(tStRt mrdrtrptN0),()(人口指数人口指数1）人口总数）人口总数2）平均年龄）平均年龄3）平均寿命）平均寿命t时刻出生的人，死亡率按时刻出生的人，死亡率按 (r,t) 计算的平均存活时间计算的平均存活时间4）老龄化指数）老龄化指数控制生育率控制生育率控制控制 N(t)不过大不过大控制控制 (t)不过高不过高精选ppt18二、二、 随机人口模型随机人口模型背景背景 一个人的出生和死亡是随机事件一个人的出生和死亡是随机事件一个国家或地区一个国

12、家或地区平均生育率平均生育率平均死亡率平均死亡率确定性模型确定性模型一个家族或村落一个家族或村落出生概率出生概率死亡概率死亡概率随机性模型随机性模型对象对象X(t) 时刻时刻 t 的人口的人口, 随机变量随机变量.Pn(t) 概率概率P(X(t)=n), n=0,1,2,研究研究Pn(t)的变化规律；得到的变化规律；得到X(t)的期望和方差的期望和方差精选ppt19若若X(t)=n, 对对t到到t+ t的出生和死亡概率作以下假设的出生和死亡概率作以下假设1)出生一人的概率与出生一人的概率与 t成正比，记成正比，记bn t ;出生二人及二人以上的概率为出生二人及二人以上的概率为o( t).2)死

13、亡一人的概率与死亡一人的概率与 t成正比，记成正比，记dn t ;死亡二人及二人以上的概率为死亡二人及二人以上的概率为o( t).3)出生和死亡是相互独立的随机事件。出生和死亡是相互独立的随机事件。 bn与与n成正比，记成正比，记bn= n , 出生概率出生概率；dn与与n成正比，记成正比，记dn= n， 死亡概率死亡概率。进一步假设进一步假设模型假设模型假设精选ppt20)()1)()()()(1111totdtbtPtdtPtbtPttPnnnnnnnn建模建模为得到为得到Pn(t) P(X(t)=n),的变化规律，的变化规律，考察考察Pn(t+ t) =P(X(t + t)=n).事件事

14、件X(t + t)=n的分解的分解X(t)=n-1, t内出生一人内出生一人X(t)=n+1, t内死亡一人内死亡一人X(t)=n, t内没有出生和死亡内没有出生和死亡其它其它(出生或死亡二人，出生或死亡二人，出生且死亡一人，出生且死亡一人， )概率概率Pn(t+ t) Pn-1(t), bn-1 t Pn+1(t), dn+1 t Pn(t), 1-bn t -dn t o( t)精选ppt21)()()() 1()() 1(11tnPtPntPndtdPnnnn)()()()(1111tPdbtPdtPbdtdPnnnnnnnn一组递推微分方程一组递推微分方程求解的困难和不必要求解的困难和

15、不必要00, 0, 1)0(nnnnPn(t=0时已知人口为时已知人口为n0)转而考察转而考察X(t)的期望和方差的期望和方差bn= n，dn= n微分方程微分方程建模建模精选ppt221)()()()(nntEtnPdtdE)()()()1()()1(121111tPntPnntPnndtdEnnnnnn1)()(nntnPtEX(t)的期望的期望求解求解)()()() 1()() 1(11tnPtPntPndtdPnnnn基本方程基本方程1)()1(kktPkkn-1=k1nndtdPndtdEn+1=k)()1(1tPkkkk精选ppt23求解求解0)0()()(nEtEdtdErtex

16、tx0)( 比较：确定性指数增长模型比较：确定性指数增长模型)()()(212tEtPntDnnX(t)的方差的方差E(t)- (t) - = r D(t) rentErt,)(0E(t)+ (t)Et0n0 , D(t) 1)()()(0tteentDX(t)大致在大致在 E(t) 2 (t) 范围内（范围内（ (t) 均方差）均方差）r 增长概率增长概率r 平均增长率平均增长率精选ppt24三、经济增长模型三、经济增长模型精选ppt25精选ppt26精选ppt27精选ppt28精选ppt29精选ppt30精选ppt31精选ppt32四、传染病模型四、传染病模型问题问题 描述传染病的传播过程

17、描述传染病的传播过程 分析受感染人数的变化规律分析受感染人数的变化规律 预报传染病高潮到来的时刻预报传染病高潮到来的时刻 预防传染病蔓延的手段预防传染病蔓延的手段 按照传播过程的一般规律，按照传播过程的一般规律，用机理分析方法建立模型用机理分析方法建立模型精选ppt33 已感染人数已感染人数 (病人病人) i(t) 每个病人每天有效接触每个病人每天有效接触(足以使人致病足以使人致病)人数为人数为 模型模型1 1假设假设ttititti)()()(若有效接触的是病人，若有效接触的是病人，则不能使病人数增加则不能使病人数增加必须区分已感染者必须区分已感染者(病病人人)和未感染者和未感染者(健康人健

18、康人)建模建模0)0(iiidtdiitteiti0)(?精选ppt34sidtdi1)()(tits模型模型2 2区分已感染者区分已感染者(病人病人)和未感染者和未感染者(健康人健康人)假设假设1）总人数）总人数N不变，病人和健康不变，病人和健康 人的人的 比例分别为比例分别为)(),(tsti 2）每个病人每天有效接触人数）每个病人每天有效接触人数为为 , 且且使接触的健康人致病使接触的健康人致病建模建模ttNitstittiN)()()()(0)0()1(iiiidtdi 日日接触率接触率SI 模型模型精选ppt35teiti1111)(00)0()1(iiiidtdi模型模型21/2t

19、mii010t11ln01itmtm传染病高潮到来时刻传染病高潮到来时刻 (日接触率日接触率) tm 1itLogistic 模型病人可以治愈！病人可以治愈！?t=tm, di/dt 最大最大精选ppt36模型模型3传染病无免疫性传染病无免疫性病人治愈成病人治愈成为健康人，健康人可再次被感染为健康人，健康人可再次被感染增加假设增加假设SIS 模型模型3）病人每天治愈的比例为）病人每天治愈的比例为 日日治愈率治愈率ttNittitNstittiN)()()()()(建模建模/ 日接触率日接触率1/ 感染期感染期 一个感染期内一个感染期内每个病人的每个病人的有效接触人数，称为有效接触人数，称为接触数接触数。0)0()1(iiiiidtdi精选ppt371,01,11)(i)11 (iidtdi模型模型3i0i0接触数接触数 =1 阈值阈值/1)(ti形曲线增长按Sti )(感染期内感染期内有效接触感染的有效接触感染的健康者人数不超过病人数健康者人数不超过病人数小01i1-1/ i0iiidtdi)1 (模型模型2(SI模型模型)如何