一种基于插值的多信源二维频率估计方法技术

本发明专利技术适用于信号处理技术改进领域，提供了一种基于插值的多信源二维频率估计方法，包括以下步骤：S1、对接收各路信号进行频率残差估计、幅度值估计进行初始化归零处理；S2、对接收各路信号进行傅里叶变换并进行迭代处理获得精确的参数估计；S3、根据步骤S2中获得的估计参数计算获得信号的频率估计。二维插值算法精确度高于UE算法和IMDF算法，且所用时间明显小于其他算法。实验也说明了多目标的DOD或DOA有相同角度的情况下，所提算法仍然能够高效精确运行。

【技术实现步骤摘要】
一种基于插值的多信源二维频率估计方法
本专利技术属于信号技术改进领域，尤其涉及一种基于插值的多信源二维频率估计方法。
技术介绍
2维频率估计在MIMO雷达(multipleinputmultipleoutput)的DOA(DirectorofArrival)和DOD(DirectorofDepature)联合估计以及雷达、声呐和无线通信中的2维DOA估计中都有着广泛的应用。关于这一问题已经有许多研究方法，这些方法可以被分为两类：一类是非参数估计方法，一类是参数估计方法。二维傅里叶变换就是使用最为广泛的一种非参数估计方法。尽管通过快速傅里叶变换，该方法的计算效率较高，但是由于快速傅里叶变换的网格特性，导致该方法分辨率较低。参数类方法，假设已知信号的方程形式然后通过求解闭式解得到信号方程准确的参数估计。最大似然估计方法利用信号的范德蒙结构能够成功估计出一维信号频率，但是应用在二维频率估计中时，拥有很高的计算复杂度。因此出现了一些基于子空间的估计算法，如ESPRIT(estimationofsignalparametersviarotationalinvariancetechnique)算法、MEMP(matrixenhancementandmatrixpencil)算法和PUMA(principal-singular-vectorutilizationformodalanalysis)算法，这类算法利用EVD(eigenvaluedecomposition)和SVD(singularvaluedecomposition)分解将接收到的信号分解为信号子空间和噪声子空间两部分。这类技术由于应用了SVD和EVD分解，计算复杂度较高的问题仍然未得到彻底解决。而且对于二维的频率估计，许多算法需要额外的匹配步骤这也增加了算法的复杂度，需要指出的是，当同一纬度有相同频率时，PUMA算法并不能成功匹配出信号频率。
技术实现思路
本专利技术的目的在于提供一种基于插值的多信源二维频率估计方法，旨在解决上述的技术问题。本专利技术是这样实现的，一种基于插值的多信源二维频率估计方法，所述多信源二维频率估计方法包括以下步骤：S1、对接收各路信号进行频率残差估计、幅度值估计进行初始化归零处理；S2、对接收各路信号进行傅里叶变换并进行迭代处理获得精确的参数估计；S3、根基步骤S2中获得的参数估计计算获得信号的频率估计。本专利技术的进一步技术方案是：所述步骤S2中还包括以下步骤：S21、对接收各路信号进行傅里叶变换求出二维频谱最大窗估计，其函数式：S22、利用函数计算各路信号的二维频率残差估计，其函数式为：S23、通过减掉其他路信号在k路信号所在频率点处产生的幅值余项来获得幅度值的估计，其函数式：其中，M1和M2分别是信号在一维和二维数据长度，m1＝1，2，...，M1，m2＝1，2，...，M2，x(m1，m2)接收到的二维信号模型，X(u1，u2)x(m1，m2)的频谱，e是自然指数，j是虚部，是第i次迭代获得的第k个目标信号的一维频率残差，是第i次迭代获得的第k个目标信号的二维频率残差，是第i次迭代获得的关于第k个目标信号c维计算数据，是第i次迭代获得的第k个目标信号的幅度，是第k个信号的前l个信号的第i次迭代结果，其中1＜＝l＜k，是第k个信号的后p个信号的第i-1次迭代结果，其中k＜p＜＝M，是第p个信号迭代i-1次后得到的幅度。本专利技术的进一步技术方案是：所述步骤S3中计算信号的频率估计的函数式其中，M1和M2分别是信号在一维和二维数据长度，是第k个目标信号的一维频谱最大窗估计，是第k个目标信号的一维频谱最大窗估计，是第k个目标信号的一维频率残差估计，是第k个目标信号的一维频率残差估计。本专利技术的进一步技术方案是：所述多信源二维频率估计方法中参数的估计通过不断迭代提高估计的精度，在迭代大于等于两次后具有收敛性。本专利技术的有益效果是：二维插值算法精确度高于UE算法和IMDF算法，且所用时间明显小于其他算法。实验也说明了多目标的DOD或DOA有相同角度的情况下，所提算法仍然能够高效精确运行。附图说明图1是本专利技术实施例提供的基于插值的多信源二维频率估计方法的流程图。图2是本专利技术实施例提供的频率估计的RMSE随迭代次数变化趋势图。图3是本专利技术实施例提供的频率估计的RMSE随SNR变化趋势图。图4是本专利技术实施例提供的有相等频率情况下RMSE随SNR变化趋势图。具体实施方式如图1所示，本专利技术提供的基于插值的多信源二维频率估计方法的流程图，其详述如下：给出接收到的二维信号模型为：其中K是信源数，M1和M2分别是信号在一维和二维数据长度，sk(m1，m2)是k路信源信号，表达式为：其中，βk是信号幅度，fk1∈(0，1)和fk2∈(0，1)是k路信源信号分别在第一维和第二维的频率。w(m1，m2)是均值为0方差为σ2/2的高斯白噪声。定义第k路信号的信噪比(SNR，SignalNoiseRate)为ρk＝|βk|/σ2。以求解第k路信源，c维频率为例，定义为第k个相关目标信号的c维频谱最大窗估计，其中c＝1，2。可以通过寻找2维频谱的峰值进行计算获得。定义δkc作为第k个相关目标信号的c维频率残差，记是δkc的估计。此外，X(k1，+0.5)和X(k2，+0.5)是k路信号分别在和的离散傅里叶系数。因此我们可以得到：其中，是第k路信号在相应c维的离散傅里叶系数。和是由于其他路信号的干扰所存在的希望去掉的余项，可以用表示为：和结合(3)到(6)，我们可以得到所期望的单路信号的离散傅里叶系数：并且定义这样我们就可以得到δkc和之间的差：从而得到δkc，c＝1，2的估计，然后得到频率估计：估计k路信号的傅里叶系数时，首先对其他路信号产生的余项进行估计，然后需要将该项减掉已获得准确的估计结果。然而其他路信号的余项是也是一个估计值，所以这个减掉的过程会引起最终结果δk1，δk2和βk的估计偏差。因此我们提出一种迭代的方法来获得精确的参数估计。在每一次迭代中，会利用前一次所有路信号相关参数的估计来求出本次迭代的估计和迭代步骤如下所示：在第i次迭代，求第k路信号的相关参数时。已获得前k-1路信号相关参数在第i次迭代后的估计以及从k到K路共K-k+1路信号相关参数在第i-1次迭代后的估计接下来，依据公式(6)，我们可以得到：将(10)带入(7)，并根据(8)，我们可以得到：一但获得和就可以开始幅度值βk的估计。可以通过减掉其他路信号在k路信号所在频率点处产生的幅值余项来获得：中，和是l和p路信号在处的离散傅里叶系数。基于插值的多信源二维频率估计算法1.初始化2.Forq＝1toQdo：Fork＝1toK，do：(1)如果q＝l，找到最大窗的大致估计和(2)利用(3)到(11)得到残差估计和(3)利用(12)得到3.最终利用(9)得到从前面的分析可知，所提算法中，参数的估计需要通过不断迭代以提高估计的精确度。因此我们给出估计结果的均方根误差(rootmeansquareerror，RMSE)随迭代次数Q从0到10变化的结果图。这里需要强调的一点是，当Q＝0时，所提算法就是简单的FFT频谱峰值点的搜索算法。本次测试的参数设定为：β1＝1，β2＝0.9，f11＝0.107，f本文档来自技高网...

【技术保护点】
1.一种基于插值的多信源二维频率估计方法，其特征在于，所述多信源二维频率估计方法包括以下步骤：S1、对接收各路信号进行频率残差估计、幅度值估计进行初始化归零处理；S2、对接收各路信号进行傅里叶变换并进行迭代处理获得精确的参数估计；S3、根据步骤S2中获得的参数估计计算获得信号的频率估计。

【技术特征摘要】
1.一种基于插值的多信源二维频率估计方法，其特征在于，所述多信源二维频率估计方法包括以下步骤：S1、对接收各路信号进行频率残差估计、幅度值估计进行初始化归零处理；S2、对接收各路信号进行傅里叶变换并进行迭代处理获得精确的参数估计；S3、根据步骤S2中获得的参数估计计算获得信号的频率估计。2.根据权利要求1所述的多信源二维频率估计方法，其特征在于，所述步骤S2中还包括以下步骤：S21、对接收各路信号进行傅里叶变换求出二维频谱最大窗估计，其函数式：S22、利用函数计算各路信号的二维频率残差估计，其函数式为：S23、通过减掉其他路信号在k路信号所在频率点处产生的幅值余项来获得幅度值的估计，其函数式：其中，M1和M2分别是信号在一维和二维数据长度，m1＝1，2，...，M1，m2＝1，2，...，M2，x(m1，m2)接收到的二维信号模型，X(u1，u2)x(m1，m2)的频谱，e是自然指数，j是虚部...

【专利技术属性】
技术研发人员：孙维泽，黄磊，包为民，王佳佳，张沛昌，黄敏，
申请(专利权)人：深圳大学，
类型：发明
国别省市：广东,44