【第3章】离散傅里叶变换(DFT).ppt

第三章离散傅里叶变换(DFT) 主要内容 离散傅立叶变换的定义 离散傅立叶变换的性质 频率域采样定理 DFT的应用 3.1.2 DFT和ZT的关系 结论: 序列x(n)的N点DFT就是x(n)的ZT在单位圆上的N点等间隔采样. 序列x(n)的N点DFT就是x(n)的FT在区间[0, 2p]上的N点等间隔采样. 3.1.3 DFT的隐含周期性 在DFT的变换对中，x(n)与X(k)均为有限长序列，由于 的周期性，使DFT和IDFT中的X(k)隐含周期性，且周期均为N. 对任意整数m，总有： 在DFT表达式中，X(k)满足： 同理在IDFT表达式中，x(n)满足： 式中： N点长序列x(n)的DFT，正好是x(n)的周期延拓x((n))N 的DFS系数的主值序列. 结论： 3.2 DFT的基本性质 1. 线性性质 如果x1(n)和x2(n)是两个长度分别为N1和N2的序列，且 y(n)=ax1(n)+bx2(n)，式中a、b为常数，取N=max[N1, N2]，则y(n)的N点DFT为：  Y(k)=DFT[y(n)]=aX1(k)+bX2(k)，0≤k≤N-1 其中X1(k)和X2(k)分别为x1(n)和x2(n)的N点DFT. (2) 时域循环移位定理 设x(n)是N点序列，y(n)为x(n)的循环移位，即y(n)=x((n+m))NRN(n)，则 其中：X(k)=DFT[x(n)]， 0≤k≤N-1 (3) 频域循环移位定理 设X(k)=DFT[x(n)]，0≤k≤N-1，Y(k)= X((k+l))NRN(n)，则 其中：x(n) =IDFT[X(k)] 上式所表示的运算为x1(n)与x2(n)的循环卷积. 循环卷积与线性卷积不同. 两个N点长序列的循环卷积长度仍为N. 【说明】 X(z)在单位圆上的N点等间隔采样X(k)的IDFT为原序列x(n)以N为周期的周期延拓序列的主值序列. 如果序列x(n)的长度为M，则只有当频域采样点数N≥M时，才有 xN(n)=IDFT[X(k)]=x(n) 即可由频域采样X(k)恢复原序列x(n)，否则产生时域混叠现象. 结论(频域采样定理) 【例题】 用频域采样X(k)表示X(z)的内插公式和内插函数 设序列x(n)长度为M，在频域0~2p之间等间隔采样N点，N≥M，则有 式中： 当z=ejw时，式子就成为x(n)的傅里叶变换X(ejw)的内插函数和内插公式。 即 化简得： 3.4 DFT的应用 3.4.1 用DFT计算线性卷积 则由时域循环卷积定理有 可见，循环卷积可以按照下图所示的计算框图在频域计算. 由于DFT有快速算法FFT，因而常用DFT(FFT)计算循环卷积. 实际中，往往需要计算两个序列的线性卷积，而DFT只能直接用来计算循环卷积，为此需要导出线卷积和循环卷积之间的关系. 若h(n)和x(n)是长度为N和M的有限长序列. 则 对照线性卷积表达式，可以看出上式中： 结论: 循环卷积yc(n)等于线性卷积yL(n)以L为周期的周期延拓序列的主值序列. yL(n)的长度为N+M-1，因此只有当循环卷积长度L≥N+M-1时，yL(n)以L为周期进行周期延拓才无混迭现象. 循环卷积等于线性卷积的条件是L≥N+M-1. 下图画出h(n)、x(n)、 h(n)*x(n)和L取6、8、10时的h(n)⊙x(n)的波形. 由于h(n)长度N=4，x(n)长度M=5，N+M-1=8，所以只有当L≥8时，h(n)⊙x(n)的波形才与h(n)*x(n)相同. 【例题】 线性卷积与循环卷积 如果取L=N+M-1，则可以用DFT来计算线性卷积，计算框图如下: 当两个序列长度相差很大时，应用上述方法时要求对短序列补充很多零点. 长序列必须全部输入后才能计算，因此要求存储容量大，计算时间长. 解决的方法是采用分段卷积的办法. 分段处理的方法有重叠相加法和重叠保留法. 【说明】 设序列h(n)长度为N，x(n)为无限长序列. 将x(n)均匀分段，每段长度取M，则 (1)重叠相加法 计算h(n)和x(n)的线性卷积时，可先进行分段线性卷积yk(n)=h(n)*xk(n)，然后把结果叠加. 每一段卷积yk(n)长度为N+M-1，因此yk(n)和yk+1(n)有N-1个点重叠，必须把重叠部分的yk(n)和yk+1(n)相加才能得到完整的卷积序列y(n) . 可以用快速算法(FFT)计算分段卷积. 结论: 重叠相加法卷积示意图 【例题】 (2)重叠保留法 【例题】 设连续信号xa(t)持续时间为Tp，最高频率为fc，xa(t)的FT为 1.用DFT对连续信号进行谱分析 3.4.2 用DFT对信号进行谱分析 对xa(t)以