题目
已知两个有限长序列:[ x_1(n) = [1,2,0,1], x_2(n) = [2,1,1,-1], n = 0,1,2,3 ]1. 直接计算 4 点循环[1]卷积 y_c(n) = x_1(n) otimes x_2(n) 。2. 分别计算 x_1(n) 和 x_2(n) 的 4 点 DFT: X_1(k) 和 X_2(k) 。3. 利用频域相乘对应时域循环卷积的性质,由 Y_c(k) = X_1(k)X_2(k) 求 y_c(n) ,验证与(1)的结果一致。4. 如果要计算 x_1(n) 与 x_2(n) 的线性卷积,DFT 的长度至少需要多少点?简要说明原因。
已知两个有限长序列:
$x_1(n) = [1,2,0,1], x_2(n) = [2,1,1,-1], n = 0,1,2,3$
1. 直接计算 4 点循环[1]卷积 $ y_c(n) = x_1(n) \otimes x_2(n) $。
2. 分别计算 $ x_1(n) $ 和 $ x_2(n) $ 的 4 点 DFT:$ X_1(k) $ 和 $ X_2(k) $。
3. 利用频域相乘对应时域循环卷积的性质,由 $ Y_c(k) = X_1(k)X_2(k) $ 求 $ y_c(n) $,验证与(1)的结果一致。
4. 如果要计算 $ x_1(n) $ 与 $ x_2(n) $ 的线性卷积,DFT 的长度至少需要多少点?简要说明原因。
题目解答
答案
-
4点循环卷积计算
根据循环卷积定义:
$y_c(n) = \sum_{m=0}^{3} x_1(m) \cdot x_2((n - m) \mod 4)$
计算得:
$y_c(n) = \{1, 6, 2, 3\}$ -
4点DFT计算
- $x_1(n) = \{1, 2, 0, 1\}$ 的DFT:
$X_1(k) = \{4, 1 - j, -2, 1 + j\}$ - $x_2(n) = \{2, 1, 1, -1\}$ 的DFT:
$X_2(k) = \{3, 1 - 2j, 3, 1 + 2j\}$
- $x_1(n) = \{1, 2, 0, 1\}$ 的DFT:
-
频域相乘后IDFT验证
- 频域相乘:
$Y_c(k) = X_1(k) \cdot X_2(k) = \{12, -1 - 3j, -6, -1 + 3j\}$ - 对 $Y_c(k)$ 做4点IDFT,还原时域序列:
$y_c(n) = \{1, 6, 2, 3\}$
与直接计算结果一致,验证成立。
- 频域相乘:
-
线性卷积所需最小DFT长度
- 线性卷积输出长度为 $N + M - 1 = 4 + 4 - 1 = 7$。
- 若DFT长度小于7,循环卷积会产生混叠,无法等效于线性卷积。
- 因此,DFT长度至少需为7点,才能保证循环卷积等于线性卷积。
解析
本题主要考察有限长序列的循环卷积、离散傅里叶变换(DFT)、频域相乘与时域循环卷积的性质以及线性卷积所需的最小DFT长度等知识。解题思路如下:
- 计算4点循环卷积:根据循环卷积的定义公式$y_c(n) = \sum_{m=0}^{N - 1} x_1(m) \cdot x_2((n - m) \mod N)$(其中$N = 4$),对$n$从$0$到$3$分别进行计算。
- 当$n = 0$时:
$\begin{align*} y_c(0) &= \sum_{m=0}^{3} x_1(m) \cdot x_2((0 - m) \mod 4)\\ &= x_1(0)x_2(0)+x_1(1)x_2(3)+x_1(2)x_2(2)+x_1(3)x_2(1)\\ &= 1\times2 + 2\times(-1)+0\times1 + 1\times1\\ &= 2 - 2 + 0 + 1\\ &= 1 \end{align*}$ - 当$n = 1$时:
$\begin{align*} y_c(1) &= \sum_{m=0}^{3} x_1(m) \cdot x_2((1 - m) \mod 4)\\ &= x_1(0)x_2(1)+x_1(1)x_2(0)+x_1(2)x_2(3)+x_1(3)x_2(2)\\ &= 1\times1 + 2\times2+0\times(-1) + 1\times1\\ &= 1 + 4 + 0 + 1\\ &= 6 \end{align*}$ - 当$n = 2$时:
$\begin{align*} y_c(2) &= \sum_{m=0}^{3} x_1(m) \cdot x_2((2 - m) \mod 4)\\ &= x_1(0)x_2(2)+x_1(1)x_2(1)+x_1(2)x_2(0)+x_1(3)x_2(3)\\ &= 1\times1 + 2\times1+0\times2 + 1\times(-1)\\ &= 1 + 2 + 0 - 1\\ &= 2 \end{align*}$ - 当$n = 3$时:
$\begin{align*} y_c(3) &= \sum_{m=0}^{3} x_1(m) \cdot x_2((3 - m) \mod 4)\\ &= x_1(0)x_2(3)+x_1(1)x_2(2)+x_1(2)x_2(1)+x_1(3)x_2(0)\\ &= 1\times(-1) + 2\times1+0\times1 + 1\times2\\ &= -1 + 2 + 0 + 2\\ &= 3 \end{align*}$
所以$y_c(n) = \{1, 6, 2, 3\}$。
- 当$n = 0$时:
- 计算4点DFT:根据DFT的定义公式$X(k)=\sum_{n = 0}^{N - 1}x(n)W_N^{kn}$(其中$W_N = e^{-j\frac{2\pi}{N}}$,$N = 4$)分别计算$X_1(k)$和$X_2(k)$。
- 对于$x_1(n) = \{1, 2, 0, 1\}$:
- 当$k = 0$时:
$\begin{align*} X_1(0) &= \sum_{n = 0}^{3}x_1(n)W_4^{0\times n}\\ &= x_1(0)+x_1(1)+x_1(2)+x_1(3)\\ &= 1 + 2 + 0 + 1\\ &= 4 \end{align*}$ - 当$k = 1$时:
$\begin{align*} X_1(1) &= \sum_{n = 0}^{3}x_1(n)W_4^{1\times n}\\ &= x_1(0)+x_1(1)W_4^1+x_1(2)W_4^2+x_1(3)W_4^3\\ &= 1 + 2\times e^{-j\frac{2\pi}{4}}+0\times e^{-j\frac{4\pi}{4}} + 1\times e^{-j\frac{6\pi}{4}}\\ &= 1 + 2\times(-j)+0 + 1\times j\\ &= 1 - j \end{align*}$ - 当$k = 2$时:
$\begin{align*} X_1(2) &= \sum_{n = 0}^{3}x_1(n)W_4^{2\times n}\\ &= x_1(0)+x_1(1)W_4^2+x_1(2)W_4^4+x_1(3)W_4^6\\ &= 1 + 2\times e^{-j\frac{4\pi}{4}}+0\times e^{-j\frac{8\pi}{4}} + 1\times e^{-j\frac{12\pi}{4}}\\ &= 1 + 2\times(-1)+0 + 1\times(-1)\\ &= -2 \end{align*}$ - 当$k = 3$时:
$\begin{align*} X_1(3) &= \sum_{n = 0}^{3}x_1(n)W_4^{3\times n}\\ &= x_1(0)+x_1(1)W_4^3+x_1(2)W_4^6+x_1(3)W_4^9\\ &= 1 + 2\times e^{-j\frac{6\pi}{4}}+0\times e^{-j\frac{12\pi}{4}} + 1\times e^{-j\frac{18\pi}{4}}\\ &= 1 + 2\times j+0 + 1\times(-j)\\ &= 1 + j \end{align*}$
所以$X_1(k) = \{4, 1 - j, -2, 1 + j\}$。
- 当$k = 0$时:
- 对于$x_2(n) = \{2, 1, 1, -1\}$:
- 当$k = 0$时:
$\begin{align*} X_2(0) &= \sum_{n = 0}^{3}x_2(n)W_4^{0\times n}\\ &= x_2(0)+x_2(1)+x_2(2)+x_2(3)\\ &= 2 + 1 + 1 - 1\\ &= 3 \end{align*}$ - 当$k = 1$时:
$\begin{align*} X_2(1) &= \sum_{n = 0}^{3}x_2(n)W_4^{1\times n}\\ &= x_2(0)+x_2(1)W_4^1+x_2(2)W_4^2+x_2(3)W_4^3\\ &= 2 + 1\times e^{-j\frac{2\pi}{4}}+1\times e^{-j\frac{4\pi}{4}} - 1\times e^{-j\frac{6\pi}{4}}\\ &= 2 + 1\times(-j)+1\times(-1) - 1\times j\\ &= 1 - 2j \end{align*}$ - 当$k = 2$时:
$\begin{align*} X_2(2) &= \sum_{n = 0}^{3}x_2(n)W_4^{2\times n}\\ &= x_2(0)+x_2(1)W_4^2+x_2(2)W_4^4+x_2(3)W_4^6\\ &= 2 + 1\times e^{-j\frac{4\pi}{4}}+1\times e^{-j\frac{8\pi}{4}} - 1\times e^{-j\frac{12\pi}{4}}\\ &= 2 + 1\times(-1)+1\times1 - 1\times1\\ &= 3 \end{align*}$ - 当$k = 3$时:
$\begin{align*} X_2(3) &= \sum_{n = 0}^{3}x_2(n)W_4^{3\times n}\\ &= x_2(0)+x_2(1)W_4^3+x_2(2)W_4^6+x_2(3)W_4^9\\ &= 2 + 1\times e^{-j\frac{6\pi}{4}}+1\times e^{-j\frac{12\pi}{4}} - 1\times e^{-j\frac{18\pi}{4}}\\ &= 2 + 1\times j+1\times1 - 1\times(-j)\\ &= 1 + 2j \end{align*}$
所以$X_2(k) = \{3, 1 - 2j, 3, 1 + 2j\}$。
- 当$k = 0$时:
- 对于$x_1(n) = \{1, 2, 0, 1\}$:
- 频域相乘后IDFT验证:先计算$Y_c(k) = X_1(k)X_2(k)$,再对$Y_c(k)$做4点IDFT。
- 计算$Y_c(k)$:
- $Y_c(0)=X_1(0)X_2(0)=4\times3 = 12$
- $Y_c(1)=X_1(1)X_2(1)=(1 - j)(1 - 2j)=1 - 2j - j + 2j^2=1 - 3j - 2=-1 - 3j$
- $Y_c(2)=X_1(2)X_2(2)=(-2)\times3=-6$
- $Y_c(3)=X_1(3)X_2(3)=(1 + j)(1 + 2j)=1 + 2j + j + 2j^2=1 + 3j - 2=-1 + 3j$
所以$Y_c(k) = \{12, -1 - 3j, -6, -1 + 3j\}$。
- 对$Y_c(k)$做4点IDFT,根据IDFT公式$y(n)=\frac{1}{N}\sum_{k = 0}^{N - 1}Y(k)W_N^{-kn}$($N = 4$):
- 当$n = 0$时:
$\begin{align*} y_c(0) &=\frac{1}{4}\sum_{k = 0}^{3}Y_c(k)W_4^{-0\times k}\\ &=\frac{1}{4}(Y_c(0)+Y_c(1)+Y_c(2)+Y_c(3))\\ &=\frac{1}{4}(12-1 - 3j-6-1 + 3j)\\ &=\frac{1}{4}\times4\\ &= 1 \end{align*}$ - 当$n = 1$时:
$\begin{align*} y_c(1) &=\frac{1}{4}\sum_{k = 0}^{3}Y_c(k)W_4^{-1\times k}\\ &=\frac{1}{4}(Y_c(0)+Y_c(1)W_4^{-1}+Y_c(2)W_4^{-2}+Y_c(3)W_4^{-3})\\ &=\frac{1}{4}(12+(-1 - 3j)j+(-6)(-1)+(-1 + 3j)(-j))\\ &=\frac{1}{4}(12 - j - 3j^2+6 + j - 3j^2)\\ &=\frac{1}{4}(12 + 3+6 + 3)\\ &= 6 \end{align*}$ - 当$n = 2$时:
$\begin{align*} y_c(2) &=\frac{1}{4}\sum_{k = 0}^{3}Y_c(k)W_4^{-2\times k}\\ &=\frac{1}{4}(Y_c(0)+Y_c(1)W_4^{-2}+Y_c(2)W_4^{-4}+Y_c(3)W_4^{-6})\\ &=\frac{1}{4}(12+(-1 - 3j)(-1)+(-6)(1)+(-1 + 3j)(-1))\\ &=\frac{1}{4}(12 + 1 + 3j-6 + 1 - 3j)\\ &=\frac{1}{4}\times8\\ &= 2 \end{align*}$ - 当$n = 3$时:
$\begin{align*} y_c(3) &=\frac{1}{4}\sum_{k = 0}^{3}Y_c(k)W_4^{-3\times k}\\ &=\frac{1}{4}(Y_c(0)+Y_c(1)W_4^{-3}+Y_c(2)W_4^{-6}+Y_c(3)W_4^{-9})\\ &=\frac{1}{4}(12+(-1 - 3j)(-j)+(-6)(-1)+(-1 + 3j)j)\\ &=\frac{1}{4}(12 + j + 3j^2+6 - j + 3j^2)\\ &=\frac{1}{4}(12 - 3+6 - 3)\\ &= 3 \end{align*}$
所以$y_c(n) = \{1, 6, 2, 3\}$,与直接计算结果一致。
- 当$n = 0$时:
- 计算$Y_c(k)$:
- 计算线性卷积所需最小DFT长度:线性卷积输出长度为$N + M - 1$(其中$N$和$M$分别是两个序列的长度),这里$N = M = 4$,所以输出长度为$4 + 4 - 1 = 7$。若DFT长度小于7,循环卷积会产生混叠,无法等效于线性卷积。因此,DFT长度至少需为7点,才能保证循环卷积等于线性卷积。