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

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

实例：杭州市2000人口布密度[根据2000年人口普查的街道数据经环带(rings)平均

计算得到的结果，数据由冯健博士处理]。下面的变换实质是一种空间自相关的析过程。

第一步，录入数据

在Excel中录入数据不赘述（见表1）。

表1原始数据序列表2补充后的数据序列

第二步，补充数据

由于Fourier变换（FT）一般是借助快速Fourier变换（FastFourierTransformation,FFT）

NN

算法，而这种算法的技术过程涉及到对称处理，故数据序列的长度必须是2（=1,2,3,…，）。

N264

如果数据序列长度不是2，就必须对数据进行补充或者裁减。现在数据长度是，介于2=16

到25=32之间，而26到32更近一些，如果裁减数据，就会损失许多信息。因此，采用补充

数据的方式。

补充的方法非常简单，在数据序列后面加0，直到序列长度为32=25为止（表2）。当然，

延续到64=26也可以，总之必须是2的整数倍。不过，补充的“虚拟数据”越多，变换结果

的误差也就越大。

1

第三步，Fourier变换的选项设置

沿着工具（Tools）→数据析（DataAnalysis）的路径打开数据析复选框（图1）。

图1数据分析（DataAnalysis）的路径

在数据析选项框中选择傅立叶析（FourierAnalysis）（图2）。

图2数据分析（DataAnalysis）

在Fourier析对话框中进行如下设置：在输入区域中输入数据序列的单元格范围

“$B$1:$B$33”；选中“标志位于第一行（L）”；将输出区域设为“$C$2”或者“$C$2:$C$33”

（图3a）。

a

2

b

图3傅立叶分析（FourierAnalysis）

注意：如果“输入区域”设为“$B$2:$B$33”，则不选“标志位于第一行（L）”（图3b）。

表3FFT的结果

3

第四步，输出FFT结果

选项设置完毕以后，确定（OK），立即得到FFT结果（表3）。

ft

显然，表3给出的都是复数（complexnumbers）。假定一个数据序列表为()，则理论

上Fourier变换的结果为



FftejtdtFft

()()

=[()],()

Fωω

表3中给出的正是相应于()的复数，这里为角频率。

第五步，计算功率谱

Excel好像不能自动计算功率谱，这需要我们利用有关