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 ]

回复此楼

» 猜你喜欢

上海工程技术大学【激光智能制造】课题组招收硕士已经有6人回复
带资进组求博导收留已经有11人回复
自荐读博已经有5人回复
求个博导看看已经有16人回复
上海工程技术大学张培磊教授团队招收博士生已经有4人回复
求助院士们，这个如何合成呀已经有4人回复
临港实验室与上科大联培博士招生1名已经有9人回复
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？已经有6人回复
最近几年招的学生写论文不引自己组发的文章已经有11人回复
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）已经有3人回复

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

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

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

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

fanbao_keke

铁虫 (正式写手)

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

大家帮帮忙，谢谢！！！

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

回复此楼

2楼2010-06-11 23:46:26

已阅回复此楼关注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);

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

赞一下(1人)

回复此楼

3楼2010-06-12 00:25:35

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

change0618

铁杆木虫 (著名写手)

方丈大师

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

F=0:0.001:10 %这个步长精度才能画出放大图的效果

赞一下

回复此楼

4楼2010-06-12 00:30:38

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

fanbao_keke

铁虫 (正式写手)

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

这两个图还是不一致，在原图中第一个区域最大的那条线在第二个区域成为第二大（如粉红线），然后交替变化，同样第一个区域最小的那条线在第二个区域为第三大，而第三条线变为最小的，然后交替变化，而我这个程序画出的图中最大的始终最大，最小的始终最小，没有交叉，这应该是在本征值的排序上处理错了，不知如何改？麻烦大家帮帮我，谢谢！！！