clear,close all

load('E:\matlab2006\work\H.mat')
Hqc(1:1397,1:7493)=H
br=127
[cb,tb]=size(Hqc)
D=Hqc(1:cb,tb-cb+1:tb)
Hqc1=gauss(Hqc)
D1=gauss(D)
rd=rank(D1)
r= rank(Hqc1)
if r==cb
    delta_r=r-rd

switch delta_r
    case 0
        D=[D,eye(cb)]
        D=gauss(D)
        
        invD=D(1:cb,cb+1:2*cb)%求得D的逆陣
        M=Hqc(1:cb,1:tb-cb)
        TEM=mod(invD*M,2)
        I=eye(tb-cb)
        Z=mod(TEM*I,2)
        P=conj(Z)'
        Gqc=[I,P]
        
        
    otherwise %如果Hqc滿(mǎn)秩D非滿(mǎn)秩，則將Hqc序列重排
        tmp=Hqc(1:cb,tb-delta_r+1:tb)
        Hqc(1:cb,1+delta_r:tb)=Hqc(1:cb,1:tb-delta_r)
        Hqc(1:cb,1:delta_r)=tmp
        per_D=Hqc(1:cb,tb-cb+1:tb)
        per_D1=gauss(per_D)
        delta_l=0
        %如果D仍然非滿(mǎn)秩，則增加列數(shù)
        
       while (r-rank(per_D1))>0
            delta_l=delta_l+1
            per_D=Hqc(1:cb,tb-cb+1-delta_l*br:tb)
            per_D1=per_D,end
            lb=delta_l*br+cb, %得到新的列數(shù)
    end

else
    delta_r=r-rd
    switch delta_r
    case 0
        per_D=D1
        per_D1=gauss(per_D(1:rank(Hqc1),1:cb))
        Q=[per_D1;per_D(cb-rank(Hqc1)+1:cb,1:cb)]
        Q=Q(1:cb,cb-rank(Hqc1)+1:cb)
    otherwise
        tmp=Hqc(1:cb,tb-delta_r+1:tb)
        Hqc(1:cb,1+delta_r:tb)=Hqc(1:cb,1:tb-delta_r)
        Hqc(1:cb,1:delta_r)=tmp
        per_D=Hqc(1:cb,tb-cb+1:tb)
        per_D1=gauss(per_D)
        delta_l=0
       while (r-rank(per_D1))>0
            delta_l=delta_l+1
            per_D=Hqc(1:cb,tb-cb+1-delta_l*br:tb)
            per_D1=per_D,end
            lb=delta_l*br+cb, %得到新的列數(shù)
    end
        per_D2=gauss(per_D1(1:r,1:cb))
        w=[per_D2;per_D1(r+1:cb,1:cb)]
        q=w(1:cb,r+1:cb)


    
end