gauss.m

ldpc编码

[M N]=size(H);

HH=fliplr(H);      %按列左右交换
HH=flipud(H);      %按行上下交换

BakH=HH;

sign=0;
i=1;
while i<=M
    while ((i<=M) && (HH(i,i)==0))
        if sign==0
            HH=[HH;HH(i,:)];
            HH(i,:)=[];
            BakH=[BakH;BakH(i,:)];
            BakH(i,:)=[];
            sign=1;
        else
            HH=[HH,HH(:,i)];
            HH(:,i)=[];
            BakH=[BakH,BakH(:,i)];
            BakH(:,i)=[];
            sign=0;
        end
        if HH(i,:)==0
            HH(i,:)=[];
            M=M-1;
        end
    end
    
    if i<M
        L=eye(M);
        L(i+1:end,i)=-HH(i+1:end,i)/HH(i,i);
        HH=L*HH;
    end
    i=i+1;
    HH=mod(HH,2);
end
HH=fliplr(HH);
HH=flipud(HH);

BakH=fliplr(BakH);
BakH=flipud(BakH);

GauseH=HH;