function de=demodulation(fftin)
n=length(fftin);
fout=zeros(1,n);
for j=1:n/8
    x=fftin(j*8-7:j*8);
    fout(j*8-7:j*8)=fft(x,8);
end
de=zeros(1,n);
for j=1:n
    if real(fout(j))>=0&&imag(fout(j))>=0
        de(j)=0;
    elseif real(fout(j))<0&&imag(fout(j))>=0
        de(j)=1;
    elseif real(fout(j))<0&&imag(fout(j))<0
        de(j)=2;
    else de(j)=3;
    end
end