资源描述：
《matlab上机实践六(计算物理班)》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、实验要求：为达到理想的实验效果，同学们务必做到：(1)实验前认真准备，要根据实验目的和实验内容，复习好实验中可能要用到的命令，越好编程的思路，做到胸*成竹，提离上机效卓。(2)实验过程中积极思考，耍深入分析侖令.程序的执彳亍结果以及各种屏幕信息的含义.出现的冰因并提出解决办法。(3)实验后认真总结，要总结本次实验希哪罢收获，还存在哪些问题，并写出实验报吿。实验报吿应包括实验目的.实验内容.流程图(校大程序).程序(命令)淸单、运行结果以及实验的收获与体会等内容。同学们在上机过程中会碰到各种各样的问题，分析问题和
2、解决问题的过程就是积累经验的过程。只要同学们按照上面3点要求去做，在学完本课程后就一定会有很大的收获。实验仪器：计算机实验时间：2018.5.实验原理：1.Matlab数值积分的实现:辛普牛方法：z=quad(@fname,a,b,tol,trace)z=quadl(@fname,a,b,tol,trace)高斯■克朗罗德法：[I,err]=quad(@fname,a,b)梯形法：z二trapz(X,Y)二重和三重积分：z=dblquad(@fname,a,b,c,d,tol)z=triplequad(@fna
3、me,a,b,c,d,e,f,tol)2.Matlab数值微分实现：DX=diff(X)DX=diff(X,n)DX=diff(X,n,dim)实验目的：1.掌握matlab数值积分方法2.掌握matlab数值微分方法实验内容:1.用数值方法求定积分。（每个题要求用3种方法计算。）(1)Ix=Jcosf+4sin⑵)力程序：»f=@(t)sqrt(cos(t)+4*sin(2*t));»[S,n]=quad(f,0,2*pi)S=4.8113+4.8113i325ln(l+x)」—axl+x2»f=@(x)lo
4、g(l+x)./(l+x.A2);»[s,n]=quad(f,04)0.272217程序：»x=0:0.01:l;»I2=log(l+x)./(l+x.A2);»trapz(x,I2)ans=0.2722⑶厶=jtcos(t)e~l2,dt»t=-pi:pi/100:pi;»I3=t.*cos(t).*exp(-i.*2.*t);»trapz(t,I3)ans=-0.0000-4.1878i程序：>>f=@(t)t・*cos(t)・*exp(-i*2*t);>>[S,n]=quad(f,-pi,pi)S二0-4
5、.1888in二161求二重和三重定积分dxdy»f=@(x,y)exp(x.A2+y.A2);»T=dblquad(f,0,2,0,1)24.0645zex^ydxdydz2.求下列函数的数值导数八x）（x范围・2到2）,然后绘制f，（x）的图形，注:利用向前差分计算（步长取0.01）ox2+lf=inline(1exp(x)・/(x•人2+1)');x=-2:0.01:2;dx=diff(f([x,2.01]))./0.01plot(xzdx,1-1);1.4
