利用Excel进行FFT和Fourier分析的基本步骤.pdfVIP

利用 Excel 进行 FFT 和 Fourier 分析的基本步骤 实例：杭州市 2000 人口分布密度 [根据 2000 年人口普查的街道数据经环带(rings)平均 计算得到的结果，数据由冯健博士处理] 。下面的变换实质是一种空间自相关的分析过程。 第一步，录入数据 在 Excel 中录入数据不赘述（见表 1）。 表 1 原始数据序列 表2 补充后的数据序列 第二步，补充数据 由于Fourier 变换（FT ）一般是借助快速Fourier 变换（Fast Fourier Transformation, FFT ） 算法，而这种算法的技术过程涉及到对称处理，故数据序列的长度必须是2N （N=1,2,3, …，）。 如果数据序列长度不是 2N 4 ，就必须对数据进行补充或者裁减。现在数据长度是 26 ，介于 2 =16 到 25=32 之间，而 26 到 32 更近一些，如果裁减数据，就会损失许多信息。因此，采用补充 数据的方式。 补充的方法非常简单，在数据序列后面加 0，直到序列长度为 32=25 为止（表 2 ）。当然， 延续到 64=26 也可以，总之必须是 2 的整数倍。不过，补充的“虚拟数据”越多，变换结果 的误差也就越大。 1 第三步，Fourier 变换的选项设置 沿着工具 （Tools ）→数据分析 （Data Analysis ）的路径打开数据分析复选框（图 1）。 图 1 数据分析（Data Analysis）的路径 在数据分析选项框中选择傅立叶分析 （Fourier Analysis ）（图2 ）。 图 2 数据分析（Data Analysis） 在 Fourier 分析对话框中进行如下设置：在输入区域中输入数据序列的单元格范围 “$B$1:$B$33”；选中“标志位于第一行（L）”；将输出区域设为“$C$2”或者“$C$2:$C$33” （图 3a）。 a 2 b 图 3 傅立叶分析（Fourier Analysis） 注意：如果“输入区域”设为“$B$2:$B$33”，则不选“标志位于第一行（L）”（图3b）。 表 3 FFT 的结果 3 第四步，输出 FFT 结果 选项设置完毕以后，确定（OK ），立即得到FFT 结果（表 3 ）。 显然，表 3 给出的都是复数（complex numbers ）。假定一个数据序列表为f (t) ，则理论 上 Fourier 变换的结果为 ∞ F (ω) ∫−∞f (t)e −j ωt dt = F[f (t)], ( −∞ω∞) 表 3 中给出的正是相应于F(ω) 的复数，这里ω为角频率。 第五步，计算功率谱 Excel 好像不能自动计算功率谱，这需要我们利用有关函数进行计算。计算公式为 1 2 1 2 2