CRC.m

gx=[1 0 0 0 0 0 1 0 0 1 1 0 0 0 0 0 1 0 0 0 1 1 1 0 1 1 0 1 1 0 1 1 1];
x=ceil(2*rand(1,8000)-1);
xs=~x(1:33);
for i=1:length(x)-33
    if xs(1)~=0
        xs=xor(xs,gx);
    end
    xs(1:32)=xs(2:33);
    xs(33)=x(i+33);
end
fcs=~xs