计算方法复习新课题.pptVIP

例求Euler法的绝对稳定区间 所以，局部截断误差 局部截断误差为O(h3 ) 所以精度为2阶的 局部截断误差主项为： 梯形法的绝对稳定区间？ 解之得： 即：梯形法是无条件稳定的（绝对稳定） 方法 公式 显式 局部截断误差 步数 Euler yi+1=yi+hf(xi,yi) yes O(h2) 单步 后退Euler yi+1=yi+hf(xi+1,yi+1) no O(h2) 单步 梯形 no O(h3) 单步 预报-校正 yes O(h3) 单步 四阶龙格库塔 （6-20) yes O(h5) 单步 Adams外插 yes O(h5) 多步（4) Adams内插 no O(h5) 多步（3) 几种方法的比较 定理 n+1个节点的求积公式 为插值型求积公式的充要条件是公式 至少具有n次代数精度。 的代数精度 可以验证, 对于f(x)=1, x时公式两端相等, 再将f(x)=x2代入公式 左端 例2 考察求积公式 两端不相等, 所以该求积公式具有 1 次代数精度. 三个节点不具有2次代数精度， 不是插值型的 右端 例3 求证 不是插值型的 证明: 设 x0 = -1, x1 =0, x2 =1, A0 =1/2, A1=1, A2=1/2 则以这三点为插值节点的Lagrange插值 基函数为 Newton-Cotes求积公式 梯形公式，精度、截断误差 抛物线（Simpson）公式，精度、截断误差 N-C求积公式的稳定性 梯形公式 梯形公式具有1阶精度 截断误差： 抛物线（Simpson）公式 抛物线公式具有3阶精度 截断误差： Newton-Cotes求积公式的稳定性 当n≥8时,误差有可能传播扩大，收敛性也没有保证。因此，牛顿 ― 柯特斯求积公式不宜采用。 复化求积公式 复化梯形公式、复化Simpson公式、截断误差 如何根据截断误差及精度确定n值，从而写出求积公式 Romberg求积公式、意义（思想） 复化梯形公式 截断误差： 复化 截断误差： 根据截断误差及精度确定n值 由余项确定截断误差绝对值的上界表示 由给定精度确定误差限 利用组成的不等式求得n值 例9 若用复合梯形公式计算积分 问积分区间要等分多少才能保证有误差不超过 写出相应的复化梯形公式(重点) 解 由余项 则当0＜x＜1时,有 又 故 即: 即: 所以: 则: 取 而用复合simpson公式n=3.706 取n=4即能达到这个精度 相应的复化梯形公式为: Romberg求积公式 Gauss型求积公式 Gauss型求积公式的定义，判断 具有2n+1次代数精度的插值型求积公式 数值微分 两点公式及其截断误差 三点公式及其截断误差 两点公式 截断误差 三点公式 三点公式 截断误差 第六章 常微分方程初值问题的数值解法 欧拉法 欧拉法、改进欧拉法、预报-校正公式 截断误差 精度 显式（隐式）、单步（多步） 绝对稳定区间的求法 龙格-库塔方法 思想、精度 显式（隐式）、单步（多步） 欧拉法 后退欧拉法 梯形公式 预报-校正公式 预报-校正公式 龙格-库塔公式 补充定义2：用单步法解模型方程 得到解满足稳定性方程 若 方法是绝对稳定的 例证明解常微分方程初值问题的梯形方法精度是二阶的，并求其局部截断误差主项 而局部截断误差按定义是： y(x i+1)-yi+1= y(x i+h)- yi+1 先将y(x i+h)在x=xi处tailor展开 分析：证明精度是二阶，即证明其局部截断误差是o(h3） 证明: 再来看 yi+1 代入上式 yi+1 (i=2,….n, j=1,....i-1) (p54 3-25) -------按行计算公式 按列计算公式： 对j=1,2,…n i=j+1,j+2,…..n 迭代方法 Jacobi迭代法 A=D-L-U Gauss-Seidel迭代法 公式及收敛条件（充分条件、充要条件） Jacobi迭代法 分量形式 矩阵形式 记 ， Jacobi迭代法 Jacobi迭代收敛条件 充分条件(1) :迭代矩阵至少存在一种矩阵范数||.||,