数值分析第六章.pptVIP

复化Simpson公式的几何意义 小抛物面积之和近似 * * 将每个小区间[xi, xi+1] 做4等分，分点分别记为 ， 在区间[xi, xi+1]采用Cotes公式: 6.3.3 复化Cotes公式 复化Cotes公式的截断误差 复化Cotes公式 * * 例6-4 将区间做4等分，然后用复化梯形公式及复化Simpson 公式计算 解 复化梯形 将积分区间[0,1]划分为4等分，则 2位有效数字，6位有效数字, 精度差别很大. 复化Simpson 1 4 , h = * * （1）使用复化梯形公式、Simpson公式，首先要确定步长h ； （2）而步长要根据余项确定，这就涉及到高阶导数的估计； （3）高阶导数的估计一般比较困难，且估计值往往偏大； 注意事项： 步长h如何选取？ * * 6.4 龙贝格求积公式/* Romberg Integration */ 6.4.1 变步长求积公式 变步长复化求积公式的思想 将积分区间逐次分半，建立递推公式计算，直到满足精度要求. 变步长复化梯形公式 逐次分半算法 直到满足精度要求为止。 * * 设将区间[a,b] 分为n等分，共有n+1个分点， 如果将求积区间再二分一次，则分点增至2n+1个， 每个 区间[xi , xi+1] 经过二分增加了一个分点 用复化梯形公式求该子区间上的积分值 * * 变步长复化梯形公式实际计算时的递推公式 直到 为止， 作为积分的近似值. * * 解 例6-5 利用变步长复化梯形公式计算 先对整个区间 使用梯形公式，即h=1 将区间二等分， h=1/2 ，求出中点的函数值 进一步二分求积区间，h=1/4 ，并计算新分点上的函数值 的近似值，精确到0.000001. 记 这样不断二分下去，计算结果见表6-1. 它表明用复化梯形公式计算积分I 要达到7位有效数字的精度需要二分区间10次，即要有分点1025个，计算量很大. * * 6.4.2 Romberg求积公式 由复化梯形公式的余项知 变化不大时 由此得到近似关系式 收敛慢！ * * 变步长复化辛普森公式 利用复化梯形公式前后两次积分近似值Tn 和T2n ，作出的 线性组合，即为复化Simpson公式，具有更高的精度。 * * 一般有： Romberg 序列 ? Romberg 算法： ? ? ? ? ? ? … … … … … … ? T1 = ) 0 ( 0 T ? T8 = ) 3 ( 0 T ? T4 = ) 2 ( 0 T ? T2 = ) 1 ( 0 T ? S1 = ) 0 ( 1 T ? R1 = ) 0 ( 3 T ? S2 = ) 1 ( 1 T ? C1 = ) 0 ( 2 T ? C2 = ) 1 ( 2 T ? S4 = ) 2 ( 1 T * * 例6-6 用龙贝格算法计算积分 解 k 0 1 2 3 4 5 3 3.1 3.1311764 3.1409416 3.1333333 3.1415687 3.1415925 3.1421177 3.1415942 3.1415858 3.1409416 3.1415927 3.1415927 3.1415927 3.1414299 3.1415927 3.1415927 3.1415927 * * 前面介绍的 n+1个节点的 Newton -Cotes求积公式，其特征是节点是等距的。这种特点使得求积公式便于构造，复合求积公式易于形成。但同时也限制了公式的精度。 n是偶数时，代数精度为n+1， n是奇数时，代数精度为n 。 我们知道 n+1个节点的插值型求积公式的代数精确度不低于n 。设想：能不能在区间[a,b]上适当选择n+1个节点 x0, x1, x2, …, xn , 使插值求积公式的代数精度高于n？ 答案是肯定的，适当选择节点，可使公式的精度最高达到2n+1，这就是本节所要介绍的高斯求积公式。 * * 数值计算方法NumericalCalculation Method 武汉学院信息系 授课人：江成顺 6.5 Gauss求积公式/* Gaussian Quadrature */ 6.5.1 Gauss公式的定义 构造具有2n+1次代数精度的求积公式 将节点 x0 … xn 以及系数 A0 … An 都作为待定系数。令 f (x)