homoque.m

来自「this program deals with homomorphic anal」· M 代码 · 共 14 行

function homoQue(x,Fs,nfft)

index1 = Fs/1000;                 
index2 = Fs/50;  
% range for searching maxima


Y=fft(x.*hamming(length(x)), nfft);
C=(ifft(log(abs(Y)),1024));
plot(linspace(1,250,nfft),C);
title('Cepstrum');xlabel('Quefrency(ms)');ylabel('Amplitude');

[junk,locs]=max(abs(C(index1:index2)));
fprintf('Pitch Frequency in Quefrency Domain = %gHz\n',Fs/(index1+locs-1));