使用matlab对振动信号进行滤波后绘制包络谱

处理振动信号时，常常需要对信号进行滤波以去除噪声等干扰。本文将介绍如何使用 MATLAB 对振动信号进行滤波，并绘制其包络谱。

首先，假设我们有一个振动信号 `x`，其采样频率为 `Fs`。我们可以使用 MATLAB 中的 `filter` 函数对其进行滤波，例如：

[b,a] = butter(4, 100/(Fs/2), 'low'); % 低通滤波器
y = filter(b,a,x); % 滤波后的信号

上述代码中，我们选择了一个 4 阶巴特沃斯低通滤波器，截止频率为 100 Hz。使用 `filter` 函数对信号进行滤波后，得到了滤波后的信号 `y`。

接下来，我们可以使用 MATLAB 中的 `hilbert` 函数对信号进行解析，得到其包络信号 `z`，例如：

z = abs(hilbert(y)); % 包络信号

最后，我们可以使用 MATLAB 中的 `pwelch` 函数对包络信号进行功率谱密度估计，并绘制出其包络谱，例如：

[P,f] = pwelch(z,[],[],[],Fs,'onesided'); % 包络信号的功率谱密度估计
loglog(f,P); % 绘制包络谱
xlabel('Frequency (Hz)'); ylabel('Power/frequency (dB/Hz)');

上述代码中，我们选择了对数坐标系来绘制包络谱。可以看到，使用 MATLAB 对振动信号进行滤波并绘制出其包络谱非常简单。