matlab画动画

一般来说，matlab制作动画有四种方式。

第一 、以质点运动轨迹的方式显示

使用comet、comet3函数，前者是二维，后者是三维

comet（y）显示质点绕向量y，comet（x，y）显示质点绕向量y与x，comet(x,y,p)，其中为轨迹尾巴的长度

以comet（x，y）为例，

显示平抛运动

vx = 40;

t = 0:0.001:10;

x = vx*t;

y = -9.8*t.^2/2;

comet(x,y)

显示导弹发射

vx = 100*cos(1/4*pi);

vy = 100*sin(1/4*pi);

t = 0:0.001:15;

x = vx*t;

y = vy*t-9.8*t.^2/2;

comet(x,y)

匀速圆周运动

sita = 0:0.0001:2*pi;

r = 10;

x=r*cos(sita);

y=r*sin(sita);

comet(x,y)

comet3与comet的用法相类似，可以在帮助文件里的例子

t = -10*pi:pi/250:10*pi;

comet3((cos(2*t).^2).*sin(t),(sin(2*t).^2).*cos(t),t)

第二、以电影播放的方式显示

保存想要产生动画的图片，存储为一系列各种类型的二维、三维图，再像放电影的方式按次序播放出来。步骤由getframe函数将当前的图片抓取为电影的画面，再由movie函数将动画显示出来。

如：

[x,y] = meshgrid([-1.05:.2:3.75]);

z = x.*exp(-x.^2-y.^2);

axis tight;

set(gca,\'nextplot\',\'replacechildren\');

for j = 1:40

surf(x*sin(pi*j/100),y*sin(pi*j/100),z*sin(-pi*j/100));

m(j) = getframe

end

movie(m)

第三、以对象方式显示

设置对象的属性EraseMode，更新对象来产生新图，drawnow（）函数进而覆盖旧图，从而使得图形不断发生变化。

例：

x = -pi:pi/30:pi;

h = plot(x,cos(x),\'o\',\'MarkerEdgeColor\',\'k\',\'MarkerFaceColor\',\'r\',\'MarkerSize\',8,\'EraseMode\',\'Xor\')

for j = 1:10000

y = 1/2*sin(3*x+0.006*j);

set(h,\'ydata\',y);

drawnow;

end

第四、以旋转颜色的方式显示

matlab中如何输出avi格式影片文件

以上文《穿越（fly through)图形动画实例》为基础，将屏幕上显示的动画保存到avi格式文件中。

clc; %清屏

clear all; %清除所有变量

load wind; %加载matlab自带数据程序

wind_speed = sqrt(u.*u + v.*v + w.*w);

%计算风速

hpatch=patch(isosurface(x,y,z,wind_speed,35)); %isosurface计算wind_speed=45的等势面

%patch命令绘制等势面

isonormals(x,y,z,wind_speed,hpatch); %isonormals计算等势面的正交向量，使画面更光滑

set(hpatch,\'FaceColor\',\'red\',\'EdgeColor\',\'none\'); %设定等势面的表面颜色为红色

daspect([1 1 1]); �ta aspect

[f verts]=reducepatch(isosurface(x,y,z,wind_speed,45),.05);

% reducepatch ,减少需要描绘等势面所需要表面数， 减少至 5%

hcone=coneplot(x,y,z,u,v,w,verts(:,1),verts(:,2),verts(:,3),2); %绘制穿越等势表面的圆锥流体

set(hcone,\'FaceColor\',\'blue\',\'EdgeColor\',\'none\');

camproj perspective %设置该图形的投影方式为透视

camva(25); %设置视角为25度

hlight = camlight(\'headlight\'); %在照相机出创建光源.

set(hpatch,\'AmbientStrength\',.1,... %设置环境光的光度为0.1，表面环境光近乎为黑色

\'SpecularStrength\',1,... %设置块对象的镜面反射强度为1

\'DiffuseStrength\',1); %设置快对象的漫反射强度为1

set(hcone,\'SpecularStrength\',1); %设置锥形流体的镜面反射强度为1

set(gcf,\'Color\',\'k\'); %将图像窗的颜色设置为黑色

lighting phong %sets the lighting to phong.

set(gcf,\'Renderer\',\'zbuffer\'); %设置光源的着色属性为zbuffer，也可以尝试openGL

hsline=streamline(x,y,z,u,v,w,80,30,11); %创建一个从点（80，30，11）开始的流线

xd = get(hsline,\'Xdata\'); %取得这条流线的x,y,z坐标

yd = get(hsline,\'Ydata\');

zd = get(hsline,\'Zdata\');

delete(hsline); %流线已无用，可删除

%以下为此动画的核心部分

aviobj=avifile(\'fly_through.avi\',\'fps\',3);

%创建avi对象，命名为fly_through.avi, 3帧每秒

for i=1:length(xd)-50 %从点（80，30，11）开始，到终点前的50个点出结束。

campos([xd(i),yd(i),zd(i)]); %设定照相机的位置为当前点位置

camtarget([xd(i+5),yd(i),zd(i)]); %设置照相机的目标位置，在照相机位置向前5个点出

camlight(hlight,\'headlight\'); %光源也移至照相机处

drawnow %绘制移动后的图形

frame=getframe(gcf); %获取当前frame

aviobj=addframe(aviobj,frame); %将当前图像加入到avi对象中

end

aviobj=close(aviobj); %结束时关闭avi对象

%注意，此方法，运行时不能移动或改变动画窗口的大小。

function makeavi(filepath,frames,format)

% This function make AVI video from images

% filepath: the location of target image

% frames: the numbles of frame that you want to use

% fromat: the format of images

%================设置相关属性======================

location=filepath;

fig=figure;

set(fig,\'DoubleBuffer\',\'on\');

set(gca,\'xlim\',[-80 80],\'ylim\',[-80 80],...

\'NextPlot\',\'replace\',\'Visible\',\'off\')

aviObj=avifile(\'examples.avi\'); % Create a new AVI file

aviObj.Quality = 100; %只对压缩格式有效，[0,100],default:75

aviObj.fps=1; %the speed of the AVI movie in frames per second (fps),default:15 fps

aviObj.Compression=\'None\'; % compression codec指定压缩编解码器

%=======设置黑白图像属性===========

%cola=0:1/255:1;

%cola=[cola;cola;cola];

%cola=cola\';

%aviobj.colormap=cola;

% =====多帧循环开始,读序列图像=================

for num=1:frames

Im_origin = imread( strcat(location,\'\\',num2str(num),\'.\',format) );

Im= im2uint8(Im_origin);

% F = getframe(gca);

aviObj = addframe(aviObj,uint8(Im));

end

aviObj=close(aviObj);

mov = aviread(\'examples.avi\');

movie(mov);

close all;

转载自：http://blog.sina.com.cn/s/blog_6b7d710b0101m39u.html