Hilbert變換簡介

希爾伯特變換是訊號處理中的一種常用手段，數學定義如下：

與卷積的概念進行對比，可以發現，上面的Hilbert變換的表示式實際上就是將原始訊號和一個訊號做卷積的結果。這個用來卷積的訊號就是

h(t)=1πth(t)=1πt

因此，Hilbert變換可以看成是將原始訊號通過一個濾波器，或者一個系統，這個系統的衝擊響應為h(t)。

對h(t)做傅立葉變換，可以得到：

H(jω)=−jsgn(ω)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)，x^(t)=Hilbert(x(t))x^(t)=Hilbert(x(t))為它的希爾伯特變換，那麼定義復隨機過程：

x~(t)=x(t)+x^(t)x~(t)=x(t)+x^(t)

為x(t)的解析過程。

這個過程有如下特點，首先，實部和虛部功率譜相同，自相關函式相同；另外，實部和虛部的互相關函式是一個奇函式。其他的還有：

RX^X(τ)=−R^X(τ)RX^X(τ)=−R^X(τ)

RXX^(τ)=R^X(τ)RXX^(τ)=R^X(τ)

RX~(τ)=2[RX(τ)+jR^X(τ)]RX~(τ)=2[RX(τ)+jR^X(τ)]

以及一個最重要的特點，就是解析訊號的功率譜只有正頻段，強度為原來的四倍。或者說是隻有正頻段且幅度值為原來的兩倍：

尤拉公式（Euler‘s formula）的啟發

尤拉一生中有公式和定理，這裡說的尤拉公式指的是這個：

這個公式說明，用復指數訊號可以表示成一個實數訊號和一個虛數訊號的和的形式。而且，這個實部和虛部是有關係的，一個是cos，一個是sin，兩者相差pi/2，看sin和cos的傅立葉變換：

可以看出，在正頻率上和負頻率上兩者的相位上的先後順序剛好相反，但是都是保持90°的差值。

看到這裡，大概可以理解Hilbert變換的用意了吧。尤拉公式實際上是一種特殊的，或者說，最簡單的Hilbert變換。

復指數訊號，就是等號左邊的那個，頻譜就是一個脈衝，而且是2πδ(ω−ω0)2πδ(ω−ω0)。只有正頻率，且是兩倍。雖然時域上是複數，但是在頻域只有正分量，實際上是一種簡化。

希爾伯特變換的意義

首先，將實數訊號變換成解析訊號的結果就是，把一個一維的訊號變成了二維複平面上的訊號，複數的模和幅角代表了訊號的幅度和相位，如圖所示：（圖片來源：https://zhuanlan.zhihu.com/p/25250010）

這樣看來，似乎複數訊號才是完整的，而實訊號只是在複平面的實軸上的一個投影。我們知道，解析訊號可以計算包絡（瞬時振幅）和瞬時相位。在上圖中可以看到，實際上我們計算的包絡就是黑色的線圍成的立體圖形的邊界在實部的投影，而計算這個邊的投影也很簡單，就是在複平面上的螺旋線中的每一個點的模值，也就是A(t) = sqrt(x^2(t) + Hilbert(x(t))^2)，而瞬時相位就是虛部（Hilbert變換後的）和實部（原始訊號）在某一時間點的比值的arctan，瞬時頻率就是它的導數。