zoukankan      html  css  js  c++  java
  • Modeling Filters and Whitening Filters

    Colored and White Process

    White Process

    White Process,又称为White Noise(白噪声),其中white来源于白光,寓意着PSD的平坦分布,white noise指的就是在PSD上具有平坦分布(常数)的随机过程。PSD是auto-correlation的傅里叶变换,PSD上为常数意味着auto-correlation是一个位于零点上的脉冲函数。

    image回顾auto-correlation的定义:

    $R_{xx}( au) = EBig{x(t+ au)x(t)Big}$

    Auto-correlation在零点上为脉冲,除此之外都是0,这表明了各个采样点之间是相互独立的,正是由于各采样点之间相互独立,才体现出了white noise的white这一特点

    理想的White noise要求PSD是一个常数,不过由于传输媒介以及有限的监测能力的关系,通常我们所得到的PSD都受限于某个带宽范围之内,因此如果一个process的PSD在监测范围(或者我们所感兴趣的频率范围)表现得平坦的话,就可以称该process为white noise。比如,如果一个音频的PSD在20~20000Hz内表现平坦的话,我们就可以称该音频为white noise。

    i.i.d. process就是一个最简单white noise的例子。i.i.d. process中样本的PDF常见的有Uniform、Poisson、Gaussian等分布方式,如果PDF为Gaussian,则该process称为Gaussian white noise。

    Colored Process

    PSD表现得平坦的就是white noise,否则就是colored process。

    Modeling Filter

    输入为white process,输出为colored process的LTI系统被称为Modeling filter或者shaping filter。如果我们需要某个PSD为$S_{xx}(e^{jOmega})$的process,可以通过设计一个LTI系统,并且对该系统输入mean为0的Bernoulli process(PSD为常数1)来得到,所设计的系统需要使得输出的process的PSD为$S_{xx}(e^{jOmega})$。根据WSS Process在LTI系统中的性质,有

    $displaystyle{ H(e^{jOmega}) = sqrt{S_{xx}(e^{jOmega})} }$

    我们只需要对$H(e^{jOmega})$进行傅里叶逆变换就可以得到所需LTI系统的脉冲响应$h[n]$。

    image

    其中也可能会出现一些特殊情况,比如说所需要的process的PSD,即$S_{xx}(e^{jOmega})$中包含脉冲,这些情况下,PSD可以表示为脉冲PSD非脉冲PSD之和,由于傅里叶变换的线性性质,LTI系统的脉冲响应最终也可以表示为脉冲PSD的脉冲响应非脉冲PSD的脉冲响应之和。非脉冲PSD可以按照上面的步骤求解,而脉冲PSD主要分为两种情况:

    • 脉冲在零点。那么脉冲PSD的脉冲响应是一个常数,该常数的大小与脉冲PSD的大小有关。
    • 脉冲在$Omega = pm Omega_0 eq 0$上。那么脉冲PSD的脉冲响应是一个余弦函数$Acos(Omega_0 n+Theta)$,其中$A$以及$Theta$都与脉冲PSD的大小有关。

    hybrid_PSD

    ※这一段只是计算时需要,不算理解时的重点,如想加深了解,请参考常用序列的DTFT

    Whitening Filter

    输入为colored process,输出为white process的LTI系统被称为Whitening filter。如下系统框图:

    image

    输入process为$x[n]$,其PSD为$displaystyle{S_{xx}(e^{jOmega}) = mathcal{F}R_{xx}[m]}$,输出为白噪声,该白噪声的variance为$sigma_{w}^2$,mean为0。

    那么根据WSS Process在LTI系统中的性质,有

    $displaystyle{S_{ww}(e^{jOmega}) = |H(e^{jOmega})|^2S_{xx}(e^{jOmega})}$

    进一步可以推导出(详情请查看Reference的第二条链接)

    $displaystyle{|H(e^{jOmega})|^2 = frac{sigma_w^2}{S_{xx}(e^{jOmega})}}$

    Reference:

    Alan V. Oppenheim: Signals, Systems and Inference, Chapter 10:Power Spectral Density

    Where does the delta function come from if we derive autocorrelation directly?

  • 相关阅读:
    java纯数字加密解密实例
    C++手稿:std::string
    java裁剪图片
    java打开windows系统的浏览器
    Android手机无线adb
    office-word
    设计模式-享元模式(13)
    设计模式-外观模式(12)
    charles工具过滤腾讯视频播放器广告
    js将json格式的list转换为按某个字段分组的map数组
  • 原文地址:https://www.cnblogs.com/TaigaCon/p/9236625.html
Copyright © 2011-2022 走看看