功率谱估计Levinson 递推法和 Burg 法

《功率谱估计Levinson 递推法和 Burg 法》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、数字信号处理实验报告姓名：学号：日期：2015.12.141.实验任务信号为两个正弦信号加高斯白噪声，各正弦信号的信噪比均为10dB，长度为，信号频率分别为和，初始相位，取，取不同的数值：0.3，0.25。为采样率。（1）分别用Levinson递推法和Burg法进行功率谱估计，并分析改变数据长度、模型阶数对谱估计结果的影响。（2）当正弦信号相位、频率、信噪比改变后，上述谱估计的结果有何变化？并作分析说明。2.原理分析2.1现代谱估计中的参数建模根据参数模型来描述随机信号的方法，我们可以知道，如果能确定

2、信号的信号模型，根据信号观测数据求出模型参数，系统函数用表示，模型输入白噪声，其方差为，信号的功率谱用下式求出：按照这种求功率谱的思路，功率谱估计可分为三个步骤：（1）选择合适的信号模型；（2）根据有限的观测数据，或者它的有限个自相关函数的估计值，估计模型的参数；（1）计算墨香的输出功率谱。其中以（1）、（2）两步最为关键。按照模型的不同，谱估计的方法有许多种，它们共同的特点是对信号观测区以外的数据不假设为0，而先根据信号观测数据估计模型参数，按照求模型输出功率的方法估计信号功率谱，回避了数据观测区以

3、外的数据假设问题。下面分析AR谱估计的两种方法：自相关法——列文森（Levenson）递推法和伯格（Burg）递推法。这两种方法均为已知信号观测数据，估计功率谱，两者共同特点是由信号观测数据求模型系数时采用信号预测误差最小的原则。对于长记录数据，这些方法的估计质量是相似的，但对于短记录数据，不同方法之间存在差别。2.2自相关法——列文森（Levenson）递推法自相关法的出发点是选择AR模型参数使预测误差功率最小，预测误差功率为假设信号的数据区在范围，有个预测系数，个数据经过冲激响应为的滤波器，输出预

4、测误差的长度为，因此应用下式计算：的长度长于数据的长度，上式中数据的两端需补充零点，相当于对无穷长的信号加窗处理，得到长度为N的数据。上式对系数的实部和虚部求微分使预测误差功率最小，得到（1）式中自相关函数采用有偏自相关估计，即对比上式，可知式（1）即为已推导出的Yule-Walker方程，因此自相关法也是基于解Yule-Walker方程的一种方法。但是直接解该方程，需要计算逆矩阵，不方便，因此，基于Yule-Walker方程中自相关矩阵的性质，导出Levinson-Durbin递推法，这是一种高效的

5、解方程的方法。Levinson-Durbin算法首先由一阶AR模型开始：一阶AR模型的Yule-Walker方程为由该方程解出然后令，以此类推，可以得到一般递推公式如下：称为反射系数，。，随着阶数增加，预测误差功率将减少或不变。由k=1开始递推，递推到k=p，依次得到各阶模型参数，AR模型的各个系数及模型输入白噪声方差求出后，信号功率谱用下式计算这种方法计算简单，但需要预先估计出信号自相关函数，实际中只能按照信号的有限个观测数据估计自相关函数。当观测数据长度较短时，估计误差较大，会出现谱峰频率偏移和谱

6、线分裂（在信号谱峰附近产生虚假谱线）；如数据很长，估计自相关函数较准确，但计算量大，应适当选择数据长度。2.3伯格（Burg）递推法Levinson-Durbin递推法需要由观测数据估计自相关函数，这是它的缺点。而伯格递推法则由信号观测数据直接计算AR模型参数。伯格递推法利用Levinson-Durbin递推公式，导出前向预测误差与后向预测误差，并按照使它们最小的原则求出，从而实现不用估计自相关函数，直接用观测数据得出结果。Burg递推法思想：借助格型预测误差滤波器，求前向、后向预测误差平均功率，选择

7、使其最小，求出。之后，再利用Levinson-Durbin递推法求模型参数和输入噪声方差。设信号的观测数据区间：,前向、后向预测误差功率分别用和表示，预测误差平均功率用表示，公式分别为前向、后向观测误差公式分别为上式中，信号项的自变量最大的是n,最小的是n-p，为了保证计算范围不超出给定的数据范围，在和计算公式中，选择求和范围为：。为求预测误差平均功率最小时的反射系数，令,将前、后向预测误差的递推公式代入得Burg递推法求AR模型参数的递推公式总结如下：(1)(2)(3)(4)(5)(6)(7)3.编

8、程思想（1）编写程序产生题目要求的信号和噪声（2）然后分别用两种方法的递推流程进行谱估计（3）改变题目中要求的变量参数，分析结果的变化4.代码Levensonclc;clearall;fs=100;%采样频率Ts=1/fs;N=2^7;%数据长度p1=20;%阶数f1=0.2*fs;f2=0.25*fs;%设置信号频率pha1=0;pha2=0;%初始相位SNR=2;%设置信噪比%产生信号w=randn(1,N);Am=sqrt(2*10^(SNR/10