倒频谱(Cepstrum)或称为倒谱、二次谱和对数功率谱
工具/原料
matlab
方法/步骤
1、意义:该方法方便提取、分析原频谱图上肉眼难以识别的周期性信号,能将原来频谱图上成族的边频带谱线简化为单根谱线,受传感器的测点位置及传输途径的影响小
2、计算方法:信号求功率谱+求对数+求傅里叶逆变换
3、倒频谱擢爻充种的求解函数为rceps(实倒频谱),在MATLAB的帮助文档中,rceps的计算公式为real(if熠硒勘唏ft(log(abs(fft(y))))),即信号→频谱→对数→傅里叶逆变换,而倒频谱的定义表述中却是信号→功率谱→对数→傅里叶逆变换。即功率谱被换成了频谱。私以为是因为功率谱为频谱值的平方,在取对数后平方会变成系数2,对后续计算影响不大,因而可以近似认为结果相同。
4、在仿真中要看出倒频谱的作用,需要手动生成一组调制信号。下列程序将高频(主频为50/龀音孵茧100/200Hz傧韭茆鳟)和低频(主频为5/10/20Hz)信号进行调制,分别画出低频、高频和调制信号的时域图和频谱图。在图ftt_y中可以看到边频带的形成。(边频带相关概念参见信号的各种频域分析方法的理解(频谱、能量谱、功率谱、倒频谱、小波分析)
5、需要注意的是横轴的计算方法,单位为时间和频率如何对应上?