MATLAB模拟伽尔顿板实验代码展示: clear;b=1;k=1;dxdy=0.81;r=0.1;ttt=0:pi/300:2*pi;xxx=r*exp(i*ttt);for n=1:80for m=10:20if xor(mod(m,2),mod(n,2))xy0(1,k)=dxdy*n;xy0(2,k)=dxdy*m;ax=xy0(1,k);ay=xy0(2,k);plot(xxx+-1+i*(-2),'r.'), plot(ax,ay,'r.','markersize',15) ,hol