新整理的MATLAB实验报(要打印的).docVIP

实验一 实验目的： 熟悉MATLAB编译环境； 熟悉MATLAB简单使用方法； 掌握MATLAB程序的调试及其技巧 实验内容： 求出的具有64位有效数字的积分值。 在时间区间[0,10]中，绘制曲线。 先运行产生矩阵Z。（1）矩阵Z中有都少个“非数数据”？（2）用指令观察所绘图形；（3）请写出绘制相应的无裂缝图形指令。 采用数值计算，画出在[0,10]区间曲线。 实验结果： 1．syms x; f=exp(-abs(x))*abs(sin(x)); int(f,x,-5*pi,1.7*pi); F=vpa(ans,64) F = 1.087849417255503701102633764498941389696991336803454392428439159 2．syms t; t=0:0.01:10; y=1-exp(-0.5*t).*cos(2*t); plot(t,y) 3．x=-3*pi:pi/15:3*pi; xx=x+(x==0)*eps; yy=xx; [X,Y]=meshgrid(xx,yy); warning off; Z=sin(X).*sin(Y)./X./Y; surf(X,Y,Z); shading interp 4．t=0:0.01:10; tt=t+(t==0)*eps; yy=sin(tt)./tt; sx=0.01*cumtrapz(yy); plot(t,sx) xlabel(x),ylabel(y); grid on [x,y]=ginput(2) 实验二 实验目的： 熟悉MATLAB编译???境； 熟悉MATLAB简单使用方法； 掌握MATLAB程序的调试及其技巧 实验内容： 1. 2.用绿实线绘制的三维曲线。 3．在区域，绘制的图。 4．在区域里，绘制图。5．n=[0.2:0.2:1.4]; 1. t=[0:0.1:18]; L=length(n); for k=1:L; m=sqrt(abs(1-n(k).^2)); if k=1k5 y=(1./m).*exp(-(n(k).*t)).*sin(m.*t); plot(t,y) xlabel(t),ylabel(y) text(0.2,0.7,\fontsize{14}\it\zeta=0.2); hold on end if k=7k5 y=(1./(2*m)).*(exp((-(n(k)-m)).*t)-exp(-(n(k)+m).*t)) ; plot(t,y,r) xlabel(t),ylabel(y) text(0.2,0.1,\fontsize{14}\it\zeta=1.4); text(10,0.7,\fontsize{14}\it\delta\zeta=0.2); hold on end if k==5 y=(1./m).*exp(-(n(k).*t)).*sin(m.*t) ; plot(t,y,k,LineWidth,3) xlabel(t),ylabel(y) hold on end end grid on 2．t=0:0.1:15; x=sin(t); y=cos(t); z=t; xlabel(x),ylabel(y),zlabel(z); plot3(x,y,z,-g); box on 3．x=-3:0.1:3; y=x; [X,Y]=meshgrid(x,y); Z=4.*X.*exp(-X.^2-Y.^2); mesh(X,Y,Z),colormap(gray); hidden off view(-32,30) 4．x=-4*pi:4*pi; y=x; [X,Y]=meshgrid(x,y); Z=sin(X+Y)./(X+Y); surf(X,Y,Z) view(33,30) shading interp;