一种基于频谱匹配的短样本低频正弦信号参数估计方法技术

本发明专利技术提供了一种基于频谱匹配的短样本低频正弦信号参数估计方法，属于信号处理技术领域。所述方法主要针对信号波形周期数小于0.8情况下的单频实正弦信号，首先计算信号线性预测自相关函数并预设信号频率估计初值；然后基于最小二乘准则构造误差函数并通过求导使误差最小化从而对信号频率估计初值进行误差校正；最后通过频谱匹配和二分法进行迭代频率估计误差校正，经一定迭代次数后求得信号频率、幅度、初相位估计值。该方法能够在波形周期数取值很小的情况下实现单频实正弦信号的全参数估计，参数估计精度明显高于最大似然方法，尤其在噪声干扰非常严重的低信噪比情况下，本发明专利技术方法表现出很好的估计性能。因此，本发明专利技术的提出和应用可以有效解决短样本低频正弦信号因波形周期数不足所导致的参数估计难题。

【技术实现步骤摘要】
一种基于频谱匹配的短样本低频正弦信号参数估计方法
本专利技术涉及信号处理
，具体涉及含噪实正弦信号的参数估计技术。
技术介绍
正弦信号参数估计是信号处理中的基础课题，短样本低频正弦信号即为信号频率偏低靠近零频且信号样本数量极少的正弦信号，常见于地震勘探、电磁波探测、工频信号检测等应用场景中。此类信号特征直观体现为信号波形周期数(CyclesinRecord,CiR)不足，由于信号样本包含的信息量有限，参数估计的难度大幅增加，许多传统算法估计精度显著下降甚至完全失效。针对低频正弦信号的参数估计问题，中国专利文献CN109782074A“一种低频正弦波快速响应全参数估计方法及装置”提出了以连续的4个采样点进行全参数估计的方法，具有响应速度快、参数少、应用要求低等优点，但该方法主要针对工程应用，由于计算样本点数少估计精度偏低，且并未考虑样本总数很少时CiR对参数估计的影响。对于CiR取值影响的问题，陈奎孚等人(陈奎孚,王建立,张森文.低频成分的频谱校正[J].振动工程学报,2008,21(1):38-42.)根据CiR取值大小指出信号频率高低是一个相对概念，但其提出的频谱校正法仍然主要针对CiR大于1的情况。文献(K.H.Li,H.W.Fung,K.C.Teh,etal.Parameterestimationofarealsingletonefromshortdatarecords[J].SignalProcessing,2004,84(3):601-617.)基于两跟或多跟傅里叶变换谱线对短样本正弦信号进行了准确的参数估计，但该方法同样难以适用于CiR小于1的情况。进一步，黄清等人(黄清,谈振辉.采样信号周期数对克拉美-罗下界的影响[J].铁道学报,2010,32(1):114-117.)通过理论分析讨论了CiR取值对参数估计性能限的影响，指出当CiR小于0.8时克拉美罗上界会随CiR减小急剧上升且对初相位非常敏感，此时无论采用何种方法均无法获得较好的估计性能。为了解决CiR小于1时的参数估计难题。闫格等人(闫格,黄翔东,刘开华.解析全相位短区间正弦波频率估计算法研究[J].信号处理,2012,28(11):1558-1564.)基于全相位傅里叶变换提出了一种通过计算相位差估计信号频率的方法，该方法具有一定工程应用价值，但在噪声影响下精度仍然较低且易受相位模糊影响。综上所述，由于短样本低频信号具有低CiR的特点,其参数估计问题仍是目前信号处理中面临的技术难题，急需有针对性地研究解决。
技术实现思路
为解决上述参数估计难题，本专利技术提出一种基于频谱匹配的短样本低频正弦信号参数估计方法。该方法可以在CiR小于0.8情况下实现单频实正弦信号的频率、幅度和初相位估计，具有比最大似然方法更好的估计性能。在CiR大于0.3小于0.8范围内，该算法具有较高的频率和幅度估计精度，尤其在噪声干扰严重的低信噪比情况下精度优势非常明显。本专利技术的应用可以有效解决低CiR情况下的正弦信号参数估计问题，为短样本低频信号相关应用场景提供可靠的参数估计结果。本专利技术的技术方案是，一种基于频谱匹配的短样本低频正弦信号参数估计方法，总体技术方案包括：首先利用含噪信号样本序列x(n)的线性预测自相关函数R(k)和预设频率估计初值ω0计算信号频率粗估计值ω1；然后将ω1作为迭代初始值代入迭代频谱匹配运算中，重复更新第i次迭代对应的信号复振幅Fi和频率估计值ωi+1；最后经过总共Q次迭代，通过ωQ+1和FQ获得到信号频率、幅度和初相位估计结果。所述频率粗估计值ω1的计算步骤具体包括：(1)记p为线性预测自相关函数R(k)序列长度，R(k)的计算式为：(2)代入R(k)和ω0计算ω1：所述第i次迭代对应的信号复振幅Fi的计算步骤具体包括：(1)优选地，令频移分量ωc＝ωi，计算频移信号(2)优选地，记整数q＝round(Nωi/π)，同时定义Fi计算式为：所述第i次迭代对应的频率估计值ωi+1的计算步骤具体包括：(1)计算两处细化频谱S+φ和S-φ，其中φ＝π/N：(2)利用二分法思想计算ωi+1：所述信号频率估、幅度和初相位估计结果和的具体计算式分别为和本专利技术提出的方法可以在CiR小于0.8情况下实现单频实正弦信号的全参数估计，在噪声干扰比较严重情况下(信噪比低于15dB)本方法于传统算法相比精度优势非常明显，此时由于噪声影响严重，多次迭代带来的性能提升效果并不明显，迭代次数Q通常可设置为3次至5次，对于高信噪比情况则可以通过大量的迭代获得更高的估计精度，但同时也将增加该方法的计算复杂度。附图说明为清楚地说明本专利技术技术方案，下面给出实施例描述所需的附图和简要说明。显然所示附图仅是本专利技术的一些实施例，本领域技术人员在不付出创造性劳动前提下还可根据所示附图创造其他类似附图。本专利技术所示附图有：图1本专利技术信号处理流程图；图2本专利技术信号的谱间干扰特征；图3本专利技术与其他方法在不同CiR取值下的频率估计方差；图4本专利技术与其他方法在不同CiR取值下的幅度估计方差；图5本专利技术与其他方法在不同CiR取值下的初相位估计方差；图6本专利技术在不同信噪比和迭代次数下的性能变化。具体实施方式以下结合附图和实施例进一步说明本专利技术的详细
技术实现思路
，本实施例以本专利技术技术方案为前提给出了详细的实施方式和计算流程，但本专利技术的保护范围不限于下述实施例，应理解实施例仅为了说明本专利技术，而不是为了限制本专利技术的保护范围。根据图1所示信号处理流程，考虑输入采样信号x(n)＝acos(ωn+θ)+z(n)，其中a、ω和θ分别为信号幅度、频率和初相位；z(n)为零均值方差为σ2的加性高斯白噪声；信号样本数量为N，当N取值较小且ω取值靠近零时即为所述短样本低频信号，此时由信号波形周期数定义式CiR＝ωN/(2π)可知，CiR值很小是所述短样本低频信号的基本特征。计算信号样本的线性预测自相关函数，由于z(n)的期望值为0，可近似地认为其中同时优选地可令p＝round(0.48N)。那么，对cos(kω)在预设频率ω0处进行一阶泰勒级数展开后即得到T(k)＝cos(kω0)-ksin(kω0)Δω0。由此即可根据最小二乘准则即可构造误差函数再由误差最小化准则，通过对误差函数求导即可求得Δω0的估计值从而对ω0进行校正得到粗频率估计值ω1：以ω1为初始值进入迭代，迭代计算的目的在于逐步消除图2所示谱间干扰，即通过迭代计算使得频移分量ωc逐步逼近ω，从而近似达到图2(b)中负频率分量仅在零频处有非零值的效果。对于第i次迭代过程：首先，优选地，令频移分量ωc＝ωi，计算频移信号此时可近似认为xm(n)的负频率谱仅在零频处有非零值基于该值可进一步求解第i次迭代的信号复振幅Fi，即利用xm(n)在零频处的谱Xm(0)值减去其正频率分量在零频处谱值后除以N后得到，其中可由两点和本文档来自技高网...

【技术保护点】
1.一种基于频谱匹配的短样本低频正弦信号参数估计方法，其特征在于，所述方法步骤为：①根据信号频率大致范围合理设定频率估计初始值ω

【技术特征摘要】
1.一种基于频谱匹配的短样本低频正弦信号参数估计方法，其特征在于，所述方法步骤为：①根据信号频率大致范围合理设定频率估计初始值ω0；②计算信号的线性预测自相关函数R(k)并求得其在ω0处的一阶泰勒展开式T(k)；③基于最小二乘准则使得R(k)与T(k)误差最小化，计算ω0的估计误差并进行误差校正后得到粗频率估计值ω1；④基于频谱匹配思想进行迭代计算，对于第i次迭代：首先基于前次迭代得到的频率估计值ωi通过频谱匹配求解信号复振幅Fi，然后根据傅里叶变换定义计算φ＝π/N时的两处细化频谱S+φ和S-φ，并通过二分法进行误差校正得到第i次迭代对应的频率估计值ωi+1；⑤经过共Q次迭代运算，利用第Q次迭代计算结果ωQ+1和FQ获得信号频率、幅度和初相位估计值。


2.根据权利要求1所述方法，其特征在于，所述步骤④中频谱匹配思想具体内...

【专利技术属性】
技术研发人员：王魁，涂亚庆，郭凌，闫华，苏永东，
申请(专利权)人：中国人民解放军陆军勤务学院，
类型：发明
国别省市：重庆;50