例题
计算以下图像的离散傅里叶变换
f = [ 1 4 4 1 2 4 4 2 2 4 4 2 1 4 4 1 ] f=\begin{bmatrix}1&4&4&1\\2&4&4&2\\2&4&4&2\\1&4&4&1\end{bmatrix} f=⎣⎢⎢⎡1221444444441221⎦⎥⎥⎤
计算方式当然不是通过原始公式硬算,而是构造权值矩阵来计算。
考虑二维DFT的分解:
F ( u , v ) = ∑ x = 0 M − 1 ∑ y = 0 N − 1 f ( x , y ) e − j 2 π ( u x / M + v y / N ) = ∑ x = 0 M − 1 e − j 2 π u x / M ∑ y = 0 N − 1 f ( x , y ) e − j 2 π v y / N = ∑ x = 0 M − 1 e − j 2 π u x / M F ( x , v ) F(u,v)=\sum^{M-1}_{x=0}\sum^{N-1}_{y=0}f(x,y)e^{-j2\pi (ux/M+vy/N)}\\=\sum^{M-1}_{x=0}e^{-j2\pi ux/M}\sum^{N-1}_{y=0}f(x,y)e^{-j2\pi vy/N}\\=\sum^{M-1}_{x=0}e^{-j2\pi ux/M}F(x,v) F(u,v)=x=0∑M−1y=0∑N−1f(x,y)e−j2π(ux/M+vy/N)=x=0∑M−1e−j2πux/My=0∑N−1f(x,y)e−j2πvy/N=x=0∑M−1e−j2πux/MF(x,v)
注意到,对每一个 x x x,都有: