Matlab中fft2（二维）与多次一维FFT计算结果位置互换问题

首先，咱们得先明确Matlab里fft函数的维度规则：默认情况下fft(A)是对矩阵的列向量（维度1）执行一维FFT。如果要对行向量（维度2）做FFT，必须显式指定第三个参数dim=2，也就是fft(A, [], 2)。

你提到的“位置多处互换”问题，大概率是因为你在对行做FFT的时候，没有指定维度参数，而是用了转置这类绕弯的操作（比如先转置矩阵再做fft，然后再转置回来），这就导致最终的频谱矩阵是原fft2结果的转置，自然会出现数值相同但位置互换的情况。

正确的先行后列2D FFT实现方式 #

咱们用代码来对比：

• 直接调用fft2：

  X = randn(4,4); % 随机测试矩阵
  output2 = fft2(X);
  

• 先行后列的一维FFT组合（正确写法）：

  % 第一步：对每一行做一维FFT（指定dim=2）
  row_fft = fft(X, [], 2);
  % 第二步：对row_fft的每一列做一维FFT（指定dim=1，也可以省略，因为默认是列）
  output1 = fft(row_fft);
  

  这时候你再对比output1和output2，会发现它们的数值完全一致，位置也不会有任何互换。

为什么顺序不影响结果？ #

从数学原理来说，2D FFT是可分离的，先对所有行做1D FFT再对列做，和先对所有列做再对行做，最终的结果是完全相同的——这也是fft2内部的实现逻辑（Matlab的fft2本质就是先列后行的一维FFT组合，但因为可分离性，先行后列的结果和它完全等价）。

你之前操作出错的可能场景 #

如果你的代码是类似这样的：

% 错误写法：用转置来实现行FFT
row_fft_wrong = fft(X').';
output1_wrong = fft(row_fft_wrong);

那output1_wrong其实等价于fft2(X')'，也就是原fft2(X)的转置矩阵，这就会出现你看到的“数值相同但位置互换”的现象——因为频谱的行和列被颠倒了。

总结一下：只要你在对行做FFT时正确指定dim=2，先行后列的结果就会和fft2完全一致，不会有位置互换的问题。

内容的提问来源于stack exchange，提问作者four_lines