１引言
现代战争中，警戒雷达遭遇最严重的威胁来自有源干扰。由于威胁目标的雷达回波可能被掩盖在干扰杂波中，导致指挥系统不能确定空空或空海导弹的发射时机，无法引导飞机在安全距离上完成攻击任务。为了打击和摧毁干扰辐射源，辐射源定位技术伴随雷达电子对抗产生。辐射源定位技术最基本的定位方法是平面定位法和空间定位法，根据定位条件的不同，又可分为单点定位和多点定位。单点定位有飞越目标定位法、方位／仰角定位法，多点定位有测向交叉定位法、多重采样相关定位法等［１～５］。本文提出的基于灰度质心法的雷达捷变频技术干扰源定位方法属于单点定位法。利用雷达捷变频技术，采用灰度质心算法对干扰源进行定位，即可掌握对方干扰掩护的方位，又可实现为飞机提供干扰源打击的概略引导信息。
２定位技术实现的原理及条件
２．１实现条件
基于灰度质心法的雷达捷变频技术干扰源定位方法需要在特定的条件下完成。１）雷达具有超低副瓣或副瓣抑制能力，且干扰机与雷达有一定的距离；２）雷达具有捷变频功能，且干扰信号基于频率瞄准或储频转发式；３）雷达具有数字化信号处理显示功能。超低副瓣的雷达天线设计技术不仅可以聚焦辐射功率，增大雷达的探测距离，也可以有效避免干扰信号从雷达的副瓣进入，提高雷达的抗干扰能力。这种技术迫使干扰信号仅能从雷达主瓣进入，见图１（ａ），从而使干扰信号携带了方向信息，为干扰源定位提供了基础条件。捷变频是一种抗干扰技术，即雷达快速改变射频频率，以迫使远距离干扰机和威胁目标自带的干扰机分散其现有的干扰功率，来覆盖增宽了的射频带宽。如果远距离干扰机坚持采用窄带瞄频或数字储频转发方式产生干扰信号，干扰信号进入雷达接收机就存在一个延迟时间，这个延迟由干扰机与雷达之间的双程信号传输距离、干扰机侦察测频处理时间、干扰机内部射频线长度等因素决定。脉间或脉组高速捷变频技术恰好利用了这个延迟，其高速的频率跳变使远距离干扰机无法紧随雷达信号实施干扰，造成相同频率的干扰信号在干扰机与雷达距离之内传播滞后于雷达回波，此距离内的雷达回波处理未受到干扰，形成如图１（ｃ）的效果［６～７］。雷达数字化显示处理使干扰杂波与目标回波以点迹形式在雷达显示，有利于对干扰源位置的数据采集，减少误差，如图１（ｂ）、１（ｃ）。
２．２实现原理
在真实战场条件中，当雷达与干扰机距离较远时，大多数情况，干扰能量仅能从雷达的主瓣和第一副瓣进入，并在雷达显示器上形成主瓣式条状干扰扇区。如果雷达采用捷变频抗干扰手段有效，往往可以消弱干扰的能量，并能去除干扰机与雷达之间的干扰亮区。这时，可以清晰观测干扰亮区的起始位置，在经过十几个扫描周期的采样后，可构画出干扰机位置分布图，之后使用定位算法可计算出干扰机位置。如果针对一部雷达进行多次干扰源定位误差校准，可进一步缩小结果误差，减少计算所需的数据采集样本，达成对运动干扰源的准实时定位［８～１０］。过程原理如图２所示。
３定位技术方法
采用灰度质心算法的雷达捷变频技术对干扰源进行定位方法［１１］，是一种以灰度为权值的加权型心法。干扰源位置以干扰噪点的起始点为准，数量Ｎ＝ｎ×Ｔ／ｔ，其中Ｔ为采集区间的时间长度，ｔ为雷达天线的扫描周期，ｎ为每个周期干扰源数量。例如，图２（ｃ）是雷达的某次扫描图像，图上有ｎ＝６个干扰源位置点。表１中列出了数据的部分样本。为了便于计算，数据由经纬度值转换为Ｘ／Ｙ坐标系，由于无高度信息，Ｚ坐标省略。根据式（１）、式（２），计算并画出干扰机位置，效果如图３所示，式（１）得到数据在Ｘ轴投影的质心，式（２）得到数据在Ｙ轴投影的质心。由于数据分布集中在两个位置，在计算之前应划分数据区域１、区域２，并且利用式（４）剔除偏离点，这样就能在图中得到两个质心Ｓ１（１６４５．８３，１０１４．９），Ｓ２（１１４２．６７，１２６２．６），以”＋”表示；干扰机真实位置Ｒ１＝（１６０８，１１０７），Ｒ２＝（１１０４，１４０１），以“＊”表示。
４定位结果分析
根据图３所示，比例尺为１∶９１．８６（单位：ｍ），区域１的质心距离实际点９１４６．４１ｍ，区域２的质心距离实际点１３２００．６８ｍ。由于无校准的定位误差较大，需要进行进一步处理才能获得理想结果［１２～１３］。
４．１随机误差
采样点相对质心的分布趋势表现数据的随机性，图４、图５显示了位置数据的均值和波动，经计算，区域１数据的标准差为７．９５ｋｍ，区域２的标准差为１０．０９ｋｍ。图３中，区域１、区域２采样点集中分布于区域底部，构成较明显的集中区，其它分散点采样来自于副瓣干扰，在质心阈值设置时可以剔除；区域２中形成了两个集中区，从雷达屏幕观察上部集中区可能是副瓣干扰，在采样时可以滤除上部的集中区域。经过阈值优化设置后，随机误差得到降低，如图６所示。经过阈值优化设置后，区域１数据的标准差为４９１８．９ｍ，区域２的标准差为１８１０．８ｍ，区域１的质心Ｓ１（１６３９．５，１０４７．４）距离实际点６１９２．６４ｍ，区域２的质心Ｓ２（１１１８．９，１３７４．８）距离实际点２７６８．７７ｍ，误差分别降低３２．３％和７９％。事实上，Ｓ１点的干扰机较Ｓ２点干扰机的干扰信号能量及能量起伏较大，对雷达副瓣的干扰不稳定，使得雷达对Ｓ１的定位不准确，因此在数据采样时可以考虑选择距离最新的干扰点，即每次扫描周期只录取干扰方向上最近的干扰点。这种录取方法操作简单，结果精度高，但对分辨主瓣干扰的能力要求较高。
４．２系统误差
由图３、图６所示，质心在同一个方向上偏离实际点，这是由于干扰机内部在转发雷达信号或测量雷达信号存在系统延迟Ｔ，这个延迟导致雷达观测到的干扰机位置比实际的位置要远一些，不同的干扰机系统延迟Ｔ是不同的，同一部干扰机工作在不同频率时Ｔ也会不同。另外一个系统误差来自雷达本身，由于雷达精度与量程、天线扫描方式、工作频率等有关，在对干扰机定位时应考虑以上因素。
５结语
数据处理结果表明雷达捷变频技术的干扰源定位方法不能实现高精度定位，该方法的实际应用价值在于可以利用现有雷达装备提供干扰威胁源位置的估计信息，定位误差小于５ｋｍ。该方法不属于高精度雷达目标跟踪功能，而是一种雷达探测能力的延伸，其最大的优势在于可以估计运动干扰源的移动线路，比如描绘携带干扰源飞机的航线，有利于对敏感空域或建筑经济期刊海域敌情的侦察预警。从经济条件和现有技术条件来看，该方法在雷达上推广易实现，技术改造需求少，是一种行之有效的雷达辅助探测功能。
作者：魏永峰 单位：９１４０４部队９２分队
