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

学科分类号 110.3420 本 科 毕 业 论 文 题 目 几种常用数值积分方法的比较 姓 名 潘晓祥 学 号 1006020540200 院 （系） 数 学 与 计 算 机 科 学 学 院 专 业 数学与应用数学 年 级 2010 级 指导教师 雍 进 军 职 称 讲 师 二一四年五月贵州师范学院本科毕业论文（设计）诚信声明本人郑重声明：所呈交的本科毕业论文（设计），是本人在指导老师的指导下，独立进行研究工作所取得的成果，成果不存在知识产权争议，除文中已经注明引用的内容外，本论文不含任何其他个人或集体已经发表或撰写过的作品成果。对本文的研究做出重要贡献的个人和集体均已在文中以明确方式标明。本人完全意识到本声明的法律结果由本人承担。 本科毕业论文作者签名： 年 月 日贵州师范学院本科毕业论文(设计)任务书毕业设计题目几种常用数值积分方法的比较作 者 姓 名潘晓祥学号1006020540200年级2010级所属学院数学与计算机科学学院专业数学与应用数学班级四班指导教师签名雍进军 讲师职称讲师开题日期2013年7月10日主要目标 1了解什么数值积分基本思想和一些常用的数值积分方法；2对各种数值积分方法的误差以及代数精度进行分析；3对各积分方法进行比较总结出优缺点。主要要求 通过对几种常用的数值积分方法进行了的分析，并用这几种方法对被积函数是普通函数做了数值积分，并在计算机上进行实验。数值积分是计算方法或数值分析理论中非常重要的内容，数值积分方法也是解决实际计算问题的重要方法，对几种常用数值积分方法的分析很必要。主要内容 本文通过对复化求积公式, NewtonCotes求积公式, Romberg求积公式,高斯型求积公式进行分析讨论并在计算机上积分实验,从代数精度,求积公式误差等角度对这些方法进行分析比较，并总结出每种求积分法的优缺点以及实用性。 贵州师范学院本科毕业论文（设计）开题报告书论 文 题 目几种常用数值积分方法的比较作 者 姓 名潘晓祥学号1006020540200年级2010级所属学院数学与计算机科学学院专业数学与应用数学班级数本（4）班指导教师姓名雍进军职称讲师预计字数5000.00字题目性质应用研究日期2013年7月05 日选题的原由：研究意义:数值积分是数学上的重要课题之一,是数值分析中的重要内容之一,也是数学的研究重点.并在实际问题及应用中有着广泛的应用.常用于科学与工程的计算中,如涉及到积分方程,工程计算,计算机图形学,金融数学等应用科学领域都有着相当重要的应用,所以研究数值积分问题有很重要的意义. 数值积分是研究如何求出一个积分的数值.这一课题的起源可追溯到古代,其中一个突出的例子是希腊人用内接与外接正多边形推算出圆面积的方法.也正是此法使阿基米德得以求出值得上界与下界,若干世纪以来,尤其是十六世纪后,已提出了多种数值积分方法,其中有矩形求积法,内插求积法,牛顿科特斯公式,复化求积公式,龙贝格求积公式,高斯型求积公式.但各种方法都有特点,在不同的情况下试用程度不同,我们将着重从求积公式的代数精度和余项等角度对这些方法进行分析比较.研究动态:这些年来,有关数值积分的研究已经成为一个很活跃的研究领域,历史上,阿基米德,牛顿,欧拉,高斯,切比雪夫等人都对此有过贡献.研究出各种各样的数值求积公式,但一个好的数值求积公式应该满足:计算简单,误差小,代数精度高.我们将对矩形求积法,内插求积法，牛顿科特斯公式,化求积公式,贝格求积公式,斯型求积公式进行比较.对数值求积公式能有进一步的了解和学习.主要内容：1 数值积分方法的基本思想2 几类常用数值积分方法的基本分析2.1 NewtonCotes求积公式2.2 复化求积公式2.3 Romberg求积公式2.4 高斯型求积公式3 几类数值积分方法的简单比较评述4利用MATLAB编程应用对几类求积算法的分析比较研究方法：本论文主要通过对相关文献和书籍的参考,合自己的见解,复化求积公式,NewtonCotes求积公式,Romberg求积公式,高斯型求积公式进行讨论并进行上机实验，从代数精度,求积公式误差等角度对这些方法进行分析比较.完成期限和采取的主要措施：本论文计划用6个月的时间完成,阶段的任务如下:(1)7月份查阅相关书籍和文献;(2)8月份完成开题报告并交老师批阅;(3)9月份完成论文初稿并交老师批阅;(4)10月份完成论文二搞并交老师批阅;(5)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 圣宝建. 关于数值积分若干问题的研究J. 南京信息工程大学. 2009.05.01. : 429 刘绪军. 几种求积公式计算精确度的比较J. 南京职业技术学院. 2009.10 史万明.吴裕树.孙新.数值分析M. 北京理工大学出版社.2010.4.指导教师意见： 签 名： 年 月 日 开 题 报 告 会 纪 要时间2013年8月26日地点宁静楼229教师办公室与会人员姓 名职务（职称）姓 名职务（职称）姓 名职务（职称）雍进军导师（讲师）邓喜才副教授李晟副教授龙林林组长会议记录摘要：指导小组针对课题二次函数性质的应用提问了以下问题以及报告人的回答：雍老师问:选择此题目的目的?潘晓祥答:随着计算机和计算方法的飞速发展,几乎所有学科都走向定量化和精确化,计算数学中的数值计算方法则是解决“计算”问题的桥梁和工具。邓老师问:对这个问题进行研究有什么实际的意义? 潘晓祥答:计算方法既有数学类课程中理论上的抽象性和严谨性,又有实用性和实验性的技术特征,计算方法是一门理论性和实践性都很强的学科.在科学研究和工程技术中都要用到各种计算方法.例如,在航天航空、地质勘探、汽车制造、桥梁设计、 天气预报和汉字字样设计中都有计算方法的踪影。李老师问:对这个问题你有什么自己的看法?潘晓祥答:随着计算机技术的迅速发展和普及,现在计算方法课程几乎已成为所有理工科学生的必修课程.我们知道,计算能力是计算工具和计算方法的效率的乘积,提高计算方法的效率与提高计算机硬件的效率同样重要.科学计算已用到科学技术和社会生活的各个领域中.所以,研究数值计算方法可以让数学的应用更大更广。会议主持人签名：记录人签名：年 月 日指导小组意见负责人签名： 年 月 日学 院 意 见负责人签名： 年 月 日贵 州 师 范 学 院 数学与计算机科学 学院指导教师指导本科毕业论文情况登记表论文（设计）题 目几种常用数值积分方法的比较用学生姓名潘晓祥学号1006020540200年级2010级所属学院数计学院专业数学与应用数学班级四班指导教师姓名雍进军职 称讲师学 历硕士指导时间指导地点指 导 内 容指导教师签名备 注2013年06月10日致远楼416论文选题，资料准备面授2013年06月22日网上确定毕业论文选题电子邮件2013年06月26日网上怎样撰写毕业论文开题报告电子邮件2013年06月28日网上指导学生撰写开题报告电子邮件2013年07月14日网上帮助学生查找有关参考文献电子邮件2013年07月17日手机如何构思自己的毕业论文手机飞信2013年08月21日手机听取学生毕业论文写作进展情况汇报手机飞信2013年08月28日网上解答学生在论文写作中遇到的疑惑电子邮件2013年9月09日网上帮助学生查找有关参考文献电子邮件2013年11月28日网上如何规划自己的论文电子邮件2013年12月04日手机怎样写好论文引言手机飞信2013年12月08日网上怎样写好论文引言电子邮件2013年12月12日网上怎样写论文摘要电子邮件2013年12月16日网上怎样选取论文关键词电子邮件2013年12月20日网上怎样编辑论文中的公式电子邮件2014年01月05日手机督促学生在寒假中写好论文的初稿电子邮件2014年02月27日宁静楼219检查学生论文完成情况面授2014年03月03日宁静楼219对学生的论文初稿提出修改时意见面授2014年03月07日宁静楼219解答学生在修改时的困惑面授2014年03月11日宁静楼219指导学生修改论文面授贵州师范学院数学与计算机科学学院本科毕业论文（设计）交叉评阅表学院（盖章）：学号1006020540200姓名潘晓祥专业数学与应用数学论文（设计）题目几种常用数值积分方法的比较班级四班指导教师意见评语： 该同学在论文撰写过程中对相关文献阅读范围广泛，方法正确，内容完整，能综合运用所学知识分析和解决实际问题。毕业论文撰写过程中态度端正，勤奋刻苦。论文研究了NewtonCotes求积公式、复化求积公式、Romberg积分、高斯积分方法，通过算例分析，得出几种常用数值积分方法是解决实际计算问题的重要方法。论文结构合理，符合逻辑，文章层次分明，语言准确，文字通顺，达到规范性要求，建议作为学士论文答辩。成绩： （满分100分） 指导教师（签名）： 年 月 日评阅教师意见评语：该同学具备较好的基础理论与专业知识，学习态度认真，阅读教师指定的参考资料、文献，较好的完成了任务书规定的工作量。论文研究了NewtonCotes求积公式、复化求积公式、Romberg积分、高斯积分方法，通过算例分析，得出几种常用数值积分方法是解决实际计算问题的重要方法。论文结构合理，符合逻辑，文章层次分明，语言准确，文字通顺，达到本科毕业论文相关要求。同意参加答辩。成绩： （满分100分） 评阅教师（签名）： 年 月 日贵州师范学院本科毕业论文答辩记录表论文题目几种常用的多项式插值方法作者姓名潘晓祥学号1006020540200年级2010级所属学院数计专业数学与应用数学班级本科（四）班指导教师姓名、职称雍进军 讲师答 辩 会 纪 要时间2014年 5 月11 日地点致远楼406答辩小组成员姓 名职务（职称）姓 名职务（职称）姓 名职务（职称）左羽教授崔忠伟副教授廖玉梅讲师答辩中提出的主要问题及回答的简要情况记录：1自己做的有哪些？答：第11页至第12页，总结进行比较。2程序运行过没有？答：运行过。320页程序代码中，if后的是什么符号？答：连接作用的符号。4解释一下什么时候用分号，什么时候不用？答：回答不清。5摘要中英文拼写有错答辩后修改答辩小组负责人签名：左羽记录人签名：梅林林2014年 5 月11 日 答辩小组意见评语：该生能在规定时间叙述论文的主要内容，对提出的问题一般能回答，无原则错误。答辩小组经过充分讨论，根据该生论文质量和答辩中的表现，同意评定论文成绩为“中等”。评定成绩：77 负责人（签名）：左羽 2014 年 5月11日欢迎下载目 录摘 要1Abstract21 前 言32 数值积分方法的基本思想33 几类常用数值积分方法的简单分析53.1 NewtonCotes求积公式53.2 复化求积公式633 Romberg求积公式734 高斯型求积公式84 几类数值积分方法的简单比较评述95 利用MATLAB编程应用对几类求积算法的分析比较10结束语13致 谢14附 录15欢迎下载摘 要 我们在求函数的积分时，往往因为原函数非常复杂以至于难以求出或用初等函数表示，这让我们计算起来非常困难，所以我们只能想办法求它的近似值，因此直接借助牛顿-莱布尼兹公式计算定积分的情况是非常少见的。这时候数值积分就是解决这种问题的一种很好很有效的方法。本文从数值积分问题的产生出发，详细介绍了一些数值积分的常用方法（NewtonCotes求积公式，复化求积公式，Romberg求积公式，高斯型求积公式）并对其进行了简要的分析，在探讨了这些数值积分算法的优缺点的理论之外，我们还将这些数值积分算法在计算机上通过matlab软件编程实现应用，并分别用各自求积公式进行运算，以此来分析比较各种求积公式的代数精度和计算误差。关键词 :数值积分;求积公式;代数精度Abstract function is very complex that it is difficult to find the elementary functions, which makes u We in the function for the integration, often because the original s very difficult to calculate, so we can only think of a way to find the approximate value, thus directly with Newton - Leibniz formula calculating definite integral situation is very rare. When numerical integration is to solve this problem in a very effective method. From the numerical integration problem, introduces some methods of numerical integration (Newton - Cotes quadrature formula, composite quadrature formulas, Longbei lattice quadrature formula, Gauss type quadrature formulas) and has carried on brief analysis, discusses the advantages and disadvantages of these numerical integration algorithm theory, we will these numerical integration algorithm in the computer by MATLAB software programming application, and separately with their respective quadrature formula for computing, in order to analyze the algebraic calculation precision and error comparison of various quadrature formulas.Keywords: Numerical integration; Calculationmeth; numerical analysis1 前 言微积分的发明是世界数学史上一项辉煌的成就。但在实际求积问题的时候，求解积分却有着非常多局限性。比如对于定积分在求某函数的定积分时，在一定条件下，虽然有牛顿-莱布里茨公式可以计算定积分的值，但在很多情况下的原函数不易求出或非常复杂。被积函数的原函数很难用初等函数表达出来，例如等；有的函数的原函数存在，但其表达式太复杂,计算量太大，有的甚至无法有解析表达式。因此能够借助牛顿-莱布尼兹公式计算定积分的情形是不多的。另外，许多实际问题中的被积函数往往是列表函数或其他形式的非连续函数对这类函数的定积分，也不能用不定积分方法求解，只能设法求其近似值。因此，探讨近似计算的数值积分方法是有明显的实际意义的即有必要研究定积分的数值计算方法，以解决定积分的近似计算。而数值积分就是解决此类问题的一种有效的方法，它的特点是利用被积函数在一些节点上的信息求出定积分的近似值。在很多实际应用中，只能知道积分函数在某些特定点的取值比如天气测量中的气温、湿度、气压等，医学测量中的血压、浓度等等。通过研究，我们将会更熟练掌握一些数值积分方法去计算一些特定条件的数值计算，以便我们得到自己想要的结果。2 数值积分方法的基本思想在数学分析中，计算连续函数f(x)在区间a,b上的积分是通过f(x)的原函数F(x),由下列定积分公式 得到的。但由于大量被积函数的原函数不能用初等函数表示，因此，很难用求原函数的公式得到积分；有些被积函数不是明显知道的，例如由数值表给出它的离散值，或者是它被定义为某个微分方程的解，而这个微分方程是不能显示解出的。这说明，按公式计算定积分是有很大局限性的。 因而常常采用在电子计算机上很有效的数值积分方法。我们从定积分的定义 出发。推导出两个简单的数值积分公式。式的几何意义，就是把整块曲线梯形的面积积分成若干个小曲边梯形面积的和，当无限细分时这个和取极限就是真正曲边梯形面积。去掉取极限这一步，用有限个小曲边梯形面积的和，代替整块的曲边梯形面积，从而求得一个近似值，这就是数值积分的基本思想。根据小区间的不同分割方法和各分点f()值的不同选择，就得到不同的数值积分公式。 数值求积公式是取上若干个点处的高度，通过加权后，再求和 从而得到积分的近似值。数值求积公式写成一般形式式中称求积节点，称求积系数，也称伴随节点的权。当积分区间确定后，求积系数仅仅与节点的选取有关，而不依赖被积函数的具体形式。记把称为求积公式的截断误差或余项。 数值求积方法的特点是直接利用积分区间上一些离散节点的函数值进行线性组合来近似计算定积分的值，从而将定积分的计算归结为函数值的计算，这就避开了牛顿-莱布尼兹公式需要寻求原函数的困难，并为计算机求积分提供了可行性。 3 几类常用数值积分方法的简单分析3.1 NewtonCotes求积公式 常用的梯形公式和Simpson公式是低阶的牛顿-柯特斯公式，牛顿-柯特斯公式是积分区间上等距节点的插值求积公式。插值求积公式在积分区间上，所取节点是等距时称为牛顿-柯特斯公式，即 其中为Cotes求积公式的系数，是n和k的函数。当n=1时，为梯形公式：梯形公式的代数精度为1，有两个积分节点。当n=2时，为Simpson公式：Simpson公式的代数精度为3，有三个积分节点。 由于只增加一个节点，其代数精度增加2，由此可知，Simpson公式比梯形公式代数精度高。当n=4时，NewtonCotes求积公式为Cotes公式：Newton-Cotes公式的代数精度为5，有5个积分节点。 所以对于Newton-Cotes积分公式，n为偶数时的代数精度要比n为奇数时的积分公式效果比较优越。但并不是n的值越大越好，当n过大时（n=8），求积公式的数值稳定性不好。3.2 复化求积公式 由于Newton-Cotes的节点n越大对应的精度就越高，但是n=8时公式的数值是不稳定的，因此就不能用增加求积节点的方法来提高精度，因此，我们常常将求积区间a,b分成若干小区间，然后在每个小区间上采用数值稳定的Cotes公式求小区间上的积分，然后把每个小区间上的结果加起来作为原定积分的近视值，这种方法构造的求积公式就叫做复化求积公式。常用的复化求积公式有：复化梯形公式： 变步长梯形公式为：复化Simpson公式：变步长复化Simpson公式：3.3 Romberg求积公式 Romberg积分方法也叫做逐次分半加速法，它是在复化梯形公式误差估计的基础上，应用线性外推的方法构造出的一种加速算法。 将积分区间分成n等分和2n等分时，求得积分近似值和，并没有误差估计式 积分近似值的误差大致等于，当用对进行修正时，与之和比更接近于真值,故是对误差的一种补偿，因此可以期望下式是一个更好的结果，即 下面说明即是分成n等分时Simpson公式的值。将复化梯形公式梯形变步长求积公式代入上式表达式得 这就是说，用梯形法二分前后两个梯形值和作线性外推，结果得到Simpson法的积分值。将误差由变为，从而提高了逼近精度。 再考察Simpson法。其截断误差与成正比，因此，若将步长折半，则误差减至，即有 由此得 不难验证，上式右端的值其实等于，就是说，用Simpson法二分前后的两个积分值与,按上式再作线性外推，结果得到柯特斯法的积分值，即有这时将误差由变为，逼近精度又一次得以提高。 同样的方法，依据柯特斯法的误差公式，可进一步导出下列龙贝格公式 逼近积分值的误差为，这样Romberg公式将误差由变为，逼近精度再次得以提高。Romberg公式有7次代数精度，这表明该公式不是牛顿-柯特斯公式。 在步长二分的过程中运用、表达式加工三次，就能将粗糙的积分值逐步加工成精度较高的Romberg值，或者说，将收敛缓慢的梯形值序列加工成收敛迅速的Romberg值序列，这种加速方法称Romberg算法。3.4 高斯型求积公式 前面介绍的个节点的 Newton -Cotes求积公式，其特征是节点是等距的。这种特点使得求积公式便于构造，复化求积公式易于形成。但同时也限制了公式的精度。是偶数时，代数精度为，是奇数时，代数精度为；我们知道个节点的插值型求积公式的代数精确度不低于。能不能在区间上适当选择个节点 使插值求积公式的代数精度高于呢？ 答案是肯定的，适当选择节点，可使公式的精度最高达到，这就是所学的高斯型求积公式。 不失一般性，将求积公式的求积区间转换成的形式。 对任意求积区间作变换 可以变换到区间上，这时其中。 高斯-勒让德求积公式在这里简称高斯公式，它是在区间上进行讨论的。4 几类数值积分方法的简单比较评述 由于我们在计算实际问题是往往要考虑到代数精度和计算量，所以不同类型的求积公式有着不同的特点： Simpson积分方法和梯形积分方法虽然计算简便，但是精度比较差，不理想。但对于光滑性较差的被积函数有时会比高精度的积分方法更为有效。特别是梯形积分方法对被积函数是周期函数的求积效果更为突出。n7时，NewtonCotes公式是不稳定的，然而复化梯形公式和复化Simpson公式不仅保留了低阶公式的优点还能够获得比较较高的精度，所以在实际计算中应用得最为广泛。 Romberg积分方法的算法简单，方便编程的实现。收敛速度快、计算精度较高，但是计算量较大。 Gauss积分方法的精度较高，数值稳定、收敛速度较快，但因为其节点不规则，计算比较麻烦。5 利用MATLAB编程应用对几类求积算法的分析比较 在简单的认识积分方法比且理论比较之后，则要进行数学实验进行验证，因此就要通过matlab软件对各种积分方法进行编程并运算，然后对其各种方法的运算结果进行分析比较，掌握和理解各方法的优缺点。 规定各个程序都以为例子进行运算。原积分的精确值为例 分别用不同的方法计算积分，并作比较。 用以上介绍的几类积分方法分别计算积分，得出误差，并进行比较：1、用Newton-Cotes公式当n=1时，即用梯形公式，用程序一 （程序见附录） 在MATLAB命令窗口中输入 NCotes(0,1,1,2) 得0.92703549240395 0.01904757796323当n=2时, 即用Simpson公式,用程序一 （程序见附录） 在MATLAB命令窗口中输入 NCotes(0,1,2,2)得0.94614588227359 0.000062811906407当n=4时, 即用科特斯公式，用程序一 （程序见附录） 在MATLAB命令窗口中输入 NCotes(0,1,4,2)得0.94608300406367 0.0000000663035132、用复化梯形公式 令h=1/8=0.125，用程序二 （程序见附录） 在MATLAB命令窗口中输入 trapr1(f,0,1,8),得=0.94569086352700.0003922068401823、用复化 Simpson 公式 令h=1/8=0.125，用程序三 （程序见附录） 在MATLAB命令窗口中输入 simpr1(f,0,1,8),得 =0.94608308538495 0.000000015017767 4、用Romberg公式 用程序四 （程序见附录） 在MATLAB命令窗口中输入 romber(f,0,1,5,0.5*(10(-8)，得0.94608307036718 0.0000000000000025、用高斯-勒让德求积公式 令，（1） 用2个节点的Gauss公式0.94604115827633（2） 用3个节点的Gauss公式，用程序五 （程序见附录） 在MATLAB命令窗口中输入 GuassLegendre (0,1,2,2)，得0.946083134078473 0.000000063711290算法比较:1.原积分的精确值为：2.由例题的各种求积算法可知：（1）对Newton-cotes公式，当n=1时只有1位有效数字，当n=2时有3位有效数字，当n=4时有7位有效数字。（4）用复化梯形公式有2位有效数字，对复化Simpson公式有7位有效数字。（5）用复化梯形公式，对积分区间二分11次用了2049个函数值，才可以得7位有效数字。（6）用Romberg公式对区间二分3次用了9个函数值，就可以得到7位有效数字；二分4次用了14个函数值，却可以得到14位有效数字。（7）用高斯-勒让德求积公式仅仅用了3个函数值，就能得到比较精确的6位有效数字。结束语 本文主要研究了常用的几类数值积分的求积算法并通过例题计算积分进行分析比较。Newton-Cotes积分方法是一种非常普遍的积分方法，然而梯形积分方法的误差最大，近似效果最差，Simpson积分方法的精度比梯形积分方法高了一个数量级；Cotes积分方法精度比Simpson积分方法高两个数量级。则Cotes代数精度比较高。由此可知一般情况下，积分公式代数精度越高，计算精度也越高。但是高阶的Cotes积分方法收敛性没有保证，因此实际应用中很少用。复化梯形积分方法比梯形积分方法精度高，同样的，复化Simpson积分方法比Simpson积分方法精度高，高了差不多7个数量级，所以复化积分方法比较优越。Romberg积分方法收敛速度快、计算精度较高，但是计算量较大。Gauss积分方法精度高、数值稳定、收敛速度较快，但是计算麻烦。经研究可以知道Newton-Cotes方法的代数精度越高，数值积分的效果越好、越精确。当积分区间比较大的时候，积分数值不稳定，这个时候可以利用复化积分方法效果会更好；Romberg积分方法可以利用变步长复化积分公式得到更为精确的数值结果，是比较好的积分方法。高斯求积方法精确度高，收敛性快，比其他积分方法优越。具有很广泛的运用。参考文献1 关治. 陆金甫. 数学分析基础（第二版） M. 北京:等教育出版社.2010.72 胡祖炽. 林源渠. 数值分析 M 北京:等教育出版社.1986.33 薛毅. 数学分析与实验 M 北京:业大学出版社 2005.34 徐士良. 数值分析与算法 M. 北京:械工业出版社2007.15 王开荣. 杨大地. 应用数值分析 M 北京:等教育出版社 2010.76 杨一都. 数值计算方法M. 北京:等教育出版社 . 2008.47 韩明. 王家宝. 李林. 数学实验（MATLAB）版M. 上海:济大学出版社 2012.18 圣宝建. 关于数值积分若干问题的研究J. 南京信息工程大学. 2009.05.01. : 42 9 刘绪军. 几种求积公式计算精确度的比较J. 南京职业技术学院.2009.10 史万明.吴裕树.孙新.数值分析M. 北京理工大学出版社.2010.4.致 谢 行文至此，我的这篇论文已接近尾声；岁月如梭，我四年的大学时光也即将敲响结束的钟声。离别在即，站在人生的又一个转折点上，心中难免思绪万千，一种感恩之情油然而生。首先感谢贵州师范学院四年来对我的培养，是博学的老师们教会了我学习的方法、锻炼了我思考的能力、指明了我未来奋斗的方向，从而使我进一步明确了人生的目标。其次，我要感谢我的指导老师雍进军老师，他的严谨细致、一丝不苟的作风一直是我工作、学习中的榜样；他的循循善诱的教导和不拘一格的思路给予我无尽的启迪。在撰写整个毕业论文的过程当中，他为我们考虑到了每一个细节，从开题报告到毕业论文的拟定修改上，雍老师更是不厌其烦的为我们做好每一步的细心指导。对此，我表示衷心地感谢。没有雍老师，我的论文也不可能这么顺利的完成。同时，我也要感谢每一位给过我帮助的老师和同学，在我撰写论文的过程当中同样给了我大量有益的建议，在此一并向他们表示真诚的感谢，感谢他们对我的支持和帮助。最后感谢这篇论文所涉及到的各位学者，本文引用了数位学者的研究文献，如果没有各位学者的研究成果带给我的的帮助和启发，我将很难完成本篇论文的写作。由于我的学术水平有限，所写论文难免有不足之处，恳请各位老师和学友批评指正。最后，衷心感谢评阅论文及参加答辩的各位老师！附 录1 NewtonCotes求积公式的MATLAB实现 先用M文件定义一个名为f1.m的函数：% i是要调用第几个被积函数g(i)，x是自变量function f=f1(i,x) g(1)=sqrt(x);if x=0 g(2)=1;elseg(2)=sin(x)/x;endg(3)=4/(1+x2);f=g(i);程序一：function C,g=NCotes(a,b,n,m)% a，b分别为积分的上下限；% n是子区间的个数；% m是调用上面第几个被积函数；% 当n=1时计算梯形公式；当n=2时计算辛浦生公式，以此类推； i=n; h=(b-a)/i; z=0;for j=0:i x(j+1)=a+j*h; s=1; if j=0 s=s; elsefor k=1:j s=s*k;endendr=1;if i-j=0 r=r;elsefor k=1:(i-j) r=r*k;endendif mod(i-j),2)=1 q=-(i*s*r);else q=i*s*r;endy=1;for k=0:i if k=j y=y*(sym(t)-k); endendl=int(y,0,i);C(j+1)=l/q; z=z+C(j+1)*f1(m,x(j+1);endg=(b-a)*z1）当输入，时，即在MATLAB命令窗口输入 NCotes(0,1,1,2)即可得用梯形公式的积分值和相应科特斯系数如图1.12）当输入，时，即在MATLAB命令窗口输入 NCotes(0,1,2,2)即可得用辛浦生公式的积分值和相应科特斯系数如图1.23）当输入，时，即在MATLAB命令窗口输入 NCotes(0,1,4,2)即可得用科特斯公式的积分值和相应科特斯系数 如图1.3 图 1.1图 1.2图1.32 复化梯形求积公式的MATLAB实现通过的个等步长节点逼近积分其中,程序二：function s=trapr1(f,a,b,n)% f是被积函数；% a,b分别为积分的上下限；% n是子区间的个数；% s是梯形总面积；h=(b-a)/n;s=0;for k=1:(n-1) x=a+h*k; s=s+feval(f,x);endformat long s=h*(feval(f,a)+feval(f,b)/2+h*s;先用M文件定义一个名为f.m的函数：function y=f(x)if