PART.5 模型应用:
一维无限深势阱的定态薛定谔方程能级的数值方法
仍是求解薛定谔方程,此时势函数由
积分法代码:
function schrodingerEq_el1
% 使用积分的形式获取矩阵元并求得特征值
% 循环计算矩阵元素
N=16;S=zeros(N,N);H=zeros(N,N);
for n=0:N
for m=0:N
S(n+1,m+1)=compute_Snm(n,m);
H(n+1,m+1)=compute_Hnm(n,m);
disp(['元素计算进度: ',num2str(n),' ',num2str(m)])
end
end
% 计算广义特征值
ERG=eig(H,S);
ERG=sort(ERG);