人口增长模型.pptVIP

微分方程模 型 华中农业大学数学建模实践基地 姓名: 侯志敏 联系方式: E-mail: houzhimin@yahoo.cn QQ: 电话: 微分方程模型 涉及“改变”、“变化”、“增加”、“减少”、“衰变”、“边际”、“速度”、 “运动”、“追赶”、“逃跑”、、、等等词语的确定性连续问题。 1、寻找改变量 一般说来微分方程问题都遵循这样的文字等式 变化率（微商）=单位增加量--单位减少量 等式通常是利用已有的原则或定律。 3、用微元法建立微分方程； 4、确定微分方程的定解条件（初边值条件）； 5、求解或讨论方程（数值解或定性理论）； 6、模型和结果的讨论与分析。 Malthus模型预测美国人口 Malthus模型预测美国人口 Malthus模型预测的优缺点 阻滞增长模型预测美国人口 阻滞增长模型预测美国人口 阻滞增长模型预测的优缺点 已感染人数 (病人) i(t) 每个病人每天有效接触(足以使人致病)人数为? 模型1 假设 若有效接触的是病人，则不能使病人数增加 必须区分已感染者(病人)和未感染者(健康人) 建模 ? 模型2 区分已感染者(病人)和未感染者(健康人) 假设 1）总人数N不变，病人和健康 人的 比例分别为 2）每个病人每天有效接触人数为?, 且使接触的健康人致病 建模 ? ~ 日 接触率 SI 模型 模型2 1/2 tm i i0 1 0 t tm~传染病高潮到来时刻 ? (日接触率)? ? tm? Logistic 模型 病人可以治愈！ ? t=tm, di/dt 最大 模型3 传染病无免疫性——病人治愈成为健康人，健康人可再次被感染 增加假设 SIS 模型 3）病人每天治愈的比例为? ? ~日治愈率 建模 ? ~ 日接触率 1/? ~感染期 ? ~ 一个感染期内每个病人的有效接触人数，称为接触数。 * * * b、微分方程建模的基本方法 a、微分方程建模的对象 （1）根据规律列方程 利用数学、力学、物理、化学等学科中的定理或经过实验检验的规律等来建立微分方程模型。 （2）微元法 利用已知的定理与规律寻找微元之间的关系式，与第一种方法不同的是对微元而不是直接对函数及其导数应用规律。 （3）模拟近似法 在生物、经济等学科的实际问题中，许多现象的规律性不很清楚，即使有所了解也是极其复杂的，建模时在不同的假设下去模拟实际的现象，建立能近似反映问题的微分方程，然后从数学上求解或分析所建方程及其解的性质，再去同实际情况对比，检验此模型能否刻画、模拟某些实际现象。 c、微分方程建模的基本规则 2、对问题中的特征进行数学刻画 导弹追踪问题 设位于坐标原点的甲舰向位于x轴上点A(1, 0)处的乙舰发射导弹，导弹头始终对准乙舰.如果乙舰以最大的速度v0(是常数)沿平行于y轴的直线行驶，导弹的速度是5v0，求导弹运行的曲线方程.又乙舰行驶多远时，导弹将它击中？ 解法一（解析法） 由(1),(2)消去t整理得模型: To Matlab(chase1) 轨迹图见程序chase1 解法二(数值解) 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，建立主程序ff6.m如下： 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）处击中乙舰 To Matlab(ff6) 令y1=y,y2=y1’，将方程（3）化为一阶微分方程组。 解法三(建立参数方程求数值解) 设时刻t乙舰的坐标为(X(t)，Y(t))，导弹的坐标为(x(t)，y(t)). 3．因乙舰以速度v0沿直线x=1运动，设v0=1，则w=5，X=1，Y=t 4. 解导弹运动轨迹的参数方程 建立m-文件eq2.m如下： function dy=eq2(t,y) dy=zeros(2,1); dy(1)=5*(1-y(1))/sqrt((1-y(1))^2+(t-y(2))^2); dy(2)=5*(t-y(2))/sqrt((1-y(1))^2+(t-y(2))^2); 取t0=0，tf=2，