功率谱估计和频率估计.pdfVIP

实验四 功率谱估计 实验内容、步骤： 实验内容包括三个： 实验一、宽带AR 过程x (n ) 是由单位方差的高斯白噪声通过滤波器 1 H (z ) −1 −2 −2 (1−0.5z +0.5z )(1+0.5z ) a. 生成x (n ) 的N 256 个样本，取p 4 并用自相关方法来计算功率谱，画出估计的 功率谱并与真实功率谱相比。 b. 重复a 中的计算20 次，分别画出20 次的重迭结果和平均结果。评论估计的方差并 说明怎样才能提高自相关方法估计功率谱的精度； c. 分别取p 6,8,12 来重复b 中的计算，描述模型阶数增加时会出现什么结果。 d. 分别采用协方差方法、修改的协方差方法来重复b,c 中计算过程，说明对宽带AR 过 程而言，哪种方法最好。 e. 把宽带AR 过程改为下列窄带AR 过程， 1 H (z ) −1 −2 −1 −2 (1−1.585z +0.96z )(1−1.152z +0.96z ) 重复a,b,c,d 中的所有分析。 y (n ) x (n ) =+w (n ) 实验二、本实验是验证最大熵方法的功率谱估计。 对随机过程 ， 2 w(n) σ x (n ) AR (2) w 是方差为 的白高斯噪声， 是 过程，由单位方差的白噪声通过如下滤波 器所获得 1 H (z ) −1 −2 1−1.585z +0.96z a. 画出x (n ) 和y (n ) 的理论功率谱。 2 σ 0.5,1,2,5 y (n ) N 100 p 2 y (n ) b. 取 w ，取 的 个样本，采用 的MEM 方法由 来 估计x (n ) 的功率谱，看看噪声对功率谱估计的精度有多大影响。 c. 改p 5 ，再重复b 中的过程，分析所观测的结果； 2 r (k ) r (k ) =+σδ(k ) d. 由于自相关序列为 y x w ，如果在计算 MEM 功率谱前从自相关值 r (0) σ2 y ω 中减去 ，用修改后的自相关序列来估计MEM 功率谱，重复c 中的过程。会提高功 - 1 - 率谱估计精度吗？ 试验三、本试验主要验证频率估计。 3 x (n ) ∑Aejn ωi =+w(n ) 令x (n ) 是谐波过程 i 1 i ，其中w(n) 是单位方差的高斯白噪声。令 A1 4e j φ1 , A2 3e j φ2 , A3 e j φ3 φi [π,−π] ， 是在 间均匀分布的不相关随机变量，取 ω 0.4πω,