[信息与通信]第6讲 谱估计2传统法.pptVIP

4.7 传统功率谱估计 4.7 传统功率谱估计 传统功率谱估计包括：1898年由Schster提出的周期图（也称样本谱或直接法）和1958年由Blackman 与Tukey提出的BT法（也称自相关法或间接法）。 4.7.1 间接法 间接法是由随机信号N个观察值 ，估计出自相关函数 ，然后求其傅里叶变换作为功率谱的估计。 4.7.2 直接法 设有限长实序列 的傅里叶变换为： 周期图的计算公式为： 周期图 的特点： 1）是 的周期为 的周期函数； 2）非负，具有真实谱的特性； 3）与复数序列傅里叶变换的相角无关； 周期图与间接法（BT法）都以FFT为基础。 4.7.3 直接法和间接法的关系 周期图 令m=k-n，则k=n+m，得 可见周期图法就是BT法中M=N-1时的功率谱估计。 4.8 传统谱估计方法的估计质量 4.8.1 均值 M=N-1时， 又对于零均值的随机过程有 所以 即直接法与间接法估计是功率谱的有偏（渐进无偏）估计。 自相关函数无偏估计 对应的功率谱为： 期望 可见自相关函数估计无偏，但功率谱估计却有偏。 功率谱估计有偏的原因： 三角窗 矩形窗 4.8.2 方差 功率谱估计方差的精确表示式非常繁冗，通常用较易理解的近似表示式。 当M≠N-1时，可以证明BT法是一致估计。 当M＝N-1时，周期图和BT法的方差： 白色高斯过程情况： 是平稳白色高斯信号的方差。 可见此时周期图和BT法不是功率谱的一致估计。 当频率 和 时（k和l是整数） 当 时，上式等于零，说明频率上彼此间隔( ) 整数倍的周期图值是不相关的。于是当N增加，这些 协方差为零的不相关频率取样将聚集一起，趋于常 数，和白色信号的真实谱相一致。 因此，记录长度N变长时，传统谱估计起伏加快。 非白色高斯过程的情况： 非白色随机序列可认为是白色噪声通过某线性系统后 的输出。于是线性系统输出的功率谱密度是输入功率 谱密度和系统幅度频率响应平方的乘积。 色噪声的传统功率谱估计： 白噪声的传统功率谱估计： 线性系统幅度频率响应的平方是 ，因此 此时谱估计方差为： 不是功率谱密度的一致估计。 随着N增加，方差正比于频谱的平方，即在谱的真 值附近摆动相当大，需要进行改进。 4.8.3 谱分辨率与“泄漏”现象 泄漏现象可能引起功率谱估计产生畸变，对正弦分量的信号检测带来有害的影响，因为弱信号的主瓣很容易被强信号泄漏到邻近的副瓣所淹没，从而造成谱的模糊与失真。 谱分辨率：窗的变换与实际信号变换的卷积意味着使信号的主瓣变宽，增加的宽度由窗的主瓣宽度决定。由于主瓣宽度一般近似等于观测时间的倒数，所以对于短观测数据序列，其功率谱估计的分辨率比较低。 谱分辨率与栅栏效应： 栅栏效应通过补零可以使谱线增密，使离散谱线的 包络逐渐接近连续频率的周期图，对克服离散谱线 的模糊不定性有所改善，但无法改善谱分辨率。这 是由于补零没有改变窗函数的宽度，因此窗函数频 谱的主瓣宽度没有减小的原因。 传统谱估计的弱点是频率分辨率的限制。 非均匀窗的使用可以降低旁瓣泄漏，但引起窗口主瓣的加宽，因而谱分辨率下降。 4.9 传统谱估计方法的改进措施 改进周期图谱估计性能的方法包括：平均和平滑。 4.9.1 多个周期图取平均估计功率谱密度 （平均周期图法） 把N点数据分成互不重叠和互相重叠的K段，分别对 每段数据求周期图再取平均。分段越多，方差越小， 估值偏差越大。 1．多个不重叠数据段周期图的平均 Bartlett提出把数据序列分成K段，每段有M个取样，N=KM，形成序列段： 对应K个周期图 如果 时， 很小，则可假设各周期图彼此 独立无关。定义平均周期图： 平均周期图法的期望值：