数值分析NA08b.pptVIP

* §2 Composite Quadrature ? 复化 Simpson 公式： 4 4 4 4 4 = Sn 注：为方便编程，可采用另一记法：令 n’ = 2n 为偶数， 这时 ，有 §2 Composite Quadrature ? 收敛速度与误差估计： 定义 　　 若一个积分公式的误差满足 且C ? 0，则称该公式是 p 阶收敛的。 ~ ~ ~ 例：计算 解： 其中 = 3.138988494 其中 = 3.141592502 运算量基本相同 §2 Composite Quadrature Q: 给定精度 ?，如何取 n ? 例如：要求 ，如何判断 n = ？ ? 上例中若要求 ，则 即：取 n = 409 通常采取将区间不断对分的方法，即取 n = 2k 上例中2k ? 409 ? k = 9 时，T512 = 3 S4 = 3.141592502 注意到区间再次对分时 可用来判断迭代 是否停止。 §3 龙贝格积分 /* Romberg Integration */ 例：计算 已知对于? = 10?6 须将区间对分 9 次，得到 T512 = 3由 来计算 I 效果是否好些？ 考察 = 3.141592502 = S4 一般有： Romberg 序列 ? Romberg 算法： ? ? ? ? ? ? … … … … … … ? T1 = ) 0 ( 0 T ? T8 = ) 3 ( 0 T ? T4 = ) 2 ( 0 T ? T2 = ) 1 ( 0 T ? S1 = ) 0 ( 1 T ? R1 = ) 0 ( 3 T ? S2 = ) 1 ( 1 T ? C1 = ) 0 ( 2 T ? C2 = ) 1 ( 2 T ? S4 = ) 2 ( 1 T §3 Romberg Integration ? 理查德森外推法 /* Richardson’s extrapolation */ 利用低阶公式产生高精度的结果。 设对于某一 h ? 0，有公式 T0(h) 近似计算某一未知值 I。由Taylor展开得到： T0(h) ? I = ?1 h + ?2 h2 + ?3 h3 + … ?i 与 h 无关 现将 h 对分，得： ( ) ( ) ( ) ... ) ( 3 2 3 2 2 2 2 1 2 0 + + + = - h h h h I T a a a Q：如何将公式精度由 O(h) 提高到 O(h2) ？ ... 4 3 2 1 1 2 ) ( ) ( 2 3 3 2 2 0 2 0 - - - = - - - h h I h T T h a a 即： §4 高斯型积分 /* Gaussian Quadrature */ 构造具有2n+1次代数精度的求积公式 将节点 x0 … xn 以及系数 A0 … An 都作为待定系数。令 f (x) = 1, x, x2, …, x2n+1 代入可求解，得到的公式具有2n+1 次代数精度。这样的节点称为Gauss 点，公式称为Gauss 型求积公式。 例：求 的 2 点 Gauss 公式。 解：设 ，应有 3 次代数精度。 ? + ? 1 0 1 1 0 0 ) ( ) ( ) ( x f A x f A dx x f x 代入 f (x) = 1, x, x2, x3 不是线性方程组，不易求解。 §4 Gaussian Quadrature 证明： “?” x0 … xn 为 Gauss 点, 则公式 至少有 2n+1 次代数精度。 对任意次数不大于n 的多项式 Pm(x)， Pm(x) w(x)的次数不大于2n+1，则代入公式应精确成立： 0 = 0 ? “?” 要证明 x0 … xn 为 Gauss 点，即要证公式对任意次数不大于2n+1 的多项式 Pm(x) 精确成立，即证明： 设 0 ? 　　　 x0 … xn 为 Gauss 点 ?