计算方法ch3.pptVIP

? 龙贝格(Romberg)算法 Step 1: 用梯形公式计算T1； Step 2: 将[a，b]二分，用半分公式求出T2，再用梯形加速公式求出S1； Step 3: 再将新区间二分，用半分公式求出T4，用梯形加速公式求出S2，再用Simpson加速公式求出C1； Step 4: 再将新区间二分，用半分公式求出T8，用梯形加速公式求出S4，再用Simpson加速公式求出C2，再同Romberg公式求出R1； * Chapter 3 数值积分与数值微分 (Newton -Leibniz) 定义 若存在实数xi?[a, b]及 Ai , i=0,1, ???, n, 且任取 f (x)?C[a, b], 都有 式(?)称为一个数值求积公式。 ----求积系数; 与 xi有关, 但与 f (x)无关; ----求积节点; ----求积余项 或求积公式(?)的截断误差. ? 插值型求积公式 设 f (x)在一组节点 上的函数值为 f (xk), k=0,1,???, n, 做n次Lagrange插值多项式 式中 为n次插值基函数。 ----插值型求积公式 令 由Lagrange插值余项有求积截断误差为: 式中 且与x有关。 1 Newton-Cotes求积公式 只与n有关。 (1) n=1 (2) n=2 ----梯形公式 ----Simpson公式(抛物线公式) 部分Cotes系数表 840 41 35 9 280 9 105 34 280 9 35 9 840 41 6 288 19 96 25 144 25 144 25 96 25 288 19 5 90 7 45 16 15 2 45 16 90 7 4 8 1 8 3 8 3 8 1 3 6 1 6 4 6 1 2 2 1 2 1 1 ) ( n k C n Cotes系数有两个基本性质: 对称性: 和为1: ----Simpson 3/8公式 ----Cotes公式 (3) n=3 (4) n=4 求积公式的代数精度 定义 若求积公式 对于xi (i=0,1,???, m)均能准确成立, 但对于x m+1不准确成立, 则称该求积公式具有m次代数精度。 例 求梯形公式 梯形公式具有1次代数精度； 同理可得Simpson公式具有3次代数精度。 的代数精度。 例 确定下列公式中的参数，使其代数精度尽量高，并指明所得公式的代数精度。 定理 对于Newton-Cotes求积公式， 当n为偶数时，至少具有n+1次代数精度。 Newton-Cotes公式的余项： 1. 梯形公式的余项 设 f ??(x)在[a, b]上连续, 则梯形公式的余项为 当n为奇数时，至少具有n次代数精度； 2. Simpson公式的余项 设 f (4)(x)在[a, b]上连续, 则Simpson公式的余项为 3. Cotes公式的余项 设 f (6)(x)在[a, b]上连续, 则Cotes公式的余项为 例 用梯形公式、Simpson公式及Cotes公式计算积分 ，并估计误差。 ? 复合求积法 1. 复合梯形公式 2. 复合Simpson公式 3. 复合Cotes公式 例 利用复合求积公式计算积分 复合求积公式的截断误差 1. 复合梯形公式的误差 设 f ??(x)在[a, b]上连续, 则复合梯形公式的余项为 2. 复合Simpson公式的误差 设 f (4)(x)在[a, b]上连续, 则复合Simpson公式的余项为 3. 复合Cotes公式的误差 设 f (6)(x)在[a, b]上连续, 则复合Cotes公式的余项为 例 考虑用复合Simpson公式计算积分 要使误差小于 ，那么求积区间[0，1]应分多少个子区间？并以此计算积分近似值。 ? 误差的事后估计方法 1. 复合梯形公式 所以若用T2n作为I的近似值，则其截断误差约为 否则将每个小区间再次分半进行计算，直到满足要求为止。 2. 复合Simpson公式 则其截断误差约为 3. 复合Cotes公式 则其截断误差约为 例 利用梯形法公式计算积分 的近似值，使误差不超过 。 ? 梯形法的步长逐次分半算法 例 用梯形法的步长逐次分半算法T16计算积分 的近似值。 * *