第四章数值积分与数值微分(第一次).ppt

* * * * * * * * * * * * 估计截断误差为 解 用梯形公式计算: =2.1835 估计截断误差为 =0.6796 用Simpson公式计算： =2. 0263． 例3 试分别使用梯形公式和Simpson公式计算积分 的近似值，并估计截断误差. =198.43 =0.06890 4、 求例3中求积公式 的余项 解 由于此求积公式的代数精度为2，故余项表达式为 . 令 ，得 ，于是有 故得 本 节 主 要 内 容 1、复化求积公式 构造思想 公式余项 2、龙贝格算法 构造思想 上机计算 3、高斯求积公式 构造过程 §4、3 复化求积公式 高次插值有Runge 现象，故采用分段低次插值 ? 分段低次合成的 Newton-Cotes 复合求积公式。 一、复化梯形公式: 在每个 上用梯形公式： =Tn /*中值定理*/ 二、复化辛普森公式: 4 4 4 4 4 = Sn 注：为方便编程，可采用另一记法：令 n’ = 2n 为偶数， 这时 ，有 三、收敛速度与误差估计： 定义 　　 若一个积分公式的误差满足 且C ? 0，则称该公式是 p 阶收敛的。 ~ ~ ~ 例4：计算 解： 其中 = 3.138988494 其中 = 3.141592502 运算量基本相同 问题: 给定精度 ?，如何取 n ? 例如：要求 ，如何判断 n = ？ ? 上述例4中若要求 , 则 即：取 n = 409 通常采取将区间不断对分的方法，即取 n = 2k 上述例4中2k ? 409 ? k = 9 时，T512 = 3例4中： S4 = 3.141592502 注意到区间再次对分时 可用来判断迭代 是否停止。 (1) (2) (3) 事后误差估计 xi 0 1/8 1/4 3/8 1/2 5/8 3/4 7/8 1 f (xi) 1 0.9973978 … … … … … … … … … 0.8414709 作业 P135： 2(1),3,6 一、梯形法的递推化——逐次分半法 上一节介绍的复化求积方法对提高精度是行之有效的，但在使用求积公式之前必须给出合适的步长，步长取得太大精度难以保证，步长太小则会导致计算量的增加，而事先给出一个恰当的步长又往往是困难的． 实际计算中常常采用变步长的计算方案，即在步长逐次分半(即步长二分)的过程中，反复利用复化求积公式进行计算，直至所求得的积分值满足精度要求为止． 设将求积区间[a，b]分成n等分，则一共有n+1个分点，按梯形公式计算积分值Tn，需要提供n+1个函数值．如果将求积区间再二分一次，则分点增至2n+1个，我们来考察二分前后两个积分值之间的联系． §4、4 龙贝格求积公式 逐次分半计算方案的实现: 注意到每个子区间[xk，xk+1]经过二分只增加了一个分点 xk+1/2＝( xk+xk+1)/2，用复化梯形公式求得该子区间上的积分值为 这里 代表二分前后的步长.将每个子区间上的积分值相加得 二、龙贝格算法 根据复化梯形公式的余项表达式 可见，利用两种步长计算的结果能估计截断误差.若将该截断误差加到计算结果中， 就得出“改进的梯形求积公式”： 事后误差估计 例：计算 已知对于? = 10?6 须将区间对分 9 次，得到 T512 = 3由 来计算 I 效果是否好些？ = 3.141592502 = S4 改进梯形求积公式的右边实际是 这就是说用梯形法二分前后的两个积分值Tn与T2n的线性组合的结果得到复化辛普森法求积公式 类似的情况，用辛普森法二分前后的两个积分值Sn与S2n的线性组合的结果可得到复化柯特斯求积公式 重复同样的手续，用柯特斯法二分前后的两个积分值Cn与C2n的线性组合的结果可得到龙贝格