24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

fanbao_keke

铁虫 (正式写手)

应助: 0 (幼儿园)
金币: 60.6
散金: 35
帖子: 388
在线: 52.3小时
虫号: 987790
注册: 2010-04-02
性别: GG
专业: 动力学与控制

[交流] 【求助】求高人帮我修改程序已有3人参与

我编的下面这个程序画出的图形与原图不一致，可能是在能量的排序上存在问题，我不知如何修改，请大家帮我找找解决的办法，万分感谢！
clear all
omega=1;
m=1;
for F=0:0.1:10
J=40;  T=2*pi/omega;
tstep=T/J;  U=eye(4,4);
for t=tstep:tstep:T
H_0=zeros(4,4);
H_0(1,2)=0.2;
H_0(2,1)=0.2;
H_0(2,3)=0.2;
H_0(3,2)=0.2;
H_0(3,4)=0.2;
H_0(4,3)=0.2;
H=zeros(4,4);
for n=1:1:4
H(n,n)=n*F*cos(omega*t);
end
U=expm(-i*(H_0+H)*tstep)*U;
end

[C,Ueig]=eig(U);
QuasiEnergy=i*log(diag(Ueig))/T;
ExpectH_0=C'*H_0*C;
[ExpSort,Index]=sort(diag(ExpectH_0));
E(:,m)=(QuasiEnergy(Index(1:4)));
m=m+1;
end
plot(0:0.1:10,sort(E))

以下是原图形：

[ Last edited by fanbao_keke on 2010-6-11 at 23:50 ]

回复此楼

» 猜你喜欢

论文终于录用啦！满足毕业条件了已经有10人回复
2025年遐想已经有4人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复
自然科学基金委宣布启动申请书“瘦身提质”行动已经有4人回复
求个博导看看已经有18人回复

» 本主题相关价值贴推荐，对您同样有帮助:

催化版上2008-04-25前沸石分子筛，ZSM-系、SAPO-系共466个话题集锦已经有110人回复

1楼 2010-06-11 23:35:30

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师

应助: 44 (小学生)
金币: 17724.5
红花: 17
帖子: 2413
在线: 546.7小时
虫号: 496517
注册: 2008-01-19
专业: 化学反应工程

★
robert2020(金币+1):多谢应助，欢迎常来仿真模拟版交流。 2010-06-12 18:43:23
fanbao_keke(金币+10):谢谢！只是现在程序中本征值在排序上还是有点问题，没法与原图一致，能否帮我一下。 2010-06-12 21:57:26
fanbao_keke(金币+20): 2010-07-22 07:44:02

QuasiEnergy=real(i*log(diag(Ueig))/T);
ExpectH_0=real(C'*H_0*C);

计算精度的原因出现了复数