【信号处理】给定信号幅度谱的包络提取（Matlab实现）

“在代码的海洋里，有无尽的知识等待你去发现。我就是那艘领航的船，带你乘风破浪，驶向代码的彼岸。

 💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码实现

💥1 概述

 给定信号幅度谱的包络提取是信号处理中的一个重要任务。其目的是从给定信号的幅度谱中获取其包络信息。 在这个过程中，首先需要对信号进行频谱分析以获得幅度谱。然后，通过特定的算法和技术来提取幅度谱的包络。 包络提取有助于了解信号的整体变化趋势和主要特征，能够去除一些细节和噪声的影响，突出信号的主要轮廓。 这对于许多应用具有重要意义，例如在通信系统中用于信号的检测和识别、音频处理中对声音特征的分析、故障诊断中对异常信号的检测等。 提取包络的方法多种多样，包括希尔伯特变换、峰值检测、平滑处理等，选择合适的方法取决于信号的特点和具体的应用需求。

📚2 运行结果

主函数部分代码：

clear, clc, close all

% load an audio file
[x, fs] = audioread('track.wav', [1 4096]);   % load an audio file
x = x(:, 1);                                  % get the first channel

% perform spectral analysis
N = length(x);   
w = hanning(N, 'periodic');
[Xamp, f] = periodogram(x, w, N, fs, 'power');
Xamp = 20*log10(sqrt(Xamp)*sqrt(2));

% spectral envelope extraction
Xenv = specenv(Xamp, f);

% plot the spectrum
plot(f, Xamp)
grid on
hold on
plot(f, Xenv, 'r', 'LineWidth', 1.5)
xlim([0 max(f)])
set(gca, 'FontName', 'Times New Roman', 'FontSize', 14)
xlabel('Frequency, Hz')
ylabel('Magnitude, dB')
title('Amplitude spectrum of the signal and its envelope')
legend('Original spectrum', 'Spectral envelope')

🎉3 参考文献

文章中一些内容引自网络，会注明出处或引用为参考文献，难免有未尽之处，如有不妥，请随时联系删除。

[1]陈炫昂,葛玖浩,胡宝旺.交流电磁场钢轨表面裂纹高速检测信号处理模块开发[J/OL].测控技术:1-8[2024-06-21].https://doi.org/10.19708/j.ckjs.2023.05.224.

[2]程锦房,谢昌奇,张伽伟,等.相干累加与AR滤波相结合的舰船轴频电场信号处理方法[J].国防科技大学学报,2024,46(03):229-236.

🌈4 Matlab代码实现