复合辛普森.docVIP

实验题目：复合梯形及复合辛普森求积公式 学生姓名： 专 业： 学 号： 完成日期： 实验题目： 数值分析之复合梯形及复合辛普森求积公式 学生姓名： 学号： 完成日期： 2011.11.20 1 实验目的 利用不同数值方法计算积分。 2 实验步骤 2.1 算法原理 复合梯形原理： 将区间划分为n等分，分点在每个子区间上采用梯形公式，则得 记 ，称为复合梯形公式，其余项可得 由于 所以使 于是复合梯形公式的余项为 可以看出误差是阶，且当 即复合梯形公式是收敛的。事实上只要设，则可得到收敛性，因为只要把 当时，上式右端括号内的两个和式均收敛到积分，所以复合梯形公式收敛。此外，的求积系数为正，则复合求积公式是稳定的。 复合辛普森原理： 将区间划分为n等分，在每个子区间上采用辛普森公式，若记则得 记 称为复合辛普森求积公式，其余项为 于是当时，与复合梯形公式相似有 易知误差阶为，收敛性是显然的，实际上，只要则可得到收敛性，即 此外，由于中求积公系数均为正数，故知辛普森公式计算稳定。 ·2.2 算法步骤 复合梯形： 首先将区间划分为n等分，分点在每个子区间上采用梯形公式，则得 复合辛普森： 首先将区间划分为n等分，在每个子区间上采用辛普森公式，若记则得 ·2.3 程序流程图 ·3 实验结果分析 4 实验结论 复合梯形及复合辛普森在求积计算积分时，复合辛普森的精度比复合梯形要高。 5 实验心得体会 此次实验，我体会到了不同的求积公式得到的结果精度不同，并且可以利用计算机解决一些比较复杂的积分。 参考文献 白峰杉.数值分析引论.北京：高等教育出版社,2004 李庆扬,王能超,易大义.数值分析.武汉：华中理工大学出版社,1986. 附录（源代码） #include stdio.h #include math.h float f(float x) { float y; y=float(sqrt(x)*log(x)); return y; } void main() { float a=0.000001; int b=1; int n; printf(输入划分数目n: ); scanf(%d,n); int k; float h; h=float((b-a)/n); float ff1,ff2,ff11,ff22; float f1=0; float f2=0; for(k=1;kn-1;k++) { f1=f1+f(a+k*h); } for(k=0;kn-1;k++) { f2=f2+f(1.0/2.0*h+k*h); } ff1=h/2*(f(a)+2*f1+f(b)); ff11=fabs(ff1-(-4/9.0)); ff2=h/6*(f(a)+4*f2+2*f1+f(b)); ff22=fabs(-4/9.0-ff2); printf(复合梯形 Tn=%f 误差%f,ff1,ff11); printf(\n); printf(复合辛普森 Sn=%f 误差%f,ff2,ff22); printf(\n); } 结束 输出 ， ， 输入数据a,b,n 开始