matlab 求一維信號的上下包絡線

clear all; clc; close all;
fs=400;                                 % 採樣頻率
N=400;                                  % 數據長度
n=0:1:N-1;
dt=1/fs;
t=n*dt;                                 % 時間序列
A=0.5;                                  % 相位調製幅值
x=(1+0.5*cos(2*pi*5*t)).*cos(2*pi*50*t+A*sin(2*pi*10*t));  % 信號序列
z=hilbert(x');                          % 希爾伯特變換
a=abs(z);                               % 包絡線
figure;
plot(x);hold on ;plot(a); hold on; plot((-1)*a,'b')