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

量子力学训练项目基于量子力学谐振子波函数的MATLAB图像实现姓名：陈万 学号专业年级：物理学2013级 指导老师：顾永建基于量子力学谐振子波函数的MATLAB图像实现陈万摘要：量子力学中定态薛定谔方程可精确求解的典型例子是线性谐振子。为了直观的理解这一模型，通过MATLAB绘制出了一维谐振子不同能级下波函数和概率密度分布。通过编程和作图加深对量子力学波函数和概率密度的认识。关键词：线性谐振子 MATLAB 引言简谐运动广泛存在于自然界中。任何体系在平衡位置附近的小振动( 例如，分子的振动、晶格的振动、原子核表面振动以及辐射场的振动等) 一般都可以看成是简谐运动即谐振动。通过对比经典谐振子可以加深对量子力学的认识。对众多实际的谐振动的求解一般是从简化的物理模型即线性谐振子出发，然后再考虑具体的物理情景。本文主要借助MATLAB 语言，给出了一维线性谐振子的可视化图形，从理论公式和可视化图形两个角度理解谐振子的波函数和概率密度分布的特性。 一维谐振子的理论研究对于一维谐振子势能函数为 (1)定态薛定谔方程为 (2)引入无量纲变量 (3)用，薛定谔方程可写为 (4)式中K是以(1/2)为单位的能量： (5)分析知ψ(x)可以写成如下形式 (6)薛定谔方程变为 (7)引入厄米多项式 (8)求得归一化的谐振子定态是 (9)物理上要求K=2n+1,也就是能量必须是 (10)能量是量子化的。其中前几个厄米多项式为 MATLAB源程序%%不同能级下波函数和概率密度分布图像clear allclcm=1; %粒子质量k=1; %劲度系数w=sqrt(k/m); %圆频率h_=1; %普朗克常量，这里为了作图方便，被放大了A=(m*w/(pi*h_))^(1/4); %归一化系数Asyms xdimV=sqrt(m*w/h_)*x; %the dimensionless variableξ%**************the first few Hermite polynormials,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);hold onezplot(psi0);ezplot(psi1+2);ezplot(psi2+4);ezplot(psi3+6);axis([-5 5 0 7]);ylabel(\psi);title(势能和\psi的图像);%在一幅图中显示ψns