《数学建模微分方程模型PPT课件.ppt》由会员分享，可在线阅读，更多相关《数学建模微分方程模型PPT课件.ppt（42页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、关于数学建模微分方程模型第一张，PPT共四十二页，创作于2022年6月动态动态模型模型 描述对象特征随时间描述对象特征随时间(空间空间)的演变过程的演变过程.分析对象特征的变化规律分析对象特征的变化规律.预报对象特征的未来性态预报对象特征的未来性态.研究控制对象特征的手段研究控制对象特征的手段.根据函数及其变化率之间的关系确定函数根据函数及其变化率之间的关系确定函数.微分微分方程方程建模建模 根据建模目的和问题分析作出简化假设根据建模目的和问题分析作出简化假设.按照内在规律或用类比法建立微分方程按照内在规律或用类比法建立微分方程.第二张，PPT共四十二页，创作于2022年6月1 目标跟踪问题目
2、标跟踪问题 设位于坐标原点的甲舰向位于设位于坐标原点的甲舰向位于x轴上点轴上点A(1,0)处的乙舰发射导弹，处的乙舰发射导弹，导弹头始终对准乙舰如果乙舰以最大的速度导弹头始终对准乙舰如果乙舰以最大的速度v0(常数常数)沿平行于沿平行于y轴轴的直线行驶，导弹的速度是的直线行驶，导弹的速度是5v0，求导弹运行的曲线方程乙舰行驶，求导弹运行的曲线方程乙舰行驶多远时，导弹将它击中？多远时，导弹将它击中？第三张，PPT共四十二页，创作于2022年6月由(1),(2)消去t,整理得模型:第四张，PPT共四十二页，创作于2022年6月解法二解法二(数值解法)1建立M文件eq1m function dy=eq
3、1(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=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）化为一阶微分方程组第五张，PPT共四十二页，创作于2022年6月结果见图导弹大致在（1，02）处击中乙舰，与前面的结论一致返 回 结论：时刻结论：时刻t
4、=021时，导弹在（时，导弹在（1，021）处击中乙舰）处击中乙舰第六张，PPT共四十二页，创作于2022年6月背景背景 年份年份 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 如何预报人口的增长如何预报人口的增长做
5、出较准确的预报做出较准确的预报 建立人口数学模型建立人口数学模型 第七张，PPT共四十二页，创作于2022年6月指数增长模型指数增长模型马尔塞斯马尔塞斯1798年年提出提出常用的计算公式常用的计算公式x(t)时刻时刻t的的人口人口基本假设基本假设:人口人口(相对相对)增长率增长率 r 是常数是常数今年人口今年人口 x0,年增长率年增长率 rk年后人口年后人口随着时间增加，人口按指数规律无限增长随着时间增加，人口按指数规律无限增长.与常用公式的一致与常用公式的一致?第八张，PPT共四十二页，创作于2022年6月指数增长模型的应用及局限性指数增长模型的应用及局限性 与与19世纪以前欧洲一些地区人口
6、统计数据吻合世纪以前欧洲一些地区人口统计数据吻合.适用于适用于19世纪后迁往加拿大的欧洲移民后代世纪后迁往加拿大的欧洲移民后代.可用于短期人口增长预测可用于短期人口增长预测.不符合不符合1919世纪后多数地区人口增长规律世纪后多数地区人口增长规律.不能预测较长期的人口增长过程不能预测较长期的人口增长过程.1919世纪后人口数据世纪后人口数据人口增长率人口增长率r不是常数不是常数(逐渐下降逐渐下降)第九张，PPT共四十二页，创作于2022年6月阻滞增长模型阻滞增长模型Logistic 模型模型人口增长到一定数量后，增长率下降的原因：人口增长到一定数量后，增长率下降的原因：资源、环境等因素对人口增
7、长的阻滞作用资源、环境等因素对人口增长的阻滞作用,且阻滞作用随人口数量增加而变大且阻滞作用随人口数量增加而变大假设假设r固有增长率固有增长率(x很小时很小时)xm人口容量（资源、环境能容纳的最大数量）人口容量（资源、环境能容纳的最大数量）r是是x的减函数的减函数第十张，PPT共四十二页，创作于2022年6月dx/dtxOxmxm/2txOx增加先快后慢增加先快后慢xmx0 xm/2阻滞增长模型阻滞增长模型(Logistic模型模型)指数增长指数增长模型模型Logistic 模型的应用模型的应用 经济领域中的增长规律经济领域中的增长规律(耐用消费品的售量耐用消费品的售量).).种群数量模型种群数
8、量模型(鱼塘中的鱼群鱼塘中的鱼群,森林中的树木森林中的树木).S形曲线形曲线第十一张，PPT共四十二页，创作于2022年6月参数估计参数估计用指数增长模型或阻滞增长模型作人口预报，用指数增长模型或阻滞增长模型作人口预报，必须先估计模型参数必须先估计模型参数 r 或或 r,xm.模型的参数估计、检验和预报模型的参数估计、检验和预报 指数增长模型指数增长模型阻滞增长模型阻滞增长模型由统计数据用由统计数据用线性最小二乘法线性最小二乘法作参数估计作参数估计例：美国人口数据例：美国人口数据(百万百万)t 1860 1870 1880 1960 1970 1980 1990 2000 x 31.4 38.
9、6 50.2 179.3 204.0 226.5 251.4 281.4 第十二张，PPT共四十二页，创作于2022年6月r=0.2022/10年，x0=6.0450 模型的参数估计、检验和预报模型的参数估计、检验和预报 指数增长模型指数增长模型阻滞增长模型阻滞增长模型r=0.2557/10年，xm=392.0886 年年实际实际人人口口计计算人口算人口(指数增指数增长长模型模型)计计算人口算人口(阻滞增阻滞增长长模型模型)17903.96.03.918005.37.45.01960179.3188.0171.31970204.0230.1196.21980226.5281.7221.2199
10、0251.4344.8245.32000422.1指数增长模型指数增长模型阻滞增长模型阻滞增长模型第十三张，PPT共四十二页，创作于2022年6月用模型计算用模型计算2000年美国人口年美国人口误差约误差约2.5%与实际数据比较与实际数据比较(2000年年281.4)=274.5模型的参数估计、检验和预报模型的参数估计、检验和预报 为作为作模型检验模型检验在参数估计时未用在参数估计时未用2000年实际数据年实际数据加入加入2000年数据重估模型参数年数据重估模型参数r=0.2490，xm=434.0 x(2010)=306.0 预报预报美国美国2010年人口年人口 美国人口普查局美国人口普查局
11、2010年年12月月21日公布：截止到日公布：截止到2010年年4月月1日美国总日美国总人口为人口为3.087亿亿.预报误差不到预报误差不到1%！第十四张，PPT共四十二页，创作于2022年6月场景场景3 如何施救药物中毒如何施救药物中毒两位家长带着孩子急匆匆来到医院急诊室两位家长带着孩子急匆匆来到医院急诊室.诉说两小时前孩子一次误吞下诉说两小时前孩子一次误吞下11片片治疗哮喘病、剂量治疗哮喘病、剂量100mg/片片的氨茶碱的氨茶碱片，已出现呕吐、头晕等不良症状片，已出现呕吐、头晕等不良症状.按照药品使用说明书，氨茶碱的每次用量成人是按照药品使用说明书，氨茶碱的每次用量成人是100200mg，
12、儿童，儿童是是35 mg/kg.过量服用可使血药浓度过量服用可使血药浓度(单位血液容积中的药量单位血液容积中的药量)过高，过高，100g/ml浓度会浓度会出现出现严重中毒严重中毒,200g/ml浓度可致命浓度可致命.医生需要判断：孩子的血药浓度会不会达到医生需要判断：孩子的血药浓度会不会达到100200 g/ml；如果会达到，；如果会达到，应采取怎样的应采取怎样的紧急施救紧急施救方案方案.第十五张，PPT共四十二页，创作于2022年6月调查与分析调查与分析转移率正转移率正比于比于x排除率正排除率正比于比于y胃肠道胃肠道血液系统血液系统口服药物口服药物体外体外认为血液系统内药物的分布，即血药浓度
13、是均匀的，可以将血液认为血液系统内药物的分布，即血药浓度是均匀的，可以将血液系统看作一个房室，建立系统看作一个房室，建立“一室模型一室模型”.药量药量x(t)药量药量y(t)血液系统对药物的吸收率血液系统对药物的吸收率(胃肠道到血液系统的转移率胃肠道到血液系统的转移率)和排除率可以和排除率可以由由半衰期半衰期确定确定.半衰期半衰期可以从药品说明书上查到可以从药品说明书上查到.第十六张，PPT共四十二页，创作于2022年6月通常，血液总量约为人体体重的通常，血液总量约为人体体重的7%8%，体重，体重5060 kg的成年人有的成年人有4000ml左右的血液左右的血液.目测这个孩子的体重约为成年人的
14、一半，可认为其目测这个孩子的体重约为成年人的一半，可认为其血液总量约为血液总量约为2000ml.调查与分析调查与分析血药浓度血药浓度=药量药量/血液总量血液总量 口服活性炭来吸附药物，可使药物的排除率增加到原来口服活性炭来吸附药物，可使药物的排除率增加到原来（人体自身）的（人体自身）的2倍倍.临床施救的办法：临床施救的办法：体外血液透析，药物排除率可增加到原来的体外血液透析，药物排除率可增加到原来的6倍，但是安全倍，但是安全性不能得到充分保证性不能得到充分保证.第十七张，PPT共四十二页，创作于2022年6月模型假设模型假设 1.胃肠道中药物向血液的转移率与胃肠道中药物向血液的转移率与x(t)
15、成正比，比例系数成正比，比例系数(0)，总剂量，总剂量1100 mg药物在药物在t=0瞬间进入胃肠道瞬间进入胃肠道.2.血液系统中药物的排除率与血液系统中药物的排除率与y(t)成正比，比例系数成正比，比例系数(0)，t=0时血液中无药物时血液中无药物.3.氨茶碱被吸收的半衰期为氨茶碱被吸收的半衰期为5 h，排除的半衰期为，排除的半衰期为6 h.4.孩子的血液总量为孩子的血液总量为2000 ml.胃肠道中药量胃肠道中药量x(t),血液系统中药量血液系统中药量y(t)，时间，时间t以孩子以孩子误服药的时刻为起点（误服药的时刻为起点（t=0）.第十八张，PPT共四十二页，创作于2022年6月模型建立
16、模型建立x(t)下降速度与下降速度与x(t)成正比成正比(比例系数比例系数),总剂量总剂量1100mg药物在药物在t=0瞬间进瞬间进入胃肠道入胃肠道.转移率转移率正比于正比于x排除率排除率正比于正比于y胃肠道胃肠道血液系统血液系统口服药物口服药物体外体外药量药量x(t)药量药量y(t)y(t)由吸收而增长的速度是由吸收而增长的速度是x，由排除而减少的速度与，由排除而减少的速度与y(t)成正比成正比(比比例系数例系数),t=0时血液中无药物时血液中无药物.第十九张，PPT共四十二页，创作于2022年6月模型模型求解求解 药物吸收的半衰期为药物吸收的半衰期为5 h 药物排除的半衰期为药物排除的半衰
17、期为6 h 只考虑血液对药物的排除只考虑血液对药物的排除第二十张，PPT共四十二页，创作于2022年6月血液总量血液总量2000ml血药浓度血药浓度200g/ml结果及分析结果及分析 胃肠道药量胃肠道药量血液系统药量血液系统药量血药浓度血药浓度100g/mly(t)=200mg严重中毒严重中毒y(t)=400mg致命致命t=1.62t=4.87t=7.89y=442孩子到达医院前已严重中毒，如不及时施救，约孩子到达医院前已严重中毒，如不及时施救，约3h3h后将致命！后将致命！y(2)=236.5 第二十一张，PPT共四十二页，创作于2022年6月施救方案施救方案 口服活性炭使药物排除率口服活性
18、炭使药物排除率增至原来的增至原来的2倍倍.孩子到达医院孩子到达医院(t=2)就开始施救，血液中药量记作就开始施救，血液中药量记作z(t)=0.1386(不变)，=0.11552=0.2310 第二十二张，PPT共四十二页，创作于2022年6月施救方案施救方案 t=5.26z=318 施救后血液中药量施救后血液中药量z(t)显著低于显著低于y(t).z(t)最大值低于致最大值低于致命水平命水平.要使要使z(t)在施救后立在施救后立即下降，可算出即下降，可算出至至少应为少应为0.4885.若采用体外血液透析，若采用体外血液透析，可增至可增至0.11556=0.693，血液，血液中药量下降更快；临床
19、上是否需要采取这种办法，当由医生中药量下降更快；临床上是否需要采取这种办法，当由医生综合考虑并征求病人家属意见后确定综合考虑并征求病人家属意见后确定.第二十三张，PPT共四十二页，创作于2022年6月偏微分方程与数学模型2022/9/19济南大学 数学科学学院24第二十四张，PPT共四十二页，创作于2022年6月偏微分方程偏微分方程偏微分方程偏微分方程（Partial Differential Equations）指在物理学、力学、工程技术以及其他自然科学、技术指在物理学、力学、工程技术以及其他自然科学、技术科学、管理科学、甚至社会科学等的研究中归纳出来的科学、管理科学、甚至社会科学等的研究中
20、归纳出来的一些含有未知函数及其偏导数的方程一些含有未知函数及其偏导数的方程2022/9/19济南大学 数学科学学院25第二十五张，PPT共四十二页，创作于2022年6月什么是偏微分方程？什么是偏微分方程？2022/9/19济南大学 数学科学学院26物理量物理量(如位移、温度等如位移、温度等)-)-时间、空间位置时间、空间位置 -物理量的变化规律物理量的变化规律(偏微分方程偏微分方程)第二十六张，PPT共四十二页，创作于2022年6月 例子 2022/9/19济南大学 数学科学学院27第二十七张，PPT共四十二页，创作于2022年6月研究内容研究内容 2022/9/19济南大学 数学科学学院28
21、一般规律一般规律+定解条件定解条件(初始条件、边界条件初始条件、边界条件)定解问题定解问题 定解问题的适定性：定解问题的适定性：存在性（存在性（ExistenceExistence）唯一性（唯一性（UniquenessUniqueness）稳定性（稳定性（StabilityStability）+附加条件附加条件方程方程第二十八张，PPT共四十二页，创作于2022年6月4 人口预测和控制人口预测和控制 年龄分布对于人口预测的重要性年龄分布对于人口预测的重要性.只考虑自然出生与死亡，不计迁移只考虑自然出生与死亡，不计迁移.人口人口发展发展方程方程第二十九张，PPT共四十二页，创作于2022年6月人
22、口发展方程人口发展方程一阶偏微分方程一阶偏微分方程第三十张，PPT共四十二页，创作于2022年6月人口发展方程人口发展方程p0(r)已知函数已知函数(人口调查人口调查)f(t)生育率生育率(控制人口手段控制人口手段)0tr第三十一张，PPT共四十二页，创作于2022年6月生育率的分解生育率的分解 总和生育率总和生育率h生育模式生育模式0第三十二张，PPT共四十二页，创作于2022年6月人口发展方程和生育率人口发展方程和生育率总和生育率总和生育率控制生育的多少控制生育的多少生育模式生育模式控制生育的早晚和疏密控制生育的早晚和疏密 正反馈系统正反馈系统 滞后作用很大滞后作用很大第三十三张，PPT共
23、四十二页，创作于2022年6月人口指数人口指数1）人口总数）人口总数2）平均年龄）平均年龄3）平均寿命）平均寿命t时刻出生的人，死亡率按时刻出生的人，死亡率按 (r,t)计算的平均存活时间计算的平均存活时间4）老龄化指数）老龄化指数控制生育率控制生育率控制控制 N(t)不过大不过大控制控制 (t)不过高不过高第三十四张，PPT共四十二页，创作于2022年6月5 烟雾的扩散与消失烟雾的扩散与消失现象现象和和问题问题 炮弹在空中爆炸，烟雾向四周扩散，形成炮弹在空中爆炸，烟雾向四周扩散，形成圆形不透光区域圆形不透光区域.不透光区域不断扩大，然后区域边界逐渐不透光区域不断扩大，然后区域边界逐渐明亮，区
24、域缩小，最后烟雾消失明亮，区域缩小，最后烟雾消失.建立模型描述烟雾扩散和消失过程，分析消失建立模型描述烟雾扩散和消失过程，分析消失时间与各因素的关系时间与各因素的关系.问题问题分析分析 无穷空间由瞬时点源导致的扩散过程，用二阶无穷空间由瞬时点源导致的扩散过程，用二阶偏微分方程描述烟雾浓度的变化偏微分方程描述烟雾浓度的变化.观察到的烟雾消失与烟雾对光线的吸收、以观察到的烟雾消失与烟雾对光线的吸收、以及仪器对明暗的灵敏程度有关及仪器对明暗的灵敏程度有关.第三十五张，PPT共四十二页，创作于2022年6月模型模型假设假设1）烟雾在无穷空间扩散，不受大地和风的）烟雾在无穷空间扩散，不受大地和风的影响；
25、扩散服从热传导定律影响；扩散服从热传导定律.2）光线穿过烟雾时光强的相对减少与烟雾浓度）光线穿过烟雾时光强的相对减少与烟雾浓度成正比；无烟雾的大气不影响光强成正比；无烟雾的大气不影响光强.3）穿过烟雾进入仪器的光线只有明暗之分，明暗）穿过烟雾进入仪器的光线只有明暗之分，明暗界限由仪器灵敏度决定界限由仪器灵敏度决定.模型模型建立建立1）烟雾浓度）烟雾浓度 的变化规律的变化规律热传导定律：单位时间通过单位法向面热传导定律：单位时间通过单位法向面积的流量与浓度梯度成正比积的流量与浓度梯度成正比.第三十六张，PPT共四十二页，创作于2022年6月曲面积分曲面积分奥氏公式奥氏公式1）烟雾浓度）烟雾浓度
26、的变化规律的变化规律的微分形式，并利用积分中值定理的微分形式，并利用积分中值定理第三十七张，PPT共四十二页，创作于2022年6月 初始条件初始条件Q炮弹释放的烟雾总量炮弹释放的烟雾总量 单位强度的点源函数单位强度的点源函数 对任意对任意t,C的等值面是球面的等值面是球面 x2+y2+z2=R2;RC 仅当仅当 t,对任意点对任意点(x,y,z),C01）烟雾浓度）烟雾浓度 的变化规律的变化规律第三十八张，PPT共四十二页，创作于2022年6月2）光强穿过烟雾时的变化规律）光强穿过烟雾时的变化规律假设假设2）光强的相对减少与烟雾浓度成正比）光强的相对减少与烟雾浓度成正比.I(l)沿沿l方向的光
27、强，方向的光强，C(l)沿沿l方向的烟雾强度方向的烟雾强度记未进入烟雾记未进入烟雾(l l0)时光强为时光强为 I(l0)=I0第三十九张，PPT共四十二页，创作于2022年6月3）仪器灵敏度与烟雾明暗界限）仪器灵敏度与烟雾明暗界限烟雾浓度连续变化烟雾浓度连续变化烟雾中光强连续变化烟雾中光强连续变化仪器仪器z-设光源在设光源在z=-,仪器在仪器在z=,则观测到的则观测到的明暗界限为明暗界限为不透光区域有扩大、不透光区域有扩大、缩小、消失的过程缩小、消失的过程穿过烟雾进入仪器的光线只有明暗之穿过烟雾进入仪器的光线只有明暗之分，明暗界限由仪器灵敏度决定分，明暗界限由仪器灵敏度决定.不透光区域边界不透光区域边界第四十张，PPT共四十二页，创作于2022年6月4）不透光区域边界的变化规律）不透光区域边界的变化规律对任意对任意t,不透光区域边界是圆周不透光区域边界是圆周不透光区域不透光区域边界半径边界半径第四十一张，PPT共四十二页，创作于2022年6月感感谢谢大大家家观观看看第四十二张，PPT共四十二页，创作于2022年6月