请教高手，帮忙看下我的求倒谱过程是否有误！

zzz417

本人实验采样故障信号频率fs=10000,共采集了10000个数据点进行分析，用MATLAB的倒谱函数求得的效果不理想，现在用倒谱定义求的程序如下，y1为故障信号函数，但还是在50HZ即20毫秒处的峰值不明显，请问我的程序是否编写有误。是否有参考程序帮帮我，谢谢！
nfft=10000;
nn=1:nfft/2;
t=1000*(0:1/fs:1-1/fs);%化成毫秒
w=hanning(nfft);      
zw=psd(y1,nfft,fs,w,nfft/2); %求功率谱密度函数
zzw=real(ifft(log(abs(zw)))); %根据定义求倒谱
figure(5);
grid on
plot(t(nn),abs(zzw(nn)));
xlabel('时间(s)');
ylabel('幅值');
grid on;