matlab求kcf算法响应图_HHT算法的若干问题及改进方法—采样周期对 HHT 算法的影响...

4.5采样周期对 HHT 算法的影响

4.5.1 整周期采样的定义

在仿真实验中我们发现，对于同一个待测信号，当采样频率与采样点数选取不同的值时，得到的HHT频谱图会有所差异。经过大量的对比分析，作者发现当整周期采样时得到的HHT频谱图效果要比非整周期采样时得到的频谱图效果好。

整周期采样最初是在傅立叶变换中根据离散信号的要求提出来的，在此我们

把这个概念引入到希尔伯特-黄变换中来。设HHT频谱分辨率为Δf，则：

其中N是采样点数，fs是采样频率。在此我们给出一个整周期采样的简单判断式，即：

f 是待测信号频率，当 K 是整数时，我们就称其为整周期采样，否则为非整周期采样。

4.5.2 整周期采样对算法的影响

我们用以下两个实例来说明采样周期对HHT频谱的影响。

例1：对余弦信号x(t)求瞬时频率。

其中f=3MHz，t=[0,1000]s，采样频率fs=100MHz，采样点数N=1000时，代入到式(4-35)中，1000/100=K/3，则K=30是个整数，得到的HHT频谱图如图4.9所示。当采样频率fs=100MHz，采样点数N=1001时，代入到式(4-35)中，1001/100= K/3，则K=30.3是非整数，相应的频谱图如图4.10所示：

图 4.9、4.10 都是对同一个信号进行分析的，只是略微调整了一下采样点数，使前一个是整周期采样，后一个是非整周期采样，但是得到的瞬时频率效果图却相差很大，非整周期采样得到的频谱图波动非常大。

例 2：两个正弦叠加信号x(t)，采样频率fs=1000MHz，采样点数N =1000其中：

当f1=80MHz，f2=50MHz时，两个信号都是整周期采样，得到的HHT频谱图如图4.11所示：

当f1=80MHz， f2=50.5MHz时，第一个信号是整周期采样，而第二个信号是非整周期采样，得到的HHT频谱图如图4.12所示：

从图4.12我们可以发现当我们把待测信号中的f2由50MHz改为50.5MHz，其它的采样频率和采样点数保持不变的情况下，对于f1=80MHz的信号来说依然是整周期采样，而对f2=50.5MHz的信号来说却是非整周期采样，两者在HHT频谱图中也发生了相应的变化。为了看得更清楚，我们将图4.12中分解出的f2=50.5MHz的频谱图进行局部放大，如图4.13：

由以上两个例子我们可以看出非周期采样对HHT频谱图有较大的影响，尤其是在端点处。分析其原因主要是因为希尔伯特变换 MATLAB中的实现，实际上是把实信号进行傅立叶变换，然后把负频率部分强制赋零，然后进行短时傅立叶变换。从这个意义上说，两端差距较大，是因为傅立叶变换以后存在能量扩散，正频率部分存在从负频率扩散而来的能量，因此即使把负频率部分强制赋零，也不能完全保证消除所有的负频率能量。从滤波的角度来说，希尔伯特变换实际上是通过一个冲击响应为1/t的滤波器，在最前端由于认为全为零，所以滤波存在误差，在实验当中我们可以人为的修改一下参数，使其为整周期采样。在实际应用中我们无法确定频率到底是多少，但是在采样的时候可以做一些措施，比如用键相信号来保证起止采样点，多出现一些倍频或者分频，这样会改善实际信号的测量效果。

《来源科技文献，经本人分析整理，以技术会友，广交天下朋友》

知乎视频​www.zhihu.com