白色雜訊

**白噪音**，是一種功率譜密度為常數的隨機訊號或隨機過程。即此訊號在各個頻段上的功率一致。由於白光是由各種頻率（顏色）的單色光混合而成，因而此訊號的平坦功率譜性質稱為「白色」，此訊號也因此得名為白噪音。相對的，其他不具有這一性質的噪音訊號則稱為有色噪音。

理想的白噪音具有無限頻寬，因而其能量是無限大，這在現實世界是不可能存在的。實際上，人常常將有限頻寬的平整訊號視為白噪音，以方便進行數學分析。

## 統計特性

術語白噪音也常用於表示在相關空間的自相關為 0 的空域噪音訊號，於是訊號在空間頻率域內就是「白色」的，對於角頻率域內的訊號也是這樣，例如夜空中向各個角度發散的訊號。右面的圖片顯示了計算機產生的一個有限長度的離散時間白噪音過程。

需要指出，相關性和概率分布是兩個不相關的概念。「白色」僅意味著訊號是不相關的，白噪音的定義除了要求均值為零外並沒有對訊號應當服從哪種概率分布作出任何假設。因此，如果某白噪音過程服從高斯分布，則它是「高斯白噪音」。類似的，還有泊松白噪音、柯西白噪音等。人們經常將高斯白噪音與白噪音相混同，這是不正確的認識。根據中心極限定理，高斯白噪音是許多現實世界過程的一個很好的近似，並且能夠生成數學上可以跟蹤的模型，這些模型用得如此頻繁以至於加性高斯白噪音成了一個標準的縮寫詞：AWGN。此外，高斯白噪音有著非常有用的統計學特性，因為高斯變量的獨立性與不相關性等價。

白噪音是維納過程或者布朗運動的廣義均方導數（generalized mean-square derivative）。

白噪音的數學期望為 0：

　
: $\mu _{n}=\mathbb {E} \{n(t)\}=0$

其自相關函數為狄拉克 δ 函數：

　
: $r_{nn}=\mathbb {E} \{n(t)n(t-\tau )\}=\delta (\tau )$

上式正是對白噪音的「白色」性質在時域的描述。由於隨機過程的功率譜密度是其自相關函數的傅立葉變換，而 _δ_ 函數的傅立葉變換為常數，因此白噪音的功率譜密度是平坦的。

## 噪音的顏色

也有其它「顏色」的噪音存在，最常用的有粉紅、棕色和藍色噪音。

## 應用

白噪音的應用領域之一是建築聲學，為了減弱內部空間中分散人注意力並且不希望出現的噪音（如人的交談），使用持續的低強度噪音作為背景聲音。一些緊急車輛的警報器也使用白噪音，因為白噪音能夠穿過如城市中交通噪音這樣的背景噪音並且不會引起反射，所以更加容易引起人們的注意。

在電子音樂中也有白噪音的應用，它被直接或者作為濾波器的輸入訊號以產生其它類型的噪音訊號，尤其是在音訊合成中，經常用來重現類似於鐃鈸這樣在頻域有很高噪音成分的打擊樂器。

白噪音也用來產生脈波響應。為了在一個演出地點保證音樂會或者其它演出的均衡效果，從 PA 系統發出一個瞬間的白噪音或者粉紅噪音，並且在不同的地方監測噪音訊號，這樣工程師就能夠建築物的聲學效應能夠自動地放大或者削減某些頻率，從而就可以調整總體的均衡效果以得到一個平衡的和聲，且根據網路調查，許多民眾提到人聲、雨聲等白噪音可以協助專心、入睡、放鬆心情

白噪音可以用於放大器或者電子濾波器的頻率響應測試，有時它與響應平坦的話筒或和自動均衡器一起使用。這個設計的思路是系統會產生白噪音，話筒接收到揚聲器產生的白噪音，然後在每個頻率段進行自動均衡從而得到一個平坦的響應。這種系統用在專業級的設備、高端的家庭立體聲系統或者一些高端的汽車收音機上。

白噪音也作為一些隨機數位生成器的基礎使用。

白噪音也可以用於審訊前使人迷惑，並且可能用於感覺剝奪技術的一部分。上市銷售的白噪音機器產品有私密性增強器、睡眠輔助器以及掩飾耳鳴。

## 數學定義

### 白色隨機向量

一個隨機向量 $\mathbf {w} $ 為一個白色隨機向量若且唯若它的平均值函數與自相關函數滿足以下條件：

　
: $\mu _{w}=\mathbb {E} \{\mathbf {w} \}=0$

　
: $R_{ww}=\mathbb {E} \{\mathbf {w} \mathbf {w} ^{T}\}=\sigma ^{2}\mathbf {I} $

意即它是一個平均值為零的隨機向量，並且它的自相關函數是單位矩陣的倍數。

### 白色隨機過程 (白噪音)

一個時間連續隨機過程 $w(t)$ ，其中 $t\in \mathbb {R} $ ， 為一個白噪音若且唯若它的平均值函數與自相關函數滿足以下條件：

　
: $\mu _{w}(t)=\mathbb {E} \{w(t)\}=0$

　
: $R_{ww}(t_{1},t_{2})=\mathbb {E} \{w(t_{1})w(t_{2})\}=(N_{0}/2)\delta (t_{1}-t_{2})$

意即它是一個對所有時間其平均值為零的隨機過程，並且它的自相關函數是狄拉克 δ 函數，有無限大的功率。

由上述自相關函數可推出以下的功率譜密度。

　
: $S_{xx}(\omega )=(N_{0}/2)\,\!$

由於 δ 函數的傅立葉變換為 1。而對於所有頻率來說，此功率譜密度是一樣的。因此這是對白噪音之「白色」性質在頻域的表述。

## 隨機向量變換

白色隨機向量的兩個理論應用是_類比_以及 _whitening_ 另外一個任意隨機向量。為了_類比_一個任意隨機向量，我們使用一個仔細選擇的矩陣對白色隨機向量進行變換。我們選擇的變換矩陣能夠是被變換的白色隨機向量的平均值和協方差矩陣與類比的任意向量的平均值和協方差矩陣相匹配。為了 _whiten_ 一個任意的隨機向量，我們使用仔細選擇的矩陣對它進行變換，這樣得到的隨機向量就是一個白色隨機向量。

這兩個思想在通訊和音訊領域中通道估計和通道均衡這樣的應用中是很關鍵的。這些思想在數據壓縮中也有應用。

### 類比隨機向量

假設隨機向量 $\mathbf {x} $ 有協方差矩陣 $K_{xx}$ ，由於這個矩陣是 共軛對稱和半正定，根據線性代數中的譜定理，我們可以用以下方法對角線或者分解矩陣，

　
: $\,\!K_{xx}=E\Lambda E^{T}$

其中 $E$ 是特徵向量的正交矩陣， $\Lambda $ 是特徵值的對角矩陣。

通過對白色向量 $\mathbf {w} $ 進行下面變換我們可以類比這個平均為 $\mathbf {\mu } $ 、協方差矩陣為 $K_{xx}$ 的隨機向量 $\mathbf {x} $ 的一階和二階矩量屬性：

　
: $\mathbf {x} =H\,\mathbf {w} +\mu $

其中

　
: $\,\!H=E\Lambda ^{1/2}$

這樣，這個變換輸出的期望是

　
: $\mathbb {E} \{\mathbf {x} \}=H\,\mathbb {E} \{\mathbf {w} \}+\mu =\mu $

協方差矩陣是

　
: $\mathbb {E} \{(\mathbf {x} -\mu )(\mathbf {x} -\mu )^{T}\}=H\,\mathbb {E} \{\mathbf {w} \mathbf {w} ^{T}\}\,H^{T}=H\,H^{T}=E\Lambda ^{1/2}\Lambda ^{1/2}E^{T}=K_{xx}$

### Whitening 隨機向量

whitening 一個平均值為 $\mathbf {\mu } $ 、協方差矩陣為 $K_{xx}$ 的向量 $\mathbf {x} $ 的方法是執行下面的計算：

　
: $\mathbf {w} =\Lambda ^{-1/2}\,E^{T}\,(\mathbf {x} -\mathbf {\mu } )$

這樣，這個變換輸出的期望是

　
: $\mathbb {E} \{\mathbf {w} \}=\Lambda ^{-1/2}\,E^{T}\,(\mathbb {E} \{\mathbf {x} \}-\mathbf {\mu } )=\Lambda ^{-1/2}\,E^{T}\,(\mu -\mu )=0$

協方差矩陣

　
: $\mathbb {E} \{\mathbf {w} \mathbf {w} ^{T}\}=\mathbb {E} \{\Lambda ^{-1/2}\,E^{T}\,(\mathbf {x} -\mathbf {\mu } )(\mathbf {x} -\mathbf {\mu } )^{T}E\,\Lambda ^{-1/2}\,\}$

　
: $=\Lambda ^{-1/2}\,E^{T}\,\mathbb {E} \{(\mathbf {x} -\mathbf {\mu } )(\mathbf {x} -\mathbf {\mu } )^{T}\}E\,\Lambda ^{-1/2}\,$

　
: $=\Lambda ^{-1/2}\,E^{T}\,K_{xx}E\,\Lambda ^{-1/2}$

對角線化 $K_{xx}$ 得到：

　
: $\Lambda ^{-1/2}\,E^{T}\,E\Lambda E^{T}E\,\Lambda ^{-1/2}=\Lambda ^{-1/2}\,\Lambda \,\Lambda ^{-1/2}=I$

這樣，通過上面的變換就可以將隨機向量 whiten 成平均值為 0、協方差矩陣是單位矩陣。

## 隨機訊號變換

我們將類比和白化這兩個概念推廣到連續時間隨機訊號或者隨機過程。我們創建一個濾波器用於類比，將白噪音注入其中，用輸出訊號類比任意隨機過程的一階和二階矩。對於白化，我們將任意隨機訊號注入所選濾波器中，濾波器輸出是白噪音。

### 類比連續時間隨機訊號

我們可以使用固定的平均值 $\mu $ 、協方差函數

　
: $K_{x}(\tau )=\mathbb {E} \left\{(x(t_{1})-\mu )(x(t_{2})-\mu )^{*}\right\}{\mbox{ where }}\tau =t_{1}-t_{2}$

和功率譜密度

　
: $S_{x}(\omega )=\int _{-\infty }^{\infty }K_{x}(\tau )\,e^{-j\omega \tau }\,d\tau $

類比任何廣義的穩定、連續時間隨機過程 $x(t):t\in \mathbb {R} \,\!$

我們可以使用頻域技術類比這個訊號。

由於 $K_{x}(\tau )$ 是個半正定的埃爾米特矩陣，所以 $S_{x}(\omega )$ 是實數並且若且唯若 $S_{x}(\omega )$ 滿足佩維維納標準（Paley-Wiener criterion）

　
: $\int _{-\infty }^{\infty }{\frac {\log(S_{x}(\omega ))}{1+\omega ^{2}}}\,d\omega <\infty $

時可以 factored 為

　
: $S_{x}(\omega )=|H(\omega )|^{2}=H(\omega )\,H^{*}(\omega )$

如果 $S_{x}(\omega )$ 是有理函數，我們可以將它分解成極點 - 零點格式

　
: $S_{x}(\omega )={\frac {\Pi _{k=1}^{N}(c_{k}-j\omega )(c_{k}^{*}+j\omega )}{\Pi _{k=1}^{D}(d_{k}-j\omega )(d_{k}^{*}+j\omega )}}$

選擇最小相位 (minimum phase) $H(\omega )$ 保證極點和零點都位於 S 面的左側，這樣我們就可以使用 $H(\omega )$ 作為濾波器的傳遞函數來類比 $x(t)$ 。

我們可以構建下面的線性、非時變 (time-invariant) 濾波器來類比 $x(t)$

　
: ${\hat {x}}(t)={\mathcal {F}}^{-1}\left\{H(\omega )\right\}*w(t)+\mu $

其中 $w(t)$ 是有如下一階和二階矩屬性的連續時間的白噪音：

　
: $\mathbb {E} \{w(t)\}=0$

　
: $\mathbb {E} \{w(t_{1})w^{*}(t_{2})\}=K_{w}(t_{1},t_{2})=\delta (t_{1}-t_{2})$

這樣，結果訊號 ${\hat {x}}(t)$ 與所期望的訊號 $x(t)$ 一樣有同樣的二階矩量屬性。

### 連續時間隨機訊號的白化

假設我們有一個廣義的穩定、連續時間隨機過程 $x(t):t\in \mathbb {R} \,\!$ ，與上面定義的訊號同樣的平均值 $\mu $ 、 協方差函數 $K_{x}(\tau )$ 和功率譜密度 $S_{x}(\omega )$ 。

我們可以使用頻域技術 **白化** 這個訊號，用上面的過程 factor 功率譜密度 $S_{x}(\omega )$ 。

選擇最小相位 $H(\omega )$ 得到極點和零點都位於 s 面左側，這樣就可以用下面的 inverse 濾波器 whiten $x(t)$

　
: $H_{inv}(\omega )={\frac {1}{H(\omega )}}$

選擇的最小相位濾波器保證逆濾波器穩定的。另外，必須保證 $H(\omega )$ 在所有 $\omega \in \mathbb {R} $ 上都嚴格為正，這樣 $H_{inv}(\omega )$ 就不會有任何 奇點。

白化過程的最終格式如下所示：

　
: $w(t)={\mathcal {F}}^{-1}\left\{H_{inv}(\omega )\right\}*(x(t)-\mu )$

這樣 $w(t)$ 就是一個白色噪音隨機過程，它的平均值為零、功率譜密度為

　
: $S_{w}(\omega )={\mathcal {F}}\left\{\mathbb {E} \{w(t_{1})w(t_{2})\}\right\}=H_{inv}(\omega )S_{x}(\omega )H_{inv}^{*}(\omega )={\frac {S_{x}(\omega )}{S_{x}(\omega )}}=1$

注意這個功率譜密度 對應於 $w(t)$ 的協方差函數的 δ 函數。

　
: $K_{w}(\tau )=\,\!\delta (\tau )$

## 參見

* 電子學
* 電子噪音
* δ 函數
* 獨立成分分析
* 噪音 (通訊學)
* 主成分分析
* 統計學
* 白噪音機（White noise machine）
* 粉紅噪音

## 外部連結

* A mathematical application of noise whitening of pictures - pdf
* White noise calculator, thermal noise - Voltage in microvolts, conversion to noise level in dBu and dBV and vice versa （頁面存檔備份，存於網際網路檔案館）
* Free hour-long white noise MP3 for download
* White Noise Machine