function C = CRC_code(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
        if(t(1))
           t=mod((t+G),2);
        end
        t=[t(2:n-k+1);C1(jj,ii)];
    end
    if(1==t(1))
        t=mod((t+G),2);
    end
    C1(k+1:n,ii)=t(2:n-k+1);%2-9,共8位
end
C=reshape(C1,1,n*len2/k);