希爾伯特轉換

在數學和信號處理中，希爾伯特變換（Hilbert transform）是一個對函數 u(t) 產生定義域相同的函數 H(u)(t) 的線性算子。

希爾伯特變換在信號處理中很重要，能夠導出信號 u(t) 的解析表示。這就意味著將實信號 u(t) 拓展到複平面，使其滿足柯西－黎曼方程。

例如，希爾伯特變換引出了傅立葉分析中給定函數的調和共軛，也就是調和分析。等價地說，它是奇異積分算子與乘子 (傅立葉分析)的一個例子。

希爾伯特變換最初只對周期函數（也就是圓上的函數）有定義，在這種情況下它就是與希爾伯特核的卷積。然而更常見的情況下，對於定義在實直線 R（上半平面的邊界）上的函數，希爾伯特變換是指與柯西核卷積。希爾伯特變換與帕利-維納定理有著密切的聯繫，帕利-維納定理是將上半平面內的全純函數與實直線上的函數的傅立葉變換相聯繫起來的另一種結果。

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

<image>

## 定義

$u$ 的希爾伯特變換可以認為是 $u(t)$ 與函數 $h(t) = \frac{1}{\pi t}:=\delta(jt)$ 的卷積。由於 $h(t)$ 是不可積的，定義卷積的積分不收斂。因而希爾伯特變換是使用柯西主值（這裡記為$p.v.$）定義的。準確說來，函數（或信號） $u(t)$ 的希爾伯特變換是：

$$H(u)(t) = \operatorname{p.v.} \int_{-\infty}^\infty u(\tau) h(t-\tau)\, d\tau = \frac{1}{\pi} \ \operatorname{p.v.} \int_{-\infty}^\infty \frac{u(\tau)}{t-\tau}\, d\tau$$

假設此積分作為主值存在。這就是 u 與緩增分布 p.v. 1/Πt 的卷積（由於；參見）。另外，通過改變變量，主值積分可以顯式地Zygmund寫為：

$$H(u)(t) = \frac{2}{\pi}\lim_{\varepsilon\rightarrow 0}\int_{\varepsilon}^\infty \frac{u(t + \tau) - u(t - \tau)}{2\tau}\,d\tau.$$

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

$$H(H(u))(t) = u(-t)$$

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

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

### 頻率響應

希爾伯特轉換之頻率響應由傅立葉變換給出：

$H(\omega ) = (-i\cdot \sgn(\omega)) \cdot \mathcal{F}\{h\}(\omega)\,$

其中

- $\mathcal{F}$是傅立葉變換。
- i (有時寫作j )是虛數單位。
- $\omega \,$是角頻率，以及
- $$\sgn(\omega) =\begin{cases} \ \ 1, & \mbox{for } \omega > 0,\\ \ \ 0, & \mbox{for } \omega = 0,\\ -1, & \mbox{for } \omega < 0, \end{cases}$$

即為符號函數。

既然：

$\mathcal{F}\{\widehat s\}(\omega) = H(\omega )\cdot \mathcal{F}\{s\}(\omega)$,

希爾伯特轉換會將負頻率成分$s(t)\,$偏移+90°，而正頻率成分偏移−90°。

### 反（逆）希爾伯特轉換

我們也注意到：$H^2(\omega ) = 1^{-}\,$。因此將上面方程式乘上$-H(\omega )\,$，可得到：

$$\mathcal{F}\{s\}(\omega) = -H(\omega )\cdot \mathcal{F}\{\widehat s\}(\omega)$$

從中，可以看出反（逆）希爾伯特轉換

$$s(t) = (h *h *h * \widehat s)(t) = \mathcal{H}^3\{\widehat s\}(t).\,$$

## 希爾伯特轉換表格

訊號 $u(t)\,$｜希爾伯特轉換 $H(u)(t)$

$\sin(t)$ ｜$-\cos(t)$

$\cos(t)$ ｜$\sin(t)\,$

$\exp \left( i t \right)$｜$- i \exp \left( i t \right)$

$\exp \left( -i t \right)$｜$i \exp \left( -i t \right)$

$1 \over t^2 + 1$｜$t \over t^2 + 1$

$e^{-t^2}$｜$2\pi^{-1/2}F(t)$ 參見道森積分

Sinc函數 $\sin(t) \over t$｜$1 - \cos(t)\over t$

矩形函數 $\sqcap(t)$｜${1 \over \pi} \log \left | {t + {1 \over 2} \over t - {1 \over 2}} \right |$

狄拉克δ函數 $\delta(t) \,$｜${1 \over \pi t}:=\delta(jt)$

指示函數 $\chi_{[a,b]}(t) \,$｜$\frac{1}{\pi}\log \left \vert \frac{t - a}{t - b}\right \vert$

Notes

常數之希爾伯特轉換為零

## 特性

### 邊界

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

$$\|Hu\|_p \le C_p\| u\|_p$$

對所有 u∈Lp(R)。這個定理由所推得；請一併參見。

最佳常數Cp可由下列算式得到:

$$C_p = \begin{cases} <pre> \tan \frac{\pi}{2p} & \text{for } 1 < p \leq 2\\ \cot \frac{\pi}{2p} & \text{for } 2 < p < \infty </pre> \end{cases}$$

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

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

$$S_R f = \int_{-R}^{R}\hat{f}({\xi})e^{2\pi i x\xi}\,d\xi$$

請參見Duoandikoetxea。

### 反自伴性

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

$$\langle Hu, v \rangle = \langle u, -Hv \rangle$$

對 u ∈ Lp(R) 且 v ∈ Lq(R) Titchmarsh.

### 逆轉換

希爾伯特轉換為一反-對合 Titchmarsh，意即

$$H(H(u)) = -u$$

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

$$H^{-1} = -H$$

### 微分

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

$$H\left(\frac{du}{dt}\right) = \frac{d}{dt}H(u)$$

此一特性亦可迭代

$$H\left(\frac{d^ku}{dt^k}\right) = \frac{d^k}{dt^k}H(u)$$

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

### 旋積

希爾伯特轉換可表示為與一緩增分布之旋積 Duistermaat

$$h(t) = \text{p.v. }\frac{1}{\pi t}$$

因此可如此表示

$$H(u) = h*u$$

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

$$H(u)(t) = \frac{d}{dt}\left(\frac{1}{\pi} (u*\log|\cdot|)(t)\right)$$

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

$$H(u*v) = H(u)*v = u*H(v)$$

若 u 和 v 為緊支撐分布，則此項論述嚴格成立，在這個狀況下

$$h*(u*v) = (h*u)*v = u*(h*v)$$

### 不變性

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

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

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

$$\displaystyle{U_{g}^{-1}f(x) = (cx + d)^{-1} f\left({ax + b \over cx + d}\right),\,\,\,g = \begin{pmatrix} a & b \\ c & d \end{pmatrix}}$$

## 希爾伯特轉換例子

注意：有些作者，例如Bracewell，將我們的$-\mathcal{H}$當作其正轉換的定義。這樣的結果就是下表右行要乘上一個負號。

## 離散希爾伯特轉換

<image>

<image>

<image>

<image>

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

$$H(u)[n] = \scriptstyle{DTFT}^{-1} \displaystyle \{U(\omega)\cdot \sigma_H(\omega)\}$$

其中

$$\sigma_H(\omega)\ \stackrel{\mathrm{def}}{=}\ \begin{cases} <pre> e^{+i\pi/2}, & -\pi < \omega < 0 \\ e^{-i\pi/2}, & 0 < \omega < \pi\\ 0, & \omega = -\pi, 0, \pi </pre> \end{cases}$$

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

$$H(u)[n] = u[n] * h[n]$$

其中

$$h[n]\ \stackrel{\mathrm{def}}{=}\ \scriptstyle{DTFT}^{-1} \big \{\displaystyle \sigma_H(\omega)\big \} = \begin{cases} <pre> 0, & \mbox{for }n\mbox{ even}\\ \frac2{\pi n} & \mbox{for }n\mbox{ odd} </pre> \end{cases}$$

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

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

$$h_N[n]\ \stackrel{\text{def}}{=}\ \sum_{m=-\infty}^{\infty} h[n - mN]$$

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

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

## 外部連結

- The Discrete Hilbert Transform; A Brief Tutorial_w236
- Derivation of the boundedness of the Hilbert transform
- Mathworld Hilbert transform — Contains a table of transforms
- Analytic Signals and Hilbert Transform Filters
- Mathias Johansson, "The Hilbert transform" a student level summary to Hilbert transformation. (via www.archive.org)
- GS256 Lecture 3: Hilbert Transformation, an entry level introduction to Hilbert transformation. (via www.archive.org)