改进Simpson公式及误差分析.pdfVIP

维普资讯 第28卷 第4期 高 师 理 科 学 刊 Vo1．28 No．4 2008庄 7月 JournalofScienceofTeachers CollegeandUniversity Ju1． 2oo8 文章编号 ：1007—9831(2008)04—0027—40 改进 Simpson公式及误差分析 杜跃鹏。，肖泽昌 (1．南阳理工学院 计算机系，河南 南阳473004；2．南阳理工学院 数学系，河南 南阳 473004) 摘要：给出了改进 Simpson公式的截断误差，分析了复化改进 Simpson公式的收敛阶．数值算例 验证 了理论分析的正确性． 关键词：数值积分；代数精度；截断误差；收敛阶 中图分类号：0241．4 文献标识码：A Simpson求积公式是有 3个等距节点的插值型求积公式 ，是计算定积分近似值的有效方法．但是 Simpson公式代数精度低，为了获得较高精度，需要更细的剖分和更多的求积节点．文献[3]给出了一些数 值求积公式的渐近性质，这些公式包括求积分的矩形法则、梯形法则和抛物线法则，获得了这些法则的改 进公式．文献[4]给出了Simpson公式误差的最佳估计．文献[5]利用代数精度的概念，对文献[3]中的左矩形 校正公式和梯形校正公式重新推导，简化了文献[3]的推理过程．文献[3]和文献[5]都未给出改进公式的截断 误差．本文利用文献[5]的思想方法，对文献[3]中改进Simpson公式重新推导，利用泰勒中值定理，给出了 改进Simpson公式的截断误差，分析了复化改进Simpson公式的收敛阶，给出了数值算例． 记，(，)=，(．文献2[]queansimps。n公式，(，)= (，(口)+4，()+，6()】一 · ， ()， ∈(口，6)．令 =a+ (6一口)，求 ∈(0，1)，使数值积分公式 = 6 (f(口)w(字)州卟一 ，4( ) (1) 具有尽可能高的代数精度． 容易验证，当，()：1， ， ：， ，， 时，式 (1)精确成立，当，( =X5时，有 X ： ， w( ]_ )= 4( f_51 (6一口)】，从而由式 (1)解得 = ，NNN~Nsimps。n公式 ( 6卜 )(] 容易验证，()= 时，I(x)≠ ()，说明改进 Simpson公式 (2)的代数精度是 5．改进 Simpson 公式的截断误差记为E (，)，即 )-Qs( b．-6a(f㈤w(]州 卜 ，4()(] 定理设，()在区间n【，纠上具有6阶连续导数，M=maxI。f()l}，则l(，)l≤笔 M． 收稿 日期：2008-03-01 作者简介：杜跃鹏 (1965一)，男，河南南阳人，副教授 ，硕士，从事计算数学研究．E-maiI：nitduyp@163．eom 维普资讯 28 高 师 理 科 学 刊 第28卷 证明设F(x)：f~f(t)dt，记。：_a+b，由泰勒中值定理知 ¨ )+ ：+ + + s+ + 其中 介于 和 ．]c。之间，故 ㈤ ，(]+ )()+ () )+ 字( ) ))