代数精度 Chapter 3 Numerical Integreation and Differentiation 对于机械求积公式fx)d=∑Af(x) i=0 R(f=I-∑Af(x) i=0 定义若上述公式对所有不超过m次的多项式Pm(x)都精确成立, 即R(Pm)=0而对某一个m+1次多项式Pm+x近似成立 即Rn(Pm+)≠0则称机械求积公式具有m次代数精度 梯形公式1(+的代数精度为1 HUST
代数精度 n b i i a i 0 f(x)dx Af(x ) = 对于机械求积公式 ∫ ≈ ∑ 定义 若上述公式对所有不超过m次的多项式Pm(x)都精确成立, 即Rn(Pm)=0,而对某一个m+1次多项式Pm+1(x)近似成立, 即Rn(Pm+1) 0.则称机械求积公式具有m次代数精度. 梯形公式 的代数精度为 梯形公式 的代数精度为1. b a (b a) f(x)dx [f(a) f(b)] 2− ≈ + ∫ = − ∑ = n I A f(x ) i i i R f) 0 ( n
判断代数精度的方法 Chapter 3 Numerical Integreation and Differentiation 当f(x)=1,Xx2,…Ⅻm时求积公式精确成立, 而f(x)=xm+时公式近似成立分求积公式的代数精度为m次 证明:必要性显然,下证充分性 对任意m次多项式Pn(X)=a+aX+a2x2+…anx(an≠0) 由于求积公式!减=∑Ax)对于f(x)=1,xx2…xm时精确成立 b ∫1dx=∑A b k xdX=∑Ax X"dx=∑AkxR a k≠0 k=a,+a×+…+21xk=a2A+a2Ax+…+a2∑A k=0 k=0 k=0 ∑A+aX+…+anx)=∑APnx) k=0 ∴求积公式对Pn(X)精确成立 但对m+1次多项式公式近似成立(R≠0),由定义知 HUST 该公式的代数精度是m次
当f(x)=1,x,x2,…,xm时,求积公式精确成立, 而f(x)= xm+1时公式近似成立, 证明: 必要性显然.下证充分性 对任意m次多项式Pm(x)=a0+a1x+ a2x2+…+ amxm.(am 0) 由于求积公式 对于 由于求积公式 对于 f(x)=1,x,x2,…,xm时精确成立 n b i i a i 0 f(x)dx Af(x ) = ∫ ≈ ∑ ≠≠ ≠ ∴= = = ∫∫ ∫ ∑∑ ∑ nn n bb b m m k kk kk aa a k0 k0 k0 1dx A , xdx A x , ... , x dx A x = + ++ ∫ ∫∫ ∫ b bb b m m 01 m a aa a P (x)dx a dx a xdx ... a x dx 判断代数精度的方法 求积公式对Pm(x)精确成立. 但对m+1次多项式,公式近似成立(R 0),由定义知. 该公式的代数精度是m次 == = = + ++ ∑∑ ∑ nn n m 0 k 1 kk m kk k0 k0 k0 a A a A x ... a A x = = + + + ∑ n m k 0 1k mk k 0 A (a a x ... a x ) = = ∑ n km k k 0 A P (x ) !求积公式的代数精度为m次
Chapter 3 Numerical Integreation 验证梯形公式的代数精度为1 and Differentiation b-a 解梯形公式x)d [f(a+f(b) 令f(x)=1左=[=b-a右 a [1+1]=b-a左=右 公式对f(x)=1精确成立。 b-a 令(x)=x×x=12,右=2+]=2,在=右 公式对f(x)=X精确成立 令f()=X2左=「x2d 右 a [a2+b2]≠左 2 公式对f(x)=×2不再精确成立 梯形公式代数精度为1 例 Simpson公式的代数精度为3 f(x)、b a tb [f(a)+4f(--)+f(b)] 2 HUST
例 验证梯形公式的代数精度为 1. 解:梯形公式 b a b a f(x)dx [f(a) f(b)] 2 − ≈ + ∫ 令f(x)=1 b a b a 1dx b a, [1 1] b a, 2 左 右 左右 − = =− = + =− = ∫ 公式对 f(x)=1精确成立. 令f(x)=x 22 22 b a b a ba b a xdx , [a b] , 22 2 右 左右 −− − = = += = ∫ 公式对 f(x)=x精确成立 令f(x)=x 2 3 3 b 2 22 a b a ba x dx , [a b ] 3 2 左 右左 − − = = = +≠ ∫ 公式对 f(x)=x 2不再精确成立 梯形公式代数精度为 1. 例 Simpson公式的代数精度为 3 b a ba ab f(x)dx [f(a) 4f( ) f(b)] 6 2 − + ≈ ++ ∫
Chapter 3 Numerical Integreation vr设有求积公式∫f(x)dx=Af(-1)+Af()+A(1) and Differentiation 试确定系数AA1A2使这个公式具有最高的代数精度 分析:要确定公式中3个待定常数A0A1A2, 可令公式对1Xx2都准确成立 解:令f()=1,XⅩ2公式都准确成立,则 A+A1+A2=2 A+A,=0 解得A=1/3A1=43A2=13 A+A ∷该求积公式为∫,fxkx=f(1)+4(0)+f( 3 易验证:f(×)=3时,求积公式也准确成立 而fx)=x时xdk= [(-1)+4×0+1] 该求积公式具有3次代数精度,它是[-1,1上的 Simpson公式 HUS
例 设有求积公式 1 0 12 1 f(x)dx A f( 1) A f(0) A f(1) − ≈ −+ + ∫ 分析: 要确定公式中 3个待定常数 A 0,A 1,A 2 可令公式对 1 ,x,x 2都准确成立. 解:令f(x)= 1,x,x 2.公式都准确成立 , 则 012 0 2 0 2 AAA2 AA0 2 A A 3 ++= − + = + = 解得 A 0=1/3, A 1=4/3, A 2=1/3 该求积公式为 1 1 1 f(x)dx [f( 1) 4f(0) f(1)] − 3 ≈ −+ + ∫ 易验证:f(x)= x 3 时, 求积公式也准确成立 而f(x)= x4 时 1 4 44 1 221 x dx [( 1) 4 0 1 ] − 533 ≈ ≠ = − +×+ ∫ 该求积公式具有 3次代数精度 ,它是[-1,1]上的Simpson公式. 试确定系数 A 0,A 1,A 2,使这个公式具有最高的代数精度
Chapter 3 Numerical Integreation and Differentiation 一般对于n+1给节点上的机械求积公式∫x)dx=∑A(x) k=0 若使其代数精度至少为n,则可确定A构造出求积公式 只需令上式对f(x)=1XX2,…,x都准确成立,则 An+A,+A2+…+A=b-a AnX+AX1+…+AnXn 2 +1 x+AX1+…+Anxn n+1 上面是关于A0,A1,A2…,An的线性方程组, 其系数行列式为范德蒙行列式其值非零, 可求得唯一解 HUST
一般,对于n+1给节点上的机械求积公式 n b k k a k 0 f(x)dx A f(x ) = ∫ ≈ ∑ 若使其代数精度至少为n,则可确定Ak,构造出求积公式. 只需令上式对f(x)=1,x,x2,…,xn都准确成立,则 012 n 2 2 00 11 nn n1 n1 nn n 00 11 nn A A A ... A b a b a A x A x ... A x 2 (3.4) ... ... b a A x A x ... A x n 1 + + + + + + = − − + ++ = − + ++ = + 上面是关于A0 , A1 , A2 ,…,An的线性方程组, 其系数行列式为范德蒙行列式,其值非零, 可求得唯一解