1. 引言
通过侦收敌方辐射源发射的电磁波信号,电子侦察设备能够检测、识别目标,判定来波方向,并能经过无源定位解算,给出目标的位置。由于其不发射能量,隐蔽性强,相比有源探测设备,作战距离远,并且在一定程度上能够侦测雷达隐身目标,因此电子侦察设备一直是现代战争中的重要武器装备 [1] [2]。其中对辐射源的无源定位能够使得指挥员更好的了解平台周围的电磁态势,同时其引导的火力打击可以做到出其不意,有效的压制敌方有源探测器,保护我方平台。目前,无源定位技术仍然是各国持续发展研究的热点,其对应的武器平台也是竞相生产、装备。
现有的无源定位方法首先可分为单站和多站定位。多站定位涉及到各站之间的同步和数据传输问题,因此其必然涉及到复杂的传输设备和交互流程,且对于窄波束雷达,分布部署的侦收站难以同时接收到辐射源信号,存在共视的问题 [3]。因此,通过单站机动,对固定辐射源进行定位,更符合作战需求 [4] [5] [6] [7]。例如,载机平台经过飞行,对地面雷达站进行无源定位。单站无源定位方法可通过对目标的角度进行交叉定位。该方法的定位精度与平台的机动距离相关。例如,为了达到一定的精度,载机平台必须经过长时间飞行,才能给出定位估计值。为了提高定位的实时性,孙仲康等提出了基于相位差变化率的定位法 [1],刘文跃等给出了相应的计算相位差变化率的一些算法 [8]。该类方法存在相位差变化率模糊跳变的问题,且涉及到高度非线性的运算,不能直接获得目标位置的解析解。传统的扩展卡尔曼滤波(EKF)无法给出较好的结果。无迹卡尔曼滤波(UKF)的性能相当于二阶高斯滤波器,它的滤波精度、稳定性要优于EKF。另外也有利用粒子滤波、鸟群算法等优化算法来解决高非线性运算的问题,但此类方法计算量较大 [9] [10]。
根据以上分析,本文采用基于UKF的滤波算法,并对定位过程进行了适当的优化。首先在定位平台短距离机动之后,即对辐射源目标进行三角交叉定位,快速给出初始定位点。然后以该值作为UKF滤波的起始点,进行滤波更新,直至收敛至真实值附近。初始定位点的引入,可有效加快滤波收敛的速度,抑制发散。同时,对于二维平面的定位运算,UKF滤波只需进行2阶矩阵的运算,对硬件的算力要求较低,适合于对体积、重量、功耗敏感的平台。
2. 单站无源定位方法
2.1. 定位模型
设定场景,假定一架飞行中的直升机,对地面雷达站进行无源定位。在飞行航路中,直升机的电子侦察系统对地面雷达站持续进行角度测量。由于直升机飞行高度相对目标的距离较小,可等效为二维平面上进行运算。如图1所示。令每次测量的角度值为θk,相应载机的位置(xk, yk),则由此估计目标的位置(xT, yT)。

Figure 1. Geometric diagram of passive positioning of single station
图1. 单站无源定位的几何示意图
2.2. 三角交叉定位
为了快速给出目标位置估计值,在传感器移动一段距离之后,即采用三角交叉定位法,对目标进行定位。若两次观测时,传感器位置为(xk-1, yk-1)、(xk, yk),所测量得到的来波到达角度为θk-1、θk,两条方向的射线可以交于一点,该点即为目标的位置估计。
根据角度定义可以得到:
(1)
写成矩阵形式为:
(2)
其中,
(3)
可以得到目标位置的解析解:
(4)
2.3. UKF滤波
在进行交叉定位时,因定位基线短,初始定位值精度较差。根据应用背景,采用基于动态系统模型的三角定位方法,每次获得新的观测值之后,不断更新的带目标位置的估计。因为假定目标是固定的,所以这里的状态转移矩阵为单位矩阵。
(5)
这种情况下,状态转移矩阵是线性的,方位角θk为观测量,其观测模型可表示为
(6)
其中,ηk表示观测噪声,这里假定为零均值加性高斯白噪声,并且误差矩阵为R。由于观测模型的高度非线性,采用扩展卡尔曼滤波不能满足稳定要求,极易发散。在此我们采用UKF无迹卡尔曼滤波进行处理。
无迹卡尔曼滤波不同于扩展卡尔曼滤波,它是概率密度分布的近似,由于没有将高阶项忽略,所以在求解非线性时精度较高。UT变换的核心思想为近似一种概率分布比近似任意一个非线性函数或非线性变换要容易。其基本原理如下:
假设n维随机向量X,其均值为
,误差协方差为PX。X通过非线性函数
变换后得到随机变量Y。通过UT变换可以比较高的精度和较低的计算复杂度求得Y的均值
和方差PY。
UT的具体过程如下:
1) 计算2n + 1个Sigma点及其权值:
(7)
其中,
;α决定Sigma点的散布程度,通常取一小的正值;
通常取0;β用来描述X的分布信息,高斯情况下,β的最优值为2。
2) 计算Sigma点通过非线性函数f()的结果:
(8)
从而得知
(9)
由于X的均值和方差都精确到二阶,计算得到Y的均值和方差也精确到二阶,比线性化模型精度更高。
2.4. 算法流程
整个单站无源定位算法实现步骤如下:
1) 启动定位功能;
2) 平台移动一定距离后,进行三角交叉定位,获取初始定位点;
3) 以初始定位点作为UKF滤波起始估计值;
4) 由状态滤波估计值产生Sigma点;
5) 根据滤波方程,迭代更新目标状态估计,即目标位置估计值。
3. 仿真分析
设定仿真场景,验证方法有效性。假定载机初始位置在原点,沿Y轴向上飞行。分析在不同约束条件下(目标初始距离R,初始方位角θ,测角精度σ,滤波周期T,飞行速度v等)的定位精度。仿真持续时间为300s,每次场景仿真100次。仿真结果见图2。
由仿真结果可得,在滤波前期,因为滤波器尚未稳定,定位精度较差,且规律不明显。随着观测时间的推移,目标估计位置向真实值收敛,定位误差变小,且大致规律为载机速度高、距离近、初始角度大、测角精度高,则定位精度高。
(a) R = 30 km, σ = 1˚, T = 0.1 s, v = 50 m/s
(b) θ = 60˚, σ = 1˚, T = 0.1 s, v = 50 m/s
(c) R = 30 km, θ = 60˚, T = 0.1 s, v = 50 m/s
(d) R = 20 km, θ = 60˚, σ = 1˚, T = 0.1 s
Figure 2. Analysis diagram of simulation
图2. 仿真结果分析图
为验证在不同约束条件下的CEP定位精度,设置如下仿真实验,获取180s时刻点的定位位置。目标估计定位位置如下图所示。其中的圆半径大小为10%R,圆心为目标真实位置。孤立点为各个估计定位位置。仿真100次。
约束条件为R = 40 km, θ = 60˚, T = 0.1 s, v = 50 m/s, σ = 1˚,96%的估计点落入圆内,满足大于50%的要求。结果见图3。

Figure 3. Estimated location of single-station target (1)
图3. 单站目标估计定位位置(一)
约束条件为R = 40 km, θ = 60˚, T = 0.1 s, v = 50 m/s, σ = 1.5˚,94%的估计点落入圆内,满足大于50%的要求。结果见图4。

Figure 4. Estimated location of single-station target (2)
图4. 单站目标估计定位位置(二)
约束条件为R = 40 km, θ = 60˚, T = 0.1 s, v = 50 m/s, σ = 2˚,93%的估计点落入圆内,满足大于50%的要求。结果见图5。

Figure 5. Estimated location of single-station target (3)
图5. 单站目标估计定位位置(三)
根据以上分析,对40 km外目标,在收敛时间不大于180s时,单站定位精度优于10%R (CEP)。
4. 结论
单站无源定位是电子侦察设备的一个重要功能。然而,对于直升机等载机平台,由于其可携行重量和体积的限制,电子侦察设备需极力做到简化,且从作战需求的角度考虑,其处理响应时间也必须短。如此,单站无源定位必须在较少的输入信息下,快速给出目标的初始定位信息,然后再进行滤波优化逐步收敛。本文方法紧贴实际需求,具有较强的工程应用价值。
参考文献