1、准备待处理数据，可将需要进行FFt变换的数据存储在TXT文件中待使用，也可直接在.m文件中列出。
2、在.m文件中编写：
clear;
clc;
close all;
s=[待处理的数据列表];%或是直接引用存储在TXT文件中的数据（a=textread('归一化-差1.txt','%s')';）;
figure
plot(s)
Fs=200;%采样频率的设置
f = fftshift(fft(s));
w = linspace(-Fs/2, Fs/2, length(s));%频率坐标，单位Hz
figure
plot(w,abs(f));
axis([0 Fs/2 min(abs(f)) max(abs(f))])