[发明专利]一种适用于复杂噪声背景的谐波信号频率估计方法

说明书

技术领域

本发明涉及信号处理领域，尤其涉及一种适用于复杂噪声背景的谐波信号频率估计方法。

背景技术

谐波信号的参数估计问题在多个信号处理领域有着广泛的应用，主要是从被噪声污染的观测信号中估计谐波信号的频率。依据噪声污染的不同情况一般可分为加性噪声和复杂背景噪声(又称为乘性和加性噪声)两种情况。

目前，复杂噪声背景中谐波信号频率的估计方法主要有循环统计量方法(李宏伟，程乾生.“乘性和加性噪声中谐波恢复的循环统计量方法”，电子学报，第26卷，第7期，1998年)和广义协方差矩阵方法(杨世永.“基于广义协方差矩阵的乘性和加性噪声中的谐波恢复”，信号处理，第28卷，第2期，2012年)。循环统计量方法是基于循环统计量，利用快速傅里叶变换和峰值搜索方法实现，由于受瑞利限的影响，循环统计量方法的估计精度和频率分辨率不高。广义协方差矩阵方法利用子空间旋转不变技术估计谐波信号的频率，但其实现过程复杂，频率估计的分辨率低。因此，有必要对上述技术问题予以解决。

发明内容

本发明的目的在于克服现有技术的不足，适应现实需要，提供一种计算精度高、可实现超分辨率估计的适用于复杂背景噪声的谐波信号频率估计方法。

为了实现本发明的目的，本发明所采用的技术方案为：

设计一种适用于复杂噪声背景的谐波信号频率估计方法，它包括以下步骤：

1)计算循环协方差c(k)；

2)构造两个循环协方差矩阵A和矩阵B；

3)对矩阵A进行特征值分解，将所得特征值的最小值记做λ_min；

4)构造F和G，其中F＝A-λ_minI，G＝B-λ_minZ；

5)对矩阵F进行奇异值分解，将所得左奇异向量矩阵记做U，奇异值矩阵记做Σ，右奇异值向量矩阵记做V；

6)将步骤5)中的左奇异向量矩阵U的前P列组成一个新的矩阵记做R，将奇异值矩阵Σ的前P行和前P列组成一个新的矩阵记做S，将右奇异值向量矩阵V的前P列组成一个新的矩阵记做T；

7)构造矩阵W＝R^HGT，其中(.)^H表示共轭转置运算；

8)构造矩阵束{S，W}，并对矩阵束{S，W}进行广义特征值分解；

9)计算频率估计值。

进一步的，所述计算循环协方差c(k)的方法为：设谐波信号的N个数据测量值为x(1),x(2),…，x(N),P为谐波分量个数,对于一个取值范围在[P+1,N/2]内的整数K，计算循环协方差c(k)，其中k＝0,1,2,…,K，则：

其中(.)^*表示取共轭运算；

进一步的，所述构造两个循环协方差矩阵A和矩阵B方法为：利用循环协方差c(0),c(1),c(2),…,c(K)构造两个K×K的矩阵A和B，即：

进一步的，在所述步骤4)中，其中：

进一步的，对矩阵束{S，W}进行广义特征值分解后，将所得P个广义特征值记做ω₁，ω₂，…，ω_P。

进一步的，所述计算频率估计值的方法为：设为谐波信号频率的估计值，则：k＝1，2，…，P…，其中∠表示取幅角运算。

本发明的有益效果在于：

本发明的方法可以实现对复杂背景噪声谐波信号进行超分辨率的频率估计，且本方法谐波信号频率估计的精度高，解决了现有方法实现过程复杂，频率估计分辨率低的技术问题。

附图说明

图1为本发明方法的流程示意图；

具体实施方式

下面结合附图和实施例对本发明进一步说明：

实施例1：一种适用于复杂噪声背景的谐波信号频率估计方法，参见图1，本方法它包括以下步骤予以实现：

1)计算循环协方差c(k)；设谐波信号的N个数据测量值为x(1),x(2),…，x(N),P为谐波分量个数,对于一个取值范围在[P+1,N/2]内的整数K，计算循环协方差c(k)，其中k＝0,1,2,…,K，则：

其中(.)^*表示取共轭运算。

2)构造两个循环协方差矩阵A和矩阵B；利用循环协方差c(0),c(1),c(2),…,c(K)构造两个K×K的矩阵A和B，即：