CN112904298B - 一种基于局部网格分裂的网格偏离空时自适应处理方法 - Google Patents

Abstract

本发明公开了一种基于局部网格分裂的网格偏离空时自适应处理方法，包括：获取杂波空间的空时自适应处理字典及杂波协方差矩阵；初始化杂波空间的正交投影矩阵和噪声投影矩阵，初始化空时导向矢量集合、杂波原子下标集合、全局迭代次数、全局迭代门限、局部迭代次数、局部迭代的误差门限、局部的最大迭代次数；在所述空时自适应处理字典中搜索杂波的空时导向矢量，将其构成更新后的杂波空间的空时导向矢量集合，并计算更新后的杂波空间的噪声投影矩阵；根据更新后的杂波空间的噪声投影矩阵获得网格失配的空时自适应处理加权矢量。该方法克服了估计精度与局部网格划分的个数有关的问题，减小了运算量，并具有更优的处理效率。

Description

一种基于局部网格分裂的网格偏离空时自适应处理方法

技术领域

本发明属于数据处理技术领域，具体涉及一种基于局部网格分裂的网格偏离空时自适应处理方法，可用于广域预警雷达。

背景技术

空时自适应处理(Space-time adaptive processing,STAP)是一种检测慢速目标的有效方法。然而，在极度非均匀杂波条件下，由于距离空变性的原因，采用相邻距离门的杂波单元作为训练样本估计杂波空间的方法将不再适用。近年来，常采用稀疏恢复算法来重建杂波的角度-多普勒谱，以获得与最优处理器接近的处理性能，即将空时二维平面离散为网格。然而，由于真实杂波点并不落于网格的中心，即造成网格偏离，导致杂波抑制的性能将会下降。同时，系统误差同样会使稀疏STAP的性能下降。

对于网格偏离问题，通常用划分网格的方法来解决。但是，过密的划分网格将会导致计算量增加，同时将会提高相邻杂波原子间的相关性。为了解决上述问题，基于稀疏贝叶斯学习的方法被提出来。然而，这类方法仅仅解决了在空域频率的网格偏离问题，不能解决在空时模型下二维频域内的网格偏离问题。

目前，解决二维频域内的网格偏离问题方法主要有以下三种：

第一种是利用连续字典的原子范数最小化方法。该方法利用原子范数的性质，把最小化原子范数的问题转化为对半正定问题的求解。然而在两维频域，该方法所用到的块托普利茨矩阵的范德蒙德分解性质并未被严格证明，仅在特定的情况下可以获得近似的解，故其实用性较差。同时，由于该方法占用了大量的内存与计算量，因而在实际的工程应用中受到了限制。

第二种是通过凸优化的方法来解决此问题。Teke O等人中提出了一种扰动正交匹配追赶法，该方法通过在每次迭代中调整所选的杂波原子来不断的降低误差以达到匹配所有杂波原子的目的。然而，该方法计算量大，并且倾向于解决随机扰动的问题，不具备普遍性。

第三种是通过在全局匹配搜索到的杂波原子周围划分更小的格子，该方法相比于前两种更具一般性。G.Bai，R.Tao等人提出了一种参数搜索的正交匹配追赶方法，该方法采用了最速下降法来进行匹配搜索，这可能导致一个非收敛的结果。李等人提出了一种基于子空间投影的降维局部搜索估计方法(Local search clutter subspace estimationSTAP，LSCSE-STAP)，该方式通过设置搜索门限来降低计算量，并且具有固定的网格划分方式。但是，该方法在局部搜索过程中仍然需要大量的计算，且在搜索最佳精度时需要进行数次实验对局部网格的大小进行合理划分，因而效率较低。

发明内容

为了解决现有技术中存在的上述问题，本发明提供了一种基于局部网格分裂的网格偏离空时自适应处理方法。本发明要解决的技术问题通过以下技术方案实现：

本发明提供了一种基于局部网格分裂的网格偏离空时自适应处理方法，包括：

步骤1：获取杂波空间的空时自适应处理字典及杂波协方差矩阵；

步骤2：初始化杂波空间的正交投影矩阵

和噪声投影矩阵

初始化空时导向矢量集合

杂波原子下标集合Ω、全局迭代次数k、全局迭代门限η、局部迭代次数p、局部迭代的误差门限σ_l、局部的最大迭代次数；

步骤3：在所述空时自适应处理字典中搜索杂波的空时导向矢量，构成更新后的杂波空间的空时导向矢量集合Φ_Ω，并计算更新后的杂波空间的噪声投影矩阵P_n；

步骤4：根据更新后的杂波空间的噪声投影矩阵P_n获得网格失配的空时自适应处理加权矢量。

在本发明的一个实施例中，所述步骤1包括：

利用机载雷达系统观测空间中的杂波散射点，计算第i个杂波散射点的空时导向矢量s(f_d,i,f_s,i)，并根据所有杂波散射点的空时导向矢量s(f_d,i,f_s,i)构造杂波空间的空时自适应处理字典Φ：

其中，第i个杂波散射点的空时导向矢量s(f_d,i,f_s,i)为：

其中，

表示Kronecker积，s_d表示杂波散射点的时域导向矢量，s_s表示杂波散射点的空域导向矢量：

s_d(f_d,i)＝[1,exp(j2πf_d,i),...,exp(j2π(M-1)f_d,i)]^T

s_s(f_s,i)＝[1,exp(j2πf_s,i),...,exp(j2π(N-1)f_s,i)]^T

其中，(·)^T表示矩阵的转置操作，N_s表示利用稀疏恢复算法后空时二维平面空域频率被离散为N_s＝ρ_sN个网格，N_d表示利用稀疏恢复算法后时域频率被离散为N_d＝ρ_dM个网格，ρ_d,ρ_s＞1，表示频域的离散程度，N为阵元个数，M为脉冲个数，f_s,i表示第i个杂波散射点的空域频率，f_d,i表示第i杂波散射点的归一化的多普勒频率；

根据杂波散射点的回波数据x_l构造杂波的协方差矩阵

其中，l表示第l个距离门。

在本发明的一个实施例中，所述步骤3包括：

步骤3.1：在所述空时自适应处理字典中进行全局搜索，通过子空间投影的方式挑选出与当前剩余杂波最相关的空时导向矢量作为进行局部迭代的初始杂波原子的空时二维平面坐标；

步骤3.2：以步骤3.1中挑选的原始杂波原子为中心进行多次网格分裂，直至杂波原子空时权矢量的响应满足局部迭代误差门限，对应的杂波原子为本次全局搜索到的杂波原子，其下标记为

步骤3.3：对搜索到的杂波原子下标集合Ω进行扩充

对搜索到的杂波原子的空时导向矢量集合

进行更新

对

进行更新杂波空间的正交投影矩阵

对杂波空间的噪声投影矩阵

进行更新

其中I为单位阵，

为矩阵的求伪逆操作，k表示全局迭代的次数；

步骤3.4：令k＝k+1，重复步骤3.1至步骤3.3，直到剩余杂波能量达到全局搜索的门限，获得此时更新的噪声投影矩阵

在本发明的一个实施例中，所述步骤3.1包括：

步骤3.11：通过子空间投影的方式，在所述空时自适应处理字典中挑选出与当前剩余杂波最相关的空时导向矢量，

表示第k次全局搜索到的原子下标，

的表达式为：

其中，m表示空时自适应字典中杂波原子的下标，s_m表示下标为m的杂波原子所对应的空时导向矢量；

步骤3.12：挑选

的表达式中惩罚项最大的一个，作为在当前的杂波子空间中搜索到的杂波原子。

在本发明的一个实施例中，所述步骤3.2包括：

步骤3.21：以步骤3.1中挑选的原始杂波原子为中心，分别在时域以1/M、在空域以1/N为边长划分一个2×2的“田”字网格，所述网格的9个顶点构成3×3的原子矩阵Φ_Ωl，计算9个杂波原子的空时权矢量的响应：

其中，p表示局部分裂的次数，l表示第l个距离门，

表示第p-1次局部分裂后更新的空时导向矢量集合里的空时导向矢量，

表示第p-1次局部分裂后更新的原子所组成的空时导向矢量集合，

表示第p次局部分裂后下标为m_l的杂波原子的空时权矢量的响应；

步骤3.22：利用所述9个杂波原子中在空时权矢量的最大响应点作为顶点，向第二与第三大响应的点做边，闭合成为下一次分裂的次网格；在次网格内重新划分2×2的网格并更新，重复步骤3.21的过程直到

其中，σ_l表示误差门限，或者达到局部的最大迭代次数；

步骤3.23：通过最后一次局部迭代挑选出最优杂波原子

在本发明的一个实施例中，步骤3.4包括：

令k＝k+1，重复步骤3.1至步骤3.3，直到终止条件成立，得到最终更新的噪声投影矩阵

所述终止条件为：

E^(k)＜ηE⁽⁰⁾，

其中，E⁽⁰⁾表示所有杂波原子的能量，E^(k)表示剩余杂波原子的能量，η为预设的门限。

在本发明的一个实施例中，所述S4包括：

根据得到的更新后的噪声投影矩阵

获得网格失配的空时自适应处理加权矢量为：

其中，s₀表示目标的空时导向矢量。

与现有技术相比，本发明的有益效果在于：

1、本发明基于局部网格分裂的网格偏离空时自适应处理方法克服了估计精度与局部网格划分的个数有关的问题，减小了运算量，节约了运算资源，并具有更优的处理效率。

2、本发明的方法提高了杂波原子的搜索精度，输出信杂噪比曲线凹口变窄，有效地提升了稀疏空时处理的性能。

以下将结合附图及实施例对本发明做进一步详细说明。

附图说明

图1是本发明实施例提供的一种基于局部网格分裂的网格偏离空时自适应处理方法的流程图；

图2是本发明实施例提供的一种基于局部网格分裂的网格偏离空时自适应处理方法的网格分裂示意图；

图3a是利用现有技术的LSCSE方法挑选的杂波原子的分布图；

图3b是利用本发明实施例的方法挑选的杂波原子的分布图；

图4是本发明实施例的方法与现有技术方法输出的信杂噪比的对比图。

具体实施方式

为了进一步阐述本发明为达成预定发明目的所采取的技术手段及功效，以下结合附图及具体实施方式，对依据本发明提出的基于局部网格分裂的网格偏离空时自适应处理方法进行详细说明。

有关本发明的前述及其他技术内容、特点及功效，在以下配合附图的具体实施方式详细说明中即可清楚地呈现。通过具体实施方式的说明，可对本发明为达成预定目的所采取的技术手段及功效进行更加深入且具体地了解，然而所附附图仅是提供参考与说明之用，并非用来对本发明的技术方案加以限制。

应当说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的物品或者设备中还存在另外的相同要素。

请参见图1，图1是本发明实施例提供的一种基于局部网格分裂的网格偏离空时自适应处理方法的流程图。本实施例的网格偏离空时自适应处理方法包括：

步骤1：获取杂波空间的空时自适应处理字典及杂波协方差矩阵。

首先，对信号模型进行简要说明。

利用机载雷达系统观测空间中静止的杂波散射点，假设一机载系统采用线型阵列天线，阵元个数为N，阵元间距为d。以重频f_r在相参积累时间内发射了M个脉冲。载机速度为V，则杂波散射点的空域频率f_s与归一化的多普勒频率f_d可分别表示为：

f_s＝dcosψ_a/λ

f_d＝Vcosψ_v/λf_r

其中，λ表示信号波长，ψ_a表示观测点的空间锥角，ψ_v表示观测点与载机运动方向的速度锥角。因此，第i个杂波散射点的空时导向矢量s(f_d,i,f_s,i)可以写为：

其中，

表示Kronecker积，s_d表示时域导向矢量，s_s表示空域导向矢量，表达式如下：

s_d(f_d,i)＝[1,exp(j2πf_d,i),...,exp(j2π(M-1)f_d,i)]^T

s_s(f_s,i)＝[1,exp(j2πf_s,i),...,exp(j2π(N-1)f_s,i)]^T

其中，(·)^T表示矩阵的转置操作，N_s表示利用稀疏恢复算法后空时二维平面空域频率被离散为N_s＝ρ_sN个网格，N_d表示利用稀疏恢复算法后时域频率被离散为N_d＝ρ_dM个网格，ρ_d,ρ_s＞1，表示频域的离散程度，N为阵元个数，M为脉冲个数。

获得杂波空间的空时自适应处理字典为：

其中，s(f_d,i,f_s,i)表示第i个杂波散射点的空时导向矢量。

随后，利用杂波散射点的回波信号x_l估计杂波第l个距离门的协方差矩阵：

步骤2：初始化杂波空间的正交投影矩阵

和噪声投影矩阵

初始化空时导向矢量集合

杂波原子下标集合Ω、全局迭代次数k、全局迭代门限η、局部迭代次数p、局部迭代的误差门限σ_l、局部的最大迭代次数。

具体地，Ω表示搜索到的杂波原子的下标集合，初始化为空集；Φ_Ω表示由Ω下标标志的在Φ中的导向矢量的集合，在迭代开始时初始化为空集；P_c和P_n分别表示向杂波空间的正交投影矩阵与噪声投影矩阵，将P_c和P_n初始化为单位阵；初始化全局迭代的次数k＝1，初始化全局迭代的门限η，初始化局部迭代的误差门限σ_l、局部迭代次数p、局部迭代的误差门限σ_l、局部的最大迭代次数。

步骤3：搜索空时导向矢量，重构杂波子空间

在空时自适应处理字典中搜索杂波的空时导向矢量，将其构成新的杂波空间的空时导向矢量集合Φ_Ω，并计算新的杂波空间的噪声投影矩阵P_n。

进一步地，步骤3具体包括：

步骤3.1：在空时自适应处理字典中进行全局搜索

在空时自适应处理字典中进行全局搜索，通过子空间投影的方式，挑选出与当前剩余杂波最相关的空时导向矢量作为进行局部迭代的初始杂波原子的空时二维平面坐标。

首先，通过子空间投影的方式，在空时自适应处理字典中挑选出与当前剩余杂波最相关的空时导向矢量。

表示在第k次全局迭代中选择的杂波原子的下标，则在第k次全局迭代中

通过下式进行挑选：

其中，m表示空时自适应字典中杂波原子的下标，s_m表示下标为m的杂波原子所对应的空时导向矢量。

挑选上式惩罚项最大的一个，作为在当前的杂波子空间中搜索到的杂波原子。

步骤3.2：进行局部迭代

以步骤3.1中挑选的杂波原子为中心，分别在时域以1/M、在空域以1/N为边长划分一个2×2的“田”字网格，该网格的9个顶点构成3×3的原子矩阵Φ_Ωl，计算这9个杂波原子空时权矢量的响应，通过下式计算：

其中，p表示局部分裂的次数，l表示第l个距离门，

表示第p-1次局部分裂后更新的空时导向矢量集合里的空时导向矢量，

表示第p-1次局部分裂后更新的原子所组成的空时导向矢量集合，

表示第p次局部分裂后下标为m_l的杂波原子的空时权矢量的响应。

请参见图2，利用9个点中在空时权矢量的最大响应点作为顶点，向第二与第三大响应的点做边，从而闭合成为下一次分裂的网格。在次网格内重新划分2×2的网格并更新，重复上述过程直到

其中，σ_l表示误差门限，或者达到局部的最大迭代次数。

最终，最优杂波原子

通过上式在最后一次局部迭代中被挑选出来。

步骤3.3：更新杂波子空间

在经过上述局部迭代过程后，第k次全局搜索得到的最优杂波原子为

此时，对搜索到的杂波原子下标集合Ω进行扩充

对搜索到的杂波原子的空时导向矢量集合

进行更新

分别以

和

进行更新，其中，I为单位阵，

为矩阵的求伪逆操作。

步骤3.4：判断搜索终止条件

令k＝k+1，重复步骤3.1至步骤3.3，直到终止条件成立，得到最终更新的噪声投影矩阵

所述终止条件为：

E^(k)＜ηE⁽⁰⁾

其中，E⁽⁰⁾表示所有杂波原子的能量，

表示剩余杂波原子的能量，η为预设的门限。

步骤4：计算空时自适应处理加权矢量矩阵。

根据得到的更新后的噪声投影矩阵

获得网格失配的空时自适应处理加权矢量为：

其中，s₀表示目标的空时导向矢量。

本发明实施例基于局部网格分裂的网格偏离空时自适应处理方法的效果可以通过下述仿真实验加以说明：

(1)试验条件：假设平台以150m/s的速度飞行，阵列方位向沿载机方位向摆放，阵元个数N＝16，脉冲数M＝16，脉冲重频6000Hz，载机高度8km，下视角为30°。设中心距离环360个杂波单元沿方位向80°到100°均匀分布。杂噪比为35dB。字典大小N_s＝5N，N_d＝5M。

(2)试验内容与结果分析：

试验1:

检验了挑选的杂波原子在空时平面上的分布情况，设置现有技术LSCSE-STAP方法中局部的字典中N_ld＝N_ls＝5，本发明实施例的方法(附图中表示为LMSSE-STAP)中的最大分裂次数为3。全局迭代门限η＝0.01E⁽⁰⁾，迭代终止次数为25。

请参见图3a和图3b，图3a是利用现有技术的LSCSE方法挑选的杂波原子的分布图；图3b是利用本发明实施例的方法挑选的杂波原子的分布图。可以看出，本发明实施例的方法比LSCSE-STAP方法具有更高的搜索精度。

试验2:

通过比较输出的信杂噪比损失作为性能检验的结果，信杂噪比损失定义如下：

其中，

表示STAP加权矢量，R表示估计的杂波矩阵。采用试验1中的实验设置，设目标的空间频率为0.06。请参见图4，图4是本发明实施例的方法与现有技术方法输出的信杂噪比的对比图，其中，横坐标为归一化多普勒频率，纵坐标为信杂噪比损失。OPT表示最优处理器。可以看出，本方法实施例的方法LMSSE-STAP比现有技术的方法LSCSE-STAP与SRCN-STAP更加接近于OPT曲线，即利用本方法实施例的方法LMSSE-STAP搜索到的杂波原子所张成的杂波子空间更加接近于真实的杂波子空间。因此，本方法实施例的方法LMSSE-STAP比现有技术的方法LSCSE-STAP与SRCN-STAP具有更好的性能。从实验结果中，可以看出本方法实施例的方法LMSSE-STAP比现有技术的方法LSCSE-STAP与SRCN-STAP在计算量较少的情况下可以达到更高的性能。

其中，SRCN-STAP(Sparse Representationbased ClutterNulling STAP)是一种基于原子挑选的杂波零陷STAP方法，该方法的基本思路为通过迭代搜索在STAP字典中挑选杂波原子构造杂波子空间，挑选准则为在之前未被选中的剩余杂波原子中，选出空时权响应最大的那一个杂波原子，并利用本次挑选的杂波原子更新噪声投影矩阵，再利用最后更新的噪声投影矩阵计算STAP权矢量。

通过实验仿真，可以证明本发明实施例的方法可以在计算较少时获得较高的精度。综合图3a和图3b的挑选原子在空时谱平面上的分布，图4的信杂噪比损失曲线，可以得出：通过局部的网格分裂，本发明实施例的方法可以快速准确地找到失配的空时网格，有效地提升了稀疏空时处理的性能，为利用阴影辅助检测动目标提供了有效的解决途径。

综上，本发明基于局部网格分裂的网格偏离空时自适应处理方法的克服了估计精度与局部网格划分的个数有关的问题，减小了运算量，节约了运算资源，并具有更优的处理效率。该方法提高了杂波原子的搜索精度，输出信杂噪比曲线凹口变窄，有效地提升了稀疏空时处理的性能。

以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演或替换，都应当视为属于本发明的保护范围。

Claims (7)

1.一种基于局部网格分裂的网格偏离空时自适应处理方法，其特征在于，包括：

步骤1：获取杂波空间的空时自适应处理字典及杂波协方差矩阵；

步骤2：初始化杂波空间的正交投影矩阵

和噪声投影矩阵

初始化空时导向矢量集合

杂波原子下标集合Ω、全局迭代次数k、全局迭代门限η、局部迭代次数p、局部迭代的误差门限σ_l、局部的最大迭代次数；

步骤3：在所述空时自适应处理字典中搜索杂波的空时导向矢量，构成更新后的杂波空间的空时导向矢量集合Φ_Ω并计算更新后的杂波空间的噪声投影矩阵P_n；

步骤4：根据更新后的杂波空间的噪声投影矩阵P_n获得网格失配的空时自适应处理加权矢量。

2.根据权利要求1所述的基于局部网格分裂的网格偏离空时自适应处理方法，其特征在于，所述步骤1包括：

利用机载雷达系统观测空间中的杂波散射点，计算第i个杂波散射点的空时导向矢量s(f_d,i,f_s,i)，并根据所有杂波散射点的空时导向矢量s(f_d,i,f_s,i)构造杂波空间的空时自适应处理字典Φ：

其中，第i个杂波散射点的空时导向矢量s(f_d,i,f_s,i)为：

其中，

表示Kronecker积，s_d表示杂波散射点的时域导向矢量，s_s表示杂波散射点的空域导向矢量：

s_d(f_d,i)＝[1,exp(j2πf_d,i),...,exp(j2π(M-1)f_d,i)]^T

s_s(f_s,i)＝[1,exp(j2πf_s,i),...,exp(j2π(N-1)f_s,i)]^T

其中，(·)^T表示矩阵的转置操作，N_s表示利用稀疏恢复算法后空时二维平面空域频率被离散为N_s＝ρ_sN个网格，N_d表示利用稀疏恢复算法后时域频率被离散为N_d＝ρ_dM个网格，ρ_d,ρ_s＞1，表示频域的离散程度，N为阵元个数，M为脉冲个数，f_s,i表示第i个杂波散射点的空域频率，f_d,i表示第i杂波散射点的归一化的多普勒频率；

根据杂波散射点的回波数据x_l构造杂波的协方差矩阵

其中，l表示第l个距离门。

3.根据权利要求2所述的基于局部网格分裂的网格偏离空时自适应处理方法，其特征在于，所述步骤3包括：

步骤3.1：在所述空时自适应处理字典中进行全局搜索，通过子空间投影的方式挑选出与当前剩余杂波最相关的空时导向矢量作为进行局部迭代的初始杂波原子的空时二维平面坐标；

步骤3.2：以步骤3.1中挑选的原始杂波原子为中心进行多次网格分裂，直至杂波原子空时权矢量的响应满足局部迭代误差门限，对应的杂波原子为本次全局搜索到的杂波原子，其下标记为

步骤3.3：对搜索到的杂波原子下标集合Ω进行扩充

对搜索到的杂波原子的空时导向矢量集合

进行更新

对杂波空间的正交投影矩阵

进行更新

对杂波空间的噪声投影矩阵

进行更新

其中I为单位阵，

为矩阵的求伪逆操作，k表示全局迭代的次数；

步骤3.4：令k＝k+1，重复步骤3.1至步骤3.3，直到剩余杂波能量达到全局搜索的门限，获得此时更新的噪声投影矩阵

4.根据权利要求3所述的基于局部网格分裂的网格偏离空时自适应处理方法，其特征在于，所述步骤3.1包括：

步骤3.11：通过子空间投影的方式，在所述空时自适应处理字典中挑选出与当前剩余杂波最相关的空时导向矢量，

表示第k次全局搜索到的原子下标，

的表达式为：

其中，m表示空时自适应字典中杂波原子的下标，s_m表示下标为m的杂波原子所对应的空时导向矢量；

步骤3.12：挑选

的表达式中惩罚项最大的一个，作为在当前的杂波子空间中搜索到的杂波原子。

5.根据权利要求4所述的基于局部网格分裂的网格偏离空时自适应处理方法，其特征在于，所述步骤3.2包括：

步骤3.21：以步骤3.1中挑选的原始杂波原子为中心，分别在时域以1/M、在空域以1/N为边长划分一个2×2的“田”字网格，所述网格的9个顶点构成3×3的原子矩阵Φ_Ωl，计算9个杂波原子的空时权矢量的响应：

其中，p表示局部分裂的次数，l表示第l个距离门，

表示第p-1次局部分裂后更新的空时导向矢量集合里的空时导向矢量，

表示第p-1次局部分裂后更新的原子所组成的空时导向矢量集合，

表示第p次局部分裂后下标为m_l的杂波原子的空时权矢量的响应；

步骤3.22：利用所述9个杂波原子中在空时权矢量的最大响应点作为顶点，向第二与第三大响应的点做边，闭合成为下一次分裂的次网格；在次网格内重新划分2×2的网格并更新，重复步骤3.21的过程直到

其中，σ_l表示误差门限，或者达到局部的最大迭代次数；

步骤3.23：通过最后一次局部迭代挑选出最优杂波原子

6.根据权利要求5所述的基于局部网格分裂的网格偏离空时自适应处理方法，其特征在于，步骤3.4包括：

令k＝k+1，重复步骤3.1至步骤3.3，直到终止条件成立，得到最终更新的噪声投影矩阵

所述终止条件为：

其中，E⁽⁰⁾表示所有杂波原子的能量，E^(k)表示剩余杂波原子的能量，η为预设的门限。

7.根据权利要求1所述的基于局部网格分裂的网格偏离空时自适应处理方法，其特征在于，所述步骤 4包括：

根据得到的更新后的噪声投影矩阵

获得网格失配的空时自适应处理加权矢量为：

其中，s₀表示目标的空时导向矢量。

Images