数学实验报告实验序号：2日期：2013年 12月 5日班级2011应数一班姓名孙婉婉学号1101114143实验名称定积分的近似计算问题背景描述： 牛顿莱布尼兹公式仅使用于被积函数的原函数能用初等函数表达出来的情形。如果这点办不到或者不容易办到，就有必要考虑近似计算的方法。在定积分的许多应用问题中，被积函数甚至没有解析表达式，可能只是一条实验记录曲线，或者是一组离散的采样值，这时应用近似方法去计算相应的定积分,以便更好地熟练地掌握定积分的近似计算. 实验目的:1 加深理解积分理论中分割、近似、求和、取极限的思想方法；2 了解定积分近似计算的矩形法、梯形法与抛物线法；3 会用MATLAB语言编写求定积分近似值的程序，会用MALAB中的命令求定积分。实验内容1 分别用梯形法与抛物线法，计算，取n=120.要求使用函数trapz( )、quad( )进行计算求解，并比较结果的差异；2 试计算定积分.（注意：可以运用trapz( )、quad( )或附录程序求解吗？为什么？）；3 学习fuluBsum.m的程序设计方法，尝试用函数sum改写附录C的程序，避免for循环。实验原理与数学模型:实验原理与数学模型： 1 矩形法：根据定积分的定义，每一个积分和都可以看作是定积分的一个近似值，即在几何意义上，这是用一系列小矩形面积近似小曲边梯形的结果，所以把这个近似计算方法称为矩形法不过，只有当积分区间被分割得很细时，矩形法才有一定的精确度针对不同的取法，计算结果会有不同。（1） 左点法：对等分区间，在区间上取左端点，即取。（2）右点法：同（1）中划分区间，在区间上取右端点，即取。（3）中点法：同（1）中划分区间，在区间上取中点，即取。2 梯形法等分区间，相应函数值为 （）曲线上相应的点为 （）将曲线的每一段弧用过点，的弦（线性函数）来代替，这使得每个上的曲边梯形成为真正的梯形，其面积为，于是各个小梯形面积之和就是曲边梯形面积的近似值，即 ，称此式为梯形公式。3 抛物线法将积分区间作等分，分点依次为，对应函数值为（），曲线上相应点为（）现把区间上的曲线段用通过三点，的抛物线来近似代替，然后求函数从到的定积分：由于，代入上式整理后得同样也有将这个积分相加即得原来所要计算的定积分的近似值：，即这就是抛物线法公式，也称为辛卜生（Simpson）公式 实验所用软件及版本：2012B主要内容： 1，分别用梯形法与抛物线法，计算，将积分区间1,2作120等分。并尝试用函数trapz(),quad()进行算求解，比较结果的差异。2，试计算定积分.(注意：可以运用trapz()、quad()、或附录程序求解吗？为什么？) 3，学习fuluBsum.m的程序设计方法，尝试用函数sum改写矩形法和抛物线法的程序，避免for循环。实验过程记录（含基本步骤、主要程序清单及异常情况记录等）：1、（1）梯形法：、format longn=120;a=1;b=2;inum=0;syms x fxfx=1/x; for i=1:n xj=a+(i-1)*(b-a)/n; xi=a+i*(b-a)/n; fxj=subs(fx,x,xj); fxi=subs(fx,x,xi); inum=inum+(fxj+fxi)*(b-a)/(2*n);endinumx=1:1/120:2;y=1./x;trapz(x,y)（2）抛物线法： format longn=120;a=1;b=2;inum=0;syms x fxfx=1/x;for i=1:n xj=a+(i-1)*(b-a)/n; xi=a+i*(b-a)/n; xk=(xi+xj)/2; fxj=subs(fx,x,xj); fxi=subs(fx,x,xi); fxk=subs(fx,x,xk); inum=inum+(fxj+4*fxk+fxi)*(b-a)/(6*n);endinum quad(1./x,1,2)2、（1）符号求积分：int(sin(x)/x,x,0,inf)（2）quad(sin(x)./x,0,inf)3 ex2zhi2tixingfa.m format long n=120;a=1;b=2;inum=0;syms x fxfx=1./x;for i=1:n xj=a+(i-1)*(b-a)/n; xi=a+i*(b-a)/n; fxj=subs(fx,x,xj); fxi=subs(fx,x,xi); inum=inum+(fxj+fxi)*(b-a)/(2*n);end inum integrate=int(fx,1,2)integrate=double(integrate)fprintf(The relative error between inum and real-value is about:%gnn,.abs(inum-integrate)/integrate)ex2zhi2paowuxianfa.mformat long n=120;a=1;b=2;m=0;p=0;syms x fxfx=1./x;for k=1:n i=2*k-1; j=2*k; xi=a+(b-a)*i/(2*n); xj=a+(b-a)*j/(2*n); fxi=subs(fx,x,xi); fxj=subs(fx,x,xj); m=m+fxi;p=p+fxj;end fxa=subs(fx,x,a); fxb=subs(fx,x,b); inum=(b-a)*(fxa+fxb+4*m+2*p)/(6*n)integrate=int(fx,1,2)integrate=double(integrate)fprintf(The relative error between inum and real-value is about:%gnn,.abs(inum-integrate)/integrate)ex2zhi2paowuxianfa.mformat long n=120;a=1;b=2;m=0;p=0;syms x fxfx=1./x;for k=1:n i=2*k-1; j=2*k; xi=a+(b-a)*i/(2*n); xj=a+(b-a)*j/(2*n); fxi=subs(fx,x,xi); fxj=subs(fx,x,xj); m=m+fxi;p=p+fxj;end fxa=subs(fx,x,a); fxb=subs(fx,x,b); inum=(b-a)*(fxa+fxb+4*m+2*p)/(6*n)integrate=int(fx,1,2)integrate=double(integrate)fprintf(The relative error between inum and real-value is about:%gnn,.abs(inum-integrate)/integrate)改写矩形法ex6gaixiejuxingfa.mformat long n=100;a=0;b=1; syms x fxfx=1/(1+x2);i=1:n;xj=a+(i-1)*(b-a)/n;xi=a+i*(b-a)/n;fxj=subs(fx,x,xj);fxi=subs(fx,x,xi);fxij=(fxi+fxj)/2;m=fxj*(b-a)/n;p=fxi*(b-a)/n;k=fxij*(b-a)/n;inum1=sum(m)inum2=sum(p)inum3=sum(k)改写抛物线法 ex6gaixiepaowuxianfa.mformat long n=100;a=0;b=1; syms x fxfx=1/(1+x2);i=0:(n-1);xj=a+(2*i)*(b-a)/(2*n);xi=a+(2*i+1)*(b-a)/(2*n);xk=a+(2*i+2)*(b-a)/(2*n);fxj=subs(fx,x,xj);fxi=subs(fx,x,xi);fxk=subs(fx,x,xk);m=(fxj+4*fxi+fxk)*(b-a)/(6*n);inum=sum(m)实验结果报告与实验总结：run(E:matlabex2zhi2tixingfa.m)inum = 0.693integrate =log(2)integrate =0.693The relative error between inum and real-value is about:6.26164e-006 0.693integrate =log(2)integrate =0.693The relative error between inum and real-value is about:6.26164e-006run(E:matlabex2zhi2paowuxianfa.m)inum = 0.69453606945825integrate =log(2)integrate = 0.693The relative error between inum and real-value is about:0.00200374run(E:matlabex2zhi2trapz.m)inum = 0.693integrate =log(2) integrate = 0.693The relative error between inum and real-value is about:6.26164e-006、（1）梯形法 inum = 0.6938 ans = 0.6938(2)抛物线法 inum = 0.6934 ans = 0.69302、