DSP第02章 离散傅里叶变换(DFT).pptVIP

* 2.3.2 循环卷积定理 若x1(n) X1(k)， x2(n) X2(k)，则 x1(n)*x2(n) X1(k)?X2(k) 或 x1(n)?x2(n) 1/N(X1(k)*X2(k)) 。 卷积定理指出了相乘与卷积运算的关系，同时给出了卷积运算的另一途径，即 由卷积定理也可以推论出x1(n)*x2(n) 的周期为N。 DFT DFT x1(n) x2(n) X1(k) X2(k) x1(n)*x2(n) X1(k)?X2(k) IDFT * 2.3.3 循环卷积与线性卷积的关系 设x1(n)与x2(n)是宽度分别为N1、N2的有限时宽序列，则 循环卷积： 线性卷积： 其中，循环卷积长度 N3=max{N1，N2} 线性卷积长度 N4=N1+N2-1 显然，N3 N4 可以证明 也即循环卷积是线性卷积以循环卷积长度进行周期化并叠加（时域混叠）的结果。 * 2.3.4 利用循环卷积求解线性卷积 实际信号处理中，特别是分析处理模拟系统时，我们经常需要获得线性卷积，而利用计算机这类数字设备只能实现循环卷积，所以需要解决用循环卷积求解线性卷积的问题。 * 两个有限长序列的线性卷积 从循环卷积与线性卷积的关系可知，如果将循环卷积的长度定义为线性卷积有效结果的长度，即取 N3≥N4=N1+N2-1 则循环卷积一个周期内的序列值与线性卷积完全相同。这样，通过周期扩展，就可以用循环卷积的计算结果表示线性卷积，称此为扩展周期法。 这可以理解为：将时域x1(n)与x2(n)补零至N4长度，导致其频谱采样间隔变小，采样精度提高，从而减小(消除)了由频谱X1(k)?X2(k)确定的时域卷积序列x1(n)*x2(n)中的混叠现象，使得循环卷积与线性卷积结果一致。 * 因此，对两个有限长序列计算线性卷积的方法为： ⅰ. 确定循环卷积的计算长度，即 N≥N1+N2-1 其中， N1与N2分别为序列x1(n)与x2(n)长度(即周期)； ⅱ. 将x1(n)与x2(n)分别补零至长度为N； ⅲ. 对x1(n)与x2(n)分别求出N点的DFT为X1(k)与X2(k)； ⅳ. 计算X1(k)?X2(k)； ⅴ. 求X1(k)?X2(k)的N点IDFT为x1(n)*x2(n)，该结果即为线性卷积。 注：算法中的DFT与IDFT均可利用FFT求解； ⅲ ~ ⅴ 也可以通过直接计算N长度的离散卷积所替代。 * 例：已知h=[1,2,-1,1]， x=[1,1,2,1,2,2,1,1]， 求y=h*x。 解： y=[1,3,3,5,3,7,4,3,3,0,1] Ly=8+4-1=11 x0 x1 x2 x3 x4 h0 h0x0 h0x1 h0x2 h0x3 h0x4 h1 h1x0 h1x1 h1x2 h1x3 h1x4 h2 h2x0 h2x1 h2x2 h2x3 h2x4 h3 h3x0 h3x1 h3x2 h3x3 h3x4 线性卷积表 h\x 1 1 2 1 2 2 1 1 1 1 1 2 1 2 2 1 1 2 2 2 4 2 4 4 2 2 -1 -1 -1 -2 -1 -2 -2 -1 -1 1 1 1 2 1 2 2 1 1 * 显然，当N确定时， [AN]与[AN*]为常数矩阵，只要给定[x(n)]或[X(k)]就可以通过矩阵计算出[X(