 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 959 | 回復(fù): 5 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
fanbao_keke鐵蟲(chóng) (正式寫(xiě)手)
|
[交流]
【求助】求高人幫我修改程序 已有3人參與
|
||
|
我編的下面這個(gè)程序畫(huà)出的圖形與原圖不一致，可能是在能量的排序上存在問(wèn)題，我不知如何修改，請(qǐng)大家?guī)臀艺艺医鉀Q的辦法，萬(wàn)分感謝！ 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 ] |
回復(fù)此樓» 猜你喜歡
352分-085602-一志愿985
已經(jīng)有3人回復(fù)
-
291求調(diào)劑
已經(jīng)有4人回復(fù)
-
330分求調(diào)劑
已經(jīng)有4人回復(fù)
-
309求調(diào)劑
已經(jīng)有4人回復(fù)
-
0856求調(diào)劑
已經(jīng)有13人回復(fù)
-
0703化學(xué)
已經(jīng)有11人回復(fù)
-
329求調(diào)劑
已經(jīng)有7人回復(fù)
-
321求調(diào)劑
已經(jīng)有8人回復(fù)
-
0856材料化工調(diào)劑 總分330
已經(jīng)有12人回復(fù)
-
070300化學(xué)354求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 催化版上2008-04-25前 沸石分子篩，ZSM-系、SAPO-系 共466個(gè)話題集錦
已經(jīng)有110人回復(fù)
change0618
鐵桿木蟲(chóng) (著名寫(xiě)手)
方丈大師
- 應(yīng)助: 44 (小學(xué)生)
- 金幣: 17724.5
- 紅花: 17
- 帖子: 2413
- 在線: 546.7小時(shí)
- 蟲(chóng)號(hào): 496517
- 注冊(cè): 2008-01-19
- 專(zhuān)業(yè): 化學(xué)反應(yīng)工程
★
robert2020(金幣+1):多謝應(yīng)助，歡迎常來(lái)仿真模擬版交流。 2010-06-12 18:43:23
fanbao_keke(金幣+10):謝謝！只是現(xiàn)在程序中本征值在排序上還是有點(diǎn)問(wèn)題，沒(méi)法與原圖一致，能否幫我一下。 2010-06-12 21:57:26
fanbao_keke(金幣+20): 2010-07-22 07:44:02
robert2020(金幣+1):多謝應(yīng)助，歡迎常來(lái)仿真模擬版交流。 2010-06-12 18:43:23
fanbao_keke(金幣+10):謝謝！只是現(xiàn)在程序中本征值在排序上還是有點(diǎn)問(wèn)題，沒(méi)法與原圖一致，能否幫我一下。 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); 計(jì)算精度的原因出現(xiàn)了復(fù)數(shù) |
3樓2010-06-12 00:25:35
fanbao_keke
鐵蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 60.6
- 散金: 35
- 帖子: 388
- 在線: 52.3小時(shí)
- 蟲(chóng)號(hào): 987790
- 注冊(cè): 2010-04-02
- 性別: GG
- 專(zhuān)業(yè): 動(dòng)力學(xué)與控制
2樓2010-06-11 23:46:26
change0618
鐵桿木蟲(chóng) (著名寫(xiě)手)
方丈大師
- 應(yīng)助: 44 (小學(xué)生)
- 金幣: 17724.5
- 紅花: 17
- 帖子: 2413
- 在線: 546.7小時(shí)
- 蟲(chóng)號(hào): 496517
- 注冊(cè): 2008-01-19
- 專(zhuān)業(yè): 化學(xué)反應(yīng)工程
4樓2010-06-12 00:30:38
fanbao_keke
鐵蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 60.6
- 散金: 35
- 帖子: 388
- 在線: 52.3小時(shí)
- 蟲(chóng)號(hào): 987790
- 注冊(cè): 2010-04-02
- 性別: GG
- 專(zhuān)業(yè): 動(dòng)力學(xué)與控制
這兩個(gè)圖還是不一致，在原圖中第一個(gè)區(qū)域最大的那條線在第二個(gè)區(qū)域成為第二大（如粉紅線），然后交替變化，同樣第一個(gè)區(qū)域最小的那條線在第二個(gè)區(qū)域?yàn)榈谌�，而第三條線變?yōu)樽钚〉模缓蠼惶孀兓�，而我這個(gè)程序畫(huà)出的圖中最大的始終最大，最小的始終最小，沒(méi)有交叉，這應(yīng)該是在本征值的排序上處理錯(cuò)了，不知如何改？麻煩大家?guī)蛶臀�，謝謝�。。�  [ Last edited by fanbao_keke on 2010-6-13 at 10:21 ] |
5樓2010-06-12 08:53:38
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 309求調(diào)劑 +4 | 誰(shuí)不是少年 2026-03-29 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 0856材料化工調(diào)劑 總分330 +12 | zhubinhao 2026-03-27 | 12/600 |
|
|
[考研] 299求調(diào)劑 +8 | 15188958825 2026-03-25 | 8/400 |
|
|
|
[考研] 332求調(diào)劑 +6 | 蕉蕉123 2026-03-28 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +3 | QiMing7 2026-03-25 | 4/200 |
|
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
|
[考研] 070300求調(diào)劑306分 +4 | 26要上岸 2026-03-27 | 4/200 |
|
|
|
[考研] 328求調(diào)劑 +7 | 嗯滴的基本都 2026-03-27 | 7/350 |
|
|
|
[考研] 315分求調(diào)劑 +7 | 26考研上岸版26 2026-03-26 | 7/350 |
|
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
|
[考研] 287求調(diào)劑 +10 | land xuxu 2026-03-26 | 10/500 |
|
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|
|
|
[考研] 一志愿211，335分，0856，求調(diào)劑院校和導(dǎo)師 +4 | 傾____蕭 2026-03-27 | 5/250 |
|
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
|
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4 | 鋼鐵大炮 2026-03-24 | 4/200 |
|
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
|
[考研] 26考研-291分-廈門(mén)大學(xué)（085601）-柔性電子學(xué)院材料工程專(zhuān)業(yè)求調(diào)劑 +3 | min3 2026-03-24 | 4/200 |
|
|
|
[考研] 材料/農(nóng)業(yè)專(zhuān)業(yè)，07/08開(kāi)頭均可，過(guò)線就行 +3 | 呵唔哦豁 2026-03-23 | 4/200 |
|
|
|
[考研] 336化工調(diào)劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|
|
|
[考研] 280分求調(diào)劑 一志愿085802 +4 | PUMPT 2026-03-22 | 7/350 |
|