希尔伯特变换(Hilbert Transform)的物理意义

[版权声明:本文主要内容归功于以下两篇文章,链接如下:]
https://blog.csdn.net/edogawachia/article/details/79366444
https://blog.csdn.net/weixin_41608328/article/details/89322214

Hilbert变换简介

希尔伯特变换是信号处理中的一种常用手段,数学定义如下:

希尔伯特变换的定义.jpg

与卷积的概念进行对比,可以发现,上面的Hilbert变换的表达式实际上就是将原始信号和一个信号做卷积的结果。这个用来卷积的信号就是
h(t)=\dfrac{1}{\pi t}​
希尔伯特变换从本质上来看是做了一次特殊的卷积积分,其中这个卷积的冲击响应为1/πt。换句话说,Hilbert变换可以看成是将原始信号通过一个滤波器或者一个系统,这个滤波器或系统的冲击响应为h(t)。

希尔伯特变换的本质建模.jpg

对h(t)做傅里叶变换,可以得到:
H(jω)=−j\,{\rm sgn}(ω)= \left\{ \begin{array}{cc} -j & \omega \ge 0\\ j & \omega < 0 \end{array} \right.
sgn()是符号函数。从频谱上来看,这个滤波器将我们的原始信号的正频率部分乘以 -j ,而负频率部分乘以了 j 。根据之前工程数学课程上的结论,Multiple a complex number by j is to leave the modulus unaltered but to increase the argument by \pi/2。换句话说,上述系统响应的实质是在保持幅度不变的条件下,将正频率部分的相位移动了 -\pi/2 ;将负频率成分移动了 \pi/2
从上述Hilbert变换可以看出,希尔伯特变换的作用上是一个90移相器,它将信号中的正频率部分相移-90°,相当于顺时针转90°;将信号中的负频率部分相移90°,相当于逆时针转90°。希尔伯特变换不会改变实信号x(t)的振幅和能量,仅仅在相位上发生了改变分而已。

下面这个示意图很直观地表示了Hilbert变换,这里画出了对原始信号做1到4次Hilbert变换的频谱示意图:

希尔伯特变换示意图.jpg

首先,可以看到,两次希尔伯特变换后,原信号相位翻转了180°,所以,Hilbert逆变换的公式显而易见,就是将正变换加一个符号即可。另外,还可以看到,Hilbert变换四次后就变回本身了。还有其它的性质,比如:

  • 如果一个信号是两个信号的卷积,即 y = conv(v,x) ,那么 Hilbert(y) = conv(Hilbert(v),x) = conv(v,Hilbert(x))

  • x(t)Hilbert(x(t)) 的能量以及平均功率相等,相关函数和功率谱相同。

Hilbert变换物理意义

解析过程概念

对于一个实随机过程 x(t)\hat{x}(t) = Hilbert(x(t)) 为它的希尔伯特变换,那么定义复随机过程:
\tilde{x}(t) = x(t) + j\hat{x}(t)
为x(t)的解析过程

这个过程有一个最重要的特点,就是解析信号的功率谱只有正频段,强度为原来的四倍。或者说是只有正频段且幅度值为原来的两倍(幅度与功率之间还有一次平方转换关系):
S_{\tilde{X}}(\omega)= \left\{ \begin{array}{cc} 4S_X(\omega) & \omega \ge 0\\ 0 & \omega < 0 \end{array} \right.
\tilde{X}(\omega)= \left\{ \begin{array}{cc} 2X(\omega) & \omega \ge 0\\ 0 & \omega < 0 \end{array} \right.

欧拉公式(Euler‘s formula)的启发

欧拉公式:e^{j\theta}=cos(\theta)+jsin(\theta)
这个公式说明,用复指数信号可以表示成一个实数信号和一个虚数信号的和的形式。

Functions Fouier Transform
cos(\omega_0t) \pi[\delta(\omega+\omega_0)+\delta(\omega -\omega_0)]
sin(\omega_0t) j\pi[\delta(\omega+\omega_0)-\delta(\omega -\omega_0)]

可以看出,在正频率上和负频率上两者的相位上的先后顺序刚好相反,但是都是保持90°的差值。欧拉公式实际上是一种特殊的 Hilbert 变换。
复指数信号 e^{j\theta},频谱就是一个脉冲,而且是 2\pi\delta(\omega - \omega_0)。只有正频率,且是两倍。虽然时域上是复数,但是在频域只有正分量。

希尔伯特变换的意义

首先,将实数信号变换成解析信号的结果就是,把一个一维的信号变成了二维复平面上的信号,复数的模和幅角代表了信号的幅度和相位,如图所示:

希尔伯特变换物理意义图

这样看来,似乎复数信号才是完整的,而实信号只是在复平面的实轴上的一个投影。我们知道,解析信号可以计算包络(瞬时振幅)和瞬时相位。在上图中可以看到,实际上我们计算的包络就是黑色的线围成的立体图形的边界在实部的投影,而计算这个边的投影也很简单,就是在复平面上的螺旋线中的每一个点的模值,也就是 A(t) = \sqrt{x^2(t) + Hilbert(x(t))^2} ,而瞬时相位就是虚部(Hilbert变换后的)和实部(原始信号)在某一时间点的比值的 arctan ,瞬时频率就是它的导数。

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 203,324评论 5 476
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 85,303评论 2 381
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 150,192评论 0 337
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,555评论 1 273
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,569评论 5 365
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,566评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,927评论 3 395
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,583评论 0 257
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,827评论 1 297
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,590评论 2 320
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,669评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,365评论 4 318
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,941评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,928评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,159评论 1 259
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 42,880评论 2 349
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,399评论 2 342