数值分析-第二章--数值积分.ppt

第二章 数值积分 （Numerical Integration） 2.1 机械求积 2.2 牛顿-柯特斯公式 2.3 龙贝格算法 2.4 高斯公式 2.5 数值微分 内容提纲 数值积分的必要性 求积公式及其代数精度 插值型求积公式  1. 数值积分的基本思想(必要性) 本章主要讨论如下形式的一元函数积分 在微积分里，按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)没有解析表达式，只有数表形式: 2. 求积公式及其代数精度 求积公式的概念： 积分值　　　　　　　　　 在几何上可解释为由x=a, x=b, y=0和 y=f (x)所围成的 曲边梯形的面积. 积分计算之所以有困难，就是因为 这个曲边梯形有一条边 y=f (x)是曲的. 　　 依据积分中值定理,对于连续函数f(x) ,在[a,b]内存在一点ξ,使得 如果简单地选取区间[a, b]的一个端点或区间 中点的高度作为平均高度, 这样建立的求积公式 分别是: 左矩形公式: I( f )≈(b-a) f(a) 右矩形公式: I( f )≈(b-a) f(b) 中矩形公式: 此外, 众所周知的梯形公式: I( f )≈(b-a)[ f(a)+f(b)]/2 和 Simpson公式: I( f )≈(b-a)[ f(a)+4f((a+b)/2)+f(b)]/6 则分别可以看作用 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(ξ), 这类求积方法称为机械求积: 或写成: 我们的任务就是确定积分系数Ai ，使 I( f )≈In( f )。具体地，就是用插值多项式近似代替被积函数 f(x) 来确定Ai 。 在微积分中，定积分是Riemann和的极限，即 数值积分就是取定积分极限中的有限项的和，即 其中，xi 称为积分节点, Ai 称为积分系数。 求积公式的代数精度 定义1 称求积公式(2)具有m次代数精度,如果它满 足如下两个条件: (i) 对所有次数≤ m次的多项式 , 有 (ii) 存在m+1次多项式 , 使得 3. 插值型求积公式 在积分区间[a,b] 上取n+1个节点xi，i=0,1,2,…,n,作f(x)的n次代数插值多项式（拉格朗日插值公式）: 则有 为插值余项. 于是有 取 称(4)式为插值型求积公式, 其中求积系数Ak由(5) 式确定. 注 对于一般的函数 f(x)，En( f )≠0, 但若 f(x) 为次数小于 n 的多项式时，因为 f (n+1) (x) = 0, ∴ En( f ) = 0。因此 n 次插值多项式形式的数值积分公式至少有 n 阶代数精度。 例1 建立[0,2]上节点为 x0=0, x1=0.5, x2=2的数值积分公式. 解 n = 2, I2( f ) = α0 f(x0) + α1 f(x1) + α2 f(x2). 因为 记 定理2 当阶数n为