计算方法 数值积分.docVIP

三、复合公式 1. 复合梯形公式 ∫ ≈ a b f ( xdx ≈ ∑ h [ f (xk + f ( x k+1] k =0 2 + 2 f ( xn?1 + f ( x n ] n?1 h [ f (x0 + 2 f ( x1 + 2 R=? 误差 (b ? a 2 h f ′′(η 12 ； 2. 复合辛浦生公式 h= b?a ,n =2m 为偶数, n f ( xdx ≈ ∑ h k =0 3 m?1 ∫ a b [ f2k +4f 2k+1+f 2k+2] ≈ h ( f0 +4f1 +2f2 +4f3 + +2f2m?2 +4f 2m?1+ f2m 3 R=? 误差 (b ? a 4 (4 h f (η 180 1 0 。 例：用复合公式计算 ∫ 4arctgx dx ，精确值 1 x=0 2 4 xarctgx a ln(1 x ? =? ? + ? ? = π ? 2ln 2 = 1.7552983 由 f ( x = 4arctgx ，将区间[0,1]分成 4 等份，n=4，h=0.25，得数表： 31 x 0 0.25 0.5 0.75 1 f (x 0 0.9799146 1.8545904 2.5740044 3.1415927 用复合梯形公式: h T4 = ( f0 + 2 f1 + 2 f2 + 2 f3 + f4 = 1.7448264 ，误差 2 R=0.0104719； 用复合抛物线公式： h S4 = ( f0 + 4 f1 + 2 f2 + 4 f3 + f4 3 四、高斯型求积公式 1. Legendre 多项式 =1.7553750， 误差 R=0.0000767。 1 dn ?( x2 ?1n ? Pn ( x = n n ? ?, 2 ? 2! dx n = 0,1,2, 对于任意次数不超过 n 的多项式 q(x，成立 ∫ 勒让德多项式 1 ?1 q( xPn+1 ( xdx = 0 。 2. 高斯—勒让德求积公式 Pn+1 ( x 的零点是高斯—勒让德求积公式的求积节点。 n k =0 ∫ b a f ( xdx ≈ ∑ Ak f ( xk 当 n=1 时 ? 1 ? ? 1 ? f ( x dx f f ≈ ? + ? ? ? ?， ∫ ?1 ? 3? ? 3? 1 其代数精确度为 3。 ? 3 ?? 1? ? 3? 当 n=2 时, ∫ ?1 f ( xdx ≈ ?5 f ? ?? 5 ? ? + 8 f ( 0 + 5 f ? ? 5? ?? 9? ? ? ?? ? ? ? 1 其代数精确度为 5。 一般地，n+1 个求积节点的高斯型求积公式的代数精确度为 2n+1。 32 对于一般的有界区间[a,b]上的定积分， 可以通过变量代换转化为区间[-1,1] 上的定积分。即在积分 ∫ b a f ( xdx 中，令 x = b?a b+a t+ 。 2 2 例：用 n=2 的高斯—勒让德公式计算 ∫ 1 0 4arctgx dx 。 1 x=0 2 4 xarctgx a ln(1 x ? ? + 精确值 = ? ? ? = 1.7552983 。 1 1 = x t + ，则 此处，a=0，b=1，作变换 2 2 1 1 1 1 = 4 arctgx dx 2 arctg ( t + dt ∫0 ∫ ?1 2 2 1 1 f ( t = 2 arctg ( t + ，故 f(0=0.9272652， 于是 2 2 f (? 3 = 0.2244562 5 ， f( 3 = 1.4515062 5 ，所以原式 ? 3 ?? 1? ? 3? ≈ ?5 f ? ?? 5 ? ? + 8 f ( 0 + 5 f ? ? 5? ?? =1.7553526，误差 R=0.0000543。 9? ? ? ?? ? ? ? 三个点的高斯型公式比五个点的复合抛物线公式还要精确，误差还要小。 33