1. 问题描述
信号处理中，对于较为负复杂的信号，我们有时候需要提取包络线分析。如下图：
mark
2. 技术背景
提取包络线，在信号处理中往往采用Hilbert变换。
3. 解决方案
h = hilbert(V);
4. 实施示例
4.1 波形参数设置
%% 基本参数
A = 1; % 激励信号幅值
N = 15; % cycle数，即激励信号波峰数
fck = 100; % 激励中心频率 kHz
fc = 100e3; % 激励信号中心频率，Hz
T0 = 1e-4; % 导波传播时间
T = 4.0*T0; % 激励持续时长
dt = 1/(20*fc)/2; % 时间步长，在最大步长基础上除以2
t = [0:dt:T]'; % 时间序列
4.2 波形构造与Hilbert变换
V = A*[heaviside(t)-heaviside(t-N/fc)].*... % 时域输入信号求解
(1-cos(2*pi*fc*t/N)).*sin(2*pi*fc*t);
h = hilbert(V); % Hilbert变换
mark
4.3 包络绘制
%% 包络绘制
figure
hold on
plot(t,V)
plot(t,abs(h),'-.')
grid on
此处注意，对变换后的变量h取范数（绝对值）。
mark
5. 常见问题
对于较复杂的波形，包络效果可能不理想。
mark
本文用时 30 m
