第4章数值积分课题.ppt

第四章 数值积分;4.0 引言 函数f(x)在区间[a,b]上连续且其原函数为F(x),则有Newton-Leibnitz公式 ;?(1) 被积函数f(x)不一定能找到初等函数形式的原函数F(x)，例如： Newton-Leibnitz公式就无能为力了;(3) 被积函数f(x)没有具体的解析表达式, 其函数 关系由表格或图形表示。 对于这些情况, 要计算积分的准确值都是十分困难的。由此可见, 通过原函数和N-L公式不能或很难解决所有的积分问题, 这时就需要用数值解法来建立积分的近似计算方法。 将积分区间细分,在每一个小区间内用简单函数代替复杂函数进行积分，这就是数值积分的思想，用代数插值多项式去代替被积函数f(x)进行积分是本章讨论数值积分的主要内容。 ;4.1 数值积分概述 4.1.1 数值积分的基本思想 积分值 在几何上可以解释为由x=a,x=b,y=0以及y=f(x)这四条边所围成的曲边梯形面积。如图4-1所示，而这个面积之所以难于计算是因为它有一条曲边y=f(x) ; 建立数值积分公式的途径比较多, 其中最常用的有两种： 方法（1）由积分中值定理可知，对于连续函数f(x)，在积分区间[a,b]内存在一点ξ，使得 即所求的曲边梯形的面积恰好等于底为(b-a)，高为 的矩形面积。但是点ξ的具体位置一般是未知的, 因而 的值也是未知的, 称 为f(x) 在区间[a,b]上的平均高度。那么只要对平均高度 提供一种算法，相应地就获得一种数值求积方法;三个求积分公式;y=f(x); 一般情况下，只要对平均高度 提供一种算法，相应地就获得一种数值求积方法，具体方法如下：;方法(2) 先用某个简单函数 近似f(x),用 的积分, 近似函数f(x)的积分，即 ;设已知f(x)在节点 有函数值 ,作n次拉格朗日插值多项式 ;其中 ;设插值求积公式的余项为 ,由插值余项定理得 ;定义 （代数精度） 设求积公式（4.1）;这是关于 的线性方程组，其系数矩阵;定理4.1 n+1个节点的求积公式 为插值型求积公式的充要条件是公式至少具有n次代数精度。 ;充分性 ： 若求积公式至少具有n次代数精度, 则对n次多项式;（1）;（2）;例4.1 设积分区间[a, b]为[0, 2]，取 时, 分别用梯形和辛卜生公式 ; f(x) 1 x x2 x3 x4 ex 准确值 2 2 2.67 4 6.40 6.389 梯形公式计算值 2 2 4 8 16 8.389 辛卜生计算值 2 2 2.67 4 6.67 6.421; f(x) 1 x x2 x3 x4 ex 准确值 2 2 2.67 4 6.40 6.389 梯形公式计算值 2 2 4 8 16 8.389 辛卜生计算值 2 2 2.67 4 6.67 6.421;取f(x)=1时， ;例4.2 试确定一个至少具有2次代数精度的公式 ;例4.3 试确定求积系数A,B,C 使 具有最高的代数精度 解:分别取f(x)=1,x,x2 使求积公式准确成立,即 得如下方程组。; 代入公式两端，左端和右端都等于(b4-a4)/4,公式两端严格相等，再将f(x)=x4代入公式两端，两端不相等，所以该求积公式具有3次代数精度。;可以验证, 对于f(x)=1, x 时公式两端相等, 再将f(x)=x2代入公式;例4.5 给定求积公式如下： ;则以这