[发明专利]一种乘性和加性噪声中谐波信号频率估计方法

摘要

本发明公开了一种乘性和加性噪声中谐波信号频率估计方法，其包括以下步骤计算循环协方差；构造循环协方差矩阵；特征值分解；构建噪声子空间矩阵；计算空间谱；计算频率估计值。本发明的乘性和加性噪声中谐波信号频率估计方法，能提高频率估计的精度和频率分辨率，且易于实现。

主权项

一种乘性和加性噪声中谐波信号频率估计方法，其特征在于，其包括以下步骤：步骤1：计算循环协方差：所述计算循环协方差的方法为：设谐波信号的N个数据测量值为x(1),x(2),…,x(N)，P为谐波分量个数，对于一个取值范围在[P+1,N/2]内的整数M，计算循环协方差c(α)，α＝0,1,2,…,M‑1：c(α)=1M-αΣn=1M-αx2(n)(x2(n+α))*-1(M-α)2Σn=1M-αx2(n)Σn=1M-α(x2(n+α))*;]]>其中(·)*表示取共轭运算；步骤2：构造循环协方差矩阵：所述构造循环协方差矩阵C的方法为：步骤3：特征值分解；步骤4：构建噪声子空间矩阵；步骤5：计算空间谱；步骤6：计算频率估计值:所述计算频率估计值的方法为：依次找出空间谱W(ωk)的P个峰值，首先找出第一大峰值将对应的ω值记为ω1，其次找出第二大峰值将对应的ω值记为ω2，由此类推，最后找出第P大的峰值将对应的ω值记为ωP；计算谐波信号频率的估计值为：m＝1,2,…,P。