简单的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)|')

上一篇:【2021/09/24】每日总结


下一篇:[学习笔记]半平面交