 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 501 | 回復(fù): 1 | |||
[求助]
H*pinv(H)為什么在非對角元出現(xiàn)非零值？
|
|
如題，代碼如下。多謝各位同學(xué)�。�！ clc; clear; N=8; M=5; A1=zeros(N,N); A2=zeros(N,N); B1=zeros(N,N); B2=zeros(N,N); C1=zeros(N,N); C2=zeros(N,N); H=zeros(M*N,M*N); Gr=zeros(M*N,M*N); t=2.8; %設(shè)定zigzag邊界系統(tǒng)哈密頓矩陣 for j=1:N-1 if(mod(j,2)==0) A1(j,j+1)=-t; A1(j+1,j)=-t; end if(mod(j,2)==1) A2(j,j+1)=-t; A2(j+1,j)=-t; end end for j=1:N B1(j,j)=-t; B2(j,j)=-t; end C2=B1'; C1=B2'; %系統(tǒng)整體哈密頓矩陣 for j=1:M if(mod(j,2)==1) H((j-1)*N+1:j*N, (j-1)*N+1:j*N)=A1; else H((j-1)*N+1:j*N, (j-1)*N+1:j*N)=A2; end end for j=1:M if(mod(j,2)==1) if j==1 H((j-1)*N+1:j*N,j*N+1: (j+1)*N)=B1; else if j<M H((j-1)*N+1:j*N,j*N+1: (j+1)*N)=B1; H((j-1)*N+1:j*N, (j-2)*N+1: (j-1)*N)=C1; else H((j-1)*N+1:j*N, (j-2)*N+1: (j-1)*N)=C1; end end else H((j-1)*N+1:j*N,j*N+1: (j+1)*N)=B2; H((j-1)*N+1:j*N, (j-2)*N+1: (j-1)*N)=C2; end end Gr=pinv(H); HH=Gr*H; 這里應(yīng)該得到單位矩陣，但是在HH非對角元得到非零值。很奇怪 �。�！求高人解答！！ [ Last edited by 田山東 on 2013-8-22 at 16:05 ] |
回復(fù)此樓» 猜你喜歡
347求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿西交288化工專碩求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿中國科學(xué)院大學(xué)265求調(diào)劑
已經(jīng)有4人回復(fù)
-
291求調(diào)劑
已經(jīng)有13人回復(fù)
-
070300化學(xué)求調(diào)劑
已經(jīng)有3人回復(fù)
-
322求調(diào)劑
已經(jīng)有7人回復(fù)
-
0703化學(xué)
已經(jīng)有15人回復(fù)
-
材料學(xué)碩333求調(diào)劑
已經(jīng)有15人回復(fù)
-
313求調(diào)劑
已經(jīng)有3人回復(fù)
-
083000學(xué)碩274求調(diào)劑
已經(jīng)有11人回復(fù)
alorsmail
鐵桿木蟲 (著名寫手)
- 應(yīng)助: 52 (初中生)
- 金幣: 7398.3
- 紅花: 41
- 帖子: 2104
- 在線: 237小時
- 蟲號: 2537936
- 注冊: 2013-07-09
- 專業(yè): 系統(tǒng)生物學(xué)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
xzhdty: 金幣+1, 感謝參與 2013-08-22 20:02:01
田山東: 金幣+5, ★★★很有幫助, 多謝回復(fù)，一會摸索下 2013-08-22 20:10:03
感謝參與，應(yīng)助指數(shù) +1
xzhdty: 金幣+1, 感謝參與 2013-08-22 20:02:01
田山東: 金幣+5, ★★★很有幫助, 多謝回復(fù)，一會摸索下 2013-08-22 20:10:03
|
這里應(yīng)該是誤差， 如果你對原 HH 做以下變換的話，就會得到單位矩陣 threshold = 0.5; HH1 = double(HH>=threshold); 當(dāng)然，如果你取 threshold =0.1; 的話，你會得到另外兩個非零 對角 研究一下為什么 threshold 會影響你的結(jié)果， |
2樓2013-08-22 17:04:07
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 田山東 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 070300化學(xué)求調(diào)劑 +3 | 小黃鴨寶 2026-03-30 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 材料科學(xué)與工程求調(diào)劑 +10 | 深V宿舍吧 2026-03-29 | 10/500 |
|
|
[考研] 370求調(diào)劑 +3 | 080700調(diào)劑 2026-03-30 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | 圖鑒212 2026-03-30 | 3/150 |
|
|
|
[考研] 一志愿中海洋材料357 +3 | 麥恩莉. 2026-03-30 | 3/150 |
|
|
|
[考研] 332求調(diào)劑 +6 | @MZB382400 2026-03-28 | 6/300 |
|
|
|
[考研] 284求調(diào)劑 +14 | junqihahaha 2026-03-26 | 15/750 |
|
|
|
[考研] 348求調(diào)劑 +6 | 小懶蟲不懶了 2026-03-28 | 6/300 |
|
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +4 | 邱gl 2026-03-26 | 7/350 |
|
|
|
[考研] 求調(diào)劑 +4 | QiMing7 2026-03-25 | 5/250 |
|
|
|
[考研] 299求調(diào)劑 +10 | 15188958825 2026-03-25 | 10/500 |
|
|
|
[考研] 298求調(diào)劑 +3 | 種圣賜 2026-03-29 | 3/150 |
|
|
|
[考研] 332求92調(diào)劑 +8 | 蕉蕉123 2026-03-28 | 8/400 |
|
|
|
[考研] 339求調(diào)劑，想調(diào)回江蘇 +6 | 烤麥芽 2026-03-27 | 8/400 |
|
|
|
[考研] 086000調(diào)劑 +3 | 7901117076 2026-03-26 | 3/150 |
|
|
|
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +4 | 哇呼哼呼哼 2026-03-26 | 4/200 |
|
|
|
[考研] 調(diào)劑推薦 +5 | 清酒714 2026-03-26 | 6/300 |
|
|
|
[考研] 292求調(diào)劑 +4 | 求求了收下我吧?/a> 2026-03-26 | 4/200 |
|
|
|
[考研] 324求調(diào)劑 +5 | hanamiko 2026-03-26 | 5/250 |
|
|
|
[考研] 334分 一志愿武理-080500 材料求調(diào)劑 +4 | 李李不服輸 2026-03-25 | 4/200 |
|