文档介绍：该【matlab实验报告--定积分的近似计算-- 】是由【幸福人生】上传分享，文档一共【20】页，该文档可以免费在线阅读，需要了解更多关于【matlab实验报告--定积分的近似计算-- 】的内容，可以使用淘豆网的站内搜索功能，选择自己适合的文档，以下文字是截取该文章内的部分文字，如需要获得完整电子版，请下载此文档到您的设备，方便您编辑和打印。matlab实验报告--定积分的近似计算--数学实验报告实验序号:2日期:2013年11月30日班级应数二班姓名丁慧娜学号1101114088实验名称定积分的近似计算实验所用软件及版本MATLABR2012b问题背景描述:利用牛顿—莱布尼兹公式虽然可以精确地计算定积分的值,,,被积函数甚至没有解析表达式,可能只是一条实验记录曲线,或者是一组离散的采样值,:本实验将主要研究定积分的三种近似计算算法:矩形法、梯形法、抛物线法。加深理解积分运算中分割、近似、求和、取极限的思想方法。,尝试用函数sum改写附录1和附录3的程序,避免for循环。实验原理与数学模型:1.?矩形法根据定积分的定义,每一个积分和都可以看作是定积分的一个近似值,即每个上的曲边梯形成为真正的梯形,其面积为,.于是各个小梯形面积之和就是曲边梯形面积的近似值,,即,称此式为梯形公式。?3.?抛物线法将积分区间作等分,分点依次为,,对应函数值为(),曲线上相应点为().现把区间上的曲线段用通过三点,,的抛物线来近似代替,然后求函数从到的定积分:由于,代入上式整理后得同样也有……将这个积分相加即得原来所要计算的定积分的近似值:,即这就是抛物线法公式,也称为辛卜生(Simpson)公式.?主要内容(要点):1.?分别用梯形法与抛物线法,计算,()、quad()进行计算求解,.?试计算定积分.(注意:可以运用trapz()、quad()或附录程序求解吗?为什么?)3.?,尝试用函数sum改写附录1和附录3的程序,避免for循环。实验过程记录(含基本步骤、主要程序清单及异常情况记录等):1:梯形法formatlongn=120;a=1;b=2;symsxfxfx=1/x;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);%所有右点值f=(fxi+fxj)/2*(b-a)/n;%梯形面积inum=sum(f)%加和梯形面积求解integrate=int(fx,1,2);integrate=double(integrate)fprintf('Therelativeerrorbetweeninumandreal-valueisabout:%g/n/n',...abs((inum-integrate)/integrate))【调试结果】>>TXFinum==-valueisabout:-06/n/n>>抛物线法:%抛物线法formatlongn=120;a=1;b=2;inum=0;symsxfxfx=1/x;fori=1:nxj=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);endinumintegrate=int(fx,1,2);integrate=double(integrate);fprintf('Therelativeerrorbetweeninumandreal-valueisabout:%g/n/n',...abs((inum-integrate)/integrate))【调试结果】>>clear>>PWXFinum=-valueisabout:-11/n/n>>使用函数trapz()x=1:1/120:2;y=1./x;trapz(x,y)【调试结果】ans=