数值剖析数值积分.pptVIP

第4章 数值积分;§1 引言 ; （1）被积函数，诸如 等等，找不到用 初等函数表示的原函数； ;就是说，底为 而高为 的矩形面积恰等于所求 ; 问题在于点ξ的具体位置一般是不知道的，因而难以 ;图4-2; 一般地，可以在区间 上适当选取某些节点 ，; 这类数值积分方法通常称为机械求积，其特点是将积 分求值问题归结为函数值的计算，这就避开了牛顿-莱布尼 茨公式需要寻求原函数的困难. ; 2. 代数精度的概念 ; 欲使求积公式 具有 次代数精度，则只要令它; 如果事先选定求积节点 ，譬如，以区间 的等 距分点作为节点，这时取 ，求解方程组即可确 定求积系数 ，而使求积公式 至少具有 次代数精度. ;; 3. 插值型的求积公式;称为是插值型的，式中求积系数 通过插值基函数 积分得出 ; 当 是次数不超过 的??项式时，插值多项式就是; 定理1; 4 . 求积公式的收敛性与稳定性;如果对任给小正数; 定理2;由定义3知，求积公式是稳定的. ;§2 牛顿-柯特斯公式 ;; 当 时，; 的牛顿-柯特斯公式称为柯特斯公式，;25; 从柯特斯系数表看到 时，柯特斯系数 出现 负值，;它表明初始数据误差将会引起计算结果误差增大，即计算 不稳定，故 的牛顿-柯特斯公式是不用的. ; 2. 偶阶求积公式的代数精度 ;另一方面，直接求积得 ; 证明;因为被积函数; 3 . 几种低阶求积公式的余项 ; ;§3 复化求积公式 ;; 其余项;; 此外， 的求积系数为正，由定理2知复化梯形公式是 稳定的. ;;; 2. 复化辛普森求积公式 ;称为复化辛普森求积公式. ; 实际上，只要; 例2;而如果将[0,1] 分为4等分，应用复化辛普森法有 S4=0.9460832 ;于是 ;对复化辛普森公式， ;;§4 Richardson外推法 ;;; §5 龙贝格求积公式 ;;; 若 (预先给定的精度)，则终止计算，; 可以证明，如果 充分光滑，那么T表每一列的元素及对角线元素均收敛到所求的积分值 ，即 ; 例4; 从表中看到用龙贝格算到 的精度与辛普森求积 精度相当. 这里 的精确值为