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

北京邮电大学数学系 第五章 微分方程模型 常微分方程的数值解 5.1 传染病模型 常微分方程的数值解及实验 在生产和科研中所处理的微分方程往往很复杂且大多得不出一般解。而在实际上对初值问题，一般是要求得到解在若干个点上满足规定精确度的近似值，或者得到一个满足精确度要求的便于计算的表达式。 高数中微分方程解法在实际中基本不会直接使用 (一)常微分方程数值解 因此，研究常微分方程的数值解法十分必要。 1、用差商代替导数 若步长h较小，则有 故有公式： 此即欧拉法(向前欧拉法). (二)建立数值解法的一些途径 对应有隐式欧拉法 2、使用数值积分 对方程y’=f(x,y), 两边由xi到xi+1积分,并利用梯形公式, 有 实际应用时，与欧拉公式结合使用 故有公式 梯形方法/*trapezoid formula*/ 此即改进的欧拉法 ? 中点欧拉公式 /* midpoint formula */ 中心差商近似导数 x0 x2 x1 4、数值公式的精度 当一个数值公式的截断误差可表示为o(hk)时(k为正整数，h为步长)，称它是一个k阶公式。 k越大，则数值公式的精度越高。 欧拉法是一阶公式，改进的欧拉法是二阶公式。 线性多步法有四阶阿达姆斯外插公式和内插公式。 [t，x]= solver( ‘f’, ts, x0, options） ode45 ode23 ode113ode15sode23s 由待解方程写成的m-文件名 ts=[t0，tf]，t0、tf为自变量的初值和终值 函数初值条件 自变量值 函数值 用于设定误差限(缺省时设定相对误差10-3, 绝对误差10-6), 命令为：options=odeset(’reltol’,rt,’abstol’,at), rt, at：分别为设定的相对误差和绝对误差. (三)用Matlab软件求常微分方程的数值解 [t,x]=ode23(‘f’, ts, x0) 3级2阶龙格-库塔公式 [t,x]=ode45(@f, ts, x0) 5级4阶龙格-库塔公式 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 传染病模型 描述传染病的传播过程 分析受感染人数的变化规律 预报传染病高潮到来的时刻 预防传染病蔓延的手段 本世纪初, 瘟疫常在世界上某地流行, 随着人类文明的不断进步, 很多疾病, 诸如天花、霍乱已经得到有效的控制.