数学建模微分方程模型讲稿.ppt

《数学建模微分方程模型讲稿.ppt》由会员分享，可在线阅读，更多相关《数学建模微分方程模型讲稿.ppt（42页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、关于数学建模微分方程模型第一页，讲稿共四十二页哦动态动态模型模型 描述对象特征随时间描述对象特征随时间(空间空间)的演变过程的演变过程.分析对象特征的变化规律分析对象特征的变化规律.预报对象特征的未来性态预报对象特征的未来性态.研究控制对象特征的手段研究控制对象特征的手段.根据函数及其变化率之间的关系确定函数根据函数及其变化率之间的关系确定函数.微分微分方程方程建模建模 根据建模目的和问题分析作出简化假设根据建模目的和问题分析作出简化假设.按照内在规律或用类比法建立微分方程按照内在规律或用类比法建立微分方程.第二页，讲稿共四十二页哦1 目标跟踪问题目标跟踪问题 设位于坐标原点的甲舰向位于设位于

2、坐标原点的甲舰向位于x轴上点轴上点A(1,0)处的乙舰发射处的乙舰发射导弹，导弹头始终对准乙舰如果乙舰以最大的速度导弹，导弹头始终对准乙舰如果乙舰以最大的速度v0(常数常数)沿平沿平行于行于y轴的直线行驶，导弹的速度是轴的直线行驶，导弹的速度是5v0，求导弹运行的曲线方程，求导弹运行的曲线方程乙舰行驶多远时，导弹将它击中？乙舰行驶多远时，导弹将它击中？第三页，讲稿共四十二页哦由(1),(2)消去t,整理得模型:第四页，讲稿共四十二页哦解法二解法二(数值解法)1建立M文件eq1m function dy=eq1(x,y)dy=zeros(2,1);dy(1)=y(2);dy(2)=1/5*sqr

3、t(1+y(1)2)/(1-x);2 取x0=0，xf=09999，建立主程序如下：x0=0，xf=09999 x,y=ode15s(eq1,x0 xf,0 0);plot(x,y(:,1),b)hold on y=0:001:2;plot(1,y,b*)结论结论:导弹大致在（导弹大致在（1，02）处击中乙舰）处击中乙舰.令y1=y,y2=y1，将方程（3）化为一阶微分方程组第五页，讲稿共四十二页哦结果见图导弹大致在（1，02）处击中乙舰，与前面的结论一致返 回 结论：时刻结论：时刻t=021时，导弹在（时，导弹在（1，021）处击中乙舰）处击中乙舰第六页，讲稿共四十二页哦背景背景 年份年份

4、1625 1830 1930 1960 1974 1987 1999人口人口(亿亿)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研究人口变化规律研究人口变化规律控制人口过快增长控制人口过快增长2 如何预报人口的增长如何预报人口的增长做出较准确的预报做出较准确的预报 建立人口数学模型建立人口数学模型 第七页，讲稿共四十二页哦指数增长模型指数增长模型马尔塞斯马尔塞斯1

5、798年年提出提出常用的计算公式常用的计算公式x(t)时刻时刻t的的人口人口基本假设基本假设:人口人口(相对相对)增长率增长率 r 是常数是常数今年人口今年人口 x0,年增长率年增长率 rk年后人口年后人口随着时间增加，人口按指数规律无限增长随着时间增加，人口按指数规律无限增长.与常用公式的一致与常用公式的一致?第八页，讲稿共四十二页哦指数增长模型的应用及局限性指数增长模型的应用及局限性 与与19世纪以前欧洲一些地区人口统计数据吻合世纪以前欧洲一些地区人口统计数据吻合.适用于适用于19世纪后迁往加拿大的欧洲移民后代世纪后迁往加拿大的欧洲移民后代.可用于短期人口增长预测可用于短期人口增长预测.不

6、符合不符合1919世纪后多数地区人口增长规律世纪后多数地区人口增长规律.不能预测较长期的人口增长过程不能预测较长期的人口增长过程.1919世纪后人口数据世纪后人口数据人口增长率人口增长率r不是常数不是常数(逐渐下降逐渐下降)第九页，讲稿共四十二页哦阻滞增长模型阻滞增长模型Logistic 模型模型人口增长到一定数量后，增长率下降的原因：人口增长到一定数量后，增长率下降的原因：资源、环境等因素对人口增长的阻滞作用资源、环境等因素对人口增长的阻滞作用,且阻滞作用随人口数量增加而变大且阻滞作用随人口数量增加而变大假设假设r固有增长率固有增长率(x很小时很小时)xm人口容量（资源、环境能容纳的最大数量

7、）人口容量（资源、环境能容纳的最大数量）r是是x的减函数的减函数第十页，讲稿共四十二页哦dx/dtxOxmxm/2txOx增加先快后慢增加先快后慢xmx0 xm/2阻滞增长模型阻滞增长模型(Logistic模型模型)指数增长指数增长模型模型Logistic 模型的应用模型的应用 经济领域中的增长规律经济领域中的增长规律(耐用消费品的售量耐用消费品的售量).).种群数量模型种群数量模型(鱼塘中的鱼群鱼塘中的鱼群,森林中的树木森林中的树木).S形曲线形曲线第十一页，讲稿共四十二页哦参数估计参数估计用指数增长模型或阻滞增长模型作人口预报，用指数增长模型或阻滞增长模型作人口预报，必须先估计模型参数必须

8、先估计模型参数 r 或或 r,xm.模型的参数估计、检验和预报模型的参数估计、检验和预报 指数增长模型指数增长模型阻滞增长模型阻滞增长模型由统计数据用由统计数据用线性最小二乘法线性最小二乘法作参数估计作参数估计例：美国人口数据例：美国人口数据(百万百万)t 1860 1870 1880 1960 1970 1980 1990 2000 x 31.4 38.6 50.2 179.3 204.0 226.5 251.4 281.4 第十二页，讲稿共四十二页哦r=0.2022/10年，x0=6.0450 模型的参数估计、检验和预报模型的参数估计、检验和预报 指数增长模型指数增长模型阻滞增长模型阻滞增

9、长模型r=0.2557/10年，xm=392.0886 年年实际实际人人口口计计算人口算人口(指数增指数增长长模型模型)计计算人口算人口(阻滞增阻滞增长长模型模型)17903.96.03.918005.37.45.01960179.3188.0171.31970204.0230.1196.21980226.5281.7221.21990251.4344.8245.32000422.1指数增长模型指数增长模型阻滞增长模型阻滞增长模型第十三页，讲稿共四十二页哦用模型计算用模型计算2000年美国人口年美国人口误差约误差约2.5%与实际数据比较与实际数据比较(2000年年281.4)=274.5模型的

10、参数估计、检验和预报模型的参数估计、检验和预报 为作为作模型检验模型检验在参数估计时未用在参数估计时未用2000年实际数据年实际数据加入加入2000年数据重估模型参数年数据重估模型参数r=0.2490，xm=434.0 x(2010)=306.0 预报预报美国美国2010年人口年人口 美国人口普查局美国人口普查局2010年年12月月21日公布：截止到日公布：截止到2010年年4月月1日美国总人日美国总人口为口为3.087亿亿.预报误差不到预报误差不到1%！第十四页，讲稿共四十二页哦场景场景3 如何施救药物中毒如何施救药物中毒两位家长带着孩子急匆匆来到医院急诊室两位家长带着孩子急匆匆来到医院急诊

11、室.诉说两小时前孩子一次误吞下诉说两小时前孩子一次误吞下11片片治疗哮喘病、剂量治疗哮喘病、剂量100mg/片片的氨的氨茶碱片，已出现呕吐、头晕等不良症状茶碱片，已出现呕吐、头晕等不良症状.按照药品使用说明书，氨茶碱的每次用量成人是按照药品使用说明书，氨茶碱的每次用量成人是100200mg，儿童是，儿童是35 mg/kg.过量服用可使血药浓度过量服用可使血药浓度(单位血液容积中的药量单位血液容积中的药量)过高，过高，100g/ml浓度会出浓度会出现现严重中毒严重中毒,200g/ml浓度可致命浓度可致命.医生需要判断：孩子的血药浓度会不会达到医生需要判断：孩子的血药浓度会不会达到100200 g

12、/ml；如果会；如果会达到，应采取怎样的达到，应采取怎样的紧急施救紧急施救方案方案.第十五页，讲稿共四十二页哦调查与分析调查与分析转移率正转移率正比于比于x排除率正排除率正比于比于y胃肠道胃肠道血液系统血液系统口服药物口服药物体外体外认为血液系统内药物的分布，即血药浓度是均匀的，可以将血液系认为血液系统内药物的分布，即血药浓度是均匀的，可以将血液系统看作一个房室，建立统看作一个房室，建立“一室模型一室模型”.药量药量x(t)药量药量y(t)血液系统对药物的吸收率血液系统对药物的吸收率(胃肠道到血液系统的转移率胃肠道到血液系统的转移率)和排除率可和排除率可以由以由半衰期半衰期确定确定.半衰期半衰

13、期可以从药品说明书上查到可以从药品说明书上查到.第十六页，讲稿共四十二页哦通常，血液总量约为人体体重的通常，血液总量约为人体体重的7%8%，体重，体重5060 kg的成年人有的成年人有4000ml左右的血液左右的血液.目测这个孩子的体重约为成年人的一半，可认为其血目测这个孩子的体重约为成年人的一半，可认为其血液总量约为液总量约为2000ml.调查与分析调查与分析血药浓度血药浓度=药量药量/血液总量血液总量 口服活性炭来吸附药物，可使药物的排除率增加到原来（人体口服活性炭来吸附药物，可使药物的排除率增加到原来（人体自身）的自身）的2倍倍.临床施救的办法：临床施救的办法：体外血液透析，药物排除率可

14、增加到原来的体外血液透析，药物排除率可增加到原来的6倍，但是安全倍，但是安全性不能得到充分保证性不能得到充分保证.第十七页，讲稿共四十二页哦模型假设模型假设 1.胃肠道中药物向血液的转移率与胃肠道中药物向血液的转移率与x(t)成正比，比例系数成正比，比例系数(0)，总剂量，总剂量1100 mg药物在药物在t=0瞬间进入胃肠道瞬间进入胃肠道.2.血液系统中药物的排除率与血液系统中药物的排除率与y(t)成正比，比例系数成正比，比例系数(0)，t=0时血液中无药物时血液中无药物.3.氨茶碱被吸收的半衰期为氨茶碱被吸收的半衰期为5 h，排除的半衰期为，排除的半衰期为6 h.4.孩子的血液总量为孩子的血

15、液总量为2000 ml.胃肠道中药量胃肠道中药量x(t),血液系统中药量血液系统中药量y(t)，时间，时间t以孩子以孩子误服药的时刻为起点（误服药的时刻为起点（t=0）.第十八页，讲稿共四十二页哦模型建立模型建立x(t)下降速度与下降速度与x(t)成正比成正比(比例系数比例系数),总剂量总剂量1100mg药物在药物在t=0瞬间进瞬间进入胃肠道入胃肠道.转移率转移率正比于正比于x排除率排除率正比于正比于y胃肠道胃肠道血液系统血液系统口服药物口服药物体外体外药量药量x(t)药量药量y(t)y(t)由吸收而增长的速度是由吸收而增长的速度是x，由排除而减少的速度与，由排除而减少的速度与y(t)成正比成

16、正比(比例比例系数系数),t=0时血液中无药物时血液中无药物.第十九页，讲稿共四十二页哦模型模型求解求解 药物吸收的半衰期为药物吸收的半衰期为5 h 药物排除的半衰期为药物排除的半衰期为6 h 只考虑血液对药物的排除只考虑血液对药物的排除第二十页，讲稿共四十二页哦血液总量血液总量2000ml血药浓度血药浓度200g/ml结果及分析结果及分析 胃肠道药量胃肠道药量血液系统药量血液系统药量血药浓度血药浓度100g/mly(t)=200mg严重中毒严重中毒y(t)=400mg致命致命t=1.62t=4.87t=7.89y=442孩子到达医院前已严重中毒，如不及时施救，约孩子到达医院前已严重中毒，如不

17、及时施救，约3h3h后将致命！后将致命！y(2)=236.5 第二十一页，讲稿共四十二页哦施救方案施救方案 口服活性炭使药物排除率口服活性炭使药物排除率增至原来的增至原来的2倍倍.孩子到达医院孩子到达医院(t=2)就开始施救，血液中药量记作就开始施救，血液中药量记作z(t)=0.1386(不变)，=0.11552=0.2310 第二十二页，讲稿共四十二页哦施救方案施救方案 t=5.26z=318 施救后血液中药量施救后血液中药量z(t)显著低于显著低于y(t).z(t)最大值低于致最大值低于致命水平命水平.要使要使z(t)在施救后在施救后立即下降，可算出立即下降，可算出至少应为至少应为0.48

18、85.若采用体外血液透析，若采用体外血液透析，可增至可增至0.11556=0.693，血液中，血液中药量下降更快；临床上是否需要采取这种办法，当由医药量下降更快；临床上是否需要采取这种办法，当由医生综合考虑并征求病人家属意见后确定生综合考虑并征求病人家属意见后确定.第二十三页，讲稿共四十二页哦偏微分方程与数学模型2022/10/15济南大学 数学科学学院24第二十四页，讲稿共四十二页哦偏微分方程偏微分方程偏微分方程偏微分方程（Partial Differential Equations）指在物理学、力学、工程技术以及其他自然科学、技术科学、指在物理学、力学、工程技术以及其他自然科学、技术科学、

19、管理科学、甚至社会科学等的研究中归纳出来的一些含有未管理科学、甚至社会科学等的研究中归纳出来的一些含有未知函数及其偏导数的方程知函数及其偏导数的方程2022/10/15济南大学 数学科学学院25第二十五页，讲稿共四十二页哦什么是偏微分方程？什么是偏微分方程？2022/10/15济南大学 数学科学学院26物理量物理量(如位移、温度等如位移、温度等)-)-时间、空间位置时间、空间位置 -物理量的变化规律物理量的变化规律(偏微分方程偏微分方程)第二十六页，讲稿共四十二页哦 例子 2022/10/15济南大学 数学科学学院27第二十七页，讲稿共四十二页哦研究内容研究内容 2022/10/15济南大学

20、数学科学学院28一般规律一般规律+定解条件定解条件(初始条件、边界条件初始条件、边界条件)定解问题定解问题 定解问题的适定性：定解问题的适定性：存在性（存在性（ExistenceExistence）唯一性（唯一性（UniquenessUniqueness）稳定性（稳定性（StabilityStability）+附加条件附加条件方程方程第二十八页，讲稿共四十二页哦4 人口预测和控制人口预测和控制 年龄分布对于人口预测的重要性年龄分布对于人口预测的重要性.只考虑自然出生与死亡，不计迁移只考虑自然出生与死亡，不计迁移.人口人口发展发展方程方程第二十九页，讲稿共四十二页哦人口发展方程人口发展方程一阶偏

21、微分方程一阶偏微分方程第三十页，讲稿共四十二页哦人口发展方程人口发展方程p0(r)已知函数已知函数(人口调查人口调查)f(t)生育率生育率(控制人口手段控制人口手段)0tr第三十一页，讲稿共四十二页哦生育率的分解生育率的分解 总和生育率总和生育率h生育模式生育模式0第三十二页，讲稿共四十二页哦人口发展方程和生育率人口发展方程和生育率总和生育率总和生育率控制生育的多少控制生育的多少生育模式生育模式控制生育的早晚和疏密控制生育的早晚和疏密 正反馈系统正反馈系统 滞后作用很大滞后作用很大第三十三页，讲稿共四十二页哦人口指数人口指数1）人口总数）人口总数2）平均年龄）平均年龄3）平均寿命）平均寿命t时

22、刻出生的人，死亡率按时刻出生的人，死亡率按 (r,t)计算的平均存活时间计算的平均存活时间4）老龄化指数）老龄化指数控制生育率控制生育率控制控制 N(t)不过大不过大控制控制 (t)不过高不过高第三十四页，讲稿共四十二页哦5 烟雾的扩散与消失烟雾的扩散与消失现象现象和和问题问题 炮弹在空中爆炸，烟雾向四周扩散，形成圆炮弹在空中爆炸，烟雾向四周扩散，形成圆形不透光区域形不透光区域.不透光区域不断扩大，然后区域边界逐不透光区域不断扩大，然后区域边界逐渐明亮，区域缩小，最后烟雾消失渐明亮，区域缩小，最后烟雾消失.建立模型描述烟雾扩散和消失过程，分建立模型描述烟雾扩散和消失过程，分析消失时间与各因素的

23、关系析消失时间与各因素的关系.问题问题分析分析 无穷空间由瞬时点源导致的扩散过程，用二无穷空间由瞬时点源导致的扩散过程，用二阶偏微分方程描述烟雾浓度的变化阶偏微分方程描述烟雾浓度的变化.观察到的烟雾消失与烟雾对光线的吸收、观察到的烟雾消失与烟雾对光线的吸收、以及仪器对明暗的灵敏程度有关以及仪器对明暗的灵敏程度有关.第三十五页，讲稿共四十二页哦模型模型假设假设1）烟雾在无穷空间扩散，不受大地和风的影）烟雾在无穷空间扩散，不受大地和风的影响；扩散服从热传导定律响；扩散服从热传导定律.2）光线穿过烟雾时光强的相对减少与烟雾浓度）光线穿过烟雾时光强的相对减少与烟雾浓度成正比；无烟雾的大气不影响光强成正

24、比；无烟雾的大气不影响光强.3）穿过烟雾进入仪器的光线只有明暗之分，明暗界）穿过烟雾进入仪器的光线只有明暗之分，明暗界限由仪器灵敏度决定限由仪器灵敏度决定.模型模型建立建立1）烟雾浓度）烟雾浓度 的变化规律的变化规律热传导定律：单位时间通过单位法向热传导定律：单位时间通过单位法向面积的流量与浓度梯度成正比面积的流量与浓度梯度成正比.第三十六页，讲稿共四十二页哦曲面积分曲面积分奥氏公式奥氏公式1）烟雾浓度）烟雾浓度 的变化规律的变化规律的微分形式，并利用积分中值定理的微分形式，并利用积分中值定理第三十七页，讲稿共四十二页哦 初始条件初始条件Q炮弹释放的烟雾总量炮弹释放的烟雾总量 单位强度的点源函

25、数单位强度的点源函数 对任意对任意t,C的等值面是球面的等值面是球面 x2+y2+z2=R2;RC 仅当仅当 t,对任意点对任意点(x,y,z),C01）烟雾浓度）烟雾浓度 的变化规律的变化规律第三十八页，讲稿共四十二页哦2）光强穿过烟雾时的变化规律）光强穿过烟雾时的变化规律假设假设2）光强的相对减少与烟雾浓度成正比）光强的相对减少与烟雾浓度成正比.I(l)沿沿l方向的光强，方向的光强，C(l)沿沿l方向的烟雾强度方向的烟雾强度记未进入烟雾记未进入烟雾(l l0)时光强为时光强为 I(l0)=I0第三十九页，讲稿共四十二页哦3）仪器灵敏度与烟雾明暗界限）仪器灵敏度与烟雾明暗界限烟雾浓度连续变化

26、烟雾浓度连续变化烟雾中光强连续变化烟雾中光强连续变化仪器仪器z-设光源在设光源在z=-,仪器在仪器在z=,则观测到的则观测到的明暗界限为明暗界限为不透光区域有扩大、不透光区域有扩大、缩小、消失的过程缩小、消失的过程穿过烟雾进入仪器的光线只有明暗之穿过烟雾进入仪器的光线只有明暗之分，明暗界限由仪器灵敏度决定分，明暗界限由仪器灵敏度决定.不透光区域边界不透光区域边界第四十页，讲稿共四十二页哦4）不透光区域边界的变化规律）不透光区域边界的变化规律对任意对任意t,不透光区域边界是圆周不透光区域边界是圆周不透光区域不透光区域边界半径边界半径第四十一页，讲稿共四十二页哦感感谢谢大大家家观观看看第四十二页，讲稿共四十二页哦