数据保存在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)|')