数值积分与数值微分课件.ppt

北京科技大学数理学院 卫宏儒 weihr168@ 数值积分与数值微分 引言 在实际问题中，往往会遇到被积函数f x 的原函数无法用初等函数来表示，或函数只能用表格表示，或有的虽然能用初等函数表示，但太复杂，所以这些情形都需要去建立定积分的近似计算公式。 在数值积分方面，最容易得到的是用f x 的代数插值函数p x 来代替它，即： 将积分 区间细分，在每小区间内用简单函数代替复杂函数，这是数值积分的基本思想。 对替代函数的要求： 1：精度要高。 2：计算量要小。 上述公式称为数值求积公式。其中 仅与 一、代数精度的定义及确定: 定义：若求积公式 ， 对一切不高于m次的多项式都准确成立，而对于m+1次多项式等号不成立，则称此公式的代数精度为m。 代数精度越高，公式越精确。 代数精度的求法： 从 依次验证求积公式是否成立，若第一个不成立的等式是 ，则该求积公式的代数精度就是 。 几个常用的Newton-Cotes公式 柯特斯系数 当n 2时，为抛物线公式 等距结点求积公式的余项问题 N-C求积公式的数值稳定性 例2：用Newton-Cotes公式计算 解：当n取不同值时，计算结果如下所示。 I准 0.座机电话号码 复合求积公式 从余项的讨论看到，积分区间越小，也可使求积公式的截断误差变小。因此，我们经常把积分区间分成若干小区间，在每个小区间上采用次数不高的插值公式，如梯形公式或抛物线公式，构造出相应的求积公式，然后再把它们加起来得到整个区间上的求积公式，这就是复合求积公式的基本思想。 复合求积公式克服了高次Newton-Cotes公式计算不稳定的问题，其运算简单且易于在计算机上实现。 常用的复合求积公式是复合梯形公式和复合抛物线公式。 复合梯形公式、复合抛物线公式 复合梯形公式 把区间[a,b] n等分，取节点xk a+kh,k 0,1,...n, h b-a /n,对每个小区间[xk,xk+1]用梯形求积公式，再累加起来得： （1） 复合抛物线Simpson公式 令n 2m,m为正整数，在每个小区间[x2k-2,x2k]上用抛物线求积公式，有： （2） 用复化梯形公式计算 令h 1/8 0.125，n 8 用复化抛物线计算 令h 1/8 0.125,m 4,n 8 Richardson外推法 假设有一个量F*，用一个步长为h的函数F1 h 去逼近，F*与h无关，既有h ? 0 时， F1 h ? F*，其与F1 h 的截断误差有估计式： R F* F*-F1 h a1h p1+a2hp2+····+a k p k+··· 1 p k p k-1 ··· p2 p1 0,a i i 1,2, ··· 都是与h无关的常数，也就是说，F1 h 逼近F*的阶是hp1 ，现在提出的问题是能否通过构造出一个新的序列，它逼近F*的阶要比hp1更高，如为hp2 。 将 1 中的h用qh来代替， q?0， 则有 F*-F1 qh a1 qh p1+a2 qh p2+···+a k qh pk+··· 现在用hp1乘 1 的两边后和上式相减,整理得 1-qp1 F*- F1 qh -qp1F1 h a2 qp2 - qp1 hp2+···+a k q p k- qp1 h p k+··· 因为 1-qp1 不等于零，用 1-qp1 除等式两边有 a2 2 , . . . ,a k 2 , . . . 都是与h无关的常数，令 F2 h 3 那么F2 h 逼近F*的误差由 2 知道为hp2. 依次做下去,计算公式为 F m+1 H ， m 1,2, . . . 4 用归纳法容易证明，由 4 得到的Fm h 逼近F*的误差为 上面的这种方法，称为Richardson外推法。 Romberg求积公式 注： 外推法不只用在积分的计算，也可用于计算其他量，只要能把该量看成f h 当 h? 0时的极限为f 0 ,而f 0 不易直接计算，要通过 hi 序列对应的函数值序列 f hi 来外推。 如果f h 有类似的展式，即 f h 则一切的算法和讨论可类似进行，更一般地设为 f h 其中 0 r1 r2 ··· rm，ri不必是整数，?I是不依赖于h的常数，?0 limf h ，可以类似地推出外推计算公式。 例：单位圆内接正n边形的面积为 S sin 设h 2/n,当n ? ?有h ?0可以类似地推出外推计算公式。 解:将S与h的关系写成 S h 。 当h ? 0时,S h 的极限是?，即单位圆的面积。将S h 作Taylor展开： ··· 即S h 有如上展式的形式，注意?是无理数，所以可按Romberg算法计算积分的过程来计算?到任意精度。 取n 6,12,24,即h 1/3,1/6,