5MATALAB2014a编程与作图合编.ppt

特殊的三维图形函数 1、空间等值线图： contour 3(x,y,z,n) 其中n表示等值线数。 例 山峰的三维和二维等值线图。 解 [x,y,z]=peaks; subplot(1,2,1) contour3(x,y,z,16,s) subplot(1,2,2) contour(x,y,z,16,s) 3、三维散点图 scatter3（X,Y,Z,S,C） 在向量X,Y和Z指定的位置上显示彩色圆圈. 向量X,Y和Z的大小必须相同. 解 输入命令: [x,y,z]=sphere(16); X=[x(:)*.5 x(:)*.75 x(:)]; Y=[y(:)*.5 y(:)*.75 y(:)]; Z=[z(:)*.5 z(:)*.75 z(:)]; S=repmat([1 .75 .5]*10,prod(size(x)),1); C=repmat([1 2 3],prod(size(x)),1); scatter3(X(:),Y(:),Z(:),S(:),C(:),filled),view(-60,60) 例 绘制三维散点图。 7 动画设计 函数m=moviein(n)用来建立一个足够大的n列的帧矩阵m，存放动画中的帧。 Getframe=m(:,i) 函数getframe对当前的图象进行快照 movie函数 movie(m,n)以每秒n幅图形的速度播放由矩阵m的列向量所组成的画面。 例 将曲面peaks做成动画。 解 [x,y,z]=peaks(30); surf(x,y,z) axis([-3 3 -3 3 -10 10]) m=moviein(15); for i=1:15 view(-37.5+24*(i-1),30) m(:,i)=getframe; end movie(m) 【例】 播放一个不断变化的眼球程序段。 m=moviein(20); for j=1:20 plot(fft(eye(j+10))) m(:,j)=getframe; end movie(m,10) 例3 设置一水波动画。 m=moviein(20); x=-1:0.05:1; for j=1:20 y=sqrt(j.^2-x.^2) plot(x,y) hold on z=-sqrt(j.^2-x.^2) plot(x,z) m(:,j)=getframe; end movie(m) Fun(x) x代表数组 * * * * * 几点说明： 1、 若X,Y都为向量,则分别表示点集的横坐标和纵坐标，例如 plot(X,Y)--画实线，S－曲线的性质（线性、粗细、颜色等）。 2、 若X向量,Y为矩阵，则将在同一坐标系下绘制m条曲线，每一行和X的关系将绘制出一条曲线。注意，要求Y矩阵的列数应该等于向量X的长度。 例：x=0:0.1:2; A=[sin(pi*x);0.5+0.5*x]; plot(x,A) %2条不同的 曲线 3、 若X,Y都为矩阵,且假设X和Y的行、列数相同，则将绘制出X矩阵每行和Y矩阵对应行之间关系的曲线。 例：A=[1 4 2;0 -3 4;0 4 3]; B=[1 0 0 ;0 1 0;0 0 1]; plot(A,B) %对应点的连线 符号函数(显函数、隐函数和参数方程)画图 (1) ezplot ezplot(‘x(t)’,’y(t)’,[tmin,tmax]) 表示在区间tminttmax绘制参数方程 x=x(t),y=y(t)的函数图 ezplot(‘f(x)’,[a,b]) 表示在axb绘制显函数f=f(x)的函数图 ezplot(‘f(x,y)’,[xmin,xmax,ymin,ymax]) 表示在区间xminxxmax和 yminyymax绘制 隐函数f(x,y)=0的函数图 例1 在[0,pi]上画y=cosx的图形 解 输入命令 ezplot(sin(x),[0,pi]) 解 输入命令 ezplot(cos(t)^3,sin(t)^3,[0,2*pi]) 解 输入命令 ezplot(exp(x)+sin(x*y)=0,[-2,0.5,0,2]) 2.3.2 图形处理 　1．subplot函数 subplot（