基于动态有序矩阵的外辐射源雷达CFAR算法

饶云华; 周健康; 万显荣; 龚子平; 柯亨玉

doi:10.11999/JEIT191024

基于动态有序矩阵的外辐射源雷达CFAR算法

doi: 10.11999/JEIT191024

饶云华^{1, 2, ,},
周健康¹,
万显荣^{1, 2},
龚子平^{1, 2},
柯亨玉^{1, 2}

1.
武汉大学电子信息学院武汉 430072
2.
武汉大学深圳研究院深圳 518063

基金项目: 国家自然科学基金(U1933135, 61271400)，国家重点研发计划(2016YFB0502403)，湖北省技术创新专项重大项目(2016AAA017)，深圳市科技计划项目(JCYJ20170818112037398)

详细信息

作者简介:
饶云华：男，1972年生，副教授，研究方向为新体制雷达、雷达系统设计、无线通信网等

周健康：男，1993年生，硕士生，研究方向为雷达电磁环境辨识、恒虚警检测

万显荣：男，1975年生，教授，研究方向为无源雷达、超视距雷达、新体制雷达系统与雷达信号处理等

龚子平：男，1977年生，讲师，研究方向为电波传播与无线电海洋遥感等

柯亨玉：男，1957年生，教授，研究方向为电磁场理论、高频雷达海洋遥感技术

通讯作者:
饶云华　ryh@whu.edu.cn

中图分类号: TN957.51
计量
- 文章访问数: 1284
- HTML全文浏览量: 628
- PDF下载量: 118
- 被引次数: 14
出版历程
- 收稿日期: 2019-12-23
- 修回日期: 2020-10-20
- 网络出版日期: 2020-12-08
- 刊出日期: 2021-04-20

CFAR for Passive Radar Based on Dynamic Ordered Matrix

Yunhua RAO^{1, 2
, ,},
Jiankang ZHOU¹,
Xianrong WAN^{1, 2},
Ziping GONG^{1, 2},
Hengyu KE^{1, 2}

1.
School of Electronic Information, Wuhan University, Wuhan 430072, China
2.
Shenzhen Research Institute, Wuhan University, Shenzhen 518063, China

Funds: The National Natural Science Foundation of China (U1933135,61271400), The National Key Research and Development Project (2016YFB0502403), The Hubei Province Technology Innovation Special Major Project (2016AAA017), The Shenzhen Science and Technology Project (JCYJ20170818112037398)

摘要

摘要: 外辐射源雷达采用不可控的第三方辐射源，其电磁传播条件复杂，尤其是在低空目标探测中，检测性能极大地受到杂波特性的影响，使得传统恒虚警算法性能明显下降。为了改善检测性能，该文提出一种基于雷达杂波空间划分的动态有序矩阵恒虚警检测算法(DOM-CFAR)。该算法将杂波空间从距离和多普勒维进行划分，构造为有序矩阵，再根据背景杂波变化进行动态极值替换、提取杂波估计中值用以计算检测阈值，从而使得检测算法阈值可动态适应杂波功率变化。仿真和实测结果表明，该算法可以在均匀杂波、多目标和杂波边缘等复杂情况下保持稳定的检测性能。
- 目标检测 /
- 杂波矩阵 /
- 杂波图 /
- 恒虚警
Abstract: Passive radar uses third party radiation source, which is uncontrollable. Due to the complicated electromagnetic propagation conditions, especially in a low altitude target detection, the detection performance of the radar is greatly affected by clutters, leading to significant degradation of the performance of traditional constant false alarm algorithm. In order to improve the detection performance, a Dynamic Ordered Matrix Constant False Alarm Rate (DOM-CFAR) algorithm based on radar clutter space partition is proposed. In this algorithm, an ordered matrix is constructed after dividing the clutter space from distance and Doppler dimension. Then the dynamic extreme value is replaced according to the background clutter change and the estimated median value of clutter is extracted so as to calculate the detection threshold. This algorithm makes the threshold value of detection can dynamically adapt to the clutter power change. The simulation and measurement results show that the algorithm can maintain excellent detection performance under the complex environment with uniform clutter, multi-target and clutter edge.
- Target detection /
- Clutter matrix /
- Clutter map /
- Constant False Alarm Rate (CFAR)

HTML全文

1. 引言

在外辐射源雷达低空小目标探测中，其背景杂波受传播环境影响而复杂多变，尤其是功率突变的杂波会在距离多普勒谱上形成点状、块状的杂波边缘，极大地影响了低空小目标的恒虚警检测性能。

当前常用的恒虚警算法根据杂波在空间与时序上的分布特性与均匀性可划分为：空域类(也称为滑窗类)和时序类。

空域类恒虚警算法中具有代表性的是单元平均(Cell Averaging)CA-CFAR算法^[1,2]，其采用临近待检测单元的杂波均值来计算检测阈值，在均匀环境下效果较好，但是在多目标及杂波边缘环境下无法消除临近目标与能量突变干扰，使得检测性能恶化。在该算法的基础上，Trunk^[3]提出了最小选择(Smallest Of)SO-CFAR以解决多目标干扰问题，Hansen等人^[4]提出最大选择(Greatest Of)GO-CFAR以解决杂波边缘环境下的能量突变问题，但是这两种算法都在均匀环境下存在一定性能损失，且并不适用于另一种非均匀环境。Smith等人^[5]在以上算法的基础上提出了变指数(Varibility Index)VI-CFAR。该算法先判决当前杂波环境的类型，然后选择与其匹配的其他恒虚警算法^[6,7]，其本质是多种算法的组合。Rohling等人^[8-11]提出有序统计(Order Statistics)OS-CFAR，虽然该算法在多目标、杂波边缘环境下性能得到部分改善，但在均匀环境下性能也有一定损失。

时序类恒虚警算法^[12-16]中，杂波图(Clutter Map)CM-CFAR是该类算法的典型代表，CM-CFAR将杂波空间划分为杂波图单元用以存储杂波数据，每个杂波图单元依靠当前数据和历史数据进行加权迭代更新。由于这类算法在迭代初期采样的数据量少，所以检测性能较差；随着迭代次数增加，检测性能会得到较大的提升。但是这类算法也具有因迭代计算而带来的记忆效应这一不可忽视的缺点。

在外辐射源雷达低空小目标探测中，由于目标回波弱、背景杂波复杂多变，空域类算法在该类应用中性能损失严重、虚警概率恶化，传统时序类算法迭代初期性能较差并且还有记忆效应等问题，为此本文提出了动态有序矩阵(Dynamic Order Matrix)DOM-CFAR算法，其将杂波空间划分成有序矩阵，利用雷达初次扫描数据的滑窗中值对有序矩阵中的有序数组进行初始化用以弥补传统时序类算法在迭代初期的性能劣势，并通过雷达扫描数据对有序数组进行动态极值替换用以获取杂波估计中值并弥补记忆效应造成的影响，最后利用杂波估计中值来计算检测阈值。

2. 检测模型

2.1 杂波与目标

在判决雷达回波 $x(t)$ 中是否存在目标时，根据二元假设检测可知

$\left. \begin{aligned} & x(t) = n(t),\qquad\;\;\;{{\rm{H}}_0} \\ & x(t) = n(t) + s(t),{{\rm{H}}_1} \end{aligned} \right\}$

(1)

其中， $s(t)$ 为目标信号， $n(t)$ 为杂波和噪声信号， ${{\rm{H}}_{\rm{0}}}$ 表示只有杂波， ${{\rm{H}}_{\rm{1}}}$ 表示存在目标。

杂波分布模型按照其概率密度函数参数个数可以分为单参和多参分布模型^[17-20]。例如瑞利分布、指数分布等只由1个参数决定其幅度分布规律的为单参分布模型；韦布尔分布、对数正态分布、莱斯分布等由多个参数决定其幅度分布规律的为多参分布模型。本文以杂波数据经平方率检波后服从指数分布为背景，其杂波概率密度函数(Probability Density Function, PDF)为

$f(x) = \frac{1}{\lambda }\exp \left( - \frac{x}{\lambda }\right),0 \le x$

(2)

其中， $\lambda$ 为每个杂波单元总功率。在 ${{\rm{H}}_{\rm{0}}}$ 假设下， $\lambda$ 为杂波和噪声的功率 $\mu$ ；在 ${{\rm{H}}_{\rm{1}}}$ 假设下， $\lambda$ 为杂波、噪声和目标的总功率 $(1 + s)\mu$ , $s$ 表示目标与杂波噪声功率 $\mu$ 的比值。由式(2)可得检测概率(Probability of Detection, PD)与虚警概率(Probability of False Alarm, PFA)为

$\begin{split} {\rm{PD}} =\,& {\rm{P}} [x > {x_{\rm{t}}}|{{\rm{H}}_1}] = \int_{{x_{\rm{t}}}}^\infty {f(x){\rm{d}} x} |{{\rm{H}}_1}\ \ \\ =\,& \exp \left( { - \frac{{{x_{\rm{t}}}}}{{(1 + s)\mu }}} \right) \end{split}$

(3)

$\begin{split} {\rm{PFA}} =\,& {\rm{P}} [x > {x_{\rm{t}}}|{{\rm{H}}_{\rm{0}}}] = \int_{{x_{\rm{t}}}}^\infty {f(x){\rm{d}} x} |{{\rm{H}}_{\rm{0}}} \\ =\,& \exp \left( { - \frac{{{x_{\rm{t}}}}}{\mu }} \right) \end{split}$

(4)

在已知参数 $\mu$ 和预设虚警概率 ${\rm{PFA}}$ 时，可由式(4)求得检测阈值T为

$T = - \mu \ln ({\rm{PFA}})$

(5)

2.2 有序矩阵

对于雷达回波的RD谱数据 ${{X}}$ ，每帧数据既有距离信息也有多普勒信息，其表达式为

${{X}} = \left[ {\begin{array}{*{20}{c}} {{x_{1,1}}}&{{x_{1,2}}}& ··· &{{x_{1,d}}}& ··· &{{x_{1,D}}}\\ {{x_{2,1}}}&{{x_{2,2}}}& ··· &{{x_{2,d}}}& ··· &{{x_{2,D}}}\\ \vdots & \vdots &{\ddots}& \vdots &{\ddots}& \vdots \\ {{x_{r,1}}}&{{x_{r,2}}}& ··· &{{x_{r,d}}}& ··· &{{x_{r,D}}}\\ \vdots & \vdots &{\ddots}& \vdots &{\ddots}& \vdots \\ {{x_{R,1}}}&{{x_{R,2}}}& ··· &{{x_{R,d}}}& ··· &{{x_{R,D}}} \end{array}} \right]$

(6)

其中， ${x_{r,d}}$ 表示RD谱 ${{X}}$ 中待检测的第 $(r,d)$ 单元数据，其中 $r$ 为第 $r$ 个距离元、 $d$ 为第 $d$ 个多普勒元。由于本文所提算法利用待检测杂波的历史数据来估计中值 ${y_{{\rm{mid}}}}$ 求解检测阈值T，为了准确地估算估计中值 ${y_{{\rm{mid}}}}$ ，这里利用有序矩阵 ${{M}}$ 来存储相关历史数据，有序矩阵 ${{M}}$ 的表达式为

${{M}} = \left[ {\begin{array}{*{20}{c}} {{{{Y}}_{1,1}}}&{{{{Y}}_{1,2}}}& ··· &{{{{Y}}_{1,d}}}& ··· &{{{{Y}}_{1,D}}}\\ {{{{Y}}_{2,1}}}&{{{{Y}}_{2,2}}}& ··· &{{{{Y}}_{2,d}}}& ··· &{{{{Y}}_{2,D}}}\\ \vdots & \vdots &{\ddots}& \vdots &{\ddots}& \vdots \\ {{{{Y}}_{r,1}}}&{{{{Y}}_{r,2}}}& ··· &{{{{Y}}_{r,d}}}& ··· &{{{{Y}}_{r,D}}}\\ \vdots & \vdots &{\ddots}& \vdots &{\ddots}& \vdots \\ {{{{Y}}_{R,1}}}&{{{{Y}}_{R,2}}}& ··· &{{{{Y}}_{R,d}}}& ··· &{{{{Y}}_{R,D}}} \end{array}} \right]$

(7)

其中， ${{{Y}}_{r,d}}$ 为有序矩阵 ${{M}}$ 中第 $(r,d)$ 单元的有序数组，其长度为 $2K + 1$ ，用于存储 ${x_{r,d}}$ 历史数据，表达式为

${{{Y}}_{r,d}} = \left[ {\begin{array}{*{20}{c}} {{y_1}}\;\;\;{{y_2}}\;\;\; ···\;\;\; {{y_k}}\;\;\; ···\;\;\; {{y_{2K + 1}}} \end{array} } \right]$

(8)

其中， ${y_k}$ 为有序数组 ${{{Y}}_{r,d}}$ 中的第 $k$ 大的数据， $K$ 为本算法的半参考窗值，K值越大，所需的存储空间越多。待检测单元 ${x_{r,d}}$ 的杂波估计中值 ${y_{{\rm{mid}}}}$ 为 ${{{Y}}_{r,d}}$ 中的 ${y_{K + 1}}$ 。

3. DOM-CFAR原理

DOM-CFAR算法实现流程如图1所示，分为初始化、迭代优化与检测3个步骤。

图 1 DOM-CFAR算法实现流程图

下载: 全尺寸图片幻灯片

3.1 初始化

初始化时，以RD谱数据 ${{X}}$ 中待检测单元 ${x_{r,d}}$ 为中心的长度为 $2K + 1$ 的距离多普勒域滑窗进行排序整理，用滑窗中值 $x{'_{r,d}}$ 对有序数组 ${{{Y}}_{r,d}}$ 中的所有数据赋初值，如式(9)所示

${y_k} = x{'_{r,d}},k \in [1,2K + 1]$

(9)

初始化后，有序数组 ${{{Y}}_{r,d}}$ 中的数据 ${y_k}$ 服从 $N({u_0} + \tau ,{\sigma ^2})$ 的高斯分布，其均值与方差为

$\qquad {u_0} + \tau = \frac{{\displaystyle\sum\limits_{k = 1}^{k = 2K + 1} {{y_k}} }}{{2K + 1}}$

(10)

$\qquad {\sigma ^2} = \frac{{\displaystyle\sum\limits_{k = 1}^{k = 2K + 1} {{{[{u_0} + \tau - {y_k}]}^2}} }}{{2K + 1}}$

(11)

其中， ${u_0}$ 为杂波理论中值， $\tau$ 为 ${y_k}$ 与 ${u_0}$ 的差值， ${\sigma ^2}$ 表示 ${y_k}$ 的方差。

由式(9)可知，初始化后有序数组 ${{{Y}}_{r,d}}$ 中 ${y_k}$ 的数值都为 $x{'_{r,d}}$ ，方差值为0，并且当K值越大，采样滑窗所包含的杂波样本越多， $x{'_{r,d}}$ 越趋近 ${u_0}$ ， $\left| \tau \right|$ 数值越小。

由于本算法初始化参考OS-CFAR的采样策略，所以在迭代初期可获得与其接近的算法性能。

3.2 迭代优化

算法迭代通过将有序数组 ${{{Y}}_{r,d}}$ 中偏离理论中值 ${u_0}$ 的数据替换，保留接近 ${u_0}$ 的数据，从而获取更加准确的估计中值 ${y_{{\rm{mid}}}}$ 。

在数值替换中对极小值 ${y_1}$ 和极大值 ${y_{2K + 1}}$ 进行动态轮询替换，以雷达扫描次数t的奇偶性为标志，当雷达扫描次数t为奇数时，将待检测单元 ${x_{r,d}}$ 赋值给 ${{{Y}}_{r,d}}$ 中的 ${y_1}$ ；当t为偶数时，将 ${x_{r,d}}$ 赋值给 ${y_{2K + 1}}$ 。重新对 ${{{Y}}_{r,d}}$ 进行排序以保持有序性，以被替换极值点为起点，未被替换极值点为终点进行单向冒泡排序。K值越大，排序耗时越多。

由于初始化时 ${{{Y}}_{r,d}}$ 中的数据都相同，所以在算法迭代中需要经过2K+1轮迭代才能将 ${{{Y}}_{r,d}}$ 中的初始数据全部替换，之后 ${{{Y}}_{r,d}}$ 中的 ${y_{{\rm{mid}}}}$ 开始变化，初始数据被全部替换也标志着算法渡过2K+1轮迭代的初期阶段。

3.3 检测阈值计算

通过杂波估计中值 ${y_{{\rm{mid}}}}$ 逆推求解杂波模型参数，即令 $\displaystyle\int_0^{{y_{{\rm{mid}}}}} {f(x){\rm{d}} x} {\rm{ = }}\frac{1}{2}$ 可得杂波参数 $\mu = \dfrac{{{y_{{\rm{mid}}}}}}{{\ln 2}}$ 。

由于预设虚警概率已知，可通过式(4)求解标称因子 $\alpha$ 和检测阈值T为

$\qquad \alpha = - \frac{{\ln ({\rm{PFA}})}}{{\ln 2}}$

(12)

$\qquad T = \alpha {y_{{\rm{mid}}}}{\rm{ = }} - \frac{{{y_{{\rm{mid}}}}\ln ({\rm{PFA}})}}{{\ln 2}}$

(13)

使用本算法进行检测处理时，若杂波数据服从单参分布模型，可直接按原流程处理；若杂波数据服从多参分布模型，则需要先保留其中1个参数为未知状态，将其它参数由参数估计方法(最大似然估计、矩估计等)估算，之后按照原流程估算这个未知参数，计算检测阈值进行检测处理。

3.4 半参考窗值K的取值

半参考窗数值K需要综合考虑内存需求、运算速度、迭代速度、期望检测性能这4个因素进行设置。为便于分析K值对迭代速度和期望检测性能的影响，这里设计如下2个评价指标

${m_{t,K}} = \frac{1}{N}\sum\limits_{i = 1}^N {|{y_{t,K,i}} - {{\bar y}_{{\rm{mid}}}}|}$

(14)

$\sigma _{t,K}^2 = \frac{1}{N}\sum\limits_{i = 1}^N {{{(|{y_{t,K,i}} - {{\bar y}_{{\rm{mid}}}}| - {m_{t,K}})}^2}}$

(15)

其中， ${y_{t,K,i}}$ 为杂波估计中值， ${\bar y_{{\rm{mid}}}}$ 为杂波理论中值， ${m_{t,K}}$ 为 ${y_{t,K,i}}$ 与 ${\bar y_{{\rm{mid}}}}$ 之差绝对值的均值， $\sigma _{t,K}^2$ 为 ${y_{t,K,i}}$ 与 ${\bar y_{{\rm{mid}}}}$ 之差绝对值的方差， $t$ 为迭代次数， $K$ 为半参考窗值， $i$ 为蒙特卡罗实验次数， $N$ 为蒙特卡罗实验总次数。均值 ${m_{t,K}}$ 越小，表明 ${y_{t,K,i}}$ 越接近 ${\bar y_{{\rm{mid}}}}$ 、准确性越高；方差 $\sigma _{t,K}^2$ 越小，表明稳定性越高；准确性与稳定性越高计算出的检测阈值更接近理论值，可获得更优的期望检测性能。

利用杂波功率 $\lambda$ 为0 dBm的指数分布杂波数据进行2×10⁴组蒙特卡罗实验，每组实验进行2000次迭代、K取值范围为[1, 100]，对获得的 ${m_{t,K}}$ 和 $\sigma _{t,K}^2$ 取自然对数绘制等高线图谱如图2、图3所示。

图 2 均值等高线

下载: 全尺寸图片幻灯片

图 3 方差等高线

下载: 全尺寸图片幻灯片

由图2、图3可知，两图左上角出现三角形异常区域，这是由算法迭代初期杂波估计中值不变暂稳态导致。K取值在[10,20]范围内时，两图的评价指标处于等高线山谷区域，均值和方差较低(异常区域数据除外)，准确性和稳定性较高，所以应尽量将K取值在[10,20]范围内。

4. 仿真分析

将DOM-CFAR与最优检测器(opt), CA-CFAR, SO-CFAR, GO-CFAR, OS-CFAR和CM-CFAR进行4×10³组蒙特卡罗仿真实验，其中背景杂波功率 $\lambda$ 为20 dBm, DOM-CFAR的半参考窗值K为15, CA-CFAR, SO-CFAR, GO-CFAR和OS-CFAR的保护单元为4、滑窗长度为30, CM-CFAR的遗忘因子为1/128，预设虚警概率PFA为10^–4。

4.1 均匀环境下的性能分析

对均匀环境下单目标分别进行10次、200次、1000次迭代优化，其中信杂比取值范围为[0 dB, 30 dB]，其检测结果如图4、图5、图6所示。

图 4 均匀杂波下迭代10次

下载: 全尺寸图片幻灯片

图 5 均匀杂波下迭代200次

下载: 全尺寸图片幻灯片

图 6 均匀杂波下迭代1000次

下载: 全尺寸图片幻灯片

由图4可知，DOM迭代初期检测性能接近OS，优于SO，较CA, GO略差。

由图5可知，经过200次迭代DOM性能已优于参照对比算法。

由图6可知，经过1000次迭代DOM性能趋近于最优检测器性能。

4.2 多目标环境下的性能分析

在目标参考前窗设置相同信杂比的第2个目标形成拥有多目标的背景杂波环境进行10次、200次迭代优化，其中信杂比取值范围为[0 dB, 30 dB]，检测结果如图7、图8所示。

图 7 多目标下迭代10次

下载: 全尺寸图片幻灯片

图 8 多目标下迭代200次

下载: 全尺寸图片幻灯片

由图7可知，DOM迭代初期拥有与OS相近的多目标检测优势。

由图8可知，DOM检测性能不受邻近目标的影响，且经过200次迭代后，其检测性能已经优于参照对比算法。

4.3 杂波边缘环境下的性能分析

设置前500个距离单元杂波功率为20 dBm，后500个距离单元杂波功率为30 dBm，形成杂波边缘环境，分别进行10次、200次迭代检测，其结果如图9、图10。

图 9 杂波边缘下迭代10次

下载: 全尺寸图片幻灯片

图 10 杂波边缘下迭代200次

下载: 全尺寸图片幻灯片

由图9可知，DOM迭代初期的虚警概率在杂波边缘的低功率区相较于CA, GO和CM更接近预设值，在杂波边缘的高功率区相较于CA, SO更接近预设值。

由图10可知，DOM经过200次迭代，虚警概率在杂波边缘的低功率区更接近预设值，在杂波边缘的高功率区相较于CA, SO和OS更接近预设值。

4.4 记忆效应分析

DOM与CM都是利用杂波历史数据计算检测阈值的时序类恒虚警算法，算法处理时每帧数据(无论是否有目标能量信息)都参与迭代并影响后续检测，从而导致算法的记忆效应。

由于DOM的参考数据为估计中值 ${x_{{\rm{DOM}}}}(t)$ , CM的参考数据为估计均值 ${x_{{\rm{CM}}}}(t)$ ，为了统一评价指标进行如下处理

$\quad {s_{{\rm{DOM}}}}(t) = \frac{{|{x_{{\rm{DOM}}}}(t) - {P_{{\rm{dB}}}} \cdot {{\bar x}_{{\rm{DOM}}}}|}}{{{P_{{\rm{dB}}}} \cdot {{\bar x}_{{\rm{DOM}}}}}}$

(16)

$\quad {s_{{\rm{CM}}}}(t) = \frac{{|{x_{{\rm{CM}}}}(t) - {P_{{\rm{dB}}}} \cdot {{\bar x}_{{\rm{CM}}}}|}}{{{P_{{\rm{dB}}}} \cdot {{\bar x}_{{\rm{CM}}}}}}$

(17)

其中， ${P_{{\rm{dB}}}}$ 为杂波功率值， ${\bar x_{{\rm{DOM}}}}$ 为理论中值， ${\bar x_{{\rm{CM}}}}$ 为理论均值。评价指标 $s$ 的数值越小，参考值越接近理论值，估算的检测阈值更准确。

先在均匀杂波无目标的环境下进行2000次迭代使DOM和CM的性能稳定，从第2001次迭代开始，连续引入50次单目标回波，仿真目标能量参与迭代对性能的影响，其中目标能量信杂比分别为4 dB, 8 dB, 16 dB，结果如图11所示。

图 11 目标干扰对性能的影响

下载: 全尺寸图片幻灯片

由图11可知，当目标干扰出现并参与迭代，对DOM的影响有一定延时，当目标干扰消失时，DOM会以更快的速度重新优化，故其拥有更优的抗目标干扰能力。

4.5 算法复杂度与耗时

由于算法中每个待检测单元都需要一组长度为2K+1的有序数组用于存储其历史数据，因此空间复杂度(使用大O表示法)为O(MK)，其中M由待检测单元总个数决定。在迭代优化中动态极值替换的逻辑计算量很小几乎可以忽略，对运算速度起决定影响的是单向冒泡排序，因此其运算复杂度为O(K)。

以大O表示法为评价指标对各个算法的空间、运算复杂度以及在硬件为i7-7700HQ CPU条件下处理10帧实测RD谱数据的串行、并行计算耗时进行分析，结果如表1所示，其中N由滑窗大小决定、K由半参考窗值决定、M由待检测单元总个数决定。

表 1 算法复杂度及耗时比较

算法名称	空间复杂度	运算复杂度	串行耗时(s)	并行耗时(s)
CA	1	MN	51.337	2.776
SO	1	MN	73.588	3.260
GO	1	MN	75.176	3.289
OS	1	MNlgN	449.651	7.632
CM	M	M	0.087	0.073
DOM	MK	MK	88.490	5.560

下载: 导出CSV

| 显示表格

由表1及前文中各算法在均匀、多目标和杂波边缘环境下的性能对比可知，本算法与现有大部分算法相比，虽然存储空间需求较大、处理速度优势不突出，但在性能上得到了较大提高。同时根据已进行的雷达低空小目标探测实验可知，其算法复杂度对于目前的GPU硬件条件来说是可行的。

5. 实测数据验证

使用武汉大学电波传播实验室的外辐射源雷达系统^[21]在有风力发电机组环境下的实测数据进行验证分析，探测条件如图12所示。

图 12 雷达探测环境

下载: 全尺寸图片幻灯片

实测数据经多径杂波和多普勒扩展杂波抑制后的RD谱中，风机叶片转动的多普勒扩展杂波残余仍然较强，形成杂波边缘环境，如图13所示。

图 13 实测RD谱

下载: 全尺寸图片幻灯片

在实测RD谱上添加100个信杂比为15 dB的随机目标则形成拥有多目标的杂波边缘环境。将DOM-CFAR与CA-CFAR, SO-CFAR, GO-CFAR, OS-CFAR和CM-CFAR进行实测验证，其中DOM-CFAR的半参考窗值K为15, CA-CFAR, SO-CFAR, GO-CFAR和OS-CFAR的保护单元为4、滑窗长度为30, CM-CFAR的遗忘因子为1/128，迭代总次数为400，预设虚警概率为10^–6，结果如表2所示。

由表2数据可知，在迭代初期，DOM相较于参照对比算法拥有较优的检测性能，而虚警性能相对较差。但是DOM通过迭代优化后，检测概率会稳步提升，虚警概率会迅速下降，从而获得更优的检测性能与虚警性能。

表 2 实测的检测概率与虚警概率

CFAR	CA	SO	GO	OS	CM迭代次数				DOM迭代次数
CFAR	CA	SO	GO	OS	10	100	200	400	10	100	200	400
检测概率(%)	72.4	51.3	78.8	86.8	1.3	62.1	76.4	93.0	88.0	98.3	98.1	98.2
虚警概率(10^–5)	2.68	9.58	1.93	22.06	0.13	2.93	4.25	4.45	20.93	5.08	2.85	1.28

下载: 导出CSV

| 显示表格

6. 结束语

本文提出了一种基于动态有序矩阵的外辐射源恒虚警检测算法，利用历史杂波数据来动态替换有序矩阵中的极值，使得杂波估计中值向理论中值收敛，通过估计中值求解检测阈值。仿真与实测结果表明，该算法在迭代初期拥有较优的检测性能，通过多次迭代可以快速达到检测概率与虚警概率的最优性能，在均匀杂波、多目标和杂波边缘的复杂环境中检测性能稳定，并且相较于同属时序类算法的CM-CFAR，拥有更优的抗目标干扰能力。

图 1 DOM-CFAR算法实现流程图

下载: 全尺寸图片幻灯片

图 2 均值等高线

下载: 全尺寸图片幻灯片

图 3 方差等高线

下载: 全尺寸图片幻灯片

图 4 均匀杂波下迭代10次

下载: 全尺寸图片幻灯片

图 5 均匀杂波下迭代200次

下载: 全尺寸图片幻灯片

图 6 均匀杂波下迭代1000次

下载: 全尺寸图片幻灯片

图 7 多目标下迭代10次

下载: 全尺寸图片幻灯片

图 8 多目标下迭代200次

下载: 全尺寸图片幻灯片

图 9 杂波边缘下迭代10次

下载: 全尺寸图片幻灯片

图 10 杂波边缘下迭代200次

下载: 全尺寸图片幻灯片

图 11 目标干扰对性能的影响

下载: 全尺寸图片幻灯片

图 12 雷达探测环境

下载: 全尺寸图片幻灯片

图 13 实测RD谱

下载: 全尺寸图片幻灯片

表 1 算法复杂度及耗时比较

算法名称空间复杂度运算复杂度串行耗时(s) 并行耗时(s)

CA 1 MN 51.337 2.776
SO 1 MN 73.588 3.260
GO 1 MN 75.176 3.289
OS 1 MNlgN 449.651 7.632
CM M M 0.087 0.073
DOM MK MK 88.490 5.560

下载: 导出CSV

表 2 实测的检测概率与虚警概率

CFAR CA SO GO OS CM迭代次数 DOM迭代次数
10 100 200 400 10 100 200 400

检测概率(%) 72.4 51.3 78.8 86.8 1.3 62.1 76.4 93.0 88.0 98.3 98.1 98.2
虚警概率(10^–5) 2.68 9.58 1.93 22.06 0.13 2.93 4.25 4.45 20.93 5.08 2.85 1.28

下载: 导出CSV

参考文献(21)

FINN H M. Adaptive detection in clutter[C]. The 5th Symposium on Adaptive Processes, New Jersey, USA, 1966: 562–567. doi: 10.1109/SAP.1966.271149.

WANG Weijiang, WANG Runyi, JIANG Rongkun, et al. Modified reference window for two-dimensional CFAR in radar target detection[J]. The Journal of Engineering, 2019, 2019(21): 7924–7927. doi: 10.1049/joe.2019.0687

TRUNK G V. Range resolution of targets using automatic detectors[J]. IEEE Transactions on Aerospace and Electronic Systems, 1978, AES-14(5): 750–755. doi: 10.1109/TAES.1978.308625

HANSEN V G and SAWYERS J H. Detectability loss due to "greatest of" selection in a cell-averaging CFAR[J]. IEEE Transactions on Aerospace and Electronic Systems, 1980, 16(1): 115–118. doi: 10.1109/TAES.1980.308885

SMITH M E and VARSHNEY P K. Intelligent CFAR processor based on data variability[J]. IEEE Transactions on Aerospace and Electronic Systems, 2000, 36(3): 837–847. doi: 10.1109/7.869503

WANG Leiou, WANG Donghui, and HAO Chengpeng. Intelligent CFAR detector based on support vector machine[J]. IEEE Access, 2017, 5: 26965–26972. doi: 10.1109/ACCESS.2017.2774262

CARRETERO M V I, HARMANNY R I A, and TROMMEL R P. Smart-CFAR, A machine learning approach to floating level detection in radar[C]. The 16th European Radar Conference (EuRAD), Paris, France, 2019: 161–164.

ROHLING H. Radar CFAR thresholding in clutter and multiple target situations[J]. IEEE Transactions on Aerospace and Electronic Systems, 1983, AES-19(4): 608–621. doi: 10.1109/TAES.1983.309350

VILLAR S A, DE PAULA M, SOLARI F J, et al. A framework for acoustic segmentation using order statistic-constant false alarm rate in two dimensions from sidescan sonar data[J]. IEEE Journal of Oceanic Engineering, 2018, 43(3): 735–748. doi: 10.1109/JOE.2017.2721058

柳向, 李东生, 胡瑞. 基于有序统计类恒虚警检测的脉冲压缩雷达移频特征消隐多载波干扰研究[J]. 兵工学报, 2017, 38(11): 2134–2142. doi: 10.3969/j.issn.1000-1093.2017.11.008

LIU Xiang, LI Dongsheng, and HU Rui. Research on blanking shift-frequency-multi-carrier jamming against pulse-compression radar based on OS-CFAR[J]. Acta Armamentarii, 2017, 38(11): 2134–2142. doi: 10.3969/j.issn.1000-1093.2017.11.008

LIN C H, LIN Y C, BAI Yue, et al. DL-CFAR: A Novel CFAR target detection method based on deep learning[C]. The 90th IEEE Vehicular Technology Conference (VTC2019-Fall), Honolulu, USA, 2019: 1–6. doi: 10.1109/VTCFall.2019.8891420.

JIN Erwen, YAN Danqing, ZHANG Zhongjin, et al. FOD Detection on Airport Runway with Clutter Map CFAR Plane Technique[M]. LIANG Qilian, WANG Wei, MU Jiasong, et al. Communications, Signal Processing, and Systems. New York: Springer, 2012: 335–342. doi: 10.1007/978-1-4614-5803-6_34.

AKÇAPINAR K and BAYKUT S. CM-CFAR parameter learning based square-law detector for foreign object debris radar[C]. The 48th European Microwave Conference, Madrid, Spain, 2018: 421–424. doi: 10.23919/EuMC.2018.8541714.

TAO Ding, ANFINSEN S N, and BREKKE C. Robust CFAR detector based on truncated statistics in multiple-target situations[J]. IEEE Transactions on Geoscience and Remote Sensing, 2016, 54(1): 117–134. doi: 10.1109/TGRS.2015.2451311

WU Fengtao, WU Nan, and WU Maosong. A fast and slow time combined CFAR detection algorithm used in through-the-wall radar[C]. 2017 IEEE Electrical Design of Advanced Packaging and Systems Symposium, Haining, China, 2017: 1–3. doi: 10.1109/EDAPS.2017.8276955.

LAYEGHY S, ODABAEE M, KHLIF M S, et al. A time frequency approach to CFAR detection[C]. 2011 IEEE International Symposium on Signal Processing and Information Technology (ISSPIT), Bilbao, Spain, 2011: 230–234. doi: 10.1109/ISSPIT.2011.6151565.

代振, 王平波, 卫红凯. 非高斯背景下基于Sigmoid函数的信号检测[J]. 电子与信息学报, 2019, 41(12): 2945–2950. doi: 10.11999/JEIT190012

DAI Zhen, WANG Pingbo, and WEI Hongkai. Signal detection based on sigmoid function in Non-Gaussian noise[J]. Journal of Electronics &Information Technology, 2019, 41(12): 2945–2950. doi: 10.11999/JEIT190012

NARASIMHAN R S, RAMAKRISHNAN K R, and VENGADARAJAN A. Robust variability index CFAR for non-homogeneous background[J]. IET Radar, Sonar & Navigation, 2019, 13(10): 1775–1786. doi: 10.1049/iet-rsn.2018.5435

ZHANG Xin, ZHANG Renli, SHENG Weixing, et al. Intelligent CFAR detector for non-homogeneous weibull clutter environment based on skewness[C]. 2018 IEEE Radar Conference (RadarConf18), Oklahoma, USA, 2018: 322–326. doi: 10.1109/RADAR.2018.8378578.

赵文静, 刘畅, 刘文龙, 等. K分布海杂波背景下基于最大特征值的雷达信号检测算法[J]. 电子与信息学报, 2018, 40(9): 2235–2241. doi: 10.11999/JEIT171092

ZHAO Wenjing, LIU Chang, LIU Wenlong, et al. Maximum eigenvalue based radar signal detection method for K distribution sea clutter environment[J]. Journal of Electronics &Information Technology, 2018, 40(9): 2235–2241. doi: 10.11999/JEIT171092

万显荣. 基于低频段数字广播电视信号的外辐射源雷达发展现状与趋势[J]. 雷达学报, 2012, 1(2): 109–123. doi: 10.3724/SP.J.1300.2012.20027

WAN Xianrong. An overview on development of passive radar based on the low frequency band digital broadcasting and TV signals[J]. Journal of Radars, 2012, 1(2): 109–123. doi: 10.3724/SP.J.1300.2012.20027

施引文献

期刊类型引用(5)

1.	刘言，刘宁波，黄勇，王中训. 利用相位特征筛选参考单元的改进CFAR方法. 烟台大学学报(自然科学与工程版). 2023(03): 371-378 . 百度学术
2.	杜改丽，封治华. 一种改进现有机载雷达探测性能的多普勒处理新技术. 现代雷达. 2022(02): 75-83 . 百度学术
3.	邹俊杰，程丰，万显荣. 外源雷达空时联合恒虚警检测分析与实验. 雷达科学与技术. 2022(04): 415-420+428 . 百度学术
4.	周子铂，王彬彬，张朝伟，刘建卫，徐颖鑫，王志会. 基于迭代对消的外辐射源雷达目标检测方法. 雷达科学与技术. 2022(05): 555-564 . 百度学术
5.	杨威，崔恒荣. 77GHz近程雷达CFAR算法研究. 通信技术. 2021(09): 2125-2131 . 百度学术