clear;
D_k=round(rand(10,100));%產(chǎn)生10*100的2進(jìn)制數(shù)
d_k1=reshape(D_k,1,sum(size(D_k)));
d_k2=reshape(D_k.',1,1000);
b=[0 0 0 0];
c=[0 0 0 0];
for n=1:length(d_k1)
    x_k(1,n)=d_k1(1,n);
    a_k1=mod(b(1,1)+b(1,2)+b(1,3)+b(1,4)+d_k1(1,n),2);
    y_k1(1,n)=mod(a_k1+b(1,4),2);
    b(1,4)=b(1,3);
    b(1,3)=b(1,2);
    b(1,2)=b(1,1);
    b(1,1)=a_k1;
end
for n=1:length(d_k2)
    a_k2=mod(c(1,1)+c(1,2)+c(1,3)+c(1,4)+d_k2(1,n),2);
    y_k2(1,n)=mod(a_k2+c(1,4),2);
    c(1,4)=c(1,3);
    c(1,3)=c(1,2);
   c(1,2)=c(1,1);
    c(1,1)=a_k2;
end   
z_k=zeros(1,2*length(d_k1));
for m=1:length(d_k1)
    k=2*m;
    z_k(1,k-1)=x_k(1,m);
    z_k(1,k)=y_k1(1,m);
end
z_k