离散时间傅里叶变换
定义
一组离散的实数或复数: x [ n ]( n 为所有整数)的离散时间傅里叶变换是产生以频率为变量的周期函数的一个傅里叶级数。当频率变量 ω 的单位是归一化的 弧度/样本 时,周期为 2π,而傅里叶级数为:
此频率域函数的性质源于 泊松求和公式 ( 英语 : Poisson summation formula ) 。令 X(f) 为任意函数 x(t) 的傅里叶变换,采样间隔为 T( 秒 ),等价于序列 x[n](或与之成正比),即 T ⋅ ⋅ --> x ( n T ) = x [ n ] {\displaystyle T\cdot x(nT)=x[n]} 。则以傅里叶级数表示的周期函数是 X(f) 的周期求和。赫兹以赫兹( 周期/秒 )为单位的频率 f {\displaystyle \textstyle f} 的话就会是:
图一. 傅立叶变换(左上)和左下的其周期求和(DTFT)的图示。右下角显示了用离散傅里叶变换(DFT)计算DTFT的采样。
整数 k 的单位为 转/样本 ,采样频率是 1/ T , f s ( 样本/秒 )。因而 X 1/ T ( f ) 含有移位 f s 倍数赫兹了的 X ( f ) 的精确副本,并加和在一起。对于足够大的 f s ,可以在区间 [− f s /2, f s /2] 有很少或没有失真(混叠)地观察到 k = 0 项。在图1中,该左上角分布的末端在左下图中被周期求和的混叠遮盖住了。
我们还注意到 e − − --> i 2 π π --> f T n {\displaystyle e^{-i2\pi fTn}} 是 δ δ --> ( t − − --> n T ) {\displaystyle \textstyle \delta (t-nT)} 的傅里叶变换。因此,DTFT的另一个定义为:
调制的 狄拉克梳状函数 ( 英语 : Dirac comb ) 是一个数学抽象,有时被称为 脉冲采样 。
频谱的周期性与混叠
频谱周期性
F D T F T ( e i ω ω --> T ) {\displaystyle F_{DTFT}(e^{i\omega T})} 具有周期性:
显然有: F D T F T ( e i ω ω --> T ) = F D T F T ( e i ( ω ω --> + 2 π π --> / T ) T ) {\displaystyle F_{DTFT}(e^{i\omega T})\,\!=F_{DTFT}(e^{i(\omega +2\pi /T)T})}
频谱混叠
根据DTFT的定义,有
即, f(nT) 的DTFT是 f(t) 的傅里叶变换以Ω为周期的延拓,这也从另一个角度证明了DTFT的周期性。很显然,如果 f(t) 的频谱带不限于Nyquist间隔([-Ω/2, Ω/2]), f(nT) 的DTFT必然发生混叠(aliasing),如右图所示。混叠使得信号的低频部分被高频部分“污染”,造成信号的失真。为避免这种情况,通常在进行进一步的数字信号处理之前要对采样序列进行抗混叠滤波(anti-aliasing filtering),这一处理通常是由低通滤波器除去高频分量实现的。
DTFT与DFT
DFT(离散傅里叶变换)是对离散周期信号的一种傅里叶变换,对于有限长信号,则相当于对其周期延拓进行变换。在频域上,DFT的离散谱是对DTFT连续谱的等间隔采样。
DTFT与DFT频率分辨率
图2 DTFT与DFT。上图为10点DFT,下图为补零到12点的DFT。
N 点序列 f(n) (n=0, ... ,N-1)的DFT离散谱对应于对 f(nT) 连续谱(即DTFT)的N点采样,因此DFT的频率分辨率 Δ Δ --> ω ω --> = 2 π π --> / N {\displaystyle \Delta \omega =2\pi /N} 。为了提高频率分辨率,可以考虑增加在DTFT频域上的采样点数,对偶在时域就是增加对时域信号 f(n) 的采样数。对于有限长信号 f(n) ,在时刻 0 至 N-1 以外的值实际上是已知的——都为 0 。因此,只要在序列 f(n) 前后 补零 就增加了在时域的采样,假设在f(n)前后补上 M-N (其中 M>N )个零,则补零之后序列的DFT的频率分辨率就相应提高到 Δ Δ --> ω ω --> = 2 π π --> / M {\displaystyle \Delta \omega =2\pi /M} 。相关证明如下:
File:DFT 20 points.jpg 图3 20点采样的DFT能够分辨 f 1 与 f 2
通过上式可以清楚的看到, f(n) 补零之后的DFT增加了在 f(nT) 连续频域谱上的采样。采样点数从N增加到M,从而提高了DFT频谱的分辨率。另一方面,补零之后在频域采样的位置发生了变化,因此可以观察到其他的频点。
如图2所示,对信号 f ( t ) = 0.5 cos --> ( 2 π π --> f 1 t ) + cos --> ( 2 π π --> f 2 t ) {\displaystyle f(t)=0.5\cos(2\pi f_{1}t)+\cos(2\pi f_{2}t)} (其中 f 1 = 300 H z {\displaystyle f_{1}=300Hz} , f 2 = 250 H z {\displaystyle f_{2}=250Hz} )按照 f s = 1000 H z {\displaystyle f_{s}=1000Hz} 采样。 f 1 {\displaystyle f_{1}} 与 f 2 {\displaystyle f_{2}} 对应的数字角频率分别为 ω ω --> 1 = 0.6 π π --> {\displaystyle \omega _{1}=0.6\pi } , ω ω --> 2 = 0.5 π π --> {\displaystyle \omega _{2}=0.5\pi } 。上图为取10个采样点作DFT,可在0.6π处看到对应的频率分量,然而由于采样点少,看不到0.5π处的分量。下图为补2个零之后的DFT离散谱,可以见到离散谱的分辨率提高到了π/6,而且能够观察到10点DFT无法看到的0.5π频率分量。另外,虚线为DTFT连续谱,可见,DFT确实是在频域对DTFT的采样。
但是我们可以看到,即使是10点DTFT的连续谱也不能分辨 f 1 {\displaystyle f_{1}} 与 f 2 {\displaystyle f_{2}} (只有一个峰)。这是因为10点DTFT的分辨率为 f s /10 = 100Hz ,大于 f 1 - f 2 = 50 Hz 。所以,只有采样的点数超过20(即分辨率小于50Hz),才能分辨出 f 1 {\displaystyle f_{1}} 与 f 2 {\displaystyle f_{2}} 这两个频率分量(如图3所示)。而前面提到的对有限长信号补零作DFT以提高频谱分辨率的说法,也只是针对在DTFT连续谱上采样而言,只有增加采样点数和提高采样频率才能真正提高离散谱的分辨率。
以DFT近似DTFT
前面提到,在时间序列前后补零之后作DFT可以增加在DTFT上的采样点数。可以想见,如果补上无穷多个零,则可以得到无穷多个DTFT连续谱上的采样点,从而以DFT逼近DTFT。即,使得离散谱的分辨率足够小,即为连续谱。
这种想法显然是错误的,首先应该注意到,DFT的周期型,虽然进行DFT的序列x(k)是有限长序列,但是DFT隐含对x(k)的周期延拓,x(k)是信号序列,周期延拓后还是原来的信号,如果对x(k)进行补零后,1.补零后序列长度不是原信号序列长度的整数倍,那么补零后的序列就不再是原信号的序列,已经破坏了原信号的信息。2.补零后序列长度是原信号序列的整数倍,原信号信息没有被破坏,仍能通过DFT后获得原信号信息,但是补零后序列的DFT虽然频率步长小了,可以观察的频谱信息更丰富了,但是这些额外的频谱并不是原信号的频谱,而是窗函数的频谱信息,也就是包络线是窗函数的形状。所以补零并没有给出关于原信号频谱的更多信息。补零并不能使补零后的DFT结果的包络线和原信号频谱的包络线一致。
DTFT与Z变换
离散时间傅里叶变换可以被看作Z变换的特例。Z变换被定义为:
如果在z平面的单位圆( z = e i ω ω --> {\displaystyle z=e^{i\omega }} )上对信号 f ( n ) {\displaystyle f(n)} 做Z变换:
此即为 f ( n ) {\displaystyle f(n)} 的离散时间傅里叶变换。因此通常用 F ( e i ω ω --> ) {\displaystyle F(e^{i\omega })} ,而不是 F ( ω ω --> ) {\displaystyle F(\omega )} 表示DTFT。
参看
离散傅里叶变换
Z变换
傅里叶级数
采样定理
参考文献
Alan V. Oppenheim and Ronald W. Schafer : Discrete-Time Signal Processing , Prentice Hall Signal Processing Series, ISBN 978-0-13-754920-7
Boaz Porat : A Course in Digital Signal Processing , pp. 27-29, pp. 104-105, John Wiley and Sons, ISBN 978-0-471-14961-3
Sophocles J. Orfanidis : Introduction to Signal Processing , Prentice Hall International, Inc., ISBN 7-32-03059-6
免责声明:以上内容版权归原作者所有,如有侵犯您的原创版权请告知,我们将尽快删除相关内容。感谢每一位辛勤著写的作者,感谢每一位的分享。
- 有价值
- 一般般
- 没价值
{{item.userName}} 举报
{{item.time}} {{item.replyListShow ? '收起' : '展开'}}评论 {{curReplyId == item.id ? '取消回复' : '回复'}}