| 查看: 1176 | 回复: 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)的曲线图 |
回复此楼» 猜你喜欢
这样的filecode谁见过
已经有13人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
霜小妹二
木虫 (正式写手)
哈哈
- 应助: 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
3楼2016-05-05 10:08:58
霜小妹二
木虫 (正式写手)
哈哈
- 应助: 154 (高中生)
- 金币: 1806.7
- 红花: 25
- 帖子: 531
- 在线: 66小时
- 虫号: 3237838
- 注册: 2014-05-27
- 专业: 机械动力学
4楼2016-05-05 11:10:20
20