%Spektrumalakítás cepstrum szűréssel
%első indítás előtt a parancs sorba: l=0;
load fon; load e1; load e2; load e3; load e4; load e5;
Fh=22050;
l=l+1;
ho=512;
figure(1)
clf
hold on
kod(l)
xx=fon(l,:)';
%xx=e1';
for n=1:ho/2:length(xx)-ho
x=xx(n:n+ho-1);

x=x-mean(x);
x=x.*hamming(ho);
X=fft(x);

a=lpc(x,24);
[h,w]=freqz(1,a);
%plot(w(1:160)/pi*256,8*abs(h(1:160)))

%plot(abs(X(1:80)))
%hold on
%			X=cbm(x);
PS=(abs(X)).^2;
LP=log(PS+eps);
%XC=rceps(x);
XC=ifft((LP));

j=60;
XC(j:ho+1-j)=zeros(1,ho+2-2*j);
%XC(j:ho)=zeros(1,ho-j+1);

xi=fft((XC));
xi=exp((xi));
%xi=sqrt(xi);
%figure(2)
plot((1:2:160)/512*11025,1*abs(xi(1:80)),'k'),grid
end