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

SHANGHAIJIAOTONGUNIVERSITY 题目名称：复合梯形公式与复合辛普森公式比较 学生姓名: 学生学号: 班级: 学院(系): 目录 1.概括3 2.问题提出4 3.算法推导5 4.算法框图6 4.1复合梯形公式算法流程图6 4.2复合辛普森公式算法流程图7 5.MATLAB源程序8 6.结论与展望9 图表目录 图4-1复合梯形公式算法流程图..................................6 图4-2复合辛普森公式算法流程图................................7 图6-1MATLAB计算结果.......................................... 9 表2-1函数计算结果表..........................................4 概括 梯形求积公式和辛普森求积公式分别是牛顿-科斯特公式中n=1和n=2时的情形。其中 梯形求积公式可表示为 其公式左端是以[a,b]区间上积分，右端为b-a为高、端点函数值为上下底的梯形的面 积值，故通称为梯形公式，具有1次代数精准度。 近似的，辛普森求积公式能够表示为 该公式一般在立体几何中用来求拟柱体的体积，由于偶数n阶牛顿-科特斯求积公式至 少具有n+1次代数精准度，所以辛普森公式实际上具有3次代数精准度。 由于牛顿-科斯特公式在n≥8时不具有稳定性，故不可能通过提高阶的方法来提高求 积精度。为了提高精度平时可把积分区间分红若干子区间（平时是平分），再在每个子区间上用低阶求积公式。这种方法称为复合求积法。 本文主要议论复合梯形公式和复合辛普森公式在同一数学问题中的应用。首先给出了 复合梯形公式和复合辛普森公式的推导过程以及其余项的表达形式，然后用流程图的形式 介绍算法思路，再运用MATLAB编写代码计算结果，最后对结果进行比较议论。 希望通过两个算法在同一个算例中的应用比较，更好的理解和掌握复合梯形公式和复 合辛普森公式的适用范围和适用条件。并且能够熟悉MATLAB编程求解问题的流程，掌握编 程化的思想方法。同时对两种方法的计算结果比较解析，议论两种求积方法的计算精度。 问题提出 ( ) = ???????? 对于函数fx ?? 1???????? 积分I=∫????。 0?? 表2-1函数计算结果表 xf(x) 01 1/8 1/4 3/8 1/2 5/8 3/4 7/8 1 算法推导 3.1复合梯形公式 根据梯形公式， 将区间[a,b]划分为n等份，分点????=??+????,h= ??-?? ??,k=0,1, ,n，在每个子区间 [??, ??+1](k=0,1,,n-1)上采用梯形公式，则得： ?? ?? 记 则T为复合梯形公式。 ?? 其他，复合梯形公式的余项可表示为 3.2复合辛普森公式 根据辛普森公式 将区间[a,b]划分为n等份，在每个子区间[????,????+1](k=0,1,,n-1)上采用辛普森公式。 若记 1 ????+1/2=????+2?? 则得 记 该公式即为复合辛普森公式。 复合辛普森公式的余项可表示为 4.算法框图 4.1复合梯形公式算法流程图 输入区中止点 求出步长h，各节点????及 各节点函数值f(xk)求和 ?? T??=2(??(??)+??(??) 图4-1复合梯形公式算法流程图 4.2复合辛普森公式算法流程图 输入区中止点a, 求出步长h，各节点????，相邻节点中点 各节??点??函+1数/2及值相f应(x的k)求函和数值sum1, 各相邻节??点中心点函数值f(x)求和 S??=6(??(??)+??(??)+2??????1 ?????? 输出积分值S 图4-2复合辛普森公式算法流程图 MATLAB源程序 %复合梯形公式及复合辛普森积分公式 clearall; formatlong; a=0; b=1; n=8; h=(b-a)/n;%步长 fori=1:n+1 x(i)=a+(i-1)*h; ifisnan(sin(x(i))/x(i)) symst; tmp=limit(sin(t)./t,t,x(i));%当被积函数在某点值不存在时，求其极限 y(i)=eval(tmp); else y(i)=sin(x(i))/x(i);%被积函数求节点的值 end end %复合梯形公式及复合辛普森积分公式 s1=0; fork=2:n s1=s1+y(k); end T8=h/2*(y(1)+2*s1++y(n+1)) %复合辛普森积分公式 s2=0; s3=0; fork=2:2:n s2=s2+y(k); end fork=3:2:n-1 s3=s3+y(k); end h1=2*h;%注：此时步长是原来的2倍