数值分析第4章1节.pptVIP

第4章 数值积分与数值微分 4.1.2 代数精度的概念 4.1.3 插值型的求积公式 4.1.4 求积公式的收敛性与稳定性 4.1 引言 4.1.1 数值求积的基本思想 一、问题 如何求积分 牛顿-莱布尼茨(Newton-Leibniz)公式： N-L公式失效的情形： 数学分析中的处理方法： （1）被积函数，诸如 等等，找不到用 初等函数表示的原函数； （2）当 是由测量或数值计算给出的一张数据表. 这时，牛顿-莱布尼茨公式也不能直接运用; 因此有必要研究积分的数值计算问题及数值积分问题. N-L公式失效的情形： （3）原函数很难求. 问题：点ξ的具体位置一般是不知道的，因而难以 准确算出 的值，怎么办？ 只要对平均高度 提供一种算法，相应地便可获得 一种数值求积方法. 由积分中值定理知，在积分区间 内存在一点ξ， 成立 二、构造数值积分公式的基本思想 （1）左矩形公式 （3）用区间中点 的“高度” 近似地取代平均 高度 ，则又可导出所谓中矩形公式 （2）右矩形公式 （4）用两端点“高度” 与 的算术平均作为平均高度 的近似值，这样导出的求积公式 是梯形公式. 一般地，可以在区间 上适当选取某些节点 ， 然后用 加权平均得到平均高度 的近似值，这样 权 仅仅与节点 的选取有关， 构造出的求积公式具有下列形式： 的具体形式. 而不依赖于被积函数 式中 称为求积节点； 称为求积系数,亦称伴随节点 的权. k A 将这种思想一般化： 用上面式子求积分近似值的特点：将积分求值问题转化为了计算函数值的问题，避开了求原函数.这类数值积分方法通常称为机械求积。 定义1 如果某个求积公式对于次数不超过 的多项式 均能准确地成立，但对于 次多项式就不准确成立， 则称该求积公式具有 次代数精度. 数值求积是近似方法，为保证精度，自然希望求积公 式对尽可能多的函数准确成立. 设给定一组节点 且已知函数 在这些节点上的值， 作插值函数 . 取 作为积分 的近似值， 这样构造出的求积公式 称为是插值型的，式中求积系数 通过插值基函数 积分得出 由插值余项定理(第2章的定理2)即知，其余项 式中ξ与变量 有关， 问题：上述插值型求积公式至少具有多少次代数精度？ 当 是次数不超过 的多项式时，插值多项式就是 函数本身， 余项 为零， 事实上， 这时求积公式对于插值基函数 应准确 成立，即有 至少具有 次代数精度. 所以这时插值型求积公式 反之， 如果一个求积公式 至少具有 次代数精度，则它必可利用插值多项式推导出来。 注意到 上式右端实际上等于 因而 成立. 这样，有下面定理. 定义2 其中 实际得的将是 ， 即 则称此求积公式是收敛的. 记 在求积公式 中，若 在求积公式 中，由于计算 可能产生误差 ， 如果对任给小正数 只要误差 充分小就有 （*） 则表明求积公式计算是稳定的， 由此给出下面定义. 定义3 对任给 若 只要 就有(*)成立，则称求积公式 是稳定的.