数字信号处理方法——谱估计实验.doc

《数字信号处理方法——谱估计实验.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、实验三随机信号的功率谱估计方法报告人：####报告时间：2013.1.2一、实验目的1．利用自相关函数法和周期图法实现对随机信号的功率谱估计。2．观察数据长度、自相关序列长度、信噪比、窗函数、平均次数等对谱估计的分辨率、稳定性、主瓣宽度和旁瓣效应的影响。3．学习使用FFT提高谱估计的运算速度。4．体会非参数化功率谱估计方法的优缺点。二、实验原理与方法假设信号x(n)为平稳随机过程，其自相关序列定义为(3-1)其中E表示取数学期望，*表示共轭运算。根据定义，x(n)的功率谱密度与自相关序列存在下面关系：(3-2)(3-3

2、)然而，实际中我们很难得到准确的自相关序列，只能通过随机信号的一段样本序列来估计信号的自相关序列，进而得到信号的功率谱估计。目前，常用的线性谱估计方法有两种，即自相关函数法和周期图法。1．自相关函数法假设我们已知随机信号x(n)的M长的自相关序列{}，利用自相关函数法可以得到x(n)的功率谱估计：(3-4)利用窗函数，上式又可表达为(3-5)其中，为矩形窗函数，定义为(3-6)因此，实际上是：真正功率谱与窗函数傅立叶变换的卷积。矩形窗函数不仅降低了谱估计的分辨率，而且使谱估计产生了旁瓣。为了降低旁瓣影响，可以采用具有较

3、小旁瓣的窗函数，如Hamming窗，它定义为(3-7)这种窗函数可以有效的抑制旁瓣，但是这种方法使主瓣宽度增大，从而降低了谱估计的分辨率，这种主瓣大小和旁瓣干扰之间的矛盾在线性谱估计方法中是无法解决的。2.周期图方法假设已知随机信号的N个样本，利用周期图方法，信号x(n)的功率谱估计为(3-8)利用上述方法得到的谱估计方差与信号的功率谱平方成正比，谱估计的方差如下：（3-9）为了减小它的方差，可以将信号序列进行分段处理，然后再求各分段结果的平均，这就是平均周期图方法，即Bartlett方法。（1）Bartlett平均周

4、期图方法将一个随机序列(0≤n≤N）分成K段，每段长度为L，各段之间互不重迭，因而N=KL，可以想到，第i段的信号序列可表示为，(3-10)对于每一段的周期图又可写成,(3-11)于是，功率谱估计定义为(3-12)因此，对于固定的记录长度来讲，分段数K增大可使谱估计的方差减小，但是由于L的减小，相应的功率谱主瓣增宽，谱分辨率降低，显然，方差和分辨率也是矛盾的。除了分辨率降低以外，分段处理还会引起序列的长度有限所带来的旁瓣效应。为减小这种影响，最有效的办法是给分段序列用适当的窗函数加权，可以得到较平滑的谱估计，当然，相应

5、的分辨率也有所下降。（2）平滑平均周期图方法这时一种改进的Bartlett周期图方法，它特别适用于FFT直接计算功率谱估值。将长度为N的平稳随机信号序列x(n)分成K段，每段长度为L，即L=N/K。但这里在计算周期图之前，先用窗函数给每段序列加权，K个修正的周期图定义为,(3-13)其中U表示窗函数序列的能量，(3-14)在这种情况下，功率谱估计可按下面表达式给出：(3-15)本实验主要是利用自相关函数法和周期图方法对下面受噪声干扰扰的正弦信号进行谱估计：(3-16)其中NS为正弦个数，，和分别为第i个正弦信号的数字频

6、率、相位和幅度，随机的分布在（0，）之间，w(n)为零均值方差等于的复高斯白噪声。三、实验内容和分析1．仔细阅读有关线性谱估计的内容，根据给出的框图编制自相关函数法谱估计的程序。运行程序，输入选择矩形窗。观察谱峰位置是否正确（注意：由于窗效应可能引起谱估计的非正定）。开始答：a.程序流图如3.1所示：（程序代码见附录一）输入参数：数据长度N，自相关函数个数M,平均次数K信号产生：输入正弦个数Ns，每个正弦信号的数字频率、相位和幅度，白噪声信号的方差，按照公式（3-16）产生复正弦加白噪声信号的N个采样图表3.1实验流程

7、框图周期图方法输入FFT点数NF，按照公式（3-11）、（3-12）、（3-13）计算NF点的功率谱。自相关函数法由N个x(n)估计出自相关序列（M长），并对此自相关序列加矩形窗或Hamming窗，利用公式（3-5）计算[0,2*pi]之间的128个功率谱抽样点结束图3.1程序流程图b.根据代码运行结果如图3.2：图3.2矩形窗进行功率谱估计的结果分析：由图可知，功率谱的峰值在0.6处，与理论计算值相同。2．观察并记录参数变化对谱估计性能的影响。（1）改变M=5，其它输入同步骤1，观察功率谱估计的主瓣宽度和旁瓣大小随自

8、相关序列长度的变化情况。（2）选择窗函数为Hamming窗，其它输入同步骤1，观察不同的窗函数对谱估计性能的影响。（3）改变，其它输入同步骤1，观察初始相位的变化对谱估计性能的影响。（4）改变，其它输入同步骤1，观察信噪比变化对谱估计性能的影响。（5）改变N=10，，其它输入同步骤1，结合(4)的内容，观察数据长度及信噪比对谱估计