计算方法课件：第6次课new 计算方法数值积分.ppt

* 计算方法3 数值积分 航空科学与工程学院 * 3 数值积分 问题的提出 各种典型的问题 数值积分 * 3.0 问题的提出 回顾 数学分析中关于积分的定义。 数值积分 * 积分函数 假设f(x)为定义在[a，b]上的可积函数，计算积分 如果F(x)为原函数，则有Newton-Leibniz公式 数值积分 * Q: (1) 有些函数的原函数不能用初等函数表现为有限的形式； (2) 原函数的形式复杂； (3) 原函数没有具体的表达式，只有离散点。——定积分的数值解法(效率+精度)。 数值积分 * 3.1 数值积分公式 一般形式 设(a，b)为有限或无限区间， 其中w(x)为权函数，满足： (1) 在(a，b)中，w(x)≥0，并且只有有限个零点； (2) 存在。 数值积分 * 用被积函数 f(x) 的若干节点xi (a≤x0 < x1< …< xn≤b)处的函数值 f(xi) 的线性组合 作为I(f )的近似值——数值积分公式，求积公式。 求积节点：xi(i=0，1，…，n) 求积系数：Ai(i=0，1，…，n)与f(x)无关，仅与求积节点有关； 余项或离散误差：En(f)=I(f)-In(f)。 数值积分 * 3.2 插值求积公式 一般思想 取简单的、便于积分且又逼近于被积函数f(x)的函数φ(x)代替f(x)来构造求积公式； 典型——插值多项式 数值积分 * Lagrange多项式： 数值积分 从而插值求积公式： 余项： * Newton-Cotes插值公式 设[a，b]为有限区间，取w(x)=1，h=(b-a)/n，等距节点xi=a+i·h(i=0，1，…，n)。 记x=a+t·h(0≤t ≤n)，则： 数值积分 * 而dx =hdt， ——n阶Newton-Cotes插值公式。 数值积分 * 通常取 Ci(i=0，1，…，n)为Cotes系数。 从而: 数值积分 * Cotes系数 Ci的特点 数值积分 * 特例 梯形公式 n=1，x0=a，x1=b，h=b-a 示意图 数值积分 * Simpson公式(抛物线插值公式) n=2，x0=a， x1=(b+a)/2， x2=b，h=(b-a)/2 数值积分 * Cotes公式 ——四阶(n=4)Newton-Cotes公式 适用性：高阶公式稳定性差 数值积分 * 例题 利用梯形公式和Simpson公式求积分 解：I1(f)=0.75; I2(f)=0.69444; I(f)=0有效数字？ 数值积分 * 代数精度 对 　　　 如果该求积公式对于一切次数≤m的多项式是准确的，但对于 m+1次多项式不准确，则称其具有m次代数精度。 描述形式? 数值积分 * 　 定理：求积公式具有n次代数精度的充要条件是：它是插值型的。 证明： 　　 (1) ≤n次的多项式的插值函数是其本身，从而插值型具有至少n阶精度； (2) 求积函数具有n次代数精度，则对其插值基函数lk(x)，有： 数值积分 * 典型插值公式的代数精度：二阶、四阶； 由插值精度构造插值公式； 数值积分 * 误差 定理： 设n为偶数，且f(x)在[a，b]上有n+2阶连续导数，则Newton-Cotes型求积公式的离散误差： 设n为奇数，且f(x)在[a，b]上有n+1阶连续导数，则Newton-Cotes型求积公式的离散误差： 数值积分 * 特例： n=1时，梯形公式的离散误差为