第九讲常微分方程模型2.pptVIP

Malthus模型预测美国人口 Malthus模型预测美国人口 Malthus模型预测的优缺点 阻滞增长模型预测美国人口 阻滞增长模型预测美国人口 阻滞增长模型预测的优缺点 参数 估计法二 % 构造非线性拟合函数 [TWOEXPS.M] function E=twoexps(a,x,y) x=x(:);y=y(:); Y=a(1)*exp(-a(3)*x)+a(2)*exp(-a(4)*x); E=sum((y-Y).^2) a0=[100 1 1 ]; options=optimset(fminsearch); options.TolX=0.01; options.Display=off; a=fminsearch(@ps,a0,options,x,y) 被传染人数的估计法二 X=fzero(‘x-1.2*log(x/0.96)-0.99’,0.5) X=0.8651 % M 函数 function dy=cwf1 (t,y) dy=zeros(2,1); dy(1)=0.1*y(1)*(1-0.001*y(1)-0.0008*y(2)); dy(2)=0.2*y(2)*(1-0.0012*y(1)-0.001*y(2)); %m程序 function dy=cwf2 (t,y) dy=zeros(2,1); dy(1)=0.1*y(1)*(1-0.001*y(1)+0.0005*y(2)); dy(2)=0.2*y(2)*(-1+0.0015*y(1)-0.001*y(2)); [T,X]=ode45(‘cwf2’,[0,100],[1600,2800]); [T,Y]=ode45(‘cwf2’,[0,100],[1200,2500]); [T,Z]=ode45(‘cwf2’,[0,100],[2500,2200]); Plot(X(:,1),X(:,2),Y(:,1),Y(:,2),Z(:,1),Z(:,2)) Text(2000,2600, 2000,2600,{\sigma_11,\sigma_21,\sigma_1\sigma_21}) 龙格—库特方法 考虑微分中值定理 2阶龙格－库特公式 4阶龙格－库特公式 一阶方程组与高阶方程 Volterra方程解曲线 与Volterra方程相轨线 利用泰勒展式得到 s i 1 0 1 D 模型4 SIR模型 相轨线 及其分析 传染病蔓延 传染病不蔓延 s(t)单调减?相轨线的方向 P1 s0 im P1: s01/? ? i(t)先升后降至0 P2: s01/? ? i(t)单调降至0 1/? ~阈值 P3 P4 P2 S0 模型4 SIR模型 预防传染病蔓延的手段 ? (日接触率)? ? 卫生水平? ?(日治愈率)? ? 医疗水平? 传染病不蔓延的条件——s01/? ? 的估计 降低 s0 提高 r0 提高阈值 1/? 降低 ?(=?/?) ? ?, ? ? 群体免疫 SIR模型 被传染人数的估计法一 记被传染人数比例 xs0 i 0 P1 ? i0 ?0, s0 ?1 ? 小, s0 ? ?1 提高阈值 1/? ?降低被传染人数比例 x s0 - 1/? = ? 6.4多种群生态数学模型 意大利生物学家D’Ancona曾致力于鱼类种群相互制约关系的研究，他从第一次世界大战期间,地中海各港口捕获的几种鱼类捕获量百分比的资料中，发现鲨鱼等的比例有明显增加（见下表），而供其捕食的食用鱼的百分比却明显下降.显然战争使捕鱼量下降，食用鱼增加，鲨鱼等也随之增加，但为何鲨鱼的比例大幅增加呢？ 他无法解释这个现象，于是求助于其岳父，著名的意大利数学家V.Volterra，希望建立一个食饵—捕食系统的数学模型，定性或定量地回答这个问题. 该模型反映了在没有人工捕获的自然环境中食饵与捕食者之间的制约关系，没有考虑食饵和捕食者自身的阻滞作用，是Volterra提出的最简单的模型. 模型（一）不考虑捕获 定理 Volterra微分方程组对应初值问题 的解 是周期函数，且解的周期平均值为 首先，建立m-文件shier.m如下： function dx=shier(t,x) dx=zeros(2,1); dx(1)=x(1)*(1-0.1*x(2)); dx(2)=x(2)*(-0.5+0.02*x(1)); 其次，建立主程序shark.m如下： [t,x]=ode45(shier,[0 1