高等数值分析.doc

数值积分方法小述 一、背景 数值积分方法发展的前提是在17世纪以牛顿和莱布尼茨为首的一批数学家发展起来的微积分。在最初的研究中，求解积分的方法便是找到求解原函数的方法，得到原函数，以此为基础解决其他问题。但是在深入的研究中，逐渐发现一些函数的原函数求解极其困难，甚至无法表示出来，是超越函数，还有的根本没有原函数，比如对于延拓函数： 无法求出它的原函数，这时要求它的积分就无法使用牛顿-莱布尼茨公式了，解决积分的问题便受到阻碍。这种情况下就需要寻求一种新的求积分的方法来解决这些问题了。数值积分方法便在数学家们的需求下发展起来。 二、发展历程 等距节点的多项式插值求积法的观点最早是1676年出现在Newton给Leibniz的一封信中。1711年，Cotes在总结了牛顿的观点后，系统归纳了小于10个节点的插值求积方法，并发表了一篇相关论文。1743年，Simpson发表他所研究的求积方法。但是从历史上看，对于辛普森的方法，数学家Cavalieri和Gregory似乎研究的更早，而且Cotes也早就得到了这种方法。1814年，数学王子Gauss在研究这个问题时，通过优化那些求积节点得到一种更高精度的数值求积分方法，随后便发表了他的第一篇关于数值求积分的论文。时间过了100多年，数学家Fejer于1933年，将Chebyshev点作为节点应用于数值求积分中，得到了一种新的方法。1960年，数学家Clenshaw和Curtis研究得到一种更为高效的数值求积公式。Kronrod在1964年发表了他自己的数值求积方法，4年后的Patterson对这种方法进行了推广，得到的方法也为世人所知。值得一提的是Richardson在1927年发现的外推法，当时并没有用来做数值积分问题。而数学家Romberg在1955年将它应用到数值积分上，取得不小的成果。 三、五种数值积分方法 对于不能求得其原函数的函数，要计算 就要研究计算定积分的近似方法，即数值积分法。 若能已知在部分点上的函数值，利用拉格朗日插值，可以构造一个多项式来逼近被积函数，而多项式为被积函数，在区间上的定积分是容易计算的，这样就可以得到计算定积分的一种数值积分方法，即 等式右边是多项式的积分，所以可以得到相应的结果 （1.0） 其中为积分后得到的相应系数。 利用这种思想可以得到插值型求积公式。 1、Newton-Cotes公式 将积分区间划分为等分，步长为，选取等距节点（）构造插值型求积公式。由于此处是变量，对公式（1.0）进行稍微变化，写成 称为Newton-Cotes公式，其中是Cotes系数，。令，则有 Cotes系数的计算都是多项式的积分，没有实质性的困难。 a、当时，解得Cotes系数为 此时 这就是梯形求积公式。在平面几何中，是以为顶的曲边梯形的面积，而是以曲线为斜边的梯形面积。因此，梯形公式就是以梯形的面积来近似代替以为顶的曲边梯形面积。 若函数在上具有连续的二阶导数，则梯形公式的截断误差为 其中。梯形公式只有1次代数精度，即对所有一次函数精确成立。 b、当时，解得Cotes系数为 ， 此时 这就是Simpson求积公式。从几何上看是以抛物线为顶点的曲边梯形面积来近似代替为顶的曲边梯形面积，所以Simpson公式也称为抛物线求积公式。 若函数在上具有连续的四阶导数，则Simpson公式的截断误差为 其中。Simpson公式有3次代数精度，对所有3次以下的函数精确成立。 c、当时的Newton-Cotes公式则特别的称为Cotes公式 其中。Cotes公式有5次代数精度。 2、复合求积公式 将积分区间划分为等分，，，，在每个子区间上分别采用梯形公式就得到了复合梯形公式 复合梯形公式的余项为 在每个子区间上分别采用Simpson公式就等到复合Simpson公式 复合Simpson公式的余项为 3、Romberg求积公式 设，从梯形公式出发，其余项为，其中只有是变量，则可设 于是 （3.1） (3.2) 用4乘（3.2）式减去（3.1）式再除以3后所得的式子记为，则有 （3.3） 于是 （3.4） 用16乘（3.4）式减去（3.3）式后除以15，将所得式子记为 事实上对于复合梯形公式、复合Simpson公式、复合Cotes公式有 根据上述思想，依次递推可以得到一个统一的公式 （3.5） 此种处理方法为Richardson外推加速方法。