且
H(k)?DFT[h(n)]LX(k)?DFT[x(n)]L 0?k?L?1
则由DFT的循环卷积定理有
Yc(k)?DFT[yc(n)]L?H(k)X(k) 0?k?L?1
4 复共轭序列的DFT
性质:设x(n)是x(n)的复共轭序列,长度为N,X(k)?DFT[x(n)]N, 则DFT[x(k)]N?X(N?k)***0?k?N?1
例:给定一16-点实序列x(n), 其 16-点 DFT 记为X(k), 已知X(13)= 2 + j3,则X*(3) = ( 2 + j3 )。 说明:DFT的性质。实序列的DFT的共轭对称性:X(k) =X*(N-k),或X(N-k) =X*(k)。
5 DFT的共轭对称性(重点)
可总结出DFT的共轭对称性质:如果序列x(n)的DFT为X(k),则x(n)的实部和虚部(包括j)的DFT分别为X(k)的共轭对称分量和共轭反对称分量;而x(n)的共轭对称分量和共轭反对称分量的DFT分别为X(k)的实部和虚部乘以j。
3.4 频域采样定理
离散傅里叶变换相当于信号傅里叶变换的等间隔采样,也就是说实现了频域的采样,便于计算机计算。那么是否任一序列都能用频域采样的方法去逼近呢?这是一个很吸引人的问题。 我们考虑一个任意的绝对可和的序列x(n),它的z变换为
如果对X(z)单位圆上进行等距离采样
现在要问,这样采样以后,信息有没有损失?或者说,采样后所获得的有限长序列xN(n)能不能代表原序列x(n)。
为了弄清这个问题,我们从周期序列
开始
第36页 共60 页
由于
所以
也即
是原非周期序列x(n)的周期延拓序列,其时域周期为频域采样点数N。在第一章我们看
到,时域的采样造成频域的周期延拓,这里又对称的看到,频域采样同样造成时域的周期延拓。 因此,如果序列x(n)不是有限长的,则时域周期延拓时,必然造成混叠现象,因而一定会产生误差。 对于长度为M的有限长序列,只有当频域采样点数N大于或等于序列长度M时,才有
即可由频域采样值X(k)恢复出原序列x(n),否则产生时域混叠现象,这就是所谓的频域采样定理。(重点)
1?z?NX(z)?N内插公式:
X(k)??k?1k?01?WNz
N?13.5 DFT的应用举例 1.用DFT计算线性卷积(重点)
用循环(周期)卷积计算有限长序列的线性卷积(重点)
对周期要求:N?N1?N2?1(N1、N2分别为两个序列的长度)(记住!!)
2.用DFT进行谱分析的误差问题(重点)
(1)混叠现象
利用DFT逼近连续时间信号的傅里叶变换,为避免混叠失真,按照抽样定理的要求,采样频率至少是信号最高频率的两倍。
解决混叠问题的唯一方法是保证采样频率足够高。 (2)截断效应
任何带限信号都是非时限的,任何时限信号都是非带限的。实际问题中遇到的离散时间序列可能是非时限的、无限长序列,在对该序列利用DFT进行处理时,由于作DFT的点数总是有限的,因此就有一个必须将该序列截断的问题。序列截断的过程相当于给该序列乘上一个矩形窗口函数RN(n)。如果原来序列的频谱为
,矩形窗函数的频谱为
,则截断后有限长序列的频谱为
第37页 共60 页
截断后序列的频谱与原序列频谱必然有差别,这种差别对谱分析的影响主要表现在如下两个方面: ①频谱泄露:由于矩形窗函数频谱的引入,使卷积后的频谱被展宽了,即
的频谱“泄露”
到其它频率处,称为频谱泄露。
在进行DFT时,由于取无限个数据是不可能的,所以序列的时域截断是必然的,泄露是难以避免的。为了尽量减少泄露的影响,截断时要根据具体的情况,选择适当形状的窗函数,如汉宁窗或汉明窗等。
②谱间干扰。 在主谱线两边形成很多旁瓣,引起不同频率分量间的干扰(简称谱间干扰),特别是强信号谱的旁瓣可能湮没弱信号的主谱线,或者把强信号谱的旁瓣误认为是另一频率的信号的谱线,从而造成假信号,这样就会使谱分析产生较大偏差。 (3)栅栏效应
由于DFT是有限长序列的频谱等间隔采样所得到的样本值,这就相当于透过一个栅栏去观察原来信号的频谱,因此必然有一些地方被栅栏所遮挡,这些被遮挡的部分就是未被采样到的部分,这种现象称为栅栏效应。由于栅栏效应总是存在的,因而可能会使信号频率中某些较大的频率分量由于被“遮挡”而无法得到反映。此时,通常在有限长序列的尾部增补若干个零值,借以改变原序列的长度。这样对加长的序列作DFT时,由于点数增加就相当于调整了原来栅栏的间隙,可以使原来得不到反映的那些较大的频率分量落在采样点上而得到反映。
产生原因说明:由傅里叶变换理论知道,若信号持续时间有限长,则其频谱无限宽;若信号的频谱有限宽,则其持续时间必然为无限长。(重点)
如果用DFT分析连续信号的频谱,在对连续信号采样时,无法满足采样定理,那么就会出现频谱混叠现象。解决混叠问题的唯一方法是保证采样频率足够高。当连续信号无限长或很长时,在对连续信号采样时,采样点数太多以致无法存储和计算,需要将信号截断,这样将导致频谱的泄漏现象。为了尽量减少泄露的影响,截断时要根据具体的情况,选择适当形状的窗函数,如汉宁窗或汉明窗等。用DFT计算连续信号的频谱只能得到采样点上的频谱,而不能看到整个频谱,这种现象称为栅栏效应。可以通过增加点数,因为点数增加就相当于调整了原来栅栏的间隙,可以使原来得不到反映的那些较大的频率分量落在采样点上而得到反映。
3.用DFT进行谱分析的参数选择问题(重点)
对模拟信号频谱的采样间隔,称之为频率分辨率。
(1)在已知信号的最高频率fc(即谱分析范围)时,为了避免频率混叠现象,要求采样频率Fs满足:Fs>2fc。
(2)采样频率Fs,采样点数N,谱分辨率F=Fs/N,如果保持采样点数N不变,要提高频谱分辨率(减小F),就必须降低采样频率,采样频率的降低会引起谱分析范围变窄和频谱混叠失真。如维持Fs不变,为提高频率分辨率可以增加采样点数N。因为NT=Tp,T=Fs-1,只有增加对信号的观察时间Tp,才能增加N。
(3)采样点数N>2fc/F (4)最小记录时间Tp≥1/F 例:用DFT对实信号进行谱分析,要求频率分辨率F?50Hz,信号最高频率为fmax?1000Hz,试确定以下参数:(1)最小记录时间;(2)最大取样间隔;(3)最少采样点数;(4)若要求频率分辨率提高一倍,求最少采样点数。(重点)
解 (1)
Tpmin?1?0.02sFmax
第38页 共60 页
(2)
Tmax?1?0.5ms2fmax
(3)
Nmin?2fmaxTpmin??40FmaxTmax
2fmax2000??80Fmax25
(4)
Nmin?
简答题:(重点)
1. 一个序列的DFT与序列的傅里叶变换之间的关系是什么?
2. 序列的DTFT和序列的z变换间的关系是什么?序列的DFT和序列的Z变换间的关系是什么? 3. 有限长序列x(n)的长度为M,对其进行频域采样,不失真的条件是什么?两个有限长序列
x1(n),0?n?M,x2(n),0?n?N,对它们进行线性卷积,结果用y(n)表示,y(n)的长
度是多少?如果进行循环卷积,那么什么时候线性卷积和循环卷积的结果相等? 4. 用DFT进行谱分析带来哪些误差问题?采取什么措施可以减少这些误差? 5. 时域采样定理的要点是什么?频域采样定理的要点是什么?
第39页 共60 页
第四章:快速傅里叶变换并不是一种新的变换,它是离散傅里叶变换的一种快速算法。
4.1 直接计算DFT的问题及改进的途径
直接计算DFT,需要次复数乘法,次复数加法。
直接计算离散傅里叶变换,由于计算量近似正比于N2,显然对于很大的N值,直接计算离散傅里叶变换要求的算术运算量非常大。(重点)我们可以利用系数WNnk的特性来改善离散傅里叶变换的计算效率。
(1)的对称性
(2)的周期性
利用的对称性和周期性,将大点数的DFT分解成若干个小点数的DFT,FFT正是基于这个
基本思路发展起来的。(重点)
说明:快速傅里叶变换FFT(Fast Fourier Transform)
分类:按时间抽取(DIT)算法和按频率抽取(DIF)算法。
第40页 共60 页