%****************************************************************
% Program 4:本程序完全从数值的角度出发,把求解薛定谔方程的过程看成是一个
% 求解哈密顿数值矩阵的特征值以及对应特征向量的一个过程
%****************************************************************
% 确定参数|求解薛定谔方程中使用到的一些参数:
hbar=1;m=1; %在原子单位体系里(hbar= 1,质量m= 1)
L = 5; % 区间长度为2L:-L < x < L
N = 1000;x = linspace(-L,L,N)'; % 总共在区间内均匀地取N个点
dx = x(2) - x(1); % 步长大小
% 选定U-x关系|在这里可以任意以多项式的形式进行组合 (这里首先假定的是一维谐振子的情况)
U = 1/2*100*x.^2; %一维谐振子满足的U-x关系
%U = 1/2*100*x.^2+1/2*x.^(4); %经过多项式组合的U-x关系(可以把上面U-x注释掉,替换成这个U-x关系)
%读者也可以自己以多项式的形式组合
%写出哈密顿的数值表达式
e = ones(N,1); Lap = spdiags([e -2*e e],[-1 0 1],N,N)/dx^2; %二阶导数的数值表达形式
H = -1/2*(hbar^2/m)*Lap + spdiags(U.*e,0,N,N); % 完整哈密顿矩阵
%计算能量本征值与能量本征函数
nmodes = 3; options.disp = 0;
[V,E] = eigs(H,nmodes,'sa',options); %利用求解特征值与特征向量的方式求出本征值与本征函数
[E,ind] = sort(diag(E));% 将E按数值大小从低到高排序
V = V(:,ind); %取出前三级的本证函数
Usc = U*max(abs(V(:)))/max(abs(U)); %为了便于观察,将U的范围按比例缩小U-x关系不变
%展示结果(如果读者认为展现的图片看不清其中的曲线,请调整显示的曲线数目)
plot(x,V,x,Usc,'--k'); %展现前三级在一定的U-x关系下本征函数V-x的关系
lgnd_str = [repmat('E=',nmodes,1),num2str(E)]; %添加标注
legend(lgnd_str)
shg %显示图形窗口
%%读者可以取消下面几行程序的注释,验证求得的能量本征函数是否满足归一化
% SumV = 0;
% for numV=1:1000
% SumV = SumV+V(numV,1)^2;
% end
% SumV
没有合适的资源?快使用搜索试试~ 我知道了~
温馨提示
【资源概览】 基于Matlab数值计算在解薛定谔方程当中的应用源码+项目资料齐全+部署说明文档 高分毕业设计.zip基于Matlab数值计算在解薛定谔方程当中的应用源码+项目资料齐全+部署说明文档 高分毕业设计.zip 【资源说明】 高分项目源码:此资源是在校高分项目的完整源代码,经过导师的悉心指导与认可,答辩评审得分高达95分,项目的质量与深度有保障。 测试运行成功:所有的项目代码在上传前都经过了严格的测试,确保在功能上完全符合预期,您可以放心下载并使用。 适用人群广泛:该项目不仅适合计算机相关专业(如软件工程、计科、区块链、人工智能、电子信息、物联网、通信工程、自动化等)的在校学生和老师,还可以作为毕业设计、课程设计、作业或项目初期立项的演示材料。对于希望进阶学习的小白来说,同样是一个极佳的学习资源。 代码灵活性高:如果您具备一定的编程基础,可以在此代码基础上进行个性化的修改,以实现更多功能。当然,直接用于毕业设计、课程设计或作业也是完全可行的。 欢迎下载,与我一起交流学习,共同进步!
资源推荐
资源详情
资源评论
收起资源包目录
基于Matlab数值计算在解薛定谔方程当中的应用源码+项目资料齐全+部署说明文档 高分毕业设计.zip (4个子文件)
部署说明文档.md 13KB
1D-harmonic-oscillator-master
harmonic_oscillator1D.m 2KB
result.png 46KB
171265889347208773632.zip 416B
共 4 条
- 1
资源评论
IT狂飙
- 粉丝: 4779
- 资源: 2658
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功