基于量子力学谐振子-波函数的MATLAB图像实现

基于量子力学谐振子波函数的MATLAB图像实现姓名：陈万学号：13020011006专业年级：物理学2013级指导老师：顾永建量子力学训练项目１基于量子力学谐振子波函数的MATLAB图像实现陈万摘要：量子力学中定态薛定谔方程可精确求解的典型例子是线性谐振子。为了直观的理解这一模型，通过MATLAB绘制出了一维谐振子不同能级下波函数和概率密度分布。通过编程和作图加深对量子力学波函数和概率密度的认识。关键词：线性谐振子MATLAB1.引言简谐运动广泛存在于自然界中。任何体系在平衡位置附近的小振动(例如，分子的振动、晶格的振动、原子核表面振动以及辐射场的振动等)一般都可以看成是简谐运动即谐振动。通过对比经典谐振子可以加深对量子力学的认识。对众多实际的谐振动的求解一般是从简化的物理模型即线性谐振子出发，然后再考虑具体的物理情景。本文主要借助MATLAB语言，给出了一维线性谐振子的可视化图形，从理论公式和可视化图形两个角度理解谐振子的波函数和概率密度分布的特性。2.一维谐振子的理论研究对于一维谐振子势能函数为2212Vxmx(1)定态薛定谔方程为22222122dmxEmdx(2)引入无量纲变量mx(3)用，薛定谔方程可写为222()dKd(4)式中K是以(1/2)为单位的能量：EK2(5)分析知ψ(x)可以写成如下形式２22()()xhe(6)薛定谔方程变为222(1)0dhdhKhdd(7)引入厄米多项式22()(1)()nnndHeed(8)求得归一化的谐振子定态是2412)(!21)()(eHnmxnn(9)物理上要求K=2n+1,也就是能量必须是,2,1,0,)21(nnnE(10)能量是量子化的。其中前几个厄米多项式为，10H,21H，2422H12833H3.MATLAB源程序%%不同能级下波函数和概率密度分布图像clearallclcm=1;%粒子质量k=1;%劲度系数３w=sqrt(k/m);%圆频率h_=1;%普朗克常量，这里为了作图方便，被放大了A=(m*w/(pi*h_))^(1/4);%归一化系数AsymsxdimV=sqrt(m*w/h_)*x;%thedimensionlessvariableξ%**************thefirstfewHermitepolynormials,Hn(ξ)H0=1;H1=2*dimV;H2=4*dimV^2-2;H3=8*dimV^3-12*dimV;%******************************************************V=0.5*k*x^2;%势能函数psi0=A*H0*exp(-(dimV^2)/2);psi1=A*(1/sqrt(2))*H1*exp(-(dimV^2)/2);psi2=A*(1/sqrt(8))*H2*exp(-dimV^2/2);psi3=A*(1/sqrt(48))*H3*exp(-dimV^2/2);%*********************ezplot(psi0);ylabel('\psi_0');title('\psi_0的图像');figureezplot(psi1);ylabel('\psi_1');title('\psi_1的图像');figureezplot(psi2);ylabel('\psi_2');title('\psi_2的图像');figureezplot(psi3);ylabel('\psi_3');title('\psi_3的图像');figure%单独的ψ(x)图像%*********************subplot(1,2,1);ezplot(V);holdonezplot(psi0);４ezplot(psi1+2);ezplot(psi2+4);ezplot(psi3+6);axis([-5507]);ylabel('\psi');title('势能和\psi的图像');%在一幅图中显示ψnsubplot(1,2,2);ezplot(V);holdonezplot(psi0^2);ezplot(psi1^2+2);ezplot(psi2^2+4);ezplot(psi3^2+6);axis([-5507]);ylabel('|\psi_n|^2');title('势能和|\psi_n|^2的图像');%在一幅图中显示|ψn|^2,即概率4.不同能级下波函数和概率密度分布图像５６5.总结本文系统地分析了一维谐振子的波函数和概率密度分布，并借助MATLAB语言画出了每种情况下的前4个能级的波函数和几率密度图像。从这些图像中，我们可以更加形象地理解比较抽象的谐振子问题，而且这些图像的绘制也锻炼了学７生的计算机应用能力，激发了学生的学习兴趣。此外在做图中，考虑到普朗克常量太小，进行了放大处理=1。还有m、等的参量也应取合适的值，才能使图像趋于完美。参考文献：[1]DavidJ.Griffiths.IntroductiontoQuantumMechanics[M].北京：机械工业出版社，2006：(51-56)[2]张迪,王丽,姜其畅,苏艳丽.量子力学中线性谐振子的可视化研究[J].运城学院学报，2015(3):(34-36)