| 查看: 445 | 回复: 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 |
回复此楼» 猜你喜欢
一志愿安徽大学0817化学工程与技术，求调剂
已经有12人回复
-
求材料调剂，一志愿郑州大学289分
已经有3人回复
-
265求调剂
已经有9人回复
-
085600材料与化工调剂
已经有20人回复
-
专硕 351 086100 也是考的材科基 本科也是材料
已经有6人回复
-
085600专硕材料与化工348分求调剂
已经有10人回复
-
085600 295分求调剂
已经有21人回复
-
285求调剂
已经有5人回复
-
一志愿0817化学工程与技术，求调剂
已经有8人回复
-
271分求调剂学校
已经有3人回复
