第四章 数值积与数值微分.pptVIP

华长生制作 4.4.1 梯形法的递推化 由上节讨论得知加密节点可以提高求积公式的精度，复化 求积方法对提高精度是行之有效的，但必须事先给出合适 的步长（即n的选取），这又是一个难题。 在电子计算机上通常采用把区间逐次二分，反复利用复化 求积公式进行计算，直到二分前后两次积分近似值之差符 合精度要求为止。 各节点为 复化梯形公式为 --------(1) 经过二分只增加了一个分点 --------(3) --------(2) 用复化梯形公式求得该子区间上的积分值为 这里h仍为二分前的步长.将每个子区间上的积分值相加得 由(1)(2)两式可 (3)式称为递推的梯形公式 递推梯形公式加上一个控制精度,即可成为自动选取步长的复化梯形公式 优点：梯形法计算简单 缺点：收敛慢，为了达到要求的精度，需要二分区间 多次，分点大量增加，计算量很大 4.4.2 龙贝格算法 根据复化梯形公式的余项表达式可知 假定 ，则有 即 依此类推 这样直接用计算结果来估计误差的方法通常称作误差的 事后估计法. 用积分近似值 的误差作为 的一种补偿，得到 例1. 原积分的精确值为 2位有效数字 3位有效数字 6位有效数字 将 用于计算 ò = 1 0 sin dx x x I 同理由复化辛普森公式的余项 可得 直接验证易知 由复化Cotes公式的余项 得 称为龙贝格公式，是一种加速技术 0.9460831 0.9460831 0.9460833 0.9456909 3 0.9460830 0.9460869 0.9445735 2 0.9461459 0.9397933 1 0.9207355 0 例3 将以上三个加速公式用于求 从表中可以看出三次加速求得R1=0.9460831每位数字都是 有效数字 4.4.2 理查德森外推加速法 定理4 设 ，则有 其中系数 与h无关 用h/2代替h，有 由(4.7)及(4.8)两式可得 比较(4.9)与 可知，这样构造的序列 就是辛普森公式序列 一般的，若记 ，则有 上述处理方法称为理查森外推加速法 设以 表示二分k次后求得的梯形值，且以 表示 的m次加速值，则依递推公式(4.10)可得 公式(4.12)也称为龙贝格求积算法 表4-5 T表 本章作业 P159 8(1)(3) 第四章 数值积分与数值微分 4.5 高斯求积公式 牛顿—柯特斯型求积公式是封闭型的（区间[a,b]的两端点a, b均是求积节点）而且要求求积节点是等距的,受此限制，牛顿—柯特斯型求积公式的代数精度只能是n(n为奇数)或n+1(n为偶数)。而如果对求积节点也适当的选取,即在求积公式中不仅Ak而且xk也加以选取,这就可以增加自由度，从而可提高求积公式的代数精度。 4.5.1 一般理论 定义4： 上式称为梯形求积公式,也称两点公式，记为 求积公式为 梯形公式的余项为 2.辛普森公式及其余项 柯特斯系数为 求积公式为 上式称为辛普森求积公式，也称三点公式或抛物线公式 记为 辛普森公式的余项为 辛普森公式具有3次代数精度 3.柯特斯公式及其余项 柯特斯系数为 求积公式为 上式称为柯特斯求积公式，也称五点公式 记为 柯特斯公式的余项为 柯特斯公式具有5次代数精度 考察柯特斯系数 因此用牛顿-柯特斯公式计算积分的舍入误差主要由 其值可以精确给定 §4.2.4 牛顿-柯特斯公式的稳定性(舍入误差) 记 而理论值为 牛顿-柯特斯公式的舍入误差只是函数值误差的 此时,公式的稳定性将无法保证 在实际应用中一般不使用高阶牛顿-柯特斯公式 例: 用n=2和n=3的牛顿-柯特斯公式 解： 求 的近似值。 1.n=2时 2.n=3时 （ 的真实值为0.7668010） 本章作业 P158 2(3)(4)(用辛普森公式计算)，3，4 第四章 数值积分与数值微分 4.3 复化求积法 4.3 复化求积法 直接使用牛顿-柯特斯公式余项将会较大 当n8时，公式的舍入误差又很难得到控制 此时，使用复化方法， 然后在每个小区间上使用低阶牛顿-柯特斯公式，最后将每个小区间上的积分的近似值相加，这种方法称为复化求积法 各节点为 记为 4.3.1 复化求积公式 由积分的区间可加性 复化求积公式 复化梯形公式分解 复化辛普森公式分解 例1. 解: 依次使用8阶复化梯形公式、4阶复化辛普森公式和2阶复化柯特斯公式 给出n=8的函数表 0 1 0.125 00.25 00.375 00.5