资源描述:
《matlab 绘图技巧》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、一、目录1.彗星图 二维彗星图 三维彗星图2.帧动画3.程序动画4.色图变换5.Voronoi图和三角剖分 Voronoi图 三角剖分6.四面体7.彩带图 彩带图 三维流彩带图8.伪彩图9.切片图 切片图 切片轮廓线图10.轮廓图 显示轮廓线 显示围裙 瀑布效果 带光照模式的阴影图11.函数绘图 轮廓线、网格图、曲面图、轮廓网格图 轮廓曲面图、二维曲线、极坐标曲线图、自定义函数12.三维图形控制 视点 灯光效果 色彩控制二、图形示例1.彗星图二维彗星图t=0:.01:2*p
2、i;x=cos(2*t).*(cos(t).^2);y=sin(2*t).*(sin(t).^2);comet(x,y);title('二维彗星轨迹图')holdonplot(x,y) 三维彗星图a=12;b=9;T0=2*pi;%T0是轨道的周期T=5*T0;dt=pi/100;t=[0:dt:T]';f=sqrt(a^2-b^2);%地球与另一焦点的距离th=12.5*pi/180;%未经轨道与x-y平面的倾角E=exp(-t/20);%轨道收缩率x=E.*(a*cos(t)-f);y=E.*(b*co
3、s(th)*sin(t));z=E.*(b*sin(th)*sin(t));plot3(x,y,z,'g')%画全程轨线holdon,sphere(20);%画地球axisofftitle('卫星返回地球示例')x1=-18*T0;x2=6*T0;y1=-12*T0;y2=12*T0;z1=-6*T0;z2=6*T0;axis([x1x2y1y2z1z2])%axis([-1510-1510-1010])axisequalcomet3(x,y,z,0.02);%画运动轨线holdoff 2.帧动画Z=pea
4、ks;surf(Z)%绘制网格表面图axistightset(gca,'nextplot','replacechildren');%设定axis覆盖重画模式title('帧动画播放示例')forj=1:20 surf(sin(2*pi*j/20)*Z,Z)%重新绘制网格表面图,这里后面一个Z当成了颜色矩阵 F(j)=getframe;%创建帧endmovie(F,20)%播放动画20次 3.程序动画t=0:pi/50:10*pii=1;h=plot3(sin(t(i)),cos(t(i)),t(i)
5、,'*','erasemode','none');%设定擦除模式gridonaxis([-22-22-110*pi])title('程序动画示例')fori=2:length(t) set(h,'xdata',sin(t(i)),'ydata',cos(t(i)),'zdata',t(i)); drawnow pause(0.01)end 4.色图变换loadspineimage(X)colormapcoolspinmap(10) 5.Voronoi图和三角剖分Voronoi图rand('sta
6、te',5)x=rand(1,10);y=rand(1,10);subplot(131)voronoi(x,y);%绘制voronoi图形axisequalaxis([-0.21.6-0.52.5])subplot(132)[vx,vy]=voronoi(x,y);plot(x,y,'r+',vx,vy,'b-');%应用返回值绘制axisequalaxis([-0.21.6-0.52.5])subplot(133)rand('state',5);x=rand(10,2);[v,c]=voronoin(x)
7、;%返回值v参数维voronoi顶点矩阵,返回值c参数为voronoi元胞数组fori=1:length(c) ifall(c{i}~=1) patch(v(c{i},1),v(c{i},2),i);%应用色图i endendaxisequalaxis([-0.21.6-0.52.5])boxon 三角剖分 [x,y]=meshgrid(1:15,1:15);tri=delaunay(x,y);z=peaks(15);trimesh(tri,x,y,z) 6.四面
8、体d=[-11];[x,y,z]=meshgrid(d,d,d);%定义一个立方体x=[x(:);0];y=[y(:);0];z=[z(:);0];%[x,y,z]分别为加上中心的立方体顶点X=[x(:)y(:)z(:)];Tes=delaunayn(X);%返回m×n的数组值tetramesh(Tes,X);%绘制四面体图camorbit(20,0);%旋转摄像目标位置 7.彩带图彩带图[x,y]=meshgr