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('单边频谱')