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

复合梯形公式与复合辛普森公式对比复合梯形公式与复合辛普森公式对比复合梯形公式与复合辛普森公式对比复合梯形公式与复合辛普森公式对比编制仅供参考审核批准生效日期地址：电话：传真:邮编:SHANGHAIJIAOTONGUNIVERSITY题目名称：复合梯形公式与复合辛普森公式对比学生姓名:学生学号:班级:学院(系):目录1. 概述 32. 问题提出 43. 算法推导 54. 算法框图 6复合梯形公式算法流程图 6复合辛普森公式算法流程图 75. MATLAB源程序 86. 结论与展望 9图表目录TOC\h\z\c"图"图41复合梯形公式算法流程图 6图42复合辛普森公式算法流程图 7图61MATLAB计算结果 9TOC\h\z\c"表"表21函数计算结果表 4概述梯形求积公式和辛普森求积公式分别是牛顿-科斯特公式中n=1和n=2时的情形。其中梯形求积公式可表示为其公式左端是以[a,b]区间上积分，右端为b-a为高、端点函数值为上下底的梯形的面积值，故通称为梯形公式，具有1次代数精确度。类似的，辛普森求积公式可以表示为该公式一般在立体几何中用来求拟柱体的体积，由于偶数n阶牛顿-科特斯求积公式至少具有n+1次代数精确度，所以辛普森公式实际上具有3次代数精确度。由于牛顿-科斯特公式在n≥8时不具有稳定性，故不可能通过提高阶的方法来提高求积精度。为了提高精度通常可把积分区间分成若干子区间（通常是等分），再在每个子区间上用低阶求积公式。这种方法称为复合求积法。本文主要讨论复合梯形公式和复合辛普森公式在同一数学问题中的应用。首先给出了复合梯形公式和复合辛普森公式的推导过程以及其余项的表达形式，然后用流程图的形式介绍算法思路，再运用MATLAB编写代码计算结果，最后对结果进行对比讨论。希望通过两个算法在同一个算例中的应用对比，更好的理解和掌握复合梯形公式和复合辛普森公式的适用范围和适用条件。并且能够熟悉MATLAB编程求解问题的流程，掌握编程化的思想方法。同时对两种方法的计算结果对比分析，讨论两种求积方法的计算精度。问题提出对于函数fx=sinx表STYLEREF1\s2SEQ表\*ARABIC\s11函数计算结果表xf(x)011/81/43/81/25/83/47/81算法推导复合梯形公式根据梯形公式，将区间[a,b]划分为n等份，分点xk=a+kh,h=b-a记则Tn另外，复合梯形公式的余项可表示为复合辛普森公式根据辛普森公式将区间[a,b]划分为n等份，在每个子区间若记x则得记该公式即为复合辛普森公式。复合辛普森公式的余项可表示为算法框图开始复合梯形公式算法流程图开始输入输入区间断点a,b及等分数n求出步长h求出步长h，各节点xk及相应的函数值f(各节点函数值各节点函数值f(xk)求和sumTT输出输出积分值T结束结束图图STYLEREF1\s4SEQ图\*ARABIC\s11复合梯形公式算法流程图复合辛普森公式算法流程图开始开始输入输入区间断点a,b及等分数n求出求出步长h，各节点xk，相邻节点中点xk+1/2及相应的函数值f各节点函数值各节点函数值f(xk各相邻节点中心点函数值f(xk+1k=1,2,..n-1SS输出输出积分值S结束结束图图STYLEREF1\s4SEQ图\*ARABIC\s12复合辛普森公式算法流程图MATLAB源程序%复合梯形公式及复合辛普森积分公式clearall;formatlong;a=0;b=1;n=8;h=(b-a)/n;%步长fori=1:n+1x(i)=a+(i-1)*h;ifisnan(sin(x(i))/x(i))symst;tmp=limit(sin(t)./t,t,x(i));%当被积函数在某点值不存在时，求其极限y(i)=eval(tmp);elsey(i)=sin(x(i))/x(i);%被积函数求节点的值endend%复合梯形公式及复合辛普森积分公式s1=0;fork=2:ns1=s1+y(k);endT8=h/2*(y(1)+2*s1++y(n+1))%复合辛普森积分公式s2=0;s3=0;fork=2:2:ns2=s2+y(k);endfork=3:2:n-1s3=s3+y(k);endh1=2*h;%注：此时步长是原来的2倍S4=h1/6*(y(1)+4*s2+2*s3+y(n+1))fprintf('梯形积分公式：%\n辛普森公式积分：%\n',T8,S4)结论与展望图STYLEREF1\s6SEQ图\*ARABIC\s11MATLAB计算结果运行MATLAB程序，得到复合梯形求积公式的积分值为，复合辛普森求积公式的积分值为（四舍五入后保留6位小数）。而实际的积分准确值保留到6位小数的结果为。通过上述结果对比可以得出，虽然复合梯形公式将区间分成了8等分而复合辛普森公式将区间分成了4等分，但两种计算方法实际都需要使用9个点上的函数值，计算量基本也相同，然而最终精度差别却很大。在保留6位小数的前提下，复合辛普森法计算结果与精确解完全一致，而复合梯形公式的计算结果却只有前两位数字与精确解相同，误差相对比较大。下面利用余项公式来估计两种算法的误差。首先需要求fx=,所以有,于是.从而复合梯形公式的误差.而复合辛普森公式的误差.从而，对比两者可得，复合辛普森公式在计算该问题时的精度远高于复合梯形公式。通过以上分析，本文所得结论如下：复合梯形公式和复合辛普森公式都可以用来作为数值积分估算的替代公式。在计算量基本相同的前提下，复合辛普森公式计算结果的计算精度要比复合梯形公式计算精度高的多。本算例也验证了辛普森公式作为偶数阶牛顿-柯特斯公式的更为精确的代数精度。关于如何开展下一步研究，提出以下构想：对多个算例进行分析，保证计算量基本相同的情况下去比较计算精度，验证复合辛普森公式具有更高精度的结论。对多个算例进行MATLAB编程分析，在要求相同计算精度的前提下去比较计算量的大小，从而分析复合梯形公式与复合辛普森公式的优劣。参考文献[1]穆耶赛尔·艾合买提,阿布都热西提·阿布都外力.改进复合梯形求积公式[J].首都师范大学学报(自然科学版),2016,06: