encode1.m

来自「除法电路编码」· M 代码 · 共 74 行

M=zeros(64,64);
rand('state',sum(100*clock));
for i=1:57;
    for j=1:57;
        M(i,j)=floor(rand(1)+0.5);
    end;
end;
G=[   1     0     0     0     0     1     1;
      1     1     0     0     0     1     1;
      1     1     1     0     0     1     1;
      1     1     1     1     0     1     1;
      1     1     1     1     1     1     1;
      1     1     1     1     1     0     1;
      0     1     1     1     1     1     1;
      1     0     1     1     1     0     1;
      0     1     0     1     1     1     1;
      1     0     1     0     1     0     1;
      0     1     0     1     0     1     1;
      1     0     1     0     1     1     1;
      1     1     0     1     0     0     1;
      0     1     1     0     1     0     1;
      0     0     1     1     0     1     1;
      1     0     0     1     1     1     1;
      1     1     0     0     1     0     1;
      0     1     1     0     0     1     1;
      1     0     1     1     0     1     1;
      1     1     0     1     1     1     1;
      1     1     1     0     1     0     1;
      0     1     1     1     0     1     1;
      1     0     1     1     1     1     1;
      1     1     0     1     1     0     1;
      0     1     1     0     1     1     1;
      1     0     1     1     0     0     1;
      0     1     0     1     1     0     1;
      0     0     1     0     1     1     1;
      1     0     0     1     0     0     1;
      0     1     0     0     1     0     1;
      0     0     1     0     0     1     1;
      1     0     0     1     0     1     1;
      1     1     0     0     1     1     1;
      1     1     1     0     0     0     1;
      0     1     1     1     0     0     1;
      0     0     1     1     1     0     1;
      0     0     0     1     1     1     1;
      1     0     0     0     1     0     1;
      0     1     0     0     0     1     1;
      1     0     1     0     0     1     1;
      1     1     0     1     0     1     1;
      1     1     1     0     1     1     1;
      1     1     1     1     0     0     1;
      0     1     1     1     1     0     1;
      0     0     1     1     1     1     1;
      1     0     0     1     1     0     1;
      0     1     0     0     1     1     1;
      1     0     1     0     0     0     1;
      0     1     0     1     0     0     1;
      0     0     1     0     1     0     1;
      0     0     0     1     0     1     1;
      1     0     0     0     1     1     1;
      1     1     0     0     0     0     1;
      0     1     1     0     0     0     1;
      0     0     1     1     0     0     1;
      0     0     0     1     1     0     1;
      0     0     0     0     1     1     1; 
    ];
c1=mod(M(1:57,1:57)*G,2);
M1=[M(1:57,1:57) c1];
M2=M1';
d=mod(M2*G,2);
M3=[M2,d];
M4=M3';
G'