Calton板的数值模拟

Calton板的数值模拟 

EmilMatthew(EmilMatthew@126.com)

摘要: 

       本文基于Matlab，实现了Calton板的数值模拟。

关键词: Calton板，二项分布

 

Numerical Simulation of Calton Board

                             EmilMatthew(EmilMatthew@126.com)

Abstract: 

       In this article , I implemented the numerical simulation of Calton Board based on the tool of matlab.

Key Words: Calton Board 

      

1前言:

       Calton板是一个常用来在物理的热力学中说明粒子分布现象的实验。

       对于N个粒子分布在ε1～εＭ上共M个能级上会有多种分布，其中热力学概率最大的分布称为最概然分布。

　　又对于热力学系统N≥10２４，N，V，E确定的系统达平衡时(即系统的热力学态)，粒子的分布方式几乎将不随时间而变化，这种分布称为平衡分布。可以证明，平衡分布即为最概然分布所能代表的那些分布。

而最终Calton板最概然分布的极限分布接近于二项分布[1]，这就为数值仿真带来了方便。

      

2实现原理:

由于Calton的极限分布接近于一个二项分布，所以，在小球即到每个钉子时，可以根据概率p来选择是向左转还是向右转。（这里，采用以概率p向左转）。假设有n排钉子，则最终小球可以落入的位置为n+1列（见图1）。有了以上认识后，进行相应的数值模拟是方便的。                                             

                                                 

                                                                图1

 

       3实验代码:

%----simulate of calton board----

 

%--global variable--

 

% the possibility of left drop direction.

  p=0.7;

% the level that the ball should pass

  level=7;

% the number of the ball generated which in one trial.

  n_balls=10000;

% array to store the num in each tunnel.

  a=zeros(1,level+1);

  aStatics=zeros(1,n_balls);

% the tunnel looks from left.

  leftPos=1;

 

% iter val

  i=0;

  iBallCount=0;

  tmpRnd=0; 

 

%--core simulatioin--

  %total ball simulate

  while iBallCount<n_balls

              %--each ball selection--

              i=0;leftPos=1;

              %drop levels selection

              while i<level

                    tmpRnd=rand(1,1);

                    %select forward direction.

                    if(tmpRnd>=p)

                         leftPos=leftPos+1;

                    end     

                    i=i+1; 

              end

              %note each simulation ball's last pos.

              a(leftPos)=a(leftPos)+1;

           aStatics(iBallCount+1)=leftPos;

          

        iBallCount=iBallCount+1;

  end

 

  hist(aStatics,level+1)

 

       4实验结果:

       图2及图3展示了测试粒子为10000粒子时，左转概率p分别为0.5及0.7的分布情况，可以看出，结果很好的验证了Calton板的极限分布服从二项分布这一推论。                                    

图2

图3

参考文献: 

[1]基于Matlab的Calton板实验仿真设计及实现，计算机仿真,2006.3

                                                                                                         程序完成日: 06/04/30

                                                                                                         文章完成日: 06/05/01

 

附录:

1测试程序下载:

http://emilmatthew.51.net/EmilPapers/06_15Calton/code.rar

若直接点击无法下载(或浏览)，请将下载(或浏览)的超链接粘接至浏览器地址栏后按回车.若不出意外，此时应能下载。

若下载中出现了问题，请参考:

http://blog.csdn.net/emilmatthew/archive/2006/04/08/655612.aspx