绘制信号包络谱的代码

要绘制信号的包络谱（也称为瞬时幅度谱），您可以使用Matlab中的`periodogram`或`spectrogram`函数，它们分别用于计算信号的一维周期ogram或多维谱图。这里是一个简单的例子，假设您有一个时间序列信号`x`：

% 假设 x 是您的时间序列数据
if nargin < 1
    t = linspace(0, 1, 1000); % 创建一个示例时间序列
    x = sin(2*pi*50*t) + sin(2*pi*120*t); % 正弦波叠加，模拟信号
else
    t = 1:length(x);
end

% 使用 periodogram 函数绘制一维包络谱
fs = 1/mean(diff(t)); % 假设采样频率为1/t步长
[Pxx, f] = periodogram(x, [], [], fs); % 计算并返回功率谱密度，f 是频率轴
env_spectrum = abs(Pxx); % 获取绝对值得到包络谱
figure;
plot(f, env_spectrum, 'LineWidth', 2); % 绘制包络谱
xlabel('Frequency (Hz)');
ylabel('Envelope Power');
title('Signal Envelope Spectrum');

% 或者使用 spectrogram 函数绘制二维谱图
window = hamming(length(x)); % 使用汉明窗进行平滑处理
H = spectrogram(x, window, [], [], fs); % 计算短时傅立叶变换
env_spectrogram = abs(H(:, :, 1)); % 取第一个频率通道即包络
imagesc(f, t, log(env_spectrogram)); % 对数显示，f 和 t 是频域和时域轴
colorbar;
xlabel('Frequency (Hz)');
ylabel('Time (samples)');
title('Signal Envelope Spectrogram');

这个例子首先计算了信号的功率谱密度（PSD），然后取其绝对值得到包络谱。如果希望看到信号随时间变化的频率响应，可以使用`spectrogram`函数。