chap7常微分的方程的数值解法.pptVIP

第7章 常微分方程（组）的数值解法;第7章 常微分方程（组）的数值解法;第7章 常微分方程（组）的数值解法;定义7.0.1 若存在常数 L 0, 使得对一切的 x∈[a , b] 及 y , , 均有 ; 我们首先给出常微分方程初值问题解的存在惟一性定理。;然后在节点上建立逼近于原初值问题的计算格式 (或差分格式), 由此计算出原问题的解 y( x ) 在节点 x1 , x2 , . . . , xN 处的近似值： y1 , y2 , . . ., yN， 称它们为常微分方程初值问题的数值解. 相邻两个节点的距离 hn= xn+1 - xn 称为步长, 通常取定步长h 0, 即节点 xn = x0 + nh , n = 0,1, … , N.;7.1 初值问题计算格式的建立;略去余项, 并以数值解 yn , yn+1 代替 y (xn) 及 y (xn+1), 则得差分方程;如果利用下列数值微分公式;类似地，可导出;设 y (x)∈C2[a , b], 由 Taylor 公式有;3. 数值积分方法;若对此式的积分采用梯形公式, ;上式称为梯形公式. 由于它关于 yn+1 成隐式形式, 故其为隐格式. 隐格式求解比较困难, 当 yn 已知时, 要求 yn+1 ，需解关于 yn+1 的非线性方程. 在实际应用时, 上式常与 Euler 公式联合使用, 构成如下计算格式:;隐式梯形公式的迭代格式;定理 7.1.1 设 f (x , y) 在区域 D 上关于 y 满足 Lipschitz 条件, 即; 为了减少计算量, 可采用预测-校正格式. 方法是先用Euler公式求得一个初始近似值 称为预测值, 再把 带入梯形公式右端计算一次求得 yn+1 称之为校正值, 即;例7.1.1：利用Euler公式与改进的Euler公式求解初值问题（步长h = 0.1）;计算结果见下表(见书P227表7.1) ;解此问题的改进的Euler公式为 ;7.1.2 截断误差与方法的精度;整体截断误差 en+1 是在没有引进舍入误差的情况下, 纯粹因为不准确的计算格式造成的, 故又称为方法误差.它不仅与 x = xn+1 这一步的计算有关, 而且和 xn , xn-1 ,. . . , x1 这几步的计算都有关系. ;如Euler公式;定义7.1.2 若某数值方法的局部截断误差为 则称该方法具有 P 阶精度, 或称其为 P 阶方法.;7.2 Runge-Kutta 方法;进一步假设设 y (x)∈Cp+1[a , b],由 Taylor 公式有;7.2 Runge-Kutta 方法;在(7.2.3)中略去余项,用yn , yn+1 代替 y (xn) 及 y (xn+1);由于局部截断误差 ;7.2.1 Runge-Kutta方法的基本思想;其一般形式为:; 7.2.2 二阶 Runge-Kutta 方??;;得到;再根据 y(xn+1) 在点 xn的一元 3 阶 Taylor 展开式;即由; 取; 7.2.3 四阶Runge-Kutta方法; Gill 公式;例7.2.1 用标准的 四阶Rung-Kutta 法解初值问题,取步长h=0.2.;7.5一阶常微分方程组与高阶方程初值问题的数值解法 ;写成向量形式： ;利用向量值函数的微积分理论, 很容易推导出一阶常微分方程组初值问题的数值解法. ;(7.5.3)的分量形式为 ;四阶标准的Runge-Kutta公式 ;四阶标准的Runge-Kutta公式的分量形式 ;例7.5.1：试写出用中点公式解下列初值问题的计算公式. ;中点公式的向量形式 ;整理得分量计算格式;7.5.2 高阶常微分方程初值问题的数值解法;等价的一阶常微分方程组初值问题：;例7.5.2写出用标准四阶Runge-Kutta公式求解;求此问题数值解的标准四阶Runge-Kutta公式为;7.5.3 高阶常微分方程组初值问题的数值解法;等价的一阶常微分方程组初值问题;本章小结