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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】求高人帮我修改程序
51/1返回列表
查看: 947  |  回复: 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

铁杆木虫 (著名写手)

方丈大师
F=0:0.001:10     %这个步长精度才能画出放大图的效果
赞一下 回复此楼
4楼2010-06-12 00:30:38
已阅   回复此楼   关注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

铁杆木虫 (著名写手)

方丈大师
★
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的回帖

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 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 318求调剂 +4 plum李子 2026-03-21 7/350 2026-03-22 14:17 by ColorlessPI
[基金申请] 山东省面上项目限额评审 +4 石瑞0426 2026-03-19 4/200 2026-03-22 08:50 by Wei_ren
[考研] 南昌大学材料专硕311分求调剂 +6 77chaselx 2026-03-20 6/300 2026-03-21 07:24 by JourneyLucky
[考研] 机械专硕299求调剂至材料 +3 kkcoco25 2026-03-16 4/200 2026-03-21 03:52 by JourneyLucky
[考研] 材料工程（专）一志愿985 初试335求调剂 +3 hiloiy 2026-03-17 4/200 2026-03-21 03:04 by JourneyLucky
[考研] 初始318分求调剂（有工作经验） +3 1911236844 2026-03-17 3/150 2026-03-21 02:33 by JourneyLucky
[考研] 华东师范大学-071000生物学-293分-求调剂 +3 研究生何瑶明 2026-03-18 3/150 2026-03-21 01:30 by JourneyLucky
[考研] 一志愿西南交大，求调剂 +5 材化逐梦人 2026-03-18 5/250 2026-03-21 00:26 by JourneyLucky
[考研] 308求调剂 +3 阿姐阿姐家啊 2026-03-18 3/150 2026-03-20 23:24 by JourneyLucky
[考研] 一志愿中海洋材料工程专硕330分求调剂 +8 小材化本科 2026-03-18 8/400 2026-03-20 23:16 by JourneyLucky
[考研] 一志愿武汉理工材料工程专硕调剂 +9 Doleres 2026-03-19 9/450 2026-03-20 22:36 by JourneyLucky
[考研] 287求调剂 +7 晨昏线与星海 2026-03-19 8/400 2026-03-20 22:19 by JourneyLucky
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
[考研] 281求调剂（0805） +14 烟汐忆海 2026-03-16 25/1250 2026-03-20 15:47 by yuncha
[考研] 298-一志愿中国农业大学-求调剂 +9 手机用户 2026-03-17 9/450 2026-03-20 14:24 by 无懈可击111
[考研] 288求调剂，一志愿华南理工大学071005 +5 ioodiiij 2026-03-17 5/250 2026-03-19 18:22 by zcl123
[考研] 一志愿福大288有机化学，求调剂 +3 小木虫200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
[考研] 0854可跨调剂，一作一项核心论文五项专利，省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 材料，纺织，生物（0856、0710），化学招生啦 +3 Eember. 2026-03-17 9/450 2026-03-18 10:28 by Eember.
信息提示
关闭
请填处理意见
关闭