一种目标运动要素纯方位解算方法

程善政，陈双，何心怡

(中国人民解放军92578 部队，北京 100161)

0 引 言

水下无人航行器是未来海军战斗力的倍增器，可在浅水区或其他潜艇、水面舰无法进入的未知复杂水域执行水下威胁目标侦察、战术海洋学调查、监视预警、远程通信中继、攻击敌目标舰艇、诱骗等军事任务，也可应用于资源及地形探测、水下光缆等设施的维护及打捞等，在国家领海防御和海洋资源开发等方面将发挥日益重要的作用[1]。攻击型水下无人航行器具备自主作战能力，将作战模式从平台级交战变为武器级交战，将以兵器攻防代替平台攻防，可大大降低人员的战损率。

潜艇、水下无人航行器等可采用基于主动探测的方位、距离法解算目标运动要素，但测量参数较多，水下主动探测方式的隐蔽性差。纯方位解算方法是一种从探测器观察到的加噪方位信息中估计目标运动要素的方法。纯方位法的特点是仅需要测量目标的方位值，测量参数少、算法比较简单、运算量较小，但解算时间长，收敛率低，但在航海、航空等领域中也可广泛采用基于被动探测的纯方位法解算目标运动要素[2]。王艳艳等[3]考虑到AUV 机动性能，采用分段轨迹优化方法，AUV 跟踪目标的初始阶段以定位的位置误差作为优化对象，提出一种基于短期预测的轨迹优化方法。董严红等[4]分析了自导鱼雷纯方位攻击的使用过程，采用大样本仿真的方法，得出不同情况下鱼雷攻击提前角的最优解。舒象兰等[5]提出一种基于匹配场处理思想改进的纯方位目标运动分析方法。

1 纯方位解算的传统算法

1.1 符号说明

1.2 数学模型

在目标等速直航的条件下，由图1 可得如下关系：

之后，传统方法是基于式（8），使用最小二乘法解参数A，B，C，这时记

图 1 纯方位解算方法示意图Fig. 1Schematic diagram of azimuth-only solution method

其中 ε为误差，未知参数 β可从如下正规方程解出

在目标和本舰等速直航的条件下，上述方程有无数的解，因此，为了得到有意义的解潜艇在测量过程中应适当机动。

传统解算方法认为解算空间是在无限的实数空间中进行，在无限空间上寻找方程的最优解，这给解算结果带来了许多不确定性，有时会造成方程的不收敛。

实际上，鱼雷攻击的目标运动要素的解算空间是有限的，将有限的连续空间范围离散化，就可以发现目标运动要素解的个数有限，在有限个目标运动要素的估计中，寻找估计方差最小的即为目标运动要素估计。

2 一种新的纯方位解算方法

设解算向量X=(D0,Vm,Qm)∈R3，根据声呐测量特性及目标的运动特性， X应属于有限的3 位实数空间。将三维实数空间离散化，可以得到有限个目标运动要素的实数解。

将D0,Vm,Qm分别划分为n1,n2,n3个水平，水平间隔在解算误差允许的范围之内，离散化后可能的解算值应该有n1×n2×n3个，在 n1×n2×n3个可能的解算值中，找出一个解D0i,Vmi,Qmi使得拟合值 yi与实际测量值 yi之间的偏差平方和达到最小。

这种方法我们称为网格法，理论上可行，但是需要进行 n1×n2×n3次计算，计算时间长，不适合工程上实时计算。

仔细分析，在解算空间中，D0,Vm,Qm之间存在着很强的相关性，找出这种相关的关系：

由此，解算值可以划分为 n2个，其值分别是(f1(Vmi),Vmi,f2(Vmi))，在 n2个解算值中推算观测值 yi的拟合值达到最小。

现代潜艇中，被动测距声呐可以测目标距离，但误差较大，但是我们由此可以判定目标距离的最大最小范围(d1,d2)，将此范围划分成 n个水平，可以得到n个初始距离的解，在 n个可能的解算值中，找出一个解(D0k,Vmk,Qmk)使得拟合值与实际测量值 yi之间的偏差平方和达到最小。

之后，传统方法是基于式（21），使用最小二乘法解参数A,B,C，这时记

根据式（26），解出 Vmj,Qmj。

对每个采样时间ti；

这种方法解算时间是传统的最小二乘滤波算法的n倍，同时也要占用一定的存储空间，但在现代的计算机条件下这点存储量和解算时间是微不足道的。为了优化解算方法，可以在解算精度和解算时间上进行折中。

3 实 例

假设目标初始距离 D0= 90 链，目标航速 Vm=22 kn，目标初始方位Qm0=10°。

图 2 初距解算误差随时间变化比较图Fig. 2Comparison diagram of initial distance solution error with time

为了验证改进的滤波算法正确性与合理性，分别对传统的最小二乘滤波法、新的滤波算和网格法进行仿真计算，解算精度比较如图2～图4 所示。结果表明，新的滤波算法在收敛时间和收敛精度上都优于传统的最小二乘滤波法，而网格法略优于新的滤波算法，但在滤波时间上占用的CPU 时间太长，不适于实时滤波。

图 3 敌速解算误差随时间变化比较图Fig. 3Comparison diagram of initial speed solution error with time

图 4 解算敌向误差随时间变化比较图Fig. 4Comparison diagram of initial course solution error with time

影响新的滤波算法的收敛精度的主要因素之一是测初始方位 f0误差，如果 f0误差很小，收敛精度高，如果 f0误差大，影响估算的精度，造成滤波误差较大。

4 结 语

本文针对传统纯方位纯方位法解算目标运动要素时存在的解算时间长，收敛率低的问题，提出一种新的纯方位解算方法。根据案例计算结果可以看出，该滤波算法在收敛时间和收敛精度上都优于传统的最小二乘滤波法，而网格法略优于新的滤波算法，但在滤波时间上占用的CPU 时间太长，不适于实时滤波。下一步将针对复杂态势下的解算情况开展深化研究，为水下无人航行器作战使用以及武器装备射击诸元计算的提供支撑。