为什么二维函数的matlab fft根据尺寸顺序给出不同的答案？

Question

Matlab的fft2是fftn的包装，但是一维fft沿尺寸的操作顺序不同。为了追赶，下面给出了一个非零的答案。为什么是这样？我找不到在线错误报告。

foo=rand(64);
err=sum(sum(fft(fft(foo,[],2),[],1)-fft(fft(foo,[],1),[],2)))