DSP-4-2014讲述.pptVIP

解： 直接由频域采样定义得 时域混叠 * * 4.4 频率域采样 【例2】已知序列 的N点离散DFT为 求序列 解：根据循环移位性质以及前例： 所以 最后 * * 4.3离散傅立叶变换（DFT）的基本性质 循环卷积定义： 设 x1(n) 和 x2(n) 都是长度为 N 的有限长序列，把它们分别拓展为周期序列 和 ，定义循环卷积为： ——周期序列卷积后取主值 4.3.3 循环卷积定理 * * 4.3离散傅立叶变换（DFT）的基本性质 因为上式的求和范围是 m 由 0 到 N-1，因此第一个序列 x1(m) 可以不作周期拓展，即 注意两个 N 点序列的线性卷积将导致一个更长的序列。而循环卷积将区间限制在 0≤n≤N-1，结果仍为 N 点序列，它与线性卷积的结构类似。不同点在于求和范围和 N 点循环移位。它与 N 有关，也叫做 N 点循环卷积。 N * * 4.3离散傅立叶变换（DFT）的基本性质 两个有限长序列循环卷积的过程： (1)上式中求和变量为m，n为参变量； (2)将x2(m) 以N为周期作周期延拓得到x2((m))N ； (3)翻转x2((m))N 形成x2((-m))N (4)对x2((-m))N进行循环移位x2((n-m))N，取主值序列，形成x2((n-m))N RN (m) ; (5) n=0,1,…N-1时， x1(m) 和x2((n-m))N R N(m)对应相乘，并对m在0～N-1区间求和。 x2 m x2(m) m x1(m) 0 0 N N m x2((0-m))N 0 N m x2((1-m))N 0 N n 0 N 注意 n=0 时的循环反转 * * 4.3离散傅立叶变换（DFT）的基本性质 循环卷积的时频映射关系 由 DTFT 的性质可知，两个序列时域上的线性卷积运算在频域上表现为两个序列 DTFT 结果的乘积。 同样的，若 则 即当在频域中进行两个 N 点 DFT 相乘时，在时域中映射为循环卷积（而不是通常的线性卷级 ！！！）。 * * 4.3离散傅立叶变换（DFT）的基本性质 证明： 直接对上式两边进行DFT 令n-m=n′， 则有 x2 * * 4.3离散傅立叶变换（DFT）的基本性质 4.3 离散傅立叶变换（DFT）的基本性质 * * 同心圆法 可用两个同心圆来表示： x1(n) ：内圆顺时针方向排列 x2(n) ：外圆逆时针方向排列 x1(0) 与 x2(0) 对齐。 * * 4.3离散傅立叶变换（DFT）的基本性质 两圆上对应数两两相乘后求和，得 x3(0) ; 将x2(n-m)移位一位，即外圆顺时针转动一位，重复（1）步骤，得x3(1) ； 依次下去，求得 * * 4.3离散傅立叶变换（DFT）的基本性质 【例】 * * 4.3离散傅立叶变换（DFT）的基本性质 例 设 x1(n) = {1,2,2}，x2(n)={1,2,3,4}，计算 4 点循环卷积 解： 注意 x1(n) 为 3 点序列，进行循环卷积之前在其尾部填一个零，使其成为 4 点序列，分别在时域和频域中求解这个问题。 * * 4.3离散傅立叶变换（DFT）的基本性质 1）时域方法(同心圆) 4 点循环卷积由下式给出 对每个 n 产生一个循环移位序列，将它的样本逐个与 x1(m) 相乘，然后求和，得此 n 值的循环卷值，在 0≤n ≤3 上重复此过程。考虑 x1(m) = {1,2,2,0} 和 x2(m) = {1,2,3,4} 注意此为 x2(0)，而不是 {4,3,2,1} * * 4.3离散傅立叶变换（DFT）的基本性质 n = 2 时， n = 3 时， 因此， * * 4.3离散傅立叶变换（DFT）的基本性质 2）频域方法： 首先计算 x1(n) 和 x2(n) 的 4 点 DFT，逐个样本相乘，取IDFT，得到循环卷积。 则（注意：X1(k) 和 X2(k) 对应位相乘） IDFT 后，得 表面看来，这样做反而更复杂，且涉及复数运算。但后面我们会看到，DFT 有快速算法 FFT，特别是当 N 很大时，效率会显著提高。 * * 4.3离散傅立叶变换（DFT）的基本性质 循环卷积也可由矩阵形式定义，如下： 则有 矩阵第1行为序列的循环反转，其它行为第1行的循环移位 * * 4.3离散傅立叶变换（DFT）的基本性质 例 序列为： x1(n) = {1,2,2}，x2(n) ={1，2，3，4}，在这个例子中，研究 N 的大小对循环卷积的影响，显然 N≥4；否则 x2(n) 将在时域