| 查看: 1188 | 回复: 3 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
急急急 matlab程序运行不出来 哪位大神可以帮忙修改一下 谢谢 谢谢
|
||
|
程序是计算一维薛定谔方程的本征值，中间薛定谔矩阵用的是有限差分法，我们老师要求的程序运行结果，是要求出波函数phi的图，和能量E的值。 下面是我写的程序，程序目前的问题主要是矩阵赋值的地方，还有求能量E的地方如何写，本征值的地方。 各位大神帮帮忙，谢谢！ clc;clear all; %量子阱参数 w=30; %势阱宽30nm N=100; %取点 h=1.054*10^-34; %普朗克常数1.054*10^-34J-s m0=9.11*10^-31; %静电子质量9.11*10^-31kg m=0.067*m0; %有效质量 V=300; %势能300mev %薛定谔矩阵 a=0;b=30; %端点坐标 deltax=(b-a)/N; %计算区间步长 p=ones(N-1,1); %创建元素为1的（N-1）*1阶矩阵 K=spdiags([-2*p,-2*p,-2*p],[-1,0,1],N-1,N-1); %生成三对角线稀疏矩阵 for i=1:N-1 y1=a+deltax*i; %计算各节点坐标值 K(i,i)=h^2/m./(deltax).^2+V; %主对角线元素赋值 K(i,i+1)=-h^2/(2*m*(deltax).^2); K(i,i-1)= -h^2/(2*m*(deltax).^2); %从对角线元素赋值 end B=zeros(N-1,1); %创建元素为0的（N-1）*1阶矩阵 B(1)=h^2/(2*m*(deltax).^2); B(N-1)=h^2/(2*m*(deltax).^2); %设定矩阵两端的元素值 phi=K\B; %矩阵左除，计算各节点函数近似值 %本征值计算 E=eig(K); %求矩阵K的全部特征值 %显示结果 spy(K); %查看稀疏矩阵 plot(y1,phi,'.') %绘制函数phi(x)的曲线图 |
回复此楼» 猜你喜欢
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
已经有48人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有29人回复
-
关于Filecode分析方法
已经有12人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有11人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有10人回复
-
阿姨
已经有4人回复
-
一作与独作在应聘高校教师时区别大吗
已经有4人回复
-
静等基金结果
已经有21人回复
-
FileCode能看出啥？
已经有28人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有6人回复
3楼2016-05-05 10:08:58
霜小妹二
木虫 (正式写手)
哈哈
- 应助: 154 (高中生)
- 金币: 1806.7
- 红花: 25
- 帖子: 531
- 在线: 66小时
- 虫号: 3237838
- 注册: 2014-05-27
- 专业: 机械动力学
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
岚泽泽zz: 金币+5, ★有帮助 2016-05-05 10:10:34
感谢参与，应助指数 +1
岚泽泽zz: 金币+5, ★有帮助 2016-05-05 10:10:34
|
%%%%%%现在能够运行处结果，但是不知道对不对。 clc;clear all; %量子阱参数 w=30; %势阱宽30nm N=100; %取点 h=1.054*10^-34; %普朗克常数1.054*10^-34J-s m0=9.11*10^-31; %静电子质量9.11*10^-31kg m=0.067*m0; %有效质量 V=300; %势能300mev %薛定谔矩阵 a=0;b=30; %端点坐标 deltax=(b-a)/N; %计算区间步长 p=ones(N-1,1); %创建元素为1的（N-1）*1阶矩阵 y1=p; K=spdiags([-2*p,-2*p,-2*p],[-1,0,1],N-1,N-1); %生成三对角线稀疏矩阵 for i=1:N-1 y1(i)=a+deltax*i; %计算各节点坐标值 K(i,i)=h^2/m./(deltax).^2+V; %主对角线元素赋值 if i<N-1 K(i,i+1)=-h^2/(2*m*(deltax).^2); end if i>1 K(i,i-1)= -h^2/(2*m*(deltax).^2); %从对角线元素赋值 end end B=zeros(N-1,1); %创建元素为0的（N-1）*1阶矩阵 B(1)=h^2/(2*m*(deltax).^2); B(N-1)=h^2/(2*m*(deltax).^2); %设定矩阵两端的元素值 phi=K\B; %矩阵左除，计算各节点函数近似值 %本征值计算 E=eig(K); %求矩阵K的全部特征值 %显示结果 spy(K); %查看稀疏矩阵 plot(y1,phi,'b-') %绘制函数phi(x)的曲线图 |
2楼2016-05-05 09:33:33
霜小妹二
木虫 (正式写手)
哈哈
- 应助: 154 (高中生)
- 金币: 1806.7
- 红花: 25
- 帖子: 531
- 在线: 66小时
- 虫号: 3237838
- 注册: 2014-05-27
- 专业: 机械动力学
4楼2016-05-05 11:10:20
