一种雷达抗主瓣干扰方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN113687314A(43)申请公布日2021.11.23(21)申请号202110984150.0(22)申请日2021.08.25(71)申请人西安电子科技大学地址710000陕西省西安市雁塔区太白南路2号(72)发明人宋万杰张璐沈佳波陈伯孝张龙涛(74)专利代理机构西安嘉思特知识产权代理事务所(普通合伙)61230代理人刘长春(51)Int.Cl.G01S7/36(2006.01)G06F17/16(2006.01)权利要求书2页说明书10页附图7页(54)发明名称一种雷达抗主瓣干扰方法(57)摘要本发明公开了一种雷达抗主瓣干扰方法，方法包括：接收雷达多通道的混合信号；针对每个通道的混合信号，对所述混合信号进行去中心化处理得到零均值预处理信号；对所述零均值预处理信号进行等变自适应盲源分离处理得到分离信号；对所述分离信号进行脉冲压缩处理；根据峰值检测方法检测每个通道的脉冲压缩后的分离信号确定目标信号，得到雷达抗主瓣干扰结果。本发明提高了抗主瓣干扰效果。CN113687314ACN113687314A权利要求书1/2页1.一种雷达抗主瓣干扰方法，其特征在于，包括：接收雷达多通道的混合信号；针对每个通道的混合信号，对所述混合信号进行去中心化处理得到零均值预处理信号；对所述零均值预处理信号进行等变自适应盲源分离处理得到分离信号；对所述分离信号进行脉冲压缩处理；根据峰值检测方法检测每个通道的脉冲压缩后的分离信号确定目标信号，得到雷达抗主瓣干扰结果。2.根据权利要求1所述的雷达抗主瓣干扰方法，其特征在于，接收的所述混合信号公式表示为：x(t)＝Hs(t)+n(t)；T其中，x(t)表示在t时刻N个通道的混合信号，s(t)＝[s1(t),s2(t),…,sM(t)]表示在t时刻的待处理目标信号和干扰信号，M表示待处理目标信号和干扰信号总数量，n(t)＝[n1T(t),n2(t),…,nM(t)]表示在t时刻的噪声，H表示N×M的混合满秩矩阵，表示为hij表示第i个通道中第j个信号源的响应因子，所述信号源为待处理目标信号的信号源或干扰信号的信号源。3.根据权利要求2所述的雷达抗主瓣干扰方法，其特征在于，对所述混合信号进行去中心化处理得到零均值预处理信号公式表示为：其中，表示在t时刻的所述零均值预处理信号，x(t)表示在t时刻的所述混合信号，E[·]求期望。4.根据权利要求3所述的雷达抗主瓣干扰方法，其特征在于，所述对所述零均值预处理信号进行等变自适应盲源分离处理得到分离信号，包括：初始化分离矩阵和步长因子，及选取非线性激活函数；其中，所述非线性激活函数依照混合信号源的概率分布特性选取的，所述混合信号源包括所述待处理目标信号的信号源和所述干扰信号的信号源；根据所述分离矩阵和所述零均值预处理信号计算观测信号；根据所述观测信号、所述步长因子和所述非线性激活函数，利用预设的分离矩阵更新策略更新分离矩阵；计算更新前的分离矩阵与更新后的分离矩阵的抖动差值；比较所述抖动差值与预设的收敛判断因子的大小，若所述抖动差值小于或等于所述收敛判断因子时，根据更新后的分离矩阵计算所述分离信号，否则判断当前迭代次数是否达到变步长迭代阈值，若是，则更新步长因子，根据更新后的分离矩阵、更新后的步长因子和所述零均值预处理信号继续更新分离矩阵，否则根据更新后的分离矩阵、初始化的步长因子和所述零均值预处理信号继续更新分离矩阵。2CN113687314A权利要求书2/2页5.根据权利要求4所述的雷达抗主瓣干扰方法，其特征在于，所述分离矩阵更新策略的公式表示为：其中，k表示第k次迭代，W(k+1)表示第k+1次迭代的分离矩阵，即更新后的分离矩阵，W(k)表示第k次迭代的分离矩阵，即更新前的分离矩阵，μ(k)表示第k次迭代的步长因子，y(k)表示第k次迭代的观测信号，g(y)＝G′(y)表示为G(y)的一阶导函数，G(y)表示非线性激活函数，(·)T表示求转置，|·|表示求绝对值。6.根据权利要求5所述的雷达抗主瓣干扰方法，其特征在于，依照混合信号源的概率分布特性选取非线性激活函数，包括：若所述混合信号源为亚高斯分布时，选取的非线性激活函数表示为：若所述混合信号源为超高斯分布或者对算法的稳健性要求很高时，选取的非线性激活函数表示为：或其中，1≤a≤2。7.根据权利要求5所述的雷达抗主瓣干扰方法，其特征在于，所述对所述分离信号进行脉冲压缩处理公式表示为：其中，表示在t时刻脉冲压缩后的分离信号，表示在t时刻的分离信号，表示抖动差值小于或等于收敛判断因子时对应的分离矩阵，表示在t时刻的所述零均值预处理信号，s0(t0‑t)表示一已知的雷达发射信号s0(t)的匹配滤波器，conv(·)表示卷积运算，(·)*表示