基于分段建模的目标运动分析方法

王国刚, 李洪瑞

(1.海军驻716研究所军代表室， 江苏 连云港　222061;2.江苏自动化研究所, 江苏连云港　222061)

纯方位目标运动分析(TMA)是一种重要的被动跟踪定位技术,在声呐目标定位跟踪中有重要应用,在观测平台性能一定的情况下,跟踪系统的性能主要依赖TMA算法。纯方位TMA系统具有两个显著特点:一是观测平台必须机动[1-3];二是跟踪系统是非线性系统。因此跟踪结果精度低、收敛速度慢,一定程度上限制了TMA的性能,也导致了算法处理上的一些困难,不容易满足实际需求。鉴于问题在实际中的重要性,理论界和工程中提出了很多方法,在观测平台机动、改善系统的可观测性、TMA建模等方面采取了一些有效的处理方法[4-5]。尤其近年来国内外学者从引入新的信息的角度开展了相关研究,如方位-多普勒频率TMA[6]、基于浅海射线声学多途结构的单水听器水下TMA[7]、基于线谱瞬时频率估计的声呐目标定位跟踪方法[8]等,这些研究成果获得了较纯方位TMA系统具有更宽容的可解性条件、更快速的收敛性以及更高的参数估计精度,以及对近距离目标运动分析的适应性等优越性能,在实际中收到了不同程度的效果。但是总体看,TMA的性能距离不断提高的实际应用需求还存在差距,其技术进展还在不断完善和深入研究中。

考虑到实际中,通常为了对态势进行进一步判断需要观测平台按固定航向和速度运动一定的时间[9],因此本文提出在观测平台不机动时求解目标相对航向、在观测平台机动后求解初始距离和其他参数的分步建模与计算算法,并进行了可解性和病态分析,提出观测平台机动优化策略,从理论上得到可解条件和克服病态方法,并有利目标运动参数估计。该算法的特点是只要求观测平台进行机动一次,所建立的最小二乘方程的维数较低。

1　问题描述

如图1所示,假定目标作匀速直线运动,T0T″是其运动轨迹。从初始时刻t0观测平台于O(0,0)测得目标方位B0并开始对目标进行定位跟踪。观测平台在t0≤ti

为减少待求参数个数,降低模型维数,目标初始方位由测量值B0取代且认为已知,因此设目标初始距离为D0、速度向量为(VTx,VTy)T,则目标t时刻位置坐标为

(1)

2　纯方位TMA模型

设在t时刻目标的测量方位为B(t)(或Bt),那么不考虑误差的测量方程为:

(2)

其中(xO(t),yO(t))为观测平台的位置坐标:

(3)

VO(t)和CO(t)是观测平台速度和航向。

当t0≤t

(4)

2.1　相对航向计算模型

令

(5)

称Cr为目标相对航向。当t0≤t

(6)

(7)

其中

At=((t-t0)cosB(t),-(t-t0)sinB(t))

(8)

2.2　距离计算模型

当t′≤t

atD0=bt(t′≤t

(9)

其中

at=sin(Bt-B0)+(t-t0)VrysinBt-(t-t0)VrxcosBt

(10)

bt=yO(t)sinBt-xO(t)cosBt+(t-t0)VO0sin(CO0-Bt)

(11)

2.3　速度计算模型

速度的计算可以利用已经计算得到的D0及(Vrx,Vry)T,从公式(5)反演得到:

(12)

2.4　最小二乘解

根据上述航向、距离、速度计算模型,由多个量测方位Bi=B(ti)(t0≤ti

(13)

其中

A=

(14)

而由式(9)及最小二乘法得到D0满足的线性方程:

(15)

其中

(16)

3　算法分析及观测平台机动策略

3.1　解的存在性

矩阵A可逆是方程(13)有唯一解的充要条件。不难计算出矩阵A的行列式为:

(17)

可见,只要在t0≤ti

对于方程(15),为简单起见,假定在t′≤t

(18)

可见,若在t′≤t

3.2　病态分析与观测机动策略

线性方程(组)(13)与(15)的病态可能导致相对航向、初距的解算精度不可靠,应用中应当避免。方程组(13)的病态性可以由方程的系数矩阵的条件数κ(A)(即系数矩阵的最大特征值与最小特征值的比值)来度量。可以计算出矩阵A的特征方程如下:

λ2-cλ+d=0

(19)

其中

(20)

通过一些计算得:

(21)

由此,可对κ(A)进行估计如下:

(22)

即

κ(A)=O(ε)

(23)

κ(A)=O(ΔB-2)

(24)

可见条件数κ(A)与ΔB-2同阶。因此ΔB过小可能导致最小二乘问题病态。这意味着在目标初距、速度一定的条件下目标舷角较小时不利于解算目标相对航向。

方程(15)的病态问题主要源自ai过小。

综合以上论述,得到观测平台机动策略如下:

1)当观测平台速度不变、仅改变航向时,则:

(25)

2)当观测平台航向不变、仅改变速度时,由于

(26)

所以速度机动量不宜过小。

3)由于需要在t0≤t

4　仿真计算

仿真算例1:目标初距2km,初始方位10°,航向165°,速度18kn;t0、t′、t″分别取为0、8min、15min;当t0≤t≤t′时观测平台速度8kn、航向345°,当t′

仿真算例2:初始距离为12000m、t′=6min,其余参数与算例1相同,态势如图3所示。采用Monte Carlo进行100次仿真计算,并与最小二乘法进行了比较。对两个算例仿真计算进行误差统计,其误差曲线分别如图4、图5所示。

由仿真计算结果可以看出,与直接最小二乘法相比,本文提出的分段最小二乘法解算目标运动参数提高了精度、缩短了收敛时间,尤其是在跟踪时间较短时,能够较快地获得可信的目标运动参数,适应于目标跟踪初始阶段使用。当解算时间足够长时,本文算法解算的运动参数存在一定偏差,与传统算法结合使用,有望满足不同阶段对目标运动参数的使用要求。此外,由于本文算法只需要观测平台机动一次,比起一般的二次机动或复杂机动方案及其参数求解算法,更有利观测平台占位。从仿真计算例进行对比也可以看出,针对近距离目标,本文提出的分段最小二乘法解算的目标运动参数精度、收敛时间性能较好,相比远距离目标明显提高,该算法对中近距离目标适应性好。

5　结束语

目标运动分析是工程中备受重视的十分复杂的问题,在许多被动定位中尤其在纯方位定位中广泛应用,它受目标运动特性、观测平台机动性能、传感器特性等多方面因素的影响。特别是它与观测平台的机动紧密联系着,而观测平台机动是一个棘手的老大难问题,它既要考虑目标跟踪定位也要考虑攻击占位,同时还要考虑观测平台隐蔽性,机动次数不宜过频。本文提出的分阶段建模目标运动参数解算算法,由于对观测平台机动次数要求较低,只需进行一次有效机动,解算模型维数低,比较适合观测平台难以多次有效机动的应用场景。由于解算模型维数较低,计算量较一般滤波算法和批处理算法少。文中提出的目标运动参数的可解条件及克服病态的方法和观测平台机动策略,有利于提高模型的稳健性和工程适用性。本文提出的方法在跟踪时间较短时,能够较快地获得可信的目标运动参数,适应于目标跟踪初始阶段使用,具有一定的工程应用价值。