数值积分4.1、4.2、4.3数值积分微分.pptVIP

华长生制作 以上这些现象,Newton-Leibniz很难发挥作 用，只能建立积分的近似计算方法。本章 将研究求积分值的数值求积公式，即形如 的求积公式，其中 叫求积节点， 叫求积系数，它们 均是与f(x)无关的常数， 叫该求积公式的截断误差。 衡量一个求积公式好坏的标准之一是所谓 的代数精度。 定义 如果上述求积公式中，当f(x)为任意 次数不超过 m的多项式时误差 ，而 当 时 ，则称该求积公式 具有m次代数精度。 易证，求积公式具有m次代数精度的充分 必要条件是：f(x)分别取 时 ，而 时 凡是至少具有零次代数精度的求积公式一 定满足 即求积系数之和等于积分区间长度，这是 求积系数的基本特性。 如果 其中 则称求积公式是收敛的。 下面考虑稳定性 设 有舍入误差 ，则 与 之间的误差为 记 ，当求积系数全为正时 说明该误差不超过最大舍入误差 的常数 倍，即计算过程是稳定的。 则复合梯形公式的余项为 由于 即有 比较三种复合公式的的余项 为此介绍收敛阶的概念 定义1. 不难知道,复合梯形、Simpson、Cotes公式的收敛阶分别为 2阶、4阶和6阶 通常情况下,定积分的结果只要满足所要求的精度即可 三、复合求积公式步长的自动选取 复合梯形公式的余项为 因此有 即 依此类推 步长自动选取的步骤: 依此类推 以上这种方法称为区间逐次分半求积法 有时也去掉 精度会更高 4.3 龙贝格求积算法 * 第四章 数值积分与数值微分 本章要点 (1) 等距节点下的:Newton-Cotes公式和Romberg公式 (2) 数值微分公式 4.1 数值积分概述 对于积分 但是在工程技术和科学研究中,常会见到以下现象: 例1. 试确定下面积分公式中的参数使其代数精确度尽量高. 解: 因此 所以该积分公式具有3次代数精确度 构造数值求积公式的方法很多,最常用的一种方法是用 插值多项式来构造数值求积公式,具体步骤如下: 不同的 插值方法 有不同的 基函数 4.2 牛顿-柯特斯公式 这就是插值型求积公式 一、Newton-Cotes数值求积公式 Newton-Cotes公式是指等距节点下使用Lagrange插值 多项式建立的数值求积公式 各节点为 其中 而 因此对于定积分 有 令 即有 n阶Newton-Cotes求积公式 Newton-Cotes公式的余项(误差) 注意是等距节点 所以Newton-Cotes公式化为 思考 使用n次Lagrange插值多项式的Newton-Cotes 公式至少具有n次代数精度,并且n为偶数时至 少具有n+1次代数精度,试以n=1,2,4为例说明 该结果 二、低阶Newton-Cotes公式及其余项 在Newton-Cotes公式中,n=1,2,4时的公式是最常用也 最重要三个公式,称为低阶公式 1.梯形(trapezia)公式及其余项 Cotes系数为 求积公式为 上式称为梯形求积公式,也称两点公式，记为 梯形公式的余项为 第二积分 中值定理 梯形(trapezia)公式具有1次代数精度 故 2.Simpson公式及其余项 Cotes系数为 求积公式为 上式称为Simpson求积公式，也称三点公式或抛物线公式 记为 Simpson公式的余项为 Simpson公式具有3次代数精度 3.Cotes公式及其余项 Cotes系数为 求积公式为 上式称为Cotes求积公式，也称五点公式 记为 Cotes公式的余项为 Cotes公式具有5次代数精度 表 4.1 989/28350 5888/28350 751/17280 -928/28350 3577/17280 41/ 840 10496/28350 1323/ 17280 216/ 840 19/288 -4540/ 28350 2989/ 17280 27/840 75/288 7/90 10496/28350 2989/ 17280 272/ 840 50/288 32/90 1/8 -928/28350 1323/17