高尔顿钉板实验的matlab代码动画演示

代码中的n，T，delay等参数均可调，但是过大的时候可能画出来的图不太好看。演示效果大致如下：

慢速的动画版代码

n=10;T=10;delay=0.2;
count=zeros(1,n);
f1=figure;f2=figure;
if mod(n,2)==1 
    left=-floor(n/2);right=floor(n/2);
else
    left=-n/2;right=n/2-1;
end
pause(1);
for k=1:T
    figure(f1);
    clf;
    for i=1:n
      text(-i+n/2+0.5,0.5,int2str(count(i)),'Color','red','FontSize',14);%计数
    end
    x=[];y=[];
    for i=n:-1:1
        x=[x,(-i:2:i)/2];
        y=[y,zeros(1,i+1)+n+1-i];%计算钉子坐标
    end
    figure(f1);hold on;
    scatter(x,y+0.15,'ok','filled');%画钉子
    r=round(rand(1,n-1)).*2-1;
    xx=0;
    for i=n-1:-1:1
        xx=xx+r(i);
        figure(f1);
        scatter(xx/2,i+0.5,160,'oy','filled');%画乒乓球
        pause(delay);
    end
    text(xx/2,0.5,int2str(count(-xx/2+(n+1)/2)),'Color','white','FontSize',14);%擦除原有的数字
    count(-xx/2+(n+1)/2)=count(-xx/2+(n+1)/2)+1; %该位置自增
    text(xx/2,0.5,int2str(count(-xx/2+(n+1)/2)),'Color','red','FontSize',14);%写更新的数字
    figure(f2);
    bar(left:right,count(n:-1:1));%更新直方图
    text(left,max(count),"t="+k,'Color','black','FontSize',14);
    pause(delay*2);
end
figure(f1);
text(left,n-1,"n="+n+" T="+T,'Color','black','FontSize',14);
hold off

快速的计算演示代码

n=15;T=10000;
count=zeros(1,n);
x=[];y=[];
for i=n:-1:1
    x=[x,(-i:2:i)/2];
    y=[y,zeros(1,i+1)+n+1-i];%计算钉子坐标
end
for k=1:T
    r=round(rand(1,n-1)).*2-1;
    xx=sum(r);
    count(-xx/2+(n+1)/2)=count(-xx/2+(n+1)/2)+1; %该位置自增
end
if n<30
    f1=figure;figure(f1);
    scatter(x,y+0.15,'ok','filled');%画钉子
    for i=1:n
      text(-i+n/2+0.1,0.5,int2str(count(i)),'Color','red','FontSize',8);%显示计数
    end
    text(-(n-1)/2,n,"n="+n+" T="+T,'Color','black','FontSize',14);%显示测试参数
end
f2=figure;figure(f2);
bar(-(n-1)/2:(n-1)/2,count(n:-1:1));%画直方图
text(-(n-1)/2,max(count),"n="+n+" T="+T,'Color','black','FontSize',14);%显示测试参数
hold off