Privacy Policy Cookie Policy Terms and Conditions 白雜訊 - Wikipedia

白雜訊

维基百科,自由的百科全书

Image:03wiki-zn-frontpage-icon.gif白雜訊正在翻译。欢迎您积极翻译与修订
目前已翻译95%,原文在en:White noise

Template:Colors of noise

白噪声频谱
白噪声频谱

白-{A|zh-cn:噪声;zh-tw:噪聲}-白雜訊,是一種功率頻譜密度為常數的隨機信號或随机过程。換句話說,此信號在各個频段上的功率是一樣的,由于白光是由各種頻率(颜色)的单色光混合而成,因而此信号的這種具有平坦功率谱的性质被称作是“白色的”,此信号也因此被称作白噪声。相对的,其他不具有这一性质的噪声信号被称为有色噪声。

理想的白噪声具有無限頻寬,因而其能量是無限大,這在现实世界是不可能存在的。实际上,我們常常將有限頻寬的平整訊號視為白噪音,因為這讓我們在數學分析上更加方便。

目录

[编辑] 統計特性

Image:Whitenoise.png
千分之四秒的白噪声
白噪声过程实现实例
白噪声过程实现实例

术语白噪声也常用于表示在相关空间的自相关为 0 的空域噪声信号,于是信号在空间频率域内就是“白色”的,对于角频率域内的信号也是这样,例如夜空中向各个角度发散的信号。右面的图片显示了计算机产生的一个有限长度、离散时间的白噪声过程。

然而,在时域不相关并不能够约束信号的可能值,任何数值分布都是可能的,尽管它必须有值为 0 的直流分量。例如,一个只能有 1 和 0 值的二进制信号如果 0 和 1 的序列统计上不相关的话它就是白色信号。如正态分布这样的连续分布的噪声当然就是白噪声。

人们通常不正确地假定高斯噪声(例如高斯幅度分布的噪声——参见正态分布)一定是白噪声,然而,二者之间并不是必然的。高斯性指的是信号值分布的方式,术语'白色'表示在两个不同时间的相关性,这并不依赖于信号的幅度分布。

Image:Noise.jpg
频谱图上显示的左边的粉红噪声和右边的白噪声

这样我们就能找到高斯白噪声,也可以找到Poisson、Cauchy等白噪声。需要注意的是分布必须是无限变化的,这样,“高斯”和“白色”这两个词经常用系统的数学模型来说明。高斯白噪声是许多现实世界状态的一个很好的近似,并且能够生成数学上可以跟踪的模型,这些模型用得如此频繁以至于叠加高斯白噪声成了一个标准的缩写词:AEGN。高斯白噪声有非常有用的统计学特性,它的值都是相互独立的(参见统计独立)。

白噪声是维纳过程或者布朗运动的generalized mean-square derivative。

白噪声的数学期望为0:

\mu_n =  \mathbb{E} \{ n(t) \} = 0

而其自相关函数δ函数:

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

由于随机过程的功率谱是其自相关的傅里叶变换,而δ函数的傅里叶变换为常数,因此这是对白噪声之“白色”性质在时域的表述。

[编辑] 噪声的颜色

主條目:Colors of noise

也有其它“颜色”的噪声存在,最常用的有粉红棕色和蓝色噪声。

[编辑] 應用

白噪声的应用领域之一是建筑声学,为了减弱内部空间中分散人注意力并且不希望出现的噪声(如人的交谈),使用持续的低强度噪声作为背景声音。一些紧急车辆的警报器也使用白噪声,因为白噪声能够穿过如城市中交通噪声这样的背景噪声并且不会引起反射,所以更加容易引起人们的注意。

电子音乐中也有白噪声的应用,它被直接或者作为滤波器的输入信号以产生其它类型的噪声信号,尤其是在音频合成中,经常用来重现类似于铙钹这样在频域有很高噪声成分的打击乐器。

白噪声也用来产生冲击响应。为了在一个演出地点保证音乐会或者其它演出的均衡效果,从 P A 系统发出一个瞬间的白噪声或者粉红噪声,并且在不同的地方监测噪声信号,这样工程师就能够建筑物的声学效应能够自动地放大或者削减某些频率,从而就可以调整总体的均衡效果以得到一个平衡的和声。

白噪声可以用于放大器或者电子滤波器的频率响应测试,有时它与响应平坦的话筒或和自动均衡器一起使用。这个设计的思路是系统会产生白噪声,话筒接收到扬声器产生的白噪声,然后在每个频率段进行自动均衡从而得到一个平坦的响应。这种系统用在专业级的设备、高端的家庭立体声系统或者一些高端的汽车收音机上。

白噪声也作为一些随机数字生成器的基础使用。

白噪声也可以用于审讯前使人迷惑,并且可能用于感觉剥夺技术的一部分。上市销售的白噪声机器产品有私密性增强器、睡眠辅助器以及掩饰耳鸣

[编辑] 數學定義

[编辑] 白色隨機向量

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

\mu_w =  \mathbb{E}\{ \mathbf{w} \} = 0
R_{ww} = \mathbb{E}\{ \mathbf{w} \mathbf{w}^T\} = \sigma^2 \mathbf{I}

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

[编辑] 白色隨機過程(白雜訊)

一個時間連續隨機過程w(t) where 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)\} = \sigma^2 \delta(t_1 - t_2)

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

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

S_{xx}(\omega) = \sigma^2 \,\!

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

[编辑] 随机向量变换

白色随机向量的两个理论应用是模拟以及whitening另外一个任意随机向量。为了模拟一个任意随机向量,我们使用一个仔细选择的矩阵对白色随机向量进行变换。我们选择的变换矩阵能够是被变换的白色随机向量的平均值和协方差矩阵与模拟的任意向量的平均值和协方差矩阵相匹配。为了whiten一个任意的随机向量,我们使用仔细选择的矩阵对它进行变换,这样得到的随机向量就是一个白色随机向量。

这两个思想在通信和音频领域中通道估计和通道均衡这样的应用中是很关键的。这些思想在数据压缩中也有应用。

[编辑] 模拟随机向量

假设随机向量 \mathbf{x}协方差矩阵 Kxx,由于这个矩阵是 Hermitian 对称(Hermitian symmetric)和半正定(en:positive semidefinite),根据线性代数中的谱定理,我们可以用以下方法对角线或者 factor 矩阵,

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

其中 E特征向量正交矩阵Λ特征值对角矩阵

通过对白色向量 \mathbf{w} 进行下面变换我们可以模拟这个平均\mathbf{\mu}、协方差矩阵为Kxx的随机向量 \mathbf{x} 的一阶和二阶moment)属性:

\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}协方差矩阵Kxx 的向量 \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}

对角线化 Kxx 得到:

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

这样,通过上面的变换就可以将随机向量 whiten 成平均值为0、协方差矩阵是单位矩阵。

[编辑] 随机信号变换

我们将模拟和whitening这两个概念推广到连续时间随机信号或者随机过程。我们创建一个滤波器用于模拟,将白噪声注入其中,用输出信号模拟任意随机过程的一阶和二阶矩。对于whitening,我们将任意随机信号注入所选滤波器中,滤波器输出是白噪声。

[编辑] 模拟连续时间随机信号

将白噪声注入线性时不变滤波器中模拟任意随机过程的一阶和二阶矩
将白噪声注入线性时不变滤波器中模拟任意随机过程的一阶和二阶矩

我们可以使用固定的平均值 μ协方差函数

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}\,\!

我们可以使用频域技术模拟这个信号。

由于 Kx(τ) Hermitian 对称(:en:Hermitian symmetric)和半正定,所以 Sx(ω)实数并且当且仅当 Sx(ω) 满足 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)

如果 Sx(ω)有理函数,我们可以将它 factor 成极-零格式

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)}

选择最小相位 H(ω) 保证极点和零点都位于S 面的左侧,这样我们就可以使用 H(ω) 作为滤波器的传递函数来模拟 x(t)

我们可以构建下面的线性、非時變滤波器来模拟 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) 输入一个线性时不变滤波器,滤波器将  x(t)白化为白噪声
任意随机过程 x(t) 输入一个线性时不变滤波器,滤波器将 x(t)白化为白噪声

假设我们有一个广义的稳定、连续时间随机过程 x(t) : t \in \mathbb{R}\,\!,与上面定义的信号同样的平均值 μ协方差函数 Kx(τ) 和功率谱密度 Sx(ω)

我们可以使用频域技术 白化 这个信号,用上面的过程 factor 功率谱密度 Sx(ω)

选择最小相位 H(ω) 得到极点和零点都位于s 面左侧,这样就可以用下面的 inverse 滤波器 whiten x(t)

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

选择的最小相位滤波器保证逆滤波器稳定的。另外,必须保证 H(ω) 在所有 \omega \in \mathbb{R} 上都严格为正,这样 Hinv(ω) 就不会有任何 奇点。

whitening 过程的最终格式如下所示:

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)协方差函数的 delta 函数。

K_w(\tau) = \,\!\delta (\tau)

[编辑] 参见

[编辑] 外部链接

THIS WEB:

aa - ab - af - ak - als - am - an - ang - ar - arc - as - ast - av - ay - az - ba - bar - bat_smg - be - bg - bh - bi - bm - bn - bo - bpy - br - bs - bug - bxr - ca - cbk_zam - cdo - ce - ceb - ch - cho - chr - chy - closed_zh_tw - co - cr - cs - csb - cu - cv - cy - da - de - diq - dv - dz - ee - el - eml - en - eo - es - et - eu - fa - ff - fi - fiu_vro - fj - fo - fr - frp - fur - fy - ga - gd - gl - glk - gn - got - gu - gv - ha - haw - he - hi - ho - hr - hsb - ht - hu - hy - hz - ia - id - ie - ig - ii - ik - ilo - io - is - it - iu - ja - jbo - jv - ka - kg - ki - kj - kk - kl - km - kn - ko - kr - ks - ksh - ku - kv - kw - ky - la - lad - lb - lbe - lg - li - lij - lmo - ln - lo - lt - lv - map_bms - mg - mh - mi - mk - ml - mn - mo - mr - ms - mt - mus - my - mzn - na - nah - nap - nds - nds_nl - ne - new - ng - nl - nn - no - nov - nrm - nv - ny - oc - om - or - os - pa - pag - pam - pap - pdc - pi - pih - pl - pms - ps - pt - qu - rm - rmy - rn - ro - roa_rup - roa_tara - ru - ru_sib - rw - sa - sc - scn - sco - sd - se - searchcom - sg - sh - si - simple - sk - sl - sm - sn - so - sq - sr - ss - st - su - sv - sw - ta - te - test - tet - tg - th - ti - tk - tl - tlh - tn - to - tokipona - tpi - tr - ts - tt - tum - tw - ty - udm - ug - uk - ur - uz - ve - vec - vi - vls - vo - wa - war - wo - wuu - xal - xh - yi - yo - za - zea - zh - zh_classical - zh_min_nan - zh_yue - zu

Static Wikipedia 2008 (no images)

aa - ab - af - ak - als - am - an - ang - ar - arc - as - ast - av - ay - az - ba - bar - bat_smg - bcl - be - be_x_old - bg - bh - bi - bm - bn - bo - bpy - br - bs - bug - bxr - ca - cbk_zam - cdo - ce - ceb - ch - cho - chr - chy - co - cr - crh - cs - csb - cu - cv - cy - da - de - diq - dsb - dv - dz - ee - el - eml - en - eo - es - et - eu - ext - fa - ff - fi - fiu_vro - fj - fo - fr - frp - fur - fy - ga - gan - gd - gl - glk - gn - got - gu - gv - ha - hak - haw - he - hi - hif - ho - hr - hsb - ht - hu - hy - hz - ia - id - ie - ig - ii - ik - ilo - io - is - it - iu - ja - jbo - jv - ka - kaa - kab - kg - ki - kj - kk - kl - km - kn - ko - kr - ks - ksh - ku - kv - kw - ky - la - lad - lb - lbe - lg - li - lij - lmo - ln - lo - lt - lv - map_bms - mdf - mg - mh - mi - mk - ml - mn - mo - mr - mt - mus - my - myv - mzn - na - nah - nap - nds - nds_nl - ne - new - ng - nl - nn - no - nov - nrm - nv - ny - oc - om - or - os - pa - pag - pam - pap - pdc - pi - pih - pl - pms - ps - pt - qu - quality - rm - rmy - rn - ro - roa_rup - roa_tara - ru - rw - sa - sah - sc - scn - sco - sd - se - sg - sh - si - simple - sk - sl - sm - sn - so - sr - srn - ss - st - stq - su - sv - sw - szl - ta - te - tet - tg - th - ti - tk - tl - tlh - tn - to - tpi - tr - ts - tt - tum - tw - ty - udm - ug - uk - ur - uz - ve - vec - vi - vls - vo - wa - war - wo - wuu - xal - xh - yi - yo - za - zea - zh - zh_classical - zh_min_nan - zh_yue - zu -

Static Wikipedia 2007:

aa - ab - af - ak - als - am - an - ang - ar - arc - as - ast - av - ay - az - ba - bar - bat_smg - be - bg - bh - bi - bm - bn - bo - bpy - br - bs - bug - bxr - ca - cbk_zam - cdo - ce - ceb - ch - cho - chr - chy - closed_zh_tw - co - cr - cs - csb - cu - cv - cy - da - de - diq - dv - dz - ee - el - eml - en - eo - es - et - eu - fa - ff - fi - fiu_vro - fj - fo - fr - frp - fur - fy - ga - gd - gl - glk - gn - got - gu - gv - ha - haw - he - hi - ho - hr - hsb - ht - hu - hy - hz - ia - id - ie - ig - ii - ik - ilo - io - is - it - iu - ja - jbo - jv - ka - kg - ki - kj - kk - kl - km - kn - ko - kr - ks - ksh - ku - kv - kw - ky - la - lad - lb - lbe - lg - li - lij - lmo - ln - lo - lt - lv - map_bms - mg - mh - mi - mk - ml - mn - mo - mr - ms - mt - mus - my - mzn - na - nah - nap - nds - nds_nl - ne - new - ng - nl - nn - no - nov - nrm - nv - ny - oc - om - or - os - pa - pag - pam - pap - pdc - pi - pih - pl - pms - ps - pt - qu - rm - rmy - rn - ro - roa_rup - roa_tara - ru - ru_sib - rw - sa - sc - scn - sco - sd - se - searchcom - sg - sh - si - simple - sk - sl - sm - sn - so - sq - sr - ss - st - su - sv - sw - ta - te - test - tet - tg - th - ti - tk - tl - tlh - tn - to - tokipona - tpi - tr - ts - tt - tum - tw - ty - udm - ug - uk - ur - uz - ve - vec - vi - vls - vo - wa - war - wo - wuu - xal - xh - yi - yo - za - zea - zh - zh_classical - zh_min_nan - zh_yue - zu

Static Wikipedia 2006:

aa - ab - af - ak - als - am - an - ang - ar - arc - as - ast - av - ay - az - ba - bar - bat_smg - be - bg - bh - bi - bm - bn - bo - bpy - br - bs - bug - bxr - ca - cbk_zam - cdo - ce - ceb - ch - cho - chr - chy - closed_zh_tw - co - cr - cs - csb - cu - cv - cy - da - de - diq - dv - dz - ee - el - eml - en - eo - es - et - eu - fa - ff - fi - fiu_vro - fj - fo - fr - frp - fur - fy - ga - gd - gl - glk - gn - got - gu - gv - ha - haw - he - hi - ho - hr - hsb - ht - hu - hy - hz - ia - id - ie - ig - ii - ik - ilo - io - is - it - iu - ja - jbo - jv - ka - kg - ki - kj - kk - kl - km - kn - ko - kr - ks - ksh - ku - kv - kw - ky - la - lad - lb - lbe - lg - li - lij - lmo - ln - lo - lt - lv - map_bms - mg - mh - mi - mk - ml - mn - mo - mr - ms - mt - mus - my - mzn - na - nah - nap - nds - nds_nl - ne - new - ng - nl - nn - no - nov - nrm - nv - ny - oc - om - or - os - pa - pag - pam - pap - pdc - pi - pih - pl - pms - ps - pt - qu - rm - rmy - rn - ro - roa_rup - roa_tara - ru - ru_sib - rw - sa - sc - scn - sco - sd - se - searchcom - sg - sh - si - simple - sk - sl - sm - sn - so - sq - sr - ss - st - su - sv - sw - ta - te - test - tet - tg - th - ti - tk - tl - tlh - tn - to - tokipona - tpi - tr - ts - tt - tum - tw - ty - udm - ug - uk - ur - uz - ve - vec - vi - vls - vo - wa - war - wo - wuu - xal - xh - yi - yo - za - zea - zh - zh_classical - zh_min_nan - zh_yue - zu