crc_code_n1.m

CRC校验码的Matlab源程序,希望对大家有帮助

function C = CRC_code_n1(M,n,G)
%CRC code (n,k) function,每k个bit编成nbit。
%Input:
%   M:(1,x*k) input data
%   n:
%   G:generating polynomial.(1,n-k+1)
%Output:
%   C:(1,n) coded data;
len1=length(G);%len1=n-k+1 已经比n-k大1了。
k=n-len1+1;
len2=length(M);
M1=reshape(M,k,len2/k);
C1=zeros(n,len2/k);
C1(1:k,:)=M1;
%t=zeros(n-k+1,1);
%CRC 部分
for ii=1:len2/k
    t=C1(1:n-k+1,ii);%8=n-k;eg.n-k+1=9;  t为 n-k+1--1的向量
    for jj=n-k+2:n
        t=mod((t+t(1)*G),2);
        t=[t(2:n-k+1);C1(jj,ii)];
    end
        t=mod((t+t(1)*G),2);
        C1(k+1:n,ii)=t(2:n-k+1);%2-9,共8位
end
C=reshape(C1,1,n*len2/k);