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