基于低秩稀疏矩阵约束优化的叶簇杂波抑制方法及系统

1.本发明属于雷达信号处理领域，具体涉及一种基于低秩稀疏矩阵约束优化的叶簇杂波抑制方法及系统。


背景技术：

2.地面战场侦察雷达是现代信息化战争中战场态势感知和情报获取的重要手段。地面战场侦察雷达在地面运动目标检测中，由于地杂波的干扰，低速目标检测是其面临的一大难题。在和平年代的非对称战争和反恐战争中，常常面临着检测处于树林遮蔽中的低速运动目标的需求，因战场侦察雷达的工作频段通常是x波段或者更高，其不具有绕射能力，致使低速运动目标直接淹没在叶簇杂波中，无法检测。因此，如何有效抑制叶簇杂波对提升地面战场侦察雷达的低速目标检测性能具有非常重要的意义。目前，在地面动目标检测中，为了减少叶簇杂波的干扰，人们提出了一系列的杂波抑制方法：对于目标多普勒都在主瓣杂波之外的情况，利用杂波与运动目标之间的频率差异提出了运动目标指示(mti)、运动目标检测(mtd)在频域对杂波进行滤波；使用子回归模型对叶簇杂波进行建模，并采用奇异值分解(svd)进行杂波剔除和目标信号重构；但以上方法对于低速目标监视场景具有很大的局限性。
3.专利申请号为cn 201610060416.1的发明专利公开了一种基于小波变换的船用雷达海杂波抑制曲线设计方法，该方法使用离散小波变换对雷达回波数据进行分解来抑制背景海杂波。该方法虽有在海杂波抑制方面取得了一定的成果，但对于地面战场侦察雷达在低速目标监视场景具有很大的局限性。


技术实现要素：

4.本发明的目的在于提出一种基于低秩稀疏矩阵约束优化的地面战场侦察雷达叶簇杂波抑制方法及系统。
5.实现本发明目的的技术解决方案为：一种基于低秩稀疏矩阵约束优化的叶簇杂波抑制方法，步骤如下：
6.步骤1、对雷达回波信号进行不同阶数的frft变换，根据二维分布在平面中找出能量最大点；
7.步骤2、由步骤1求出的能量最大点确定出最优变换阶次ρu，同时对雷达回波进行ρu阶frft变换；
8.步骤3、在frft域将雷达回波矩阵进行稀疏分解建立优化模型，并对优化模型进行凸松弛；
9.步骤4、基于改进增广拉格朗日乘子法求解目标函数，初始化相关值后不停迭代直至收敛得到目标函数矩阵；
10.步骤5、将步骤4中重构的信号逆fft变换到时域，得到抑制杂波后的数据。
11.一种基于低秩稀疏矩阵约束优化的叶簇杂波抑制系统，包括：
12.第一模块，用于对雷达回波信号进行不同阶数的frft变换，根据二维分布在平面中找出能量最大点；
13.第二模块，通过第一模块找出的能量最大点确定出最优变换阶次ρu，同时对雷达回波进行ρu阶frft变换；
14.第三模块，在frft域将雷达回波矩阵进行稀疏分解建立优化模型，并对优化模型进行凸松弛；
15.第四模块，基于改进增广拉格朗日乘子法求解目标函数，初始化相关值后不停迭代直至收敛得到目标函数矩阵；
16.第五模块，用于将第四模块重构的信号逆fft变换到时域，得到抑制杂波后的数据。
17.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现上述的基于低秩稀疏矩阵约束优化的叶簇杂波抑制方法。
18.一种计算机可读存储介质，其上存储有计算机程序，程序被处理器执行时实现上述的基于低秩稀疏矩阵约束优化的叶簇杂波抑制方法。
19.本发明与现有技术相比，其显著优点在于：
20.(1)基于sa算法改进增广拉格朗日乘子法保证每一轮迭代中都能获取最优的μ
sa
，保证了算法的自适应能力，提高目标信号重构的完整度。
21.(2)提出在频域进行低秩稀疏约束重构，从而实现叶簇杂波抑制，提高低速目标回波信号的信杂比，有效提升了叶簇杂波中低速目标的检测性能。
附图说明
22.图1本发明基于低秩稀疏矩阵约束优化的地面战场侦察雷达叶簇杂波抑制方法流程图。
23.图2为本发明提出的基于模拟退火优化(simulated annealing，sa)改进增广拉格朗日乘子法算法流程图。
24.图3为雷达回波功率谱示意图。
25.图4为采用低秩稀疏联合约束对叶簇杂波抑制结果示意图。
具体实施方式
26.本发明提出一种基于低秩稀疏矩阵约束优化的叶簇杂波抑制方法，该方法首先对雷达回波进行分数阶傅里叶变换并根据二维分布找出平面内能量最大点，确定最佳变换阶数ρ，再对雷达回波进行ρ阶frft变换，使得雷达接收信号在frft域中目标和杂波的能量差取到最大。接着在frft域构造低秩稀疏约束模型并优化目标函数，通过基于模拟退火优化(simulated annealing，sa)改进增广拉格朗日乘子法迭代求解，最后将重构的信号逆fft变换到时域，得到抑制杂波后的数据。本发明针对叶簇杂波环境中地面侦察雷达检测低速运动目标的难题，在频域进行低秩稀疏约束重构，从而实现叶簇杂波抑制，提高低速目标回波信号的信杂比，有效提升了树林环境中低速目标的检测性能。
27.如图1所示，本发明基于低秩稀疏矩阵约束优化的地面战场侦察雷达叶簇杂波抑
制方法包括以下步骤：
28.步骤1、对不同阶数ρ下的雷达接收数据根据公式进行frft变换,得到所有阶次下雷达信号frft域的幅度分布。并根据二维分布在平面中找出能量最大点坐标。
29.其中，为旋转角度，ρ为frft的变换阶数。
30.步骤2、由步骤1求出的能量最大点确定出最优变换阶次ρu，对雷达接收回波进行ρu阶frft变换，使得在frft域中目标和杂波的能量差异取到最大。
31.步骤3、基于rpca在frft域将雷达回波矩阵进行稀疏分解建立优化模型，叶簇杂波的抑制可表述为低秩稀疏约束优化问题为：
[0032][0033]
其中，z
t
为目标子空间；zf为杂波和噪声子空间；rank(
·
)为矩阵的秩；||
·
||0为矩阵的0阶范数；λ为正则化参数，与分解精度相关。引用凸松弛求解上述优化问题，将范数松弛为范数，同时以核范数替换rank(
·
)进行低秩约束：
[0034][0035]
式中，||
·
||
*
为矩阵的核范数；||
·
||1为矩阵的1阶范数；λ为正则化参数与分解精度相关。
[0036]
步骤4、如图2所示，引入增广拉格朗日乘子法(alm)进行低秩求解：
[0037][0038]
式中，l(zf,z
t
,y)为增广拉格朗日函数；λ为正则化参数；μ＞0为迭代系数，用来平衡迭代次数和终止条件；y为拉格朗日乘子；||
·
||f为矩阵的frobenius范数；
[0039]
由于迭代系数μ自适应能力较弱，因而引入sa算法，以此在每轮的迭代中都能获取最优的μ，表示为μ
sa
，因此可以得到:
[0040][0041]
式中，μ
sa
为采用sa算法求出的每轮最优参数，在矩阵每次迭代时自适应更新。
[0042]
具体更新方法为：设定μ的取值范围是(其中z
fk
为第k次迭代的目标矩阵，σ
max
(z
fk
)为z
fk
的最大奇异值)，取适应度函数为(其中σi(z
fk
)为z
fk
的第i个奇异值)，基于sa算法得到最优的惩罚系数μ
sa
。
[0043]
假设最大迭代次数为k，则第k次运算过程分为两步：首先求出本次迭代的最优参数再将作为本次迭代参数更新，具体步骤为：
[0044]
定义迭代算子s
τ
(x)＝sgn(x)max(|x|-τ,0)，并带入可得：
[0045][0046]
式中，x为复数；τ为实数。
[0047]
再次定义奇异值迭代算子d
τ
(zf)＝us
τ
(∑)v
*
带入上式可得：
[0048][0049]
先固定z、z
t
、y、μ
sa
，则可得到zf的迭代更新公式为：
[0050][0051]
类似的固定z、zf、y、μ
sa
，则可得到z
t
的迭代更新公式为：
[0052][0053]
最后固定z、zf、z
t
、μ
sa
求得y的迭代更新公式为：
[0054][0055]
多轮迭代后，直至满足收敛条件(ε通常取10-4
)或迭代次数达到k
max
时，此时得到目标信号矩阵。
[0056]
步骤5、将步骤4中重构的信号使用逆频fft变换到时域，从得到所需要的抑制杂波的数据。
[0057]
本发明还提供一种基于低秩稀疏矩阵约束优化的叶簇杂波抑制系统，包括：
[0058]
第一模块，用于对雷达回波信号进行不同阶数的frft变换，根据二维分布在平面中找出能量最大点；
[0059]
第二模块，通过第一模块找出的能量最大点确定出最优变换阶次ρu，同时对雷达回波进行ρu阶frft变换；
[0060]
第三模块，在frft域将雷达回波矩阵进行稀疏分解建立优化模型，并对优化模型进行凸松弛；
[0061]
第四模块，基于改进增广拉格朗日乘子法求解目标函数，初始化相关值后不停迭代直至收敛得到目标函数矩阵；
[0062]
第五模块，用于将第四模块重构的信号逆fft变换到时域，得到抑制杂波后的数据。
[0063]
上述第一～第五模块的具体实现方法与前述的基于低秩稀疏矩阵约束优化的叶簇杂波抑制方法的各步骤具体方法相同，此处不再赘述。
[0064]
下面进行仿真验证本发明的杂波抑制算法
[0065]
雷达仿真参数设置为：雷达工作频率10ghz，信号波形采用线性调频信号，信号带宽为30mhz，时宽为0.5us，重复周期为25us，采样频率为240mhz。目标位于第200距离门，以0.5m/s的速度匀速临近雷达。回波信号积累周期为512。杂波信号采用单摆多径杂波模型进行仿真。
[0066]
采用本发明基础低秩稀疏联合约束方法对上述叶簇杂波进行抑制的结果如图3和图4所示。从图4中可以看出，经过本发明基于低秩稀疏矩阵约束优化的杂波抑制方法处理后，目标信号可以实现较好的重构，最大杂波功率谱相对于目标谱衰减达15db左右，可以有效提升战场侦查雷达叶簇环境下低速目标的检测性能。
[0067]
以上实施例仅为说明本发明的技术思想，不能以此限定本发明的保护范围，凡是按照本发明提出的技术思想，在技术方案基础上所做的任何改动，均落入本发明保护范围
之内；本发明未涉及的技术均可通过现有技术加以实现。