随机信号的功率谱估计方法.pdf

《随机信号的功率谱估计方法.pdf》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、精品文档数字信号处理II——随机信号的功率谱估计方法1欢迎下载。精品文档一、实验目的1．利用自相关函数法和周期图法实现对随机信号的功率谱估计。2．观察数据长度、自相关序列长度、信噪比、窗函数、平均次数等对谱估计的分辨率、稳定性、主瓣宽度和旁瓣效应的影响。3．学习使用FFT提高谱估计的运算速度。4．体会非参数化功率谱估计方法的优缺点。二、实验原理与方法假设信号x(n)为平稳随机过程，其自相关序列定义为:(m)Ex*(n)x(nm)(0.1)其中E表示取数学期望，*表示取共轭。根据定义，x(n)的功率谱密度P(w)与自

2、相关序列(m)存在如下关系：P()(m)ejm(0.2)m1(m)P()ejmd(0.3)2然而，实际中我们很难得到准确的自相关序列(m)，只能通过随机信号的一段样本序列来估计信号的自相关序列，进而得到信号的功率谱估计。目前常用的线性谱估计方法有两种：自相关函数法和周期图方法，本实验将对这两种方法分别予以讨论。1．自相关函数法假设已知随机信号x(n)的N个观测样本，则其自相关序列可以用下式进行估计:(m)1N

3、m

4、1x*(n)x(nm)

5、m

6、N1(0.4)N

7、m

8、n0当

9、仅使用长度为2M-1的自相关序列时，对其进行傅立叶变换即可得到功率谱估计如下:2欢迎下载。精品文档M1Pˆ()ˆ(m)ejm(0.5)mM1其中M为加窗长度，WRec(m)为矩形窗函数，定义如下:M1,

10、m

11、MWRec(m)(0.6)M0,

12、m

13、M因此，Pˆ(w)在一定程度上可以看作是“真正的功率谱P(w)”与窗函数傅立叶变换的卷积。矩形窗函数不仅降低了谱估计的分辨率，而且使谱估计产生了旁瓣，旁瓣效应使那些处于旁瓣附近功率较小的频率分量被淹没掉。为了降低旁瓣影响，可以采用具有较小旁瓣的窗函数，如Hamm

14、ing窗，它定义为：m0.540.46cos,

15、m

16、MWHam(m)M(0.7)M0,

17、m

18、M这种窗函数可以有效的抑制旁瓣，但此时主瓣宽度增大，从而降低了谱估计的分辨率。这种主瓣和旁瓣之间的矛盾在非参数化功率谱估计方法中是无法解决的。2．周期图方法假设已知随机信号x(n)的N个样本，利用周期图方法，信号x(n)的功率谱估计为:ˆ(1N1P)

19、x(n)ejm

20、2(0.8)Nn0利用上述方法得到的谱估计方差与信号x(n)的功率谱平方[P(w)]2成正比，为了减小它的方差，可以将信号序列进行分段处理，然后再

21、求各分段结果的平均，这就是平均周期图方法，即Bartlett方法,它特别适用于FFT直接计算功率谱估值。（1）Bartlett平均周期图方法将一个随机序列x(n)(0nN)分成K段，每段长度为L，各段之间互不重叠，因而NLK。可以想到，第i段的信号序列可表示为：3欢迎下载。精品文档x(n)x(n(i1)L)0n

22、x(n)ejm

23、21iK(0.10)iLin0于是，功率谱估计定义为：1KPˆ()Iˆ()(0.11)Kii1对于固定

24、的记录长度来讲，分段数K增大可使谱估计的方差减小，但是由于L的减小，相应的功率谱主瓣增宽，谱分辨率降低，显然，方差和分辨率也是矛盾的。除了分辨率降低以外，分段处理还会引起序列的长度有限所带来的旁瓣效应。为减小这种影响，最有效的办法是给分段序列用适当的窗函数加权，可以得到较平滑的谱估计，当然，相应的分辨率也有所下降。（2）平滑平均周期图方法这是一种改进的Bartlett周期图方法，将长度为N的平稳随机信号序列x(n)分成K段，每段长度为L，即LN/K。这里在计算周期图之前，先用窗函数W(n)给每段序列x(n)加权，K个修正的周期图定义

25、为：Li21L-1IˆW()=x(n)W(n)e-jm1iK(0.12)iLUiLn=0其中U表示窗函数序列每个样本的平均能量，即：21L-1U=W(n)(0.13)LLn=0在这种情况下，功率谱估计可按下面表达式给出：1KPˆ()=IˆW()(0.14)Kii=1本实验主要是利用自相关函数法和周期图方法对下面受噪声干扰扰的正弦信号进行谱估计：Nsax(n)=ej(n+)+w(n)(0.15)iiii=1其中N为正弦个数，，和a分别为第i个正弦信号的数字频率、相位siii4欢迎下载。精品文档和幅度，随机分

26、布在[0,2)之间，w(n)为零均值，方差等于2的复高斯白iw噪声。三、实验内容及步骤1、仔细阅读有关功率谱估计的内容，按照图1给出的框图编制自相关函数法谱估计程序。运行程序，输入N=100,M=10,N=1,=0