微分方程数值解法59954.pptVIP

高阶微分方程的求解 用MATLAB解常微分方程-实例1 SI 用MATLAB解常微分方程-实例2 SIR 课堂练习 常微分方程数值解 1. 两个最常用的数值解法: 欧拉（Euler）方法 龙格-库塔（Runge-Kutta）方法 主要内容 2. 龙格-库塔方法的MATLAB实现 区分已感染者(病人)和未感染者(健康人) SI 模型 传染病模型 Logistic 模型 传染病有免疫性——病人治愈后即移出感染系统，称移出者 SIR模型 传染病模型 无法求出 的解析解 相平面分析方法（Phase Plane Analysis) 时间t的相（像） 的轨迹称为相轨线 相平面分析是在相平面上分析相轨线随时间变化的变化情况，是一种定性分析的方法。它也是应用中非常重要的一种分析方法。 “常微分方程初值问题数值解”的提法 ? ，而在一系列离散点 不求解析解 通常取等步长h y1 y2 yn 欧 拉 方 法 基本思路 P1 P2 P3 各种欧拉公式 y P0 x0 x1 x2 x3 x y=y(x) y0 取不同的点 向前欧拉公式 x取左端点xn 近似：yn?y(xn) ,yn+1 ? y(xn+1) 误差分析 局部截断误差主项为 假设到第n步公式右端yn没有误差,即yn=y(xn),从xn到xn+1一步的计算值yn+1与精确值y(xn+1) 之差,称为局部截断误差