高尔顿钉板实验的matlab代码动画演示
代码中的n,T,delay等参数均可调,但是过大的时候可能画出来的图不太好看。演示效果大致如下:慢速的动画版代码n=10;T=10;delay=0.2;count=zeros(1,n);f1=figure;f2=figure;if mod(n,2)==1left=-floor(n/2);right=floor(n/2);elseleft=-n/2;right=n/2-1;endpause(1);f
·
代码中的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
更多推荐
已为社区贡献1条内容
所有评论(0)