计算方法电子教案(第三章).pptVIP

cotes系数的计算： 当 n=1时，k=0,1 例3.6 用变步长的复合Simpson公式计算定积分 再将步长折半，h=0.25,由例3.4可得。 二、 Romberg算法 Romberg算法是利用复合梯形公式，在对积分区间的步长逐次折半的过程中，求得积分 的近似值序列 （其中k表示对求积区间[a,b]的对分次数,k=0,1,2, …),并且利用误差补偿的方法，逐步将 加工(递推)成具有高精度的积分近似值。 定理3.2 Gauss型求积公式的系数Ak0(k=0, 1, …, n). 证明: 由于Gauss型求积公式对任何不大于(2n+1)次的多项式精确成立,若取 §3.3.1 Gauss-Legendre求积公式 0.3478548451 0.6521451549 ±0.8611363116 ±0.3399810436 3 0.2369268851 0.4786286705 0.5688888889 ±0.9061798459 ±0.5384693101 0 4 0.1012285363 0.2223810345 0.3137066459 0.3626837834 ±0.9602898565 ±0.7966664774 ±0.5255324099 ±0.1834346425 7 0.1294849662 0.2797053915 0.3818300505 0.4179591837 ±0.9491079123 ±0.7415311856 ±0.4058451514 0 6 0.5555555556 0.8888888889 ±0.7745966692 0 2 1 ±0.5773502692 1 0.1713244924 03607615730 04679139346 ±0.9324695142 ±0.6612093865 ±0.2386191861 5 2 0 0 Ak xk n Ak xk n 表3.2 例3.4 用n=8的复化梯形公式及n=4的复化Simpson公式,计算积分 解: 只要将区间[0, 1]分为8等分,用复化梯形公式时n=8,h=0.125,对复合Simpson公式取n=4,h=0.25.计算各分点xk的函数值f (xk).由公式(3.11)及式(3.13)得下表: 0.9588510 0.5 4 0.8414709 1 8 0.9767267 0.375 3 0.8771925 0.825 7 0.9896158 0.25 2 0.9088516 0.75 6 0.9973978 0.125 1 0.9361556 0.625 5 1.0000000 0 0 f (xk) xk k f (xk) xk k 由式(3.11)及式(3.13)得 例3.5 若用复化梯形公式求 的近似值,问要将积分区间[0,1]分成多少份才能保证计算结果有4位有效数字?若用复化Simpson求积公式呢? 解: (1) 用复化梯形求积公式 即若用复化梯形公式求I(f)的近似值,需将[0,1]41等分才能保证计算结果有4位有效数字. (2) 用复化Simpson 求积公式 §3.2.3 逐次分半梯形积分公式与龙贝格积分公式 一、逐次分半梯形积分公式 1、基本步骤: 解：利用例3.4给出的结果，先取h=b-a=1,则： 将步长折半，h=0.5,则 由3.16式 T—数表 T3(0) T2(0) T2(1) T1(0) T1(1) T1(2) T0(0) T0(1) T0(2) T0(3) 0 1 2 3 T3(k-3) T2(k-2) T1(k-1) T0(k) k 梯形公式 辛普生公式 Cotes公式 龙贝格公式 0.946 083 07 0.946 083 07 0.946 083 31 0.945 690 86 3 0.946 083 00 0.946 086 93 0.944 513 52 2 0.946 145 88 0.939 793 28 1 0.920 735 49 0 T3(k-3) T2(k-2) T1(k-1) T0(k) k §3.3 Gauss型求积公式 一. 问题：在节点数目固定为n+1的条件下,能否恰当地选择节点位置和相应的系数,使求积公式 具有最大的代数精度. 首先分析一下对于n+1个节点,求积公式(3.1)可以达到的最大代数精度是多少? 设上述求积公式对所有m (m待定)次多项式 是准确的,于是有 由于a0,