concoder.m

%Random bits
mx=ceil(rand(1,500)-0.5);
%Initialize delay line
dline=zeros(1,7);
for i=1:500
    dline(2:7)=dline(1:6);
    dline(1)=mx(i);
    xy((i-1)*2+1)=xor(xor(xor(dline(1),dline(2)),xor(dline(3),dline(4))),dline(7));
    xy(i*2)=xor(xor(xor(dline(1),dline(3)),xor(dline(4),dline(6))),dline(7));
end
%Punctuation
for i=1:166
    xyp((i-1)*4+1)=xy((i-1)*6+1);
    xyp((i-1)*4+2)=xy((i-1)*6+2);
    xyp((i-1)*4+3)=xy((i-1)*6+4);
    xyp(i*4)=xy((i-1)*6+5);
end