function cepstrum(y)

Yk=fft(y);
figure;
plot(abs(Yk));
title('Yk');

logYk=log10(abs(Yk)+0.0001);
figure;
plot(logYk);
title('log(Yk)');

logC=ifft(logYk);
figure;
plot(abs(logC));
title('Cepstrum');

win=zeros(length(logC),1);
winN=50;
w=boxcar(winN);
win(1:winN/2)=w(winN/2+1:end);
win(length(logC)-winN/2+1:end)=w(1:winN/2);

winlogYkk=win.*logC;
%winlogYkk=logYkk;
hold on;
plot(abs(winlogYkk),'r');

winlogYk=ifft(winlogYkk);
figure;
plot(real(winlogYk));
title('log(Yk) windowed');

winYk=10.^winlogYk;
figure;
plot(abs(winYk));
title('Yk windowed');

winYkc=abs(winYk)+j*imag(Yk);
figure;
spect(winYkc,44100,'r');
spect(Yk);
%spect(fft(gl),44100,'g');
%spect(fft(fi),44100,'k');
title('Ergebnisse: bl:eingang, gn:glottis, k:filter, rot:Ausgang');

out=real(ifft(winYkc));
figure;
hold on;
plot(y,'b');
%plot(gl,'g');
%plot(fi,'k');
plot(out,'r');
title('Ergebnisse: bl:eingang, gn:glottis, k:filter, rot:Ausgang');
