数值分析--第4章 数值积分与数值微分.ppt

第4章 数值积分与数值微分;4.1 数值积分概论 ; （1）被积函数，诸如 等，找不到用 初等函数表示的原函数，或者即使能求得原函数但原函数的 表达式非常复杂，计算困难； ;就是说，底为 而高为 的矩形面积恰等于所求 ; 问题在于点ξ的具体位置一般是不知道的，因而难以 ;图4-2; 一般地，可以在区间 上适当选取某些节点 ，; 这类数值积分方法通常称为机械求积，其特点是将积 分求值问题归结为函数值的计算，这就避开了牛顿-莱布尼 兹公式需要寻求原函数的困难. ; 4.1.2 代数精度的概念 ; 欲使求积公式(1.3)具有 次代数精度，则只要令它; 如果事先选定求积节点 ，譬如，以区间 的等 距分点作为节点，这时取 ，求解方程组(1.4)即可确 定求积系数 ，而使求积公式(1.3)至少具有 次代数精度. ; 例如 时，取 ，求积公式为; 当 时（1.4）式的第3个式子不成立，因为; 如 ，此时求积公式为; 再令 ，代入（1.4）的第3式有; 例1 给定形如 的 求积公式，试确定系数 ，使公式具有尽可能高的 代数精度.; 当 时，得; 4.1.3 插值型的求积公式;（1.6）; 当 是次数不超过 的多项式时，插值多项式就是;注意到 上式右端实际上即等于 ，因而; 若求积公式(1.3)的代数精度为 ，则有求积公式余项的表达式(1.7)可以证明余项形如;（1.9）; 对中矩形公式(1.2)，其代数精度为1，可以证明; 例2 求例1中求积公式; 4.1.5 求积公式的收敛性与稳定性;如果对任给小正数; 定理2 若求积公式(1.3)中系数; 定理2表明，只要求积系数 ，就能保证计算的稳 定性. ;4.2 牛顿-柯特斯公式 ;（2.2）; 当 时，按(2.2)式，; 的牛顿-柯特斯公式称为柯特斯公式，其形式是 ; 从柯特斯系数表看到 时，柯特斯系数 出现 负值，;它表明初始数据误差将会引起计算结果误差增大，即计算 不稳定，故 的牛顿-柯特斯公式是不用的. ; 4.2.2 偶阶求积公式的代数精度 ;这时有 ，即辛普森公式对次数不超过三次的多项式 均能准确成立，而它对 通常是不准确的，因此， 辛普森公式实际上具有三次代数精度.; 证明 我们只要验证，当 为偶数时，牛顿-柯特斯公式对  的余项为零. ;因为被积函数; 4.2.3 辛普森公式的余项 ;从而可得辛普森公式的余项为;4.3 复合求积公式 ;（3.1）; 由(1.10) ，余项;（3.3）; 此外， 的求积系数为正，由定理2知复合梯形公式是 稳定的. ; 4.3.2 复合辛普森求积公式 ; 由(2.5)，其余项; 例3 对于函数 ，给出 的函数表(见表4-2)，试用复合梯形公式(3.2)及复合辛普森公式(3.5)计算积分 ; 同积分的准确值 比较，复合梯形法的结果 只有两位有效数字，而复合辛普森法的结果 却有6位有效数字.;于是 ; 例4 计算积分 ，若用复合梯形公，问区间 应分多少等份才能使误差不超过 ，若改用复合辛普森公式，要达到同样的精度，区间 应分多少等份？;因此有 ，可取 ，即将区间 213等份，即可使误差不超过; 从这个例子可以看出，为达到同样的精度，复合辛普森 公式只需计算9个函数值，而复合梯形公式则需214个函数值， 工作量相差近24倍.;4.4 龙贝格求积公式 ;用复合梯形公式求得该子区间上的积分值为 ;从而利用式(3.2)可导出下列递推公式 ; 例5;利用递推公式(4.1)，有 ; 它表明用复化梯形公式计算积分 要达到7位有效数 字的精度需要二分区间10次，即要有分点1025个，计算量 很大. ; 4.4.2 外推技巧 ; 定理4 设