[]第四章 数值积分.ppt

4.1.2 代数精度的概念 4.1.4 求积公式的余项 4.2.2 偶阶求积公式的代数精度 function I=T_quad(x,y) n=length(x);m=length(y); h=(x(n)-x(1))/(n-1); a=[1 2*ones(1,n-2) 1]; I=h/2*sum(a.*y); x=-1:0.1:1; y=exp(-x.^2); I=T_quad(x,y) 4.4.3 理查森外推加速法 4.5 自适应积分方法 函数变化有急有缓，为了照顾变化剧烈部分的误 差，我们需要加密格点。对于变化缓慢的部分，加 密格点会造成计算的浪费。以此我们介绍一种算 法，可以自动在变化剧烈的地方加密格点计算，而 变化缓慢的地方，则取稀疏的格点。 针对这类问题的算法技巧是在不同区间上预测被 积函数变化的剧烈程度确定相应步长，该法称为自 适应积分。 以复合辛普森公式为例说明该算法的思想。 n xk(n) Ak(n) Rn 1 0 2 ±0.5773503 1 3 ± 0.7745967 5/9=0.5555556 0 8/9=0.8888889 4 ± 0.8611363 0.3478548 ± 0.3399810 0.6521452 5 ± 0.9061799 0.2369269 ± 0.5384693 0.4786287 0 0.5688889 §8 数值微分 定理4 设 则有 （4.7） 其中系数 与 无关. 定理4表明 是 阶， 代替 ，有 （4.8） 若用4乘(4.8)式，减去(4.7)式再除3记之为 在(4.7)中，若用 则得 这里 以及后面将出现的 均为与 无关的系数， 这样构造的 与积分值 近似的阶为 . （4.9） 比较(4.9)与(4.4)可知，这样构造的序列 就是辛普森公式序列 又根据(4.9)，有 若令 则又可进一步从余项展开式中消去 项，而有 这样构造出的 ，其实就是柯特斯公式序列. 它与积分值 的逼近阶为 如此继续下去，每加速一次，误差的量级便提高2阶. 一般地，若记 则有 （4.10） 经过 次加速后， （4.11） 上述处理方法通常称为理查森外推加速方法. 设以 表示二分 次后求得的梯形值，且以 表示 序列 的 次加速值，则依递推公式(4.10)可得 余项便取下列形式： 计算过程：  (1) 取 令 ( 记区间 的二分次数). (2) 求梯形值 即按递推公式(4.1)计算 (3) 求加速值，按公式(4.12)逐个求出T表(见表4-5)的 第 行其余各元素 公式(4.12)也称为龙贝格求积算法. （4.12） 求 (4) 若 (预先给定的精度)，则终止计算， 否则令 转(2)继续计算. 并取 可以证明，如果 充分光滑，那么T表每一列的元素及对角线元素均收敛到所求的积分值 ，即 对于 不充分光滑的函数也可用龙贝格算法计算， 只是收敛慢一些，这时也可以直接使用复化辛普森公式计算. 例4 解 在 上仅是一次连续可微， 用龙贝格算法计算积分 用龙贝格算法计算结果见表4-6. 从表中看到用龙贝格算到 的精度与辛普森求积 精度相当. 这里 的精确值为 4.6 高斯（Gauss）型求积公式 4.6.1 高斯积分问题的提出 在前面建立牛顿-柯特斯公式时，为了简化计算，对插值公式中的节点限定为