1、第 卷 第 期（总第 期）年 月火控雷达技术 .（）.收稿日期：作者简介：李哲（），男，硕士研究生。研究方向为雷达总体设计。基于中值对消 算法的机载雷达杂波抑制研究李 哲 王 伟 赵 娟 党腾飞 卢珊珊（.西安电子工程研究所 西安；.西安现代控制技术研究所 西安）摘 要：杂波抑制干涉（，）是一种常见的雷达信号处理方法，其高效抑制雷达接收回波中的地杂波干扰，可大大提高机载合成孔径雷达地面动目标指示（）的信杂噪比。然而，由于杂波抑制算法要实时进行大量滤波权系数的计算，因此其在工程中的实现较为困难。基于此研究背景，本文对机载雷达地杂波进行了建模，分析了、（）和中值对消 三种杂波抑制算法的原理和性能，
2、并进行了工程实现。由工程实验结果可知，三种算法均可提高机载雷达信杂噪比，从而提高动目标识别能力，其中中值对消 算法虚警概率最低，性能提升最大。关键词：机载雷达；杂波抑制干涉；中值对消 算法中图分类号：文献标志码：文章编号：（）引用格式：李哲，王伟，赵娟，党腾飞，卢珊珊.基于中值对消 算法的机载雷达杂波抑制研究.火控雷达技术，（）：.：.，（，；，）：，（），（），（），：；引言合成孔径雷达（）是一种空间分辨率极佳的微波成像雷达，可全天时全天候工作，通常应用于对地或对海的静止场景成像，在民用观测和军事领域都有着举足轻重的作用。然而在工程实际中雷达往往工作在较为复杂的环境下，因此接收回波中会存火
3、控 雷 达 技 术第 卷在大量的杂波和噪声分量，直接影响了雷达的信杂噪比，从而影响到雷达的动目标分辨性能。因此，为了得到更高精度的动目标检测结果，地杂波抑制的重要性便不言而喻。常见的杂波抑制技术有时间平均相干机载雷达（）技术以及偏置相位中心天线（）技术等。然而，由于杂波抑制算法要实时进行大量滤波权系数的计算，因此其较难应用在实际雷达工程信号处理中。故本文着重选择了、中值对消 三种算法加以研究，首先对机载雷达地杂波进行了建模，然后分析了上述三种算法的原理和理论性能，最后通过 将其应用在实际工程中，以工程实际数据来验证并比较几种算法的优越性。以此为思路，首先，我们需要对机载雷达地杂波的特征加以研究
4、。机载雷达地杂波模型如图 所示，为机载雷达回波数据采集几何关系。接收中心设置在由均线阵（）接收回波信息的最大速率、平台高度 的水平面上进行匀速直线运动的载机上，它的阵元数为，阵元间隔为，并设雷达波长为，脉冲重复周期（）为，在单个 内以 个脉冲发送和接收。图 机载雷达回波数据几何关系如果把全部杂波按方位角等分成个杂波源，并认为每个杂波源都彼此独立同分布，在时域中是彼此相关的。则对于没有距离模糊的系统来说，其第 个杂波源的 维空间时间导引向量，可以描述为（，）（）（）其中，（），（），（），（）为杂波的时域和空域导引向量；表示向量的 积；上标 为矩阵求转置，、则是对第 个杂波源的归一化多普勒、空间
5、频率，由其俯仰角，方位角和偏航角度确定，即 （）|（）经过距离脉压后，此距离内杂波的回波信号可用所有杂波源的线性组合来表示为 （，）（）其中（，），（，），（，），而系数矩阵 ，为第 个杂波源的对应的复幅，服从零均值复高斯分布，方差为。所以，在此距离门的杂波协方差矩阵，可以描述为 （）其中 （，），（）表示对角线元素由矢量 构成的对角阵。杂波抑制干涉（）.双通道 信号模型图 所示为机载双通道 数据采集的几何和空间模型。确定载机沿 轴航行，雷达天线在正侧视方向上运行，波长为，相位中心为，再由子口径、按航迹位置依次排列接收信息，从子口径相位中央到 的距离都是。假定某一发散点在场地上为，则其方位角和
6、俯仰角在发射孔相位中央所对应的位置，依次为、。从图中可以发现，由于接收位置在航迹上并不相同，导致当 点回波达到孔径、处时，就会形成一条路径差 ，而因为双极化天线孔径间距 远低于成像间距，图中 和 两线也近乎平行，因此 也可表达为 （）这样，点的回波信号传输到、点时将会有不同的相位差为（）第 期李 哲等：基于中值对消 算法的机载雷达杂波抑制研究图 双通道 回波信息接收空间关系图 假定 是杂波，其由于载机运行引起的相对运动形成的普勒频率就可以描述为（）即当 点为具有一定径向速度的运动目标时，其多普勒频率由目标本身和载机运动联合确定，可以描述为式（）所示。（）由此可见，同一个位置的杂波与不同运动时刻
7、的目标之间存在着截然不同的多普勒频率，这样在经过 成像后，会处在同一图像中不同的多普勒方位单元内。同理，实际位置处于不同方位的运动目标和杂波，虽然其可能均处于图像中同一像素单元上，但二者由于不同通道接收造成的相位差也不同。多通道的 系统正是由于这一点，才可进行空、时二域的协同处理，这样一来，便可使临近杂波空间的慢动目标的检测更好地进行。在进行双通道 处理时，对天线、点回波的成像中心点分别进行成像运动补偿处理。因此，、接收孔径得到的 图象中，与相应像素单元中静止目标之间具有公式（）中的相位差，所以空域的消权系数可以通过一个一个的图像单位来计算，这样就完成了杂波控制的工作。图 中提供了双通道 的信
8、号处理示意图。图 双通道 信号处理示意图 、接收孔径得到的 图像中相应像素单元的复信号分别用（，），（，）来表示，其中 为距离单元序号，为方位单元序号，该像素单元对应的空域对消权系数用（，）来表示为（，）（）即每个像素点单元对消权系数可根据系统参数和成像几何关系来计算，并借此逐像素对 通道图像进行补偿，然后对补偿后两通道图像进行相减计算：（，）（，）（，）（，）（）火 控 雷 达 技 术第 卷通道中，设某独立像素单元的杂波及动目标信号分别为、，所处实际方位位置对应的干涉相位分别为 、。那么 通道对应像素单元中杂波信号为，目标信号为，则两种分量的输出在经过上述处理后可分别表示为 （）（）（）|（
9、）也就是说杂波矢量可以全部消除，所以在对消过程中由于动目标的径向速率引入附加的多普勒频率，故能量得以保留。.自适应 技术由公式（）提供了理想情形下空域内对消权系数的准确测算，即利用系统参量和成像空间几何关系对其加以推算，这也是 有效检测的必要前提。而在实际数据处理中，由于信道幅相误差、惯导系统误差等误差的出现，简单按照公式方法（）测算所得到的对消权系数往往并不准确。此时，我们就可以利用接收数据对空域对消权系数来进行自适应的方式作出估算，我们可把在这个方式下的 算法称作 算法。.算法权系数的计算在 处理中，由邻近的像素单元信号自适应估算可得出任一像素单元空域对消权系数。在无目标出现的情况下，相邻
10、信道中与 图像对应的像素单元的复信号可由（，）（，）（，）（，）来表示，其中（，）、（，）分别是杂波和噪声分量，且彼此独立。而为了将噪杂的输出功率最小化，可按照最小均方误差原则（），选取的待测量单元的空域对消权系数符合如下要求（，）（，）（，）（，）（）求解可得（，）（，）（，）（，）（，）（）其中，上标 为共轭运算。实际处理时，由于互相关项 （，）（，）以 及 自 相 关 项 （，）（，）未知，当（，），（，）为零均值复高斯变量，且 图像中与待检测单元（，）距离相邻的训练样本（，），（，）与待检测单元杂噪分量满足独立同分布时，可用最小均方值 来代替真值，从而计算各个像素单元的空域对消权系数见
11、式（）所示。（，），（，）（，），（，）（，）（）.中值对消 算法在自适应算法中，如果想要获取空域对消权系数，就必须先在训练样本中对互相的有关参数（，）（，）及 其 自 相 关 函 数 （，）（，）作出正确估算，但在训练样本中出现干扰数据后，估算偏差将明显上升，进而造成输出 下降。由此，可将中值估计法带入 数据处理中，用式（）的最小均方法对空域相对消权关系作出了估算。但在双信道 数据处理时，由于像素单元个数大大地超过了其计算自由度，因而可以把该办法带入空域的对消处理过程中，以大大降低在扰动任务中对杂波控制的效果。根据这一点，还给出了一个大值对消 杂波的方法，在该方法中，权系数的运算可以由式（）
12、实现。（，），（，）（，）|，（，）（，）|（）其中，为取中值。可以说明，当样本数量趋于无穷时，就可以收敛为式（）中的最优权系数。为便于表达，以下过程中同样省略了式（）中像素单元序号。令：（）表达最高权系数所相应的输出信息时，其与不关联，又因，是零平均数复高斯变数，故线性形式也是零平均数复高斯变数，所以，与互相独立性。对进行求解即可得到 （）定义变量 ，则|（）第 期李 哲等：基于中值对消 算法的机载雷达杂波抑制研究对式（）中 的实部分和虚部分别取统计的数值，可得式（）。，|，|（）其中，表示对实部与虚部分别统计中值，再有|（）其中，复高斯变量，的标准差分别为和。这样一来，的分数线上下均为方差
13、为 的相互独立的零均值复高斯变量。因此，均可这样表示 实虚部的概率密度函数为，（，），（）可以看出，式（）中的函数，为，的偶函数，所以可知 的均值为零，则有，（）故，即当样本数量趋于无穷时，收敛于最优权矢量。杂波抑制算法的工程实现.信号处理系统设计由于雷达在进行回波数据处理时要接收大量的数据，故工程上雷达的信号处理机多采用 的双芯片模式，以保证数据处理的速度以及可嵌入式的程序改动。为满足算法需求，本系统选用一片 系列的，和两片飞腾公司的 内核的 处理器。图 为雷达信号处理机的工程实现方案，其特点是：采用了 的信号处理架构；芯片间的通信互联采用了 总线；外接了 存储器扩展数据存储；网口通信、系统
14、控制和数据传输采用了 协议。本系统利用了 的流水性能和多核 的数据并行处理能力，将算法合理分配到不同的处理器中。总体来说，雷达收到回波信号时，信处系统接收从前端传来的数据，通道 对原始数据进行采样和模数转换，之后传至 进行数字下变频和距离向脉冲压缩处理。脉压后的基带信号以脉冲重复频率通过 接口实时将每个脉冲的回波数据发送至 进行其他数据处理。图 和图 分别为雷达工程实现所用 板卡和 程序的调试界面。.实测数据处理结果为了经由雷达 成像和动目标显示功能验证上一章节中三种算法的性能，试验选取陕西省某地进行 成像，所用雷达系统参数如表 所示，另有三辆车辆作为配合动目标在成像区域省道上以 左右的时速匀
15、速行驶。图 雷达信号处理机实现方案火 控 雷 达 技 术第 卷图 工程实现所用的 板卡图 调试界面表 试验所用系统参数列表参数类型取值（范围）雷达工作频率成像分辨率.系统带宽脉冲重复频率载机地速 载机飞行高度雷达作用距离主波束方位宽度接收通道数图 为三种算法对同一目标场景 成像图中的数据处理结果。其中，图（）为原始 图像，三个配合目标车辆已在图中标注出；图（）为 算法下的处理结果，如图（）所示，三个配合目标只识别出了两个，虽然识别到了一个地面非配合目标车辆，但是有高达 个虚警目标；图（）为 算法下的处理结果，如图（）所示，三个配合目标都识别了出来，且还有一个非配合目标车辆，但仍有两个虚警目标；
16、图（）为中值对消 算法下的处理结果，如图（）所示，三个配合目标均有识别，虚警目标也只有一个，而且识别到多达三个地面非配合目标车辆，由试验现场确认，三个非配合目标均为实际行驶在道路上的民用车辆。由试验结果可得：中值对消 算法比起另两种算法，其动目标信杂噪比提高更大；动目标输出功率提高也更大；表现在 图像上可清晰观测并识别到动目标强点的反射信号，且虚警率大幅降低，系统更易通过 检测。因此可得，在非均匀复杂环境中，中值对消 算法比起另外两种算法，其杂波抑制效果明显更好，动目标显示的性能也更强。第 期李 哲等：基于中值对消 算法的机载雷达杂波抑制研究图 三种算法在同一目标场景中的数据处理结果 结束语本文主要对、和中值对消 三种算法在机载 系统中的应用进行了研究，首先分析了算法的原理，然后通过 进行了工程实现。由试验结果可知，中值对消 算法比起另两种算法，其动目标信杂噪比和动目标输出功率均有提高，且虚警率大幅降低，使得杂波抑制在非均匀复杂环境中具有更好的适应性，可以认为是一种更符合工程实际的多通道 杂波抑制方法。参考文献：孙智，陈海旭，蒋兴涛，等 基于 的机载雷达海面慢速小目标信号积累与杂波抑制
