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

1、SHANGHAI JIAO TONG UNIVERSITY题目名称：复合梯形公式与复合辛普森公式对比学生：_ 学生学号: 班 级:学院（系）:目录1. 概述3.2. 问题提出4.3. 算法推导6.4. 算法框图&4.1复合梯形公式算法流程图84.2复合辛普森公式算法流程图85. MATLAB 源程序106. 结论与展望12图表目录图4-1复合梯形公式算法流程图 8.图4-2复合辛普森公式算法流程图 1.0图6-1 MATLAB计算结果 121.概述梯形求积公式和辛普森求积公式分别是牛顿 -科斯特公式中n=1和n=2时的情形。其中梯形求积公式可表示为bb aaf(X)dX 丁(住)f(b)其公式

2、左端是以a,b区间上积分，右端为b-a为高、端点函数值为上下底的 梯形的面积值，故通称为梯形公式，具有 1次代数精确度。类似的，辛普森求积公式可以表示为baa bS 害f(a) 4f (畔)f(b)6 2该公式一般在立体几何中用来求 拟柱体的体积，由于偶数n阶牛顿-科特斯求 积公式至少具有n+1次代数精确度，所以辛普森公式实际上具有 3次代数精确 度。由于牛顿-科斯特公式在n8时不具有稳定性，故不可能通过提高阶的方法 来提高求积精度。为了提高精度通常可把积分区间分成若干子区间(通常是等 分),再在每个子区间上用低阶求积公式。这种方法称为复合求积法。本文主要讨论复合梯形公式和复合辛普森公式在同一

3、数学问题中的应用。首先给出了复合梯形公式和复合辛普森公式的推导过程以及其余项的表达形式，然后用流程图的形式介绍算法思路，再运用 MATLAB编写代码计算结果，最后对 结果进行对比讨论。希望通过两个算法在同一个算例中的应用对比，更好的理解和掌握复合梯形公式和复合辛普森公式的适用围和适用条件。并且能够熟悉MATLAB编程求解问题的流程，掌握编程化的思想方法。 同时对两种方法的计算结果对比分析，讨论两种求积方法的计算精度。2问题提出对于函数给出的函数表如下，试用复合梯形公式和复合辛普森公式计算积分0表2-1函数计算结果表xf(x)011/80.18221/40.80923/80.80921/20.8

4、4065/80.47403/40.11127/80.403110.7897ba f(x)dxb a2 (f(a) f(b)3.1复合梯形公式根据梯形公式,3.算法推导将区间厲工划分为n等份，分点呛二=乳，在每个子区间l+1|(fc = O.lTi-l)上采用梯形公式，则得：f(x)dxf(x)dxk 0(f(Xk)0f(Xki) Rn(f )h n 1Tn - f(Xk)2 k of(xk1) -f(a)n 12f(xQk 1f(b)则为复合梯形公式另外，复合梯形公式的余项可表示为Rn(f )if ()3.2复合辛普森公式根据辛普森公式S /戸心）4f（a b）f（b）6 2将区间厲”划分为n

5、等份，在每个子区间上采用辛普 森公式。若记则得f(x)dxaxk 1f(x)dxhn6kf(Xk)4f(Xki/2) f(Xki) Rn(f )Sn -f(a)n 14 f (xk 1/2 )k 0n 12f(xQk 1f(b)该公式即为复合辛普森公式。复合辛普森公式的余项可表示为Rn(f)需(”()4. 算法框图4.1复合梯形公式算法流程图开始求出步长h，各节点耳及相应的各节点函数值求和 sum.输出积分值结束图4-1复合梯形公式算法流程图4.2复合辛普森公式算法流程图开始求出步长h，各节点相邻节点中点 池+1/二及相应的函数值卜汽，弋：Dk=1,2,.n-1各节点函数值求和suml,各相邻

6、节点中心点函数值（知+ V2）求和sum2+ fW + 2sumA + 45m2)输出积分值卜结束图4-2复合辛普森公式算法流程图5. MATLAB程序%K合梯形公式及复合辛普森积分公式clear all;format Io ng;a=0;b=1;n=8;h=(b-a)/n;% 步长for i=1: n+1x(i)=a+(i-1)*h;if isnan (si n(x(i)/x(i)syms t;tmp=limit(si n( t)./t,t,x(i);%y(i)=eval(tmp);elsey(i)=si n(x(i)/x(i);%endend当被积函数在某点值不存在时，求其极限被积函数求节

7、点的值%复合梯形公式及复合辛普森积分公式s1=0;for k=2:ns仁 s1+y(k);endT8=h/2*(y(1)+2*s1+y( n+1)%复合辛普森积分公式s2=0;s3=0;for k=2:2:ns2=s2+y(k);endfor k=3:2: n-1s3=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. 结论与展望 aZOlDllirfuheTixinTLfuhexinpusengongshi t弟形和分公式1 0.

8、 945691辛昔森公式积分：0. 946083图6-1 MATLAB计算结果运行MATLABS序，得到复合梯形求积公式的积分值为 0.945691，复合辛普 森求积公式的积分值为0.946083 （四舍五入后保留6位小数）。而实际的积分准 确值保留到6位小数的结果为0.946083。通过上述结果对比可以得出，虽然复合梯形公式将区间分成了 8等分而复合 辛普森公式将区间分成了 4等分，但两种计算方法实际都需要使用 9个点上的函 数值，计算量基本也相同，然而最终精度差别却很大。在保留6位小数的前提下， 复合辛普森法计算结果与精确解完全一致，而复合梯形公式的计算结果却只有前 两位数字与精确解相同，

9、误差相对比较大。F面利用余项公式来估计两种算法的误差。首先需要求的高阶导f(x)sin x1cos(xt)dt,所以有kfk(x)i dki kkk(cosxt)dtt cos(xt )dt,dx02于是1/k 、fk(x)0cos(xt )2max0 x 1tkdt1tkdt0从而复合梯形公式的误差max f (x)12 0 x 1E(8)230.43410 3而复合辛普森公式的误差11R8(f) 2880(245 0271106.从而，对比两者可得，复合辛普森公式在计算该问题时的精度远高于复合梯形公 式。通过以上分析，本文所得结论如下：1. 复合梯形公式和复合辛普森公式都可以用来作为数值积分估算的替代公式。2. 在计算量基本相同的前提下，复合辛普森公式计算结果的计算精度要比复合 梯形公式计算精度高的多。3. 本算例也验证了辛普森公式作为偶数阶牛顿-柯特斯公式的更为精确的代数 精度。关于如何开展下一步研究，提出以下构想：1. 对多个算例进行分析，保证计算量基本相同的情况下去比较计算精度，验证复合辛普森公式具有更高精度的结论。2. 对多个算例进行MATLAB编程分析，在要求相同计算精度的前提下去比较计算量的大小，从而分析复合梯形公式与复合辛普森公式的优劣参考文献1 穆耶赛尔艾合买提，阿布都热西提阿布都外力改进复合梯形求积公式J. 首都师大学学报（自然科学版）,2016