基于Matlab实现现代功率谱估计[1]要点

2011年 8月 15日第 34卷第 16期

现代电子技术

M odern Electro nics T echnique

A ug. 2011V ol. 34N o. 16

基于 Matlab 实现现代功率谱估计

王春兴

(山东师范大学物理与电子科学学院 , 山东济南 250014

摘要 :功率谱估计可以分为经典谱估计和现代谱估计。现代谱的估计可建立 A R 模型对离散信号进行谱估计、建立 M A 模型和 A RM A 模型进行谱估计。基于 M atlab 对三种模型进行仿真 , 并对结果进行了分析。结果显示 , 三种模型对现代谱的获得是有效的 , 并得到较好的谱估计。

关键词 :P SE; 现代功率谱估计 ; AR 模型法 ; A RM A

中图分类号 :T N911-34; G202 文献标识码 :A 文章编号 :1004-373X (2011 16-0065-03

Modern Power Spectrum Estimation Based on Matlab

W AN G Chun -x ing

(Colleg e o f Physics and Elect ro nics, Shando ng No rm al U niversity , Jinan 250014, Chi na

Abstract :Po wer spectr um estimation can be divided into classical spectr al estimat ion and modern spectr al estimation. M odern spectr al estimation model can establish AR mo del, M A mo del and ARM A model fo r discr ete sig nals to per for m spec -t ral

estimatio n. T hese t hr ee models can be simulated based o n M atlab, and the r esults ar e analy zed. T he r esult s sho w that the three models of mo der n spect rum are valid, and can get better spectrum estimatio n.

Keywords :PSE; mo der n pow er spect rum est imatio n; A R model method; A RM A

收稿日期 :2011-03-26

基金项目 :国家自然科学基金项目资助 (10874103

随机信号在时域上是无限长的 , 在测量样本上也是无穷多的 , 因此随机信号的能量是无限的 , 应用功率信号来描述。然而 , 功率信号不满足傅里叶变换的狄里赫莉绝对可积的条件 , 因此严格意义上随机信号的傅里叶变换是不存在的。因此 , 要实现随机信号的频域分析 , 不能简单从频谱的概念出发进行研究 , 而是功率谱 [1-2]。

信号的功率谱密度描述随机信号的功率在频域随频率的分布。利用给定的 N 个样本数据估计一个平稳随机信号的功率谱密度叫做谱估计。谱估计方法分为两大类 :经典谱估计和现代谱估计。作为经典谱估计 , 其主要缺陷 [3]是描述功率谱波动的数字特征方差性能较差 , 频率分辨率低 ; 方差性能差的原因是无法获得按功率谱密度定义中求均值和求极限的运算。分辨率低的原因是在周期图法中 , 假定工作区域以外的数据全为零 , 而对相关函数法中 , 假定延迟窗以外的自相关函数全为零。这是不符合实际情况的 , 因而产生了较差的频

率分辨率。而现代谱估计的目标都是旨在改善谱估计的分辨率。

1 现代功率谱估计的参数模型

如果将平稳的随机过程 x (n 看成是一个输入序列

u(n (白噪声过程 , 其方差为 2

激励线性系统 H (z 的

输出 , 那么由已知的 x (n 或其自相关函数 r(m 就可估计 H (z 的参数 , 再由 H (z 的参数估计 x (n 的功率

谱。即当输入一线性系统 H (z 时 , 其输出 x (n 的功率谱可用系统 H (z 的参数来估计。

按照该思路 , x (n 的功率谱可表示为 :P X (e j

= 2

1+

q

k=1

b k e

-j k

2

/

1+

p

k=1

a k e

-j k

2

(1

式中 : 2是激励白噪声的方差 ; P X (e j

为功率谱密度 , a k 和 b k 为模型参数。式 (1 也称为 ARM A 模型。

如果式 (1 参数 b 1, b 2, , b q 全为 0, 就变为 AR 模型 , 它是一个全极点模型 [4], 其含义是该模型现在的输出是现在的输入和过去 p 个输出的加权和。该情况下的功率谱密度为 :

P X (e j

=

2

1+

k=1

a

k

e

-j k

2

(2

如果 ARMA 模型参数 a 1, a 2, a p 全为 0, 则变为 M A 模型 , 它是一个全零点模型。

P X (e j

= 2

1+

q

k=1

b k e -j k

2

(3

AR 模型的正则方程是一组线性方程 , 而 M A 和

ARM A 模型是非线性方程。 2 基于 AR 模型的功率谱

AR 模型又称为自回归模型 , 它是一个全极点模

型 , 它当前输出是现在输入和过去输入的加权和 :

x (n =-

p

k=1

a k x (n -k +u(n

(4

式中 :u(n 为白噪声信号 ; p 为 AR 模型的阶数。对式 (4 求 z 变换 , 便得到系统函数 :