数学建模微分方程模型解读.pptVIP

  1. 1、本文档共42页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  5. 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  6. 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  7. 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  8. 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
数学建模微分方程模型解读

生育率的分解 ?~总和生育率 h~生育模式 0 人口发展方程和生育率 ~总和生育率——控制生育的多少 ~生育模式——控制生育的早晚和疏密 正反馈系统 滞后作用很大 人口指数 1)人口总数 2)平均年龄 3)平均寿命 t时刻出生的人,死亡率按 ?(r,t) 计算的平均存活时间 4)老龄化指数 控制生育率 控制 N(t)不过大 控制 ?(t)不过高 5 烟雾的扩散与消失 现象和 问题 炮弹在空中爆炸,烟雾向四周扩散,形成圆形不透光区域. 不透光区域不断扩大,然后区域边界逐渐明亮,区域缩小,最后烟雾消失. 建立模型描述烟雾扩散和消失过程,分析消失时间与各因素的关系. 问题分析 无穷空间由瞬时点源导致的扩散过程,用二阶偏微分方程描述烟雾浓度的变化. 观察到的烟雾消失与烟雾对光线的吸收、以及仪器对明暗的灵敏程度有关. 模型假设 1)烟雾在无穷空间扩散,不受大地和风的影响;扩散服从热传导定律. 2)光线穿过烟雾时光强的相对减少与烟雾浓度成正比;无烟雾的大气不影响光强. 3)穿过烟雾进入仪器的光线只有明暗之分,明暗界限由仪器灵敏度决定. 模型建立 1)烟雾浓度 的变化规律 热传导定律:单位时间通过单位法向面积的流量与浓度梯度成正比. 曲面积分奥氏公式 1)烟雾浓度 的变化规律 的微分形式,并利用积分中值定理 初始条件 Q~炮弹释放的烟雾总量 ? ~单位强度的点源函数 对任意t, C的等值面是球面 x2+y2+z2=R2; R??C? 仅当 t??, 对任意点(x,y,z), C?0 1)烟雾浓度 的变化规律 2)光强穿过烟雾时的变化规律 假设2)光强的相对减少与烟雾浓度成正比. I(l) ~ 沿l方向的光强, C(l) ~ 沿l方向的烟雾强度 记未进入烟雾(l?l0)时光强为 I(l0)=I0 3)仪器灵敏度与烟雾明暗界限 烟雾浓度连续变化 烟雾中光强连续变化 仪器 z -? ? 设光源在z=-?, 仪器在z=?,则观测到的明暗界限为 不透光区域有扩大、缩小、消失的过程 穿过烟雾进入仪器的光线只有明暗之分,明暗界限由仪器灵敏度决定. ~不透光区域边界 4)不透光区域边界的变化规律 对任意t, 不透光区域边界是圆周 不透光区域边界半径 r(t) rm 0 t1 t2 t 结果分析 观测到不透光区域边界达到最大的时刻t1,可以预报烟雾消失的时刻t2 微分方程模型 2 如何预报人口的增长 3 如何施救药物中毒 4 人口预测和控制模型 1 目标跟踪问题 动态模型 描述对象特征随时间(空间)的演变过程. 分析对象特征的变化规律. 预报对象特征的未来性态. 研究控制对象特征的手段. 根据函数及其变化率之间的关系确定函数. 微分方程建模 根据建模目的和问题分析作出简化假设. 按照内在规律或用类比法建立微分方程. 1 目标跟踪问题 设位于坐标原点的甲舰向位于x轴上点A(1, 0)处的乙舰发射导弹,导弹头始终对准乙舰.如果乙舰以最大的速度v0(常数)沿平行于y轴的直线行驶,导弹的速度是5v0,求导弹运行的曲线方程.乙舰行驶多远时,导弹将它击中? 由(1),(2)消去t, 整理得模型: 解法二(数值解法) 1.建立M文件eq1.m 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=0.9999,建立主程序如下: x0=0,xf=0.9999 [x,y]=ode15s(eq1,[x0 xf],[0 0]); plot(x,y(:,1),’b.) hold on y=0:0.01:2; plot(1,y,’b*) 结论: 导弹大致在(1,0.2)处击中乙舰. 令y1=y, y2=y1’,将方程(3)化为一阶微分方程组. 结果见图 导弹大致在(1,0.2)处击中乙舰,与前面的结论一致. 返 回 结论:时刻t=0.21时,导弹在(1,0.21)处击中乙舰. 背景 年份 1625 1830 1930 1960 1974 1987 1999 人口(亿) 5 10 20

文档评论(0)

22ffbqq + 关注
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档