Gauss型求积公式及其Matlab程序

为什么要引入数值积分

由于Gauss型求积公式属于数值积分的内容，学东西总要知道它的来龙去脉，下面我简单介绍一下为什么要引入数值积分

给定函数$f(x)\in C[a,b]$,考虑积分

$$I(f)={\int }_a^{b}f(x)dx$$
的计算问题，从数学分析中知道，当已知$f(x)$的原函数为$F(x)$时，由牛顿-莱布尼兹公式，有
$$I(f)={\int }_a^{b}f(x)dx=F(b)-F(a)$$
然而，在实际计算中，被积函数的$f(x)$的原函数经常无法用初等函数表示，但过于复杂。还有时，$f(x)$只在一些离散点上给出。在这样的情况下，就有必要借助数值方法来求$I(f)$的近似值。

Gauss型求积公式的定义

从正交多项式理论可知，在区间[a,b]上，对给定的权函数$\rho (x)$,存在正交多项式系${{\omega }_k(x)}_{k=0}^{\infty }$并且可以将其构造出来。进一步，已经证明了${\omega }_k(x)$在[a,b]上恰有k个相异的根，取${\omega }_{n+1}(x)$的n+1个根为求积结点，构造插值型求积公式，将得到具有2n+1次代数精度的求积公式，称如此构造出来的求积公式为Gauss型求积公式即：

如果求积结点$x_0,x_1,....,x_n$,使插值型求积公式${\int }_{-1}^{1}f(x)dx\approx$