[发明专利]一种二维谐波信号频率估计方法

摘要

本发明涉及信号处理领域，尤其涉及一种二维谐波信号频率估计方法，包括以下步骤：计算协方差c(s,t)；构造矩阵F；对矩阵F进行特征值分解，将所得的特征值按大小顺序记做λ1,λ2,…,λK，相应的特征值向量记做e1,e2,…,eK；构造矩阵U；计算多项式R(z)；计算第一个谐波频率分量的估计值；构造矩阵G；对矩阵G进行特征值分解，将所得的特征值按从大到小的顺序记做γ1,γ2,…,γL，相应的特征值向量记做d1,d2,…,dL；构造矩阵V；计算多项式S(z)；计算第二个谐波频率分量的估计值；构造矩阵H；对矩阵H进行特征值分解，将所得的特征值按从大到小的顺序记做η1,η2,…,ηKL，相应的特征值向量记做b1,b2,…,bKL；构造矩阵W；计算估计量Qk,l；计算二维谐波频率对的估计值。本方法可实现对二维谐波信号进行超分辨率的频率估计，具有估计精度高，计算复杂度低的特点。

主权项

一种二维谐波信号频率估计方法，其特征在于：包括以下步骤：1)计算协方差c(s,t)；2)构造矩阵F；3)对矩阵F进行特征值分解，将所得的特征值按从大到小的顺序记做λ1,λ2,…,λK，相应的特征值向量记做e1,e2,…,eK；4)构造矩阵U；5)计算多项式R(z)；6)计算第一个谐波频率分量的估计值；7)构造矩阵G；8)对矩阵G进行特征值分解，将所得的特征值按从大到小的顺序记做γ1,γ2,…,γL，相应的特征值向量记做d1,d2,…,dL；9)构造矩阵V；10)计算多项式S(z)；11)计算第二个谐波频率分量的估计值；12)构造矩阵H；13)对矩阵H进行特征值分解，将所得的特征值按从大到小的顺序记做η1,η2,…,ηKL，相应的特征值向量记做b1,b2,…,bKL；14)构造矩阵W；15)计算估计量Qk,l；16)计算二维谐波频率的估计值。