适用于机械扫描雷达的有源干扰侦察方法
【专利摘要】本发明公开了一种适应于机械扫描雷达的有源干扰侦察方法，主要解决现有技术侦察准确性低和可靠性差的问题，其实现过程为：1)获得机械扫描雷达的回波数据；2)确定干扰侦察通道回波数据中干扰回波所在的脉组及功率值；3)利用干扰侦察通道的回波数据确定干扰强度、干扰频点、干扰指向线及干扰类型；4)利用主辅通道的回波数据消除3)中的虚假指向线。本发明增强了侦察结果的准确性和可靠性，利用侦察得到的干扰源的相关信息，可以采取对应的抗干扰措施，增强雷达的抗干扰能力，可用于机械扫描体制雷达。
【专利说明】
适用于机械扫描雷达的有源干扰侦察方法
技术领域
[0001] 本发明属于雷达信号处理领域，具体的说是一种有源干扰侦察方法，可用于机械 扫描体制雷达。
【背景技术】
[0002] 随着现代战争的发展，在高技术局部战争日益发展、日趋成熟的情况下，现代雷达 的工作环境也越来越复杂，对雷达系统性能的要求也不断提高。一方面，随着电子技术的发 展和电子技术在武器装备的不断应用，军事与民用电磁波设备使用越来越频繁，设备涉及 波段也在不断增宽，造成雷达工作环境中的电磁污染愈发严重。另一方面，在现代战争中， 敌我双方在摧毁对方雷达时，还要保证自己雷达可以正常探测目标和收集情报，所以双方 将会采取各种措施，在陆上、海上乃至空中等多维空间展开干扰与反干扰的激烈斗争，以期 在战争中取得主动权。因此军用雷达在未来战争中，将工作于各种复杂的电磁环境下，要提 高雷达的生存能力，雷达必须具有对各种恶劣环境的良好适应能力，同时也需要对复杂的 电子战环境具有良好适应能力。
[0003] 电子干扰即电子战中卓为有效的一种方式，它主要是一种通过主动发射电磁波， 或转换发射、反射特定的电磁波信号以破坏雷达对电磁信息的获取和利用的手段。随着雷 达信号处理算法研究的快速发展，干扰侦察技术已经成为现代雷达中重要的一部分。现代 雷达不仅要求信号处理要不断提高目标检测概率，而且还要求在复杂的背景，如杂波、噪 声、干扰下提取出目标信息，并进一步获取更多的目标信息，同时还要求雷达能够根据干扰 回波准确的侦察出干扰特征，如干扰频点、干扰方位、干扰强度等。正确判断出这些信息便 于后面的抗干扰处理，操作员能根据干扰信息更有针对性的采取抗干扰措施，减弱干扰信 号，提尚目标检测概率。
[0004] 传统机械扫描雷达由于受自身体制的影响，通常采取在主通道之外再增加一路辅 助通道来进行干扰侦查检测，用这种方法侦查出的干扰指向线中往往会存在虚假指向线。 而且，由于辅助通道和主通道工作在同一频点下，所以只有当干扰和主通道在同一频点时 才能够侦查到，否则便无法检测出干扰。
【发明内容】
[0005] 本发明的目的在于提出一种适用于机械扫描雷达的有源干扰侦察方法，以解决上 述现有技术的不足，提高干扰的检测准确度。
[0006] 本发明的技术方案是这样实现的：
[0007] -、技术原理
[0008] 本发明在现有机械扫描雷达两个独立的回波接收通道上增加一路干扰侦察通道， 雷达系统采用三个独立的回波接收通道，分别为主天线接收通道、辅助天线接收通道和干 扰侦察通道。其中，主天线接收通道为雷达正常检测通道;干扰侦察通道具有与主天线接收 通道一致的方向性，而且干扰侦察通道通过不断的变换中心频率提供雷达整个工作频带的 回波信息，再根据天线方向图调制特性采用自适应门限算法和质心法，通过对回波信号的 幅度信息分析处理得到干扰源所在方位的干扰指向线及干扰源强度、频点和干扰信号带宽 等信息;辅助天线接收通道没有方向性，且与主辅接收通道工作在同一频点下，辅助天线接 收通道的增益介于主接收天线主瓣增益与主接收天线副瓣增益之间，基于此来消除虚假指 向线，增强侦察结果的准确性和可靠性。
[0009] 二、技术方案
[0010] 基于上述原理，本发明的实现步骤如下：
[0011] (1)获得机械扫描雷达的回波数据
[0012] (la)从雷达的主天线接收通道分出一路作为干扰侦察通道，并根据雷达发射脉冲 的顺序分别接收三路通道中每个脉冲的回波数据，即主天线接收通道回波数据m，辅助天线 接收通道回波数据u，干扰侦察通道原始回波数据g;
[0013] (lb)舍去干扰侦察原始回波数据g中近区回波数据并选取当前工作频点之后的某 一个频点X作为起始频点，再使g从起始位置到结束位置对应的频点依次增加，增加到最大 频点后再从最小频点1开始增加直到频点X-1后结束，得到干扰侦察通道回波数据g'；
[0014] (2)确定干扰侦察通道回波数据g'中干扰回波所在的脉组及功率值：
[0015] (2a)计算干扰侦察通道回波数据g'中每个脉组的频点功率；
[0016] (2b)将雷达系统要检测的最小干扰功率作为干扰侦察门限；
[0017] (2c)将干扰侦察通道回波数据g'中脉组的每个频点功率与干扰侦察门限比较，若 频点功率大于干扰侦察门限，则保存该干扰脉组和频点功率;若小于或等于干扰侦察门限 则丢弃该脉组；
[0018] (2d)对g'中所有的脉组均做（2c)中的处理，当所有的脉组均处理完后，得到所有 的干扰脉组W;
[0019] (2e)判断干扰脉组W是否在所要求的范围之内，若在，则执行(2f);若不在，则根据 干扰脉组中的频点功率重新确定干扰侦察门限，并返回（2c)直至干扰脉组W满足所要求的 范围；
[0020] (2f)计算干扰脉组W中每个脉组收到的最大干扰功率Pi，…，Pw以及Pi，…，Pw对应 的频点fi，…，fV;
[0021] (3)确定干扰强度、干扰频点及带宽信息：
[0022] (3a)比较(2f)中每个脉组收到的最大干扰功率Pi，…，Pw，将其中的最大值作为干 扰强度，该最大值对应的频点作为干扰频点；
[0023] (3b)求出干扰脉组W对应的方位码C1，…，cw;
[0024] (3c)根据(2f)中Pi,…，Pw和(3b)中的ci，…，cw，由质心法求出干扰指向线1;
[0025] (3d)根据雷达系统要侦察的最小干扰设定宽带门限，统计频点功率大于宽带门限 的个数N，当N达到宽带干扰的频带范围时，则判定该干扰为宽带干扰，否则，则认为该干扰 为窄带干扰。
[0026] 本发明与现有技术相比具有如下优点：
[0027] 1、干扰侦察的准确性高
[0028]现有的技术只有当干扰频点与主通道的频点相同时才可以侦察出干扰，而本发明 通过不断改变接收机混频的中心频率使干扰侦察通道提供了雷达整个工作频带范围内所 有频点的回波信息，突破了干扰侦察的局限性，同时本发明还采用动态频点分配及舍去近 区数据的方法避免了地杂波的影响，提高了干扰侦察的准确性；
[0029] 2、侦察结果的可靠性高
[0030] 干扰侦察通道侦察出的干扰指向线可能会存在虚假指向线，基于辅助天线接收通 道的增益介于主接收天线主瓣增益与主接收天线副瓣增益之间，本发明通过比较主辅天线 通道的脉组功率来消除虚假指向线，提高了侦察结果的可靠性。
【附图说明】
[0031] 图1是本发明的实现流程图；
[0032] 图2是本发明中的干扰侦察通道频点分布图；
[0033] 图3是本发明中确定干扰回波所在脉组及功率值的子流程图；
[0034]图4是用本发明对噪声压制式干扰的干扰强度仿真图；
[0035]图5是用本发明对噪声压制式干扰的干扰频点仿真图；
[0036]图6是未消除虚假指向线时噪声压制式干扰的干扰指向线图；
[0037] 图7是消除虚假指向线后噪声压制式干扰的干扰指向线图。
【具体实施方式】
[0038] 为了方便阐述本发明中有源干扰侦察方法及验证本发明的有效性，下面结合具体 实施例及验证结果进一步说明本发明。
[0039]参照图1，本发明的实现步骤如下：
[0040] 步骤1:获得机械扫描雷达的回波数据。
[0041] (la)获取主天线接收通道回波数据m和辅助天线接收通道回波数据u:
[0042] 根据雷达发射脉冲的顺序分别接收主天线接收通道回波数据m、辅助天线接收通 道回波数据u，其中，主天线接收通道提供了雷达正常检测通道的回波信息;辅助天线接收 通道的增益介于主天线接收通道的主瓣增益和第一副瓣增益之间，而且与主天线接收通道 工作在同一频点。由于辅助天线接收通道没有方向性，所以它提供了全方位性的回波信息。 [0043] (lb)获取干扰侦察通道回波数据g' ：
[0044] 从雷达的主天线接收通道分出一路作为干扰侦察通道，根据雷达发射脉冲的顺序 接收干扰侦察通道回波数据g，由于近区常常有地杂波且接收机中频滤波器带外抑制有限， 所以干扰侦察通道回波数据的近区常常会受地杂波的影响，为了避免地杂波对侦察结果的 影响，所以将干扰侦察通道每个脉冲回波数据的前端丢掉而选取后端数据做处理；同时，在 本发明中，对干扰侦察通道除了将近区回波数据丢掉以外，还将当前工作频点与远距离门 数据相对应，得到干扰侦察通道回波数据g' ；
[0045] 如图2所示，本发明以当前工作频点之后的某一个频点X为例进行说明：将X作为起 始频点，从处理的回波数据的起始到回波结束对应的频点依次增加，增加到最大频点后再 从最小频点1开始增加直到频点X-1后结束，通过这种方式使当前工作频点对应后面距离门 的数据。
[0046] 步骤2:确定干扰侦察通道回波数据中干扰回波所在的脉组及功率值。
[0047]参照图3,本步骤的具体实现如下：
[0048] (2a)计算脉组中每个频点的功率：
[0049] 干扰侦察通道中每个脉冲的回波数据都包含了雷达工作频点范围内所有频点的 信息，取出每个频点的数据，对每个频点的数据求模平方和，再将模平方和除以数据个数， 便得到脉冲中每个频点的功率，但由于雷达单个脉冲的回波能量有限，所以常常将多个脉 冲的频点功率做非相参积累后得到脉组中每个频点的功率，并将雷达系统要检测的最小干 扰功率作为干扰侦察门限；
[0050] (2b)将脉组中频点功率和干扰侦察门限做比较：
[0051 ]将脉组中每个频点的功率和干扰侦察门限做比较，若频点功率大于干扰侦察门 限，则保存该频点功率和该脉组;若小于或等于干扰侦察门限则丢弃该脉组，当将所有脉组 的频点功率和干扰侦察门限做完比较后，统计连续收到干扰的脉组W;
[0052] (2c)确定干扰脉组的最大功率值Pi…Pi…Pw和频点f r ? ? f i…:
[0053] (2cl)在干扰侦察通道回波数据g'中提取W个脉组中每个频点的数据；
[0054] (2c2)计算第一个脉组中每个频点的数据模平方和的均值，得到脉组中每个频点 的功率，比较所有频点的功率，得到其中的最大值Pl，最大值丹对应的频点记为f1;再计算第 二个脉组中每个频点的数据模平方和的均值，得到脉组中每个频点的功率，比较所有频点 的功率，得到其中的最大值P2，该最大值? 2对应的频点记为f2，依次对剩余的脉组计算每个 频点的数据模平方和的均值，比较得频点功率的最大值以及该最大值对应的频点，对所有 的脉组均做完处理后，得到PfPi - Pw和斤4^1，其中上表示干扰脉组1中第1个脉组的 最大频点功率，fi表示Pi对应的频点，i的取值范围为1~W，W表示干扰脉组的个数；
[0055] (2d)保存符合条件的干扰脉组W:
[0056] 由于受雷达天线方向图的调制，所以干扰主瓣常常介于一个范围之内，设干扰主 瓣达到最大角度时干扰的脉组个数为K1，干扰主瓣达到最小角度时干扰的脉组个数为K2; [0057] 将干扰脉组W和K1、K2比较，若W小于K2则直接丢掉该干扰脉组W;若W介于K1和K2之 间时，则保存每个脉组频点功率的最大值Pr ? ? Pi…Pw及Pr ? ? Pi…Pw对应的频点f r ? ? f i…f w;若 W大于Kl，则根据Pr-Pr-Pw中的最大值重新确定干扰侦察门限，并返回（2b)。
[0058] 步骤3:利用干扰侦察通道的回波信息确定干扰强度、干扰频点、干扰指向线及干 扰类型。
[0059] (3a)确定干扰频点和干扰强度：
[0060] 比较每个脉组中频点功率的最大值Pr ? ? Pi…Pw，将Pr ? ? Pi…Pw中的最大值作为干扰 强度，该最大值对应的频点即为干扰频点；
[0061] (3b)根据方位码和干扰功率由质心法计算干扰指向线：
[0062]天线扫描一圈为360°，为了提高侦察干扰方位的准确性，将一圈360°划分为K等份 以提高精度，即将K等份依次记为方位码…K-1，每个脉组对应一个方位码，j的取值范 围为0~K-1;
[0063]将步骤2得到的干扰脉组W所对应的方位码依次记为cnw，求干扰指向线1的 方法有:平均方位法、最大方位法和质心法，虽然质心法较前两种计算复杂，但是质心法的 误差较小，所以本发明中采用质心法求干扰指向线1:
[0065] Cl表示干扰脉组W中第i个脉组的方位码;表示干扰脉组W中第i个脉组的频点 功率最大值Pl和第i个脉组的方位码^相乘；
[0066] (3c)根据干扰的频带宽度判断干扰类型：
[0067]根据宽带干扰和窄带干扰频带宽度的差异可以判断干扰为宽带干扰还是窄带干 扰，而干扰频带宽度的差异最直接的反映就是干扰频点个数的多少，当干扰的带宽越宽干 扰的频点个数就会越多。
[0068] 本步骤是根据雷达系统的工作频带设定宽带门限，统计干扰脉组W中频点功率大 于宽带门限的个数N:
[0069] 若N值达到宽带干扰的频带范围，则判定该干扰为宽带干扰；
[0070] 否则，判定该干扰为窄带干扰。
[0071 ]步骤4:消除步骤(3b)中的虚假指向线。
[0072] 由于雷达的工作环境较为复杂，同时受天线方向图副瓣的影响，在步骤3中侦察出 的干扰指向线可能有一些并不是真实的干扰源所在的位置，所以需要消除这些虚假指向 线。
[0073] 同时由于辅助天线接收通道与主天线接收通道具有相同的频率特性并且工作在 同一频点下，所以辅助天线接收通道和主天线接收通道之间有很好的相关性。
[0074] 此外由于辅助天线接收通道的增益介于主天线接收通道主瓣增益和第一副瓣增 益之间，所以若干扰从主天线接收通道的主瓣进入时，主通道接收的干扰功率应该大于辅 助通道接收的干扰功率，则保存该干扰指向线;若干扰从主天线副瓣进入，则主通道接收的 干扰功率会小于辅助通道接收的干扰功率，此时需丢弃该干扰指向线。
[0075]本步骤的具体实现如下：
[0076] (4a)根据在步骤(3b)中已经得到干扰脉组W对应的方位码cni - cw，在主辅通道 中分别取出相对应的脉组数据；
[0077] (4b)分别计算主道中每个脉组的干扰功率PV" PV" Pwm和辅通道中每个脉组的干 扰功率，其中，下标m表示主通道中每个脉组的干扰功率，下标u表示辅助通道 中每个脉组的干扰功率，P1^表示主通道第i个脉组的最大频点功率，表示辅助通道的第i 个脉组最大频点功率；
[0078] (4c)将主通道中脉组的干扰功率PV"PV"Pwm与辅通道中脉组的干扰功率PV" pv"pwu进行比较:若pW^PVPVhPVPV则保存该干扰指向线；否则，认为该干扰指 向线为虚假指向线将其丢弃。
[0079] 下面结合仿真实例叙述本发明方法的正确性。
[0080] 1 ?仿真条件
[00811设定宽带门限为7.0 X108,频点个数为30,一圈方位码共有4000,当有20个频点功 率大于宽带门限时则判断干扰为宽带干扰，否则，判断为窄带干扰。
[0082] 2.仿真内容
[0083]仿真1，仿真噪声压制式干扰的频点和强度。
[0084]提取噪声压制式干扰的回波数据，并按步骤2中的方法确定干扰脉组和干扰频点， 结果如图4所示。由图4可知有8个脉组收到干扰，每个脉组的最大干扰功率为3.8*10~8， 5.8*10~8,7.7*10~8,3.5*10~8,6.5*10~8,3.2*10~8，1.4*10~8，1.7*10~8,其中第 3 个脉组 收到的干扰功率最大为7.7*10~8,所以最大干扰强度为7.7*10~8,最大干扰功率7.7*10~8 对应的频点为5;
[0085]仿真2,仿真噪声压制式干扰的类型。
[0086] 在仿真1中得到了8个脉组收到干扰，再根据步骤(3c)中的方法确定8个干扰脉组 中频点功率大于宽带门限的个数，结果如图5所示。从图5可以看出，只有频点为5的干扰功 率7.7*10~8大于宽带门限，故可知此时只有一个频点的干扰功率大于宽带门限，其数量小 于20,所以判断为窄带干扰。
[0087] 仿真3,仿真噪声压制式干扰的干扰指向线。
[0088] 根据步骤3中的方法得到8个干扰脉组的方位码2329，2330，2331，2332，2333， 2334，2335，2336并根据仿真1中得到8个干扰脉组的干扰功率3.8*10~ 8，5.8*10~ 8，7.7*10~ 8,3.5*10~8,6.5*10~8,3.2*10~8，1.4*10~8，1.7*10~8,通过质心法得到干扰指向线，结果 如图6所示。由图6可知，除了在干扰源方位有指向线外，在195°和106°也有干扰指向线，经 过消除虚假指向线后只在210°有干扰指向线，如图7所示。其中，图6和图7中实线代表主通 道回波数据的干扰分布情况，虚线代表辅助通道回波数据的干扰分布情况。
[0089]综上可知，本发明中的有源干扰侦察方法可以比较准确的侦察出干扰指向、干扰 频点、干扰强度以及干扰类型。实现了对机械扫描雷达有源干扰的侦察，增强了雷达的抗干 扰能力，具有广泛的应用前景。
【主权项】
1. 一种适应于机械扫描雷达的有源干扰侦察方法，包括： (1) 获得机械扫描雷达的回波数据 (la) 从雷达的主天线接收通道分出一路作为干扰侦察通道，并根据雷达发射脉冲的顺 序分别接收三路通道中每个脉冲的回波数据，即主天线接收通道回波数据m，辅助天线接收 通道回波数据u，干扰侦察通道原始回波数据g; (lb) 舍去干扰侦察原始回波数据g中近区回波数据并选取当前工作频点之后的某一个 频点X作为起始频点，再使g从起始位置到结束位置对应的频点依次增加，增加到最大频点 后再从最小频点1开始增加直到频点X-I后结束，得到干扰侦察通道回波数据g' ； (2) 确定干扰侦察通道回波数据g'中干扰回波所在的脉组及功率值： (2a)计算干扰侦察通道回波数据g'中每个脉组的频点功率； (2b)将雷达系统要检测的最小干扰功率作为干扰侦察门限； (2c)将干扰侦察通道回波数据g'中脉组的每个频点功率与干扰侦察门限比较，若频点 功率大于干扰侦察门限，则保存该脉组和该频点功率;若小于或等于干扰侦察门限则丢弃 该脉组； (2d)对g'中所有的脉组均做（2c)中的处理，当所有的脉组均处理完后，得到连续收到 干扰的脉组W; (2e)判断干扰脉组W是否在所要求的范围之内，若在，则执行(2f);若不在，则根据干扰 脉组中的频点功率重新确定干扰侦察门限，并返回（2c)直至干扰脉组W满足所要求的范围； (2f)计算干扰脉组W中每个脉组收到的最大干扰功率Pr^Pi - Pw以及Pr^Pi - Pw对应的 频点fr ·· fi…f w，其中，Pi表示干扰脉组W中第i个脉组的最大频点功率，fi表示Pi对应的频 点，i的取值范围为1~w，w表示干扰脉组的个数； (3) 确定干扰强度、干扰频点、干扰指向线及干扰类型： (3a)比较(2f)中每个脉组收到的最大干扰功率Pr^Pi…Pw，将其中的最大值作为干扰 强度，该最大值对应的频点作为干扰频点； (3b)求出干扰脉组W对应的方位码(^???(^???(^。其中，(^表示干扰脉组1中第1个脉组的方 位码 (3c)根据(2f)中Ρι···Ργ··Ρ?和(3b)中的cr"cr"cw，由质心法求出干扰指向线1; (3d)根据雷达系统的工作带宽设定宽带门限，统计频点功率大于宽带门限的个数N，当 N达到宽带干扰的频带范围时，则判定该干扰为宽带干扰，否则，则认为该干扰为窄带干扰； (4) 消除(3c)中的虚假指向线： (4a)在主通道回波数据m和辅助通道回波数据u中分别取出与干扰侦察通道回波数据 g'中干扰脉组W相对应的脉组数据； (4b)分别计算主通道回波数据m和辅助通道回波数据u中W个脉组中每个脉组的干扰功 率，分别用PWvPwm和ρν··ρν··ρ\·表示，其中，Pim表示主通道第i个脉组的最大频点功 率，P 1u表示辅助通道的第i个脉组最大频点功率； (4c)比较PV"PV"Pwm与 扰指向线1，否则，丢弃(3c)中的干扰指向线1。2. 根据权利要求1适应于机械扫描雷达的有源干扰侦察方法，其特征在于，（2f)中Pr·· Pr"Pw和fr"fr"fw的计算过程如下： (2Π )在干扰侦察通道回波数据g'中提取W个脉组中每个频点的数据； (2f2)计算第一个脉组中每个频点的数据模平方和的均值，得到脉组中每个频点的功 率，比较所有频点的功率，得到其中的最大值P1，最大值应的频点记为h;再计算第二个 脉组中每个频点的数据模平方和的均值，得到脉组中每个频点的功率，比较所有频点的功 率，得到其中的最大值P 2，该最大值P2对应的频点记为f2，依次对剩余的脉组计算每个频点 的数据模平方和的均值，比较得频点功率的最大值以及该最大值对应的频点，对所有的脉 组均做完处理后，得到Pr"Pr"Pw和fr"fV"f w。3. 根据权利要求1适应于机械扫描雷达的有源干扰侦察方法，其特征在于，（3b)中计算 方位码(^，…，⑶，是先用天线扫描一圈即360°，再将一圈360°划分为K等份，将K等份依次记 为方位码0…j…K-I，j的取值范围为0~K-I，每个脉组都对应一个方位码，将W个脉组对应 的方位码记为C1· **Ci***Cw〇4. 根据权利要求1适应于机械扫描雷达的有源干扰侦察方法，其特征在于，（3c)中由质 心法求出干扰指向线1，按加下公式计笪，其中，Pr^Pi - Pw是W个脉组中每个脉组收到的最大干扰功率，cr^ci…Cw是W个脉组对 应的方位码。
【文档编号】G01S7/41GK105891799SQ201610473461
【公开日】2016年8月24日
【申请日】2016年6月24日
【发明人】何学辉, 翟小娟, 杨琨, 赵永波, 苏涛, 刘宏伟
【申请人】西安电子科技大学