第6章 数值积分解析.ppt

第六章 数值积分 引言 数值积分是数值分析的重要内容，也是函数插值的直接应用。在工程计算中，由于许多函数的不定积分无法用简单函数解析地表达出来，甚至被积函数本身都无法详尽地描述，而只能以简单的表格形式给出一些离散点上的函数值,或者定义为某个无法用显示形式表达的微分方程的解。在上述这些情况下，我们只能用数值方法求函数的定积分。 例如，在土地丈量中会遇到各种各样不规则地块，由于我们无法知道其边缘曲线满足何种函数关系，因此计算它的面积存在一定的困难。 又如，求定积分 数不能用初等函数解析地表示出来，因此这些定积分都不能用Newton-Leibniz 公式求出，只能用数值方法求 其近似值 。 § 6.1 求积公式及其代数精度 代数精度的概念 截断误差公式往往与被积函数及其导数信息有关，因而无法直接用来衡量求积公式本身的好坏。通常我们用求积公式的代数精度来说明求积公式的好坏。 例1 (2) 求积公式的构造方法 数值积分的基本问题 针对某类函数，选择合适的求积结点和求积系数，使得求积公式(6.1.1) 具有尽可能小的截断误差或尽可能高的代数精度。 寻找被积函数 f (x)在区间[a ,b]上的一个容易求积的近似函数(如某个多项式)p(x)，然后用p(x)在区间[a ,b]上的积分值去近似f (x)在区间[a ,b]上的积分值。最常用的近似函数p(x) 是f (x)在区间[a ,b]上插值多项式函数。 § 6.2 插值型求积公式及其代数精度 方法： 截断误差估计 例2 插值型求积公式的代数精度 §6.3 Newton-Cotes求积公式 截断误差为 Cotes系数 表6-1 Newton-Cotes求积公式的代数精度 梯形公式： Simpson公式： Simpson3/8公式： Cotes公式 例3 §6.4 Newton-Cotes求积公式的收敛性数值稳定性