一种基于子空间投影的外辐射源雷达对目标压缩检测方法与流程

本发明属于信号处理技术领域，涉及一种基于子空间投影的外辐射源雷达对目标压缩检测方法。

背景技术：

近年来，以广播信号、电视信号、卫星信号等作为照射源的外辐射源雷达得到了广泛的研究，其优点在于成本低、只接收信号不发射电磁波具有很好的隐蔽特性、照射源多样性等。随着雷达信号处理技术的发展，对雷达信号处理正沿着高宽带到分辨率的方向努力，但同时也面临着要处理更多的数据量。为了减少数据运算量，压缩感知理论的提出证明了在明显减少了对采样数据量的分析的情况下，基本不丢失压缩之前信号的特征。目前雷达系统中采用的压缩感知技术，通常需要利用重构算法对信号进行重构，进而实现目标检测与参数估计。但是，重构算法一般要求较高的信噪比，才能得到满足目标检测与参数估计的重构信号。对于像外辐射源雷达由于其本质是双基地非合作雷达，目标回波功率低，通常不满足信号重构的信噪比要求。因此采用压缩感知技术实现微弱信号检测需要新方法。

直达波对消技术是解决外辐射源雷达目标检测的前提。在研究直达波对消的问题上，2013年赵艳等人提出了基于自适应对消的直达波抑制技术，2008年陈多芳等人利用了子空间投影方法实现了直达波对消并完成目标检测。而基于压缩感知的信号检测的研究主要有以下几种方法。2014年王康等人提出稀疏系数幅度比较算法，其核心思想是对用两个观测测量矩阵同时接收待检测信号，利用omp算法得到的幅度和位置信息联合设计判决门限，文中提出门限系数的设置是根据经验取值，没有给出具体确定计算方法；2015年alirezahariri等人提出最大似然比值算法，其核心思想是没有对待检测信号重构，得到两种情况的似然函数，然后做比值，通过似然比函数概率密度特性得到门限公式，进行判决。但是，上述三种检测算法在信噪比为-5db以下就不能实现目标检测。这就限制了检测器的应用范围，为了探测更大的范围，就需要解决在更低信噪比下对待检测信号的检测。2017年马俊虎等人根据目标在距离维上的稀疏性，设计了压缩测量矩阵，直接对目标完成检测，无需信号重构，但是其是基于回波中的直达波对消之后而实现对目标的检测。本发明在此基础上考虑了回波信号中的直达波和多径影响，设计新的测量矩阵完成对目标的检测，无需信号重构，减少了数据量。

技术实现要素：

本发明所要解决的技术问题是对压缩后的外辐射源雷达回波信号在没有信号重构的条件下直接对目标进行检测。

为了便于理解，对本发明采用的技术做如下说明：

由于接收的目标回波信号中含有很强的直达波和多径杂波，并且直达波所在空间多普勒频移为零，而多径杂波认为服从威布尔分布。将目标回波投影到直达波正交补子空间，由于存在多普勒频移，在该空间中信号被保存了下来，进而可以实现直达波的对消。进一步，要构造外辐射源雷达系统中直达波等无多普勒频移的正交子空间信号矩阵，需要利用辅助信道中接收的参考信号sr∈r^n×1，则直达波扰矩阵可如下表示：

j＝[j1j2...jl]＝[srdsr...d^l-2sr]∈r^n×l

其中，n为信号长度，d表示的是时延矩阵，表达式如下所示：

那么，主天线通道回波信号中的直达波子空间的投影矩阵：

p0＝j(j^hj)^-1j^h∈r^n×n

并且可以得到主天线中的直达波正交补子空间的投影矩阵为：

p0^⊥＝i-p0＝i-j(j^hj)^-1j^h∈r^n×n

所以将回波信号投影到p0^⊥上可以实现直达波的对消。目标在雷达探测范围内的距离维和速度维上可认为是稀疏的。设在距离维度上其稀疏字典基：

ψ＝[sr(t-τ1)sr(t-τ2)…sr(t-τn)]∈r^n×n

其中τi为距离上每个距离单元对应的时延。因此可以设计测量矩阵：

其中压缩比为正整数。对压缩测量的回波信号取模平方输出判决目标是否存在。

为解决上述问题，本发明的技术方案为：

一种基于压缩感知的外辐射源雷达对目标检测方法，该方法用于对压缩后的信号在没有信号重构的条件下直接对信号进行检测，其特征在于，包括以下步骤：

确定检测门限，具体包括：

s1、首先，根据参考信道接收的参考信号sr(t)构造距离维上能稀疏表示目标的字典基矩阵ψ＝[[sr(t-τ1)sr(t-τ2)…sr(t-τn)]∈r^n×n，其中τi为距离上每个距离单元对应的时延；

s2、取ψ的前l列组成子空间j＝[j1j2...jl]，计算得到：p0^⊥＝i-j(j^hj)-¹j^h，设计测量矩阵

s3、对回波信号压缩测量y＝φsh(t)，并取模平方得：l＝|y|²；

s4、把数据长度为m的l分成n个向量。如果li>γ，则判为有信号；否则判为信号不存在，其中i＝1,…,n，γ为检测门限，

本发明的有益效果为，首先根据目标在距离维和速度维上具有稀疏性，然后根据参考信号构造稀疏表示目标回波的字典基，再结合子空间投影思想，设计了确定性测量矩阵，实现了对回波信号的直达波对消并直接对信号进行检测，无需信号重构，并且减少了数据运算量。在待检测探测方向具有很好的参考和实际应用。

附图说明

图1本发明方法对待检测信号检测流程图；

图2检测器测量效果

图3不同压缩比在不同信噪比下检测性能。

具体实施方式

下面将结合附图和实施例，对本发明的技术方案进行进一步说明。

实施例

本例在接收待检测信号中直达波干扰dir＝20db，信杂比sir＝-10db，信噪比snr＝-15db到-6db步进1db时对信号的存在与否进行检测：

实施例的检测实施方法如附图1、附图2所示。具体步骤如下所示：

确定检测判决门限：

试验中选择本地雷达信号为：回波信号表达式为：

其中，a1＝0.1，τi为时延，fd为多普勒频移，c(t)为杂波，本实验考虑杂波为韦布尔杂波，m(t)为高斯白噪声。仿真参数设置：设雷达探测距离rmin＝10000m，rmax＝15000m，中心频率f0＝100mhz，信号带宽b＝30mhz，时宽t＝10μs，采样频率fs＝90mhz，调频斜率为：μ＝b/t＝3×10¹²hz/s。本次实验设置1个目标，在r1＝14000m，多普勒频移fd＝200hz，回波信号的稀疏字典基ψ＝[sr(t-τ1)sr(t-τ2)…sr(t-τn)]^t，其中τi为距离上每个距离单元对应的时延。

步骤一：首先，根据参考信道接收的参考信号sr(t)构造距离维上能稀疏表示目标的字典基矩阵ψ＝[sr(t-τ1)sr(t-τ2)…sr(t-τn)]，其中τi为距离上每个距离单元对应的时延；

步骤二：取ψ的前l列组成子空间j＝[j1j2...jl]，计算得到：p0^⊥＝i-j(j^hj)^-1j^h，设计测量矩阵

步骤三：对回波信号压缩测量y＝φsh(t)，并取模平方得：l＝|y|²；

步骤四：把数据长度为m的l分成32个向量。如果li>γ，则判为有信号；否则判为信号不存在，其中i＝1,…,n，γ为检测门限，

图2是在直达波干扰dir＝20db，信杂比sir＝-10db，信噪比snr＝-10db，m/n＝0.5时，压缩信号与门限比较结果。从图2看出该算法能够在低信噪比时，没有信号重构的情况下能够很好的消除直达波的干扰，并且明显观察出目标所在位置。图2说明了本算法具有较好的降噪性能。

图3是在直达波干扰dir＝20db，信杂比sir＝-10db时，不同压缩比下检测器工作的性能曲线，结果可以看出当压缩比m/n＝0.5，信噪比为-15db时检测概率为97％，可以很精确的判决出待检测信号的待检测是否存在。当压缩比m/n＝0.2，信噪比为-9db时检测概率为95％，可以很精确的判决出待检测信号的待检测是否存在。