[2018年最新整理]ch4数值积分与数值微分.pptVIP

第四章 数值积分与数值微分4 . 1引言 很多函数的原函数不能用初等函数表示，或者难于求得。这就很难用N-L公式计算积分的精确值。因此，寻求计算积分即使只是计算积分的近似值的方法，就非常必要了。数值积分是计算积分近似值的一种方法。 1.数值求积的基本思想 数值求积的基本思想，用积分区间内被积函数值的加权平均值作为积分的近似值。 值的方法常称为机械求积。 2.代数精度的概念 3.插值型求积公式 4.求积公式的收敛性与稳定性 定义 则称该求积公式是收敛的(其中： ) 定义 有 ，则称该求积公式是 稳定的。 定理 代精，则此求积公式是稳定的。 4.2 牛顿-柯特斯公式(N-C公式)1. 柯特斯系数 将 分为n等分,步长 2．偶阶求积公式的代精 定理 证明 3．几种低阶求积公式的余项 类似 4．3 复化求积 （n=1即前述公式） 2.复化辛普生公式 3.复化柯特斯公式 4．4龙贝格方法1.龙贝格求积公式 （1）事后误差估计 （2）用误差补偿 （2）用误差补偿 （2）用误差补偿 2. 梯形法的逐次分半算法 龙贝格算法：用龙贝格求积公式加工梯形值序列，从而得到精度较高的积分值。 3．理查森外推加速法 以上不断消去 ，从而不断提高积分公式精度的方法称为理查森外推加速法。 4.5 高斯求积公式1．高斯点与高斯求积公式 定理 带权的求积公式,其代精最高不超过2n+1次。 定义 若取 使求积公式具有2n+1次代精，则称 为高斯点，相应的求积公式称为高斯求积公式。 2．高斯求积公式的特征 定理 插值型求积公式中求积节点 为高斯点的充要条件是以这些节点为根的多项式 与任意不超过n次的多 项式 带权 正交即 ， 。 求 的方法 (1)据代精需满足的关系建立线性方程组并解之； (2) 为插值基函数， 为权。 3．高斯求积公式的性质 定理 高斯求积公式是稳定的。 4．高斯-勒让德求积公式 取 ， 区间，以勒让德多项式的0点作高斯点，进而得到的求积公式称为高斯-勒让德求积公式。 5．高斯-切比雪夫求积公式 取 ，区间 ，以切比雪夫多项式的0点作高斯点，进而得到的求积公式称为高斯-切比雪夫求积公式。 4 . 6 数值微分1.插值型求导公式 n=1: n=2 : 2．理查森外推算法 3．其它数值微分法 用样条函数求导近似代替； 用数值积分得线性方程组反推出来。 取 ， 新分点 （有继承性的算法） 其中 与 无关， （消去 ） ， ， 消去 一般， 例 用龙贝格方法计算定积分 的近似值， 。 要求误差不超过 解 构造表格如下： 所以， 解毕 考虑一般的带权的求积公式 证明 ，（2n+2次） ， 证毕 例 求 使 为高斯求积公式。 解 令 代入公式，使之准确成立： 代入 ， 是高斯求积公式。 其中， 是高斯点，代精为3。 解毕 证明 证明 证毕 据此定理,高斯点是n+1次正交多项式的n+1个根。 高斯点一旦确定，高斯求积公式也就随之确定。 定理 高斯求积公式是收敛的 。 证明 证毕 （证略） ， 进而得到两点G-L公式，有3次代精。 进而得到三点G-L公式，有5次代精。 ， 例 构造三点G-L求积公式。 解1 设 其中 令 所以， 解2 设 其中 所以， 解毕 一般的 ，区间 作代换 将 变换为 。 ， 例 用四点G-L求积公式，计算 。 解 解毕 两点公式： （向前差商） 余项 （向后差商） 余项 选择 三点等距。 （中心差商） 三点公式： 河海大学理学院《数值分析》 数 值 分 析 Computational Method Chapter 4 数值积分与数值微分 定义 称 为数值求积公式，其中： 称为求积节点， 称为求积系数，又称为 的权，权 仅与 有关，而不依赖于 的具体形式， 称为求积公式的截断误差，又称为求积公式的 余项。 数值求积公式中，关健要确定求积节点 和求积系数 ，用数值求积公式求得积分近似 定义 若某求积公式对次数为 的任意