版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

田山東

捐助貴賓 (著名寫手)

應(yīng)助: 5 (幼兒園)
金幣: 77
散金: 739
紅花: 8
帖子: 1577
在線: 138.1小時
蟲號: 800307
注冊: 2009-06-28
性別: GG
專業(yè): 基礎(chǔ)物理學(xué)

[求助] 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ù)此樓

» 猜你喜歡

085601 329分調(diào)劑已經(jīng)有3人回復(fù)
本2一志愿C9-333分，材料科學(xué)與工程，求調(diào)劑已經(jīng)有3人回復(fù)
一志愿大連理工大學(xué)材料求調(diào)劑已經(jīng)有5人回復(fù)
材料與化工調(diào)劑一志愿大連海事085600，349 已經(jīng)有6人回復(fù)
一志愿中海洋320化學(xué)工程與技術(shù)學(xué)碩求調(diào)劑已經(jīng)有4人回復(fù)
材料科學(xué)與工程調(diào)劑已經(jīng)有11人回復(fù)
311求調(diào)劑一志愿合肥工業(yè)大學(xué) 已經(jīng)有4人回復(fù)
一志愿食品科學(xué)與工程083200求調(diào)劑已經(jīng)有4人回復(fù)
材料科學(xué)與工程求調(diào)劑已經(jīng)有10人回復(fù)
一志愿北化085600材料專碩275|有文章專利｜求調(diào)劑已經(jīng)有11人回復(fù)

everythinghasitsseason.enjoyyourlife.

1樓 2013-08-22 16:03:48

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

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)該是誤差，
如果你對原 HH 做以下變換的話，就會得到單位矩陣
threshold = 0.5;
HH1 = double(HH>=threshold);

當(dāng)然，如果你取 threshold =0.1; 的話，你會得到另外兩個非零對角
研究一下為什么 threshold 會影響你的結(jié)果，

贊一下(1人)

回復(fù)此樓

2樓2013-08-22 17:04:07