ts = 1000;%采樣率
t = 0:1/ts:1;
x = sin(2*pi*2*t)+2;%原始信號，頻率2Hz
x2 = sin(2*pi*50*t);%載波，頻率50Hz
y = x.*x2;%調制
plot(y);%調制信號
y = abs(hilbert(y));%hilbert變換，求絕對值得到瞬時幅度
hold on;
plot(y,'r')%包洛
plot(x+0.1,'g')%原始信號加0.1，與檢測到的包洛比較。