完整word版谈谈高斯 勒让德公式推导过程word文档良心出品.doc

4章 数值积分与数值微分 4.1 引言 4.1.1 数值求积的基本思想 实际问题当中常常需要计算积分．有些数值方法，如微分方程和积分方程的求解，也都和积分计算相联系． 依据人们所熟知的微积分基本定理，对于积分 ． 的原函数，便有下列牛顿－莱布尼兹(Newton-Leibniz)只要找到被积函数公式： 但实际使用这种积分方法往往有困难，因为大量的被积函数，诸如等等，我们 是由测量或数值计算给出的一张数据表时，当找不到用初等函数表示的原函数；另外，牛顿－莱布尼兹公式也不能直接使用．因此有必要研究积分的数值计算问题． 内存在一点，成立积分中值定理告诉我们，在积分区间 的矩形面积恰等于所求曲边梯形的面积（图４－１）而高为．问就是说，底为 的具体位置一般是不知道的，因而难准确算出称为区题在于点的值．我们将 上的平均高度．这样，只要对平均高度提供一种算法，相应地便获得一种数间值求积方法． 的算术平均平均作为平均高度高度“和”的近似值，这样如果我们用两端点 导出的求积公式 （4.1.1） ．而如果改用区间中点的便是我们所熟悉的梯形公式（几何意义参看图４－２） ，则又可导出所谓中矩形公式“高度”（今后简称矩形公式）近似地取代平均高度 （4.1.2） 加权平均得到平均然后用，上适当选取某些节点我们可以在区间更一般地， 高度的近似值，这样构造出的求积公式具有下列形式： （4.1.3） 仅仅与节点的权．权称为求积系数，亦称伴随节点式中的选称为求积节点； 取有关，而不依赖于被积函数 的具体形式．这类数值积分方法通常称为机械求积，其特点是将积分求值问题归结为函数值的计算，这就避开了牛顿－莱布尼兹公式需要求原函数的困难． 4.1.2 代数精度的概念 数值求积方法是近似方法，为要保证精度，我们自然希望求积公式能对“尽可能多”的函数准确地成立，这就提出了所谓代数精度的概念． 但对于次多如果某个求积公式对于次数不超过的多项式均能准确地成立，定义１ 项式就不准确成立，则称该求积公式具有次代数精度． 不难验证，梯形公式（4.1.1）的矩形公式（4.1.2）均具有一次代数精度． 次代数精度，只要令它对于３）具有都一般地，欲使求积公式（4.１. 能精确成立，这就要求 (4.1.4) 为简洁起见，这里省略了符号中的上下标． 的等距分点作为节点，这时取如果我们事先选定求积节点求，臂如，以区间 至少具有(4.1.3)次代数精度．本章第解方程组(4.1.4)即可确定求积系数，而使求积公式２节介绍这样一类求积公式，梯形公式是其中的一个特例． 的代数问题．的求积公式，原则上是一个确定参数和为了构造出形如(4.1.3) 插值型的求积公式4.1.3 设给定一组节点 由于代数多项式．式）(2.9)（参见第２章作插值函数在节点上的值，且已知函数． 的原函数是容易求出的，我们取 作为积分的近似值，这样构造出的求积公式 (4.1.5) 通过插值基函数的积分得出 称为是插值型的，式中求积系数 (4.1.6) 由插值余项定理（第２章的定理２）即知，对于插值型的求积公式(4.1.5)，其余项 (4.1.7) 有关，．式中 与变量 的多项式是插值型的，按式(4.1.7)，对于次数不超过，其余项如果求积公式(4.1.5) 等于零，因而这时求积公式至少具有次代数精度． 至少具有次代数精度，则它必定是插值型的．事实上，这时公反之，如果求积公式(4.1.5) 对于插值基函数应准确成立，即有(4.1.5) 式 ，上式右端实际上即等于，因而式(4.1.6)成立．注意到 综上所述，我们的结论是： 的求积公式至少具有次代数精度的充分必要条件是，它是插值型的． 形如(4.1.5)定理１4.1.4 求积公式的收敛性与稳定性 定义２ 在求积公式(4.1.3)中，若 ． 其中，则称求积公式(4.1.3)是收敛的． 算计由中式积在即得实误产可能生差，际到，求公(4.1.3)，于 ．记． ． ，只要误差充分小就有如果对任给小正数 ， (4.1.8) 它表明求积公式(4.1.3)计算是稳定的，由此给出： ，只要就有在任给(4.1.8)，若定义３ 成立，则称求积公式(4.1.3)是稳定的． 中系数，则此求积公式是稳定的． 若求积公式(4.1.3)定理２ 都有，则有，对证明 对任给，若取 由定义３可知求积公式(4.1.3)