数值分析第4章.pptVIP

数值积分 　　例１　验证求积公式　　　　　　　　　　　　　　　　　　　　　　　　　　 　 　具有3次代数精确度． 　　解：　当　　　　　　　　　　　　　 　　而　　　　　　　　　　　　 　　　　　　　　　　　　　　　　　　 　　有　　　 （1）当　　　　　　　　　　　　　　　 　　　　　　　　　　　　　　　　　　 （2）当 （3）当　　　　　 （1）当　　　　　　　　　　　　　　　 　　　　　　　　　　　　　　　　　　 故求积公式具有三次代数精确度． 插值型求积公式 这一节所讨论的求积公式，都是用在区间[a, b]上对被积函数f(x)作插值所得插值多项式Pn(x)代替被积函数f(x)导出的公式．这一类求积公式的求积节点xk，就是对f(x)作插值时的插值节点，所以这类求积公式称为插值型求积公式． 　　　 科特斯系数表 几个常见公式 举例（一） 复合求积公式 复合梯形公式 复合simpson公式 复合科特斯公式 举例（二） h 很小时的误差 例：分别用梯形公式和simpson公式计算积分 解： a＝0, b＝1, f (x) = e -x ， 由 simpson 公式可得 由梯形公式可得 与精确值 0.6321 相比得误差分别为 0.0518 和 0.0002。 提高积分计算精度的常用两种方法 用 复合公式 用 非等距节点 复合求积公式：将积分区间分割成多个小区间，然后在每个小区间上使用低次牛顿－科特斯求积公式。 将[a, b] 分成 n 等分 [xi , xi+1] ，其中节点 (i = 0, 1, …, n) 数值积分公式与多项式插值有很大的关系。因此Runge现象的存在，使得我们不能用太多的积分点计算。采用与插值时候类似，我们采用分段、低阶的方法 复合梯形公式: Tn 余项： ，? ?(a, b) 复合simpson公式: Sn 余项： ，? ?(a, b) 复合cotes公式: Cn 余项： ，? ?(a, b) 解： 例：设 ，利用下表中的数据分别用复合梯形公式和复合simpson公式计算积分 0.936 5/8 0.909 6/8 0.990 2/8 0.977 3/8 0.841 0.877 0.954 0.997 1 f (xi ) 1.0 7/8 1/8 0 4/8 xi ?i ?(xi, xi+1 ) (h 0) 定积分定义 即 同理 4.4 龙贝格求积公式 一、梯形法的递推化 二、龙贝格公式 用梯形公式得到的积分 近似值 的误差大致是 ,因此人们期望, 如果用这个误差作为对 的一种补偿,则得到的求积公式 (1) 的代数精度会有所提高。 通过直接验证可知 (2) 也就是说 , 用梯形公式二分前后的两个积分值 与 按照公式 (1) 做线形组合,其结果正好是用抛物线公式得到的积分值 。 用抛物线公式得到的积分近似值 的误差大致是 ,因此对抛物线公式进行修正, 得到 (3) 也就是说 , 用抛物线公式二分前后的积分值 与 按照公式(3) 作线形组合, 其结果正好是用柯特斯公式得到的积分值 。 通过直接验证可知 (4) 用柯特斯公式得到的积分近似值 的误差大致是 ,因此, 对柯特斯公式进行修改,得到求积公式 (5) 为此,构造求积公式 (6) 称(6)式为龙贝格(Romberg)公式。 龙贝格公式是一种计算积分的方法。在变步长的求积过程中, 运用(2) , (4) , (6) 式可以将精度低的梯形值 逐步加工成精度较高的抛物线 ,柯特斯值 与龙贝格值 。 计算f(a) , f(b) , 算出 。 (2) 把[ a , b] 2等分, 计算 , 算出 与 。 (3) 把[ a , b] 4等分,计算 算出 与 。 龙贝格求积的计算步骤如下： (5) 把[ a , b] 16等分, 计算 (4) 把[ a , b] 8 等分, 计算 算出 与 与