76定积分的数值计算.ppt

定积分的数值计算 数值积分 * * 定积分的计算——Newton-Leibniz公式 关键：找到 f (x)的一个原函数 原函数能用初等函数表示的仅占很小一部分，很多函数原函数虽存在，却不能用初等函数表出，如 在实际问题中，通常不能得到函数的解析表达式，仅能得到函数在若干离散点上的函数值。 寻找可积函数 p，满足 p的原函数能用初等函数表示 最简单的可积函数——多项式 由中值定理 y a b o x 0次多项式 左矩形公式 右矩形公式 可行吗？ a b o x y 一次多项式 梯形公式 a b o x y 二次多项式 抛物线公式(Simpson公式) 一般的，将区间[a,b] n 等分， 以 xi 为节点作 n 次 Lagrange 插值多项式 Newton- Cotes公式 Cotes系数 Cotes系数与区间[a,b]及f (x)均无关！ (1)对称性 (2)规范性 Cotes公式 Cotes系数的特点 梯形公式 抛物线 公式 3.1421 Cotes公式 3.1333 Simpson公式 3.0000 梯形公式 计算结果 方法 衡量数值积分方法精度的一个指标——代数精度 定义：若一个数值求积公式在被积函数是任意不高于n 次的多项式时都精确成立，而存在着一个n+1次多项式使公式不能精确成立，则称该求积公式具有n次代数精度。 可以证明： n为奇数，n步Newton-Cotes公式的代数精度为n n为偶数，n步Newton-Cotes公式的代数精度为n+1 命题： n步Newton-Cotes公式的代数精度至少为n。 是否代数精度越高的数值积分方法越好？ 例 -5 -4 -3 -2 -1 0 1 2 3 4 5 -1.5 -1 -0.5 0 0.5 1 1.5 2 n=2 n=4 n=6 n=8 n=10 f(x)=1/(1+x2) 不同次数的Lagrange插值多项式的比较图 Runge现象 复化求积 解决办法： 将积分区间分成若干等份； 在每一个小区间上使用低步数的Newton-Cotes公式； 再将各个小区间上的积分近似值加起来。 ? 复化梯形公式： 在每个 上用梯形公式： 复化梯形公式示意图 ? 复化Simpson公式： 在每个 上用Simpson公式： 记 为区间 的中点， *