简单的FFT分析代码

数据保存在txt文件里，以列的型式存在，MATLAB代码如下：

clc;

X = load('38-3.txt');  %打开数据存放的文件
Fs=4000;                 %采样频率，可更改为所采样的频率如3000.......
Tload=12.5;             %用X文件里的数据列数/4000得到
Len=Tload*Fs;
Va=X(1:Len,2);        %只需更改2这一值，代表所需分析的数据在第二列
T = 1/Fs;                  % Sampling period
Y = fft(Va); 
L =length(Y);            % Length of signal 
t = (0:L-1)*T;             % Time vector

P2 = abs(Y/L);
P1 = P2(1:L/2+1);
P1(2:end-1) = 2*P1(2:end-1);

f = Fs*(0:(L/2))/L;
plot(f,P1)

xlabel('f (Hz)')
ylabel('|P1(f)|')