DFT和FFT功率（波数）谱分析
参考帖子https://www.mathworks.com/matlabcentral/answers/24965-fft2-function-in-matlab一维傅里叶变换Y(omega)=【f(t)exp(-iomegat)dt】【】表示负无穷到正无穷积分对时间变换得到的是频率谱，omega=2pif，为圆频率，f为频率，单位s-1，即赫兹，单位时间内的cycle数对空间变换得到的是波数谱，omega=2pik，k为波数，单位为m-1，单位空间内的cycle数编程时积分改为累加，即离散傅
