MATLAB编写的fft2程序更新版

根据快速傅立叶变换的信号流图可知，可将整个过程中所有的数据组成一个二维数组data(N,M+1)，数组共有N行，M+1列（傅立叶变换分为M=log2(N)级，再加上第一级倒序数组输入，则共有M+1列）。除第一列单独赋值外，其余列则按照共同的规律来赋值。这里则详细说明其的规律性。并给出源程序！