随机游动的Matlab小程序

    前一段时间需要做大量的蒙特卡罗仿真，做了一个随机游动的Matlab小程序。

close all;
clear all;
s=0.04; %max step size
x=rand-0.5; %x initial value
y=rand-0.5; %y initial value
hx=animatedline('Color','r','Marker','x');
xborder=[-1 1];yborder=[-1 1];  %x and y axis boundary
axis([xborder yborder]);
box on;
while 1
    dsx=x+s*randn;
    dsy=y+s*randn;
    if(dsx>xborder(2))  %dicide if x exceed the boundary
        x=xborder(2)*2-dsx;
    elseif(dsx<xborder(1))
        x=xborder(1)*2-dsx;
    else
        x=dsx;
    end
    if(dsy>yborder(2))  %dicide if y exceed the boundary
        y=yborder(2)*2-dsy;
    elseif(dsy<yborder(1))
        y=yborder(1)*2-dsy;
    else
        y=dsy;
    end
%     plot(x,y,'or');axis([xborder yborder]);
    addpoints(hx,x,y);
    drawnow update
    pause(0.05)
end

    下面是一个运行了约2分钟效果图。