本文介绍了如何通过两次一维离散傅里叶变换（DFT）的叠加来实现二维DFT，从而降低计算成本。在MATLAB中，利用矩阵运算的方式对256x256图像的每一行和每一列进行DFT处理，构建了一个特定的矩阵（EMatrix），与图像矩阵相乘得到结果。这种方法在处理图像时效率较高，文中还给出了具体的MATLAB代码示例。
使用两次一维DFT的叠加方式实现二维的DFT运算可以有效的减少计算成本。
一维DFT公式
X [ k ] = ∑ n = 0 N − 1 x [ n ] e − j 2 π N k n X[k]=\sum_{n=0}^{N-1}x[n]e^{-j\frac{2\pi}{N}kn} X[k]=n=0∑N−1x[n]e−jN2πkn
将DFT转化为矩阵运算
在MATLAB中，使用矩阵运算比循环更加有效率。下面展示对每一行的处理：（以256x256大小图片举例）
对第i行进行DFT处理，则第i行第一个元素的卷积为：
x [ i , 1 ] e − j 2 π N 0 ∗ 0 + x [ i , 1 ] e − j 2 π N 0 ∗ 1 + . . . + x [ i , 1 ] e − j 2 π N 0 ∗ 255 x[i,1]e^{-j\frac{2\pi}{N}0*0}+x[i,1]e^{-j\frac{2\pi}{N}0*1}+...+x[i,1]e^{-j\frac{2\pi}{N}0*255} x[i,1]e−jN2π0∗0+x[i,1]e−jN2π0∗1+...+x[i,1]e−jN2π0∗255
所以将每行的每个元素都进行此运算之后得到：
X [ i ， 1 ] = x [ i , 1 ] e − j 2 π N 0 ∗ 0 + x [ i , 1 ] e − j 2 π N 0 ∗ 1 + . . . + x [ i , 1 ] e − j 2 π N 0 ∗ 255 X [ i ， 2 ] = x [ i , 2 ] e − j 2 π N 1 ∗ 0 + x [ i , 2 ] e − j 2 π N 1 ∗ 1 + . . . + x [ i , 2 ] e − j 2 π N 1 ∗ 2 . . . . . . X [ i ， 256 ] = x [ i , 256 ] e − j 2 π N 255
最低0.47元/天 解锁文章
727
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
