实验二matlab实验报告

实验序号 日期 年11 19计算，Matlab有专门函数可用。矩形法．但是，只有当积分区间被分割得很细时，矩形法才有一定的精确度．针对不同的取法，有三种办法：左点法、右点法、中点法，不同的办法计算成果会有不同。

，即，sum(a)aformatlong15（注：由于本实验要比较近似解法和精确求解间的误差，需要更高的精度double()ASCII码；若输入的是整型数值则转化为对应的实型数quad()quad(fun,a,b)fun是函数，并且为数值形式的，因此使用*、/、^等运算时要在其前加上小数点，即.*、./、.^等．x为带有步长的积分区间；y为数值形式的运算（dblquad()格式：dblquad(fun,xmin,xmax,ymin,ymax)，funinline定义，也能够通过某个函数文献的句柄fprintf（文献地址，格式，写入的变量）：subs(f，'x'，a)afxsubs(f)syms12重要内容（要点 ．并尝试直接使用函数trapz()、quad()进行计算试计算定积 （fulu2sum.msum13forformatlongsymsxfxfori=1:nfprintf('Therelativeerrorbetweeninumandreal_valueisabout:%g\n\n',abs((inum-trapz(x,fx))/trapz(x,fx)))办法二：（抛物线法formatlongsymsxfxfori=1:n about:%g\n\n',abs((inum-quad('1./x',1,2))/quad('1./x',1,2)))办法一（办法二（formatlongsymsxfxfori=1:nfprintf('Therelativeerrorbetweeninum1andreal-valueisabout:%g\n\n',...fprintf('Therelativeerrorbetweeninum2andreal-valueisabout:%g\n\n',...fprintf('Therelativeerrorbetweeninum3andreal-valueisabout:%g\n\n',...formatlongsymsxfxf1=(fxi)*(b-fprintf('Therelativeerrorbetweeninumandreal_valueisabout:%g\n\n',abs((inum-办法二（sum函数——抛物线法formatlongsymsxfxquad('1./x',1,2,1e- about:%g\n\n',abs((inum-quad('1./x',1,2,1e-16))/quad('1./x',1,2,1e-16)))办法三（sum函数——梯形法）formatlongsymsxfxfprintf('Therelativeerrorbetweeninumandreal_valueisabout:%g\n\n',abs((inum-fulu2sum.m的程序设计办法。于是学1/2*pi。2bbinfb改为较大的1001000来替代无穷大，即得到成果。但误差较大，错误成果为：inumintegrate=Therelativeerrorbetweeninum1andreal-valueisabout:boinumfTherelativeerrorbetweeninumandreal_valueisabout:0inumans3办法（1）ans办法（2）ans办法（3）inumintegrate=Therelativeerrorbetweeninum1andreal-valueisabout:2.38752e-Therelativeerrorbetweeninumandreal_valueisabout:2.78348e-0083办法一inumfTherelativeerrorbetweeninumandreal_valueisaboutinumansTherelativeerrorbetweeninumandreal_valueisabout:1.35886e-inum1inum2inum3fTherelativeerrorbetweeninumandreal_valueisabout:0Therelativeerrorbetweeninumandreal_valueisaboutTherelativeerrorbetweeninumandreal_valueisa