一维有限深方势阱

MATLAB源代码，一维有限深方势阱，量子力学专业课程，求解能级与波函数

clc;clear;
a=50*1e-10;%m
me=9.1*1e-31;%kg
V0=1;%eV
V_J=1.602e-19;
h_bar=1.0545*1e-34;%J*s
k0=sqrt（2*me*V0*V_J/（h_bar^2））;
z0=k0*a/2;
E=-1:0.001:0;
k=sqrt（2*me.*（E+V0）*V_J/h_bar^2）;
z=k*a/2;
y1=tan（z）;
y2=sqrt（（z0./z）.^2-1）;
figure（1）
plot（Ey1‘r--‘Ey2‘g.-‘）
legend（‘tan（z）‘‘sqrt（（z0./z）.^2-1）‘）
grid on
xlabel（‘E（eV）‘）
axis（[-1 0 -18 18]）
y=y1-y2;
[n1m1]=find（abs（y（1:100））<1）;
E1=E（m1）;
[n2m2]=find（abs（y（100:200））<0.1）;
E2=E（100:200）;
E2=E2（m2）;
[n2m3]=find（abs（y（300:400））<0.01）;
E3=E（300:400）;
E3=E3（m3）;
EE=[E1E2E3];
figure（2）
plot（[-50 -25 -25 25 25 50][0 0 -1 -1 0 0]）;hold on
for i=1:3
plot（[-25 25][EE（i） EE（i）]）;hold on
end
axis（[-50 50 -1.2 0.2]）
text（-10-0.05‘阱宽50A‘）;%文本标注
text（-15-1.1‘阱深1eV‘）;%文本标注
text（-15E1‘基态-0.9870eV‘）;%文本标

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1656  2015-04-08 14:04  Untitled.m