用希尔伯特黄变换HHT求时频谱和边际谱.docxVIP

用希尔伯特黄变换（HHT ）求时频谱和边际谱 1.什么是HHT? HHT就是先将信号进行经验模态分解（EMD分解），然后将分解后的每个IMF分量进行 Hilbert变换，得到信号的时频属性的一种时频分析方法。 2.EMD分解的步骤。 EMD分解 Cl)找出信号中所有局部极大值井用三次样条函数连接成上包络； Cl) 同理，利用三次样条插值函数连接所有局部极小值构成下包络； C2)如果咿足IMF C2) 如果咿足IMF的条件，那么图就是求得的第一个IMF分量; C3） 否则棚作为原始信号进行Cl） ~ C2）的步骤，直到第k次迭代后差值 成为一个IMF,记为： q (t)=E) Matlab论坛 iLoveMatlab.cn EMD分解的流程图如下: Matlab论坛 iLoveMatlab-cii 3.实例演示。 给定频率分别为10Hz和35Hz的两个正弦信号相叠加的复合信号，采样频率fs=2048Hz的 信号，表达式如下：y=5sin(2*pi*10t)+5*sin(2*pi*35t) 为了对比，先用fft对求上述信号的幅频和相频曲线 代码： function fftfenxi clear;clc; N=2048; %fft默认计算的信号是从0开始的 t=linspace(1,2,N);deta=t(2)-t(1);1/deta x=5*sin(2*pi*10*t)+5*sin(2