复合辛普森求积

《复合辛普森求积》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、第三次实验实验名称复合辛普森求积计算积分实验时间2012.05.06姓名班级数应二班学号成绩一、实验目的，内容1.目的：学习并理解复合辛普森求积计算积分的matlab实现。 2.内容：用matlab计算积分（精确值-0.4444），并求出达到。二．代码1.functionSn=ComSimpson(a,b,n)%复合辛普森求积%f表示被积函数，本题中由f.m调用%a,b分别表示积分上下限%n表示区间分割次数%sn表示该方法计算所返回的值h=(b-a)/n;fork=0:nx(k+1)=a+k*h;x_k(k+1)=x(k+1)+1/2*h;if(

2、x(k+1)==0)

3、(x_k(k+1)==0)x(k+1)=10^(-10);x_k(k+1)=10^(-10)；%误差精度endendS1=h/6*(f(x(1))+f(x(n+1)));%复合辛普森求积公式第一部分fori=2:nF_1(i)=h/3*f(x(i));endforj=1:nF_2(j)=2*h/3*f(x_k(j));endS2=sum(F_1)+sum(F_2);%复合辛普森求积公式第二部分Sn=S1+S2;%算得值f的表达式由f.m（见附）文件调用。附：functiony=f(x)y=sqrt(x).*log(x);三.

4、数值结果:在命令窗口输入如下指令时，出现如下结果Sn=ComSimpson(0,1,300)Sn=-0.4438Sn=ComSimpson(0,1,700)Sn=-0.4442四．计算结果的分析:当步长取值很小的时候，误差较大。当步长取得越细，即区间分割的越小时，实验计算值的精度越高，即越趋近精确值。五.计算中出现的问题，解决方法及体会：本实验过程中，了解了复合求积公式的几个构成，以及在matlab中的实现，深化了对于该问题的理解。一的操作中，本来个个很简单的问题，哪怕就是一个小小符号的不注意，也不会得到结果。实验的严谨性，细腻性有待进一步提升。

5、对于实验，网上的程序有不少，关键是要看懂，弄懂，在实际的问题中能运用。借鉴吸收的同时，学会编程。多练习，才会有效果。教师评语指导教师：年月日