白雜訊功率譜
白雜訊[1],或稱白噪音,是一種功率譜密度為常數的隨機訊號或隨機過程。即此訊號在各個頻段上的功率一致。由於白光是由各種頻率(顏色)的單色光混合而成,因而此訊號的平坦功率譜性質稱為「白色」,此訊號也因此得名為白雜訊。相對的,其他不具有這一性質的雜訊訊號則稱為有色雜訊。
理想的白雜訊具有無限頻寬,因而其能量是無限大,這在現實世界是不可能存在的。實際上,人常常將有限頻寬的平整訊號視為白雜訊,以方便進行數學分析。
統計特性
白雜訊過程現實實例
術語白雜訊也常用於表示在相關空間的自相關為0的空域雜訊訊號,於是訊號在空間頻率域內就是「白色」的,對於角頻率域內的訊號也是這樣,例如夜空中向各個角度發散的訊號。右面的圖片顯示了計算機產生的一個有限長度的離散時間白雜訊過程。
需要指出,相關性和機率分布是兩個不相關的概念。「白色」僅意味著訊號是不相關的,白雜訊的定義除了要求均值為零外並沒有對訊號應當服從哪種機率分布作出任何假設。因此,如果某白雜訊過程服從高斯分布,則它是「高斯白雜訊」。類似的,還有卜瓦松白雜訊、柯西白雜訊等。人們經常將高斯白雜訊與白雜訊相混同,這是不正確的認識。根據中心極限定理,高斯白雜訊是許多現實世界過程的一個很好的近似,並且能夠生成數學上可以跟蹤的模型,這些模型用得如此頻繁以至於加性高斯白雜訊成了一個標準的縮寫詞:AWGN。此外,高斯白雜訊有著非常有用的統計學特性,因為高斯變量的獨立性與不相關性等價。
白雜訊是維納過程或者布朗運動的廣義均方導數(generalized mean-square derivative)。
白雜訊的數學期望為0:
![{\displaystyle \mu _{n}=\mathbb {E} \{n(t)\}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3efdc6a7bf9b2ab855f4577d63bf181a218e4572)
其自相關函數為狄拉克δ函數:
![{\displaystyle r_{nn}=\mathbb {E} \{n(t)n(t-\tau )\}=\delta (\tau )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d7dafa18f96c8533c015ab3a1c7504ed5062a2f2)
上式正是對白雜訊的「白色」性質在時域的描述。由於隨機過程的功率譜密度是其自相關函數的傅立葉變換,而δ函數的傅立葉變換為常數,因此白雜訊的功率譜密度是平坦的。
雜訊的顏色
頻譜圖上顯示的左邊的粉紅雜訊和右邊的白雜訊
也有其它「顏色」的雜訊存在,最常用的有粉紅、棕色和藍色雜訊。
應用
白雜訊的應用領域之一是建築聲學,為了減弱內部空間中分散人注意力並且不希望出現的雜訊(如人的交談),使用持續的低強度雜訊作為背景聲音。一些緊急車輛的警報器也使用白雜訊,因為白雜訊能夠穿過如城市中交通雜訊這樣的背景雜訊並且不會引起反射,所以更加容易引起人們的注意。
在電子音樂中也有白雜訊的應用,它被直接或者作為濾波器的輸入訊號以產生其它類型的雜訊訊號,尤其是在音訊合成中,經常用來重現類似於鐃鈸這樣在頻域有很高雜訊成分的打擊樂器。
白雜訊也用來產生脈波響應。為了在一個演出地點保證音樂會或者其它演出的均衡效果,從PA系統發出一個瞬間的白雜訊或者粉紅雜訊,並且在不同的地方監測雜訊訊號,這樣工程師就能夠建築物的聲學效應能夠自動地放大或者削減某些頻率,從而就可以調整總體的均衡效果以得到一個平衡的和聲,且根據網路調查,許多民眾提到人聲、雨聲等白噪音可以協助專心、入睡、放鬆心情[2]
白雜訊可以用於放大器或者電子濾波器的頻率響應測試,有時它與響應平坦的話筒或和自動均衡器一起使用。這個設計的思路是系統會產生白雜訊,話筒接收到揚聲器產生的白雜訊,然後在每個頻率段進行自動均衡從而得到一個平坦的響應。這種系統用在專業級的設備、高端的家庭立體聲系統或者一些高端的汽車收音機上。
白雜訊也作為一些隨機數位生成器的基礎使用。
白雜訊也可以用於審訊前使人迷惑,並且可能用於感覺剝奪技術的一部分。上市銷售的白雜訊機器產品有私密性增強器、睡眠輔助器以及掩飾耳鳴。
數學定義
白色隨機向量
一個隨機向量
為一個白色隨機向量若且唯若它的平均值函數與自相關函數滿足以下條件:
![{\displaystyle \mu _{w}=\mathbb {E} \{\mathbf {w} \}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fb20f54c25061b67beb8a28dd8c8e2ff0750f043)
![{\displaystyle R_{ww}=\mathbb {E} \{\mathbf {w} \mathbf {w} ^{T}\}=\sigma ^{2}\mathbf {I} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/fd1e852bd0fc1cf42c44fa82c260d7c2b7b16f45)
意即它是一個平均值為零的隨機向量,並且它的自相關函數是單位矩陣的倍數。
白色隨機過程(白雜訊)
一個時間連續隨機過程
,其中
, 為一個白雜訊若且唯若它的平均值函數與自相關函數滿足以下條件:
![{\displaystyle \mu _{w}(t)=\mathbb {E} \{w(t)\}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/584f462decfdb64d102f1728c1fe440d7cc7e156)
![{\displaystyle R_{ww}(t_{1},t_{2})=\mathbb {E} \{w(t_{1})w(t_{2})\}=(N_{0}/2)\delta (t_{1}-t_{2})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b19561828a9761136033f4c2249d4049d6ab6c0b)
意即它是一個對所有時間其平均值為零的隨機過程,並且它的自相關函數是狄拉克δ函數,有無限大的功率。
由上述自相關函數可推出以下的功率譜密度。
![{\displaystyle S_{xx}(\omega )=(N_{0}/2)\,\!}](https://wikimedia.org/api/rest_v1/media/math/render/svg/067d7834b6557497cf51bb0a45997ddfd8719dc9)
由於δ函數的傅立葉變換為1。而對於所有頻率來說,此功率譜密度是一樣的。因此這是對白雜訊之「白色」性質在頻域的表述。
隨機向量變換
白色隨機向量的兩個理論應用是類比以及白化另外一個任意隨機向量。為了類比一個任意隨機向量,我們使用一個仔細選擇的矩陣對白色隨機向量進行變換。我們選擇的變換矩陣能夠是被變換的白色隨機向量的平均值和協方差矩陣與類比的任意向量的平均值和協方差矩陣相匹配。為了白化一個任意的隨機向量,我們使用仔細選擇的矩陣對它進行變換,這樣得到的隨機向量就是一個白色隨機向量。
這兩個思想在通訊和音訊領域中通道估計和通道均衡這樣的應用中是很關鍵的。這些思想在數據壓縮中也有應用。
類比隨機向量
假設隨機向量
有協方差矩陣
,由於這個矩陣是 共軛對稱和半正定,根據線性代數中的譜定理,我們可以用以下方法對角線或者分解矩陣,
![{\displaystyle \,\!K_{xx}=E\Lambda E^{T}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/93ee0e451ccc4362ffa54fbc3350f171f002e158)
其中
是特徵向量的正交矩陣,
是特徵值的對角矩陣。
通過對白色向量
進行下面變換我們可以類比這個平均為
、協方差矩陣為
的隨機向量
的一階和二階矩量屬性:
![{\displaystyle \mathbf {x} =H\,\mathbf {w} +\mu }](https://wikimedia.org/api/rest_v1/media/math/render/svg/11af3eca93cab3870e69ba4f2575b4df44275e85)
其中
![{\displaystyle \,\!H=E\Lambda ^{1/2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8ce234c646e3d100ae3d5c6132531ad547a5d652)
這樣,這個變換輸出的期望是
![{\displaystyle \mathbb {E} \{\mathbf {x} \}=H\,\mathbb {E} \{\mathbf {w} \}+\mu =\mu }](https://wikimedia.org/api/rest_v1/media/math/render/svg/513d9080b8adcbc0b909062f72df695e3924f297)
協方差矩陣是
![{\displaystyle \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}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/244ea43387f2d88c0cc0a7999a4331b9666129ed)
白化隨機向量
白化一個平均值為
、協方差矩陣為
的向量
的方法是執行下面的計算:
![{\displaystyle \mathbf {w} =\Lambda ^{-1/2}\,E^{T}\,(\mathbf {x} -\mathbf {\mu } )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e5ca41690d4d817d32aa5dd14fbdecf5603addce)
這樣,這個變換輸出的期望是
![{\displaystyle \mathbb {E} \{\mathbf {w} \}=\Lambda ^{-1/2}\,E^{T}\,(\mathbb {E} \{\mathbf {x} \}-\mathbf {\mu } )=\Lambda ^{-1/2}\,E^{T}\,(\mu -\mu )=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fa98cc005a496cfca66fac69cabbedb2799da309)
協方差矩陣
![{\displaystyle \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}\,\}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c80d64030d72ad6dd0241a164e26958472a4aa87)
![{\displaystyle =\Lambda ^{-1/2}\,E^{T}\,\mathbb {E} \{(\mathbf {x} -\mathbf {\mu } )(\mathbf {x} -\mathbf {\mu } )^{T}\}E\,\Lambda ^{-1/2}\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9ac912caccdc956e448631e306c4544a208bfd97)
![{\displaystyle =\Lambda ^{-1/2}\,E^{T}\,K_{xx}E\,\Lambda ^{-1/2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bcb784b1c5743d3448be0d7a0a89b620f638780f)
對角線化
得到:
![{\displaystyle \Lambda ^{-1/2}\,E^{T}\,E\Lambda E^{T}E\,\Lambda ^{-1/2}=\Lambda ^{-1/2}\,\Lambda \,\Lambda ^{-1/2}=I}](https://wikimedia.org/api/rest_v1/media/math/render/svg/26815e55cdc8c7d0905af4913fc2e6c7357b0a8d)
這樣,通過上面的變換就可以將隨機向量白化為平均值為0、協方差矩陣是單位矩陣的隨機向量。
隨機訊號變換
我們將類比和白化這兩個概念推廣到連續時間隨機訊號或者隨機過程。我們創建一個濾波器用於類比,將白雜訊注入其中,用輸出訊號類比任意隨機過程的一階和二階矩。對於白化,我們將任意隨機訊號注入所選濾波器中,濾波器輸出是白雜訊。
類比連續時間隨機訊號
將白雜訊注入線性時不變濾波器中類比任意隨機過程的一階和二階矩
我們可以使用固定的平均值
、協方差函數
![{\displaystyle K_{x}(\tau )=\mathbb {E} \left\{(x(t_{1})-\mu )(x(t_{2})-\mu )^{*}\right\}{\mbox{ where }}\tau =t_{1}-t_{2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d6660d58d3a9e9da67450fa75d37a3f271282916)
和功率譜密度
![{\displaystyle S_{x}(\omega )=\int _{-\infty }^{\infty }K_{x}(\tau )\,e^{-j\omega \tau }\,d\tau }](https://wikimedia.org/api/rest_v1/media/math/render/svg/658623f1778001a97507a905ea5463ba7ba3dbe6)
類比任何廣義的穩定、連續時間隨機過程
我們可以使用頻域技術類比這個訊號。
由於
是個半正定的埃爾米特矩陣,所以
是實數並且若且唯若
滿足佩維維納標準(Paley-Wiener criterion)
![{\displaystyle \int _{-\infty }^{\infty }{\frac {\log(S_{x}(\omega ))}{1+\omega ^{2}}}\,d\omega <\infty }](https://wikimedia.org/api/rest_v1/media/math/render/svg/58594c1d2e481a82835c03aa100b1779a01304cc)
時可以 factored 為
![{\displaystyle S_{x}(\omega )=|H(\omega )|^{2}=H(\omega )\,H^{*}(\omega )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/937361cf759828d9d0ff3c5db6ba8b89fdf148ca)
如果
是有理函數,我們可以將它分解成極點-零點格式
![{\displaystyle 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 )}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a7b6819848332c29a45174dc7c7d90d4666d1df7)
選擇最小相位 (minimum phase)
保證極點和零點都位於S 面的左側,這樣我們就可以使用
作為濾波器的傳遞函數來類比
。
我們可以構建下面的線性、非時變 (time-invariant) 濾波器來類比
![{\displaystyle {\hat {x}}(t)={\mathcal {F}}^{-1}\left\{H(\omega )\right\}*w(t)+\mu }](https://wikimedia.org/api/rest_v1/media/math/render/svg/71d57f69b1585734cadf734372894a9c90c2e607)
其中
是有如下一階和二階矩屬性的連續時間的白雜訊:
![{\displaystyle \mathbb {E} \{w(t)\}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f0c7fe41c90979cf8d9d546f1d9bb8081375daf6)
![{\displaystyle \mathbb {E} \{w(t_{1})w^{*}(t_{2})\}=K_{w}(t_{1},t_{2})=\delta (t_{1}-t_{2})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/210c77cd1ce104bdf865a08e59f5798219330001)
這樣,結果訊號
與所期望的訊號
一樣有同樣的二階矩量屬性。
連續時間隨機訊號的白化
任意隨機過程 x(t) 輸入一個線性時不變濾波器,濾波器將 x(t)白化為白雜訊
假設我們有一個廣義的穩定、連續時間隨機過程
,與上面定義的訊號同樣的平均值
、 協方差函數
和功率譜密度
。
我們可以使用頻域技術 白化 這個訊號,用上面的過程 factor 功率譜密度
。
選擇最小相位
得到極點和零點都位於s 面左側,這樣就可以用下面的反濾波器白化
![{\displaystyle H_{inv}(\omega )={\frac {1}{H(\omega )}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4f63bb40600765b540683864f1d521fc1e7a4455)
選擇的最小相位濾波器保證逆濾波器穩定的。另外,必須保證
在所有
上都嚴格為正,這樣
就不會有任何 奇異點。
白化過程的最終格式如下所示:
![{\displaystyle w(t)={\mathcal {F}}^{-1}\left\{H_{inv}(\omega )\right\}*(x(t)-\mu )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6cd5eea98aa1559281d02838321452c5511282bb)
這樣
就是一個白色雜訊隨機過程,它的平均值為零、功率譜密度為
![{\displaystyle 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}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2b42d716cd9396d3fc990be8bd4993d821ba67b3)
注意這個功率譜密度 對應於
的協方差函數的 δ函數。
![{\displaystyle K_{w}(\tau )=\,\!\delta (\tau )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/42cf8782c2bebac76ab8f291a3918485bfc5e4e5)
參見
參考資料
外部連結