例子:网格50*50,网格距为1,每个网格里有均匀分布的100个粒子,粒子电量为1,使用粒子云分室法计算格点的电荷量。 结果图 MATLAB程序 function =electron_density() q=zeros(51,51); for i=1:50 for j=1:50 for m=1:10 for n=1:10 q(i,j)=q(i,j)+(11-m)*(11-n)/(11*11); q(i+1,j)=q(i+1,j)+m*(11-n)/(11*11); q(i,j+1)=q(i,j+1)+(11-m)*n/(11*11); q(i+1,j+1)=q(i+1,j+1)+m*n/(11*11); end end end end raw=zeros(1,51);column=zeros(52,1) q_raw= ;q_column= figure(1);pcolor(q_column);colorbar; end