希爾伯特轉換

數學訊號處理中,希爾伯特轉換(英語:Hilbert transform)是一個對函數 u(t) 產生定義域相同的函數 H(u)(t) 的線性算子

希爾伯特轉換在訊號處理中很重要,能夠導出訊號 u(t) 的解析表示。這就意味着將實訊號 u(t) 拓展到複數平面,使其滿足柯西-黎曼方程。 例如,希爾伯特轉換引出了傅利葉分析中給定函數的調和共軛,也就是調和分析。等價地說,它是奇異積分算子傅利葉乘子英語Multiplier (Fourier analysis)的一個例子。

希爾伯特轉換最初只對週期函數(也就是上的函數)有定義,在這種情況下它就是與希爾伯特核摺積。然而更常見的情況下,對於定義在實直線 R上半平面邊界)上的函數,希爾伯特轉換是指與柯西核摺積。希爾伯特轉換與帕利-維納定理英語Paley–Wiener theorem有着密切的聯繫,帕利-維納定理是將上半平面內的全純函數與實直線上的函數的傅利葉轉換相聯繫起來的另一種結果。

希爾伯特轉換是以大衛·希爾伯特來命名的,他首先引入了該算子來解決全純函數黎曼–希爾伯特問題的一個特殊情況。

希爾伯特轉換結果(紅色)與原來的訊號——方波(藍色)

定義

 希爾伯特轉換可以認為是   與函數  摺積。由於   是不可積的,定義摺積的積分不收斂。因而希爾伯特轉換是使用柯西主值(這裏記為 )定義的。準確說來,函數(或訊號)   的希爾伯特轉換是:

 

假設此積分作為主值存在。這就是 u緩增分佈 p.v. 1/πt 的摺積(由於Schwartz (1950);參見Pandey (1996,Chapter 3))。另外,通過改變變量,主值積分可以顯式地(Zygmund 1968,§XVI.1)寫為:

 

若希爾伯特轉換接連用在函數 u 上兩次,結果就是負 u

 

假設定義兩次迭代的積分都收斂。特別地,反轉換是 −H。可以通過考慮 u(t) 的傅利葉轉換的希爾伯特轉換效應看出這一事實(參見下面的與傅利葉轉換的關係)。

上半平面解析函數,希爾伯特轉換描述了邊界值的實部與虛部之間的關係。也就是說,如果 f(z) 是在 Im z > 0 平面內的解析函數,而 u(t) = Re f(t + 0·i ),假設希爾伯特轉換存在,則 Im f(t + 0·i ) = H(u)(t) 取決於一個相加性常數。

頻率響應

希爾伯特轉換之頻率響應傅利葉變換給出:

   

其中

  •  是傅利葉變換,
  • i (有時寫作j )是虛數單位
  •  角頻率,以及
  •  

即為符號函數

既然:

 ,

希爾伯特轉換會將負頻率成分 偏移+90°,而正頻率成分偏移−90°。

反(逆)希爾伯特轉換

我們也注意到: 。因此將上面方程式乘上 ,可得到:

 

從中,可以看出反(逆)希爾伯特轉換

 

希爾伯特轉換表格

訊號
 
希爾伯特轉換[fn 1]
 
  [fn 2]  
  [fn 2]  
   
   
   
    參見道森積分
Sinc函數
 
 
矩形函數
 
 
狄拉克δ函數
 
 
指示函數
 
 
Notes
  1. ^ Some authors (e.g., Bracewell) use our −H as their definition of the forward transform. A consequence is that the right column of this table would be negated.
  2. ^ 2.0 2.1 The Hilbert transform of the sin and cos functions can be defined in a distributional sense, if there is a concern that the integral defining them is otherwise conditionally convergent. In the periodic setting this result holds without any difficulty.

常數之希爾伯特轉換為零

特性

邊界

若 1<p<∞,則 Lp(R)之希爾伯特轉換為一有界算子,表示存在一常數Cp使得

 

對所有 uLp(R)。這個定理由Riesz (1928,VII)所推得;請一併參見Titchmarsh (1948,Theorem 101)。 最佳常數Cp可由下列算式得到:

 

這個結果由(Pichorides 1972)所推得;請一併參見Grafakos (2004,Remark 4.1.8)。上述最佳常數計算方式應用在週期性希爾伯特轉換一樣成立。

希爾伯特轉換的邊界指的是 Lp(R) 對稱級數運算子對於在 Lp(R) 之中 f 的收斂

 

請參見(Duoandikoetxea 2000,第59頁)。

反自伴性

希爾伯特轉換為一反自伴算子,連結 Lp(R) 與其對偶空間 Lq(R),其中 pq赫爾德共軛且 1 < p,q < ∞. 以符號表示

 

u ∈ Lp(R) 且 v ∈ Lq(R) (Titchmarsh 1948,Theorem 102).

逆轉換

希爾伯特轉換為一反-對合 (Titchmarsh 1948,第120頁),意即

 

假定每一轉換皆完整定義過。由於 H 保存了 Lp(R)空間,這特別代表希爾伯特轉換在 Lp(R) 上是可逆的,且

 

微分

正式上,一個式子其希爾伯特轉換的微分即為其微分的希爾伯特轉換,意即這兩者是可以交換的線性算子

 

此一特性亦可迭代

 

給定 u 以及其前k次微分皆屬於Lp(R) (Pandey 1996,§3.3)空間,此項論述為嚴格成立。在頻域上可以輕易驗證這件事情,由於微分在頻域上即為與 ω 之乘積。

旋積

希爾伯特轉換可表示為與一緩增分佈旋積 (Duistermaat & Kolk 2010,第211頁)

 

因此可如此表示

 

然而,事前此特性可能只有對緊支撐之分佈 u定義。由於緊支撐函數在 Lp 上是稠密的,因此此項特性可能嚴格成立。另一角度來看,也可使用 h(t) 其微分之特性來證明

 

在大部分的用途,希爾伯特轉換可被視為是一旋積。舉例而言,旋積與希爾伯特轉換具備下列可交換的特性

 

uv 為緊支撐分佈,則此項論述嚴格成立,在這個狀況下

 

不變性

希爾伯特轉換在空間 L2(R) 上有下列特性

  • 可與算子 Taƒ(x) = ƒ(x + a) 交換,對所有實數 a
  • 可與算子 Mλƒ(x) = ƒ(λx) 交換,對所有 λ > 0
  • 可與鏡射 Rƒ(x) = ƒ(−x) 反交換

實際上,有更大一部分的算子可與希爾伯特轉換交換。群組 SL(2,R) 由么正算符 Ug 可在空間 L2(R) 上由以下式子表示

 

希爾伯特轉換例子

注意:有些作者,例如Bracewell,將我們的 當作其正轉換的定義。這樣的結果就是下表右行要乘上一個負號。

離散希爾伯特轉換

 
圖 1: 頻寬被限制在95%奈奎斯特頻率之濾波器頻率響應
 
圖 2: 高通頻率響應之希爾伯特轉換濾波器
 
圖 3.
 
圖 4. cos(wt)函數之希爾伯特轉換為 sin(wt)。此圖顯示了sin(wt)函數與一個利用MATLAB函式庫 hilbert(·)計算之近似希爾伯特轉換的差異

對於一離散函數 u[n],以及其 離散傅利葉轉換 函數 U(ω),可推得其希爾伯特轉換為:

 

其中

 

此外,根據摺積定律,另一個相等的方程式為:

 

其中

 

當摺積經由數值運算後,一FIR 近似將取代h[n],如 圖 1所示,可以見到頻率響應在通帶之兩端(0與奈奎斯特頻率)的陡降,形成一帶通濾波器。其中高頻部分可藉由一FIR濾波器回復,如 圖 2所示。然而實際上,一個經過適當取樣的 u[n] 序列在高頻部分已經不具有可用的分量。當脈衝響應持續越久,低頻部分也可以被回復。

用FIR近似h[n]的時候,交疊儲存法是一個對於很長的u[n] 序列做摺積運算的有效方法。有時候陣列FFT{h[n]}會被σH(ω)相對應之取樣序列所取代。如此將會有與週期疊加函數做摺積之效果:

 

圖 3比較了hN[n]之半週期與一相同長度分量之h[n]。兩者之間之差異與兩者之長度皆不短於區段長度(N)之現象為失真的來源,且失真可經由增加區段長度與交疊參數來有效減少。

MATLAB中有一函數 hilbert(u,N),此函數會回傳一複數序列,其中虛部序列為 u[n]之離散希爾伯特轉換近似,實部序列為原本輸入之序列,所以這樣的複數輸出等於是 u[n]的分析訊號。與前述類似, hilbert(u, N) 只使用來自 sgn(ω)分佈的取樣,因此是與 hN[n] 的摺積。如前段所述,失真可藉由選擇比實際之u[n]序列更大的N與捨棄適當數量的輸出取樣來有效減少。圖 4為這種失真的一個例子。

相關條目

參考文獻

外部連結