资源描述：
《matlab在量子力学中地的应用》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、《MATLAB语言》课程论文MATLAB在量子力学中的应用姓名：魏祎学号：12011243989专业：通信工程班级：2011级通信工程指导老师：汤全武学院：物理电气信息学院完成日期：2012-12-119MATLAB在量子力学中的应用（魏祎120112439892011级通信班）[摘要]量子力学的应用和成就是多方面的，迄今仍保持有旺盛的生命力，硕果颇传。虽然《大学物理》中介绍的量子力学只是一些最基本的概念，但之中涉及了许多复杂的数值计算问题，解微分方程的问题，图像显示问题，例如一维无限深势阱问题，一维运动粒子的波函数曲线问题，对其手工求解较为复杂，而MATLA
2、B语言正是处理这些复杂问题的很好工具，既能进行数值求解，又能绘制有关曲线，非常方便实用，另外利用其可以减少工作量，节约时间，加深理解对量子力学的理解，同时可以培养应用知识的能力。[关键词]量子力学MATLAB语言一维无限深方势阱波函数概率密度一、问题的提出MATLAB语言是当今国际上科学界（尤其是自动控制理论）最具影响力、也是最有活力的软件。它提供了强大的科学运算、灵活的程序设计流程、高质量的图形可视化与界面设计、便捷的与其他程序和语言接口的功能。MATLAB语言在各国高校与研究单位起着重大的作用，它是一种集数值运算、符号运算、可视化建模、仿真和图形处理等多种
3、功能。在量子力学中，可以利用其帮助初学者理解量子力学与经典力学截然不同的思维方式和观念，理解微观粒子的波粒二象性。下面将以一维势阱问题，波函数和概率密度曲线问题为例讲述MATLAB在量子力学中的应用。二、用MATLAB语言求解一维无限深势阱问题如图1所示，设想一粒子处在势能为的力场中，并沿x轴作一维运动。粒子的势能满足下述边界条件：(1)当粒子在范围内时，=0；(2)当及时，。这就是说粒子只能在宽度为a的两个无限高势垒壁之间自由运动，就像一小球被限制在无限深的平底深谷中运动那样，我们理想化了得势阱曲线叫无限深方形势阱。因为粒子只限于沿x轴方向运动，故这个势阱为
4、一维无限深的方形势阱，简称一维方势阱。有上述边界条件已知，粒子在势阱中得势能（x）与时间无关，且=0.因此，由一般的薛定谔方程（1），粒子在无限深方势阱中得定态薛定谔方程为（1）式中m为粒子的质量，E为粒子的总能量。如令k为（2）9则上式可写成（3)根据边界条件，x=0时，，则可以利用MATLAB求解微分方程。利用MATLAB语言求解此方程程序如下：y=dsolve('D2y+k^2*y','y(0)=0','x')%%求方程（16）0ax图1一维无限深方势阱中得粒子运行结果：y=C1*sin(k*x)又根据边界条件，x=a时，，此时式（16）的解为（4）一般
5、说来，A可不为零，故，有n=1,2,3…，上式也可写成将上式与式（4）相比较，可得势中粒子可能的能量值为（5）式中n为量子数，表明粒子的能量只能取离散的值，当n=19时，势阱中粒子的能量为，n=2,3,4…时，4，9，16。这就是说，一维无限深方势阱中粒子的能量是量子化的。下面在来确定常数，由于粒子被限制在和的势阱中，因此，按归一化条件，粒子在此区间内出现的概率总和为1，即（6）或（7）令,则上式左侧积分为（8）于是，可得这样，式（7）所表现得波函数即为，（9）由此可得，能量为E所表示的粒子在势阱中得概率密度为下面用MATLAB语言求解电子的各能级能量、波函数
6、曲线和概率密度曲线。程序如下：functionE=shor(m,a)%建立函数文件n=1:10;%量子数nh=6.63*1e-34;%普朗克常量E=n.^2*h^2/(8*m*a^2);%n能级的能量值x=0:1.0*1e-12:a;%x的值subplot(4,2,1)%分割绘图区域，第一个子图y1=sqrt(2/a)*sin(pi*x/a);%n=1的波函数plot(x,y1);%绘制n=1的波函数图象title('n=1');%给n=1的破函数曲线加标题subplot(4,2,3)%分割绘图区域，第三个子图y1=sqrt(2/a)*sin(2*pi*x/a
7、);%n=2的波函数plot(x,y1);%绘制n=2的波函数图象9title('n=2');%给n=2的破函数曲线加标题subplot(4,2,5)%分割绘图区域，第五个子图y1=sqrt(2/a)*sin(3*pi*x/a);%n=3的波函数plot(x,y1);%绘制n=3的波函数图象title('n=3');%给n=3的破函数曲线加标题subplot(4,2,7)%分割绘图区域，第七个子图y1=sqrt(2/a)*sin(4*pi*x/a);%n=4的波函数plot(x,y1)%绘制n=4的波函数图象title('n=4');%给n=4的波函数曲线加标
8、题subplot(4,2,2)%分割绘
