第4节6节高斯求积公式.ppt

1 § 6 高斯求积公式 2 1. 一般理论 求积公式 ? ? ? ? n k k k b a x f A x x f 0 ) ( d ) ( 含有 个待定参数 2 2 ? n ). , , 1 , 0 ( , n k A x k k ? ? 当 为等距节点时得到的插值求积公式其代数精度至少 为 次 . k x n 3 为具有一般性，研究带权积分 , d ) ( ) ( ? ? b a x x x f I ? 这里 为权函数， ) ( x ? 求积公式为 , ) ( d ) ( ) ( 0 ? ? ? ? n k k k b a x f A x x x f ? 为不依赖于 的求积系数 . ) , , 1 , 0 ( n k A k ? ? ) ( x f ) , , 1 , 0 ( n k x k ? ? 为求积节点， k k A x 及 适当选取 使其具有最高 次代数精度 4 定理 n+1 个节点的插值型求积公式的代数精度 不超过 2n+1. 证：令 g(x)=(x-x 0 ) 2 …(x -x n ) 2 次多项式不精确成立。 对 2 2 ) 0 ) ( ( 0 ) ( 0 d ) ( ) ( ) ( 0 ? ? ? ? ? ? ? ? ? n x g x g A I x x g x g I k n k k k n b a ? ? ? 5 定义 若 n+1 个互异节点的插值型求积公式 的代数精度达到 2n+1 次，则称此 n+1 个互 异节点为高斯点， 此求积公式为高斯型 求积公式。 6 根据定义要使 求积公式 具有 次代数精度，只要对 1 2 ? n ), 1 2 , , 1 , 0 ( ) ( ? ? ? n m x x f m ? . 1 2 , , 1 , 0 , d ) ( 0 ? ? ? ? ? ? n m x x x x A n k b a m m k k ? ? 当给定权函数 ，求出右端积分，则可解得 ) ( x ? ). , , 1 , 0 ( n k A x k k ? ? 及 令 精确成立， , ) ( d ) ( ) ( 0 ? ? ? ? n k k k b a x f A x x x f ? 7 ? ? ? ? ? ? ? ? ? ? ? ? ? . 9 2 1 3 1 0 3 0 A x A x 例 5 ). ( ) ( d ) ( 1 1 0 0 1 0 x f A x f A x x f x ? ? ? 解 令公式 (5.3) 对于 准确成立， 3 2 , , , 1 ) ( x x x x f ? ; 3 2 1 0 ? ? A A 试构造下列积分的高斯求积公式： 得 ; 5 2 0 0 0 0 ? ? A x A x ; 7 2 1 2 1 0 2 0 ? ? A x A x 8 由此解出 , 9 10 , 21 5 1 0 1 0 ? ? ? x x x x 从而 9 . 277556 . 0 , 389111 . 0 ; 289949 . 0 , 821162 . 0 1 0 1 0 ? ? ? ? A A x x 这样，高斯公式是 ) 821162 . 0 ( 389111 . 0 d ) ( 1 0 f x x f x ? ? ). 289949 . 0 ( 277556 . 0 f ? 由于非线性方程组较复杂， 通常 就很难求解 . 2 ? n 故一般不通过解方程 求 ， ) , , 1 , 0 ( n k A x k k ? ? 及 而从分析高斯点的特性来构造高斯求积公式 . . 1 2 , , 1 , 0 , d ) ( 0 ? ? ? ? ? ? n m x x x x A n k b a m m k k ? ? ). ( ) ( d ) ( 1 1 0 0 1 0 x f A x f A x x f x ? ? ? 10 定理 5 b x x x a n ? ? ? ? ? ? 1 0 是高斯点的充分必要条件是以这些节点为零点的多项式 ) ( ) )( ( ) ( 1 0 1 n n x x x x x x x ? ? ? ? ? ? ? 与任何次数不超过 的多项式 带权 正交， n ) ( x P ) ( x ? . 0 d ) ( ) ( ) ( 1 ? ? ? b a n x x x x P ? ? 证明 即 插值型求积公式的节点 必要性 . , ) ( n H x P ? 设 , ) ( ) ( 1 2 1 ? ? ? n n H x x P ? 则 , ) ( d ) ( ) ( 0 ? ? ? ? n k k k b a x f A x x x f ? 11 . d ) ( ) ( d ) ( ) ( ? ? ? b a b a x x x q x