梯形多步法和辛普森积分.ppt

华南师范大学数学科学学院 谢骊玲 例7.7、7.8 分别用组合梯形公式、组合辛普森公式求函数 在区间[1,6]上的定积分 考察积分区间等距划分为10,20,40,80,160个子区间的情况 M h T(f,h) ET(f,h)=O(h2) S(f,h) Es(f,h) =O(h4) 10 0.5 8-08020 0.25 8-08040 0.125 8-08080 0.0625 8-080160 0.03125 8-080积分真实值 I =8.18347920766273… 由上表可见，当 h 减半时，组合梯形公式的误差项序列ET( f ,h)的衰减因子约为1/4；而组合辛普森公式的误差项序列ES( f ,h)的衰减因子约为1/16，这验证了推论7.2和推论7.3中关于两个公式的误差阶分别为O(h2)和O(h4)的结论 华南师范大学数学科学学院 谢骊玲 利用误差阶确定区间划分个数 例7.9、7.10 计算 M 和步长 h，使得组合梯形公式和组合辛普森公式对逼近定积分 的误差ET( f ,h)和ES( f ,h)小于5×10-9 组合梯形公式：M=22822，步长h=5/22822=0.000219086846 组合辛普森公式：M=113，步长h=5/113=0.02212389381 可见，使用227次f(x)求值的组合辛普森公式与使用22823次f(x)求值的组合梯形公式得到同样的精度。前者的函数求值次数只是后者的1% 华南师范大学数学科学学院 谢骊玲 数值积分精度与函数求值、区间划分的关系 从各阶闭型N－C公式来看，函数求值的次数越多，则逼近的精度越高 但从理论上可证明M≥8的N－C公式不稳定，不能用来求解积分近似值 从组合梯形公式和组合辛普森公式计算的积分近似值来看，划分的子区间数越多，则逼近的精度越高 但从例7.9、7.10可发现，仅仅通过区间划分的方法提高精度的速度很慢 华南师范大学数学科学学院 谢骊玲 区间的划分方法 假如采用划分子区间的方式来提高精度。如何选择子区间的数目？ 采用二分区间的方法：开始时是一个区间，对分成2个子区间，再将2个子区间各自二分得到4个子区间，…，不断试验直至得到想要的精度 这个过程生成一个梯形公式的序列{T(J)} 华南师范大学数学科学学院 谢骊玲 T(0)为20=1个梯形的面积 T(1)为21=2个梯形的面积 T(2)为22=4个梯形的面积 T(3)为23=8个梯形的面积 华南师范大学数学科学学院 谢骊玲 连续梯形公式 定理7.4 设J≥1，点{xk=a+kh}将[a,b]划分为2J=2M个宽度为(b-a)/2J的子区间。梯形公式T( f ,h)和T( f ,2h)满足如下关系： 梯形公式序列 记T(0)=(h/2)(f(a)+f(b))，它是步长为h=b-a的梯形公式。对于所有J≥1，记T(J)=T(f,h)，其中T(f,h)是步长为h=(b-a)/2J的梯形公式 华南师范大学数学科学学院 谢骊玲 递推梯形公式 推论7.4 由T(0)=(h/2)(f(a)+f(b))开始，梯形公式序列{T(J)}可由以下递推公式生成： 其中J=1,2,…，h=(b-a)/2J，{xk=a+kh} 例7.11 华南师范大学数学科学学院 谢骊玲 例7.11 用连续梯形公式计算如下积分 的逼近T(0), T(1), T(2)和T(3). 只需计算9个点的函数值，而且是根据递推需要逐渐增加计算的 华南师范大学数学科学学院 谢骊玲 递推辛普森公式 定理7.5 设{T(J)}为由推论7.4产生的梯形公式序列，如果J≥1，且S(J)为区间[a,b]的2J个辛普森公式，则S(J)和T(J-1)与T(J)满足关系式： 例7.12 华南师范大学数学科学学院 谢骊玲 例7.12 用连续辛普森公式计算如下积分 的逼近S(1), S(2)和S(3). 利用连续梯形的计算结果进行组合，得到连续辛普森公式，误差阶提高了二阶，由O(h2)提高到O(h4) 华南师范大学数学科学学院 谢骊玲 华南师范大学数学科学学院 谢骊玲 华南师范大学数学科学学院 谢骊玲 第7章 数