generatega.m

来自「这是国标数字电视地面传输中生成矩阵的生成程序」· M 代码 · 共 65 行

%before  start this program,we should fill the GA cell with 24*35 hexdigit
%line in wordfile GflineA(修改).doc
GAfinal=zeros(3048,4445);
for m1=1:24
    for n1=1:35
        %transform hex to bin
           hex=GA{m1,n1};
           bin=zeros(1,128);%由于是十六进制转换为二进制，应该是128bit,而在标准          
           for k=1:32       %中是把最低位的“0”丢弃了，剩下127bit
                 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;
         GAfinal(i2:x1,j2:y1)=subG;
    end
end