常微分方程数值解法.ppt

第九章 常微分方程的数值解法 ;主 要 内 容;现实世界中大多数事物;如何利用数值方法求解微分方程（组）的问题。;寻找解析解的过程称为求解微分方程组。 ;4.什么是微分方程的数值解?;在大量的实际方程中出现的函数起码的连续性都 无法保证，更何况要求阶的导数;常微分方程的数值解法常用来求近似解;本章主要讨论一阶常微方程的初值问题;则称 f (x,y) 对y 满足李普希兹条件，L 称为Lipschitz常数.;我们以下的讨论，都在满足上述条件下进行.;写成向量形式为;也就解决了高阶方程的定解问题.;简单的数值方法与基本概念 ;1. 简单的欧拉(Euler)方法;2. 欧拉方法的导出;对微分方程(1.1)两端从;右端积分用左矩形数值求积公式;x0;3.欧拉公式有明显的几何意义;也称欧拉折线法.;从上述几何意义上得知，由Euler法所得的 折线明显偏离了积分曲线，可见此方法 非常粗糙。;定义;如果单步差分公式的局部截断误差为O(hp+1), 则称该公式为p阶方法.这里p为非负整数. 显然,阶数越高,方法的精度越高.;5. 欧拉公式的改进：;6.梯形公式 /* trapezoid formula */;方 法;? 改进欧拉法 /* modified Euler’s method */;§3 龙格 - 库塔法 /* Runge-Kutta Method */;首先希望能确定系数 ?1、?2、p，使得到的算法格式有2阶精度，即在 的前提假设下，使得 ;Step 3: 将 yi+1 与 y( xi+1 ) 在 xi 点的泰勒展开作比较;其中?i ( i = 1, …, m )，?i ( i = 2, …, m ) 和 ?ij ( i = 2, …, m; j = 1, …, i?1 ) 均为待定系数，确定这些系数的步骤与前面相似。 ;注： ? 龙格-库塔法的主要运算在于计算 Ki 的值，即计算 f 的值。Butcher 于1965年给出了计算量与可达到的最高精度阶数的关系：;§4 收敛性与稳定性 /* Convergency and Stability */;? 2.稳定性 /* Stability */;定义;例：考察显式欧拉法;例：隐式龙格-库塔法;例1 用欧拉方法，隐式欧拉方法和欧拉中点公式计算;本题精确解为y=e-x,计算结果见表9-1;梯形公式只校正一次的格式为;1. Runge-Kutta 法的一般形式 2. 2阶Runge-Kutta 方法 3. 经典Runge-Kutta 方法 4．R-K-Fehhlberg 方法 5. 隐式R-K方法 6. 变步长方法 ;1.Runge-Kutta 法的一般形式;式(9.11) 称L级p阶Runge-Kutta方法(简称R-K法)。 ??L＝1就是欧拉法，当L＝2时为改进的欧拉法。 ;其局部截断误差是：;3. 经典Runge-Kutta方法;4．R-K-Fehhlberg 方法;;5. 隐式R-K方法;通常，同级的隐式公式获得比显式公式更高的阶。常用的隐式 R-K法有:;6.变步长方法;例 分别用改进的欧拉格式和四阶龙格—库塔格式解初值问题（取步长h=0.2）：;节点 改进欧拉法 四阶龙格—库塔法 准确解 0 1 1 1 0.2 1.186667 1.183229 1.183216 0.4 1.348312 1.341667 1.341641 0.6 1.493704 1.483281 1.483240 0.8 1.627861 1.612514 1.612452 1 1.754205 1.732142 1.732051;单步法的相容性、收敛性和稳定性 ;单步法的相容性;单步法的收敛性;单步法的稳定性;相容性和收敛性的关系 ;相容性和方法阶的关系 ; 稳定性和收敛性的关系;绝对稳定性和绝对稳定域 ;绝对稳定性和绝对稳定域2;则Euler 方法的绝对稳定域是以 为半径的圆的内部。为了保证稳定性步长有所控制。 假如当 时h应满足 ，当 时， h 应满 足