 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 954 | 回復(fù): 5 | |||
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ù)此樓» 猜你喜歡
328求調(diào)劑
已經(jīng)有7人回復(fù)
-
340求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿南昌大學(xué)324求調(diào)劑
已經(jīng)有4人回復(fù)
-
求調(diào)劑推薦 材料 304
已經(jīng)有15人回復(fù)
-
291求調(diào)劑
已經(jīng)有19人回復(fù)
-
331環(huán)境科學(xué)與工程求調(diào)劑
已經(jīng)有3人回復(fù)
-
289求調(diào)劑
已經(jīng)有8人回復(fù)
-
315分求調(diào)劑
已經(jīng)有7人回復(fù)
-
295求調(diào)劑
已經(jīng)有6人回復(fù)
-
291求調(diào)劑
已經(jīng)有7人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 催化版上2008-04-25前 沸石分子篩，ZSM-系、SAPO-系 共466個(gè)話題集錦
已經(jīng)有110人回復(fù)
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)工程
★
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
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ī)蛶臀遥x謝�。�！  [ Last edited by fanbao_keke on 2010-6-13 at 10:21 ] |
5樓2010-06-12 08:53:38
凡人大五
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1820.2
- 散金: 509
- 紅花: 4
- 帖子: 328
- 在線: 207.4小時(shí)
- 蟲(chóng)號(hào): 1265461
- 注冊(cè): 2011-04-14
- 專(zhuān)業(yè): 光學(xué)
6樓2017-07-18 14:35:19
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 fanbao_keke 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料求調(diào)劑 一志愿哈工大總分298分，前三科223分 +4 | dongfang59 2026-03-27 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 求調(diào)劑 +8 | 張zz111 2026-03-27 | 9/450 |
|
|
[考研] 317求調(diào)劑 +5 | 十閑wx 2026-03-24 | 5/250 |
|
|
|
[考研] 292求調(diào)劑 +4 | 求求了收下我吧?/a> 2026-03-26 | 4/200 |
|
|
|
[考研] 286求調(diào)劑 +4 | lim0922 2026-03-26 | 4/200 |
|
|
|
[考研] 276求調(diào)劑。有半年電池和半年高分子實(shí)習(xí)經(jīng)歷 +10 | 材料學(xué)257求調(diào)劑 2026-03-23 | 11/550 |
|
|
|
[考研] 317求調(diào)劑 +7 | 蛋黃咸肉粽 2026-03-26 | 7/350 |
|
|
|
[考研] 327求調(diào)劑 +7 | prayer13 2026-03-23 | 7/350 |
|
|
|
[考研] 352求調(diào)劑 +4 | 大米飯！ 2026-03-22 | 4/200 |
|
|
|
[考研] 機(jī)械學(xué)碩310分，數(shù)一英一，一志愿211本科雙非找調(diào)劑信息 +3 | @357 2026-03-25 | 3/150 |
|
|
|
[考研]
 085602 289分求調(diào)劑
+8
|
WWW西西弗斯 2026-03-24 | 8/400 |
|
|
|
[考研] 材料科學(xué)與工程 317求調(diào)劑 +4 | JKSOIID 2026-03-26 | 4/200 |
|
|
|
[考研] 材料277求調(diào)劑 +5 | min3 2026-03-24 | 5/250 |
|
|
|
[考研] 生物技術(shù)與工程 +3 | 1294608413 2026-03-25 | 4/200 |
|
|
|
[考研] 化學(xué)調(diào)劑 +6 | yzysaa 2026-03-21 | 6/300 |
|
|
|
[考研] 289材料與化工（085600）B區(qū)求調(diào)劑 +4 | 這么名字咋樣 2026-03-22 | 5/250 |
|
|
|
[考研] 一志愿武理085500機(jī)械專(zhuān)業(yè)總分300求調(diào)劑 +3 | an10101 2026-03-24 | 7/350 |
|
|
|
[考研] 284求調(diào)劑 +3 | yanzhixue111 2026-03-23 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +5 | Zhangbod 2026-03-21 | 7/350 |
|
|
|
[考研] 求助 +5 | 夢(mèng)里的無(wú)言 2026-03-21 | 6/300 |
|