| 查看: 876 | 回复: 5 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
fanbao_keke铁虫 (正式写手)
|
[交流]
【求助】求高人帮我修改程序 已有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人回复
» 本主题相关价值贴推荐，对您同样有帮助:
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
3楼2010-06-12 00:25:35
fanbao_keke
铁虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 60.6
- 散金: 35
- 帖子: 388
- 在线: 52.3小时
- 虫号: 987790
- 注册: 2010-04-02
- 性别: GG
- 专业: 动力学与控制
2楼2010-06-11 23:46:26
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
4楼2010-06-12 00:30:38
fanbao_keke
铁虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 60.6
- 散金: 35
- 帖子: 388
- 在线: 52.3小时
- 虫号: 987790
- 注册: 2010-04-02
- 性别: GG
- 专业: 动力学与控制
5楼2010-06-12 08:53:38