参考:http://www.360doc.com/content/17/0330/15/908538_641424541.shtml
由于理论和实际之间的差异,单边、双边频谱间的相互转换是信号处理中的一个基本运算。本文档中首先说明单边频谱和双边频谱相互转换的一般过程,之后列出其中若干容易出错的环节,最后以一个简单的例子进行说明。
1. 转换的基本过程
(1) 双边谱到单边谱
对于长度为N的1维实序列
其离散傅立叶变换(DFT)定义为:
其中,k=0,1,…,N-1。设采样频率为fs,则序列{Y}的频率分辨率(间隔)为:
对应于序列{y}的DFT序列{Y}长度也为N,表示为:
根据傅立叶变换的理论,这个序列中的后半部分实际上表示的是负频率信息。实序列的傅立叶变换的元素间存在关系(共轭add by lijil):
其中,常数s为:
这样,当N是奇数时,可以将序列{Y}表示为:
单边谱是长度为(N+1)/2的序列
当N是偶数时,序列{Y}可以表示为:
单边谱是长度为(N/2)+1的序列
其中,式(8)和式(10)中的w1和w2是不为0的权系数,关于不同的频谱类型在不同的参考文献中取值不一。
(2) 从单边谱到双边谱
由从双边谱到单边谱的转换过程可以看到,单边谱序列长度是奇数还是偶数与双边谱序列长度的奇偶属性没有确定的联系,所以在没有关于双边谱长度的先验知识的情况下,认为从单边谱得到的奇数长度的双边谱和偶数长度的双边谱都是正确的。
假设现有长度为N的单边谱序列{X}
可按如下方式得到长度为2N-1(奇数)的双边谱{Y}
按如下方式得到长度为2(N-1)(偶数)的双边谱{Y}
其中,式(12)和式(13)中的权系数与式(8)和式(10)中权系数相对应。
理论上,将式(12)或式(13)中的序列{Y}进行逆傅立叶变换,就可以近似得到式(1)中的实序列{y}。