第5章 福大科学工程与计算-数值积分.pptVIP

令其对f(x)=1, x 精确成立, 有 解出 A1＝Ａ2＝1, 从而得到代数精度为3的两点高斯－勒让德求积公式 三点高斯－勒让德求积公式的形式是 其它的高斯－勒让德求积公式的形式根据常用高斯－勒让德求积公式的节点和系数表，自己可以写出来. 常用高斯－勒让德求积公式的节点和系数表 n 节点数 求积节点xk 求积系数Ａk 0 1 0 2 1 2 ±0.5773503 1 2 3 ± 0.7745966692 0.0000000000 0,5555555556 0.8888888889 3 4 ± 0.8611363116 ± 0.3399810436 0.3478548451 0.6521451549 4 5 ± 0.9061798459 ± 0.5384693101 0.0000000000 0.2369268851 0.4786286705 0.5688888889 公式(5.9)的余项由(5.8)得到 这里是 最高项系数为1的勒让德多项式，由第三章的一些性质得到 当n=1时，有 它比辛普森公式 余项(区间[a, b]) 还小，且比 辛普森公式少算一个函数值. 在区间[a, b]上的高斯－勒让德求积公式 作变量替换 然后再用[-1,1]上的高斯－勒让德求积公式． 得 这就是中矩形公式． 例如, 用一点高斯－勒让德公式有 用两点高斯－勒让德公式有 例6 利用两点Gauss-Legendre 求积公式计算 解:因为 为偶函数 Simpson公式及其余项 Cotes系数为 求积公式为 上式称为Simpson求积公式，也称三点公式或抛物线公式 记为 Simpson公式的余项为 Simpson公式具有3次代数精度 Cotes公式及其余项 Cotes系数为 求积公式为 上式称为Cotes求积公式，也称五点公式 记为 Cotes公式的余项为 Cotes公式具有5次代数精度 使用n次Lagrange插值多项式的Newton-Cotes 公式至少具有n次代数精度,并且n为偶数时至 少具有n+1次代数精度. 5.2.2 复化求积公式 直接使用Newton-Cotes公式的余项将会较大 公式的舍入误差又很难得到控制 为了提高公式的精度,又使算法简单易行,往往使用复化方法 然后在每个小区间上使用低阶Newton-Cotes公式 最后将每个小区间上的积分的近似值相加 ? 复化 Simpson 公式： 4 4 4 4 4 = Sn 注：为方便编程，可采用另一记法：令 n’ = 2n 为偶数， 这时 ，有 复化求积公式的余项和收敛的阶 我们知道,三个求积公式的余项分别为 单纯的求积公式 复化求积公式的每个小区间 则复合梯形公式的余项为 由于 即有 例1 对于函数f(x)=sinx/x，给出n=8的函数表，试用复化梯形公式和复化辛普森公式 计算积分 x f(x) 0 1/8 1/4 3/8 1/2 5/8 3/4 7/8 1 1 0.9973978 0.9896158 0.9767267 0.9588510 0.9361556 0.9088516 0.8771925 0.8414709 解 将积分区间[0,1]划分为8等分，用复化梯形公式求得 而将积分区间[0, 1]划分为2×4等分，用复化辛普森公式求得 比较上面两个计算结果T8与S4，它们都需要提供9个点上的函数值，然而精度却差别很大，同积分准确值I=0.9460831比较，应用复化梯形公式计算的结果T8=0.9456909只有2位有效数字，而应用复化辛普森公式计算的结果S4= 0.9460832却有6位有效数字. 为了利用余项公式估计误差，要求f(x)=sinx/x的高阶导数，由于 所以有 于是 复化梯形公式误差为 复化辛普森公式误差为 例2 利用复化梯形公式计算 使其误差限为0.5*10-6，应将区间[0, 1]几等分? 解 利用例1的结果 取n=236可满足要求. 由复化梯形公式的余项得 6 因此只需将区间[0, 1]四等分，即取9个点 n3.4 取n=4 使用复合Cotes公式，只需要将区间1等分，即四小等分，五个点。 前面用复化梯形公式计算此题，满足相同的精度需要将区间[0, 1]划分17等分，可见复化辛普森公式的精度的确比复化梯形公式精度高同样也可用 | S4m-S2m |ε