CN112051552B - 一种基于多站的主瓣抗干扰方法及装置 - Google Patents

Abstract

本发明提出一种基于多站的主瓣抗干扰方法及装置，属于雷达抗干扰技术领域。该装置包括回波数据存储模块、干扰源位置初始化模块、迭代算法模块、筛选模块和干扰重构模块。该方法首先对场景建模，当所有雷达站点采样到回波信号后，对各站点采集到的回波信号通过多层迭代，输出包含冗余信息的干扰源位置集合估计值和对应临时复变量集合估计值并进行筛选，得到不包含冗余信息的干扰源位置集合估计值和对应临时复变量集合的估计值；然后对干扰信号进行重构，在回波信号中减去重构后的干扰信即可实现多站主瓣抗干扰。本发明可对从主瓣进入的随机噪声干扰进行有效消除，解决了现有技术依赖于多站干扰信号相参的问题，实现更为有效的干扰抑制。

Description

一种基于多站的主瓣抗干扰方法及装置

技术领域

本发明属于雷达抗干扰技术领域，特别提出一种基于多站的主瓣抗干扰方法及装置。

背景技术

随着电子科学技术的快速发展，电子对抗在现代战争中的作用日益增大，对雷达造成显著威胁。对于雷达系统而言，雷达抗干扰事关雷达在对抗条件下的生存能力，具有重要的地位。

干扰的形式多种多样，按照干扰源是否主动辐射电磁波可以分为有源干扰和无源干扰。而有源干扰又包括欺骗式干扰和压制式干扰。其中，有源压制式干扰通过向雷达接收天线注入强干扰信号淹没目标回波信号，从而使得雷达检测、跟踪性能下降，对雷达探测目标构成显著威胁。一般而言，压制性干扰采用随机性信号作为干扰波形(如高斯白噪声)，而如何对抗强高斯白噪声干扰成为雷达抗干扰的一个重要研究问题。

对于有源干扰而言，根据干扰能量从雷达天线旁瓣还是主瓣进入雷达接收机，又可分为旁瓣干扰和主瓣干扰两类。对于旁瓣干扰，目前的抗干扰研究工作较多，如设计超低副瓣天线、旁瓣对消、旁瓣匿影、波束成形等方法，可以实现对从旁瓣进入的干扰进行有效地抑制。但对于主瓣干扰，传统的旁瓣抗干扰算法不能够直接应用。例如，如果使用波束自适应成形算法做主瓣抗干扰，由于干扰和目标回波都从主瓣进入雷达接收天线，在对干扰方向形成零陷的同时雷达主瓣也会产生严重的变形，从而影响雷达检测、跟踪性能。目前，主瓣抗干扰的方法主要有和差波束主瓣对消方法和盲源分离方法，前者需要天线主波束对准目标，而这在干扰环境下通常难以实现；后者在某些情况可以分离出目标信号，但分离后通道之间关系是未知的，给目标测角带来很大的困难。

以上提到的方法通常是基于一个雷达站点进行干扰抑制。近些年来，随着有关多站联合信号处理的研究日益深入，使用多雷达站点实现主瓣抗干扰的方法逐渐引起研究人员的关注。由于多站雷达相比单站雷达有更高的角度分辨力和更好的空间多样性，多站主瓣抗干扰很有可能获得比单站方法更好的抗干扰效果。在目前已有的多站主瓣抗干扰的研究中，一种经典的思路是通过各站回波信号构建协方差矩阵，应用矩阵奇异值分解估计干扰所在的子空间，通过子空间投影的方法去除干扰，但这样的方法要求各站接收到的信号是相参的，即在对于不同采样时刻，各站回波信号的相位差固定。而实际情况中，当干扰机使用随机旁瓣调制时这样的条件并不一定满足，从而导致该抗干扰算法失效。

发明内容

本发明的目的是为克服已有技术的不足之处，提出一种基于多站的主瓣抗干扰方法及装置。本发明不依赖于多站雷达之间所接受到的干扰信号相参这一假设，能够实现对从主瓣进入的随机噪声干扰进行有效消除，解决了现有技术依赖于多站干扰信号相参的问题，实现更为有效的干扰抑制。

本发明提出一种基于多站的主瓣抗干扰方法，其特征在于，包括以下步骤：

1)场景建模；

设在一个二维的平面上，共包括K个雷达站点，K≥2，第k个雷达站点的平面坐标表示为u_k＝[u_k，x，u_k，y]^T，k＝1，...，K，其中u_k，x表示第k个雷达站点的横坐标，u_k，y表示第k个雷达站点的纵坐标；每个雷达站点装备有一个均匀线性的接收天线阵列，其中第k个雷达站点的接收天线阵列有N_k个阵元，N_k≥2；

在雷达探测的目标附近，存在Q个干扰源，Q≥1，第q个干扰源的平面坐标表示为v_q＝[v_q，x，v_q，y]^T，q＝1，...，Q，其中v_q，x表示第q个干扰源的横坐标，v_q，y表示第q个干扰源的纵坐标；每个干扰源向空间持续辐射噪声干扰信号，该干扰信号和目标回波信号一同作为回波信号从各雷达站点接收天线阵列的主瓣进入；

2)数据采集；

将第k个雷达站点第n_k个阵元在T_k时刻进行采样得到回波信号记为 则T_k时刻第k个雷达站点采样得到的回波信号构成向量s_k：

3)抗干扰信号处理；

当所有雷达站点采样到一次回波信号后，对各站点采集到的回波信号进行抗干扰信号处理，具体步骤如下：

3-1)算法初始化；

3-1-1)将Q个干扰源位置集合初始化为空集,V⁰＝{}；

3-1-2)用Φ表示每个雷达站点对于每个干扰源的临时复变量的集合，其中第k个雷达站点对于第q个干扰源的临时复变量记为φ_q，k，则Φ＝{φ_q，k，k＝1，...，K，q＝1，2，...}；

3-1-3)假设已知每个干扰源位置在空间中的大致范围，定义为Ω；

3-1-4)定义第k个雷达站点的天线阵列的第n_k个阵元对于第q个干扰源的导引量为其中，d_k表示第k个阵列的阵元间距，λ表示回波信号的波长，θ_k(v_q)表示第q个干扰源对第k个雷达站点的天线阵列的方向角，对应正弦值表示为:

则第k个雷达站点天线阵列的各阵元对于第q个干扰源的导引量组成的向量表示为：

3-1-5)令第k个雷达站点对应的残差变量e_k＝s_k，k＝1，...，K；

3-2)执行迭代1：

3-2-1)用l表示迭代1的迭代次数，令初始迭代次数l＝1，设置迭代1的迭代次数上限为l_max，l_max＞Q；

3-2-2)将Ω均匀划分为I个格点，每个格点对应空间中的记为坐标c_i，i＝1，…,I，用Ω_I表示所有格点对应坐标的集合，Ω_I＝{c₁，...，c_I}；

3-2-3)令第l次迭代的参数v_l为：

利用v_l更新得到第l次迭代后的干扰源位置集合，使得V^l＝{V^l-1，v_l}，V^l为第l次迭代后的干扰源位置集合；

3-2-4)根据步骤3-2-3)得到的集合V^l，令

其中，B_k，l＝[g_k(v₁)，...，g_k(v_l)]，B_k，l为g_k(v₁)，...，g_k(v_l)构成的临时矩阵变量，{Φ^l}_k表示集合Φ^l的第k个向量元素，Φ^l表示第l次迭代的临时复变量集合，其中集合Φ^l的第k个向量元素{Φ^l}_k中的第q个元素表示为

3-2-5)通过迭代2对V^l进行精确估计；

用p表示迭代2的迭代次数，令初始迭代次数p＝1，设置迭代2的迭代次数上限为p_max；

将步骤3-2-3)和3-2-4)得到的V^l和Φ^l分别作为迭代2的初值和/>即/> 则：

其中，为迭代2第p次迭代的v_q值，ζ^p为迭代2第p次迭代的步长，/>表示对复数取实部的运算符，q＝1，...，l；

迭代2第p-1次迭代参数∈^p-1定义为：

其中，⊙表示Hadamard积，n_k＝[1，...，N_k]^T；

迭代2第p-1次迭代参数α_p-1，β_p-1分别定义为：

当到达迭代2的次数上限p_max时，通过迭代2得到

3-2-6)通过迭代3对Φ进行精确估计；

用z表示迭代2的迭代次数，令初始迭代次数z＝1，设置迭代3的迭代次数上限为z_max；

其中，χ^z为迭代3第z次迭代的步长，η为权重因子，q＝1，...，l，k＝1，...，K；

迭代3第z-1次迭代参数γ_z-1定义为：

当到达迭代3的次数上限z_max时，通过迭代3得到

3-2-7)更新

3-2-8)令l＝l+1，然后重新返回步骤3-2-3)，直至到达迭代1的迭代次数上限l_max时，通过迭代1得到和/>其中，/>为最终的干扰源位置集合，/>为最终的每个雷达站点对于每个干扰源的临时复变量的集合；

3-3)筛选干扰源位置；

从步骤3-2)得到的最终的干扰源位置集合中移除对应/>且满足下式的元素v_q，q＝1，...l_max：

其中，T_h表示预设的一个界；

移除完毕后，集合剩余元素数目即为干扰源数目Q的估计值/>集合/>剩余元素组成对应干扰源位置V的估计值/>对应的临时复变量集合的估计值为/>

3-4)干扰信号重构和去除；

根据步骤3-3)得到的估计值和/>对每个雷达站点T_k时刻的干扰信号进行重构：

其中，h_k表示对第k个雷达站点接收到的干扰信号的重构，k＝1，...，K；

将重构后的干扰信号h_k从回波信号s_k中减去，实现对第k个雷达站点T_k时刻采样信号的主瓣干扰的抑制。

本发明还提出一种基于上述方法的基于多站的主瓣抗干扰装置，其特征在于，包括回波数据存储模块、干扰源位置初始化模块、迭代算法模块、筛选模块、干扰重构模块；

所述回波数据存储模块的输入端连接各雷达站点，回波数据存储模块的输出端分别连接干扰源位置初始化模块的输入端、迭代算法模块的输入端和干扰重构模块的输入端，干扰源位置初始化模块的输出端连接迭代算法模块的输入端，迭代算法模块的输出端连接筛选模块的输入端，筛选模块的输出端连接干扰重构模块的输入端；

所述回波数据存储模块用于获取各雷达站点的回波信号，并将接收到的各站点的回波信号数据进行存储，然后将存储的数据分别发送给干扰源位置初始化模块、迭代算法模块和干扰重构模块；

所述干扰源位置初始化模块利用从回波数据存储模块接收到的各站点的回波信号数据，计算干扰源位置集合和对应临时复变量集合的初值估计，然后发送给迭代算法模块；

所述迭代算法模块用于将从回波数据存储模块接收到的各站点的回波信号数据、从干扰源位置初始化模块接收到的干扰源位置集合和对应临时复变量集合的初值估计代入到迭代算法模型中，输出包含冗余信息的干扰源位置集合的估计值和对应临时复变量集合的估计值，然后将估计值发送给筛选模块；

所述筛选模块用于对从迭代算法模块接收到的包含冗余信息的干扰源位置集合的估计值和对应临时复变量集合的估计值进行筛选，得到不包含冗余信息的干扰源位置集合估计值和对应临时复变量集合的估计值并发送给干扰重构模块；

所述干扰重构模块利用从筛选模块接收到的不包含冗余信息的干扰源位置集合估计值和对应临时复变量集合的估计值进行干扰信号重构，得到干扰信号的精确形式，然后再从回波数据存储模块接收到的各站点的雷达回波信号数据中减去对应的干扰信号，输出去除干扰后的信号，实现多站主瓣抗干扰。

本发明的特点及有益效果在于：

本发明使用各站阵列测角能力对干扰信号进行重构，其带来的效益体现在不需要各站接收到的信号相参的假设，扩大了主瓣抗干扰的应用场景。

本发明提出一种基于多站的主瓣抗干扰方法及装置，基于各站雷达回波数据，获取干扰源位置集合和对应复系数集合的初值，并根据迭代计算和筛选，输出干扰源位置集合和对应复系数集合的精确估计，无需各站接收到的干扰信号相参的假设，通过干扰源位置集合和对应复系数集合的精确估计即可重构干扰信号，从而将重构的干扰信号从回波信号中减去，即可实现更为有效的干扰抑制。本发明不依赖于上述多站雷达之间所接受到的干扰信号相参这一假设，能够实现对从主瓣进入的随机噪声干扰进行有效消除，解决了现有技术依赖于多站干扰信号相参的问题，实现更为有效的干扰抑制。

附图说明

图1是本发明方法的整体流程图；

图2是本发明装置的结构示意图。

具体实施方式

本发明提出一种基于多站的主瓣抗干扰方法及装置，下面结合附图和具体实施例对本发明进一步详细说明如下。

本发明提出一种基于多站的主瓣抗干扰方法，整体流程如图1所示，包括以下步骤：

1)场景建模；

本发明考虑在一个二维的平面上，共包括K个雷达站点，K≥2，第k个雷达站点的平面坐标表示为其中u_k，x表示第k个雷达站点的横坐标，u_k，y表示第k个雷达站点的纵坐标。每个雷达站点装备有一个均匀线性的接收天线阵列，可以实现对来波方向角度的测量，其中第k个雷达站点的接收天线阵列有N_k个阵元，N_k≥2(每个雷达站点接收天线阵列的阵元数目可不相同)。

在雷达探测的目标附近，存在Q个干扰源，Q≥1(Q为未知)，第q个干扰源的平面坐标表示为v_q＝[v_q，x，v_q，y]^T，q＝1，...，Q，其中v_q，x表示第q个干扰源的横坐标，v_q，y表示第q个干扰源的纵坐标。每个干扰源向空间持续辐射噪声干扰信号，该干扰信号和目标回波信号一同作为回波信号从各雷达站点接收天线阵列的主瓣进入。

2)数据采集；

将t时刻第k个雷达站点的天线阵列中第n个阵元接收到的回波信号表示为s_n，k(t)。对于各雷达站点接收的回波信号，假设第k个雷达站点是在T_k时刻进行采样得到该时刻的单快拍数据，将第k个雷达站点第n_k个阵元在T_k时刻进行采样得到回波信号记为则T_k时刻第k个雷达站点采样得到的回波信号构成向量s_k：

3)抗干扰信号处理；

当所有雷达站点采样到一次回波信号后，对各站点采集到的回波信号进行抗干扰信号处理，具体步骤如下：

3-1)算法初始化：

3-1-1)将Q个干扰源位置集合初始化为空集,V⁰＝{}；

3-1-2)用Φ表示每个雷达站点对于每个干扰源的临时复变量的集合，其中第k个雷达站点对于第q个干扰源的临时复变量记为φ_q，k，则Φ＝{φ_q，k，k＝1，...，K，q＝1，2，...}。

3-1-3)假设已知每个干扰源位置在空间中的大致范围，定义为Ω。

3-1-4)定义第k个雷达站点的天线阵列的第n_k个阵元对于第q个干扰源的导引量为其中，d_k表示第k个阵列的阵元间距，λ表示回波信号的波长，θ_k(v_q)表示第q个干扰源对第k个雷达站点的天线阵列的方向角，对应正弦值表示为：

则第k个雷达站点天线阵列的各阵元对于第q个干扰源的导引量组成的向量表示为：

3-1-5)令第k个雷达站点对应的残差变量e_k＝s_k，k＝1，...，K；

3-2)执行迭代1(用l表示迭代次数，从1到l_max，l_max＞Q)：

3-2-1)令初始迭代次数l＝1；

3-2-2)将Ω均匀划分为I个格点(通常可以以雷达距离分辨力作为间隔来划分Ω)，每个格点对应空间中的记为坐标c_i，i＝1，…,I，用Ω_I表示所有格点对应坐标的集合，即划分I个格点后所有格点的坐标集合为Ω_I＝{c₁，...，c_I}；

3-2-3)令第l次迭代的参数v_l为：

利用v_l更新得到第l次迭代后的干扰源位置集合，使得V^l＝{V^l-1，v_l}，V^l为第l次迭代后的干扰源位置集合。

3-2-4)根据步骤3-2-3)得到的集合V^l，令

其中，B_k，l＝[g_k(v₁)，...，g_k(v_l)]，B_k，l为g_k(v₁)，...，g_k(v_l)构成的临时矩阵变量，{Φ^l}_k表示集合Φ^l的第k个向量元素，Φ^l表示第l次迭代的临时复变量集合，其中集合Φ^l的第k个向量元素{Φ^l}_k中的第q个元素表示为

3-2-5)通过迭代2对V^l进行精确估计：执行迭代2(用p表示迭代次数，从1到p_max，p_max需设置得足够大使得迭代2能够收敛。由于迭代2为迭代1中的子迭代步骤，为清晰起见，迭代2中的参数形式上省略表示迭代1迭代次数的符号l)，其中将步骤3-2-3)和3-2-4)得到的V^l和Φ^l分别作为迭代2的初值的初值和/>即/> 则：

其中，为迭代2第p次迭代的v_q值，/>为迭代2第p次迭代的步长，可以由经典梯度下降算法中的回溯直线搜索算法确定，/>表示对复数取实部的运算符，q＝1，...，l；

迭代2第p-1次迭代参数∈^p-1定义为：

其中，⊙表示Hadamard积，n_k＝[1，...，N_k]^T；

迭代2第p-1次迭代参数α_p-1，β_p-1分别定义为：

当到达迭代2的次数上限p_max时，通过迭代2得到

3-2-6)通过迭代3对Φ进行精确估计：执行迭代3(用z表示迭代次数，从1到z_max，z_max需设置得足够大使得迭代3能够收敛)：

其中，χ^z为迭代3第z次迭代的步长，可以由经典梯度下降算法中的回溯直线搜索算法确定，η为根据经验设置的权重因子，q＝1，...，l，k＝1，...，K，迭代3第z-1次迭代参数γ_z-1定义为：

当到达迭代3的次数上限z_max时，通过迭代3得到

3-2-7)更新

3-2-8)令l＝l+1，然后重新返回步骤3-2-3)，直至到达迭代1的迭代次数上限l_max时，通过迭代1得到和/>其中，/>为最终的干扰源位置集合，/>为最终的每个雷达站点对于每个干扰源的临时复变量的集合。

3-3)筛选干扰源位置；

从步骤3-2)得到的最终的干扰源位置集合中移除对应/>且满足下式的元素v_q，q＝1，...l_max：

其中，T_h表示预设的一个界，其设置原则是用以将集合中信号能量不足的元素滤除，由于干扰信号的能量通常远高于目标信号能量，因此可以从C_q的计算结果中区分干扰和目标的T_h的设置范围是广泛的。通过筛选后，集合/>剩余元素数目即为干扰源数目Q的估计值/>集合/>剩余元素组成对应干扰源位置V的估计值/>对应的临时复变量集合的估计值为/>

3-4)干扰信号重构和去除；

根据步骤3-3)估计值和/>对每个雷达站点T_k时刻的干扰信号进行重构：

其中，h_k表示对第k个雷达站点接收到的干扰信号的重构，k＝1，...，K。将重构后的干扰信号h_k从回波信号s_k中减去，即可以实现对第k个雷达站点T_k时刻采样信号的主瓣干扰的抑制。

总结本发明方法有如下特点：

1)该方法不需要各站干扰信号相参的假设；

2)该方法可以得到干扰源的坐标；

3)该方法可以实现主瓣干扰消除。

本发明还提出一种基于上述方法的基于多站的主瓣抗干扰装置，结构如图2所示，包括回波数据存储模块、干扰源位置初始化模块、迭代算法模块、筛选模块、干扰重构模块。

所述回波数据存储模块的输入端通过与各雷达站点通信传输获得各站回波信号并进行存储，回波数据存储模块的输出端通过物理链路分别连接到干扰源位置初始化模块的输入端、迭代算法模块的输入端和干扰重构模块的输入端，干扰源位置初始化模块的输出端通过物理链路连接到迭代算法模块的输入端，迭代算法模块的输出端通过物理链路连接到筛选模块的输入端，筛选模块的输出端通过物理链路连接到干扰重构模块的输入端，最后干扰重构模块的输出端输出去除干扰后的信号。

所述回波数据存储模块用于获取各雷达站点的回波信号，并将接收到的各站点的回波信号数据进行存储，然后将存储的数据分别发送给干扰源位置初始化模块、迭代算法模块和干扰重构模块。

所述干扰源位置初始化模块使用从回波数据存储模块接收到的各站点的回波信号数据，计算干扰源位置集合和对应临时复变量集合的初值估计，然后发送给迭代算法模块。

所述迭代算法模块用于将从回波数据存储模块接收到的各站点的回波信号数据、从干扰源位置初始化模块接收到的干扰源位置集合和对应临时复变量集合的初值估计代入到迭代算法模型中，输出包含冗余信息的干扰源位置集合的估计值和对应临时复变量集合的估计值，然后将估计值发送给筛选模块。

所述筛选模块用于对从迭代算法模块接收到的包含冗余信息的干扰源位置集合的估计值和对应临时复变量集合的估计值进行筛选，得到不包含冗余信息的干扰源位置集合估计值和对应临时复变量集合的估计值并发送给干扰重构模块。

所述干扰重构模块利用从筛选模块接收到的不包含冗余信息的干扰源位置集合估计值和对应临时复变量集合的估计值进行干扰信号重构，得到干扰信号的精确形式，然后在从回波数据存储模块接收到的各站点的雷达回波信号数据中减去对应的干扰信号，输出去除干扰后的信号，实现多站主瓣抗干扰。

Claims (2)

1.一种基于多站的主瓣抗干扰方法，其特征在于，包括以下步骤：

1)场景建模；

设在一个二维的平面上，共包括K个雷达站点，K≥2，第k个雷达站点的平面坐标表示为u_k＝[u_k，x，u_k，y]^T，k＝1，...，K，其中u_k，x表示第k个雷达站点的横坐标，u_k，y表示第k个雷达站点的纵坐标；每个雷达站点装备有一个均匀线性的接收天线阵列，其中第k个雷达站点的接收天线阵列有N_k个阵元，N_k≥2；

在雷达探测的目标附近，存在Q个干扰源，Q≥1，第q个干扰源的平面坐标表示为v_q＝[v_q，x，v_q，y]^T，q＝1，…，Q，其中v_q，x表示第q个干扰源的横坐标，v_q，y表示第q个干扰源的纵坐标；每个干扰源向空间持续辐射噪声干扰信号，该干扰信号和目标回波信号一同作为回波信号从各雷达站点接收天线阵列的主瓣进入；

2)数据采集；

将第k个雷达站点第n_k个阵元在T_k时刻进行采样得到回波信号记为 则T_k时刻第k个雷达站点采样得到的回波信号构成向量s_k：

3)抗干扰信号处理；

当所有雷达站点采样到一次回波信号后，对各站点采集到的回波信号进行抗干扰信号处理，具体步骤如下：

3-1)算法初始化；

3-1-1)将Q个干扰源位置集合初始化为空集，V⁰＝{}；

3-1-2)用Φ表示每个雷达站点对于每个干扰源的临时复变量的集合，其中第k个雷达站点对于第q个干扰源的临时复变量记为φ_q，k，则Φ＝{φ_q，k，k＝1，...，K，q＝1，2，...}；

3-1-3)假设已知每个干扰源位置在空间中的大致范围，定义为Ω；

3-1-4)定义第k个雷达站点的天线阵列的第n_k个阵元对于第q个干扰源的导引量为其中，d_k表示第k个阵列的阵元间距，λ表示回波信号的波长，θ_k(v_q)表示第q个干扰源对第k个雷达站点的天线阵列的方向角，对应正弦值表示为：

则第k个雷达站点天线阵列的各阵元对于第q个干扰源的导引量组成的向量表示为：

3-1-5)令第k个雷达站点对应的残差变量e_k＝s_k,k＝1，…，K；

3-2)执行迭代1：

3-2-1)用l表示迭代1的迭代次数，令初始迭代次数l＝1，设置迭代1的迭代次数上限为l_max，l_max＞Q；

3-2-2)将Ω均匀划分为I个格点，每个格点对应空间中的记为坐标c_i，i＝1，…，I，用Ω_I表示所有格点对应坐标的集合，Ω_I＝{c₁，...，c_I}；

3-2-3)令第l次迭代的参数v_l为：

利用v_l更新得到第l次迭代后的干扰源位置集合，使得V^l＝{V^l-1,v_l}，V^l为第l次迭代后的干扰源位置集合；

3-2-4)根据步骤3-2-3)得到的集合V^l，令

其中，B_k,l＝[g_k(v₁)，...，g_k(v_l)]，B_k,l为g_k(v₁)，...，g_k(v_l)构成的临时矩阵变量，{Φ^l}_k表示集合Φ^l的第k个向量元素，Φ^l表示第l次迭代的临时复变量集合，其中集合Φ^l的第k个向量元素{Φ^l}_k中的第q个元素表示为

3-2-5)通过迭代2对V^l进行精确估计；

用p表示迭代2的迭代次数，令初始迭代次数p＝1，设置迭代2的迭代次数上限为p_max；

将步骤3-2-3)和3-2-4)得到的V^l和Φ^l分别作为迭代2的初值和/>即/>v_q,v_q∈V^l，/>则：

其中，为迭代2第p次迭代的v_q值，ζ^p为迭代2第p次迭代的步长，/>表示对复数取实部的运算符，q＝1，…，l；

迭代2第p-1次迭代参数∈^p-1定义为：

其中，⊙表示Hadamard积，n_k＝[1，...，N_k]^T；

迭代2第p-1次迭代参数α_p-1，β_p-1分别定义为：

当到达迭代2的次数上限p_max时，通过迭代2得到

3-2-6)通过迭代3对Φ进行精确估计；

用z表示迭代2的迭代次数，令初始迭代次数z＝1，设置迭代3的迭代次数上限为z_max；

其中，χ^z为迭代3第z次迭代的步长，η为权重因子，q＝1，...，l，k＝1，...，K；

迭代3第z-1次迭代参数γ_z-1定义为：

当到达迭代3的次数上限z_max时，通过迭代3得到

3-2-7)更新

3-2-8)令l＝l+1，然后重新返回步骤3-2-3)，直至到达迭代1的迭代次数上限l_max时，通过迭代1得到和/>其中，/>为最终的干扰源位置集合，/>为最终的每个雷达站点对于每个干扰源的临时复变量的集合；

3-3)筛选干扰源位置；

从步骤3-2)得到的最终的干扰源位置集合中移除对应/>且满足下式的元素v_q,q＝1，...l_max：

其中，T_h表示预设的一个界；

移除完毕后，集合剩余元素数目即为干扰源数目Q的估计值/>集合/>剩余元素组成对应干扰源位置V的估计值/> 对应的临时复变量集合的估计值为/>

3-4)干扰信号重构和去除；

根据步骤3-3)得到的估计值和/>对每个雷达站点T_k时刻的干扰信号进行重构：

其中，h_k表示对第k个雷达站点接收到的干扰信号的重构，k＝1，...，K；

将重构后的干扰信号h_k从回波信号s_k中减去，实现对第k个雷达站点T_k时刻采样信号的主瓣干扰的抑制。

2.一种基于如权利要求1所述方法的基于多站的主瓣抗干扰装置，其特征在于，包括回波数据存储模块、干扰源位置初始化模块、迭代算法模块、筛选模块、干扰重构模块；

所述回波数据存储模块的输入端连接各雷达站点，回波数据存储模块的输出端分别连接干扰源位置初始化模块的输入端、迭代算法模块的输入端和干扰重构模块的输入端，干扰源位置初始化模块的输出端连接迭代算法模块的输入端，迭代算法模块的输出端连接筛选模块的输入端，筛选模块的输出端连接干扰重构模块的输入端；

所述回波数据存储模块用于获取各雷达站点的回波信号，并将接收到的各站点的回波信号数据进行存储，然后将存储的数据分别发送给干扰源位置初始化模块、迭代算法模块和干扰重构模块；

所述干扰源位置初始化模块利用从回波数据存储模块接收到的各站点的回波信号数据，计算干扰源位置集合和对应临时复变量集合的初值估计，然后发送给迭代算法模块；

所述迭代算法模块用于将从回波数据存储模块接收到的各站点的回波信号数据、从干扰源位置初始化模块接收到的干扰源位置集合和对应临时复变量集合的初值估计代入到迭代算法模型中，输出包含冗余信息的干扰源位置集合的估计值和对应临时复变量集合的估计值，然后将估计值发送给筛选模块；

所述筛选模块用于对从迭代算法模块接收到的包含冗余信息的干扰源位置集合的估计值和对应临时复变量集合的估计值进行筛选，得到不包含冗余信息的干扰源位置集合估计值和对应临时复变量集合的估计值并发送给干扰重构模块；

所述干扰重构模块利用从筛选模块接收到的不包含冗余信息的干扰源位置集合估计值和对应临时复变量集合的估计值进行干扰信号重构，得到干扰信号的精确形式，然后再从回波数据存储模块接收到的各站点的雷达回波信号数据中减去对应的干扰信号，输出去除干扰后的信号，实现多站主瓣抗干扰。