MATLAB绘制粒子沿平面上任何角度的布朗运动

clear all;
close all;
z=0;%设置初始值，这里利用一个复数来表示位置
rand('state',6);%初始化随机数的状态
ph=plot(z);%绘图
hold on;
T=1;%计算时间的参数
z1=1;%位置的中间变量
ti=title('time=0','FontSize',12);%在图题处实时显示时间
while abs(z1)>0.1; %循环多次进行计算
    A = rand*pi*2;%生成一个随机角度
    z1 = z+exp(i*A);%更新位置数值
    plot(real([z,z1]),imag([z,z1]));%连接t时刻和t+1时刻两点绘制线段
    z=z1;%更新z的数值
    set(ti,'String',['time',num2str(T)]);%更新图题处的时间数值
    T=T+1;%时间值更新
    pause(0.1);%暂停一下，显示动画效果
end

详解MATLAB图形绘制技术—228页张德丰，雷小平主编；周燕副主编；［Ｍ］北京：电子工业出版社