以下m文件可以绘制小球沿着螺旋线运动的动画，其实就是两个步骤：
- 画出螺旋线背景
- 改变小球的位置，不断刷新图像，看起来就是动画了
完整的m文件代码如下
function f = ballw( K,ki )
%ballw.m演示红色小球沿着一条封闭螺旋线运动的实时动画
%仅演示实时动画的调试格式为ballw(K)
%既演示实时动画又拍摄照片的调试格式为f = ballw(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];
%data = [x',y',z'] %查看封闭曲线的坐标数值
plot3(x,y,z,'r','Linewidth',4) %绘制曲线
axis off; %不画坐标轴
%定义“线”色、“点”型（点）、点的大小（40）、擦除方式（xor)
h = line('Color',[0.67 0 1],'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));
%bw = [x(i),y(i),z(i)] %查看小球位置
drawnow; %刷新屏幕
pause(0.01) %控制球速
i = i+1;
if nargin == 2 && nargout == 1 %当输入变量个数为2并且输出变量1个时才拍摄照片
if (i == ki && j == 1)
f = getframe(gcf); %拍摄i = ki时的照片
end
end
if i > n
i = 1;
j = j+1;
if j > K
break;
end
end
end
代码注解：
- getframe(gcf) 获取整个窗口内容的图像，gcf是get current figure
- nargin是number of input arguments的缩写，意为输入参数的个数
- nargout是number of output arguments的缩写，意为输出参数的个数
- drawnow 刷新屏幕，当需要实时看到程序绘图变化，需要此函数
- pause(n) 暂停程序运行n秒
代码中注释掉的两行代码可以用来打印出坐标点，方便调试。
由此引出Matlab的m文件的调试方法：
-
直接调试法
- 打印重点怀疑命令行的结果（去掉行尾的;），或打印重点变量
- 在适当位置添加keyboard指令，中断程序，查看内存空间
- 在原函数文件首行之前加上百分号，使一个中间变量难于观察的M函数文件变为一个所有变量都保留在基空间中的M文件
- 用echo指令，使M文件在运行时在屏幕上逐行显示文件内容
-
工具调试法
- Matlab集成专门的调试工具，常用的就是设置断点，查看内存空间，跟踪变量值的变化。
可见Matlab的m文件调试方法和其他语言完全类似。
将上述的m文件保存为ballw.m，然后在Matlab命令窗口中执行如下代码：
>> f = ballw(1,200)
f =
cdata: [420x560x3 uint8]
colormap: []
>> image(f.cdata) %查看获得的照片
拍摄的小球照片
本来想将小球运动的gif动图导出贴在这儿的，可惜我不会，也懒得弄了，感兴趣的朋友可以学习一下，matlab的什么函数可以导出gif或者avi的动画。
