2MATLAB作图全解.ppt

MATLAB作图是通过描点、连线来实现的，故在画一个曲线图形之前，必须先取得该图形上的一系列的点的坐标（即横坐标和纵坐标），然后将该点集的坐标传给MATLAB函数画图. 例 在[0,2 ]用红线画sin x,用绿圈画cos x. x=linspace(0,2*pi,30); y=sin(x); z=cos(x); plot(x,y,r,x,z, g0) plot3(x,y,z,s) plot3(x,y,z) 空 间 曲 面 1. 在图形上加格栅、图例和标注 （1）GRID ON: 加格栅在当前图上 GRID OFF: 删除格栅 （3） hh = gtext(‘string’) 命令gtext(‘string’)用鼠标放置标注在现有的图上. 运行命令gtext(‘string’)时，屏幕上出现当前图形，在 图形上出现一个交叉的十字，该十字随鼠标的移动移动， 当按下鼠标左键时，该标注string放在当前十交叉的位 置. x=1200:400:4000;y=1200:400:3600;[xx,yy]=meshgrid(x,y);zz=[1130,1250,1280,1230,1040,900,500,700;1320,1450,1420,1400,1300,700,900,850;1390,1500,1500,1400,900,1100,1060,950;1500,1200,1100,1350,1450,1200,1150,1010;1500,1200,1100,1550,1600,1550,1380,1070;1500,1550,1600,1550,1600,1600,1600,1550;1480,1500,1550,1510,1430,1300,1200,980;]; surf(xx,yy,zz);hold on;[c,h]=contour(xx,yy,zz,6);clabel(c,h) %问题1. clear ax1=subplot(1,2,1); t=0:.01*pi:2*pi; f1=polar(t,1+cos(t)); title(心形线); legend(f1,r=1+cost) ax2=subplot(1,2,2); x1=linspace(-40,20,40); y1=linspace(-40,20,40); [xx1,yy1]=meshgrid(x1,y1); zz1=xx1.*yy1; f2=surf(xx1,yy1,zz1); title(马鞍面); legend(f2,z=x^2-y^2,Location,NorthWest) xlabel(x轴) ylabel(y轴) zlabel(z轴) axis([-40 20 -40 20 -1000 1500]) setAllowAxesRotate(rotate3d,ax1,false); grid on %问题2(给你四个程序，自己挑吧)：%1)clear[x,y,z]=sphere(200);i=find(x.^2+y.^2-x=0x.^2+y.^2+z.^21);z(i)=NaN;mesh(x,y,z)hold ont=0:pi/20:2*pi; x=(cos(t)+1)/2; y=sin(t)/2; z=linspace(-1.2,1.2,length(t)); X=meshgrid(x); Y=meshgrid(y); Z=[meshgrid(z)]; surf(X,Y,Z) xlabel(x),ylabel(y),zlabel(z) title(球面 X^2+Y^2+Z^2=1和圆柱面X^2+Y^2=X所围成的区域)rotate3d on %2)clear[X,Y,Z]=sphere(20);mesh(X,Y,Z)hidden onhold onaxis equalt=0:pi/20:2*pi; x=(cos(t)+1)/2; y=sin(t)/2; z=linspace(-1.2,1.2,length(t)); X=meshgrid(x); Y=meshgrid(y); Z=[meshgrid(z)]; surf(X,Y,Z) xlabel(x),ylabel(y),zlabel(z) rotate3d on %3)[x,y,z]=sphere(200);i=find(x.^2+y.^2-x=0x.^2+y.^2+z.^2=1);z(i)=NaN;mesh(x,y,z)hidden onhold onezplot(x^2+y^2=x)