9传染病模型与微分方程数值解.ppt

北京邮电大学数学系 第五章 微分方程模型 常微分方程的数值解 5.1 传染病模型 常微分方程的数值解及实验 在生产和科研中所处理的微分方程往往很复杂且大多得不出一般解。而在实际上对初值问题，一般是要求得到解在若干个点上满足规定精确度的近似值，或者得到一个满足精确度要求的便于计算的表达式。 高数中微分方程解法在实际中基本不会直接使用 (一)常微分方程数值解 因此，研究常微分方程的数值解法十分必要。 1、用差商代替导数 若步长h较小，则有 故有公式： 此即欧拉法(向前欧拉法). (二)建立数值解法的一些途径 对应有隐式欧拉法 2、使用数值积分 对方程y’=f(x,y), 两边由xi到xi+1积分,并利用梯形公式, 有 实际应用时，与欧拉公式结合使用 故有公式 梯形方法/*trapezoid formula*/ 此即改进的欧拉法 ? 中点欧拉公式 /* midpoint formula */ 中心差商近似导数 x0 x2 x1 3、使用泰勒公式 以此方法为基础，有龙格-库塔法、线性多步法等方法 库塔三阶方法 四阶龙格-库塔公式 1、在解n个未知函数的方程组时，x0和x均为n维向量，m-文件中的待解方程组应以x的分量形式写成. 2、使用Matlab软件求数值解时，高阶微分方程必须等价地变换成一阶微分方程组. 注意: 设取步长 ，从 到 用四阶龙格-库塔方法 微分方程求解实例 求解初值问题 h=0.2; ts=0:h:1; y0=1; [t, x]=ode45(dfun1,ts,y0); [t, x], plot(t,x) function dx=dfun1(x,y) dx=y-2*x/y; 0 1.0000 0.2000 1.1832 0.4000 1.3416 0.6000 1.4832 0.8000 1.6125 1.0000 1.7321 建立m-文件dfun1.m如下 输入命令 3、结果如图 解: 令 y1=x，y2=y1’, 1、建立m-文件dfun2.m如下： function dx=dfun2(t,y) dx=[y(2); (1-y(1)^2)*y(2)-y(1)]; 例 则微分方程变为一阶微分方程组： [t,y]=ode45(@dfun2,[0,20],[2,0]); [t,y] plot(t,y(:,1),r-,t,y(:,2),b.-); hold on plot(y(:,1),y(:,2),co); hold off legend(t ~ x, t ~ x`, x ~ y); 2、取t0=0, tf=20, 输入命令： 3、结果如图 [t,x]=ode45(@dfun2,[0,20],[2,0]) plot(t,x(:,1),r-,t,x(:,2),b.-); hold on plot(x(:,1),x(:,2),co); hold off legend(t ~ x,t ~ x`,x ~ y); 描述对象特征随时间(空间)的演变过程 分析对象特征的变化规律 预报对象特征的未来性态 研究控制对象特征的手段 根据函数及其变化率之间的关系确定函数 微分方程建模 根据建模目的和问题分析作出简化假设 按照内在规律或用类比法建立微分方程 动态模型 5.1 传染病模型 描述传染病的传播过程 分析受感染人数的变化规律 预报传染病高潮到来的时刻 预防传染病蔓延的手段 本世纪初, 瘟疫常在世界上某地流行, 随着人类文明的不断进步, 很多疾病, 诸如天花、霍乱已经得到有效的控制. 然而, 即使在今天, 一些贫穷的发展中国家，仍出现传染病流行的现象,医疗卫生部门的官员与专家所关注的问题是： 问题提出 感染疾病的人数与哪些因素有关? 问题分析 不同类型传染病的传播过程有不同的特点, 故不从医学的角度对各种传染病的传播过程一一进行分析，而是按一般的传播机理建立模型． 由于传染病在传播的过程涉及因素较多，在分析问题的过程中，不可能通过一次假设建立完善的数学模型. 思路： 针对结果中的不合理之处，逐步修改假设，最终得出较好的模型。 先做出最简单的假设，对得出的结果进行分析， 模型一 模型假设： （1）一人得病后，久治不愈，人在传染期内不会死亡。 (2) 假设每个病人每天有效接触(足以使人致病)人数为? 设已感染人数 (病人) x(t), 假设是连续可微函数 建模 ? 举个实例 x=0:0.1:10; y=exp(x); plot(x,y,b-); 最初只有1个病人，1