计算方法（潘林强）NA03b.pptVIP

* §2 复合求积 /* Composite Quadrature */ Haven’t we had enough formulae? What’s up now? Oh come on, you don’t seriously consider h=(b?a)/2 acceptable, do you? Why can’t you simply refine the partition if you have to be so picky? Don’t you forget the oscillatory nature of high- degree polynomials! Uh-oh 高次插值有Runge 现象，故采用分段低次插值 ? 分段低次合成的 Newton-Cotes 复合求积公式。 ? 复合梯形公式： 在每个 上用梯形公式： = Tn /*中值定理*/ §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 注意到区间再次对分时 可用来判断迭代 是否停止。 HW: p.87 #8 §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 即： HW: p.87 #10 §4 高斯型积分 /* Gaussian Quadrature */ 构造具有2n+1次代数精度的求积公式 将节点 x0 … xn 以及系数 A0 … An 都作为待定系数。令 f (x) = 1, x, x2, …, x2n+1 代入可求解，得到的公式具有2n+1 次代数精度。这样的节点称为Gauss 点，公式称为Gauss 型求积公式。 例：求