Matlab 振动信号处理-单边频谱

clc
clear
close all

%% 主要参数

Fs=100;%采样频率
T=1/Fs;%采样时间间隔
L=4000;%选取进行FFT分析的信号长度(点数)

%% 生成进行FFT分析的信号

t=(0:L-1)/Fs;%时间向量
X=sin(2*pi*13*t)+2*sin(2*pi*27*t)+randn(size(t));%正弦信号＋零均值随机噪声

%% 单边频谱

Y=fft(X);
P1=abs(Y)/L;%双边频谱P1
P2=P1(1:L/2+1);
P2(2:end-1)=2*P2(2:end-1);%单边频谱P2
f=Fs*(0:(L/2))/L;%单边频率向量

%% 绘图

figure(1)
subplot(2,1,1)
plot(t,X)
xlabel('t [s]')
ylabel('X(t) [g]')
title('时域信号')
subplot(2,1,2)
plot(f,P2)
xlabel('f [Hz]')
ylabel('|P2(f)| [g]')
title('单边频谱')

Result