数学建模微分方程模型分析.pptVIP

生育率的分解 ?~总和生育率 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