数值积分算法与MATLAB实现陈悦5133201.doc

东北大学秦皇岛分校数值计算课程设计报告数值积分算法及MATLAB实现学 院数学与统计学院专 业信息与计算科学学 号5133201姓 名陈悦指导教师姜玉山 张建波成 绩教师评语：指导教师签字： 2015年07月14日 数学与统计学院课程设计（实习）报告 第 10 页1 绪论数值分析是计算数学的一个主要部分，计算数学是数学科学的一个分支，它研究用计算机求解各种数学问题的数值检索方其理论与软件的实现.而数值分析主要研究数值计算.现科学技术的发展与进步提出了越来越多的复杂的数值计算问题，这些问题的圆满解决已远人工手算所能胜任，必须依靠电子计算机快速准确的数据处理能力.这种用计算机处理数值问题的方法，成为科学计算.今天，科学计算的应用范围非常广泛，天气预报、工程设计、流体计算、经济规划和预测以及国防尖端的一些科研项目，如核武器的研制、导弹和火箭的发射等，始终是科学计算最为活跃的领域.1.1 数值积分介绍数值积分是数值分析的重要环节，实际问题当中常常需要计算积分，有些数值方法，如微分方程和积分方程的求解，也都和积分计算相联系.求某函数的定积分时，在多数情况下，被积函数的原函数很难用初等函数表达出来，因此能够借助微积分学的牛顿-莱布尼兹公式计算定积分的机会是不多的.另外，许多实际问题中的被积函数往往是列表函数或其他形式的非连续函数，对这类函数的定积分，也不能用不定积分方法求解.由于以上原因，数值积分的理论与方法一直是计算数学研究的基本课题.对微积分学做出杰出贡献的数学大师，如I.牛顿、L.欧拉、C.F.高斯、拉格朗日等人都在数值积分这个领域作出了各自的贡献，并奠定了这个分支的理论基础.构造数值积分公式最通常的方法是用积分区间上的n 次插值多项式代替被积函数，由此导出的求积公式称为插值型求积公式.特别在节点分布等距的情形称为牛顿-科特斯公式，例如梯形公式(Trapezoidal Approximations)与抛物线公式(Approximations Using Parabolas)就是最基本的近似公式.但它们的精度较差.龙贝格算法是在区间逐次分半过程中，对梯形公式的近似值进行加权平均获得准确程度较高的积分近似值的一种方法，它具有公式简练、计算结果准确、使用方便、稳定性好等优点，因此在等距情形宜采用龙贝格求积公式(Rhomberg Integration).当用不等距节点进行计算时，常用高斯型求积公式计算，它在节点数目相同情况下，准确程度较高，稳定性好，而且还可以计算无穷积分.数值积分还是微分方程数值解法的重要依据.许多重要公式都可以用数值积分方程导出.现探讨数值积分算法以及运用MATLAB软件的具体实现1.2 MATLAB软件MATLAB是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分.MATLAB是matrix&laboratory两个词的组合，意为矩阵工厂（矩阵实验室）.是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境.它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平.MATLAB和Mathematica、Maple并称为三大数学软件.它在数学类科技应用软件中在数值计算方面首屈一指.MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域.MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且MATLAB也吸收了像Maple等软件的优点，使MATLAB成为一个强大的数学软件.在新的版本中也加入了对C，FORTRAN，C+，JAVA的支持.2 数值积分的基本概念 一般的，我们可以在区间上适当选取某些节点，然后用的加权平均得到平均高度的近似值，这样构造出的求积公式具有下列形式：，式中称为求积节点；称为求积系数，亦称伴随节点的权.权仅仅与节点的选取有关，而不依赖于被积函数的具体形式.2.1 代数精度的概念如果某个求积公式对于次数不超过的多项式均能准确的成立，但对于次多项式就不准确成立，则称该求积公式具有次代数精度（或代数精确值）一般地，欲使求积公式具有次代数精度，只要令它对于都能准确成立，这就要求:2.2 求积公式的余项令求积公式的余项为，其中.；区间可以是有限的或无限的.构造求积公式的问题就是确定和使得在某种意义下尽可能地小.3 数值积分方法及MATLAB实现3.1 复合辛普森公式3.1.1 插值型求积公式用插值多项式替换积分中的被积函数，然后计算作为积分的近似值，这样建立的求积公式称为插值型求积公式.用插值多项式的表达式，代入得，其中: .其余项为:.3.1.2 牛顿-科特斯公式介绍 取等距节点，把积分区间剖分成等分.令步长，并记，则个节点为，代入得:=.这种等距节点的插值型求积公式通常称为牛顿-科特斯公式. 3.1.3 辛普森公式利用牛顿-科特斯公式，取=2，此时为，即为辛普森公式，其余项为.3.1.4 复合辛普森公式将积分区间分成等分，分点为，其中，记区间的中点为，在每个小区间上用辛普森公式，则得到所谓的复合辛普森公式：.余项为，.3.1.5 复合辛普森公式的MATLAB实现代码如下：function s=xinpusen(fun,a,b,n)h=(b-a)/ns1=0;s2=0;for k=0:(n-1) x=a+h*k; s1=s1+feval(fun,x);endfor k=0:(n-1) x=a+h*(k+1/2); s2=s2+feval(fun,x);ends=h/6*(feval(fun,a)+feval(fun,b)+2*s1+4*s2);3.2 龙贝格公式3.2.1 梯形法的递推化将区间分成等份，共有个分点，如果将求积区间再二分一次，则分点增至个，用复合梯形公式求得该子区间上的积分值为，将每个子区间上的积分值相加得=.得到递推公式：.3.2.2 龙贝格算法公式当在上充分光滑时， 可证用逼近的截断误差是：=按理查森外推法：，其中，为满足的适当正数 .取序列：.用来逼近的误差为，这种算法就是龙贝格算法.3.2.3 龙贝格算法MATLAB实现代码如下：function s=longbeige(fun,a,b,tol)if nargintol) i=i+1;h=h/2; T(i+1,1)=T(I,1)/2+sum(feval(fun,a+h/2:h:b-h/2)*h/2; for j=1:i T(i+1,j+1)=(4j*T(i+1,j)-T(I,j)/(4j-1); endendTs=T(i+1,j+1);3.3 自适应法自适应积分法是一种比较经济而且快速的求积分的方法.他能自动地在被积函数变化剧烈的区域增多节点，而在被积函数变化平缓的地方减少节点.因此它是一种不均匀区间的积分方法.按照子区间上的积分方式它可以分为自适应辛普森积分法和自适应梯形积分法.通常是采用自适应辛普森积分法作为子区间的积分方式. 自适应积分法的基本步骤如下：(1) 将积分区间分成两个相等的1级子区间和，且；(2) 在上述两个1级子区间上用辛普森积分得到积分和;(3) 将子区间分成两个相等的2级子区间和;(4) ;(5) 比较和， 如果| - | fun=inline(4./(1+x.2); xinpusen(fun,0,1,10)ans = 3.2749259863031184.2 龙贝格算法求解longbeige(inline(4./(1+x.2),0,1,1e-6)T = 3.0000 3.1000 3.1333 3.1312 3.1416 3.1421 3.1390 3.1416 3.1416 3.1416 3.1409 3.1416 3.1416 3.1416 3.1416 3.1414 3.1416 3.1416 3.1416 3.1416 3.1416ans =3.1415926536382444.3 高斯算法求解 gaosi(inline(4./(1+x.2),0,1,2,3)ans = 3.141591222382834 gaosi(inline(4./(1+x.2),0,1,4,4)ans = 3.1415956115587354.4 三种方法比较分析结果显示每一个算法都接近真实值，但龙贝格算法相比较复合辛普森算法，高斯算法来说更加接近.对于代数精度来说，复合辛普森的代数精度为11，龙贝格代数精度为11，高斯代数精度为11.可见代数精度相同时，龙贝格的求积精度最小，所以相同条件下龙贝格求积公式最能接近准确值.总结随着数学实验的兴起， 对整个数学课程教学改革起到了积极的推动作用， 我们要熟悉的运用各种数学软件，解决数学运算中繁琐的问题，实现学习的简单，快捷化.同时意识到用MATLAB编程时，要实现代码的层次性，做到有规有矩，那样才能把MATLAB运用自如.这次课程设计，用MATLAB实验对数值积分进行了实现，简介了5种不同的数值积分的方法，并且实现了其中的3中方法，实现过程中发现了各种方法之间的区别和联系.并且在实验过程中，使自己对数值积分和MATLAB更加的熟悉.做到了学习和实践相联系.参考文献1 戈慈水.数值分析课程