使用直接法求取瞬时频率的步骤

---

1 有限差分求导公式

---

2 求取瞬时频率的【直接法步骤】

已知数据是t与x_t，二者均为离散序列；y_t是对x_t进行希尔伯特变换后的结果，即y_t=hilbert(x_t)。求取瞬时频率的具体步骤如下：

1、根据下述公式求解瞬时相位：

2、对相位展开得到单调递增函数（需要对相位展开，是因为瞬时频率应该是正的且是连续的，因此相位应该是单调递增的）。具体步骤如下两张图：


3、使用步骤2中展开后的相位对时间求导，得到瞬时频率，其中使用的公式如下：

其中φ(t)是相位，ω(t)是频率，其单位为Hz。具体来说，上式中的一阶求导公式使用的是【1 有限差分求导公式】一节的第一个公式。

4、对上一步得到的瞬时频率做平滑处理。为了去除计算中出现的畸点，可以对瞬时频率做平滑或滤波。本文采用的是滑动平均法，即每一个瞬时频率的计算结果都是 10 个采样点的平均值。

END

---

3 参考文献

1、有限差分法求导

2、任达千,杨世锡,吴昭同,严拱标,孟庆波.信号瞬时频率直接计算法与Hilbert变换及Teager能量法比较[J].机械工程学报,2013,49(09):42-48.