对振动信号进行包洛谱分析提取特征值Matlab代码

以下是用Matlab实现包络谱分析的示例代码：

% 读取信号
load('vibration_signal.mat');
signal = vibration_signal;

% 设置参数
fs = 1000;   % 采样频率
nfft = 1024; % fft点数
window_length = 0.05; % 窗口长度，单位秒
overlap_length = 0.025; % 窗口重叠长度，单位秒

% 计算包络谱
t = (0:length(signal)-1)/fs;
wlen = round(fs*window_length);
noverlap = round(fs*overlap_length);
[S,F,T,P] = spectrogram(signal,hann(wlen),noverlap,nfft,fs,'centered');
env = abs(hilbert(P));
envelope_spectrum = mean(env,2);

% 绘制包络谱
figure;
plot(F,envelope_spectrum);
xlabel('Frequency (Hz)');
ylabel('Amplitude');
title('Envelop Spectrum');

说明：

- `vibration_signal.mat` 是一个包含振动信号的Matlab数据文件，可以替换成自己的信号。
- `fs` 是采样频率，单位Hz。
- `nfft` 是fft点数，通常选择2的次幂。
- `window_length` 是窗口长度，单位秒。
- `overlap_length` 是窗口重叠长度，单位秒。
- `spectrogram` 函数可以将信号分帧，并计算每一帧的短时傅里叶变换。
- `hilbert` 函数可以计算信号的解析信号，即包络信号。
- `abs` 函数可以取包络信号的振幅。
- `mean` 函数可以计算所有帧的包络信号的平均值，得到包络谱。