% 计算包络信号 x_envelope = abs(hilbert(x)); % 绘制包络信号 t = (0:length(x)-1)/fs; % 时间轴 plot(t, x_envelope); % 计算包络谱 N = length(x_envelope); f = (0:N-1)*(fs/N); % 频率轴 envelope_spectrum = abs(fft(x_envelope)).^2 / N; % 绘制包络谱 figure; plot(f, envelope_spectrum); xlabel('Frequency (Hz)'); ylabel('Power');

这段代码计算了振动信号的包络信号，并绘制了包络信号和包络谱。

首先，使用 `hilbert` 函数计算了原始振动信号 `x` 的包络信号 `x_envelope`，并使用 `abs` 函数获取其幅值。

然后，通过计算时间轴 `t`，使用 `plot` 函数绘制了包络信号的时域波形。

接下来，计算了包络信号的长度 `N` 和频率轴 `f`，使用 `fft` 函数计算了包络信号的傅里叶变换，并通过除以 `N` 和取绝对值的平方得到包络谱 `envelope_spectrum`。

最后，创建了一个新的图形窗口，并使用 `plot` 函数绘制了包络谱，横轴为频率（单位为 Hz），纵轴为功率。

需要注意的是，这段代码中缺少了一行 `subplot(1, 2, 1)`，它用于创建一个1行2列的子图中的第1个子图，并且在这个子图中绘制包络信号。如果想要在同一个图形窗口中同时显示包络信号和包络谱，可以将这行代码添加到绘制包络信号之前。