计算方法ch5.pptVIP

第五章 插值型数值微分与数值积分;5.1 插值型数值微分公式 ; 5.1.1 常用的数值微分公式 ;若给定三点上的函数值 则由 ;进一步由 可得计算公式 ; 5.1.2 数值微分公式的误差分析 ;三点公式的截断误差为 ; 例 5.1 为计算 在 x=2 处的一阶导数值，我们可选用中点公式 ;为估计二阶导数数值微分公式的误差，可设 f (x) 四阶连续可微，故得 ;5.2 插值型数值积分 ; 5.2.1 牛顿柯特斯公式 ;特别地 ;5.2.2 复合求积公式 ;复合求积的方法： ; 例 5.2 试利用表5-3的函数表，分别用复合梯形公式、复合Simpson公式和复合Cotes公式计算定积分 ;5.2.3 插值型求积公式的误差分析与步长减半算法 ;定理 5.1 设 In 为由 N-C 公式（5-10）计算生成，则当 n 为奇数时，In 的代数精度为 n ；当 n 为偶数时， In 的代数精度为n+1，且当 在区间 [a, b] 上连续时，我们有如下误差估计式 ;特别地， ;从而可得 ; 为便于估计误差，实际计算时常常采用步长逐次减半的算法，下面介绍其思想。 ;因此，可先用 计算出T1，并把步长减半算出T2 ，若 则T2 即为所求的近似值，否则再把步长减半，算出T4， 根据式 (5-18a) 进行事后误差估计 ， 如此递推计算，直到某个n 满足 为止 ，取 为所求的近似值，这就是梯形公式的步长逐次减半算法。;为减少计算量，需建立递推公式，现对复合梯形公式推导之。 ;因此可建立梯形公式的步长逐次减半递推公式： ; 例 5.3 试用梯形公式的步长逐次减半算法计算定积分 使误差小于 。 ; 5.2.4 龙贝格积分法 ;类似地， ; 例 5.4 试用龙贝格积分法求解例5.3的定积分 使误差小于