54 Gauss求积公式 对于机械求积公式 (x)k=∑4/(x)+/(541) k=0 略去余式R/],由定理52知,它如果是插值型求 积公式,则至少有n次代数精度
5.4 Gauss求积公式 对于机械求积公式 (5.4.1) 略去余式 R f [ ],由定理 5.1.2 知,它如果是插值型求 积公式,则至少有 n 次代数精度。 ( ) ( ) [ ] 0 f x dx A f x R f k b a n k = k + =
从式(541)可以看出,待定参数 4,x(k=0,1…n)共2n+2个,如果令 ∫(x)=,x2…x2使求积公式(41)精确成 立,即R[/0,从理论上看,建立了求 Ax2(k=0,1,…m)的2m+2个代数方程组
从 式 (5.4.1) 可 以 看 出 , 待 定 参 数 , 0,1, ( ) A x k n k k = 共 2 2 n+ 个,如果令 ( ) 2 2 1 1, , , , n f x x x x + = 使求积公式(5.4.1)精确成 立,即 R f [ ]= 0 ,从理论上看,建立了求 ( ) , 0,1, A x k n k k = 的 2 2 n+ 个代数方程组
+1(b11-a“)(=0,1,2,n)(542) 如果从式(542)解得A,x(k=0,1…n),则求积 公式(541)将有2n+1次代数精度 下面我们就一个例子来看:
(5.4.2) 如果从式(5.4.2)解得 , 0,1, ( ) A x k n k k = ,则求积 公式(5.4.1)将有2 1 n+ 次代数精度。 下面我们就一个例子来看: ( ) ( 0,1,2,... ) 1 1 1 1 0 b a j n j A x j j n k j k k − = + = + + =
例541试决定参数A,A2x02x,使求积公式 f(x)dx= Aof(xo)+Af(1) 具有3次代数精度。 解令f(x)=1x,x2,x3,代入求积公式得 A+41=2 Axo+AX=0 (2) Aox+ Ax (3) (4)
例 5.4.1 试决定参数 0 1 0 1 A A x x , , , ,使求积公式 具有 3 次代数精度。 ( ) ( ) ( ) 0 0 1 1 1 1 f x dx = A f x + A f x − 解 令 ( ) 2 3 f x x x x = 1, , , ,代入求积公式得 + = + = + = + = 0 (4) (3) 3 2 0 (2) 2 (1) 3 1 3 0 2 1 2 0 0 0 1 1 0 1 0 1 0 1 A x A x A x A x A x A x A A
由式(1)及式(2)可解得 2. 由式(3)及式(4)可解得 2x Ao=3x(x1-xo 2
由式(1)及式(2)可解得 由式(3)及式(4)可解得 (5) 2 2 1 0 0 1 1 0 1 0 − − = − = x x x A x x x A (6) 3 ( ) 2 3 ( ) 2 1 0 2 0 1 1 0 2 1 0 1 0 − − = − = x x x x A x x x x A