有代码功率谱估计Levinson递推法Burg递推法随机信号处理学案.docVIP

功率谱估计 随机信号处理 学号： 姓名： 实验三 功率谱估计 1实验内容 信号为两个正弦信号加高斯白噪声，各正弦信号的信噪比均为10dB，长度为N，信号频率分别为和，初始相位，取，取不同数值：0.3，0.25。为采样频率。分别用Levinson递推法和Burg法进行功率谱估计，并分析改变数据长度、模型阶数对谱估计结果的影响。 2实验原理 2.1 Levinson递推法： 自相关法——列文森（Lenvison）递推法是已知信号观测数据，估计功率谱。它的出发点是选择AR模型参数使预测误差功率最小。 假设信号的数据区在范围，有P个预测系数，N个数据经过冲激响应为的滤波器，输出预测误差的长度为，因此有预测误差功率为 的长度长于数据的长度，上式中数据在以外补充零点，相当于对无穷长的信号加窗处理，会引入误差。 上式对系数的实部和虚部求微分使预测误差功率最小，得 （Yule-Walker方程） 式中自相关函数采用有偏自相关估计，即 Levinson-Durbin算法：使一种按阶次递推的算法。它以和模型参数作为初始条件，计算模型参数；再用模型参数计算模型参数，k阶模型参数由k-1阶模型参数计算得到。一直计算出模型参数为止。 一阶AR模型的Yule-Walker方程为 由该方程解出 然后令，以此类推，可以