6. 辐射源定位系统的精度¶
几乎所有军事资产,从大型固定设施到单个飞机、车辆或小单位,都必须发射某种信号以完成其任务。尽管这些资产可能在夜间、雾霾或烟雾中、伪装状态下或超出视距范围内不可见,但其发射机的位置对应着其物理位置。对这些位置发射的信号进行分析通常可以识别出资产的类型(武器、军事单位、飞机、舰艇)。对资产进行定位和识别能够支持以下军事活动:
- 临战警告:通过定位处于打击范围内的敌方武器平台(飞机、舰艇、炮兵单位),可以判断可能发生的攻击类型以及武器接近我们的方向。
- 威胁规避:如果知道与发射机相关的威胁系统的位置,可以避开这些区域,或至少在进入这些区域时有所准备。
- 电子对抗措施的选择与实施:威胁的定位与识别将决定哪些对抗措施有效,以及应在何时启动。
- 电子战序列(electronic order of battle)的构建:不同类型的单位具有不同类型的发射机(不同频率与调制),因此辐射源在一定精度下可以识别敌方单位的类型。通过了解敌方单位的位置、识别信息以及其近期的行动历史,可以估计敌军的兵力结构,甚至预测其行动。(例如,进攻与防御需要的部队类型和数量是不同的。)
- 目标指示:如果具备非常精确的辐射源位置信息,就有可能对超视距的重要敌方资产实施武器打击。通过与辐射源关联的资产识别,也能比单纯依靠光学方法获得更准确的目标识别。
- 窄视场侦察手段的提示:光学传感器通常在视场与分辨率之间存在显著权衡,因此常被形象地称为“吸管(soda straw)”式传感器(即好像通过吸管观察物体)。这意味着在这种窄视场中寻找敌方资产可能需要很长时间。电子辐射源定位系统可以通过提供目标大致位置来缩小搜索区域。
本章首先回顾辐射源定位技术,但仅限于支撑本章重点,即辐射源定位系统精度的规定与测量方式。关于辐射源定位系统的更详细讨论可见于《EW 101》的第8章。
精度通常以系统的均方根误差(RMS error)表述,这是一种被广泛接受的系统有效精度定义。均方根误差通过收集大量角度误差或位置误差数据(覆盖多个频率、多个角度,或两者兼有)来确定。每个数据值平方后,计算其平方均值的平方根。(关于均方根误差的更多内容见第6.4.1节。)
6.1 基本的辐射源定位方法¶
如图6.1所示,辐射源定位有三种基本方法。该图展示的是二维平面(例如地球表面上的)定位方法。当然,所有方法都可以扩展到三维。第一种方法[图6.1(a)]是三角测量,即从两个(或多个)已知位置引出方位线到辐射源,交点即为辐射源位置。第二种方法[图6.1(b)]是测定一条方位线与一条距离。第三种方法[图6.1(c)]用于精密定位系统,其通过确定两条数学描述的曲线,二者交点即为辐射源位置。

图6.1 辐射源定位的三种基本方法分别为:(a) 三角测量,(b) 角度与距离测量,(c) 数学曲线交点。
虽然也可以通过测定辐射源到两个或多个已知位置的距离来实现定位,但由于实际限制,这种方法通常只适用于合作型辐射源。由于电子战系统的目标是非合作的敌方辐射源,因此这里忽略这种方法。
6.2 角度测量技术¶
在前两种方法中,必须从已知位置测定指向辐射源位置的方位线。这是通过测量信号的到达方向(DOA, direction-of-arrival)实现的,通常称为测向(DF, direction finding)。主要的测向技术包括:
- 旋转定向天线;
- 多天线幅度比较;
- Watson-Watt;
- 多普勒;
- 干涉仪。
6.2.1 旋转定向天线¶
如图6.2所示,旋转天线的增益图形是天线波束轴线角度的函数。通过让天线旋转经过辐射源,可以根据天线方位随时间的变化确定信号的到达方向。需要注意的是,天线增益曲线的形状是有详细记录的,因此只需在主波束内有两个或更多截获点,就能确定天线轴线方向,使信号位于波束中心。窄波束的大天线能够实现极高的到达方向测量精度(约为波束宽度的十分之一)。这种方法在海军电子战系统中十分常见,因为那里可以使用相对较大的天线。

图6.2 接收天线的增益图随波束轴线角度变化。
6.2.2 多天线幅度比较¶
如图6.3所示,两副朝向不同的天线在同时截获同一信号时会产生输出信号幅度比。通过该幅度比可以计算信号的到达方向。这种技术广泛应用于飞机和小型舰艇上的雷达告警接收机,因为它不需要大天线,并且能够快速确定单个脉冲的到达方向。但其精度通常较低(约在 \(5^{\circ}\) 至 \(15^{\circ}\) 之间)。

图6.3 两副天线的增益图为两个接收机提供功率比。
6.2.3 Watson-Watt 技术¶
Watson-Watt 技术由著名的雷达先驱 Robert Watson-Watt 爵士提出,它使用直线排列的三副天线。中间的天线为感应天线,两侧的天线间隔约为四分之一波长。通过将两侧天线信号的和与差(以感应天线归一化)进行处理,可以得到随角度变化的心形图案,如图6.4所示。如果有多个对称的外侧天线对,通过切换不同天线对可以旋转该心形图案,从而计算信号的到达方向。这种技术在辐射源定位系统中应用广泛,能提供中等精度的DOA(约 \(2.5^{\circ}\) RMS)。
6.2.4 多普勒技术¶
该方法在两副天线中测量接收信号的频率,其中一副天线绕另一副旋转,如图6.5所示。由于运动天线存在与传输距离变化率成正比的多普勒频移,因此其接收到的信号频率与发射频率不同,而静止天线接收到的则是发射频率。运动天线 \((A)\) 的圆周运动会在其接收的频率与天线B接收的频率之间产生正弦变化的多普勒频移。辐射源的到达方向就是多普勒频移从正变负的角度。在实际系统中,旋转天线通常被顺序切换到接收机的环形天线阵列所取代。通过比较相邻天线切换时的信号相位,可以计算频率。这种技术广泛应用于民用水面船只的无线电测向系统,通常能提供 \(3^{\circ}\) 或更大的RMS精度。

图6.4 在 Watson-Watt 测向系统中,两侧天线的和与差图形产生心形图案。

图6.5 在多普勒测向系统中,天线A绕天线B旋转,产生一个相对于天线 \(B\) 接收频率的正弦频率偏移。
6.2.5 距离测量技术¶
如果已知发射功率和接收功率水平,就可以计算信号传播的距离——如图6.6所示。由于该技术仅用于不需要高精度距离测量的电子战系统,因此通常会忽略除扩散损耗(空间损耗)之外的所有因素。扩散损耗由以下公式给出:
其中
\(L_{S}\) 为以分贝表示的扩散损耗;
\(F\) 为以兆赫兹表示的发射频率;
\(d\) 为以千米表示的传播路径长度。

图6.6 接收天线处的功率相对于有效辐射功率(ERP)按频率和距离因子衰减。
该方程可对 \(d\) 进行求解:
其中 \(\operatorname{antilog}(fn)\) 实际上为 \(10^{fn}\)。
例如,如果一部工作在 10 GHz 的雷达,其有效辐射功率为 +100 dBm,而在接收天线处接收到 -50 dBm,则扩散损耗为 150 dB。将这些数值代入公式,可得距离约为 76 km。
在实际系统中,尤其是在飞机上,这种测量的精度通常不会好于测量距离的 \(25 \%\)。
一种更精确的技术是测量传播时间。信号的传播速度接近光速(\(3 \times 10^{8} \mathrm{~m} / \mathrm{s}\)),大约为 \(1 \mathrm{ft} / \mathrm{ns}\)。因此,如果已知信号从发射天线发出和到达接收天线的时间,就能通过以下公式精确确定传播距离:
其中
\(d=\) 传播距离(米);
\(t=\) 传播时间(秒);
\(c=\) 光速 \(\left(3 \times 10^{8} \mathrm{~m} / \mathrm{s}\right)\)。
例如,如果传播时间为 1 ms,则距离为 300 km。
这就是雷达测距的方式——因为发射机和接收机通常共址,所以很容易实现。然而,在单向通信中则要困难得多。问题在于如何精确确定发射时间和到达时间。到达时间的问题已基本通过高精度GPS时钟解决,但发射时间只能在合作系统(如GPS)中确定。
正如将在第6.3节讨论的那样,一种重要的精密辐射源定位技术(针对敌方发射源)就是基于测量信号在两个接收站的到达时间差。
6.2.6 干涉仪测向¶
当测向系统的指标为约 \(1^{\circ}\) RMS 精度时,通常使用干涉仪技术。该技术测量两副天线接收到的信号相位,并根据两者的相位差推导信号的到达方向。干涉仪测向技术的基础最好通过图6.7中的干涉三角形来说明。
两副天线形成一条基线。假设系统已知这两副天线的位置,因此可以精确计算它们的间距和朝向。考虑到达信号的“波前”。波前在自然界中并不存在,但它是一个有用的概念。波前是一条垂直于信号到达方向的线。
设发射信号为正弦波,并以光速传播。一个完整周期的长度(波长)包含 \(360^{\circ}\) 的相位。沿波前的任意点信号的相位相同,因此波前可以看作一条恒定相位线(例如正向过零点)。波长和频率的关系由以下公式给出:
其中
\(\mathrm{c}=\) 光速 \((3 \times 10^{8} \mathrm{~m} / \mathrm{s})\);
\(\lambda=\) 波长(米);
\(F=\) 信号频率(赫兹)。

图6.7 干涉三角形展示了由两副天线形成的基线几何关系,以及如何根据两副天线接收信号的相位差确定到达角。
在干涉三角形中,波前接触到一副天线,与另一副天线的距离为 \(D\)。这是由基线、波前和 \(D\) 构成的直角三角形。\(D\) 与波长的比值等于相位差与 \(360^{\circ}\) 的比值。因此,两副天线接收信号的相位差等于 \(D\) 与波长的比值(可根据接收信号的测得频率计算)。
\(D\) 与基线长度的比值等于角 \(A\) 的正弦,角 \(A\) 等于角 \(B\)。由于基线的垂线被认为是干涉仪的“零角”,所以角 \(B\) 是测得的到达角。站点的到达角还包括基线的朝向。
在“单基线”干涉仪系统(即一次使用一条基线)中,基线长度通常为 \(0.1\) 至 \(0.5 \lambda\)。基线小于 \(0.1 \lambda\) 时精度不足,而大于 \(0.5 \lambda\) 时会产生歧义。
如果天线具有 \(360^{\circ}\) 覆盖,还会产生所谓的前后歧义。来自镜像方向的信号会在天线间产生相同的相位差。解决方法是使用前后比很高的天线或使用多条基线。图6.8展示了干涉仪系统中常用的四偶极子阵列。从俯视图可见,该阵列中有六对天线(即六条基线)。正确的到达角会在不同基线的数据中得到一致的结果。
还有所谓的相关干涉仪,它使用大于 \(0.5 \lambda\) 的基线,通过相关多个基线数据解决歧义。多基线精密干涉仪则同时收集三条或更多基线的数据,这些基线的长度为多个半波长,其中两条基线的长度相差 \(0.5 \lambda\)。这种方式能够在数学上解决大量歧义。

图6.8 干涉仪测向系统中常用的四偶极子阵列,用于形成六条基线。
6.3 精密辐射源定位技术¶
精密辐射源定位技术通常指那些精度足够用于目标打击的技术,定位精度可达几十米。常用的两种精密辐射源定位技术是:到达时间差(TDOA)和到达频率差(FDOA)。它们常常结合使用,并通常与低精度定位系统配合。下面我们先关注TDOA方法。
6.3.1 到达时间差方法¶
第6.2.5节讨论了根据信号传播时间计算传播距离。信号以光速传播,因此如果我们知道信号从发射机离开和到达接收机的时间,就能知道路径长度。在处理合作信号(如GPS)或己方数据链时,信号编码可以用来确定发射时间。然而在处理敌方辐射源时,我们无法得知信号何时发射。唯一能测量的信息是信号的到达时间。但通过确定两个站点的到达时间差,可以知道发射机位于一条双曲线上。如果到达时间差测量得非常精确,辐射源位置就会非常接近这条线上,但由于双曲线是无限曲线,因此定位问题尚未完全解决。
图6.9展示了两个接收站接收来自单个发射机信号的情况。这两个站点形成一条基线。不确定区域就是可能包含目标辐射源的区域。注意,两条距离的差值决定了到达时间差。图6.10展示了无限多条双曲线中的几条。每条曲线表示一个特定的到达时间差,称为“等时线”(isochrone)。

图6.9 两个接收站形成一条基线,通过传播时间差计算辐射源的相对距离。
6.3.1.1 脉冲辐射源定位¶
如果发射信号是脉冲,确定到达时间差相对容易——前提是每个接收站都有非常精确的时钟,并且测得的时间能够传送到公共处理位置。如图6.11所示,脉冲的前沿提供了一个可靠的时间测量事件。一个问题是两个接收机必须测量到同一个脉冲。由于只需测量少量脉冲即可确定辐射源位置,因此将到达时间数据传输到处理位置(用于计算TDOA)所需的数据链带宽非常小。

图6.10 等时线是一条双曲线,它包含了所有发射机可能位于的位置,这些位置到两个站点的传播路径长度差固定,从而导致信号的到达时间差固定。

图6.11 如果在每个接收机用精确时钟测量脉冲的到达时间,则计算到达时间差非常直接。
6.3.1.2 模拟信号辐射源定位¶
现在考虑带有模拟调制的信号。这类信号具有连续的载波(在发射频率上),信息通过载波的频率、幅度或相位调制来承载。载波每个波长重复一次(通常小于1米),因此我们能够在两个接收机间相关以确定到达时间的,只有调制部分。到达时间差的确定方法是在一个接收机中以不同延迟多次采样接收信号。该时间延迟必须在足够范围内变化,以覆盖辐射源可能位于区域内的最小到最大时间差。采样信号经过数字化、加时间码后发送到公共处理点,在那里计算两个信号样本之间的相关性。
如图6.12所示,相关性随着差分延迟而变化。当差分延迟值等于到达时间差时,相关性达到峰值。需要注意的是,该相关曲线顶部相对平滑,但通常能确定到延迟步进的十分之一量级。
对于模拟信号,TDOA过程相对较慢,因为必须获取大量采样;同时需要较大的数据传输带宽,因为每个采样需要较多比特才能保证足够的定位精度。

图6.12 两个接收机接收的信号之间的相关性作为其中一个接收机信号延迟的函数,在延迟值等于到达时间差时达到峰值。对于模拟信号,该曲线呈现柔和的峰值。
6.3.1.3 定位¶
确定辐射源的实际位置需要第三个接收站,以便形成至少两条基线。如图6.13所示,每条基线形成一条双曲线等时线,这两条双曲线在辐射源位置交汇。但存在定位歧义,因为两条双曲线可能在两个位置相交。然而,只有一个位置预计位于不确定区域内(见图6.9)。
为了提供精确的辐射源位置,必须确保接收站位置准确已知。随着GPS的可用性,小型车辆甚至单兵操作员都能获得准确位置。如果接收机在移动,当然必须在绘制等时线和进行辐射源定位计算时考虑接收机的瞬时位置。
6.3.2 基于FDOA的精密辐射源定位¶
FDOA是实现精密辐射源定位的一种技术。它涉及测量来自单个发射机(通常是静止的)到两个移动接收机的接收频率差。由于接收频率差是由多普勒频移差引起的,FDOA也称为差分多普勒(DD)。
6.3.2.1 到达频率差方法¶
首先考虑从固定发射机接收信号时,若接收机在运动,接收信号的频率情况。如图6.14所示,接收信号的频率取决于发射频率、接收机速度以及发射机与接收机速度矢量之间的真实球面角。接收信号频率由下式给出:

图6.13 三个接收站可形成第二条基线。这两条基线的等时线在辐射源位置相交。

图6.14 移动接收机接收来自固定辐射源的信号,由于速度和角度 \(\theta\) 的作用,信号发生多普勒频移。
其中
\(F_{R}=\) 接收频率;
\(F_{T}=\) 发射频率;
\(V_{R}=\) 接收机速度;
\(\theta=\) 从接收机速度矢量指向发射机的角度;
\(\mathrm{c}=\) 光速。
现在考虑两个运动中的接收机从不同位置接收同一信号,如图6.15所示。两个接收机的瞬时位置形成一条基线。两个接收机接收频率的差异是 \(\theta_{1}\) 和 \(\theta_{2}\) 之间的差异以及接收机速度矢量的函数。两个接收机接收频率差由下式给出:

图6.15 两个运动中的接收机根据其速度矢量和截获几何关系接收到不同的频率。
其中
\(\Delta F=\) 频率差;
\(F_{T}=\) 发射机频率;
\(V_{1}=\) 接收机1的速度;
\(V_{2}=\) 接收机2的速度;
\(\theta_{1}=\) 从接收机1速度矢量指向发射机的真实球面角;
\(\theta_{2}=\) 从接收机2速度矢量指向发射机的真实球面角;
\(\mathrm{c}=\) 光速。
在三维空间中存在一个曲面,它定义了在现有条件下会产生所测频率差的所有可能发射机位置。如果我们将该曲面与一个平面(例如地球表面)相交,得到的曲线通常称为“等频线”(isofreq)。两个接收机可以以不同方向、不同速度运动,系统计算机可以针对每种速度/几何/频率差条件绘制正确的等频线。然而,为了便于人类直观理解,图6.16展示了一组等频线,这些等频线对应于两个以相同速度、相同方向运动的接收机(但不一定是追逐场景)下的不同频率差。需要注意的是,这些曲线充满了整个空间,看起来像高中物理课本中条形磁铁两端的磁通线。
与TDOA类似,两个接收机的频率差测量并不能直接给出发射机位置,而只能确定可能位置的曲线(即等频线)。然而,如果频率差测量足够精确,发射机的位置将非常接近该等频线(量级约为 \(50 \mathrm{~m}\))。如果有第三个运动接收机,就会有三条测量基线,每条基线都能收集FDOA数据并计算等频线。此时可以通过两条或多条基线等频线的交点来确定发射机的位置。
6.3.3 针对运动发射机的FDOA¶
在利用FDOA对运动发射机(使用运动接收机)进行定位时存在一个显著问题。所测频率差来自于接收机速度矢量已知的多普勒频移。如果发射机也在运动,它会产生与接收机运动同量级的多普勒频移——但发射机的速度矢量未知。这为辐射源定位计算引入了另一个变量。虽然数学上可以求解,但所需的计算(即计算机能力和时间)要复杂得多。因此,FDOA通常被认为只适用于从运动的空中接收机定位固定或运动非常缓慢的发射机。

图6.16 使用两个接收机时,FDOA系统确定一条通过辐射源位置的曲线(等频线)。
6.3.4 FDOA与TDOA的联合应用¶
由于频率测量和时间测量都需要高精度频率基准,因此从同一对接收机同时执行这两种功能是合乎逻辑的。这在许多精密定位系统中得以实现。考虑图6.17,这是由两个接收机组成的一条基线计算出的等时线(TDOA)和等频线(FDOA)集合。可以注意到,发射机位置位于一条等时线与一条等频线的交点处。因此,精密辐射源定位可以由两个接收机的一条基线来确定。
在实际中,定位系统通常使用三个或更多平台,因此可以分别通过TDOA、FDOA或TDOA与FDOA的联合来计算多个解。这种多解性使得系统能够在广泛的作战条件下提供更精确的最终结果。

图6.17 辐射源位置可由同一对运动接收机的TDOA和FDOA确定。
6.4 辐射源定位——定位精度报告¶
在比较测向系统的作战价值时,一个关键参数是有效精度。对于到达角(AOA)系统,这种精度通常以均方根(RMS)角误差表示。在完整的辐射源定位系统(例如多个AOA测向站)中,辐射源定位精度通常用CEP或EEP表示。
当辐射源位置报告给决策者时(例如,指挥官试图确定敌方资产的位置),CEP或EEP定义了辐射源定位系统所报告测量结果的不确定性。在地图显示上绘制的圆或椭圆可以对决策过程起到辅助作用。
6.4.1 均方根(RMS)误差¶
在任何到达角或辐射源定位测量中总是存在一定误差,但我们需要能够评估和描述系统的“有效”误差。如果要求每一次测量都小于某一特定角误差,就会采用峰值误差。然而,在实际的测向系统中,特别是那些能瞬时覆盖 \(360^{\circ}\) 的系统,可能会在某些特定角度和频率下出现明显大于平均值的测量误差。这几乎总是在现场测试中发生,因为低水平干扰信号或点反射体可能导致峰值误差,而这些误差与系统本身性能无关。如果这些峰值误差(无论来自系统内部还是外部)只出现在少数角度/频率组合上,通常不会被认为是对系统作战实用性的公平检验。因此,我们通常认为系统的有效精度更适合用均方根(RMS)误差来表示。
为了确定RMS误差,需要在系统全角范围(通常为 \(360^{\circ}\))的多个不同角度以及系统全频率范围的规律增量上采集数据。每个测得的到达角与实际到达角比较,得到误差角。实际角度可由系统停放在转台上的位置、系统所在飞机或舰艇的导航系统,或其他独立角度基准确定。然后对每个误差平方,取其平方平均值,再开平方根。公式为:
RMS误差通常被认为由两部分组成:均值和标准差。均值误差只是所有误差测量的平均值——它可以在所有输出数据中加以修正。而标准差是当均值误差从每个数据点中减去后计算的RMS误差(即不含均值误差分量的RMS误差)。RMS误差、均值误差与标准差的关系为:
其中
\(\mu=\) 数据点的均值;
\(\sigma=\) 标准差。
如果尝试套用该公式,会发现数列1、4、6、8和12的均值为6.2,标准差为3.7,RMS为7.22。
6.4.2 圆概率误差(CEP)¶
CEP是炮兵与轰炸领域的术语。如果将若干炮弹或炸弹瞄准某个目标点,并测量每个命中点与该点的距离,则圆概率误差为能包含一半命中点的圆的半径。
在评估辐射源定位系统性能时,我们使用CEP来量化定位精度。在这种情况下,CEP是以计算出的辐射源位置为圆心的圆的半径,该圆有50%的概率包含实际的辐射源位置。这也称为 \(50 \%\) CEP;如果该圆有 \(90 \%\) 的概率包含辐射源,则称为 \(90 \%\) CEP。
考虑图6.18所示的辐射源定位几何。这是一种理想几何,因为两个测向站与辐射源的距离相等,并且在辐射源视角下相隔90度。
假设每个测向接收机的误差均为零均值的正态分布,则定义RMS误差的线将包含沿中心线的 \(68.2 \%\) 的测向读数。这是因为正态分布曲线在标准差点内的面积为0.341。因此,在这种几何下,正好内切RMS线的圆应有 \(46.5 \%\) 的概率包含实际辐射源位置(\(0.682 \times 0.682=0.465\))。然而,如果将该圆半径增加到原来的1.037倍,则包含辐射源的概率为 \(50 \%\)。这是因为正态分布曲线在 \(1.037\sigma\) 处的面积为 \(35.36 \%\),如图6.19所示 \(\left([0.3536 \times 2]^{2}=0.5\right)\)。

图6.18 在 \(90^{\circ}\) 夹角条件下,计算位置为圆心的 \(1\sigma\) 圆有 \(46.5 \%\) 的概率包含辐射源位置。
为进一步澄清前文概念,考虑到从每个测向站观察,实际辐射源位置落在 \(\pm 1.037\sigma\) 误差范围内的概率为 \(70.7 \%\)。对于两个最佳方位布设的测向站,这一概率为 \(0.707^2=0.5\)。
6.4.3 椭圆概率误差(EEP)¶
现在考虑图6.20所示的较不理想的辐射源定位几何。在这种情况下,尺寸为刚好内切RMS误差线的椭圆的103.6\%的椭圆将定义EEP。与CEP类似,EEP也可以定义为 \(50 \%\) 或 \(90 \%\) 概率。
6.5 辐射源定位——误差预算¶
辐射源定位系统价值的最重要衡量指标通常被认为是定位精度。在系统规范中,必须考虑所有引入误差的因素。这称为误差预算。一些因素在多种辐射源定位方法中通用,但许多因素只与某一种方法相关。

图6.19 对于均值为零的正态分布误差,高斯曲线下的面积定义了实际位置在给定半径内的概率。

图6.20 在这种非理想几何情况下,发射机位于 \(1\sigma\) 椭圆内的概率为 \(46.5 \%\)。
6.5.1 误差因素的组合¶
辐射源定位误差有多个来源;有些是随机的,有些是固定的。通常情况下,如果误差来源随机且相互独立,则它们可以统计组合。总误差为各组成部分平方和的平方根,如下所示:
其中有 \(n\) 个独立随机误差来源。
然而,如果误差来源不是随机的,就必须直接相加。
如前所述,RMS误差由均值误差和标准差组成。在某些情况下,例如在仪器化试验场中,可以对系统误差进行非常准确和完整的测量,此时可以通过将所有定位测量或到达角测量修正均值误差来抵消偏差。这样,系统的RMS误差就等于测得误差数据的标准差。需要注意的是,这假设没有显著的场地误差——主要误差来源与平台相关。装在空中平台上的测向系统通常具有这一特性,因为机体反射会导致显著的到达角误差,而远处的多径反射造成的测量误差则较小。
6.5.2 反射对AOA误差的影响¶
靠近目标辐射源路径的反射体会通过产生多径导致误差。AOA站测量的是直达路径分量和到达天线的所有多径分量的矢量和。如图6.21所示,靠近目标辐射源的反射体会使多径信号以相对较小的偏移角到达,从而造成相对较小的误差(典型于空中系统)。然而,靠近AOA站的反射体会导致多径信号以较大角度到达,从而引起相对较大的误差。地面AOA系统受到近地形影响尤为显著。但所有AOA系统都会受到其所安装的载体(空中或地面)多径误差的影响。来自侧面的反射(相对于信号AOA的近端或对端)会造成最显著的误差。

图6.21 靠近目标辐射源的反射体会引起较小的AOA测量误差,而靠近AOA系统的反射体会引起较大的AOA误差。
6.5.3 测站位置精度¶
对于任何一种辐射源定位方法,测站位置定义的误差必须直接叠加到辐射源定位误差中,如图6.22所示。如果使用空基、舰基或地面机动测量系统,其测站位置来自平台上的惯性导航系统(INS)。直到几年前,飞机的INS位置精度会随离开固定位置(如机场或航母甲板)的时间推移而降低。然而,现代INS单元使用GPS基准进行连续位置校准。这在长时间任务中大大提高了精度。由于舰艇具有出色的导航能力,舰载测站的位置误差极小。

图6.22 测站位置的误差会直接传递到目标辐射源的位置误差中。
固定地面测站的位置经过严格测量,其位置误差主要来自塔架的静态或动态弯曲。在GPS可用之前,地面机动测站的位置精度较低,因为“低价值平台”不适合使用INS。必须停下并在已知地图位置进行架设,才能获得准确的站点位置。然而,GPS接收机体积小且定位精度极高,因此即便是卡车上的测站或单兵携带的测站,也能实现米级精度。
6.5.4 AOA辐射源定位方法的误差预算项¶
AOA系统测量的是信号到达测站方向与某一参考方向之间的角度(方位角和/或俯仰角)。其误差预算包括角度测量精度和方向参考精度,如图6.23所示。由于仰角是相对于局部竖直或水平,参考容易获得且精度很高。方位参考(通常是真北)则更具挑战性。在高价值平台(如舰艇和飞机)中,角度参考来自INS。INS具有会漂移的方向陀螺,但现代系统会利用连续GPS定位修正之间的方位向量来更新角度参考,以保证长期精度。目前已经有小型轻量的INS用于低价值平台,它们使用光纤陀螺作为角度参考,并利用GPS进行定位。虽然精度不如大型INS单元,但仍能提供高质量的方向参考和位置。

图6.23 当使用一种到达角测量技术时,信号实际的到达方向将等于测量误差与参考方向误差(如真北)之和。
在早期的地面机动平台中,必须使用磁力计作为北向参考。该设备能够感应地球磁场,可以看作是一种能电子读取的罗盘。磁力计安装在测量系统天线阵列上,因此能够测量阵列的实际朝向,即使它因风力或拉索张力而移动。然而,磁力计存在与人工罗盘相同的问题(例如磁偏角随位置变化),其精度相对较低(约 \(1.5^{\circ}\) RMS)。
AOA系统的另一个误差来源是天线阵列相对于角度参考的方向。除非角度参考装置安装在阵列上,否则每次阵列架设(如果在桅杆上)或安装在空中平台时都会引入偏差误差。
6.5.5 与信噪比相关的误差¶
通常情况下,系统定位精度是在接收强信号时规定的,但系统必须能够处理更弱的信号。评估测向系统灵敏度的一种方法是以不同接收信号强度逐步进行一系列测量(通常为5到10次)。对于强信号,所有AOA测量结果非常接近(通常相同)。当信号强度减弱时,降低的SNR会导致测得的AOA产生波动。系统灵敏度通常表述为:当一组测量的标准差等于 \(1^{\circ}\) 时的接收信号强度。也可以根据指定的SNR计算由其引起的RMS角误差分量,但这取决于具体系统配置。
6.5.6 校准误差¶
所有高精度AOA系统都需通过校准来消除来自天线安装几何、载体反射和信号处理的固定误差。校准过程包括在某种精确测试场测量AOA,并在后续操作中修正测得数据以去除校准时测得的误差。校准数据的精度本身也是角误差的一个附加来源。
6.5.7 AOA系统误差的组合¶
通常可以认为除测站位置外的所有前述误差因素都是独立且随机的。因此,它们可以统计组合以确定总体误差预算。然而,在某些情况下,某些误差必须被视为相加。例如,如果存在未修正的均值误差,就应将其视为相加误差。在某些系统中,还存在作为到达角函数的计算误差。这些误差不是随机的,如果在处理过程中未修正,也应视为相加(相对于角度)。
6.6 AOA误差到位置误差的转换¶
应用辐射源定位系统最重要的问题是它能够以何种精度对辐射源进行定位。对于任何类型的定位系统,这既依赖于测量精度,也依赖于交会几何关系。这里我们只考虑AOA系统。
6.6.1 测量精度¶
标准差是指统计意义上误差相对于均值误差的偏离角度。统计学上,这意味着 \(34 \%\) 的误差角读数与均值误差的偏离小于该角度。然而,由于误差可能出现在均值的左侧或右侧,因此两个标准差角度范围包含 \(68 \%\) 的测量误差数据。
由于均值误差主要受到携带AOA系统的平台(空中、地面或海上)的反射影响,系统校准可以消除其对最终系统RMS误差的大部分影响。因此,在已校准的AOA系统中,RMS误差主要反映的是修正数据的标准差。
如图6.24所示,单个AOA站只能将辐射源定位到一个“扇形”区域。如果将两条RMS线视为标准差线,则认为有 \(68 \%\) 的概率辐射源位于该角区域内。如果RMS误差较小,则扇形较窄,角度精度较高。然而,线性误差也取决于与AOA站的距离。定位辐射源附近误差区域的宽度可由下式计算:
其中
\(W=\) 从真实角矢量到RMS误差矢量的距离(任意单位);
\(D=\) 从站点到定位辐射源的距离(相同单位);
\(\theta=\) RMS角误差。
当使用两个AOA站时,辐射源的位置由三角测量确定。图6.25显示了两个站点角度区域的交集。理想情况下,从辐射源视角看,两个站点相隔 \(90^{\circ}\),这会使位置不确定区域最小。图中展示的是两个站点未处于理想位置的一般情况。此时,辐射源位于两个扇形交叠的菱形区域内的概率为 \(68 \%\)。

图6.24 两条RMS误差线之间的角度区域包含68\%的解,这些解来源于均值误差已消除的正态分布误差的多次计算。

图6.25 两组RMS误差线形成的菱形区域有68\%的概率包含辐射源位置。
如果对该辐射源位置进行多次计算机仿真,从两个站点随机选择误差值,得到的位置分布将形成一个椭圆区域,其中心位置密度较高,离中心越远密度逐渐减小。图6.26中的椭圆包含 \(46.5 \%\) 的仿真解。

图6.26 多次带有正态分布误差的定位仿真结果形成椭圆状散点分布,其中50\%的解位于EEP椭圆内。
6.6.2 圆概率误差(CEP)¶
CEP已在第6.4.2节定义。对于辐射源定位系统(即多个AOA传感器和必要的三角测量处理),我们可以根据图6.26的椭圆确定CEP。
第一步,需要将椭圆从 \(46.5 \%\) 包含率调整为 \(50 \%\) 包含率。为此,将椭圆的长轴和短轴同时乘以1.036。这使得椭圆有 \(50 \%\) 的概率包含辐射源。需要注意,这个椭圆可以根据AOA站的指标和交会几何计算出来,通常称为“椭圆概率误差”。
第二步是根据下式计算该椭圆的长半轴和短半轴矢量和的大小(如图6.27所示):
其中
该CEP值(据称在 \(10 \%\) 范围内估计真实CEP)出自L.H. Wegner于1971年6月在兰德公司发布的报告《R-722-PR:关于被动定位电磁辐射源的空基技术精度分析》,该报告被广泛引用。

图6.27 CEP是以测得辐射源位置为圆心的圆的半径,该圆有50%的概率包含实际辐射源。
在辐射源定位系统中,计算出的椭圆概率误差会在地图显示上绘制在测得辐射源位置周围。此时假设辐射源有 \(50 \%\) 的概率位于该椭圆内。这一信息对军事分析员非常有价值,因为他们可以结合其他战术态势和近期历史信息,做出精确度合适的战术判断。
另一方面,CEP在比较评估不同系统和战术的实际几何场景时最有用。
6.7 精密定位系统中的定位误差¶
TDOA和FDOA是第6.3节讨论的精密辐射源定位方法。
TDOA和FDOA系统计算辐射源位置的精度通常用CEP(如前所述)来表示,但CEP是根据等时线和等频线的计算精度确定的。与AOA系统中讨论的定位精度类似,TDOA和FDOA的CEP依赖于测量精度和几何条件。精度计算基于一个假设:进行了大量独立测量。因此,精度是统计意义上的。
6.7.1 TDOA系统精度¶
图6.28展示了辐射源和两个TDOA接收站的位置。等时线是平面上所有可能的位置曲线,这些位置会产生所测得的到达时间差(信号以光速传播)。位置误差预算包括站点位置的精度和时间测量的精度。以下方程计算了绘制该双曲线精度的标准差,假设站点位置是精确的。因此,它仅依赖于交会几何和时间测量精度。时间测量误差假设为高斯分布。这有时被描述为双曲线的“厚度”。为了避免将其写成一个庞大的公式,我们分几步进行。
如图6.28所示,有两个站点,沿坐标系的 \(x\) 轴相距 \(B_{1}\)。辐射源位于 \(X, Y\) 位置(相对于位于站点1的坐标原点,以千米为单位)。

图6.28 到达1号与2号站点的信号的TDOA定义了一条通过辐射源的双曲线。一个 \(1 \sigma\) 的TDOA误差会使该双曲线移动距离 \(E_{1}\)。
首先,计算两条信号路径的长度:
然后,计算该三角形边长和的一半 \(\left(S_{1}\right)\):
接下来,计算从辐射源看站点1和2之间夹角的一半正弦值 \(\left(\theta_{1}\right)\):
现在可以写出偏移量(以千米计)的标准差(\(1\sigma\))公式,即误差双曲线与真实通过辐射源的双曲线之间的偏移量。在信号位置的偏移距离称为 \(E_{1}\)。其中 \(\Delta t\) 项是TDOA测量的 \(1\sigma\) 时间误差。
6.7.1.1 TDOA的CEP¶
如第6.3.1节所述,双曲线只是一条通过辐射源的线,必须从另一条基线绘制第二条双曲线,才能通过交点确定位置。图6.29展示了引入第三个截获站的情形。我们将计算第二条基线的双曲线偏移,并绘制由此产生的误差区域。为简化处理,假设第三个站点位于\(x\)轴上,距离站点2为\(B_{2}\)。同样,计算分步进行,首先计算第三条信号路径长度:
然后,计算第二个三角形边长和的一半 \(\left(S_{2}\right)\):
接着,计算从辐射源看站点2和站点3之间夹角一半的正弦值 \(\left(\theta_{2}\right)\):
现在可以写出偏移量(以千米为单位)的标准差(\(1\sigma\))公式,即误差双曲线与通过辐射源的真实双曲线之间的偏移量。在信号位置的偏移距离称为 \(E_{2}\)。其中 \(\Delta t\) 项是TDOA测量的 \(1\sigma\) 时间误差。

图6.29 引入第三个站点后定义了第二条双曲线。两条双曲线在辐射源位置相交。可由基线 \(B_{1}\) 计算出一组 \(\pm 1\sigma\) 的双曲线偏移。
在辐射源位置,来自每条双曲线的\(\pm 1\sigma\) 误差线相交所围成的区域是图6.30所示的平行四边形。与AOA的CEP计算类似,进行随机TDOA测量误差(高斯分布)的计算机仿真将生成如图6.31所示的椭圆状点密度散布图。包含50\%解的椭圆即为EEP,而CEP可以按照第6.6.2节所述的方法由EEP计算得出。

图6.30 来自两条TDOA基线的 \(\pm 1\sigma\) 误差线在辐射源观察下相交角度为两半角之和,形成平行四边形。

图6.31 多次带有正态分布误差的TDOA测量仿真结果形成椭圆状分布。包含50\%解的椭圆即为EEP。
6.7.2 FDOA辐射源定位系统中的定位误差¶
与之前讨论的TDOA系统类似,FDOA系统的精度通常以定位的CEP或EEP来衡量。与AOA和TDOA一样,精度计算的假设是进行了大量独立测量。
这里给出的等频线(isofreq)精度公式基于Paul Chestnut博士发表于1982年3月《IEEE Transactions on Aerospace and Electronic Systems》上的一篇经典文章。此处的公式假设发射源距离传感器平台足够远,以至于可以忽略接收平台观测发射源的俯角。同时假设每个接收平台的位置和速度矢量是精确已知的。最后,为简化问题,还假设接收平台以相同速度沿水平基线直线飞行。(需要注意的是,Chestnut博士的公式更复杂,因为它们不做这些简化假设。)
6.7.2.1 等频线精度¶
图6.32展示了辐射源和两个FDOA接收站的位置。等频线是平面上所有可能位置的曲线,这些位置会在该截获情形下产生所测得的到达频率差,如第6.3.2节所述。以下方程计算了绘制该曲线的精度标准差(基于前述假设)。因此,它仅依赖于交会几何和频率测量精度。频率测量误差假设为高斯分布。这有时被描述为等频线的“厚度”。需要注意的是,等频线在真实发射源位置处垂直于夹角 \(\theta_{1}\) 的平分线,因此测量误差会导致曲线沿平分线方向偏移。为了避免写成一个庞大的公式,我们分几步进行。
如图6.32所示,有两个站点,沿坐标系\(x\)轴相距 \(B_{1}\)。辐射源位于相对于站点1(坐标原点)的\(X, Y\)位置(以千米计)。
首先,使用第6.7.1节TDOA案例中的相同公式计算两条信号路径长度:

图6.32 FDOA系统中,来自站点1和站点2的信号定义了一条通过辐射源的曲线(称为等频线)。\(\pm 1\sigma\) 的FDOA误差会使等频线沿角度 \(\theta_{1}\) 的平分线方向偏移距离 \(E_{1}\)。
接着,计算辐射源方向与各接收平台速度矢量之间的夹角:
然后,计算从辐射源到各移动接收平台连线的角速度(\(\hat{\alpha}\)):
现在可以写出\(E_{1}\)的公式,即错误等频线相对于真实通过辐射源的等频线的偏移量标准差(\(1\sigma\),以千米计),其中\(F\)为发射频率(赫兹),\(\Delta F\)为FDOA测量的\(1\sigma\)频率误差(赫兹):
注意,真实等频线和错误等频线都垂直于角度 \(\theta_{1}\) 的平分线。
6.7.2.2 FDOA定位的CEP¶
由于定位需要第三个测站来绘制第二条等频线,使其与第一条等频线在辐射源位置相交,因此考虑图6.33。
我们需要按前述步骤计算 \(D_{3}, A_{3}, \hat{\alpha}_{3}\) 和 \(\theta_{2}\):
现在,可以通过下式计算第二条等频线的误差:

图6.33 引入第三个站点后定义了第二条等频线。两条等频线在辐射源位置相交。可由基线 \(B_{2}\) 计算出一组 \(\pm 1\sigma\) 曲线偏移。
图6.34中的平行四边形展示了来自两条FDOA基线的\(\pm 1\sigma\) 错误等频线。由于误差(\(E_{1}\) 和 \(E_{2}\))为高斯分布,因此测得值有 \(46.5 \%\) 的概率位于\(\pm 1\sigma\)误差线之间。如果将 \(E_{1}\) 和 \(E_{2}\) 的数值乘以1.036,所得平行四边形将包含一半的数据点。与AOA和TDOA系统类似,对高斯分布FDOA测量进行计算机仿真将生成如图6.35所示的椭圆状点密度散布图。包含50\%解的椭圆即为EEP,而CEP可按第6.4节AOA情形的方法计算。

图6.34 来自两条FDOA基线的 \(\pm 1\sigma\) 误差线在辐射源观察下相交角度为两半角之和,形成平行四边形。

图6.35 多次带有正态分布误差的FDOA测量仿真结果形成椭圆状分布。包含50\%解的椭圆即为EEP。