Continuous function – FT & Conv
有两个时域函数 g(t)g(t)g(t) 和 h(t)h(t)h(t), 则有时域卷积等价于频域相乘 (反之亦然, 仅多一个系数), 即
结论:FT(g∗h)(g\ast h)(g∗h)=FT(g)×(g)\times(g)×FT(h)(h)(h)
证明:两者卷积后的函数 p(t)p(t)p(t) 可写作
p(t)=g(t)∗h(t)=∫−∞∞g(τ)h(t−τ)dτp(t)=g(t)\ast h(t)=\int_{-\infty}^\infty g(\tau)h(t-\tau)d\taup(t)=g(t)∗h(t)=∫−∞∞g(τ)h(t−τ)dτ
进一步, p(t)p(t)p(t) 的傅里叶变换为
P(f)=∫−∞∞p(t)ej2πftdt=∫−∞∞∫−∞∞g(τ)h(t−τ)dτe−j2πftdt=∫−∞∞h(t−τ)e−j2πftdt∫−∞∞g(τ)dτ=∫−∞∞h(t′)e−j2πft′dt′e−j2πfτ∫−∞∞g(τ)dτ=H(f)∫−∞∞g(τ)e−j2πfτdτ=G(f)H(f)P(f)=\int_{-\infty}^\infty p(t)e^{j2\pi f t}dt=\int_{-\infty}^\infty \int_{-\infty}^\infty g(\tau)h(t-\tau)d\tau e^{-j2\pi f t}dt \\ =\int_{-\infty}^\infty h(t-\tau)e^{-j2\pi f t}dt \int_{-\infty}^\infty g(\tau)d\tau =\int_{-\infty}^\infty h(t')e^{-j2\pi f t'}dt' e^{-j2\pi f \tau}\int_{-\infty}^\infty g(\tau)d\tau\\ =H(f) \int_{-\infty}^\infty g(\tau)e^{-j2\pi f \tau} d\tau= G(f)H(f) P(f)=∫−∞∞p(t)ej2πftdt=∫−∞∞∫−∞∞g(τ)h(t−τ)dτe−j2πftdt=∫−∞∞h(t−τ)e−j2πftdt∫−∞∞g(τ)dτ=∫−∞∞h(t′)e−j2πft′dt′e−j2πfτ∫−∞∞g(τ)dτ=H(f)∫−∞∞g(τ)e−j2πfτdτ=G(f)H(f)
Discrete sequence – DFT & Conv
有两个定义在整数上的函数 g[n]g[n]g[n] 和 h[n]h[n]h[n], 它们的卷积可以写作
p[n]=∑i=0∞g[i]h[n−i]p[n]=\sum_{i=0}^{\infty}g[i]h[n-i]p[n]=i=0∑∞g[i]h[n−i]
但是当 g[n]g[n]g[n] 和 h[n]h[n]h[n] 有限长 NNN, MMM时, 他们的卷积可以写作
p[n]=∑i=0N−1g[i]h[n−i],n=0,1,2,...,N+M−1p[n]=\sum_{i=0}^{N-1}g[i]h[n-i], n=0,1,2,...,N+M-1p[n]=i=0∑N−1g[i]h[n−i],n=0,1,2,...,N+M−1
结论:DFT(g∗h)(g\ast h)(g∗h) = DFT(g)×(g)\times(g)×DFT(h)(h)(h) 成立当且仅当K点DFT变换的长度 K>N+M−1K>N+M-1K>N+M−1.
这个很显然,因为 g∗hg\ast hg∗h 卷积出来的有 N+M−1N+M-1N+M−1 点,DFT点数不能小于这个值。
OFDM symbol – Cyclic Conv
在OFDM信号处理中,g[n]g[n]g[n] 一般是对IDFT后一个OFDM符号的 NNN 点采样. h[n]h[n]h[n] 一般是多径信道 (其长度 MMM 一般会小于 NNN)。 在接收端,我们的 DFT 只会做 NNN 点的,这就需要 CP 把信道的卷积变成循环卷积。
设g[n]g[n]g[n] 和 h[n]h[n]h[n] 长度为 NNN, MMM, 且 M≤NM\leq NM≤N. 他们的循环卷积可以写作
p[n]=g[n]⊗h[n]=∑i=1−NN−1g[i]h[n−i],n=0,1,2,...,N−1p[n]=g[n] \otimes h[n]=\sum_{i=1-N}^{N-1}g[i]h[n-i], n=0,1,2,...,N-1p[n]=g[n]⊗h[n]=i=1−N∑N−1g[i]h[n−i],n=0,1,2,...,N−1
公式里面的 NNN 其实都应该是 max(N,M)\max(N,M)max(N,M),若 M>NM> NM>N 就把 g[n]g[n]g[n] 和 h[n]h[n]h[n] 倒反就行了。
循环卷积等价于把其中长的那个序列(图中是 g[n]g[n]g[n])无限循环延拓,在用另一个短序列滑动相关。这样得出的序列 p[n]p[n]p[n] 是一个周期序列,且周期为长序列的长度 (图中是NNN)。每个p[n]p[n]p[n]都是 min(N,M)\min(N,M)min(N,M) 个点相关的结果