一种强弱信号的直接定位与信源数估计方法技术

本发明专利技术公开了一种强弱信号的直接定位与信源数估计方法，首先对所有基站接收的数据进行分段傅里叶变换，估计基站接收信号的协方差矩阵；然后对得到的协方差矩阵进行特征值分解，再利用分解得到的特征向量计算特征波束的空间谱，通过比较空间谱的最大值与除最大值之外的平均值之差和所设定的门限值估计各个信号的位置及信源数目。本发明专利技术方法简单实用，在有多个强弱信号并存时，能精确估计强弱信号的位置和信源数，可用于雷达、通信、导航、测控和电子侦察众多领域中位置的定位与干扰的抑制。

【技术实现步骤摘要】
一种强弱信号的直接定位与信源数估计方法
本专利技术属于信号处理
，具体涉及一种定位与信源数估计方法。
技术介绍
直接定位(direction-Position-Determination,DPD)技术是信号处理领域的一个重要研究方向。近十年来，直接定位估计技术因其优越的性能，在雷达、声纳以及通信等领域得到了广泛应用。随着现代无线通信系统的迅速发展，对于定位的精度和分辨率需要进一步的提升。因此，针对直接定位算法在信号处理领域受到了越来越多的关注，而目前的间接定位如基于信号到达角(angleofarrival，AOA)和信号到达时间(timeofarrival，TOA)的定位技术要比一步直接定位技术损失更多的信息量，其精度也不如直接定位技术。然而，目前大多数超分辨DPD算法都是基于信号幅度大致相等的情况下进行的定位，但在实际工程应用中，由于信号存在多径衰落或者雷达探测不同RCS(radarcrosssection)大小的目标，这样接收端接收到的信号幅度相差很大，被称为强弱信号。目前存在的一些DPD算法，如最大似然(ML:MaximumLikelihood)、多重信号分类(MUSIC:MultipleSignalClassification)、最小无畸变方差(MVDR:MinimumVarianceDistortionlessResponse)等方法，由于具有超分辨的多目标直接定位性能引起很多学者的关注。但是，同时存在信号幅度相差很大的强弱信号时，由于强信号的影响，利用上述方法难以准确估计出弱信号的位置，从而导致前述超分辨DPD算法性能退化，甚至失效。针对超分辨DPD算法估计问题，目前已有的方法主要有：1.AnthonyJ.Weiss发表于2004年的IEEESIGNALPROCESSINGLETTERS上的《DirectPositionDeterminationofNarrowbandRadioFrequencyTransmitters》一文中，提出用最大似然估计的方法来估计目标的位置。2.AnthonyJ.Weiss等发表于2005年的EURASIPJournalonAppliedSignalProcessing上的《DirectPositionDeterminationofMultipleRadioSignals》一文中，提出了用MUSIC的方法来估计目标的位置。3.TomTirer等发表于2016年的IEEESIGNALPROCESSINGLETTERS上的《HighResolutionDirectPositionDeterminationofRadioFrequencySources》一文中，提出了用MVDR的方法来估计目标的位置，并且性能要优于1中2005年提出的最大似然估计方法。上述的超分辨DPD算法都是针对信号源输出功率大致相等的情况下进行的信号源直接定位，目前并没有提出针对强弱信号背景下的超分辨DPD算法，上述算法对强弱信号背景下的弱信号的位置不能正确的估计，甚至完全失效，因此也不能对信号源数目进行很好地正确估计。对弱目标的漏估与过估，将会影响后续目标信息的正确提取。
技术实现思路
为了克服现有技术的不足，本专利技术提供了一种强弱信号的直接定位与信源数估计方法，首先对所有基站接收的数据进行分段傅里叶变换，估计基站接收信号的协方差矩阵；然后对得到的协方差矩阵进行特征值分解，再利用分解得到的特征向量计算特征波束的空间谱，通过比较空间谱的最大值与除最大值之外的平均值之差和所设定的门限值估计各个信号的位置及信源数目。本专利技术方法简单实用，在有多个强弱信号并存时，能精确估计强弱信号的位置和信源数，可用于雷达、通信、导航、测控和电子侦察众多领域中位置的定位与干扰的抑制。本专利技术解决其技术问题所采用的技术方案包括以下步骤：步骤1：假设空间中具有Q个信号源，L个基站；第l个基站接收到的信号为：式中，αlq是路径衰落导致的复常数，al(pq)是第l个基站对第q个信号源的导向矢量，l＝1,...,L，q＝1,...,Q，pq为第q个信号源的位置坐标；sq(t-τl(pq))是信号经过延时τl(pq)后的结果，nl(t)代表零均值的复高斯白噪声，t为时间，T为观测周期；假定信号源和基站在观测时间内都静止，将第l个基站接收到的信号在时间段[0,T]分为J段，假设每段信号有K个频率点；对每个时间段的信号做傅里叶变换，第l个基站第j个时间段的信号在第k个频率点的傅里叶变换为：式中，fk是傅里叶变换第k个频率点，rl(k,j)、sq(k,j)、nl(k,j)分别是rl(t)、sq(t)、nl(t)在第j个时间段的傅里叶变换,j＝1,...,J,k＝1,...,K；定义得到：将L个基站的信号合并，定义：s(k,j)＝[s1(k,j),...,sQ(k,j)]T(8)p＝[p1,...,pQ](9)进一步得到：r(k,j)＝A(k,p)s(k,j)+n(k,j)(10)计算L个基站接收信号的协方差矩阵为：其中：A(k,p)＝Λk(p)Bα(12)α＝[α1,...,αL]T(15)IL代表L×L的单位矩阵，1M代表M×1的全一向量，代表克罗内克积，规定||α||＝1；步骤2：对协方差矩阵进行特征值分解得到：式中，Γ为特征值从大到小排列组成的对角矩阵，U＝[u1...uLM]为特征向量组成的矩阵，ui为分解的特征向量，每个特征向量对应一个特征波束，i＝1...LM，LM为L个基站包含的所有阵元数目；步骤3：令i等于1，定义IL为信号门限值；步骤4：用第i个特征向量计算第i个特征波束的空间谱Pi(p)：其中，为的逆矩阵，λmax{·}代表取最大特征值操作；步骤5：将第i个特征波束的空间谱Pi(p)化为单位dB的形式Pi_dB(p)，计算Pi_dB(p)的最大值Pi_dBmax(p)及对应的pimax；计算除最大值位置Pi_dBmax(p)外的其余位置的平均值若则第i个信号源的位置为pimax，令i+1，转到步骤4；否则，若则信号源数目为i-1；最终得到信号源数量和每个信号源的位置。优选地，所述IL＝13.4dB。由于采用了本专利技术的一种强弱信号的直接定位与信源数估计方法，取得了以下有益效果：1.本专利技术通过比较每个特征波束空间谱的最大值与除了最大值位置之外其他位置空间谱的平均值，能同时估计强弱信号的位置与信源数，不需要信源数已知。2.相较于现有技术，本专利技术所提方法的优势是在能够对强信号进行很好地定位同时，也能实现对于弱信号位置的精确估计，同时对于当多个强弱信号相隔很近的时候仍能实现很高的分辨率；本方法充分的利用了所有基站的联合信息，能够实现对于非信号源的位置的能量很好地抑制，抗干本文档来自技高网...

【技术保护点】
1.一种强弱信号的直接定位与信源数估计方法，其特征在于，包括以下步骤：/n步骤1：假设空间中具有Q个信号源，L个基站；/n第l个基站接收到的信号为：/n

【技术特征摘要】
1.一种强弱信号的直接定位与信源数估计方法，其特征在于，包括以下步骤：
步骤1：假设空间中具有Q个信号源，L个基站；
第l个基站接收到的信号为：



式中，αlq是路径衰落导致的复常数，al(pq)是第l个基站对第q个信号源的导向矢量，l＝1,...,L，q＝1,...,Q，pq为第q个信号源的位置坐标；sq(t-τl(pq))是信号经过延时τl(pq)后的结果，nl(t)代表零均值的复高斯白噪声，t为时间，T为观测周期；
假定信号源和基站在观测时间内都静止，将第l个基站接收到的信号在时间段[0,T]分为J段，假设每段信号有K个频率点；对每个时间段的信号做傅里叶变换，第l个基站第j个时间段的信号在第k个频率点的傅里叶变换为：



式中，fk是傅里叶变换第k个频率点，rl(k,j)、sq(k,j)、nl(k,j)分别是rl(t)、sq(t)、nl(t)在第j个时间段的傅里叶变换,j＝1,...,J,k＝1,...,K；
定义得到：



将L个基站的信号合并，定义：












s(k,j)＝[s1(k,j),...,sQ(k,j)]T(8)
p＝[p1,...,pQ](9)
进一步得到：
r(k,j)＝A(k,p)s(k,j)+n(k,j)(1...

【专利技术属性】
技术研发人员：谢坚，邓均午，况梅东，陈清浪，丁立超，陶明亮，张兆林，王伶，范一飞，粟嘉，杨欣，韩闯，宫延云，
申请(专利权)人：西北工业大学，中国船舶重工集团公司第七一五研究所，
类型：发明
国别省市：陕西;61