基于循环对消的数字电视无源双基地雷达目标检测方法与流程

本发明涉及基于循环对消的数字电视无源双基地雷达目标检测方法，属于数字电视无源雷达技术领域。

背景技术：

无源雷达是指利用来自机会照射源的非合作信号，如调频广播(fm)、数字电视、模拟电视、手机信号、全球移动通信系统(gsm)等，以完全被动的方式探测移动目标的一种新体制雷达。相比于常用的fm、gsm等窄带信号照射源，数字电视照射源信号距离多普勒分辨率更高，发射功率更稳定，低空覆盖更广。

无源双基地雷达覆盖范围广，会面临大量片状运动杂波，如云雨杂波，洋流和涌波等。该类杂波与地杂波不同，不仅能量较强，可以遮盖低空小目标，而且由于多普勒效应的存在，归一化最小均方(nlms)、采样矩阵求逆(dmi)等简单的地杂波抑制算法不能抑制该类片状运动杂波。

技术实现要素：

针对现有技术存在的不足，本发明目的是提供一种基于循环对消的数字电视无源双基地雷达目标检测方法，在距离多普勒互相关处理之前有效地抑制直达波和片状运动杂波信号，从而实现弱目标检测。

为了实现上述目的，本发明是通过如下的技术方案来实现：

本发明的基于循环对消的数字电视无源双基地雷达目标检测方法，包括以下步骤：

(1)利用直达波抑制算法抑制直达波，估计强杂波的时延τ和多普勒频移f；重构各强杂波的回波信号；

(2)在监测通道中，通过时延多普勒滑动滤波器组循环对消算法，实现各强杂波的抑制；

(3)将对消结果与参考通道恢复的直达波信号再次进行距离多普勒互相关计算，从而实现弱目标的检测。

步骤(1)中，所述利用直达波抑制算法抑制直达波具体方法如下：

(1-1)在参考通道中，利用数字波束形成恢复直达波信号p(k)，并作为匹配滤波的参考信号，其中，k表示离散时间；

(1-2)在监测通道中，将监测通道接收到的信号xs(k)与恢复的直达波信号对消，将对消后的信号y(k)与参考通道恢复的直达波信号p(k)，进行距离与多普勒互相关计算。

步骤(1-1)中，所述恢复直达波信号p(k)的形成方法如下：

利用数字电视发射台位置信息，采用数字波束形成dbf实现发射站直达波恢复，波束输出公式如下：

式中：nc表示线阵天线单元数，d表示阵元间距，λ数字电视信号的波长，表示来波方向与阵列法线夹角；采用数字波束形成指向发射站，将dbf结果用作直达波信号。

步骤(1-2)中，所述对消结果y(k)具体的计算方法如下：

dmi滤波器的最优权计算公式为：

式中：rpp为直达波信号p(k)的协方差矩阵，[]^-1表示求逆运算，ppx为p(k)与监测通道接收信号xs(k)的互相关矩阵，rpp和ppx表达式如下：

其中，n1和n2分别为采样数据的起始点和终点，[·]^h表示共轭转置运算；

对消结果如下所示：

y(k)＝xs(k)-f^h(k)p(k)(8)

其中，xs(k)表示监测通道接收到的信号，p(k)表示dbf恢复的直达波信号。

步骤(1-2)中，对消后的信号y(k)与参考通道恢复的直达波信号p(k)，进行距离与多普勒互相关计算，具体方法如下：

其中，n是相关数据的长度，τ是距离延迟，f是多普勒频移；设置门限γ，当|ψ(τ,f)|＞γ时，记录τ和f的信息，此即估计的强杂波的时延和多普勒频移。

上述重构各强杂波的回波信号具体方法如下：

通过片状运动杂波检测算法获取由片状运动杂波分割的强杂波的个数na，即|ψ(τ,f)|超出门限时记录的多普勒的个数，根据f和对应该多普勒频移的最小时延τimin，重构各强杂波的回波信号，公式如下所示：

其中，fs表示采样率。

步骤(2)中，所述时延多普勒滑动滤波器组循环对消算法具体方法如下：

步骤1：初始化对消次数j＝0，输入对消前数据z0(k)，即y(k)，以及三个重构信号x1(k)，x2(k)和x3(k)；

步骤2：利用dmi计算最优权：

式中，rxixi为重构信号xi(k)的协方差矩阵，[]^-1表示求逆运算，pxiy为xi(k)与y(k)的互相关矩阵；

步骤3：计算加权后辅助通道数据：

步骤4：计算对消剩余结果：

zj(k)＝zj-1(k)-qj(k)(13)

步骤5：令j＝j+1，重复上述步骤直到j≥na/3，第j次对消结果为zj(k)，最终对消结果为zl(k)。

步骤(3)中，将对消结果zl(k)与参考通道恢复的直达波信号p(k)再次进行距离多普勒互相关处理，实现弱目标的检测，公式如下所示：

其中，设置门限γ，当|ψ(τ,f)|＞γ时，记录τ和f的信息，此即弱目标的时延和多普勒参数。

本发明在距离多普勒互相关处理之前有效地抑制直达波和片状运动杂波信号；先利用目标检测算法抑制直达波并估计强杂波的时延和多普勒频移(强杂波个数即片状运动杂波分割的个数)，重构各强杂波的回波信号；然后在监测通道中通过时延多普勒滑动滤波器组实现片状运动杂波的抑制；最后将对消后的信号与恢复的直达波信号进行距离多普勒互相关，从而实现弱目标检测。

附图说明

图1为新型目标检测框图；

图2为直达波抑制框图；

图3一次时延多普勒滑动滤波器组对消框图；

图4多级循环对消框图；

图5直达波对消前的距离多普勒图；

图6(a)为直达波对消后的距离多普勒立体图；

图6(b)为直达波对消后的多普勒剖面图；

图6(c)为直达波对消后的距离剖面图；

图7(a)为循环对消后的距离多普勒立体图；

图7(b)为循环对消后的多普勒剖面图；

图7(c)为循环对消后的距离剖面图。

具体实施方式

为使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解，下面结合具体实施方式，进一步阐述本发明。

1.数字电视外辐射源雷达杂波模型

假设数字电视无源双基地雷达的发射信号为p(k)，一个微小运动强杂波的多普勒频移为f，强杂波回波的衰减因子和时延分别为α0和τ0，那么该强杂波的回波s0(k)表示如下：

式中，k表示离散时间，θ0表示强杂波的入射角，fs表示采样率，a(θ0)表示该强杂波回波的阵列流形，可写成

其中，d表示阵元间距，λ表示数字电视信号的波长。

片状运动杂波可看作在时延和多普勒域连续的微小强杂波的集合，因此片状运动杂波的回波s(k)具有如下形式：

式中，n为s(k)在距离维上分割的个数，m为s(k)在多普勒维上分割的个数，θ表示入射角，a(θ)表示强杂波回波的阵列流形，τij表示强杂波回波的时延，fj表示多普勒频移。

2.基于循环对消的新型目标检测算法

本发明提出一种基于循环对消的新型目标检测算法，算法框图如图1所示。该算法先利用传统的目标检测算法抑制直达波从而检测片状运动杂波，估计强杂波的时延和多普勒频移(强杂波个数即片状运动杂波分割的个数)，重构各强杂波的回波信号。然后通过时延多普勒滑动滤波器组实现各强杂波的抑制，最后将对消后的信号与恢复的直达波信号进行距离多普勒互相关，从而获取弱目标的时延和多普勒信息。

(1)直达波信号恢复

利用已知的数字电视发射台位置信息，采用数字波束形成技术(dbf)实现发射站直达波恢复，波束输出公式如下：

式中：nc表示线阵天线单元数，d表示阵元间距，λ数字电视信号的波长，表示来波方向与阵列法线夹角。采用数字波束形成指向发射站，将dbf结果用作直达波信号。

(2)直达波抑制算法

数字电视无源双基地雷达系统需要具备两个通道：实现直达波恢复的参考通道与实现直达波抑制与片状运动杂波检测的监测通道。如图2所示，在参考通道中，利用数字波束形成恢复直达波信号p(k)，用作后续匹配滤波的参考信号。在监测通道中，将接收到的信号xs(k)(包括直达波信号、多径信号和噪声)与恢复的直达波信号对消，对消后的信号y(k)通过与直达波信号p(k)做距离多普勒互相关处理，提取出片状运动杂波的时延和多普勒信息，从而实现片状运动杂波的检测。

选择采样矩阵求逆算法(dmi)完成直达波的消除。dmi滤波器的最优权计算公式为：

式中：rpp为直达波信号p(k)的协方差矩阵，[]^-1表示求逆运算，ppx为p(k)与监测通道接收信号xs(k)的互相关矩阵，rpp和ppx表达式如下：

其中，n1和n2分别为采样数据的起始点和终点，[]^h表示共轭转置运算。

对消结果如下所示：

y(k)＝xs(k)-f^h(k)p(k).(8)

其中，xs(k)表示监测通道接收到的信号，p(k)表示dbf恢复的直达波信号。

(3)基于距离多普勒互相关的时延与多普勒参数估计

针对对消后的信号与参考通道恢复的直达波信号，进行距离与多普勒互相关计算，计算式如下

其中，y(k)、p(k)分别为对消剩余和参考通道恢复的直达波信号，τ是距离延迟，f是多普勒频移。设置门限γ，当|ψ(τ,f)|＞γ时，记录τ和f的信息，此即估计的强杂波的时延和多普勒参数。

(4)强杂波重构

通过片状运动杂波检测算法获取由片状运动杂波分割的强杂波的个数na(|ψ(τ,f)|超出门限时记录的多普勒的个数)，各强杂波的时延τ以及多普勒频移f，根据f和对应该多普勒频移的最小时延τimin，重构各强杂波的回波信号，公式如下所示：

(5)时延多普勒滑动滤波器组循环对消算法

考虑到对消算法的复杂度和收敛性，每三个重构信号作为一组辅助通道，通过时延多普勒滑动滤波器与上一次对消结果进行对消，抑制强杂波信号。一次时延多普勒滑动滤波器组对消框图如图3所示。循环时延多普勒滑动对消na/3次，多级循环对消框图见图4所示。对消步骤如下所示：

步骤1：初始化对消次数j＝0，输入对消前数据z0(k)，即y(k)，以及三个重构信号x1(k)，x2(k)和x3(k)。

步骤2：利用dmi计算最优权：

式中，为重构信号xi(k)的协方差矩阵，[]-1表示求逆运算，pxiy为xi(k)与y(k)的互相关矩阵。

步骤3：计算加权后辅助通道数据：

步骤4：计算对消剩余结果：

zj(k)＝zj-1(k)-qj(k).(13)

步骤5：令j＝j+1，重复上述步骤直到j≥na/3，第j次对消结果为zj(k)，最终对消结果为zl(k)。

(6)基于距离多普勒互相关的弱目标检测

最终对消结果zl(k)与参考通道恢复的直达波信号p(k)再次进行距离多普勒互相关处理，实现弱目标的检测，公式如下所示：

其中，τ是距离延迟，f是多普勒频移。设置门限γ，当|ψ(τ,f)|＞γ时，记录τ和f的信息，此即弱目标的时延和多普勒参数。

通过计算机仿真验证循环对消器的性能。仿真中设置采样率为300khz，产生1个片状运动杂波和1个弱目标，弱目标的距离和多普勒分别是102km和140hz，片状运动杂波的距离在60～65km区间内变化，多普勒在20～80hz区间内变化。直达波入射角为15°，弱目标入射角为40°，接收阵列采用25单元的半波长等间隔线阵，指向角为30°。直达波对消前的距离多普勒图如图5所示，目标被片状运动杂波和直达波的副瓣所掩盖。

直达波抑制后的距离多普勒图如图6(a)-(c)所示，直达波被完全抑制，噪声基底由图5的163.5db降低为151.7db，可以更好地检测出片状运动杂波，而弱目标被片状运动杂波的旁瓣掩盖，无法被检测到。图6(b)中幅度与门限180db相比，参数估计算法测得的片状运动杂波的多普勒为21～80hz，共60种不同的多普勒频率，后续经过20(60/3)次循环对消。对消后的距离多普勒图如图7(a)-(c)所示，在对片状运动杂波进行抑制后，弱目标的信噪比由图6的22.3db提升为36.3db，检测到弱目标的距离和多普勒分别为102km和140hz，与设置一致，可以较好地实现弱目标检测。仿真实验结果表明了该方法的有效性。

以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。