分类:
2012-07-18 18:32:10
连续时间信号x(t) 以及对应的都是连续函数。由于数字系统只能处理有限长的,因此必须将x和都离散化,并且建立对应的傅里叶变换。
假设x(t)时限于[0, L],再通过时域采样将x(t)离散化,就可以得到有限长离散信号,记为xdiscrete(t)。设采样周期为T,则时域采样点数N=L/T。
它的傅里叶变换为
这就是x(t)在时域采样后的连续傅里叶变换,也就是,它在频域依然是连续的。
下面将频域信号转化为有限长离散信号。与对时域信号的处理类似,假设频域信号是带限的,再经过离散化,即可得到有限长离散信号。依据,时域采样若要能完全重建原信号,频域信号应当带限于(0,1/T)。由于时域信号时限于[0, L],由采样定理以及时频对偶的关系,频域的采样间隔应为1/L。故,频域采样点数为:
即频域采样的点数和时域采样同为N,频域采样点为 在DTFT频域上采样:
令T=1,将其归一化,就得到前面定义的离散傅里叶变换。因此,DFT就是先将信号在时域离散化,求其连续傅里叶变换后,再在频域离散化的结果。
DFT与CT
下面考察离散傅里叶变换与连续傅里叶变换(CT,Continuous Fourier Transform)的关系。连续傅里叶变换
的采样为:
将这个积分以黎曼和的形式近似,有:
这一结论指出离散傅里叶变换确实是连续傅里叶变换在某种意义上的近似。不过必须注意以下两点:
为此,通常对连续信号的时域采样再做一次加窗处理(Windowing),这样就得到带限的有限长离散信号。
离散时间傅里叶变换(DTFT)是在时域上对的采样。DFT则是在频域上对DTFT的均匀采样。离散信号x[n](n=0,...,N-1)的DTFT为:
对在离散的频点上采样
即为x 的DFT。由于DTFT在频域是周期的,所以在DTFT频域上的均匀采样也应是周期的。实际上是这个周期序列的主值序列。
N 点序列的DFT只能在有限的N个频点上观察频谱,这相当于从的缝隙中观察景色,对于了解信号在整个频域上的特性是不够的。为了观察到其他频率的信息,需要对原信号x[n]做一些处理,以便在不同的频点上采样。
将原来在DTFT频域上的采样点数增加到M 点,这样采样点位置变为。则对应的DFT成为
若在序列x[n] 之后补上M-N个零,设为x'[n],则上式变为
因此将x[n]补零再做DFT就可以得到x[n]的DTFT在其他频率上的值,相当于移动了栅栏,因而能够从其他位置进行观察。
N 点DFT的频谱分辨率是2π / N。一节指出可以通过补零观察到更多的频点,但是这并不意味着补零能够提高真正的频谱分辨率。这是因为x[n] 实际上是x(t) 采样的主值序列,而将x[n]补零得到的x'[n] 周期延拓之后与原来的序列并不相同,也不是x(t) 的采样。因此与是不同离散信号的频谱。对于补零至M点的x'的DFT,只能说它的分辨率2π / M仅具有计算上的意义,并不是真正的、物理意义上的频谱。频谱分辨率的提高只能在满足采样定理的条件下增加时域采样长度来实现。
周期为N的离散信号构成一个N 维。在这一空间上两个信号x 和y 的定义为
下面给出上的一组:
将信号x 在这组正交基上分解,得
令
此即为离散傅里叶变换。又
则有
此即为离散傅里叶变换的逆变换。
由此可知,在正交分解的角度上说,离散周期信号x的离散傅里叶变换实质上是x在正交基{ek}上的分量。而从的角度上说,{ek}是的,则是对应的。
根据,离散信号x与y的对偶于频域上x与y离散傅里叶变换的乘积。下面揭示DFT与圆周卷积的内在关系。
对长为N的离散信号与,如果要计算它们的卷积,就必须定义与在0≤n
如此,周期为N的圆周卷积:
卷积的结果仍然是以N为周期的离散信号。
前面指出,ek是DFT的特征向量,实际上它也是圆周卷积的特征向量。定义x与y的圆周卷积算子:
则ek与y的圆周卷积为
显然,y的DFT
就是圆周卷积的特征值。
离散傅里叶变换是可逆的
其中C表示。即,任意N-维复向量都存在DFT和IDFT,而且其DFT和IDFT也是N-维复向量。
向量组exp(2πi kn/N)是N-维复数空间上的一组正交基:
其中 δkn是。
时域信号序列xn的相位移动exp(2πinm / N)(其中m为整数)在频域反映为“循环移位”,即:频域信号序列Xk变为,其中下标是指将k-m 对N 。类似的,由对偶性,时域信号序列的循环移位对应于频域信号序列的相移:
上文中一节已经证明,离散序列的傅里叶变换是周期的。有限长序列xn的离散傅里叶变换Xk可以被看作是它的周期延拓序列的离散时间傅里叶变换。由时频对偶性可知Xk也可以被看作它的周期延拓的主值。
上一节的移位定理隐含着DFT的周期性,因为DFT的幅度 | Xk | 不受输入信号的循环移位的影响,因为时移在频域对偶于相移,循环移位仅仅使DFT的相位产生平移。周期性的边界条件在DFT的许多应用中有重要作用。解时,就假设边界条件是满足周期性的,这是一个很有用的性质(参见)。
如果 Xk 和 Yk 分别是 xn 和 yn 的离散傅立叶变换,那么就有 :
其中星号表示复共扼。 是普朗歇尔定理的一个特例:
DFT在诸多多领域中有着重要应用,下面仅是颉取的几个例子。需要指出的是,所有DFT的实际应用都依赖于计算离散傅里叶变换及其逆变换的快速算法,即。
快速傅里叶变换(即FFT)是计算离散傅里叶变换及其逆变换的快速算法。按照DFT的定义计算一个长为n的序列的DFT需要的计算复杂度达到了,而同样长度FFT的计算复杂度仅为。由于DFT的逆变换可以由DFT表示,所以DFT逆变换的计算同样可以由FFT完成。FFT算法的提出,使DFT得到了广泛的实际应用。
前面指出,。因此可以对连续信号x(t)均匀采样并截断以得到有限长的离散序列,对这一序列作离散傅里叶变换,可以分析连续信号x(t)频谱的性质。前面还提到DFT应用于频谱分析需要注意的两个问题:即采样可能导致信号混叠和截断信号引起的频谱泄漏。可以通过选择适当的采样频率(见)消减。选择适当的序列长度并加窗可以抑制频谱泄漏。
由于人类感官的分辨能力存在极限,因此很多有损压缩算法利用这一点将语音、音频、图像、视频等信号的高频部分除去。高频信号对应于信号的细节,滤除高频信号可以在人类感官可以接受的范围内获得很高的压缩比。这一去除高频分量的处理就是通过离散傅里叶变换完成的。将时域或空域的信号转换到频域,仅储存或传输较低频率上的系数,在解压缩端采用逆变换即可重建信号。
离散傅里叶变换及其多维形式在偏微分方程的求解中也有应用。此时DFT被看作的近似。傅里叶级数将函数在复指数einx上展开,这正是微分算子的特征方程:d/dx einx = in einx。因此,通过傅里叶级数的形式,线性常微分方程被转换为代数方程,而后者是很容易求解的。此时得到的结果是偏微分方程解的级数表示,只要通过DFT逆变换即可得到其一般表示。这种方法被称作谱方法或级数解法。
目前长或最快速的算法是基于离散傅里叶变换的。由于整数(或多项式)乘法是逐位(或逐项)乘累加的形式,因此整数(或多项式)乘积的数字(或系数)可以用乘数数字(或乘式系数)的表示。利用,只要将数字(或系数)序列通过离散傅里叶变换变到频域,就可以将逐个乘累加的卷积变为对位的乘法,从而减少计算量,再以一次逆变换便可以得到乘法结果。需要注意整数乘法还有的问题。下面以多项式乘法为例说明这一应用:
假设需要计算多项式乘法c(x) = a(x) · b(x),这一乘积结果的各项系数的表达式实际上是线性卷积的形式。由于离散傅里叶变换对应于圆周卷积,因此需要将这两个乘式的系数按照次数升序排列,序列后“补零”,使它们的长度d 大于两式最高项次数的和:d > deg(a(x)) + deg(b(x))。然后作圆周卷积:
其中c 就是c(x) 系数的向量。由于圆周卷积的DFT是乘积,有:
则
利用快速傅里叶变换,这一算法的运算复杂度为 。
OFDM(正交频分复用)在宽带中有重要的应用。这种技术将带宽分为N个等间隔的子载波,可以证明这些子载波相互正交。尤其重要的是,OFDM调制可以由IDFT实现,而解调可以由DFT实现。OFDM还利用DFT的移位性质,在每个帧头部加上循环前缀(Cyclic Prefix),使得只要信道延时小于循环前缀的长度,就能消除信道延时对传输的影响。