数据积分-牛顿科茨法与高斯勒让德法对比及示例

在结点x _i上插值f的次数最多是n次的多项式

p(x)=∑f(x _i)l _i(x)

拉格朗日插值多项式:

l _i(x)= ∏(x-x _j)/(x _i-x _j)

 

基于插值的数值积分：

∫f(x)dx≈∫p(x)dx=∑f(x _i)  ∫l_i(x) dx =∑A_if(x_i)

当 结点等距即为牛顿 -科茨公式

 

例：设插值多项式是2次，积分区间[0,1],设f(x) ≈p(x)=c ₀+c ₁x+c ₂x ²

求牛顿-科茨积分公式。

n=2;结点(0,1/2,1)

对结点0,1/2,1的三个基本多项式是:

l ₀=(x-1/2)(x-1)/((0-1/2)(0-1))=2(x-1/2)(x-1)

l ₁=-4 x(x-1)

l ₂=2 x(x-1/2)

则A ₀= ∫l₀(x)=1/6；

A ₁= ∫l₁(x)=2/3；

A ₂= ∫l₂(x)=1/6

所以：∫f(x)dx≈ 1/6 f(0)+2/3 f(1/2)+1/6 f(1)

 

也可以用 待定系数法求解

设∫f(x)dx≈ A₀ f(0)+ A₁ f(1/2)+ A₂ f(1)

它对于所有次数 ≤2 的多项式是精确成立的。

依次把多项式f(x)=1,x,x ²作为试用函数，得到

1=∫f(x)dx=∫1dx=  A₀+ A₁ + A₂

1/2=∫f(x)dx=∫xdx=  1/2 A₁ + A₂

1/3=∫f(x)dx=∫x ²dx=  1/4 A₁ + A₂

这三个方程的联立方程组的解是：

A ₀=1/6；

A ₁=2/3；

A ₂=1/6

因为公式是线性的，所以对任何二次多项式

f(x) =c₀+c₁x+c₂x² ，它将产生积分的精确值。

 

高斯求积公式

ω(x) 为给定的正的权函数的积分公式：

∫f(x)ω(x)dx≈∑ A_if(x_i)

A_i= ∫ω(x)l_i(x) dx

 = ∫ω(x) ∏(x-x_j)/(x_i-x_j) dx

如果预先限制结点x _i，如上述拉格朗日插值积分及牛顿科茨积分，

可以找到次数 ≤n的多项式是精确成立的求积公式；

如果 不限制结点，因为有n+1个系数A _i和n+1个结点x _i，

则可以找到次数 ≤2n+1的多项式是精确成立的求积公式；

高斯求积定理：

设 ω是正的权函数，q是一个n+1次非零多项式并且与∏ _n是 ω 正交的，也就是对 任意 pє∏_n （ n 阶多项式空间），都有 ∫q(x) p(x)ω(x) dx=0

若x ₀,x ₁,…,x _n是q的零点，则具有 A_i= ∫ω(x) ∏(x-x_j)/(x_i-x_j) dx

公式中给定系数 A_i的求积公式 ∫f(x)ω(x)dx≈∑ A_if(x_i)

对于所有 fє∏_2n+1是精确成立的。

 

也就是说：与 ∏_n关于 ω正交的n+1次非零多项式的零点x ₀,x ₁,…,x _n作为积分结点的话，可以找到次数 ≤2n+1的多项式是精确成立的求积公式  ∫f(x)ω(x)dx≈∑A_if(x_i)

f(x)= c ₀+c ₁x+c ₂x ²+...+c _2n+1x ²ⁿ⁺¹

积分结点x ₀,x ₁,…,x _n叫做 高斯点。

 

权函数 ω(x)=0 是一种特别重要的情况，这时，如果积分区间取[-1,1]，则高斯积分公式变为

∫_-1 ¹ f(x) dx≈∑₀ ⁿ A_if(x_i)

以高斯点为零点的n+1次多项式q _n+1(x)=(x- x ₀) (x- x ₁)… (x- x _n)

称为 勒让德(Legendre)多项式

 

高斯-勒让德积分，对2个高斯点，2阶勒让德多项式来说，q ₂(x)=1/2(3x ²-1);  n=1

2个高斯点为：√3/3，-√3/3

所以

∫_-1 ¹ f(x) dx≈∑₀ ⁿ A_if(x_i)= A₀f(√3/3 )+ A₁f(-√3/3)

 

例; 设插值多项式是2n+1=3次，积分区间 [-1,1],设f(x) ≈p(x)=c ₀+c ₁x+c ₂x ²+c ₃x ³

求高斯-勒让德积分公式。

采用 待定系数法求解

设∫f(x)dx≈ A₀ f(-√3/3)+  A₁ f(√3/3)

它对于所有次数 ≤3 的多项式是精确成立的。

依次把多项式f(x)=1,x作为试用函数，得到

2=∫f(x)dx=∫1dx=  A₀+ A₁

0=∫f(x)dx=∫xdx=  -√3/3  A₀ +√3/3  A₀

这三个方程的联立方程组的解是：

A ₀=1；

A ₁=1

所以两点高斯-勒让德积分公式为：

∫_-1 ¹ f(x)dx≈f(-√3/3)+  f(√3/3)

 

分别令f(x)=x ², x ³, x ⁴

f(x)=x ² 时， ∫_-1 ¹ f(x)dx=2/3；f(-√3/3)+  f(√3/3)=2/3  精确成立

f(x)=x ³ 时， ∫_-1 ¹ f(x)dx=0；f(-√3/3)+  f(√3/3)=0      精确成立

f(x)=x ⁴ 时， ∫_-1 ¹ f(x)dx=2/5；f(-√3/3)+  f(√3/3)=2/9  不精确成立

验证了两点高斯-勒让德积分公式对 次数 ≤3 的多项式是精确成立的。