第五章 数值积分..pptVIP

* 八、数值导数 * * 复合求积公式的截断误差为 * * 3．外推法 * * * * * 对比例4，同样用9个点上的函数值计算的数值，这个结果已经相当好了。作为算法，还需要在计算中估计误差，为此，我们再讨论误差。 * 4．实时误差估计 其它的复合求积公式有类似的结论。 * 5．Romberg算法 * * * * * * 六、 高斯（Guass）求积公式 * * 1．高斯点 * * * 显然，高斯点就是n+1次正交多项式的零点。 * 2．高斯-勒让德求积公式 * * * * 3.区间的变化 * * 4．高斯-勒让德求积公式的余项 注：在实际运用时一般不用高阶高斯公式，两点、三点和五点高斯公式是最常用的。 * 七 、其它的类型 * * * * 第五章 数值积分和微分 * 一、一般求积公式及代数精度 * * * * * 二、 插值型求积公式 * * * * 一般插值型求积公式的构造方法如下： * * * 三 、牛顿—柯特斯（Newton-Cotes）公式 定义3 等距节点下的插值型求积公式称为牛顿—柯特斯公式。 * * * * 柯特斯系数具有下列性质： * * 注 * 几个牛顿-柯特斯公式： * * 注 * 四、余项 * * 五 、复合求积公式与龙贝格（Romberg）算法 * 1. 复合求积公式构造 复合梯形公式 * 复合Simpson公式 * * * 2．复合求积公式的截断误差 设函数在区间上具有一定的光滑度，用数值 积分方法计算 的最一般的方法是用在节点 的函数值的线性组合来近似地表示： （1） 其中是与无关的常数称为求积系数， （1）称为求积公式，称为求积节点。 注： 1. 如果去掉定义一中的第二个条件，则称公式（1）至少具有 m次代数精度。 2. 若求积公式（1）对都精确成立，那么（1）对 亦精确成立。 定义1 如果（1）对于任何次数不高于次的代数多项式都精 确成立,对于不能精确成立,则称(1)具有次代数精度。 实际上，如果且， 则令：，那么： 3. 验证（1）式具有次代数精度，只需验证（1）式对 精确成立，而对于不精确 成立即可。 因为，若（1）式精确成立， 则对有 。 设,称为区间的一个分割。 以为插值节点作插值多项式，用 近似是一个自然的想法。由于： （2） 其中。 定义2 对于给定节点，如果求积系数 由（2）式确定，则称此时的求积公式为插值型求积公式。 定理1 是插值型求积公式的充分必要 条件是它的代数精度。 (2)求积系数只与积分区间及节点有关，而与被积函数 无关。因此，可以一次求出系数，多次使用。 (3)个节点的插值型求积公式至少具有次代数精度。实际上，如果是次数不大于的多项式，那么，所以（2）式准确成立． (4)求积系数之和等于区间的长度。 因为，插值型求积公式（2）至少具有一次代数精度， 令，这是零次多项式，（2）精确成立，所以： 。 在积分区间上选取节点，； 解关于的线性方程组： 求出，。这样就得到求积公式： 。 例1 确定插值型求积公式中的系数与节点。 把区间分成等分，每个小区间的长度为， 等分点为，作变换 ， 取，则， ． 这里， 与被积函数以及积分区间均无关，称为柯特斯系数。 公式： 称为牛顿—柯特斯公式。 常用的柯特斯系数： 1 2 3 4 5 ；这是因为。 对称性：。 应用中必须考虑数值稳定性：设函数值计算产生误差为：，并记，则在牛顿—柯特斯公式计算中：，误差是： 如果，则 ， 这公式是数值稳定的。 在牛顿——柯特斯公式系列中，时，。 开始出现负系数，公式稳定性下降，一般不再使用。 (1)、梯形公式： (2)、Simpson公式: (3)、柯特斯公式 例2 利用梯形公式和Simpson公式计算积分。 解：，， 由梯形公式得 ； 由Simpson公式得 。 对比真值，可见更精确一些。 解：利用梯形公式：； 利用Simpson公式：。 例3中我们发现实际上，这说明Simpson公式公式具有三次代数精度！这不是偶然的，稍后，我们可以看到牛顿——柯特斯公式代数精度当为奇数是次，但当为偶数时却是次。 例3 利用梯形公式和Simpson公式计算积分。 定理2 设，则 ； 设，则 ； 设，则 。 我们只就梯形公式证明。 由于： 从余项不难看出，由于次多项式的阶导数为零，故：和 分别具有三次和五次代数精度，而梯形公式是一次代数精度的。 求积公式的误差与求积区间的长度有关，积分区间越小，误差一般也小，但与此同时，会增加插值多项式的次数使误差迅速增大。因此将积分区间分成若干小区间，在每个小区间上用低阶求积公式，然后将这些小区间