几种常用数值积分方法的比较汇总

1、学科分类号110.3420GUIZHOUNORMALCOLLEGE本科毕业论文题目一几种常用数值积分方法的比较姓名潘晓祥学号1006020540200院（系）数学与计算机科学学院专业数学与应用数学年级2010级指导教师雍进军职称二O一四年五月贵州师范学院本科毕业论文（设计）诚信声明本人郑重声明：所呈交的本科毕业论文（设计），是本人在指导老师的指导下，独立进行研究工作所取得的成果，成果不存在知识产权争议，除文中已经注明引用的内容外，本论文不含任何其他个人或集体已经发表或撰写过的作品成果。对本文的研究做出重要贡献的个人和集体均已在文中以明确方式标明。本人完全意识到本声明的法律结果由本人承担。本科毕

2、业论文作者签名：年月日贵州师范学院本科毕业论文（设计）任务书毕业设计题目几种常用数值积分方法的比较作者姓名潘晓祥学号1006020540200年级2010级所属学院数学与计算机科学学院专业数学与应用数学班级四班指导教师签名雍进军讲师职称讲师开题日期2013年7月10日主要目标1了解什么数值积分基本思想和一些常用的数值积分方法；2对各种数值积分方法的误差以及代数精度进行分析；3对各积分方法进行比较总结出优缺点。主要要求通过对几种常用的数值积分方法进行了的分析，并用这几种方法对被积函数是普通函数做了数值积分，并在计算机上进行实验。数值积分是计算方法或数值分析理论中非常重要的内容，数值积分方法也是解

3、决实际计算问题的重要方法，对几种常用数值积分方法的分析很必要。贵州师范学院本科毕业论文（设计）开题报告书主要内容本文通过对复化求积公式，Newton-Cotes求积公式，Romberg求积公式，图斯型求积公式进行分析讨论并在计算机上积分实验，从代数精度，求积公式误差等角度对这些方法进行分析比较，并总结出每种求积分法的优缺点以及实用性。贵州师范学院本科毕业论文(设计)开题报告书论文题目几种常用数值积分方法的比较作者姓名潘晓祥学号1006020540200年级2010级所属学院数学与计算机科学学院专业数学与应用数学班级数本(4)班指导教师姓名雍进军职称讲师预计字数5000.00字题印性质应用研究日

4、期2013年7月05S选题的原由：研究意义：数值积分是数学上的重要课题之一，是数值分析中的重要内容之一，也是数学的研究重点并在实际问题及应用中有着广泛的应用常用于科学与工程的计算中，如涉及到积分方程，工程计算，计算机图形学，金融数学等应用科学领域都有着相当重要的应用，所以研究数值积分问题有很重要的意义数值积分是研究如何求出一个积分的数值这一课题的起源可追溯到古代，其中一个突出的例子是希腊人用内接与外接正多边形推算出照面积的方法也正是此法使阿基米德得以求出n值得上界与下界，若干世纪以来，尤其是十六世纪后，已提出了多种数值积分方法，其中有矩形求积法，内插求积法，牛顿科特斯公式，复化求积公式，龙贝格

5、求积公式，高斯型求积公式但各种方法都有特点，在不问的情况下试用程度不同，我们将着重从求积公式的代数精度和余项等角度对这些方法进行分析比较研究动态：这些年来，有关数值积分的研究已经成为一个很活跃的研究领域，历史上，阿基米德，牛顿，欧拉，高斯，切比雪夫等人都对此有过贡献研究出各种各样的数值求积公式，但一个好的数值求积公式应该满足：计算简单，误差小，代数精度高我们将对矩形求积法，内插求积法，牛顿科特斯公式，化求积公式，贝格求积公式，斯型求积公式进行比较对数值求积公式能有进一步的了解和学习-主要内容：1数值积分方法的基本思想2几类常用数值积分方法的基本分析2.1 NewtonCotes求积公式2.2

6、复化求积公式2.3 Romberg求积公式2.4 高斯型求积公式3几类数值积分方法的简单比较评述4利用MATLAB编程应用对几类求积算法的分析比较研究方法：本论文主要通过对相关文献和书籍的参考，合自己的见解，纪化求积公式，Newton-Cotes求积公式Romberg求积公式，高斯型求积公式进行讨论并进行上机实验，从代数精度，求积公式误差等角度对这些方法进行分析比较完成期限和采取的主要措施：本论文计划用6个月的时间完成，阶段的任务如下：(1) 7月份查阅相关书籍和文献；(2) 8月份完成开题报告并交老师批阅；(3) 9月份完成论文初稿并交老师批阅；(4) 10月份完成论文二搞并交老师批阅；(5

7、) 11月份完成论文三稿；(6) 12月份定稿.主要措施：考相关书籍和文献，合自己的见解，老师的指导下和同学的帮助下完成主要参考文献及资料名称：1吴治陆金甫数学分析基础(第二版)M.北京：等教育出版社.2010.72胡祖炽.林源渠.数值分析M北京：等教育出版社.1986.33薛毅.数学分析与实验M北京：业大学出版社2005.34徐士良.数值分析与算法M.北京：械工业出版社2007.15王开荣.杨大地.应用数值分析M北京:等教育出版社2010.76杨一都.数值计算方法M.北京:等教育出版社.2008.47韩明.王家宝.李林.数学实验(MATLAB版M.上海：济大学出版社2012.18圣宝建.关于

8、数值积分若干问题的研究J.南京信息工程大学.2009.05.01:429刘绪军.几种求积公式计算精确度的比较J.南京职业技术学院.2009.10史万明.吴裕树.孙新.数值分析M.北京理工大学出版社.2010.4.指导教师苣见:签名：开题报告会纪要时间2013年8月26日地点宁静楼229教师办公室与会人员姓名职务（职称）姓名职务（职称）姓名职务（职称）蕤进军导师（讲师）邓喜才副教授李晟副教授组长会议记录摘要：指导小组针对课题二次函数性质的应用提问了以下问题以及报告人的回答：雍老师问：扰择此题目的目的？潘晓祥答：随着计算机和计算方法的飞速发展，几乎所有学科都走向定量化和精确化，计算数学中的数值计算

9、方法则是解决“计算”问题的桥梁和工具。邓老师何：对这个问题进行研究有什么实际的意义？潘晓祥答：计算方法既有数学类课程中理论上的抽象性和严谨性，又有实用性和实验性的技术特征，计算方法是一门理论性和实践性都很强的学科在科学研究和工程技术中都要用到各种计算方法.例如，在航天航空、地质勘探、汽车制造、桥梁设计、天气预报和汉字字样设计中都有计算方法的踪影。李老师何：对这个问题你有什么自己的看法？潘晓祥答：随着计算机技术的迅速发展和普及，现在计算方法课程几乎已成为所有理工科学生的必修课程我们知道，计算能力是计算工具和计算方法的效率的乘积，提高计算方法的效率与提高计算机硬件的效率向样重要科学计算已用到科学技

10、术和社会生活的各个领域中所以，研究数值计算方法可以让数学的应用更大更广。会议主持人签名：记录人签名：在日口指导小组意见负责人签名：年月日学院意见负责人签名：年月日贵州师范学院数学与计算机科学学院指导教师指导本科毕业论文情况登记表论文（设计）题目几种常用数值积分方法的比较用学生姓名潘晓祥学号1006020540200年级2010级所属学院数计学院专业数学与应用数学班级四班指导教师姓名雍进军职称讲师学历硕士指导时间指导地点指导内容指导教师签名备注2013年06月10日致远楼416论文扰题，资料准备面授2013年06月22日网上确定毕业论文扰题电子邮件2013年06月26日网上怎样撰写毕业论文开题报

11、告电子邮件2013年06月28日网上指导学生撰写开题报告电子邮件2013年07月14日网上帮助学生查找有关参考文献电子邮件2013年07月17日手机女M可构思自己的毕业论文手机飞信2013年08月21日手机听取学生毕业论文写作进展情况汇报手机飞信2013年08月28日网上解答学生在论文写作中遇到的疑惑电子邮件2013年9月09日网上帮助学生查找有关参考文献电子邮件2013年11月28日网上女M可规划自己的论文电子邮件2013年12月04日手机怎样写好论文引言手机飞信2013年12月08日网上怎样写好论文引言电子邮件2013年12月12日网上怎样写论文摘要电子邮件2013年12月16日网上怎样疑

12、取论文关键词电子邮件2013年12月20日网上怎样编辑论文中的公式电子邮件2014年01月05日手机督促学生在寒假中写好论文的初稿电子邮件2014年02月27日宁静楼219检查学生论文完成情况面授2014年03月03日宁静楼219对学生的论文初稿提出修改时意见面授2014年03月07日宁静楼219解答学生在修改时的困惑面授2014年03月11日宁静楼219指导学生修改论文面授贵州师范学院数学与计算机科学学院本科毕业论文（设计）交叉评阅表学院（盖章）：学号1006020540200姓名潘晓祥专业数学与应用数学论文（设计）题目几种常用数值积分方法的比较班级四班指导教师意见评语：该向学在论文撰写过程

13、中对相关文献阅读范围广泛，方法正确，内容完整能综合运用所学知识分析和解决实际问题。毕业论文撰与过程中态度端正，勤奋刻苦。论文研究了Newton-Cotes求枳公式、复化求积公式、Romberg积分、高斯积分方法，通过算例分析，得出几种常用数值积分方法是解决头际计算问题的重要方法。论文结构合理，符合逻辑，文章层次分明，语言准确，文字通顺达到规范性要求，建议作为学士论文答辩。成绩：（满分100分）指导教师（签名）：年月日评阅教师意见评语：该向学具备较好的基础理论与专业知识，学习态度认真，阅读教师指定的参考资料、文献，较好的完成了任务书规定的工作量。论文研究了Newton-Cotes求积公式、复化求

14、积公式、Romberg积分、图斯积分方法通过算例分析，得出几种常用数值积分方法是解决实际计算问题的重要方法。轮文结构合理，符合逻辑，文章层次分明，语言准确，文字通顺达到本科毕业论文相关要求。同意参加答辩。成绩：（满分100分）评阅教师（签名）：年月日贵州师范学院本科毕业论文答辩记录表论文题目几种常用的多项式插值方法作者姓名潘晓祥学号1006020540200所属学院数计专业数学与应用数学年级2010级班级本科（四）班雍进军 讲师指导教师姓名、职称答辩会纪要时间2014年5月11日地点致远楼406答辩小组成员姓名职务（职称）姓名职务（职称）姓名职务（职称）左羽教授在忠伟副教授廖玉梅讲师答辩中提出

15、的主要问题及回答的简要情况记录：1自己做的有哪些？答：第11页至第12页，总结进行比较。2 .程序运行过没有？答：运行过。3 .20页程序代码中，if后的是什么符号？答：连接作用的符号。4解释一下什么时候用分号，什么时候不用？答：回答不清。5.摘要中央文拼与有错答辩后修改答辩小组负责人签名：左羽记录人签名：梅林林2014年5月11日答j_L辩小组意见评语：该生能在规定时间叙述论文的主要内容，对提出的问题一般能回答，无原则错误。答辩小组经过充分讨论，根据该生论文质量和答辩中的表现，向意评定论文成绩为“中等”。评定成绩：77负责人（签名）：左羽2014年5月11日Abstract21前言32数值积

16、分方法的基本思想33几类常用数值积分方法的简单分析51 .1NewtonCotes求积公式52 .2复化求积公式63 .3Romberg求积公式74 .4高斯型求积公式84几类数值积分方法的简单比较评述95利用MATLAB程应用对几类求积算法的分析比较10结束语错误！未定义书签。致谢14附录16贡州师范学院毕业论文（设计）我们在求函数的积分时，往往因为原函数非常复杂以至于难以求出或用初等函数表示，这让我们计算起来非常困难，所以我们只能想办法求它的近似值，因此直接借助牛顿一莱布尼兹公式计算定积分的情况是非常少见的。这时候数值积分就是解决这种问题的一种很好很有效的方法。本文从数值积分问题的产生出发

17、，详细介绍了一些数值积分的常用方法（NewtonCotes求积公式，复化求积公式，Romberg求积公式，高斯型求积公式）并对其进行了简要的分析，在探讨了这些数值积分算法的优缺点的理论之外，我们还将这些数值积分算法在计算机上通过matlab软件编程实现应用，并分别用各自求积公式进行运算，以此来分析比较各种求积公式的代数精度和计算误差。关键词：数值积分；求积公式；代数精度贡州师范学院毕业论文（设计）Abstractfunctionisverycomplexthatitisdifficulttofindtheelementaryfunctions,whichmakesuWeinthefunctio

18、nfortheintegration,oftenbecausetheoriginalsverydifficulttocalculate,sowecanonlythinkofawaytofindtheapproximatevalue,thusdirectlywithNewton-Leibnizformulacalculatingdefiniteintegralsituationisveryrare.Whennumericalintegrationistosolvethisprobleminaveryeffectivemethod.Fromthenumericalintegrationproble

19、m,introducessomemethodsofnumericalintegration(Newton-Cotesquadratureformula,compositequadratureformulas,Longbeilatticequadratureformula,Gausstypequadratureformulas)andhascarriedonbriefanalysis,discussestheadvantagesanddisadvantagesofthesenumericalintegrationalgorithmtheory,wewillthesenumericalintegr

20、ationalgorithminthecomputerbyMATLABsoftwareprogrammingapplication,andseparatelywiththeirrespectivequadratureformulaforcomputing,inordertoanalyzethealgebraiccalculationprecisionanderrorcomparisonofvariousquadratureformulas.Keywords:Numericalintegration;Calculationmeth;numericalanalysis#员州加范学院毕业论文(设计)

21、、八1月I言微积分的发明是世界数学史上一项辉煌的成就。但在实际求积问题的时候，求解积分却有着非常多局限性。比如对于定积分Ff(x)dx在求某函数的定积分时，在一定条件下，虽然有牛顿莱布里茨公式仁阡(x)dx二F(b)一F(a)可弋a以计算定积分的值，但在很多情况下f(x)的原函数不易求出或非常复杂。被积函数f(x)的原函数很难用初等函数表达出来，例如()二沁(，等；有的X函数f(x)的原函数F(x)存在，但其表达式太复杂，计算量太大，有的甚至无法有解析表达式。因此能够借助牛顿莱布尼兹公式计算定积分的情形是不多的。另外，许多实际问题中的被积函数往往是列表函数或其他形式的非连续函数对这类函数的定积

22、分，也不能用不定积分方法求解，只能设法求其近似值。因此，探讨近似计算的数值积分方法是有明显的实际意义的即有必要研究定积分的数值计算方法，以解决定积分的近似计算。而数值积分就是解决此类问题的一种有效的方法，它的特点是利用被积函数f(x)在一些节点上的信息求出定积分的近似值。在很多实际应用中，只能知道积分函数在某些特定点的取值比如天气测量中的气温、湿度、气压等，医学测量中的血压、浓度等等。通过研究，我们将会更熟练掌握一些数值积分方法去计算一些特定条件的数值计算，以便我们得到自己想要的结果。2数值积分方法的基本思想在数学分析中，计算连续函数f(x)在区间a,b上的积分是通过f(x)的原函数F(x),

23、由下列定积分公式baf(x)dx二F(b)-F得到的。但由于大量被积函数的原函数不能用初等函数表示，因此，很难用求原函数的公式4(x)dx=F(b)-F(a)得到积分；有些被积函数f(x)不是明显知道的，例如由数值表给出它的离散值，或者是它被定义为某个微分方程的解，而这个微分方程是不能显示解出的。这说明按f(x)dx=F(b)一F(a)公式计算La定积分是有很大局限性的。因而常常采用在电子计算机上很有效的数值积分方法。我们从定积分的定义bnaf(x)dx:(b-a)C(n,k)f(xk)k=0bn出发。推导出两个简单的数值积分公式。Jaf(x)dx=n(kE)AXk式的几何意义，-ka就是把整

24、块曲线梯形的面积积分成若干个小曲边梯形面积的和，当无限细分时这个和取极限就是真正曲边梯形面积。去掉取极限这一步，用有限个小曲边梯形面积的和，代替整块的曲边梯形面积，从而求得一个近似值，这就是数值积分的基本思想。根据小区间的不同分割方法和各分点f(z)值的不同选择，就得到不同的数值积分公式。数值求积公式是取lab1上若干个点Xk处的高度f(Xk)，通过加权人后，再n求和、Akf(Xk)k=0从而得到积分的近似值。数值求积公式写成一般形式nf(x)dxAAJ(Xk)k=0式中Xk称求积节点，Ak称求积系数，也称伴随节点Xk的权。当积分区间a,bl确定后，求积系数，仅仅与节点Xk的选取有关，而不依赖

25、被积函数f(x)的具体形式。记#贵州师范学院毕业的文（设计）Rff(x)dxA.f(xk)-k=0#员州师范学院毕业论文(设计)把RIf1称为求积公式的截断误差或余项。数值求积方法的特点是直接利用积分区间l.a,b1上一些离散节点的函数值进行线性组合来近似计算定积分的值，从而将定积分的计算归结为函数值的计算，这就避开了牛顿莱布尼兹公式需要寻求原函数的困难，并为计算机求积分提供了可行性。3几类常用数值积分方法的简单分析3.1Newton-Cotes求积公式常用的梯形公式和Simpsor公式是低阶的牛顿柯特斯公式，牛顿柯特斯公式是积分区间上等距节点的插值求积公式。插值求积公式在积分区间上，所取节点

26、是等距时称为牛顿柯特斯公式，即nf(x)dx:(b-a)二C(n,k)f(xk)k=0其中c(n,k)为Cotes求积公式的系数，是n和k的函数。当n=1时，为梯形公式：f(x)dx9ff(b)2梯形公式的代数精度为1，有两个积分节点。当n=2时，为Simpson公式：b(b-a)af(x)dx定1f(a)4f(；f(b)6Simpson公式的代数精度为3，有三个积分节点由于只增加一个节点，其代数精度增加2，由此可知，Simpson公式比梯形公式代数精度高。当n=4时，NewtonCotes求积公式为Cotes公式：:f(x)dx烤7f(a)32f(专)12f(宁)32f(宁)7f(b)New

27、ton-Cotes公式的代数精度为5,有5个积分节点。所以对于Newton-Cotes积分公式，n为偶数时的代数精度要比n为奇数时的积分公式效果比较优越。但并不是n的值越大越好，当n过大时(n=8)，求积公式的数值稳定性不好。3.2 复化求积公式由于Newton-Cotes的节点n越大对应的精度就越高，但是n=8时公式的数值是不稳定的，因此就不能用增加求积节点的方法来提高精度，因此，我们常常将求积区间a,b分成若干小区间，然后在每个小区间上采用数值稳定的Cotes公式求小区间上的积分，然后把每个小区间上的结果加起来作为原定积分的近视值，这种方法构造的求积公式就叫做复化求积公式。常用的复化求积公

28、式有：复化梯形公式：b,a1nj1变步长梯形公式为：=+f(Xo)”)=猪12Mn二宁1f(X2jj)2nj#复化Simpson公式：b-a人盲“心)忖眸)变步长复化Simpson公式：nb_an心S2n八Ikf(X。)fgn)4、f(X2j2f(X2k)k：i6nkdkT3.3 Romberg求积公式Romberg积分方法也叫做逐次分半加速法，它是在复化梯形公式误差估计的基础上，应用线性外推的方法构造出的一种加速算法。将积分区间分成n等分和2n等分时，求得积分近似值和T.，并没有误差估计式1|.T2n3(T2n-Tn)3积分近似值T2n的误差大致等于3仃2人)，当用-)对T?.进行修正时，3

29、3-(T2n-Tn)与T2n之和比T?n更接近于真值I,故陶-)是对T?n误差的一种补偿，33因此可以期望下式是一个更好的结果，即T二T2n.(J-)二今n2-；Tn333下面说明T即是分成n等分时Simpson公式的值5。将复化梯形公式f(a)2、f(xOf(b)梯形变步长求积公式T2ndin4f(Xf22心k2代人上式T表达式得-hn4n4T=;Jf(a)+4送f(x-)+2送f(xj+f(b)6-心迁心这就是说，用梯形法二分前后两个梯形值和丁加作线性外推，结果得到Simpson法的积分值&。将误差由。(忤)变为。(币)，从而提高了逼近精度。再考察Simpson法。其截断误差与h成正比，因

30、此，若将步长折半,则误差减至，即有16I-En. -Sn16一一-7贵州师范学院毕业的文（设计）WWI：16S2-Sn1515不难验证，上式右端的值其实等于Cn，就是说，用Simpson法二分前后的两个积分值&与s2n,按上式再作线性外推，结果得到柯特斯法的积分值G，即有161CnS2nSn1515这时将误差由。（的变为。（的，逼近精度又一次得以提高。同样的方法，依据柯特斯法的误差公式，可进一步导出下列龙贝格公式641“nC2nCn6363用逼近积分值的误差为0（M，这样Romberg公式将误差由。（忖变为o（h），逼近精度再次得以提高。Romberg公式有7次代数精度，这表明该公式不是牛顿柯

31、特斯公式。在步长二分的过程中运用Sn、Cn、Rn表达式加工三次，就能将粗糙的积分值Tn逐步加工成精度较高的Romberg值R”，或者说，将收敛缓慢的梯形值序列入加工成收敛迅速的Romberg值序列R”，这种加速方法称Romberg算法。3.4高斯型求积公式前面介绍的n-1个节点的Newton-Cotes求积公式，其特征是节点是等距的。这种特点使得求积公式便于构造，复化求积公式易于形成。但同时也限制了公式的精度。n是偶数时，代数精度为n1，n是奇数时，代数精度为n；我们知道n1个节点的插值型求积公式的代数精确度不低于no能不能在区间a.b1上适当选择个节点Xo,X1,X2.x使插值求积公式的代数

32、精度高于n呢？答案是肯定的，适当选择节点，可使公式的精度最高达到2n1，这就是所学的高斯型求积公式。不失一般性，将求积公式/0f（x）dx八Axf（Xk）的求积区间La,b转换成T,1ak=0的形式。对任意求积区间la,b作变换D-a,a夕ttb a _L a b *b a _L a b j 22 b 8 ?(t) dt1-1,1上，这时可以变换到区间af(X)dX=乜b-a2-baf出二2其中-仕）二f（batabt）。22高斯-勒让德求积公式在这里简称高斯公式，它是在区间1-1,1上进行讨论的。4几类数值积分方法的简单比较评述由于我们在计算实际问题是往往要考虑到代数精度和计算量，所以不同类

33、型的求积公式有着不同的特点：Simpson积分方法和梯形积分方法虽然计算简便，但是精度比较差，不理想。但对于光滑性较差的被积函数有时会比高精度的积分方法更为有效。特别是梯形积分方法对被积函数是周期函数的求积效果更为突出。n7时，Newton-Cotes公式是不稳定的，然而复化梯形公式和复化Simpson公式不仅保留了低阶公式的优点还能够获得比较较高的精度，所以在实际计算中应用得最为广泛。Romberg积分方法的算法简单，方便编程的实现。收敛速度快、计算精度较高，但是计算量较大。Gauss积分方法的精度较高，数值稳定、收敛速度较快，但因为其节点不规则，计算比较麻烦5利用MATLAB编程应用对几类

34、求积算法的分析比较在简单的认识积分方法比且理论比较之后，则要进行数学实验进行验证，因此就要通9员州师范学院毕业论文(设计)过matlab软件对各种积分方法进行编程并运算，然后对其各种方法的运算结果进行分析比较，掌握和理解各方法的优缺点。规定各个程序都以I=0皿dx为例子进行运算。原积分的精确值为电XisinxoXdx”946083070367183例分别用不同的方法计算积分I二，沁dx，并作比较。x用以上介绍的几类积分方法分别计算积分，得出误差，并进行比较：1、用Newton-Cotes公式当n=1时，即用梯形公式，用程序一(程序见附录)在MATLAB命令窗口中输入NCotes(0,1,1,2

35、)得I.92703549240395Ir=0.01904757796323当n=2时，即用Simpson公式，用程序一(程序见附录)在MATLAB命令窗口中输入NCotes(0,1,2,2)得I肝0.94614588227359Ir=0.000062811906407当n=4时，即用科特斯公式，用程序一(程序见附录)在MATLAB命令窗口中输入NCotes(0,1,4,2)得I-0.94608300406367R=0.0000000663035132、用复化梯形公式令h=1/8=0.125,用程序二(程序见附录)在MATLAB命令窗口中输入trapri(0,1,8)，得J;叱dxf(0)+2f

36、(h)+.+f(7h)+f(l)=0.9456908635270x2R=0.0003922068401823、用复化Simpson公式令h=1/8=0.125,用程序三(程序见附录)在MATLAB命令窗口中输入simpr1(f,0,1,8),得1SnXdx：h(f(0)4f(h)f(7h)b：；,2f(2h)f(6h)IHf(10x3=0.94608308538495R=0.0000000150177674、用Romberg公式用程序四(程序见附录)在MATLAB命令窗口中输入romber(T,0,1,5,0.5*(10/(-8)得(沁dx4.94608307036718|R=0.000000

37、000000002X5、用高斯勒让德求积公式令x=(t+1)/2，|=fR49Zdt=t+1(1)用2个节点的Gauss公式10.94604115827633(2)用3个节点的Gauss公式，用程序五(程序见附录)在MATLAB命令窗口中输入GuassLegendre(0,1,2,2)得1-0.946083134078473|R=0.000000063711290算法比较：1 .原积分的精确值为：1 sinx1 =dx=0.946083070367183x2 .由例题的各种求积算法可知：(1)对Newton-cotes公式，当n=1时只有1位有效数字,当n=2时有3位有效数字，当n=4时有7位

38、有效数字。(4)用复化梯形公式有2位有效数字，对复化Simpson公式有7位有效数字。(5)用复化梯形公式，对积分区间0,1二分11次用了2049个函数值，才可以得7位有效数字。(6)用Romberg公式对区间0,1二分3次用了9个函数值，就可以得到7位有效数字；二分4次用了14个函数值，却可以得到14位有效数字。(7)用高斯勒让德求积公式仅仅用了3个函数值，就能得到比较精确的6位有效数字。11贡州师范学院毕业论文（设计）结束语本文主要研究了常用的几类数值积分的求积算法并通过例题计算积分进行分析比较。Newton-Cotes积分方法是一种非常普遍的积分方法，然而梯形积分方法的误差最大，近似效果

39、最差，Simpson积分方法的精度比梯形积分方法高了一个数量级；Cotes积分方法精度比Simpson积分方法高两个数量级。则Cotes代数精度比较高。由此可知一般情况下，积分公式代数精度越高，计算精度也越高。但是高阶的Cotes积分方法收敛性没有保证，因此实际应用中很少用。复化梯形积分方法比梯形积分方法精度高，同样的，复化Simpson积分方法比Simpson积分方法精度高，高了差不多7个数量级，所以复化积分方法比较优越。Romberg积分方法收敛速度快、计算精度较高，但是计算量较大。Gauss积分方法精度高、数值稳定、收敛速度较快，但是计算麻烦。经研究可以知道Newton-Cotes方法的

40、代数精度越高，数值积分的效果越好、越精确。当积分区间比较大的时候，积分数值不稳定，这个时候可以利用复化积分方法效果会更好；Romberg积分方法可以利用变步长复化积分公式得到更为精确的数值结果，是比较好的积分方法。高斯求积方法精确度高，收敛性快，比其他积分方法优越。具有很广泛15的运用。参考文献1关治.陆金甫.数学分析基础（第二版）M.北京：等教育出版社.2010.72胡祖炽.林源渠.数值分析M北京:等教育出版社.1986.33薛毅.数学分析与实验M北京:业大学出版社2005.34徐士良.数值分析与算法M.北京版工业出版社2007.15王开荣.杨大地.应用数值分析M北京:等教育出版社2010.

41、76杨一都.数值计算方法M.北京:等教育出版社.2008.47韩明.王家宝.李林.数学实验（MATLAB）版M.上海:济大学出版社2012.18圣宝建.关于数值积分若干问题的研究J.南京信息工程大学.2009.05.01.:429刘绪军.几种求积公式计算精确度的比较J.南京职业技术学院.2009.10史万明.吴裕树.孙新.数值分析M.北京理工大学出版社.2010.4.行文至此，我的这篇论文已接近尾声；岁月如梭，我四年的大学时光也即将敲响结束的钟声。离别在即，站在人生的又一个转折点上，心中难免思绪万千，一种感恩之情油然而生。首先感谢贵州师范学院四年来对我的培养，是博学的老师们教会了我学习的方法、

42、锻炼了我思考的能力、指明了我未来奋斗的方向，从而使我进一步明确了人生的目标。其次，我要感谢我的指导老师一雍进军老师，他的严谨细致、一丝不苟的作风一直是我工作、学习中的榜样；他的循循善诱的教导和不拘一格的思路给予我无尽的启迪。在撰写整个毕业论文的过程当中，他为我们考虑到了每一个细节，从开题报告到毕业论文的拟定修改上，雍老师更是不厌其烦的为我们做好每一步的细心指导。对此，我表示衷心地感谢。没有雍老师，我的论文也不可能这么顺利的完成。同时，我也要感谢每一位给过我帮助的老师和同学，在我撰写论文的过程当中同样给了我大量有益的建议，在此一并向他们表示真诚的感谢，感谢他们对我的支持和帮助。最后感谢这篇论文所

43、涉及到的各位学者，本文引用了数位学者的研究文献，如果没有各位学者的研究成果带给我的的帮助和启发，我将很难完成本篇论文的写作。由于我的学术水平有限，所写论文难免有不足之处，恳请各位老师和学友批评指正。最后，衷心感谢评阅论文及参加答辩的各位老师！员州师范学院毕业论文(设计)1 NewtonCotes求积公式的MATLAB实现先用M文件定义一个名为fl.m的函数：%i是要调用第几个被积函数g(i)，x是自变量functionf=f1(i,x)g(1)=sqrt(x);ifx=0g(2)=1；elseg(2)=sin(x)/x;endg(3)=4/(1+xA2);仁g(i);程序一：functionC

44、,g=NCotes(a,b,n,m)%a，b分别为积分的上下限；%n是子区间的个数；%m是调用上面第几个被积函数；%当n=1时计算梯形公式；当n=2时计算辛浦生公式，以此类推;i=n；h=(b-a)/i;z=0;forj=0:ix(j+1)=a+j*h;s=1;ifj=Os=s;elsefork=1:js=s*k;endendr=1;ifi-j=0r=r;elsefork=1:(i-j)r=r*k;endendifmod(i-j),2)=1q=-(i*s*r);elseq=i*s*r;endy=i;fork=0:iifk=jy=y*(sym(t)-k);endendl=int(y,0,i);C

45、(j+1尸l/q；z=z+C(j+1)*f1(m,x(j+1)；endg=(b-a)*z1)当输入a=0,b=1,n=1A=2时，即在MATLAB命令窗口输入NCotes(0,1,1,2)即可得用梯形公式的积分值和相应科特斯系数如图1.12)当输入a=0,b=1，n=2,m=2时，即在MATLAB命令窗口输入NCotes(0,1,2,2)即可得用辛浦生公式的积分值和相应科特斯系数如图1.23)当输入a=0,b=1，n=4,m=2时，即在MATLAB命令窗口输入NCotes(0,1,4,2)即可得用科特斯公式的积分值和相应科特斯系数如图1.317贵州师范学院毕业论文（设计）CommandWind

46、owHCotesLL2)8293248040994423/9007199254740992ans=(1/2,1/2)8293248040994423/90071992547409920.92073549240395CommandWindowNCotes(01.2.2)8522124485690909/9007199254740992ans=(1/6,2/3,1/68522124485690909/9007199254740992ans=0.94614588227359图1.219员州师范学院毕业论文(设计)CommandWindowNCotes(Qb42)g=38347011581064500

47、9/405323966463344640axis=7/90,16/45,2/15,16/45,7/90383470115810645009/405323966463344640an，=0.94608300406367图1.32复化梯形求积公式的MATLAB实现通过f(x)的n1个等步长节点逼近积分:/(幼办怎(/()+/(*)+方/Im)其中,xtkh,Xo=a.程序二:functions=trapr1(f,a,b,n)%f是被积函数；%a,b分别为积分的上下限；%n是子区间的个数；%S是梯形总面积；h=(b-a)/n;s=0;fork=1:(n-1)x=a+h*k;s=s+feval(f,x);endformatlongs=h*(feval(f,a)+feval(,f,b)/2+h*s;先用M文件定义一个名为f.m的函数：functiony=f(x)ifx=0y=i;elsey=sin(x)/x;end若取子区间的个数在MATLAB命令窗口中输入trapri(