anim_zzy1.m

学matlab还有比这个更好的吗?我千兴万苦才找到的,这是一个学习matlab的视频,功能强大

function f=anim_zzy1(K,ki)
% anim_zzy1.m		演示红色小球沿一条封闭螺线运动的实时动画
% 仅演示实时动画的调用格式为  					anim_zzy1(K)			
% 既演示实时动画又拍摄照片的调用格式为		f=anim_zzy1(K,ki)		
% K						红球运动的循环数（不小于1）
% ki				指定拍摄照片的瞬间，取 1 到 1034 间的任意整数。
% f						存储拍摄的照片数据，可用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,'b'),	axis off		% 绘制曲线
% 定义"线"色、"点"型（点）、点的大小（40）、擦除方式（xor)
h=line('Color',[1 0 0],'Marker','.','MarkerSize',40,'EraseMode','xor');
% 使小球运动
n=length(x);i=1;j=1;
while 1													% 无穷循环
	set(h,'xdata',x(i),'ydata',y(i),'zdata',z(i));	% 小球位置
   drawnow;												% 刷新屏幕				<21>
   pause(0.0005)											% 控制球速				<22>
   i=i+1;
   if nargin==2 & nargout==1	% 仅当输入宗量为2、输出宗量为1时，才拍摄照片
      if(i==ki&j==1);f=getframe(gcf);end				% 拍摄i=ki时的照片	<25>
   end
   if i>n 
      i=1;j=j+1;
      if j>K;break;end
   end
end