1.一种波束域补偿的匀加速机载雷达杂波抑制方法，其特征在于，具体按照以下步骤实施：步骤1，确定机载雷达，且机载雷达检测范围内存在目标，设机载雷达向其检测范围内发射信号并经目标反射后接收到的回波信号，记为雷达原始回波数据矩阵X，雷达原始回波数据矩阵X为N×K×L维的阵元‑脉冲‑距离域数据，N表示机载雷达天线的阵元数，K表示机载雷达在一个相干积累时间内发射的脉冲数，L表示机载雷达接收到的回波数据包含的距离单元总个数；

步骤2，根据阵元加权FFT的切比雪夫加权值，确定阵元FFT加权系数矢量β；

步骤3，令l＝1，l∈{1,2,…,L}，l表示第l个距离单元；

步骤4，对第l个距离单元的N×K维阵元‑脉冲维数据Xl，根据机载雷达的距离分辨率，确定第l个距离单元的斜距Rl，根据第l个距离单元的斜距Rl，确定第l个距离单元上的散射点的俯仰角 根据第l个距离单元上的散射点的俯仰角 和多波束形成的角度范围，确定对第l个距离单元的N×K维阵元‑脉冲维数据Xl进行阵元加权FFT的Q个归一化空间频率值，Q的大小和机载雷达天线的阵元数N相等；

步骤5，根据步骤4得到的归一化空间频率值 和步骤2得到的阵元FFT加权系数矢量β，确定第l个距离单元的N×K维阵元‑脉冲维数据的阵元FFT加权变换矩阵Fl；使用阵元FFT加权变换矩阵Fl对所述第l个距离单元的阵元‑脉冲维数据Xl进行阵元域加权FFT变换，得到第l个距离单元波束‑脉冲维回波数据Yl；

步骤6，令q＝1，q∈{1,2,…,Q}，q表示第q个波束；

步骤7，根据机载雷达飞行的加速度a及其它雷达参数，确定对第l个距离单元、第q个波束的K×1维杂波数据yql的K×1维补偿矢量tql，并根据补偿矢量tql对第l个距离单元、第q个波束的脉冲维回波数据yql进行补偿，得到补偿后的脉冲维数据步骤8，令q的值加1，并重复执行步骤8，得到第l个距离单元的所有Q个波束的补偿后的脉冲维数据步骤9，令l的值加1，并重复执行步骤5到步骤9，得到所有L个距离单元的所有Q个波束的补偿后的脉冲维数据 从而根据所有补偿后的脉冲维数据 构建波束域补偿后的N×K×L维的波束‑脉冲‑距离域数据步骤10，对步骤9构建的波束域补偿后波束‑脉冲‑距离域数据 在脉冲域进行FFT变换，得到波束‑多普勒‑距离域数据步骤11，对步骤10得到的波束‑多普勒‑距离域数据 进行空时降维自适应处理，得到空时降维自适应处理后的剩余杂波数据Z。

2.根据权利要求1所述的一种波束域补偿的匀加速机载雷达杂波抑制方法，其特征在于，所述步骤2中，阵元FFT加权系数矢量β的表达式为：T

β＝[β0 β1 … βN‑1] (1)

T

式(1)中，βn,n＝1,…,N表示阵元FFT加权系数的每一项，(·) 表示转置操作。

3.根据权利要求1所述的一种波束域补偿的匀加速机载雷达杂波抑制方法，其特征在于，步骤4的具体过程为：步骤4.1，对第l个距离单元的N×K维阵元‑脉冲维数据Xl，根据机载雷达的距离分辨率，确定第l个距离单元的斜距Rl；

机载雷达的距离分辨率的表达式为：

式(2)中，c表示光速，B表示机载雷达发射信号的带宽；

则第l个距离单元的斜距Rl的表达式为：

Rl＝(l‑1)ΔR (3)

步骤4.2，根据第l个距离单元的斜距Rl，确定第l个距离单元上的散射点的俯仰角第l个距离单元上的散射点的俯仰角的表达式为：‑1

式(4)中，sin (·)表示反正弦函数，H表示载机飞行高度，Re＝4/3re表示地球等效半径，re为地球真实半径；

步骤4.3，根据第l个距离单元上的散射点的俯仰角 和多波束形成的角度范围，确定对第l个距离单元的N×K维阵元‑脉冲维数据Xl进行阵元加权FFT的Q个归一化空间频率值，Q的大小和机载雷达天线的阵元数N相等；

假设多波束形成的角度范围为[π‑θmax,θmax]，其中θmax为形成多波束的波束中心指向的最大方位角；

对第l个距离单元的N×K维阵元‑脉冲维数据Xl进行阵元加权FFT的Q个归一化空间频率值的表达式为：

4.根据权利要求1所述的一种波束域补偿的匀加速机载雷达杂波抑制方法，其特征在于，步骤5的具体过程为：步骤5.1，根据步骤4得到的归一化空间频率值 和步骤2得到的阵元FFT加权系数矢量β，确定第l个距离单元的N×K维阵元‑脉冲维数据的阵元FFT加权变换矩阵Fl，其表达式为：步骤5.2，使用阵元FFT加权变换矩阵Fl对第l个距离单元的阵元‑脉冲维数据Xl进行阵元域加权FFT变换，得到第l个距离单元波束‑脉冲维回波数据Yl，表达式为：H

式(7)中，(·) 表示共轭转置操作。

5.根据权利要求1所述的一种波束域补偿的匀加速机载雷达杂波抑制方法，其特征在于，步骤7的具体过程为：步骤7.1，根据机载雷达飞行的加速度a及其它雷达参数，得到散射点的时域导向矢量bq，通过散射点的时域导向矢量bq，确定对第l个距离单元、第q个波束的K×1维杂波数据yql的K×1维补偿矢量tql；

其中，散射点的时域导向矢量bq的表达式为：

bq＝bv,q⊙ba,q (8)

式(8)中，⊙表示Hadamard积；

其中，

式(9)‑(11)中，v0表示载机飞行的初始速度，a表示载机飞行的加速度，λc表示机载雷达发射信号的波长，Tr表示机载雷达发射信号的脉冲重复周期，α表示机载雷达天线的安装角，θq表示第q个波束的波束中心指向的方位角；

对俯仰角为 方位角为θq的散射点的时域导向矢量bq中的二次项增量进行补偿的K×

1维补偿矢量tql

步骤7.2，根据补偿矢量tql对第l个距离单元、第q个波束的脉冲维回波数据yql进行补偿，得到补偿后的脉冲维数据 表达式为：

6.根据权利要求1所述的一种波束域补偿的匀加速机载雷达杂波抑制方法，其特征在于，步骤9的具体过程为：步骤9.1，根据第l个距离单元的所有Q个波束的补偿后的脉冲维数据 得到第l个距离单元的补偿后的波束‑脉冲维数据

步骤9.2，根据第l个距离单元的补偿后的波束‑脉冲维数据 可以构建波束域补偿后的Q×K×L维的波束‑脉冲‑距离域数据 数据 为数据 在距离维的堆叠。

7.根据权利要求1所述的一种波束域补偿的匀加速机载雷达杂波抑制方法，其特征在于，步骤10的具体过程为：对第l个距离单元的波束域补偿后波束‑脉冲维数据 进行脉冲域加权FFT变换，可以得到第l个距离单元波束‑多普勒域回波数据 为其中，脉冲域FFT变换矩阵为：

式(16)中，M是多普勒通道数目，K是相干积累脉冲数，ai,i＝0,…,K‑1表示FFT加权系数，表示第m个多普勒通道的归一化多普勒频率，表示为表示第M个多普勒通道的归一化多普勒频率；

对所有距离单元的波束‑多普勒域回波数据 均进行脉冲域加权FFT变换，可以得到波束‑多普勒‑距离域数据

8.根据权利要求1所述的一种波束域补偿的匀加速机载雷达杂波抑制方法，其特征在于，步骤11的具体过程为：步骤11.1，构造第m个多普勒通道的空时降维矩阵，即：式(17)中，Bm＝[Γm‑1 Γm Γm+1]表示时域降维矩阵；IN表示N×N维的单位阵，是空域降维矩阵；Γm表示脉冲域FFT变换矩阵Γ的第m列；

对于第l个距离单元的波束‑多普勒域回波数据 可将其整理为NK×1维矢量，即：用邻近第l个距离单元的P个训练单元的数据来估计第l个距离单元的杂波加噪声协方差矩阵，即：式(19)中， 表示针对第l个距离单元的P个训练单元的数据中的第p个训练单元的数据；

步骤11.2，构造目标导向矢量，即：

式(20)中， 表示Kronecker积， 为目标回波的归一化多普勒频率， 为目标回波的空间频率， 表示目标空域导向矢量， 为目标时域导向矢量，分别表示为：式(21)‑(24)中，v表示载机飞行速度，θ0为目标方位角， 为目标俯仰角，vt为目标径向速度，λc为雷达发射信号波长，fr为脉冲重复频率，d为阵元间距；

步骤11.3，对于第l个距离单元，第m个多普勒通道的数据可以进行空时降维自适应处理，得到降维后的搜索空时导向矢量，即：降维后的第l个距离单元，第m个多普勒通道的回波数据为：降维后的杂波加噪声协方差矩阵为：

式(27)中， 表示针对第l个距离单元的P个训练单元中第p个训练单元、第m个多普勒通道的回波数据；

在保证目标信号增益不变的前提下，根据线性约束最小准则，计算第l个距离单元，第m个多普勒通道的空时自适应滤波器，表达式为：对第l个距离单元，第m个多普勒通道的回波数据 可以降维自适应滤波输出为对波束‑多普勒‑距离域数据 的所有距离单元和多普勒通道数据进行空时降维自适应处理，得到空时降维自适应处理后的剩余杂波数据Z。