MATLAB-图像处理之三维绘图.ppt

《MATLAB-图像处理之三维绘图.ppt》由会员分享，可在线阅读，更多相关《MATLAB-图像处理之三维绘图.ppt（72页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、上次课内容回顾上次课内容回顾MATLAB二维绘图二维绘图离散函数和数据的可视化；离散函数和数据的可视化；二维曲线和图形二维曲线和图形plot基本调用格式、基本调用格式、衍生调用格式、衍生调用格式、可控调用格式；可控调用格式；坐标控制和图形标识（坐标轴、分格线、坐标坐标控制和图形标识（坐标轴、分格线、坐标框、图形标识、精细控制）框、图形标识、精细控制）多次叠绘、双纵坐标和多子图；多次叠绘、双纵坐标和多子图；1.获取二维图形数据的指令获取二维图形数据的指令ginput。3.3三维曲线和曲面三维曲线和曲面内容提要内容提要三维线图指令三维线图指令plot3；三维曲面三维曲面/网线图；网线图；颜色的使用

2、；颜色的使用；曲面曲面/网线图的精细修饰；网线图的精细修饰；透视、镂空和裁切透视、镂空和裁切高维可视化高维可视化1. 图形窗功能图形窗功能1. 三维线图指令三维线图指令plot3 plot3(X,Y,Z,s) （单参数三维曲线）（单参数三维曲线） 用用s 指定的点形线型色彩绘制曲线指定的点形线型色彩绘制曲线 plot3(X1,Y1,Z1,s1,X2,Y2,Z2,s2, . )用用s1, s2 指定的点形线型色彩绘制指定的点形线型色彩绘制多类曲线多类曲线 说明说明：X、Y、Z为同维向量、为同维向量、矩阵、矩阵、s的意义、的意义、四元组之间没有约束关系四元组之间没有约束关系例例5.3-1利用三维曲

3、线的参数方程利用三维曲线的参数方程绘制三维图形绘制三维图形 t=(0:0.02:2)*pi;%参数采样参数采样 x=sin(t);y=cos(t);z=cos(2*t); plot3(x,y,z,b-,x,y,z,bd) view(-82,58),box on xlabel(x),ylabel(y),zlabel(z) legend(链链,宝石宝石) 2.三维曲面三维曲面/网线图网线图三维图形较二维图形复杂：表现在：三维图形较二维图形复杂：表现在：数据的准备数据的准备三维图形色彩的使用三维图形色彩的使用明暗处理、明暗处理、光照处理、光照处理、1. 视点处理等视点处理等（1）数据准备（很重要）数

4、据准备（很重要） 画函数画函数 所代表的三维空间曲面，需所代表的三维空间曲面，需要做以下数据准备：要做以下数据准备： ： （1）确定自变量的取值范围和取值间隔；）确定自变量的取值范围和取值间隔； （2）构成自变量）构成自变量x,y的自变量的自变量“格点格点”矩阵；矩阵； （3）计算在自变量采样）计算在自变量采样“格点格点”上的函数上的函数值。值。),(yxfz （2）绘制曲面）绘制曲面/网线图的基本指令网线图的基本指令surf(Z) 以以Z矩阵列、行下标为轴自变量，矩阵列、行下标为轴自变量，画曲面图。画曲面图。surf(X,Y,Z) 最常用的曲面图调用格式最常用的曲面图调用格式surf(X,Y

5、,Z,C) 最完整调用格式最完整调用格式，画由，画由C指定用色的曲面图。指定用色的曲面图。mesh(Z) 以以Z矩阵列、行下标为轴自变矩阵列、行下标为轴自变量，画网线图。量，画网线图。mesh(X,Y,Z) 最常用的网线图调用格式最常用的网线图调用格式mesh(X,Y,Z,C) 最完整调用格式，最完整调用格式，画由画由C指定用色的网线图。指定用色的网线图。 说明说明对于最完整调用格式对于最完整调用格式surf(X,Y,Z,C) 和和mesh(X,Y,Z,C) ：四个输入宗量都是维数相同的矩阵；四个输入宗量都是维数相同的矩阵；X、Y是自变量是自变量“格点格点”矩阵；矩阵；Z是格点上的函数矩阵；是

6、格点上的函数矩阵；C是指定各点用色的矩阵。是指定各点用色的矩阵。C可以空缺，空缺时，可以空缺，空缺时，默认的用色矩阵默认的用色矩阵C=Z。单宗量输入绘图时，单宗量输入绘图时，Z矩阵的列下标当作矩阵的列下标当作x轴坐标的轴坐标的“自变量自变量”，把把Z的行下标当作的行下标当作y轴坐标的轴坐标的“自变量自变量”。例例5.3-2用曲面图表现函数用曲面图表现函数 clf x=-4:4;y=x; X,Y=meshgrid(x,y); Z=X.2+Y.2; surf(X,Y,Z); colormap(hot) hold on stem3(X,Y,Z,bo) hold off xlabel(x),ylabe

7、l(y),zlabel(z) axis(-5,5,-5,5,0,inf) view(-84,21) 22yxz单宗量输入时：单宗量输入时： clf x=-4:4;y=x; X,Y=meshgrid(x,y); Z=X.2+Y.2; surf(Z); colormap(hot) hold on mesh(X,Y,Z) hold off xlabel(x),ylabel(y),zlabel(z) view(-84,21) 3、颜色的使用、颜色的使用 但是很多时候，一个简单的二维或三维图但是很多时候，一个简单的二维或三维图形不能一次显示出想要提供的全部信息。形不能一次显示出想要提供的全部信息。 这时

8、，颜色可以对图形提供一个附加的维这时，颜色可以对图形提供一个附加的维数。数。 许多绘图函数都可以接受一个可用的颜色许多绘图函数都可以接受一个可用的颜色参量，来增加这附加的维数。参量，来增加这附加的维数。 本小节以研究颜色映象开始：本小节以研究颜色映象开始：如何使用、显示、修改和如何创建用户自如何使用、显示、修改和如何创建用户自己的颜色映象。己的颜色映象。然后，阐述在一个图形窗口中仿真多个颜然后，阐述在一个图形窗口中仿真多个颜色映象的技术或只使用颜色映象的一部分色映象的技术或只使用颜色映象的一部分的技术。的技术。最后，讨论照明模型并提供例子。最后，讨论照明模型并提供例子。 （1）颜色映像简介）颜

9、色映像简介 MATLAB有一个叫颜色映象的数据结构来有一个叫颜色映象的数据结构来代表颜色值。代表颜色值。 颜色映象定义为一个有三列和若干行的矩颜色映象定义为一个有三列和若干行的矩阵。阵。 利用利用0到到1之间的数，矩阵的每一行都代表之间的数，矩阵的每一行都代表了一种色彩。了一种色彩。 任一行的数字都指定了一个任一行的数字都指定了一个RGB值，即红、值，即红、黄、蓝三种颜色的强度，形成一种特定的黄、蓝三种颜色的强度，形成一种特定的颜色。颜色。 简单颜色映像表简单颜色映像表预定的颜色映象预定的颜色映象分别由十个分别由十个MATLAB函数产生函数产生 上面所列的各个颜色映象产生一个上面所列的各个颜色

10、映象产生一个643的的矩阵，指定了矩阵，指定了64种颜色种颜色RGB的描述。的描述。 这些函数都接受一个参量来指定所产生矩这些函数都接受一个参量来指定所产生矩阵的行数。比如阵的行数。比如hot(m)产生一个产生一个m3的矩的矩阵，它包含的阵，它包含的RGB颜色值的范围从黑经过颜色值的范围从黑经过红、橘红和黄，到白。红、橘红和黄，到白。 大多数计算机在一个大多数计算机在一个8位的硬件查色表中一次可位的硬件查色表中一次可以显示以显示256种颜色，当然有些计算机的显示卡可种颜色，当然有些计算机的显示卡可以同时显示更多的颜色。以同时显示更多的颜色。 这就意味着在不同的图中，一般一次可以用三或这就意味着

11、在不同的图中，一般一次可以用三或四个四个643的颜色映象。如果使用了更多的颜色的颜色映象。如果使用了更多的颜色映象输入项，计算机必须经常在它的硬件查色表映象输入项，计算机必须经常在它的硬件查色表中调出输入项。比如，当在画中调出输入项。比如，当在画MATLAB图形时背图形时背景图案发生了变化，就是发生了这种情况。景图案发生了变化，就是发生了这种情况。 所以，除非计算机有一次显示更多种颜色的显示所以，除非计算机有一次显示更多种颜色的显示卡，最好任何一次所用的颜色映象输入项数都小卡，最好任何一次所用的颜色映象输入项数都小于于256。 （2）颜色映象使用）颜色映象使用 语句语句colormap(M）将

12、矩阵）将矩阵M作为当前图形作为当前图形窗口所用的颜色映象。窗口所用的颜色映象。 例如，例如，colormap（cool）装入了一个有）装入了一个有64个输入项的个输入项的cool颜色映象。颜色映象。 colormap default装入了缺省的颜色映象装入了缺省的颜色映象（hsv）。）。 函数函数plot、plot3、contour（等高线图）（等高线图）和和contour3不使用颜色映象，它们使用列不使用颜色映象，它们使用列在在plot颜色和线形表中的颜色。而大多数其颜色和线形表中的颜色。而大多数其它绘图函数，比如它绘图函数，比如mesh、surf、fill、pcolor和它们的各种变形函数

13、，使用当前和它们的各种变形函数，使用当前的颜色映象。的颜色映象。 接受颜色参量的绘图函数中的颜色参量通常采用接受颜色参量的绘图函数中的颜色参量通常采用以下三种形式之一：以下三种形式之一： （1）字符串。代表）字符串。代表plot颜色或线型表中的一种颜颜色或线型表中的一种颜色，比如，色，比如， r 代表红色；代表红色； （2）三个输入的行向量。它代表一个单独的）三个输入的行向量。它代表一个单独的RGB值，比如值，比如.25 .50 .75； （3）矩阵。如果颜色参量是一个矩阵，其元素）矩阵。如果颜色参量是一个矩阵，其元素作了调整，并把它们用作当前颜色映象的下标。作了调整，并把它们用作当前颜色映象

14、的下标。最后一种形式会在以后作更多讨论。最后一种形式会在以后作更多讨论。 （3）颜色映象显示）颜色映象显示 可以用多种途径来显示一个颜色映象。其中一个方法是观可以用多种途径来显示一个颜色映象。其中一个方法是观察颜色映象矩阵的元素。察颜色映象矩阵的元素。 例如：例如： hot(8) ans =0.3333 0 00.6667 0 01.0000 0 01.0000 0.3333 01.0000 0.6667 01.0000 1.0000 01.0000 1.0000 0.50001.0000 1.0000 1.0000上面的数据显示出第一行是上面的数据显示出第一行是1/3红色，而最后一行是白色。

15、红色，而最后一行是白色。 另外，函数另外，函数pcolor可以用来显示一个颜色可以用来显示一个颜色映象映象。 n=16; colormap(jet(n) pcolor(1:n+1;1:n+1 ) title( Using Pcolor to Display a Color Map )用伪彩色来显示颜色映象用伪彩色来显示颜色映象因为上面这一段程序很有用处，它已经装入精通因为上面这一段程序很有用处，它已经装入精通MATLAB工具箱中的函数工具箱中的函数mmshow中。中。help mmshowMMSHOW PCOLOR Colormap DisplayMMSHOW uses pccolor to

16、display the current colormap.MMSHOW(MAP) displays the colormap MAP.MMSHOW(MAP(N) displays the colormap MAP having N elements.Examples:MMSHOW(hot)MMSHOW(pink(30)帮助信息：帮助信息：MMSHOW 显示显示PCOLOR颜色映象颜色映象MMSHOW 使用使用pcolor来显示当前颜色映象来显示当前颜色映象MMSHOW(MAP) 显示显示MAP颜色映象颜色映象MMSHOW（MAP(N) 显示一个有显示一个有N个元素的个元素的MAP颜色映象颜色映

17、象例子：例子：MMSHOW（hot） MMSHOW（pink（30） 函数函数mmshow取和取和colormap同样的输入参同样的输入参量，但在这种情况下它用自己的伪彩色显量，但在这种情况下它用自己的伪彩色显示而不是把颜色映象施加到当前图形。示而不是把颜色映象施加到当前图形。 另一种途径是使用另一种途径是使用MATLAB的函数的函数rgbplot，它可以把颜色映象的各列分别画成红、绿它可以把颜色映象的各列分别画成红、绿和蓝色。和蓝色。例如例如rgbplot(hot) 图中显示红色分量首先增加，然后是绿色，图中显示红色分量首先增加，然后是绿色，最后是蓝色。最后是蓝色。 rgbplot（gray

18、）表示所有三列数据均匀线）表示所有三列数据均匀线性地增加（三条线重叠）。性地增加（三条线重叠）。4、曲面、曲面/网线图的精细修饰网线图的精细修饰 （1）视角控制）视角控制viewview(az,el) 通过方位角、俯视角设置视点通过方位角、俯视角设置视点view(vx,vy,vz) 通过直角坐标设置视点通过直角坐标设置视点（2）色图）色图colormap colormap(CM) 设置当前图形窗的着色色图为设置当前图形窗的着色色图为CM MATLAB的预定义色图矩阵的预定义色图矩阵CMCM含含 义义CM含含 义义autumn红、黄浓淡色jet蓝头红尾饱和值色bone蓝色调浓淡色lines采用p

19、lot绘线色colorcube三浓淡多彩交错色pink淡粉红色图cool青、品红浓淡色prism光谱交错色copper纯铜色调线性浓淡色SPRING青、黄浓淡色flag红-白-蓝-黑交错色summer绿、黄浓淡色gray灰色调线性浓淡色winter蓝、绿浓淡色hot黑、红、黄、白浓淡色white全白色hsv两端为红的饱和值色注意：jet时默认色图。 函数函数colorbar在当前的图形窗口中增加水平或垂在当前的图形窗口中增加水平或垂直的颜色标尺以显示当前坐标轴的颜色映象。直的颜色标尺以显示当前坐标轴的颜色映象。 colorbar( horiz ) 在当前的图形下面放一个水在当前的图形下面放一个

20、水平的颜色条。平的颜色条。 colorbar( vert ) 在当前的图形右边放一个垂直在当前的图形右边放一个垂直的颜色条。的颜色条。 对无参量的对无参量的colorbar，如果当前没有颜色条就加，如果当前没有颜色条就加一个垂直的颜色条，或者更新现有的颜色条。一个垂直的颜色条，或者更新现有的颜色条。 x,y,z=peaks; mesh(x,y,z); colormap(hsv) axis(-3 3 -3 3 -6 8) colorbar（3）浓淡处理）浓淡处理shadingshading options图形对象着色的浓淡处理图形对象着色的浓淡处理mesh、surf、pcolor、fill和和f

21、ill3所创建的所创建的图形非数据点处的着色由图形非数据点处的着色由shading指令决指令决定。定。flatinterptfaceted1. shading是设置当前轴上是设置当前轴上“面面”对象的对象的Edgecolor和和Facecolor属性的高层指令。属性的高层指令。【例例5.3-3】三种浓淡处理方式比较。三种浓淡处理方式比较。 clf x=-4:4;y=x; X,Y=meshgrid(x,y); Z=X.2+Y.2; surf(X,Y,Z) colormap(jet) subplot(1,3,1),surf(Z),axis off subplot(1,3,2),surf(Z),ax

22、is off,shading flat subplot(1,3,3),surf(Z),axis off,shading interp set(gcf,Color,w) （4）透明度控制）透明度控制 alpha alpha(v) 对面、块、象三种图形对象的透明度加对面、块、象三种图形对象的透明度加以控制；以控制；0-1之间，之间，0为完全透明，为完全透明，1为不透明。为不透明。 【例例5.3-4】半透明的表面图。半透明的表面图。 clf surf(peaks) shading interp alpha(0.5) colormap(summer) （5）灯光设置）灯光设置lightlight(co

23、lor,option1,style,option2,position,option3)灯光设置灯光设置指令使用前后，原为强度相等的漫射光，后为面等子对指令使用前后，原为强度相等的漫射光，后为面等子对象有关光的属性都被激活。象有关光的属性都被激活。该指令可以不包含任何宗量，采用默认设置：白光、无该指令可以不包含任何宗量，采用默认设置：白光、无穷远、穿过穷远、穿过1 0 1射向坐标原点。射向坐标原点。任何一个宗量都可以缺省；任何一个宗量都可以缺省；option1可以采用可以采用RGB三元组或相应的色彩字符；三元组或相应的色彩字符；option2可以取可以取infinte和和local分别表示无穷远

24、光分别表示无穷远光和近光。和近光。option3为直角坐标三元组形式，对远光，表示光线穿为直角坐标三元组形式，对远光，表示光线穿过该点射向原点；对近光，表示光源所在位置。过该点射向原点；对近光，表示光源所在位置。（6）照明模式）照明模式lightinglighting options设置照明模式设置照明模式flat与与faceted配合使用，默认格式，入射配合使用，默认格式，入射光均匀洒在图形对象的每个面上。光均匀洒在图形对象的每个面上。gouraud先对顶点颜色插补，再对顶点勾先对顶点颜色插补，再对顶点勾画的面色进行插补，用于曲面表现；画的面色进行插补，用于曲面表现；phong对顶点处法线差

25、值，再计算各像素对顶点处法线差值，再计算各像素的反光。的反光。1. none使所有光源关闭。使所有光源关闭。（7）控制光反射的材质指令）控制光反射的材质指令materialmaterial options使用预定义反射模式使用预定义反射模式 shiny使对象比较明亮，镜反射份额较大，反使对象比较明亮，镜反射份额较大，反射光颜色仅取决于光源颜色。射光颜色仅取决于光源颜色。dull适用对象比较暗淡，漫反射份额较大，没有适用对象比较暗淡，漫反射份额较大，没有镜面亮点，反射光颜色仅取决于光源颜色。镜面亮点，反射光颜色仅取决于光源颜色。metal使对象带金属光泽，镜反射份额较大，背使对象带金属光泽，镜反

26、射份额较大，背景光和漫反射份额很小，反射光颜色取决于光景光和漫反射份额很小，反射光颜色取决于光源和图形表面两者的颜色。源和图形表面两者的颜色。1.Default返回默认设置模式。返回默认设置模式。【例例5.3-5】灯光、照明、材质指令灯光、照明、材质指令所表现的图形。所表现的图形。 clf; X,Y,Z=sphere(80); colormap(jet) subplot(1,2,1),surf(X,Y,Z),axis equal off,shading interp light (position,0 -10 1.5,style,infinite) lighting phong materia

27、l shiny subplot(1,2,2),surf(X,Y,Z,-Z),axis equal off,shading flat light;lighting flat light(position,-1,-1,-2,color,y) light(position,-1,0.5,1,style,local,color,w) set(gcf,Color,w)4、透视、镂空和裁切、透视、镂空和裁切 （1）图形的透视）图形的透视 hidden off透视被叠压的图形透视被叠压的图形 hidden on消隐被叠压的图形消隐被叠压的图形【例例5.3-6】透视演透视演 X0,Y0,Z0=sphere(3

28、0); X=2*X0;Y=2*Y0;Z=2*Z0; surf(X0,Y0,Z0); shading interp hold on,mesh(X,Y,Z),colormap(hot) hold off hidden off axis equal,axis off （2）图形的镂空）图形的镂空 【例例5.3-7】演示：如何利用演示：如何利用“非数非数”NaN，对图，对图形进行镂空处理。形进行镂空处理。 P=peaks(30); P(18:20,9:15)=NaN; surfc(P); colormap(hot) light(position,50,-10,5) material(0.9,0.9,0

29、.6,15,0.4) grid off,box on 镂空实例镂空实例2 clf; t=linspace(0,2*pi,100); r=1-exp(-t/2).*cos(4*t); % 旋转旋转母线母线X,Y,Z=cylinder(r,60); % 产生旋转柱面数据产生旋转柱面数据ii=find(X0&Y6|abs(Y)6); ZZ(ii)=zeros(size(ii); surf(X,Y,ZZ),shading interp;colormap(copper) light(position,0,-15,1);lighting phong material(0.8,0.8,0.5,10,0.5)

30、 5、高维可视化、高维可视化 （1）二维半图指令）二维半图指令pcolor, contour, contourf 【例例5.4-1】“二维半二维半”指令指令pcolor, contour, contourf的应用的应用 clf;clear;X,Y,Z=peaks(40);n=6;subplot(1,2,1),pcolor(X,Y,Z)shading interpzmax=max(max(Z);zmin=min(min(Z);caxis(zmin,zmax)colorbarhold onC=contour(X,Y,Z,n,k:);clabel(C)hold offsubplot(1,2,2)C,

31、h=contourf(X,Y,Z,n,k:);clabel(C,h)colormap(cool)set(gcf,Color,w) “二维半二维半”指令的演示指令的演示 （2）四维表现）四维表现人对自然界的理解和思维是多维的。人的人对自然界的理解和思维是多维的。人的感官不仅善于接受一维、二维、三维的几感官不仅善于接受一维、二维、三维的几何信息，而且对几何物体的运动，对颜色、何信息，而且对几何物体的运动，对颜色、声音、气味、触感等反应灵敏。从这意义声音、气味、触感等反应灵敏。从这意义上讲，上讲，MATLAB色彩控制、动画等指令为色彩控制、动画等指令为四维或更高维表现提供了手段。四维或更高维表现提供

32、了手段。 A.准四维表现准四维表现【例例5.4-2】用颜色表现函数的其他特征（如用颜色表现函数的其他特征（如梯度、曲率等）梯度、曲率等） clfx=3*pi*(-1:1/15:1);y=x;X,Y=meshgrid(x,y);R=sqrt(X.2+Y.2)+eps;Z=sin(R)./R;dzdx,dzdy=gradient(Z);dzdr=sqrt(dzdx.2+dzdy.2);dz2=del2(Z); subplot(1,2,1),surf(X,Y,Z,abs(dzdr)shading faceted;colorbar(SouthOutside)brighten(0.6);colormap

33、 hsvtitle(No. 1 surf(X,Y,Z,abs(dzdr)subplot(1,2,2);surf(X,Y,Z,abs(dz2)shading facetedcolorbar(NorthOutside)title(No. 2 surf(X,Y,Z,abs(dz2) B.切片图切片图 X,Y,Z=meshgrid(x,y,z) 由采样向量产生由采样向量产生三维自变量三维自变量“格点格点”数组数组 slice(X,Y,Z,V,sx,sy,sz) 三元函数切片图三元函数切片图 【例5.4-3】图形表现。 clf x,y,z = meshgrid(-2:.2:2,-2:.25:2,-2:.

34、16:2); v =x.*exp(-x.2-y.2-z.2); xs = -0.7,0.7; ys = 0; zs =0; slice(x,y,z,v,xs,ys,zs) colorbar shading interp colormap hsv xlabel(x),ylabel(y),zlabel(z) title(The color-to-v(x,y,z) mapping) view(-22,39) alpha(0.3) 222zyxxev（3）动态图形动态图形在在MATLAB的的“上层上层”图形指令中的彗图形指令中的彗星轨线指令、色图变幻指令、影片动画指星轨线指令、色图变幻指令、影片动画指

35、令，能很方便地使图形及色彩产生动态变令，能很方便地使图形及色彩产生动态变化效果。化效果。由于在由于在Notebook 和硬拷贝下，这种动和硬拷贝下，这种动态变化效果都无法表现，因此本节所有例态变化效果都无法表现，因此本节所有例题都不提供图形，而只给出有关指令。题都不提供图形，而只给出有关指令。当在当在MATLAB指令窗中运作这些指令后，指令窗中运作这些指令后，便可在图形窗中看到相应的动态图形。便可在图形窗中看到相应的动态图形。 A.彗星状轨迹图彗星状轨迹图comet(x,y,p) 二维彗星轨线二维彗星轨线comet3(x,y,z,p) 三维彗星轨线三维彗星轨线【例例 5.4-4】简单二维彗星示

36、例。简单二维彗星示例。 shg;n=2;t=n*pi*(0:0.000005:1);x=sin(t);y=cos(t);plot(x,y,g);axis squarehold on comet(x,y,0.0001)hold off B.色图的变幻色图的变幻spinmapspinmap(t)spinmap(inf)spinmap(t,inc) 【例例 5.4-5】指令指令spinmap的应用。的应用。 ezsurf(x*y,circ);shading flat;view(-18,28)C=summer;CC=C;flipud( C );colormap(CC)spinmap(30,4) C.影

37、片动画影片动画M(i)=getframemovie(M,k)【例5.4-6】三维图形的影片动画。clf;shg,x=3*pi*(-1:0.05:1);y=x;X,Y=meshgrid(x,y);R=sqrt(X.2+Y.2)+eps;Z=sin(R)./R;h=surf(X,Y,Z);colormap(jet);axis offn=12;mmm=moviein(n);for i=1:nrotate(h,0 0 1,25);mmm(:,i)=getframe;endmovie(mmm,5,10) D.实时动画实时动画 【例 5.4-7】制作红色小球沿一条带封闭路径的下旋螺线运动的实时动画（图5.

38、4-5）。（1）anim_zzy1.mfunction f=anim_zzy1(K,ki)% anim_zzy1.m演示红色小球沿一条封闭螺线运动的实时动画演示红色小球沿一条封闭螺线运动的实时动画% 仅演示实时动画的调用格式为仅演示实时动画的调用格式为 anim_zzy1(K)% 既演示实时动画又拍摄照片的调用格式为既演示实时动画又拍摄照片的调用格式为f=anim_zzy1(K,ki)% K红球运动的循环数（不小于红球运动的循环数（不小于1）% ki指定拍摄照片的瞬间，取指定拍摄照片的瞬间，取 1 到到 1034 间的任意整数。间的任意整数。% f存储拍摄的照片数据，可用存储拍摄的照片数据，可

39、用image(f.cdata)观察照片。观察照片。% t1=(0:1000)/1000*10*pi;x1=cos(t1);y1=sin(t1);z1=-t1;t2=(0:10)/10;x2=x1(end)*(1-t2);y2=y1(end)*(1-t2);z2=z1(end)*ones(size(x2);t3=t2;z3=(1-t3)*z1(end);x3=zeros(size(z3);y3=x3;t4=t2;x4=t4;y4=zeros(size(x4);z4=y4;x=x1 x2 x3 x4;y=y1 y2 y3 y4;z=z1 z2 z3 z4;plot3(x,y,z,Color,1,0

40、.6,0.4,LineWidth,2.5)axis off% h=line(xdata,x(1),ydata,y(1),zdata,z(1),Color,1 0 0,Marker, ., MarkerSize,40,EraseMode,xor);% n=length(x);i=2;j=1;while 1set(h,xdata,x(i),ydata,y(i),zdata,z(i); drawnow;% pause(0.0005)% i=i+1; if nargin=2 & nargout=1 if(i=ki&j=1);f=getframe(gcf);end% end if in i=1;j=j+1; if jK;break;end endend （2）f=anim_zzy1(2,450);（3）image(f.cdata),axis off