第六章数值积分.pptVIP

梯形公式 辛蒲生公式 柯特斯公式 牛顿-柯特斯公式： 例1：分别利用梯形公式、 辛蒲生公式、 柯特斯公式 计算积分 的近似值。 解： * * * * * 第六章 数值积分 本章主要内容： 1、牛顿-柯特斯求积公式 2、复化求积公式 3、龙贝格求积公式  数值积分的必要性 本章主要讨论如下形式的一元函数积分 在微积分里，按Newton-Leibniz公式求定积分 要求被积函数f(x) ? 有解析表达式； ? f(x)的原函数F(x)为初等函数． 实际问题 1. f(x)的原函数F(x)不能用初等函数表示 考虑一个实际问题: 建筑上用的一种铝制波纹瓦是用一种机器将一块平整的铝板压制而成的. 例如函数: 假若要求波纹瓦长4英尺,每个波纹的高度(从中 心线)为1英寸,且每个波纹以近似2π英寸为一个周期. 求制做一块波纹瓦所需铝板的长度L. 这个问题就是要求由函数 f(x)=sin x 给定的曲线, 从x=0到x=48英寸间的弧长L. 由微积分学我们知道,所求的弧长可表示为: 上述积分称为第二类椭圆积分,它不能用普通方法来计算. 2. 有些被积函数其原函数虽然可以用初等函数表示成有限形式,但表达式相当复杂,计算极不方便. 并不复杂, 但它的原函数却十分复杂: 例如函数  3. f(x)没有解析表达式，只有数表形式: x 1 2 3 4 5 f(x) 4 4.5 6 8 8.5 这些都说明,通过原函数来计算积分有它的局限性,因而,研究关于积分的数值方法具有很重要的实际意义. 6.1 求积公式及其代数精度 积分值　　　　　　　　　 在几何上可解释为由 x=a, x=b, y=0 和 y=f(x) 所围 成的曲边梯形的面积. 积分计算之所以有困难，就 是因为这个曲边梯形有一条边 y=f(x) 是曲的. 一.求积公式的概念 依据积分中值定理, 对于连续函数f(x) ,在[a,b]内 存在一点 ξ, 使得 称f(ξ)为区间[a, b]的平均高度. 问题在于点 ξ 的具体位置一般是不知道的.这样,只要对平均高度f(ξ)提供一种算法, 相应地便获得一种数值求积方法. 如果简单地选取区间[a, b]的一个端点或区间中点的高度作为平均高度, 这样建立的求积公式分别是: 左矩形求积公式； 右矩形求积公式； 中矩形求积公式； 此外,众所周知的梯形公式: 和 Simpson公式: 则分别可以看作用 a, b, c=(a+b)/2, 三点高度的加权 平均值 [f(a)+f(b)]/2 和 [f(a)+4f(c)+f(b)]/6 作为平均高度f(ξ)的近似值. 更一般地, 取区间[a,b]内n+1个点 {xi}(i=0,1,2,…n) 处的高度{f(xi)} (i=0,1,…,n), 通过加权平均的方法近 似地得出平均高度 f(ξ), 称为机械求积公式. 即取 得到近似求积公式的一般形式 (*) ? 一般把积分公式表示成如下形式： 其中 求积节点 求积系数 仅与求积节点有关 求积公式的截断误差或余项： 求积系数的特征： 代数精度的判别方法 ? 二. 求积公式的代数精度 如果求积公式 求积公式 具有m次 代数精度的充要条件是 为 时求积公式精确成立，而 为 时求积公式 不能成为等式。 定义： 定理： 对一切不高于m 次的多项式都恒精确成立，而对于某 个m+1次多项式不能精确成立，则称该求积公式具有 m 次代数精度。 ?求积公式的收敛性和稳定性 若 则称求积公式（*）是收敛的。 设 有舍入误差 ，实际计算的求积公式为： 两者的误差为 其中 求积系数全为正时，公式是稳定的 三、插值型求积公式 思想 用被积函数 在区间 上的插值多项式近似代替 f (x) 计算. 作n次拉格朗日插值多项式: 设已知函数 在节点 上的函数值 其中 插值型求积公式： 余项 形如 的求积公式至少有 证明： 充分性 设它是插值型求积公式 当 时， 即它对所有不超过 n 次的多项式精确成立，故至少有n 次代数精度。 定理6.1 n 次代数精度的充要条件是它是插值型求积公式。 则对所有不超过n次的多项式求积公式精确成立 取 因此求积公式 是插值型的。 必要性 设求积公式至少有