求助关于使用hilbert变换对调频、调相信号进行包络解调的MATLAB仿真问题

zmts

我现在要对一个调频信号进行包络解调，使用MATLAB进行仿真，信号的表达式在MATLAB中表示为：
fs=16384;%设定采样频率
fc=6000;
f=40;
N=16384; %采样点数
n=0:N-1;
t=(0:1/fs:(N-1)/fs);
x = cos(2*pi*fc*t+30*sin(2*pi*f*t));
x是调相信号。对其进行解调，步骤如下：
1、现用MATLAB中自带的函数hilbert(x),得到解析信号S；
2、再对S使用函数angle(imag(S),real(S)),得到瞬时的相位信号；
3、对相位求导，恢复出调制信号。
这个方法肯定有问题，我在运行时始终得不到正确的结果。请高手看看问题在哪？或者有满意其他的方法可以恢复出调制信号的。