项目说明

使用 Matlab 绘制一段信号的频谱与包络谱。

代码实现

clear;
clc;

%% 原始信号
dt = 1/12000;  % 两个点之间的时间采样间隔
t = 0:dt:10;
data = sin(2*pi*t);

figure
subplot(311);
plot(t, data);
title('振动信号')

%% 频域分析
Fs = 12000;  % 采样频率必须是dt的倒数
Hz = 20*18;
N = length(data);  % 信号重采样的长度
f = (0:N-1)*Fs/N;  % 频率
data_f = abs(fft(data-mean(data)))/(N/2);  % 时频转换

subplot(312);
plot(f(1:Hz),data_f(1:Hz));
title('频域信号')

%% 包络谱分析
HV = abs(hilbert(data));
HV_m = HV - mean(HV);
HV_f = abs(fft(HV_m))/(N/2);

subplot(313);
plot(f(1:Hz), HV_f(1:Hz))
title('信号包络谱');