24小时热门版块排行榜

>论坛更新日志 (2803)
>虫友互识 (400)
>文献求助 (121)
>导师招生 (84)
>休闲灌水 (53)
>论文投稿 (50)
>考研 (33)
>考博 (32)
>博后之家 (31)
>基金申请 (28)
>硕博家园 (26)
>公派出国 (21)
>论文道贺祈福 (18)
>教师之家 (15)
>招聘信息布告栏 (14)
>药学 (13)

返回列表

田山东

捐助贵宾 (著名写手)

应助: 5 (幼儿园)
金币: 77
散金: 739
红花: 8
帖子: 1577
在线: 138.1小时
虫号: 800307
注册: 2009-06-28
性别: GG
专业: 基础物理学

[求助] H*pinv(H)为什么在非对角元出现非零值？

如题，代码如下。多谢各位同学啊！！
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;

%设定zigzag边界系统哈密顿矩阵
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';

%系统整体哈密顿矩阵
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;

这里应该得到单位矩阵，但是在HH非对角元得到非零值。很奇怪啊！！求高人解答！！

[ Last edited by 田山东 on 2013-8-22 at 16:05 ]

回复此楼