希尔伯特变换(HilbertTransform)简介及其物理意义
希尔伯特变换(Hilbert Transform)简介及其物理意义
Hilbert变换简介
希尔伯特变换是信号处理中的⼀种常⽤⼿段,数学定义如下:
与卷积的概念进⾏对⽐,可以发现,上⾯的Hilbert变换的表达式实际上就是将原始信号和⼀个信号做卷积的结果。这个⽤来卷积的信号就是
1
h(t)=πt
因此,Hilbert变换可以看成是将原始信号通过⼀个滤波器,或者⼀个系统,这个系统的冲击响应为h(t)。
对h(t)做傅⾥叶变换,可以得到:
H(jω)=−jsgn(ω)
或者写成:
sgn()是符号函数。从频谱上来看,这个滤波器将我们的原始信号的正频率部分乘以-j,也就是说,保持幅度不变的条件下,将相位移动了-pi/2,⽽对于负频率成分,移动了pi/2。
下⾯这个⽰意图很直观地表⽰了Hilbert变换,在这⾥我画出了对原始信号做1到4次Hilbert变换的频谱⽰意图,是为了说明Hilbert变换的
⼏个性质:
⾸先,可以看到,两次希尔伯特变换后,原信号相位翻转了180°,所以,Hilbert逆变换的公式显⽽易见,就是将正变换加⼀个符号即可。另外,还可以看到,Hilbert变换四次后就变回本⾝了。还有其它的性质,⽐如:
如果⼀个信号是两个信号的卷积,即 y = conv(v,x) ,那么Hilbert(y) = conv(Hilbert(v),x) = conv(v,Hilbert(x))这个性质,只要意识到Hilbert变换本质上是卷积就可以明⽩。
x(t) 和 Hilbert(x(t))的能量以及平均功率相等,相关函数和功率谱相同。
Hilbert 变换物理意义
解析过程概念
对于⼀个实随机过程,为它的希尔伯特变换,那么定义复随机过程:
为x(t)的解析过程。
这个过程有如下特点,⾸先,实部和虚部功率谱相同,⾃相关函数相同;另外,实部和虚部的互相关函数是⼀个奇函数。其他的还有:
x (t )(t )=x ^Hilbert (x (t ))(t )=x ~x (t )+(t )
x ^R (τ)=X X
^−(τ)R ^X R (τ)=XX
^(τ)R ^X R (τ)=X ~2[R (τ)+X
j (τ)]R ^X
以及⼀个最重要的特点,就是解析信号的功率谱只有正频段,强度为原来的四倍。或者说是只有正频段且幅度值为原来的两倍:
欧拉公式(Euler‘s formula)的启发
欧拉⼀⽣中有公式和定理,这⾥说的欧拉公式指的是这个:
这个公式说明,⽤复指数信号可以表⽰成⼀个实数信号和⼀个虚数信号的和的形式。⽽且,这个实部和虚部是有关系的,⼀个是cos,⼀个是sin,两者相差pi/2,看sin和cos的傅⾥叶变换:
![这⾥写图⽚描述](imgconvert.csdnimg/aHR0cDovL2ltZy5ibG9nLmNzZG4ubmV0LzIwMTgwMjI1MDEyOTAxMzY?
x-oss-process=image/format,png)
可以看出,在正频率上和负频率上两者的相位上的先后顺序刚好相反,但是都是保持90°的差值。
看到这⾥,⼤概可以理解Hilbert变换的⽤意了吧。欧拉公式实际上是⼀种特殊的,或者说,最简单的Hilbert变换。
2πδ(ω−ω)0
复指数信号,就是等号左边的那个,频谱就是⼀个脉冲,⽽且是。只有正频率,且是两倍。虽然时域上是复数,但是在频域只
有正分量,实际上是⼀种简化。
希尔伯特变换的意义transform和convert的区别
这样看来,似乎复数信号才是完整的,⽽实信号只是在复平⾯的实轴上的⼀个投影。我们知道,解析信号可以计算包络(瞬时振幅)和瞬时相位。在上图中可以看到,实际上我们计算的包络就是⿊⾊的线围成的⽴体图形的边界在实部的投影,⽽计算这个边的投影也很简单,就是在复平⾯上的螺旋线中
的每⼀个点的模值,也就是A(t) = sqrt(x^2(t) + Hilbert(x(t))^2),⽽瞬时相位就是虚部(Hilbert变换后的)和实部(原始信号)在某⼀时间点的⽐值的arctan,瞬时频率就是它的导数。
2018年02⽉25⽇01:29:30
版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198,我们将在24小时内删除。
发表评论