复合梯形公式与复合辛普森公式对比

wordSHANGHAI JIAO TONG UNIVERSITY题目名称：复合梯形公式与复合辛普森公式比照学生某某: 学生学号: 班 级: 学院(系): 目录1.概述32.问题提出43.算法推导54.算法框图64.1复合梯形公式算法流程图64.2 复合辛普森公式算法流程图75.MATLAB源程序86.结论与展望9图表目录图 41 复合梯形公式算法流程图6图 42 复合辛普森公式算法流程图7图 61 MATLAB计算结果9表 21函数计算结果表41. 概述梯形求积公式和辛普森求积公式分别是牛顿-科斯特公式中n=1和n=2时的情形。其中梯形求积公式可表示为其公式左端是以a,b区间上积分，右端为b-a为高、端点函数值为上下底的梯形的面积值，故通称为梯形公式，具有1次代数准确度。类似的，辛普森求积公式可以表示为该公式一般在立体几何中用来求拟柱体的体积，由于偶数n阶牛顿-科特斯求积公式至少具有n+1次代数准确度，所以辛普森公式实际上具有3次代数准确度。由于牛顿-科斯特公式在n8时不具有稳定性，故不可能通过提高阶的方法来提高求积精度。为了提高精度通常可把积分区间分成假如干子区间通常是等分，再在每个子区间上用低阶求积公式。这种方法称为复合求积法。本文主要讨论复合梯形公式和复合辛普森公式在同一数学问题中的应用。首先给出了复合梯形公式和复合辛普森公式的推导过程以与其余项的表达形式，然后用流程图的形式介绍算法思路，再运用MATLAB编写代码计算结果，最后对结果进展比照讨论。希望通过两个算法在同一个算例中的应用比照，更好的理解和掌握复合梯形公式和复合辛普森公式的适用X围和适用条件。并且能够熟悉MATLAB编程求解问题的流程，掌握编程化的思想方法。同时对两种方法的计算结果比照分析，讨论两种求积方法的计算精度。2. 问题提出对于函数fx=sinxx给出的函数表如下，试用复合梯形公式和复合辛普森公式计算积分I=01sinxxdx。表 21函数计算结果表xf(x)011/80.9973978670818221/40.9896158370180923/80.9896158370180921/20.9588510772084065/80.9361556367047403/40.9088516800311127/80.87719257398403110.8414709848078973. 算法推导3.1复合梯形公式根据梯形公式，将区间a,b划分为n等份，分点xk=a+kh,h=b-an,k=0,1,n，在每个子区间xk,xk+1(k=0,1,n-1)上采用梯形公式，如此得：记如此Tn为复合梯形公式。另外，复合梯形公式的余项可表示为3.2 复合辛普森公式根据辛普森公式将区间a,b划分为n等份，在每个子区间xk,xk+1(k=0,1,n-1)上采用辛普森公式。假如记xk+1/2=xk+12h如此得记该公式即为复合辛普森公式。复合辛普森公式的余项可表示为4. 算法框图开始4.1复合梯形公式算法流程图输入区连续点a,b与等分数n求出步长h，各节点xk与相应的函数值f(xk),k=1,2,.n-1各节点函数值f(xk)求和sum, k=1,2,.n-1Tn=h2(fa+fb+sum)输出积分值Tn完毕图 41 复合梯形公式算法流程图4.2 复合辛普森公式算法流程图开始输入区连续点a, b与等分数n求出步长h，各节点xk，相邻节点中点xk+1/2与相应的函数值f(xk), f(xk+1/2)k=1,2,.n-1各节点函数值f(xk)求和sum1,各相邻节点中心点函数值f(xk+1/2)求和sum2 k=1,2,.n-1Sn=h6(fa+fb+2sum1+4sum2)输出积分值Sn完毕图 42 复合辛普森公式算法流程图5. MATLAB源程序%复合梯形公式与复合辛普森积分公式clear all;format long;a=0;b=1;n=8;h=(b-a)/n;%步长for i=1:n+1 x(i)=a+(i-1)*h; if isnan(sin(x(i)/x(i) syms t; tmp=limit(sin(t)./t,t,x(i);%当被积函数在某点值不存在时，求其极限 y(i)=eval(tmp); else y(i)=sin(x(i)/x(i);%被积函数求节点的值 end end%复合梯形公式与复合辛普森积分公式s1=0;for k=2:n s1=s1+y(k);endT8=h/2*(y(1)+2*s1+y(n+1)%复合辛普森积分公式s2=0;s3=0;for k=2:2:n s2=s2+y(k);endfor k=3:2:n-1 s3=s3+y(k);endh1=2*h;%注：此时步长是原来的2倍S4=h1/6*(y(1)+4*s2+2*s3+y(n+1)fprintf(梯形积分公式：%6.6fn辛普森公式积分：%6.6fn,T8,S4)6. 结论与展望图 61 MATLAB计算结果运行MATLAB程序，得到复合梯形求积公式的积分值为0.945691，复合辛普森求积公式的积分值为0.946083四舍五入后保存6位小数。而实际的积分准确值保存到6位小数的结果为0.946083。通过上述结果比照可以得出，虽然复合梯形公式将区间分成了8等分而复合辛普森公式将区间分成了4等分，但两种计算方法实际都需要使用9个点上的函数值，计算量根本也一样，然而最终精度差异却很大。在保存6位小数的前提下，复合辛普森法计算结果与准确解完全一致，而复合梯形公式的计算结果却只有前两位数字与准确解一样，误差相比照拟大。下面利用余项公式来估计两种算法的误差。首先需要求fx=sinxx的高阶导数。由于,所以有,于是.从而复合梯形公式的误差.而复合辛普森公式的误差.从而，比照两者可得，复合辛普森公式在计算该问题时的精度远高于复合梯形公式。通过以上分析，本文所得结论如下：1. 复合梯形公式和复合辛普森公式都可以用来作为数值积分估算的替代公式。2. 在计算量根本一样的前提下，复合辛普森公式计算结果的计算精度要比复合梯形公式计算精度高的多。3. 本算例也验证了辛普森公式作为偶数阶牛顿-柯特斯公式的更为准确的代数精度。关于如何开展下一步研究，提出以下构想：1. 对多个算例进展分析，保证计算量根本一样的情况下去比拟计算精度，验证复合辛普森公式具有更高精度的结论。2. 对多个算例进展MATLAB编程分析，在要求一样计算精度的前提下去比拟计算量的大小，从而分析复合梯形公式与复合辛普森公式的优劣。参考文献1穆耶赛尔艾合买提,阿布都热西提阿布都外力. 改良复合梯形求积公式J. 首都师X大学学报(自然科学版),2016,06:1-4.2X智颖,王向公,任威龙,龙耀萍. 辛普森公式的推广形式与应用J. 某某理工大学学报(自然科学版),2014,01:29-31.3陈凤忠. 辛普森公式的统一美J. 某某理工大学学报(综合版),2005,03:59-60+80.4李庆扬, 王能超, 易大义. 数值分析.第5版M. 清华大学, 2008,12:106-109.14 / 14