zf=zeros(N,1);
nf=2^6;           % Nombre de fenêtres
lf=round(N/nf);   % Longueur des fenêtres
k=0;
for f=1:nf,
    zf(k+(1:lf))=fft(y(k+(1:lf)));
    k=k+lf;
end,
zs=sort(abs(zf));
seuil=zs(round(N*(1-1/r)));
zf(abs(zf)<seuil)=0;
k=0;
for f=1:nf,
    ycf(k+(1:lf))=ifft(zf(k+(1:lf)));
    k=k+lf;
end,