电子双缝衍射实验的MATLAB演示我的毕业论文.doc

1、电子双缝衍射实验的 MATLAB 演示xxx （xxxxxx 理学院 09 级物理学 1 班 云南 玉溪 653100）指导老师：xxxx摘要：依据电子衍射的几率密度函数，运用蒙特卡洛随机模拟方法，借助 MATLAB编程和计算机可视化技术、绘图技术、构建电子双缝衍射的动态随机过程、清晰的展示电子双缝衍射的全过程。关键字：电子双缝衍射；蒙特卡洛；几率密度函数；MATLAB;计算机模拟；波函数1.引言通过光的衍射现象和实验来进一步揭示光的波动性，说明衍射现象是光在空间或物质中传播的基本方式，同时说明衍射现象是讨论现代光学问题的基础 1 。1924 年，法国物理学家德布罗意（L.V.de Brogl

2、ie)提出物质波概念，并指出可以用晶体对电子的衍射实验加以证明。1961 约恩孙（C.JOnsson)得到了电子单缝、双缝衍射，说明电子具有波动性 2。M.Born(1926)提出把粒子性与波动性统一起来，更确切地说，把微观粒子的“原子性”与波的“相干叠加性”统一起来。他是在薛定谔方程来处理散射问题时为了解释散射粒子的角度分布提出来的。他认为量子力学中的波函数所描述的，并不像经典波那样代表什么实在的物理波动性，只不过是刻画粒子在空间的概率分布的概率波（probability wave） 3。电子双缝衍射 4-7对实验仪器的稳定性要求很高，实验平台要求防震，环境的温度湿度都对实验有一定的影响。由

3、美国 Mathworks公司推出的 MATLAB8-11 是数值分析中较强的软件，具有强大的数值计算、数据可视化与编程、仿真和图形处理等功能，克服普通实验不利的因素，从而直观的展示电子双缝衍射的相关图像，加深对原理概念的理解。对电子双缝衍射几率密度的计算推导，采用 MATLAB进行编程，分别画出 e_number为 100、1000、10000、20000 个电子的双缝衍射图像与实验得出的进行比对，总结分析比对结果。2.电子双缝衍射几率密度函数2.1 电子双缝衍射几率密度函数的理论根据量子力学原理，任何微观粒子都具有波粒二象性。微观粒子的波动性并不是由于粒子大量聚集在空间中形成的，即波动性是单

4、个微观粒子所具有的特性。若已知微观粒子的波函数 12,13，则其模的平方 表示的是微观)(x2)(x粒子出现在 z点附近的几率大小，与此相对应的就是波动光学中所说的光波强度，这被称为波函数的统计解释。设入射电子流很微弱，电子几乎是一个一个经过双缝然后在感光底片上被记录下来。起初，当感光时间较短时，底片上出现一些点子，它们看起来没什么规律。当感光时间足够长时，底片上出现的点子愈来愈多，就会发现有的地方点子很密，有的地方则几乎没有点子，最后在底片上呈现出有律明亮布的条文 1。2.2 电子双缝衍射几率密度函数的计算设质量为 m的电子经电场加速后具有速度 v，根据德布罗意的物质波理论，电子的动量和波长

5、的关系为：mhp 1当 ( 为光速)时，电子波的波长为：cveUh2 2式中U为加速电压。按照波动理论，电子穿过双缝应产生衍射现象。图1是电子双缝衍射原理图。图 1 电子双缝衍射原理图电子的双缝衍射强度分布函数为：sincossinii220 baaIx 3式中， 为 处的衍射强度， 为衍射角， 为缝宽， 为双缝间距。当30I时( 为缝与屏的间距 )下式成立：DxDxsin 4把式 代入式 得：43xbaxDaIx220cossi 5由式 知，衍射强度最大值在 处，且 。令5 0maxIDbxIK220cossin 6则 为 处相对于 的衍射强度，显然 。x0I1,0K图l是电子双缝衍射的示意

6、图。衍射屏位于 平面，观测屏位于 平yoxxoy面，缝 和 的宽度均为 ，两缝中心间距为 ，衍射屏到观测屏的距离1s2aba为 。设动量为 、能量为E的自由电子沿 轴正方向入射到双缝，其波函数为：DpZ)(Etzioe 7设在 时刻，波前到达双缝处 ，由式 可知，在双缝处波函数0t 0Z7为一常数。设在衍射屏上的波函数为：8根据费曼的路径积分理论，电子经过双缝而在t时刻到达观测屏上P点的衍射波函数为： Etziebakaktx sn2cosin2),(0 9式中， ，所以，电子经过双缝在观测屏上P点出现的几率密度为：9/2k sin2cossin2ii4),() 2202 bakatxsin2

7、cossinii20 baa 10设 ，则有：baDDxsi 1代入 有：10 BxAxxbaxDax 220220 cossincossin)( 其中， ， 。AabB2.3 电子双缝衍射动态过程的实现根据几率密度函数分布的随机数系列，采用蒙特卡罗随机抽样的 Von Neumannn的舍选发。舍选发的具体内容为：计算机在一定的范围内随机地选取观测屏一坐标点 ，并计算1 iyx,，其中， 是几率密度函数 在点 的值， 是几max/iHiximax率密度函数 的最大值；计算机产生一个0到1之间均匀分布的随机数 ； 2 M将 与 进行比较，若 ，则选取该点 ，若 ,则舍去该3MHiyx,H点，并重

8、复 至 ，重新选择坐标点。通过以上方法，即可以得到一个按13分布的随机数序列。这个随机数序列就是实现电子衍射随机运动过程的依x据。3 电子双缝衍射实验 MATLAB 的可视化3.1 电子双缝衍射的 MATLAB 程序 h=6.62559e-34;m=9.10908e-31;e=1.6021e-19;U=1000;a=2e-7;b=1e-6;D=0.25; rd=h/sqrt(2*m*e*U); %计算电子波长 rd A=(pi*a)/(rd*D);B=pi*(a+b)/(rd*D); %计算 A和 B axis(-5e-5,5e-5,-4e-5,4e-5); %设置坐标轴范围 set(gca,

9、color,0.122,0.012,0.62); %设置背景色 i=1; while i=Mhold on %保留当前坐标系中已经存在的图形对象h_point=plot(x,y,.r,EraseMode,none); %用红色显示符合条件的坐标点，并赋句柄值给 h_pointi=i+1; pause(1) %暂停 1秒set(h_point,color,w) %将当前坐标点由红色变为白色endend注释：e_number 只需更改电子数目几率密度函数： h=6.62559e-34;m=9.10908e-31;e=1.6021e-19;U=1000;a=2e-7;b=1e-6;D=0.25; r

10、d=h/sqrt(2*m*e*U); %计算电子波长 rd A=(pi*a)/(rd*D);B=pi*(a+b)/(rd*D); %计算 A和 B axis(-5e-5,5e-5,-4e-5,4e-5) x=-3e-5:1e-8:3e-5; W=(sin(A*x).2)./(A*x).2).*(cos(B*x).2); %衍射几率分布函数 Warning: Divide by zero. plot(x,W,EraseMode,none); %将当前坐标点由红色变为白色 3.2 电子双缝衍射图像通过更改 e_number 分别为 100、1000、10000、20000 在 MATLAB上画出对

11、应的图像：e_number=100 的电子双缝衍射图像e_number=1000 的电子双缝衍射图像e_number=10000 的电子双缝衍射图像e_number=20000 的电子双缝衍射图像几率密度函数图像4.结论通过改变 e_number，对比 4个图像，变化十分明显。在普通的光学实验中，由于对仪器精度、实验环境的的要求非常高，在改实验变参数时引起条文变化不明显，难以观察，难以充分展示实验的特征。电子双缝衍射图像与几率密度函数图像完全吻合。致谢我衷心地感谢任 xxxx在我写论文中给予的认真指导和提供的诸多帮助。参考文献1 姚启均.光学教程 .（第四版）出版时间：2008年6月2 普通物

12、理实验3光学部分.（第四版）出版时间：2007（12）3 曾谨言.量子力学教程.（第二版）出版时间： 2003年2月4 罗礼进. 基于蒙特卡罗的电子双缝衍射的计算机模拟J. 南通大学学报(自然科学版). 2005(04)5 蒋练军,金辉霞. 电子双缝及多缝衍射的量子分析J. 云南大学学报(自然科学版). 2005(S2)6 王筠,吉紫鹃,肖鹏程. 杨氏双缝实验中的干涉与衍射现象J. 湖北第师范学院学报. 2012(08)7 林连昶. 谈双缝干涉与双缝衍射J. 宁波教育学院学报. 2001(04)8 张峥、杨文平、石博强、李海鹏. MATLAB 程序设计与实例应用.出版时间：2003（11）9

13、张亮、郭仕剑、王宝顺、贺兴华.MTHLAB7 .X系统建模与仿真.出版时间：2006（11）10 贺兴华、周媛媛、王继阳、周辉.MTHLAB7 .X 系统建模与仿真.出版时间：2006（11）11 陈凤翔,汪礼胜,谭改娟,许文英. 量子力学虚拟实验的 MATLAB演示J. 大学物理. 2012(05)12 张德伟. 浅谈量子力学中的概率波J. 科教文汇(中旬刊). 2009(11)13 张树林. 基于 MATLAB的能级、波函数及几率密度图形的绘制J. 大学物理实验. 2008(03)Electronic double-slit diffraction experiment MATLAB de

14、moxxx(Yuxi Teachers College of Science 09 Physics classes Yuxi 653100)Directed by Prof. xxxxxAbstract: Based on the probability density function of the electron diffraction, the use of stochastic simulation Monte Carlo method, with MATLAB programming and computer visualization technology, graphics technology, build dynamic random process electronic double slit diffraction clear display electronic double-slit diffraction process.Keywords: double-slit diffraction; Monte Carlo; probability density function; MATLAB; computer simulation, the wave function.