数值计算方法第2版 第6章 数值积分和数值微分.pptVIP

第6章 数值积分和数值微分  ; 本章的问题： 计算定积分∫abf(x)dx的近似值。 必要性: 如果f(x)的原函数是F(x)，则 ;;;6.1.2 代数精度;;;;代数精度与节点数的关系;;6.1.3 插值求积公式;;;6.1.4 构造插值求积公式的步骤;;;用待定系数法构造插值求积公式;;;;;;;;;;;;6.2 牛顿-柯特斯求积公式6.2.1 公式的导出;2 柯特斯系数的求取;;;n ;柯特斯系数的性质 ;梯形公式 ;;抛物线(辛卜生)公式;;辛卜生公式的几何意义 ;柯特斯公式;;;6.2.2 牛顿－柯特斯公式的代数精度; 定理 当n是偶数时，牛顿－柯特斯求积公式具有n+1次代数精确度。 梯形公式， n=1( 2个节点）,有1次代数精度，应用梯形公式不是因为其代数精度高，而是因为其简单。 辛卜生（抛物线）公式，n=2 ( 3个节点）,有3次代数精度， 柯氏公式，n=4 ( 5个节点）,有5次代数精度。因为其代数精度高，所以常采用。 当n=3 (4个节点）,因为n=3不是偶数，只有3次代数精度，所以该公式不采用。 ;;证 ; 对梯形公式余项的说明 负号 f(x)的2阶导数，有1次代数精度。 3 和区间的3次方成正比。 ;例 证梯形公式的代数精度为1。 ; 定理 （辛卜生公式的余项）设f(x)在[a,b]上具有连续的四阶导数，则辛卜生公式的余项 ; 对辛卜生公式余项的说明 负号 f(x)的4阶导数，有3次代数精度。 3 和区间的5次方成正比。 ;例 证明辛卜生公式的代数精度为3。 ; 对科特斯公式余项的说明 负号 f(x)的6阶导数，有5次代数精度。 3 和区间的7次方成正比。 ;;; 梯形公式、辛卜生公式和柯特斯公式在区间不大时，用来计算定积分是简单实用的。但当区间比较大时，由余项可以看出精度差（梯形公式、辛卜生公式和柯特斯公式的余项分别和区间长度的3，5，7次方成正比），为减小因区间过大造成的误差过大，将积分区间等分成 n 等份，对每等份（每个小区间）分别用低阶的牛顿－柯特斯公式（如梯形公式、辛卜生公式或柯特斯公式）求积，然后将其结果加起来，得到积分的近似值。;; 将积分区间等分成n个小区间，在每个小区间上分别用梯形求积公式求积，然后再将其结果加起来。 设f(x) 在[a,b]上有连??的二阶导数，n是正整数. 将[a,b]等分成n个小区间 ;在[a,b] 上的误差;例 用n=6的复化梯形公式计算积分 ;2 复化辛卜生（抛物线）公式及其误差;;使绝对误差小于10 – 6。 ;3 复化柯特斯公式及其误差;;;;;;;;;6.3 变步长求积和龙贝格算法; 将积分区间等分成n个子区间，则有n+1个分点 ; 比较;;;;;;若将区间再分半，为 则有;当 时，T2n即为所求的近似值。 是T2n 的修正项，它与T2n 之和比T2n 、Tn更接近与真值，即它是一种补偿。 取;下面说明 将Tn ,T2n的表达式;;; 2 辛卜生法的加速 当把区间分成n等份，用复化辛卜生公式计算积分的近似值为 ，截断误差为;设 连续且变化不大时，有 ,可得近似式; 3 龙贝格公式 （柯特斯法的加速） 当把区间分成 n 等份，用复化柯特斯公式计算积分的近似值为 ，截断误差为;设 连续且变化不大时，有 ,可得近似式;龙贝格积分法可以按下面表的顺序进行：; 当对角线上最后两个相邻项满足;例 用龙贝格积分法计算积分 ;;;;;6.4.2 高斯-勒让德（Gauss-Legendre）求积公式; 从定理可以看出，当 给定，节点 就确定了。 ;;;;;;本题的精确解 ,求积公式具有4为有效数字。 ;6.5 数值微分;; 两点公式;