基于压缩感知的频率捷变雷达切片干扰抑制技术研究_刘尚.pdf

第 卷 第 期（总第 期）年 月火控雷达技术 .（）.收稿日期：基金项目：国防基础科研项目（）作者简介：刘尚（），男，硕士研究生，高级工程师。研究方向为雷达信号处理、武器系统抗干扰技术。基于压缩感知的频率捷变雷达切片干扰抑制技术研究刘 尚 蒋金华 李麒麟（江南机电设计研究所 贵阳）摘 要：近些年各军事强国在战略、战术和技术上对电子战高度重视，干扰手段和干扰技术飞速发展，我方雷达探测装备在实战中面临着日益复杂的电磁干扰环境。本文主要针对雷达抗干扰的设计需求，将稀疏恢复技术与频率捷变信息处理技术相结合，基于压缩感知理论稀疏恢复算法开展技术研究。首先，通过建立时 频域捷变目标回波模型，将观测区域内的目标散射点的距离和速度离散化为多个二维分辨单元，将分辨单元里的散射点强度作为原始信号，而将观测采样矩阵的每一列作为对应单位强度散射点的回波向量。然后，利用贪婪追踪算法完成稀疏恢复。最后，通过仿真结果表明基于稀疏恢复的捷变雷达信息处理方法具有较好的目标分辨能力，可在较高的干信比下完成切片干扰的抑制，验证了该方法的可行性及有效性。关键词：压缩感知；稀疏恢复；正交匹配追踪算法；频率捷变；切片转发干扰中图分类号：文献标志码：文章编号：（）引用格式：刘尚，蒋金华，李麒麟.基于压缩感知的频率捷变雷达切片干扰抑制技术研究.火控雷达技术，（）：.：.，（，）：，.，.，.，.，.，.，.，.：；火 控 雷 达 技 术第 卷 引言随着高速采样、超大规模数字处理、射频直接存储等技术的进步，电子干扰技术也快速发展。干扰机产生的干扰波形能够实现与雷达探测波形相参，形成对雷达回波信号在多个域上的全面覆盖。不仅抵消了雷达匹配滤波带来的抗干扰增益，而且大幅压缩了雷达的探测距离。因此，日益复杂的强对抗电磁环境对雷达的抗干扰技术能力提出了更高的作战要求。近年来，信号处理领域里“压缩感知”理论的稀疏恢复技术取得了新的突破，改变了经典雷达信号处理基于“匹配滤波”的传统思想，将“稀疏恢复”思想与频率捷变信息处理技术相结合，通过发射时 频维快速、随机捷变的雷达波形，并以干扰条件下目标信息最优重建为目的进行回波信号处理，在接收处理过程中去除杂波和干扰后，通过基于非完整观测的目标信息“最优重建”，解决目标在复杂电磁干扰环境下的探测问题，改变雷达在电子对抗中的被动局面。为此，本文从压缩感知稀疏恢复算法的角度，开展了时 频维上频率捷变处理技术相关方面的研究。压缩感知稀疏恢复理论一直以来，奈奎斯特采样定理是指导如何进行信号采样的重要理论基础。该定理指出，采样速率应大于等于信号带宽，原始信号才能由采样信号不失真地精确重构。近年来压缩感知理论的研究则突破了此采样定理对信号采集与处理的限制。由、及华裔科学家 等人提出的压缩感知理论表示：如果信号通过某种变换（如傅立叶变换，小波变换等）后，是可稀疏表示或可压缩的，则可以设计一个与变换基不相关的测量矩阵测量信号，对得到的测量值求解优化，实现信号的精确或近似重构。压缩感知是针对稀疏信号或可压缩信号提出的基于数据采集和数据压缩为一体的新型信号处理方法，推动了信息领域向前发展。它与传统的信号采样、编解码的方式有着很大的不同，传统的信号采样、编码、解码的过程如图 所示。图 传统的信号采样、编解码理论框图 压缩感知理论对信号的采样、压缩编码发生在同一个步骤，利用信号的稀疏性，以远低于奈奎斯特采样率的速率对信号进行非自适应的测量编码，如图 所示。测量值并非信号本身，而是从高维到低维的投影值，从数学角度看，每个测量值是传统理论下的每个样本信号的组合函数，即一个测量值已经包括了所有样本信号的少量信息。解码过程不是编码的简单逆过程，而是利用信号稀疏分解中已有的重建方法在概率意义上实现信号的精确重建或者一定误差下的近似精确重建。解码所需测量值的数目远小于传统理论下的测量数。图 基于压缩感知理论的信号采样、编解码理论框图压缩感知理论主要包括三部分：一是信号的稀疏表示；二是设计测量矩阵对信号进行观测采样，要在降低维数的同时保证原始信号 的信息损失最小；三是稀疏恢复算法，利用 个观测值无失真地恢复出信号长度为 的原始信号。其中，第三部分信号的稀疏恢复是从 维的测量信号 恢复 维原始信号，如式（）所示。（）其中，是信号 在 域的变换向量，又称为感知矩阵。由于是稀疏的，且 ，所以求出稀疏解，再根据正交变换的可逆性，就可以很容易根据 恢复出原信号。该稀疏恢复问题可以写为式（）中的优化问题。.（）目前常用的用于求解式（）中的优化问题的算法可以简略归纳为 类：最小 范数法、贪婪算法、迭代阈值法。其中贪婪算法主要思想是通过迭代计第 期刘 尚等：基于压缩感知的频率捷变雷达切片干扰抑制技术研究算 的支撑集获得重构信号，具有低复杂度和简单几何运算的优点。文献中描述了常用的贪婪算法，包括匹配追踪（，）算法，正交匹配追踪（，）算法，压缩采样匹配追踪（，）算法，子空间追踪（）算法以及它们的改进方法。基于稀疏恢复的频率捷变处理技术“捷变”一直是雷达反侦察抗干扰的有效手段。在快速发展的干扰威胁刺激下，捷变的应用在新的信号处理理论推动下，正朝着多维、综合、全相参处理的方向进一步发展。频率捷变技术的出现改变了雷达发射单一不变波形的局面。频率捷变是指雷达发射相邻脉冲或脉冲组的载频在一定范围内快速变化，在整个相参捷变带宽内实现多个频率的同时发射，作战过程中极大改善了雷达的反侦察和抗干扰性能，成为现代雷达普遍使用的抗干扰技术手段。信号稀疏恢复技术由于利用了观测对象的稀疏特征，可以用远低于经典奈奎斯特频率的采样率获取观测对象的完整信息。在雷达信号处理方面，稀疏恢复技术首先被应用到探地、穿墙和高分遥感成像等领域，有效降低了获取目标高分辨信息所需的带宽和采样速率。此外，在多维信号处理方面，稀疏恢复技术有效降低了杂波谱估计所需的样本个数，提高了非均匀杂波中的目标检测性能。还有学者将稀疏恢复技术应用于随机步进频脉冲的高分辨成像处理中，证明了稀疏恢复技术在频率捷变波形的信号处理中的适用性和应用潜力。本文将稀疏恢复与频率捷变处理技术相结合，通过雷达发射波形在时 频维的快速随机捷变，使干扰机难以获取对波形的多维完整观测，难以在所有波形维度内形成对目标回波完整、精准破坏。在此基础上，雷达利用自身发射波形已知的信息优势，通过稀疏恢复技术利用对目标的非完整观测重建目标信息，完成目标检测跟踪和抗干扰的作战任务。.时 频域频率捷变目标回波模型传统理论认为，当发射波形捷变（尤其是频率捷变）时，目标散射特性会存在起伏，即所谓捷变“去相关”现象。分析表明，对于理想点目标基本不存在捷变引起的目标起伏。而当目标是复杂目标时，频率捷变引起的目标起伏是由于在一个分辨单元内存在多个散射中心，各中心回波的矢量叠加效果会随电磁波波长的变化而变化。因此，基于稀疏恢复的捷变雷达信号处理应结合该波形的高分辨特性，采用点散射类模型，将目标回波建模成各个散射中心回波的矢量叠加。捷变频雷达发射波形为时 频维快速、随机捷变的复杂波形，雷达以脉组为单元进行信号发射，每一个脉组包含 个脉冲，称为一个处理周期。在一个处理周期内，不同脉冲具有不同的发射载频，捷变频雷达发射的探测波形示意如图 所示。图 捷变频雷达发射探测波形示意图设原始信号 ，考虑线性采样模型，即每一次采样都是用一个线性算子作用于上，采样数据 （，）可以使用公式（）的稀疏方程表示，它描述了信号空间中的线性采样过程。以时 频维捷变为例，若将相参处理周期内的 个脉冲的回波信号写成一个 维向量，将所观测区域内的目标散射点的距离和速度离散化为 个分辨单元，假设第（，）个分辨单元上存在一强度为，的散射点，则根据脉冲信号模型可知，其 对 应 散 射 点 回 波 向 量 可 以 写 作，因此，目标多散射点回波建模仿真可以通过稀疏方程公式（）表示。（）（）（），|（）（）（）|（）其中对应散射点回波向量的元素，可表示为式（）。火 控 雷 达 技 术第 卷，（）|，（，），（，）（）（），（，），（）稀疏恢复的测量矩阵包含了散射点的距离和速度信息、信号频率调制信息等，以此构造的稀疏测量矩阵可将高维信号投影到低维空间上，尽可能多地将散射点重要的信息特征保留并完成大概率的恢复。公式（）中，对不同脉冲的发射载频信息进行了公式推导，对相关散射点的距离和速度分辨单元变量进行了分解，其中 为信号发射基带载频频率；为一个相参处理周期中脉冲个数；为一个相参处理周期中捷变频率变化量；为雷达与对应散射点之间的距离；为脉冲回波信号的序号数；为脉冲重复时间（）；为目标散射点相对雷达的速度离散化单元的速度值。公式（）中 的行向量对应的是由相参处理周期内的 个脉冲的回波信号组成 维向量，中的元素为观测场景中目标散射点的 个速度 距离分辨单元中散射点强度。中的每一列即为对应单位强度散射点的回波向量。对于捷变频雷达，发射的波形为时 频维快速、随机捷变的复杂波形，对应 中的元素产生的分辨单元信号的维数远远大于采样数据维数，利用 的稀疏性，测量矩阵 的稀疏系数中只有 个非零稀疏，且满足，因此可以通过求解稀疏方程公式（）的逆问题得到信号。.切片转发干扰信号模型切片干扰的原理是利用脉冲压缩雷达信号处理中匹配滤波带来的相干增益来提高切片干扰的压制或者欺骗效果，其实现途径是基于数字射频存储器、直接数字合成器对雷达信号进行不失真采样、存储、信号处理和形成模拟信号后转发，这样转发的干扰信号同雷达信号高度相参，干扰信号通过匹配滤波后可以形成单个假目标或多个假目标群。假设雷达信号为线性调频信号（），脉宽为，带宽为，调制斜率 为 ，则有如式（）所示。（）()（）其中，()，其它。假设真目标距离为，则雷达至目标的单程延时为 ，其中 为光速。由此可得雷达的目标回波信号（）为（）（）|（）（）其中表示目标回波幅度值，目标回波（）通过匹配滤波器（）的响应（）为（）（）（）（）（）（）|（）（）（）式（）中，为卷积运算，（）()，（）表示距离为 的目标回波通过匹配滤波器（）的响应。为了简化模型，切片干扰的数学模型不考虑多普勒频移，只在时域维增加转发延时处理，可简化为将一个矩形包络脉冲串（也称为间歇采样信号）同雷达信号相乘，然后根据不同的组合方式将其转发。间歇采样信号的数学模型可表示为（）()（）（）重复转发切片式干扰的工作示意图如图 所示。图 重复转发切片干扰工作示意图因此，切片干扰信号可表示为（）（）（）()（）（）（）第 期刘 尚等：基于压缩感知的频率捷变雷达切片干扰抑制技术研究而雷达接收到的切片干扰的回波就可表示为（）（）（）假设雷达脉冲宽度为，带宽为，系统采样频率为。切片干扰的采样周期为，占空比为，转发延时为，重复转发次数为 次，干扰调制频移为，信干比为 。则对目标信号加切片干扰回波脉冲压缩后的仿真结果如图 所示。图 目标信号加切片干扰脉压后干扰效果经过干扰调制频移后的重复式切片干扰经脉冲压缩后造成的假目标在时域上 和 处各形成了一个假目标群，这两个假目标群是干扰重复转发两次所致。由此可知，切片转发干扰经过雷达脉冲压缩后可以获得众多幅度不一的假目标，这些假目标相对于真目标幅度有强有弱，再通过调频的方式，可以实现假目标超前或滞后的欺骗干扰，对雷达实战化能力非常不利。.稀疏恢复重建算法第二章中介绍了几种常用的重建算法，贪婪追踪算法是压缩感知理论中最常用的稀疏信号重建算法，其主要思想是通过迭代计算 的支撑，即通过每次迭代时进行局部最优化寻找各个非零系数来获取信号矢量的逼近。该类算法主要是采用某种准则在每次迭代过程中从冗余字典选取一个或几个与观测值余量最大相关的原子，最终寻找一组与观测值匹配的最稀疏的原子组合，从而实现信号的重建。其算法的优点是对于维数较低的小尺度信号问题运算速度很快，计算复杂度相对较低。为了完成对时 频维捷变多散射点回波稀疏重构研究中稀疏方程的求解，即求解以下的优化问题：，.（）式（）中，模范数定义为（）得到的解 满足如式（）误差关系为 （）即是 的最佳 阶稀疏逼近，即 ，.（）基于以上稀疏方程求解的思路，本文选择相对收敛速度较快且重建算法中使用最多的正交匹配