matlab双谱估计,AR双谱估计的Matlab实现

收稿日期：2012-09-03 作者简介：何莹(1974-)，女，江西崇义人，韶关学院物理与机电工程学院讲师，硕士，主要从事随机信号分析与处理方面的研究. AR 双谱估计的 Matlab 实现 何 莹 (韶关学院 物理与机电工程学院， 广东 韶关 512005) 韶关学院学报·自然科学 Journal of Shaoguan University·NaturalScience 2012 年 10 月 第 33 卷 第 10 期 Otc.2012 Vol.33 No.10 摘要：根据 AR 双谱估计的基本原理，对其 Matlab 实现方法进行研究，给出了 AR 模型双谱估计的 Matlab 实现，最后 给出一个仿真测试实例，仿真结果表明了该方法的有效性. 关键词：双谱；AR 模型；Matlab 中图分类号：TN911 文献标识码：A 文章编号：1007-5348(2012)10-0040-04 双谱估计在通信、雷达、机械学、生物医学、电磁学、地球物理学、流体动力学、故障诊断等领域具有广泛的应用［1］.双谱提供了比功率谱更好的特性［2-3］，它可以检验信号偏离高斯分布的程度，能保留信号的相位信息，并能自动抑制加性高斯噪声的影响等，是分析非高斯信号、非因果、非最小相位系统的主要工具. Matlab 是一种以矩阵运算为基础的交互式的可视化软件系统环境，集数值分析、矩阵运算、信号处理和图形显示于一体，为科学研究、工程设计及众多学科领域提供了一种简捷、高效的编程工具.用 MATLAB 对信号进行双谱估计和分析，其实现方法非常方便、简单. 1 双谱基本理论 设{x(n)}为零均值、实平稳随机信号，其三阶累积量 c3x(m1,m2)定义为： c3x(m1,m2)=E［x(n)x(n+m1)x(n+m2)］ (1)其双谱定义为三阶累积量的二维傅里叶变换，即有： Bx(ω1, ω2)= ∞ m1=-∞ Σ ∞ m2=-∞ Σc3x(m1,m2)exp［-j(ω1m1+ω2m2］ (2) 其中，|ω1|燮π, |ω2|燮π, |ω1+ω2|燮π［1］. 在实际应用中由于可供处理的观测信号一般都是具有有限长度的，因此无法精确应用式(2)求解随机过程的双谱， 而只能对双谱进行估计.AR 模型法在观测数据较短的情况下就可获得高分辨率的双谱估计，是双谱估计的常用方法之一. 设零均值采样数据{x(n)}(n=1,2,…,M)由以下 p 阶 AR 模型表示： x(n)+ p k=1 Σakx(n-k)=w(n) (3) 式中 w(n)为零均值、具有有限的非零累积量 γ3e 的高阶白噪声.对(3)式求三阶自相关，可得如下的正则方程［4］： c3x(m,n)+ p k=1 Σakc3x(m-k,n)=γ3eδ(m,n) 第 10 期 m>0, 坌n 或 (m,n)=(0,0) (4) 上式可用矩阵表示为： Ra=b (5) a=［1,a1,…,ap］T b=［γ3e, 0, …, 0］T c3x(0,0) c3x(-1, 0) … c3x(-p,0) c3x(1,0) c3x(0, 0) … c3x(-p+1,0) … … … … c3x(p,p) c3x(p-1, p) … c3x(0,p 坌坌坌坌坌坌坌坌坌坌坌坌坌坌坌 坌坌 坌坌坌坌坌坌坌坌坌坌坌坌坌坌坌 坌 坌) 矩阵 R 是把(4)式的各个方程按 m(m=0,…,p)从小到大、而在 m 相同时 n(n=0,…,m)由小到大排列而形成的. 由于 R 不是方阵，因此方程组(4)只能在最小二乘意义下用伪逆运算求解. 其中三阶累积量