8第八章数值积分.pptVIP

第八章 数值积分 近似计算 8.1 插值型求积公式 节点 f (x) 插值型积分公式 误差 8.2 复化求积公式 如果积分区间比较大，直接地使用上述求积公式， 精度难以保证。 高次插值有Runge 现象，故采用分段低次插值 ? 分段低次合成的 Newton-Cotes 复合求积公式。 （1）等分求积区间，比如取步长 ，分[a, b]为n等分， 分点为 k = 0, 1, 2,…, n （2）在区间 [xk, xk+1]上使用以上求积公式求得Ik （3）取和值 ，作为整个区间上的积分近似值。 ? 复化梯形公式： = Tn /*积分中值定理*/ ? 复化 Simpson 公式： = Sn 例 8．1：利用数据表 计算积分 这个问题有明显的答案 取n = 8用复化梯形公式 取n=4,用辛卜生公式 8.3 变步长梯形方法 8.4 求积公式的误差 舍入误差： 取f (x) ? 1，则 若f (xk)的舍入误差小于? ，则 1．梯形公式的截断误差 2．辛卜生公式的截断误差 8.5 龙贝格求积公式 龙贝格积分法是在计算梯形和序列的基础上应用了 线性外推的加速方法，由此构成的一种具有超线性 收敛的自动积分法 方法思路 : 1.按照区间逐次分半的方法，计算梯形和序列 由此生成序列 T0, T1, …, Tn ,… o h2 T Sn I T Tn+1 Tn h2 设Tn为梯形和，I为积分真值，由复化梯形公式 f(x) 2.加速 由解析几何 令h = 0， 则此直线在T 轴上的截距为 用类似方法可推得： 柯特斯序列 龙贝格序列 由此法，可得如下三角形数表 计算方法的实现： 首先构造T数表： 计算步骤： 2．对k = 1, 2,… 计算下列各步 3．对n = 0, 1, 2,…, k = n – 1, n – 2, … 4．收敛控制 8.6 高斯型求积公式 问题：在节点个数一定的情况下，是否可以在[a, b]上自由 选择节点的位置，使求积公式的精度提得更高 ？ 代数精确度：称： 为一般求积公式。这里Ak为不依赖f (x)的常数若（8.9） 对任意不高于m次的多项式精确成立，而对于xm+1不能 精确成立，就说（8.9）式具有m次代数精确度。 例 8.2：求形如 的两点求积公式。 （1）用梯形公式（即以x0 = -1，x1 = 1为节点的插值型 求积公式）立即可得 一次代数精确度。 （2）若对求积公式中的四个待定系数A0, A1, x0, x1适当选取， 使求积公式对f (x) = 1，x，x2，x3都准确成立 o x y a b A B f(x)