应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】求高人帮我修改程序
51/1返回列表
查看: 951  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-06-11 23:35:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师
★
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的回帖
查看全部 6 个回答

fanbao_keke

铁虫 (正式写手)
大家帮帮忙，谢谢！！！

[ Last edited by fanbao_keke on 2010-6-11 at 23:48 ]
回复此楼
2楼2010-06-11 23:46:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师
F=0:0.001:10     %这个步长精度才能画出放大图的效果
赞一下 回复此楼
4楼2010-06-12 00:30:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanbao_keke

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



[ Last edited by fanbao_keke on 2010-6-13 at 10:21 ]
赞一下 回复此楼
5楼2010-06-12 08:53:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 284求调剂 +3 Zhao anqi 2026-03-22 3/150 2026-03-22 15:27 by ms629
[考研] 一志愿中南大学化学学硕0703总分337求调剂 +3 niko- 2026-03-22 3/150 2026-03-22 15:15 by 杨杨杨紫
[考研] 材料与化工085600，总分304，本科有两篇sci参与，求调剂 +3 幸运的酱酱 2026-03-22 4/200 2026-03-22 15:06 by allen-yin
[考研] 生物学071000 329分求调剂 +4 我爱生物生物爱� 2026-03-17 4/200 2026-03-22 08:34 by hxsm
[考研] 333求调剂 +5 87639 2026-03-21 7/350 2026-03-21 19:31 by ColorlessPI
[考研] 0805 316求调剂 +3 大雪深藏 2026-03-18 3/150 2026-03-21 18:55 by 学员8dgXkO
[考研] 一志愿深大，0703化学，总分302，求调剂 +4 七月-七七 2026-03-21 4/200 2026-03-21 18:20 by 学员8dgXkO
[考研] 311求调剂 +3 勇敢的小吴 2026-03-20 3/150 2026-03-21 17:40 by ColorlessPI
[考研] 一志愿天津大学化学工艺专业（081702）315分求调剂 +12 yangfz 2026-03-17 12/600 2026-03-21 03:30 by JourneyLucky
[考研] 一志愿中国石油大学（华东） 本科齐鲁工业大学 +3 石能伟 2026-03-17 3/150 2026-03-21 02:22 by JourneyLucky
[考研] 华东师范大学-071000生物学-293分-求调剂 +3 研究生何瑶明 2026-03-18 3/150 2026-03-21 01:30 by JourneyLucky
[考研] 一志愿 西北大学 ，070300化学学硕，总分287，双非一本，求调剂。 +3 晨昏线与星海 2026-03-18 3/150 2026-03-21 00:46 by JourneyLucky
[考研] 一志愿西南交大，求调剂 +5 材化逐梦人 2026-03-18 5/250 2026-03-21 00:26 by JourneyLucky
[考研] 324求调剂 +5 lucky呀呀呀鸭 2026-03-20 5/250 2026-03-20 22:30 by 促天成
[考研] 中南大学化学学硕337求调剂 +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[考研] 求调剂一志愿南京航空航天大学289分 +3 @taotao 2026-03-19 3/150 2026-03-20 21:34 by JourneyLucky
[考研] 求调剂 +3 暗涌afhb 2026-03-16 3/150 2026-03-20 00:28 by 河南大学校友
[考研] 288求调剂，一志愿华南理工大学071005 +5 ioodiiij 2026-03-17 5/250 2026-03-19 18:22 by zcl123
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
信息提示
关闭
请填处理意见
关闭