本发明属于目标检测跟踪领域,涉及一种基于信息论的预测快速协方差交互融合无源协同定位方法(information-theoreticcriterionbasedpredictedfastcovarianceintersectionalgorithminpassivecoherentlocation,it-pfci-pcl)。
背景技术:
无源协同定位(passivecoherentlocation,pcl)指的是雷达本身不发射电磁波,借助非合作外辐射源(如手机通信基站,数字电视信号基站等)发射的电磁波来检测跟踪目标。与传统的有源雷达相比,pcl系统部署灵活,低成本,抗干扰能力强,具有较强的生存能力。除此外,pcl系统利用双、多基站的空间分布性大幅提高了系统对低空和隐身目标的探测性能,受到了国内外学者的广泛关注。
与此同时,为了更加精确的实现跟踪定位,基于无源协同定位系统的航迹融合技术显得尤为重要。因此,寻找高效、实时可实现的基于无源协同定位系统的航迹融合算法具有较高的技术价值。集中式融合(centralizedtrackfusion,ctf)计算复杂度高,对通信带宽的要求高。而信息矩阵融合(informationmatrixfusion,imf)方法的跟踪精度有待提高。为了解决pcl系统中的航迹融合问题,本发明针对局部估计误差之间的互协方差未知的情况,提出了基于信息论的预测快速协方差交互融合无源协同定位方法(information-theoreticcriterionbasedpredictedfastcovarianceintersectionalgorithminpassivecoherentlocation,it-pfci-pcl)。本发明根据一种新的信息理论优化准则,得到权值有一个封闭解,并且将原先的融合状态估计和误差协方差考虑在内,通过预测协方差交互融合算法进行融合,得到目标状态。因为它不涉及非线性代价函数的优化,节省了大量的计算,而且本发明能够很好的应用于无源协同定位系统中。
技术实现要素:
本发明的目的是针对无源协同定位系统的航迹融合问题,考虑局部估计误差之间的互协方差未知的情况,提出一种基于信息论的预测快速协方差交互融合无源协同定位方法,包括以下步骤:
步骤1:构建两个接收站、一个外辐射源的无源协同定位系统的数学模型。
步骤2:两个接收站分别通过扩展卡尔曼滤波算法对目标状态进行预测更新。
步骤3:根据信息理论优化准则,将原先的融合状态估计和误差协方差考虑在内,得到各个误差协方差之间的k-l散度(kullback-leiblerdivergence),从而得到预测快速协方差交互融合算法权值的闭式解。
步骤4:通过预测快速协方差交互融合算法得到目标的状态。
在pcl系统中,集中式融合计算复杂度高,对通信带宽的要求高。而信息矩阵融合算法的跟踪精度有待提高。本发明能够在多接收站一个外辐射源的无源协同定位环境下,改善融合效果,提高目标跟踪精度,算法复杂度低。
附图说明
图1为本发明的流程图。
图2为pcl系统示意图。
图3为位置估计的均方根误差。
图4为速度估计的均方根误差。
具体实施方式
下面结合图1对本发明做进一步的分析。
具体算法步骤如下:
步骤1具体建立了二维空间中的两个接收站一个外辐射源的无源协同定位系统的数学模型。
考虑如图2所示的两个接收站和一个外辐射源的pcl系统,r1和r2分别表示接收站1和接收站2,t表示外辐射源,o表示目标,ror1表示o与r1间的距离,ror2表示o与r2间的距离,rot表示o与t间的距离,rr1t表示r1与t间的距离,rr2t表示r2与t间的距离。
首先建立目标运动模型,记目标在第k帧的状态可以为:
x(k)=fx(k-1)+v(k)(1)
式中
pcl系统中两个接收站的量测模型分别为:
其中:
式中,[xr1(k),yr1(k)],[xr2(k),yr2(k)]和[xt(k),yt(k)]分别表示r1,r2和t在第k帧的位置。量测噪声w1(k),w2(k)分别服从协方差为
步骤2具体为两个接收站分别通过扩展卡尔曼滤波算法对目标状态进行预测更新。
在扩展卡尔曼滤波算法中,以接收站r1举例说明,只要给定初始状态
1)状态一步预测方程:
其中
2)量测一步预测方程:
其中
3)新息和误差协方差预测:
其中ψ(k)表示在第k帧的新息,
4)雅可比矩阵:
其中h(k)表示在第k帧的雅可比矩阵,
5)新息协方差:
其中s(k)表示在第k帧的新息协方差,h(k)t表示雅可比矩阵h(k)的转置。
6)滤波增益:
其中w(k)表示在第k帧的滤波增益。
7)状态更新:
其中
8)误差协方差更新:
其中
接收站r2的扩展卡尔曼滤波同公式(6)-(14),最终得到两个接收站分别对目标的状态更新
步骤3具体考虑原先融合航迹的状态估计和误差协方差,根据信息理论优化准则,得到各个误差协方差之间的k-l散度(kullback-leiblerdivergence),从而得到预测快速协方差交互融合算法权值的闭式解。
由步骤2得到的两个接收站在第k帧的状态估计分别为
其中:
式中d(pa,pb)表示两个误差协方差pa,pb之间的k-l散度:
其中:tr(·)表示求秩,
而:
其中
步骤(4)、通过预测协方差交互融合算法得到目标的状态。
将融合航迹中预测融合状态估计和预测误差协方差考虑在内,得到融合状态估计
ω1(k)+ω2(k)+ωf(k|k-1)=1(27)
对本发明所提方法在无源协同定位系统中进行了仿真,如图3所示,可以看出三种方法的目标位置估计rmse都收敛,且it-pfci-pcl融合方法的目标位置估计的均方根误差rmse较信息矩阵融合与集中式融合低。如图4所示,it-pfci-pcl融合方法的目标速度估计的均方根误差rmse较信息矩阵融合与集中式融合低。即在pcl系统中,it-pfci-pcl融合方法比信息矩阵融合与集中式融合精确度更高。