一种稀疏贝叶斯学习的直接定位方法技术

本发明专利技术提供了一种稀疏贝叶斯学习的直接定位方法，利用广义双帕累托先验分布比传统分布更强的稀疏特性，保证算法不受相关信号影响的同时，还加强了信号的稀疏性，最终表现为分辨率和定位精度的提高以及低信噪比、小快拍数情况下的稳健性。本发明专利技术将稀疏贝叶斯学习引入直接定位系统，从稀疏类方法入手，解决了子空间原理类直接算法无法定位相关辐射源的问题；同时性能上的提升解决了现有最大似然原理类直接定位算法所需采样快拍数较大、分辨率不够、未知辐射源数目时无法准确定位等问题。未知辐射源数目时无法准确定位等问题。未知辐射源数目时无法准确定位等问题。

【技术实现步骤摘要】
一种稀疏贝叶斯学习的直接定位方法


[0001]本专利技术涉及无源定位领域，具体涉及一种信号的直接定位方法。

技术介绍

[0002]无源定位系统指利用目标辐射源自身发射的无线电信号进行定位的系统，在电子侦察、雷达探测等军事领域被广泛运用。直接定位作为无源定位领域一种新的定位技术，无需中间参数估计步骤，在低信噪比下相比传统两步定位方法具有估计精度高、参数自动关联等优点。该技术从2004年被提出之后，先后出现了最大似然原理类算法和子空间类算法：最大似然原理类算法具有定位性能接近性能界克拉美罗下界的优势，子空间原理类算法具有较高的定位分辨率，两类方法均被广泛研究。
[0003]然而，随着分布式集群作战概念的提出，无人机蜂群等密集目标辐射源场景对定位系统的精度及分辨力提出了更为严苛的要求。目前的直接定位算法仍存在如下不足：(1)最大似然原理类算法计算复杂度高、所需采样快拍数较大、分辨率低、未知辐射源数目时无法准确定位；(2)子空间原理类算法在辐射源之间存在互相关情况时定位效果不佳，甚至失效。

技术实现思路

[0004]为了克服现有技术的不足，本专利技术提供一种稀疏贝叶斯学习的直接定位方法。基于稀疏贝叶斯学习的参数估计方法因具备超分辨、高精度、强稳健等优越性能而被广泛应用于空间谱估计中。为了解决现有直接定位技术存在的上述不足，本专利技术以直接定位技术为核心，提出了一种稀疏贝叶斯学习的直接定位方法，利用广义双帕累托(GDP)先验分布比传统分布更强的稀疏特性，保证算法不受相关信号影响的同时，还加强了信号的稀疏性，最终表现为分辨率和定位精度的提高以及低信噪比、小快拍数情况下的稳健性。
[0005]本专利技术解决其技术问题所采用的技术方案包括以下步骤：
[0006]步骤1.构建直接定位接收数据的空间稀疏模型；
[0007](1)阵列接收数据模型建立：二维平面内有N个接收站，在接收站的远场处有Q个固定的信号辐射源，每个信号辐射源的位置为p
q
＝[x
q
,y
q
]T
,q＝1,2,
…
,Q，p
q
为第q个辐射源的位置，x
q
为辐射源横坐标，y
q
为辐射源纵坐标，各辐射源辐射的信号为互相关的窄带信号，N个接收站接收辐射源辐射的信号并将数据传送到中心处理站实现对辐射源位置的估计；
[0008]第n个接收站接收到的K个快拍的数据可建模表示为：
[0009]r
n
＝A
n
s
n
+w
n
ꢀꢀ
(1)
[0010]每个接收站阵元数目为M，采样快拍数为K，为接收站接收到的K快拍数据，为阵列流型矩阵，a
n
(p
q
)为第q个辐射源到达第n个观测站时的导向矢量，主要由目标辐射源与接收站之间的相对位置关系决定，为接收站接收到的信号包络，为第n个接收站的加性噪声；
[0011]在直接定位问题中，辐射源数目Q往往未知，定位目标就是直接利用观测到的数据r
n
定位出Q个辐射源的空间位置，n＝1,2,
…
,N；
[0012](2)构建接收数据的空间稀疏模型：为了将源定位问题转化为一个稀疏表示的问题，受稀疏重构基础理论的启发，本专利技术考虑将接收数据建立为空域稀疏的模型，将感兴趣的空间范围划分网格点，共G个网格，每个网格表示为潜在的辐射源的位置，划分的网格足够小，真实辐射源的位置在格点上；
[0013]将空间谱相应样本排列在G
×
1的向量中，因为G远大于信号源数Q，空间信号向量是稀疏的，理想情况下，它的大多数元素都接近于0，只有Q个元素与零元素有较大差异，因此采用稀疏恢复的方法来获得辐射源的位置估计；
[0014]辐射源位于其中的Q个网格点上，单采样快拍下，接收信号矢量式(1)表示为：
[0015][0016]其中，为已知的稀疏表示的过完备阵列流型矩阵，G为划分网格点的格点数，为格点上的位置矢量，并假设M<<G、Q<<G，因为目标辐射源仅位于其中的几个格点上，所以表示为空间稀疏信号，仅在对应的有信号源的网格点对应的位置上才有数值，其他元素均为0；
[0017]由稀疏信号的概念，可知每一列都是是具有稀疏度为Q的稀疏向量，其中仅有Q个非零元，其他元素均为零元素，其中，Q个非零元对应的是真实辐射源的位置；
[0018]步骤2.信号的GDP稀疏先验分布假设：
[0019](1)噪声的先验假设；
[0020]每个接收站每个阵元接收噪声均为方差的独立的零均值高斯分布，将w
n
概率密度函数表示为所有采样快拍下噪声的联合概率密度函数：
[0021][0022]噪声方差未知，为简化分析，噪声服从伽马分布：
[0023]p(α0|c,d)＝Gamma(α0|c,d)
ꢀꢀ
(4)
[0024](2)接收数据先验假设；
[0025]由噪声功率的假设得到，接收数据的先验分布表示为：
[0026][0027](3)为了对空间稀疏信号建模为GDP先验分布。
[0028]建立三层先验分布来实现；首先第一层先验，假设的每一块拍数据服从独立的复高斯分布：
[0029][0030]协方差矩阵位Λ＝diag(α)，且表示每个网格点信号的功
率；
[0031]接下来第二层先验对信号的方差α中的α
g
建模，为了得到GDP先验分布，此处对α建立超参数为ξ的伽马分布，具体概率密度函数表示为：
[0032][0033]其中，ξ＝[ξ1,ξ2,
…
,ξ
G
]T
为超参数，为了得到GDP分布先验，第三层先验假设超参数ξ为伽马分布的随机变量，具备如下的概率分布密度函数：
[0034][0035]其中，超参数h为接近于0的正常数；
[0036]综合式(5)至(7)三层先验，得到依赖于α和ξ的的概率分布密度函数表达式如下式：
[0037][0038]其中，h为可以调节稀疏度的超参数。
[0039]步骤3.基于GDP先验分布的稀疏贝叶斯学习参数推导；
[0040]基于稀疏贝叶斯学习的直接定位算法主要就是根据接收到的数据矩阵和假设的概率分布依据贝叶斯规则推导出后验概率密度，通过最大化后验概率密度函数得到目标辐射源的位置；
[0041]通过前面的先验概率假设，得到后验概率密度函数为：
[0042][0043]使得公式(10)中目标函数最大化的稀疏信号中非0元的位置即为目标辐射源的位置。
[0044]所述步骤3中，为了最大化式后验概率密度函数，采用EM算法，E步骤主要计算完全似然对数的期望；M步骤使用算法最大化期望值；
[0045](1)E步骤表达为本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种稀疏贝叶斯学习的直接定位方法，其特征在于包括下述步骤：步骤1.构建直接定位接收数据的空间稀疏模型；(1)阵列接收数据模型建立：二维平面内有N个接收站，在接收站的远场处有Q个固定的信号辐射源，每个信号辐射源的位置为p
q
＝[x
q
,y
q
]
T
,q＝1,2,
…
,Q，p
q
为第q个辐射源的位置，x
q
为辐射源横坐标，y
q
为辐射源纵坐标，各辐射源辐射的信号为互相关的窄带信号，N个接收站接收辐射源辐射的信号并将数据传送到中心处理站实现对辐射源位置的估计；第n个接收站接收到的K个快拍的数据可建模表示为：r
n
＝A
n
s
n
+w
n
ꢀꢀꢀꢀ
(1)每个接收站阵元数目为M，采样快拍数为K，为接收站接收到的K快拍数据，为阵列流型矩阵，a
n
(p
q
)为第q个辐射源到达第n个观测站时的导向矢量，主要由目标辐射源与接收站之间的相对位置关系决定，为接收站接收到的信号包络，为第n个接收站的加性噪声；在直接定位问题中，辐射源数目Q往往未知，定位目标就是直接利用观测到的数据r
n
定位出Q个辐射源的空间位置，n＝1,2,
…
,N；(2)构建接收数据的空间稀疏模型：将接收数据建立为空域稀疏的模型，将感兴趣的空间范围划分网格点，共G个网格，每个网格表示为潜在的辐射源的位置，划分的网格足够小，真实辐射源的位置在格点上；将空间谱相应样本排列在G
×
1的向量中，因为G远大于信号源数Q，空间信号向量是稀疏的，理想情况下，它的大多数元素都接近于0，只有Q个元素与零元素有较大差异，因此采用稀疏恢复的方法来获得辐射源的位置估计；辐射源位于其中的Q个网格点上，单采样快拍下，接收信号矢量式(1)表示为：其中，为已知的稀疏表示的过完备阵列流型矩阵，G为划分网格点的格点数，为格点上的位置矢量，并假设M<<G、Q<<G，因为目标辐射源仅位于其中的几个格点上，所以表示为空间稀疏信号，仅在对应的有信号源的网格点对应的位置上才有数值，其他元素均为0；由稀疏信号的概念，可知每一列都是是具有稀疏度为Q的稀疏向量，其中仅有Q个非零元，其他元素均为零元素，其中，Q个非零元对应的是真实辐射源的位置；步骤2.信号的GDP稀疏先验分布假设：(1)噪声的先验假设；每个接收站每个阵元接收噪声均为方差的独立的零均值高斯分布，将w
n
概率密度函数表示为所有采样快拍下噪声的联合概率密度函数：
噪声方差未知，为简化分析，...

【专利技术属性】
技术研发人员：谢坚，王秋平，王伶，张兆林，陶明亮，宫延云，韩闯，粟嘉，
申请(专利权)人：西北工业大学，
类型：发明
国别省市：