generategb.m

数字电视标准矩阵及其生成矩阵源代码（LDPC)

 GBfinal=zeros(4572,7493);
for m1=1:36
    for n1=1:23
        %transform hex to bin
           hex=GB{m1,n1};
           bin=zeros(1,128);
           for k=1:32
                 i1=4.*(k-1)+1;
                 j1=4.*k;
                 switch hex(1,k)
                 case '0' 
                 bin(1,i1:j1)=[0 0 0 0];
                 case '1' 
                 bin(1,i1:j1)=[0 0 0 1];
                 case '2' 
                 bin(1,i1:j1)=[0 0 1 0];
                 case '3' 
                 bin(1,i1:j1)=[0 0 1 1];
                 case '4' 
                 bin(1,i1:j1)=[0 1 0 0];
                 case '5' 
                 bin(1,i1:j1)=[0 1 0 1];
                 case '6' 
                 bin(1,i1:j1)=[0 1 1 0];
                 case '7' 
                 bin(1,i1:j1)=[0 1 1 1];
                 case '8' 
                 bin(1,i1:j1)=[1 0 0 0];
                 case '9' 
                 bin(1,i1:j1)=[1 0 0 1];
                 case 'A' 
                 bin(1,i1:j1)=[1 0 1 0];
                 case 'B' 
                 bin(1,i1:j1)=[1 0 1 1];
                 case 'C' 
                 bin(1,i1:j1)=[1 1 0 0];
                 case 'D' 
                 bin(1,i1:j1)=[1 1 0 1];
                 case 'E' 
                 bin(1,i1:j1)=[1 1 1 0];
                 otherwise
                 bin(1,i1:j1)=[1 1 1 1];
                 end
             end
          subG=zeros(127,127);
          s=bin(2:128);
          subG(1,1:127)=s;
             for m2=2:127
                 s=circshift(s,[0,1]);
                 subG(m2,1:127)=s;
             end
          i2=127.*(m1-1)+1;
          j2=127.*(n1-1)+1;
          x1=127.*m1;
          y1=127.*n1;
         GBfinal(i2:x1,j2:y1)=subG;
    end
end
I=eye(127,127);
for m3=1:36
          n2=m3+23;
          i3=127.*(m3-1)+1;
          j3=127.*(n2-1)+1;
          x2=127.*m3;
          y2=127.*n2;
         GBfinal(i3:x2,j3:y2)=I;
end