LAB05_复化求积实验.docVIP

Lab5．Matlab语言编写按复化梯形公式和复化Simpson公式进行数值积分的程序； 2．掌握用复化梯形公式和复化Simpson公式进行数值积分的方法。 3．对，用所编写的程序计算T8和S4，通过结果的比较和理论分析，以了解复化梯形公式和复化Simpson公式精度。 【实验内容】 1．根据Matlab语言特点，描述复化梯形公式和复化Simpson公式。 2．用Matlab语言编写按复化梯形公式和复化Simpson公式求积分的程序。 3．对于，用所编写的程序计算T8和S4。并比较计算结果的精度。 【实验仪器与软件】 1．CPU主频在1GHz以上，内存在128Mb以上的PC； 2．Matlab 6.0及以上版本。 实验讲评： 实验成绩： 评阅教师： 200 年 月 日Lab5．impson公式 将区间[a,b]划分为N等分，分点Xk=a+(k-1)h,h=(b-a)/n,k=1,2,…,n,n+1,在每个子区间[Xk,Xk+1](k=1,…n)上采用梯形公式，则得I=把称为复化梯形公式，其余项为 将区间划分为N等分，且N=2M，在每个子区间上采用辛普森公式， 若记则得I= 记称为复化辛普森求积公式，其余项为： 复化梯形和复化Simpson的程序 x=a; T=0; for i=1:(n-1) x=x+h; T=T+f(x); end s=(f(a)+2*T+f(b))*h/2; %复化Simpson程序 function s=fhxps(a,b,n)%[a,b] %为积分区间，n为等分数 h=(b-a)/n; x=a; x2=a+h/2; T=0; S=f(x2); for i=1:(n-1) x2=x2+h; x=x+h; S=S+f(x2); T=T+f(x); end s=(f(a)+4*S+2*T+f(b))*h/6; 数值计算实验 function f=fun(x) f=x./(4+x.^2); %利用梯形求积公式： a=0; b=1; n=8; T=ftxing(a,b,n) 运行计算得：T=0.1114 代数精度为1次 %利用辛普森求积公式： a=0; b=1; n=4; S=fxpin(a,b,n) 运行计算得：S=0.1116 代数精度为3次 四、实验总结 答：通过实验发现，利用辛普森公式进行数值计算所得结果的代数精度为三次。而梯形公式计算的记过精度没有辛普森公式计算的精度高，在数值计算中，我们应该尽可能的使计算结果的精度要。 姓名：余万龙 学号：2010051113 班级：10数信