| 查看: 451 | 回复: 0 | ||
fanbao_keke铁虫 (正式写手)
|
[求助]
求高手帮我修改一个求本征函数的程序
|
|
求高手帮我修改一个求本征函数的程序： clear all m=1 tt=1; ttt=0.1; w=10 for F=0:0.02:4 J=100; T=2*pi/w; tstep=T/J; U=eye(4,4); for t=tstep:tstep:T H=zeros(4,4); H_0=zeros(4,4); H_0(1,2)=tt; H_0(2,1)=tt; H_0(2,3)=tt; H_0(3,2)=tt; H_0(3,4)=tt; H_0(4,3)=tt; H_0(1,3)=ttt; H_0(3,1)=ttt; H_0(2,4)=ttt; H_0(4,2)=ttt; H(1,1)=F*w*sin(w*t); U=expm(-1i*(H_0+H)*tstep)*U; end [C,Ueig]=eig(U); QuasiEnergy=((1i*log(diag(Ueig))/T)); ExpectH_0=(C'*(H_0)*C); [ExpSort,Index]=sort(real(diag(ExpectH_0))); [CC,Ueig]=(eig(U)); % yy=CC(:,Index(2)) %这个命令下图2与图3的图混错 yy=CC(:,2) %这个命令下在横坐标2.4附近有几个点不连续，有不对的值。 yy1(m)=yy(1) yy2(m)=yy(2) yy3(m)=yy(3) yy4(m)=yy(4) m=m+1; end plot(0:0.02:4,abs(yy1).^2,'k.','MarkerSize',1.5) hold on plot(0:0.02:4,abs(yy2).^2,'k.','MarkerSize',1.5) hold on plot(0:0.02:4,abs(yy3).^2,'k.','MarkerSize',1.5) hold on plot(0:0.02:4,abs(yy4).^2,'k.','MarkerSize',1.5) 在计算本征函数时用 yy=CC(:,Index(2)) %这个命令下图2与图3的图混错 用yy=CC(:,2) %这个命令下在横坐标2.4附近有几个点不连续，有不对的值。 求助如何处理本征函数的排序问题？ 正确的图已给出在附件里。  新建 Microsoft Word Document (4).jpg |
回复此楼» 猜你喜欢
化学专业申博
已经有4人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
评审有感
已经有27人回复
-
论文撤稿了
已经有7人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
河北省自然科学基金
已经有9人回复
10