希尔伯特变换

在數學和訊號處理中，**希爾伯特轉換**（英語：**Hilbert transform**）是一個對函數 _u_(_t_) 產生定義域相同的函數 _H_(_u_)(_t_) 的線性算子。

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

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

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

## 定義

 $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_ 的摺積（由於 Schwartz (1950)；參見 Pandey (1996，Chapter 3)）。另外，通過改變變數，主值積分可以顯式地（Zygmund 1968，§XVI.1）寫為：

　
: $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 \operatorname {sgn}(\omega ))\cdot {\mathcal {F}}\{h\}(\omega )\,$

其中

 * ${\mathcal {F}}$ 是傅立葉變換，
* _i_ (有時寫作 _j_ ) 是虛數單位，
* $\omega \,$ 是角頻率，以及
* $\operatorname {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).\,$

## 希爾伯特轉換表格

**Notes**

常數之希爾伯特轉換為零

## 特性

### 邊界

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

　
: $\|Hu\|_{p}\leq C_{p}\|u\|_{p}$

對所有 _u_∈_L_p(**R**)。這個定理由 Riesz (1928，VII) 所推得；請一併參見 Titchmarsh (1948，Theorem 101)。
最佳常數 _Cp_ 可由下列算式得到:

　
: $C_{p}={\begin{cases}\tan {\frac {\pi }{2p}}&{\text{for }}1<p\leq 2\\\cot {\frac {\pi }{2p}}&{\text{for }}2<p<\infty \end{cases}}$

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

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

　
: $S_{R}f=\int _{-R}^{R}{\hat {f}}({\xi })e^{2\pi ix\xi }\,d\xi $

請參見（Duoandikoetxea 2000，p.59）。

### 反自伴性

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

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

對 _u_ ∈ _Lp_(**R**) 且 _v_ ∈ _L_q(**R**) （Titchmarsh 1948，Theorem 102）.

### 逆轉換

希爾伯特轉換為一反 - 對合 （Titchmarsh 1948，p.120），意即

　
: $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^{k}u}{dt^{k}}}\right)={\frac {d^{k}}{dt^{k}}}H(u)$

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

### 旋積

希爾伯特轉換可表示為與一緩增分布之旋積 （Duistermaat & Kolk 2010，p.211）

　
: $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)$

### 不變性

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

* 可與算子 _T_aƒ(_x_) = ƒ(_x_ + _a_) 交換，對所有實數 _a_
* 可與算子 _M_λƒ(_x_) = ƒ(λ_x_) 交換，對所有 λ > 0
* 可與鏡射 _R_ƒ(_x_) = ƒ(−x) 反交換

實際上，有更大一部分的算子可與希爾伯特轉換交換。群組 SL (2,**R**) 由么正算符 _U_g 可在空間 _L_2(**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}}$ 當作其正轉換的定義。這樣的結果就是下表右行要乘上一個負號。

## 離散希爾伯特轉換

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

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

其中

　
: $\sigma _{H}(\omega )\ {\stackrel {\mathrm {def} }{=}}\ {\begin{cases}e^{+i\pi /2},&-\pi <\omega <0\\e^{-i\pi /2},&0<\omega <\pi \\0,&\omega =-\pi ,0,\pi \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}0,&{\mbox{for }}n{\mbox{ even}}\\{\frac {2}{\pi n}}&{\mbox{for }}n{\mbox{ odd}}\end{cases}}$

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

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

　
: $h_{N}[n]\ {\stackrel {\text{def}}{=}}\ \sum _{m=-\infty }^{\infty }h[n-mN]$

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

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

## 相關條目

* 摺積
* 希爾伯特 - 黃轉換

## 參考文獻

* Bargmann, V., Irreducible unitary representations of the Lorentz group, Ann. of Math., 1947, **48** (3): 568–640, JSTOR 1969129, doi:10.2307/1969129
* Bedrosian, E., A Product Theorem for Hilbert Transforms (PDF), Rand Corporation Memorandum, December 1962, (RM-3439-PR) [2016-08-05], （原始內容存檔 (PDF) 於 2021-02-25）
* Benedetto, John J. Harmonic analysis and applications. Boca Raton, FL: CRC Press. 1996. ISBN 0849378796.
* Bitsadze, A.V., Boundary value problems of analytic function theory, Hazewinkel, Michiel (編), 数学百科全书，Springer, 2001, ISBN 978-1-55608-010-4 .
* Bracewell, R., The Fourier Transform and Its Applications 3rd, McGraw–Hill, 2000, ISBN 0-07-116043-4 .
* Calderón, A.P.; Zygmund, A., On the existence of certain singular integrals, Acta Mathematica, 1952, **88** (1): 85–139, doi:10.1007/BF02392130 .
* Carlson, Crilly, and Rutledge, Communication Systems 4th, 2002, ISBN 0-07-011127-8 .
* Duoandikoetxea, J., Fourier Analysis, American Mathematical Society, 2000, ISBN 0-8218-2172-5 .
* Duistermaat, J.J.; Kolk, J.A.C. Kolk, Distributions, Birkhäuser, 2010, ISBN 978-0-8176-4672-1, doi:10.1007/978-0-8176-4675-2 .
* Duren, P., Theory of $H^{p}$ -Spaces, New York: Academic Press, 1970 .
* Fefferman, C., Characterizations of bounded mean oscillation, Bull. Amer. Math. Soc., 1971, **77** (4): 587–588, MR 0280994, doi:10.1090/S0002-9904-1971-12763-5 .
* Fefferman, C.; Stein, E.M., Hp spaces of several variables, Acta Math., 1972, **129**: 137–193, MR 0447953, doi:10.1007/BF02392215 .
* Gel'fand, I.M.; Shilov, G.E., Generalized Functions, Vol. 2, Academic Press, 1967 .
* Grafakos, Loukas, An Elementary Proof of the Square Summability of the Discrete Hilbert Transform, American Mathematical Monthly (Mathematical Association of America), 1994, **101** (5): 456–458, JSTOR 2974910, doi:10.2307/2974910 .
* Grafakos, Loukas, Classical and Modern Fourier Analysis, Pearson Education, Inc.: 253–257, 2004, ISBN 0-13-035399-X .
* Hardy, G. H.; Littlewood, J. E.; Polya, G., Inequalities, Cambridge: Cambridge University Press, 1952, ISBN 0-521-35880-9 .
* Hilbert, David, Grundzüge einer allgemeinen Theorie der linearen Integralgleichungen, Chelsea Pub. Co., 1953
* Kak, Subhash, The discrete Hilbert transform, Proc. IEEE, 1970, **58**: 585–586 .
* Kak, Subhash, Number theoretic Hilbert transform, Circuits Systems Signal Processing, 2014, **33**: 2539–2548 .
* Khvedelidze, B.V., Hilbert transform, Hazewinkel, Michiel (編), 数学百科全书，Springer, 2001, ISBN 978-1-55608-010-4 .
* King, Frederick W., Hilbert Transforms **2**, Cambridge: Cambridge University Press: 453, 2009, ISBN 978-0-521-51720-1 .
* Kress, Rainer, Linear Integral Equations, New York: Springer-Verlag: 91, 1989, ISBN 3-540-50616-0 .
* Lang, Serge, SL(2,**R**), Graduate Texts in Mathematics **105**, Springer-Verlag, 1985, ISBN 0-387-96198-4
* Pandey, J.N., The Hilbert transform of Schwartz distributions and applications, Wiley-Interscience, 1996, ISBN 0-471-03373-1
* Pichorides, S., On the best value of the constants in the theorems of Riesz, Zygmund, and Kolmogorov, Studia Mathematica, 1972, **44**: 165–179
* Riesz, Marcel, Sur les fonctions conjuguées, Mathematische Zeitschrift, 1928, **27** (1): 218–244, doi:10.1007/BF01171098
* Rosenblum, Marvin; Rovnyak, James, Hardy classes and operator theory, Dover, 1997, ISBN 0-486-69536-0
* Schwartz, Laurent, Théorie des distributions, Paris: Hermann, 1950 .
* Schreier, P.; Scharf, L., Statistical signal processing of complex-valued data: the theory of improper and noncircular signals, Cambridge University Press, 2010
* Stein, Elias, Singular integrals and differentiability properties of functions, Princeton University Press, 1970, ISBN 0-691-08079-8 .
* Stein, Elias; Weiss, Guido, Introduction to Fourier Analysis on Euclidean Spaces, Princeton University Press, 1971, ISBN 0-691-08078-X .
* Sugiura, Mitsuo, Unitary Representations and Harmonic Analysis: An Introduction, North-Holland Mathematical Library **44** 2nd, Elsevier, 1990, ISBN 0444885935
* Titchmarsh, E, Reciprocal formulae involving series and integrals, Mathematische Zeitschrift, 1926, **25** (1): 321–347, doi:10.1007/BF01283842 .
* Titchmarsh, E, Introduction to the theory of Fourier integrals 2nd, Oxford University: Clarendon Press, 19481986, ISBN 978-0-8284-0324-5 .
* Zygmund, Antoni, Trigonometric series 2nd, Cambridge University Press, 19681988, ISBN 978-0-521-35885-9 .

## 外部連結

* 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
* 埃里克・韋斯坦因. Titchmarsh theorem. MathWorld.
* 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)