希爾伯特轉換

数学信号处理中,希尔伯特变换(英語: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為這種失真的一個例子。

相關條目

參考文獻

外部連結