基于波束保形的阵列雷达自适应抗干扰方法
【专利摘要】本发明公开了一种基于波束保形的阵列雷达自适应抗干扰方法，主要解决现有技术形成的波束主瓣不平坦以及不适用于多点约束的问题。其实现步骤是：1.采集阵列雷达检测目标的采样数据；2.设计优化的三个约束点的导向矢量；3.根据步骤2中的导向矢量构建数学模型；4.设计数学模型的约束矢量，并根据该约束矢量求解数学模型的自适应权矢量；5.用自适应权矢量对干扰信号进行抑制，得到输出数据。本发明通过在主波束多个点内采用线性约束，通过设置相位响应的值实现对目标信号相位的匹配，避免了干扰对检测性能的影响；同时通过设置主波束幅度响应为1，使主波束能均匀覆盖观测区域，提高了阵列雷达的目标检测性能，可用于多目标的检测。
【专利说明】
基于波束保形的阵列雷达自适应抗干扰方法
技术领域
[0001 ]本发明属于雷达技术领域，更进一步涉及阵列自适应波束形成方法，可用于目标 检测。
【背景技术】
[0002] 阵列自适应波束形成方法被广泛应用于许多领域，例如雷达，声呐，无线通信，医 疗成像等等。其中，最小方差无失真响应MVDR波束形成器是最典型的一种。传统设计波束形 成器的方法都假定感兴趣信号SOI的导向矢量是准确已知的。但是，在实际应用中，目标信 号污染训练样本、样本数不足、期望信号导向矢量误差是导致目标检测性能下降的主要原 因。传统的MVDR波束形成器对这些误差十分敏感。因此，实际应用急需稳健的自适应波束形 成方法来克服上述因素造成目标检测性能下降的问题。
[0003] 目前，许多学者提出了一些稳健的自适应波束形成方法，其中，
[0004] 第一种是：线性约束最小方差LCMV方法，该方法通过在目标附近的多个点上进行 线性约束，使其幅度响应均为1，同时最小化阵列的输出总功率，以在抑制干扰的同时展宽 波束的主瓣。但是，该方法由于没有考虑相位响应，因而使得波束的主瓣区域不平坦，造成 目标检测性能的下降；
[0005] 第二种是：两点的二次约束波束形成方法，该方法通过在靠近目标的两点上进行 二次约束，使其幅度响应模值的平方为1，同时最小化阵列的输出总功率，以提高对目标到 达角估计误差的稳健性，它可以看作是一种响应矢量进一步优化的LCMV方法。但是，该方法 只适用于临近两点的约束，不适用于多点约束的情况，当干扰在约束点以外时会严重影响 目标的检测性能。
【发明内容】
[0006] 本发明的目的在于针对上述已有技术的不足，提出一种基于波束保形的阵列雷达 自适应抗干扰方法，以提高阵列雷达的目标检测性能。
[0007] 本发明的基本思路是:通过在主波束多个点内采用线性约束，通过设置相位响应 的值实现对目标信号相位的匹配，有效实现主瓣波束的保形，通过设置主波束幅度响应为 1，使得主波束均匀覆盖观测区域，其实现方案包括如下：
[0008] 1)在阵列信号处理机中输入阵列雷达检测目标的采样数据X，该数据X包含目标和 干扰信号；
[0009] 2)设计优化的三个约束点的导向矢量：
[0010] 2a)设计三个约束点的角度，即第一个约束点的角度0:是估计的目标到达角，第二 个约束点的角度02 = 1-3°，第三个约束点的角度03 = 003° ;
[0011] 2b)根据步骤2a)中三个约束点的角度0:、02、03，计算优化的三个约束点的导向矢 量分别为a(9i)，a(%)，a(03);
[0012] 3)根据步骤2b)中三个约束点的导向矢量apihapshapj构建数学模型，即 min \v " Rw 使得WHC = uT， w
[0013]其中，w为待求的自适应权矢量，（?）H表示共辄转置操作，（?）T表示转置操作，皂 为采样数据的协方差矩阵，CztaOd a(02) a(03)]为约束矩阵，= 为约束 矢量，？>2、％为待求的两个不同约束相位，j表示虚数；
[0014] 4)设计约束矢量：
[0015] 4a)根据步骤3)中的约束矩阵C和采样协方差矩阵食算得到一个中间矩阵 an 〇126^2 a[3ejt,}
[0016] (C并R/C)-、512e-九 a22 a23eAs aueAi a2ie,bl' a i3
[0017] 其中，（?）-1 表示求逆操作，31:1、.兔.2.、a13、心、323s 5 12.、..石3、.. 523为上述 中间矩阵中的九个数值不同的中间数；
[0018] 4b)根据步骤4a)中的第七个中间数&和第九个中间数1得到第二个约束点的约 束相位熟==7r-Al2和第三个约束点的约束相位_=及-523+炉2=2;rH ;
[0019] 4c)根据步骤4b)中的第二个约束点的约束相位於和第三个约束点的约束相位奶 得到约束矢量为《= 1 e/Ml:) '
[0020] 5)求解自适应权矢量：
[0021 ]根据步骤3)中的约束矩阵C以及步骤4c)中的约束矢量u得到自适应权矢量为： [0022] ；
[0023] 6)根据步骤5)中的自适应权矢量w对干扰信号进行抑制，得到抑制干扰后的输出 信号Y=w TX，其中，（?）T表示转置操作。
[0024] 本发明与现有技术相比具有以下优点：
[0025] 本发明通过在主波束多个点内采用线性约束，通过设置相位响应的值实现对目标 信号相位的匹配，有效实现主瓣波束的保形，同时通过设置主波束幅度响应为1，使得主波 束均匀覆盖观测区域，避免了干扰对检测性能的影响，实现了阵列自适应波束形成方法的 稳健性，提高了阵列雷达的目标检测性能，可用于多目标的检测。
【附图说明】
[0026] 图1是本发明的使用场景图；
[0027] 图2是本发明的实现流程图；
[0028] 图3是用发明方法检测目标的波束形成图；
[0029]图4是用本发明方法检测目标的性能改善因子与最优性能改善因子的对比曲线 图。
【具体实施方式】
[0030] 下面结合附图对本发明实施例及效果作进一步的详细描述。
[0031] 参照图1，本发明的使用场景中，采用阵列平台作为雷达平台，雷达配置方式为一 维等距线阵，阵元数为M，阵元间距d为半波长，信号波达角为0。
[0032]参照图2,本发明的实现步骤如下：
[0033] 步骤1，采集阵列雷达检测目标的采样数据X。
[0034] la)阵列雷达在目标运动的过程中实时采集检测目标的采样数据X，该数据包含目 标和干扰信号；
[0035] lb)将采样数据X输入到阵列信号处理机中。
[0036]步骤2,在阵列信号处理机中设计优化的三个约束点的导向矢量。
[0037] 2a)设计三个约束点的角度
[0038] 设第一个约束点的角度0:是估计的目标到达角，第二个约束点的角度02 = 0^3°， 第三个约束点的角度03 = 003° ;
[0039] 当干扰靠近目标信号时，阵列雷达检测目标的性能会严重下降，所以在目标信号 及其附近两点进行约束，从而可以使波束得到保形。
[0040] 2b)计算优化的三个约束点的导向矢量：
[0041]根据步骤2a)中的0L0L03得到优化的三个约束点的导向矢量分别为8(00^(02)， a(03)即：
[0042] 第一个约束点的导向矢量：
[0043] 第二个约束点的导向矢量：
[0044] 第三个约束点的导向矢量：
[0045]其中，j表示虚数，d为阵元间距，A为发射信号的波长，M为阵列阵元数，（?）T表示 转置操作。
[0046] 步骤3,构建数学模型。
[0047] 约束目标及其附近两点的幅度及相位响应，使它们跟目标信号相匹配，同时最小 化阵列的输出总功率，在抑制干扰的同时实现波束主瓣的保形。
[0048] 3a)计算采样数据的协方差矩阵Ry
[0050] 其中，L为采样的快拍数，x(i)表示阵列雷达第i次快拍的接收数据；
[0051] 3b)根据步骤2b)中的三个约束点的导向矢量计算约束矩阵C =
[aPO a(02) a(03)]和约束矢量《=〇 其中，约、热为待求的约束相位， (? )T表示转置操作；
[0052] 3c)根据步骤3a)中的采样数据的协方差矩阵食dP3b)中的约束矩阵C和约束矢量u 得到如下数学模型：
[0053] mjn w良W.使得使wHc = uT，
[0054] 其中，w为待求的自适应权矢量，（?）H表示共辄转置操作。
[0055] 步骤4,设计约束矢量。
[0056] 4a)根据步骤3a)中的采样协方差矩阵和步骤3b)中的约束矩阵C计算得到一个 中间矩阵
[0058] 其中，（?）-1 表示求逆操作，5lt、.5I2.、513、tt、t},;、lj 2、. Ij3、. 为上述 中间矩阵中的九个数值不同的中间数；
[0059] 4b)根据步骤4a)中的第七个中间数&和第九个中间数&得到第二个约束点的约 束相位=贫-^! 2:和第三个约束点的约束相位約=贫-^23.+=.2_；r-/)12
[0060] 4c)根据步骤4b)中第二个约束点的约束相位约和第三个约束点的约束相位鈣得 到约束矢量为：
[0061] ?: = [1 々）/"釗|,
[0062] 该约束矢量中1表示主波束幅度响应，使主波束均匀覆盖观测区域，并且其相位响 应与目标信号相匹配，因此在检测目标的性能上损失较小。
[0063]步骤5,求解自适应权矢量。
[0064] 5a)定义步骤3中数学模型的拉格朗日函数+ 其中，y为拉 格朗日乘子，拉格朗日函数1是关于最优权矢量w的函数；
[0065] 5b)对上述拉格朗日函数1进行最小化处理，利用该函数求得步骤3中数学模型的 角军w，即为自适应权矢量：
[0066] w =旋/史(Cf't-尤)-1 M。
[0067]其中，C为约束矩阵，u为约束矢量。
[0068]步骤6,获得输出数据。
[0069]信号处理机对阵列雷达检测目标的采样数据X进行处理，即根据步骤5中自适应权 矢量w对干扰信号进行抑制，得到抑制干扰后信号处理机的输出数据:Y=wtX，
[0070]其中，（?）T表示转置操作，此时Y中只包含目标信号，不再含干扰信号，完成对目 标信号的检测。
[0071 ]下面通过仿真实验对本发明的效果做进一步说明。
[0072] 1 ?仿真参数：
[0073]设阵列雷达为阵元间距为半波长的等距线阵，阵元数M= 12，脉冲数K= 100，感兴 趣信号的到达角为5°，波束照射方向为8°，两个干扰的方位角为-25°和30°，三个约束点的 方位角分别为5°、8°、13°，信号噪声比为10dB，杂波噪声比为40dB。
[0074] 2.仿真内容：
[0075] 仿真1，在上述仿真参数下，用发明方法进行稳健的波束形成，结果如图3所示。
[0076] 由图3可见，本发明方法的波束形成方向图在主瓣区域非常平坦，没有发生畸变， 证明了本发明的有效性。
[0077]仿真2,在上述仿真参数下，用本发明方法检测目标的阵列输出信干噪比随输入信 噪比变化曲线，结果如图4所示。
[0078]由图4可见，本发明方法的改善因子只比最优性能下降3dB左右，可见本发明方法 的性能较好。
【主权项】
1. 基于波束保形的阵列雷达自适应抗干扰方法，包括： 1) 在阵列信号处理机中输入阵列雷达检测目标的采样数据X，该数据X包含目标和干扰 信号； 2) 设计优化的三个约束点的导向矢量： 2a)设计三个约束点的角度，即第一个约束点的角度Θ:是估计的目标到达角，第二个约 束点的角度02 = 0^3°，第三个约束点的角度θ3 = θ1+3° ; 2b)根据步骤2a)中三个约束点的角度θ:、θ2、θ3，计算优化的三个约束点的导向矢量分 别为a(0〇，a(02)，a(03); 3) 根据步骤2b)中三个约束点的导向矢量构建数学模型，即其中，w为待求的自适应权矢量，（·）H表示共辄转置操作，（·）T表示转置操作，食t为采 样数据的协方差矩阵，C= [3(0:)3(02)8(03)]为约束矩阵，= 1 e*]为约束矢量， %、於为待求的两个不同约束相位，j表示虚数； 4) 设计约束矢量： 4a)根据步骤3)中的约束矩阵C和采样协方差矩阵算得到一个中间矩阵其中，（·）-1表示求逆操作，3U、"'r、成3、矣 2.、.523、、吞3、石23为上述中间 矩阵中的九个数值不同的中间数； 4b)根据步骤4a)中的第七个中间数各12和第九个中间数&23得到第二个约束点的约束相 位灼=；T-S12和第三个约束点的约束相位约=疋-43 +? 表3 ; 4c)根据步骤4b)中的第二个约束点的约束相位於和第三个约束点的约束相位朽得到5) 求解自适应权矢量： 根据步骤3)中的约束矩阵C以及步骤4c)中的约束矢量u得到自适应权矢量为：6) 根据步骤5)中的自适应权矢量w对干扰信号进行抑制，得到抑制干扰后的输出信号Y =wTX，其中，（· "表示转置操作。2. 根据权利要求1所述的方法，其中步骤2b)中计算优化的三个约束点的导向矢量a (00,8(02)4(03),按如下公式进行：其中，j表示虚数，d为阵元间距，λ为发射信号的波长，Μ为阵列阵元数，（·）τ表示转置操 作。3.根据权利要求1所述的方法，其中，步骤3)中的采样协方差矩阵<，按如下公式计算：其中，L为采样的快拍数，x(i)表示阵列雷达第i次快拍的接收数据。
【文档编号】G01S7/41GK105929374SQ201610471771
【公开日】2016年9月7日
【申请日】2016年6月24日
【发明人】许京伟, 廖桂生, 冯阳, 王成浩, 兰岚
【申请人】西安电子科技大学