顾新锋 简 涛 何 友
(海军航空工程学院信息融合技术研究所 烟台 264001)
距离扩展目标的双门限恒虚警检测器及性能分析
顾新锋 简 涛 何 友*
(海军航空工程学院信息融合技术研究所 烟台 264001)
针对普通积累检测器检测稀疏散射点目标存在“坍塌损失”以及依赖于散射点密度的二元积累(SDD-BI)检测器检测密集散射点目标存在信杂比损失的问题,该文充分利用目标的散射点密度及幅度信息,提出了一种距离扩展目标的双门限恒虚警检测器(DT-CFAR),在高斯杂波背景下推导了DT-CFAR虚警概率和检测概率与检测门限关系的解析表达式,并给出了最优第1门限的确定方法。仿真结果表明,DT-CFAR对不同散射点密度目标都具有较优的检测性能,并且对散射点密度估计失配具有鲁棒性。
目标检测;距离扩展目标;双门限;恒虚警率
早期的雷达由于距离量化精度较低,目标尺寸远小于雷达的距离量化单元宽度,利用雷达对目标进行检测时,通常把目标看成是一个点,采用点目标模型进行建模,并形成了一系列的点目标检测方法[1]。随着技术的进步,现代雷达具有较高的距离分辨率,其分辨单元尺寸已远小于目标尺寸,这种雷达称为高分辨率雷达,其目标回波分布在多个距离单元中,形成距离扩展目标[2]。这时,如果仍采用传统的点目标检测方法,一方面目标回波能量没有充分利用,高分辨雷达的优势得不到发挥,另一方面,由于目标的部分能量分散到参考单元中,反而降低了目标的检测概率[3]。提高雷达分辨率,并采用合适的方法对目标所占据的多个距离单元回波能量进行积累,可以显著提高目标的检测概率[4−11],其中,文献[4-8]利用所有距离单元回波能量进行积累,这种积累检测器对密集散射点目标具有较好的检测性能,而对于稀疏散射点目标,由于累积了不含强散射点距离单元能量,存在一定的“坍塌损失”。文献[9]采用了二元积累检测方法,这种检测器具有结构简单,易于实现的特点,由于没有考虑目标散射点的密度信息,对于稀疏散射点目标,存在严重的检测性能损失;针对这一问题,文献[10]利用散射点的展布特性,研究了最优二元检测方法,而文献[11]利用散射点的密度信息,提出了依赖于散射点密度的二元积累检测器(Scatterer Density Dependent-Binary Integrator, SDD-BI),这两种方法对不同散射点密度目标都能进行有效的检测,但需要已知目标散射点密度的先验信息,并且第2门限检测没有利用散射点的幅度信息,存在一定的信噪比损失。
为增强高分辨率雷达对不同散射点密度目标检测的鲁棒性,以及进一步提高目标的检测概率,本文充分利用目标散射点的密度和幅度信息提出了一种双门限恒虚警(Double Threshold-Constant False Alarm Rate, DT-CFAR) 检测器,其第1门限检测和第2门限检测都采用了恒虚警检测,并得到了检测器在高斯杂波背景中虚警概率与检测门限关系的解析表达式。文章结构安排如下:第2节是检测器的设计,第3节对检测器进行理论分析,第4节针对不同散射点密度的目标模型和实际目标1维距离像,分析检测器的检测性能,最后是总结。
假设雷达只发射1个脉冲。在接收端,对这个脉冲进行匹配滤波(或是脉冲压缩),这样目标的每个散射中心都能够被分解[4]。目标占据J个距离单元,每个散射点占据1个距离单元,各散射点回波及杂波回波相互独立,对各散射点回波采用复随机变量建模,则对距离扩展目标的检测可采用如下二元假设检验来表示:
式(1)中,H0表示目标不存在,H1表示目标存在,zr表示第r个距离单元的接收信号,sr表示第r个距离单元的目标回波信号,cr表示第r个距离单元的杂波回波信号。对每个距离单元做平方率检波并归一化得到dr=|zr|2/σ2,σ2为雷达I/Q通道中杂波的方差。当信号幅度信息完全未知时,可以采用广义似然比检验判别目标有无,其等价的检测器为积累检测器,即
式(2)中,d0为检测门限。对于密集散射点目标,积累检测器有最佳的检测性能。而对于实际目标,只少数的几个距离单元内存在强散射点回波,即实际目标为稀疏散射点目标,这时如仍用所有单元回波进行积累,由于积累了不含目标散射点的距离单元回波,使得平均信杂比降低,从而产生一定的信杂比损失,这种损失称之为“坍塌损失”[12]。
如果只利用含目标强散射点的距离单元回波进行积累,就可以有效地克服“坍塌损失”,但对于不同的目标,或目标相对于雷达的姿态角不同,其强散射点所在的距离单元位置也就不同,基于此,本文先设置第1门限,进行强散射点的检测,然后,将强散射点的回波减去第1门限后,进行积累,与第2门限比较判别目标有无,其第1门限和第2门限的设定都采用恒定的虚警概率,因此,称之为双门限恒虚警检测器,具体过程如下:
第1步 强散射点的CFAR检测:
对每个距离单元采用如下二元假设检验。
式(3)中,h0表示距离单元内不存在强散射点,h1表示距离单元内存在强散射点。对于每个距离单元,强散射点检测可以表示为
式(6)中,u(⋅)表示单位阶跃函数。图 1为检测器方框图。
特别地,当第1检测门限l0=0时,式(6)与式(2)等价,检测器退化为普通的积累检测器,即积累检测器是本文提出检测器在l0=0时的特例;当第2检测门限L0(K)=0时,这时J次检测中只要有一次超过第1门限,检测器就会做出H1判决,即认为目标存在,这时检测器就等价为J中取1检测器(第2门限设置为1的二元积累检测器)。
图1 双门限恒虚警检测器方框图
便于理论分析,假设目标的散射点密度为α(每个距离单元存在强散射点的概率),各强散射点的回波幅度为瑞利分布,且强度相等;杂波采用高斯杂波模型。经平方率检波和归一化后,第r个距离单元的回波dr是服从参数为μ的指数分布,其中,μ=2(1+λ),λ为距离单元内目标强散射点的信杂比。当距离单元内不存在目标强散射点时λ=0,μ=2。
当第2检测门限L0(K)=0时,DT-CFAR检测器退化为J中取1检测器。将L0(K)=0代入式(9)得P2fa(K)=1。再将P2fa=1代入式(12)可得第1门限
由式(19)~式(25)可知检测器的检测概率Pd与强散射点的信杂比λ、目标占据单元数J、目标散射点密度α、第1门限检测虚警概率、超过第1门限的回波个数K以及相应的第2检测门限L0(K)有关。由式(17),式(7)和可知L0(K)与检测器虚警概率Pfa、第1检测门限l0有关。目标信杂比定义为J个距离单元内所有强散射点的功率与单个距离单元内杂波功率的比,即
式(33)表示参数SCR,α,Pfa,J不变,改变参数l0,取使得g(S CR,l0,α,Pfa,J)达到最大值时的l0。例如,当SCR=14时,l0取4.767可使得Pd达到最大值。根据式(33)可以计算得到参数SCR,α,Pfa,J取不同值时对应的最优l0。表1 给出了Pfa= 10−4, SCR=14时目标所占据不同距离单元数及α取不同值时对应的最优l0取值。由表1可以看出,最优l0取值随着α增大而减小,当α>0.5时,l0取0时检测器有最佳的检测性能,这表明当目标所占据的单元中一般以上的距离单元存在强散射点时,采用积累检测器就能得到最佳的检测性能。进一步比较分析图2和表1,可以得出目标散射点密度α是影响最佳l0取值的关键因素,其次是目标占据单元数,而信杂比对最佳l0取值的影响较弱。
表1 Pfa= 10-4, SCR=14时最优l0取值
图 3~图 5给出了Pfa=410−,J=40,α分别为0.05(极稀疏散射点目标),0.2(稀疏散射点目标)和0.5(中等密集散射点目标)时,本文检测方法(DTCFAR),积累检测器(Integrator)和SDD-BI的检测性能比较曲线。DT-CFAR检测器的第1门限查表1得到, SDD-BI检测器的第2门限分别取1, 2和5。从图 3可以看出,对于极稀疏散射点目标,DTCFAR检测器略优于SDD-BI检测器而明显优于积累检测器,这是因为DT-CFAR和SDD-BI利用了散射点密度的先验信息,而积累检测器由于积累了不含目标强散射点的距离单元回波,存在一定的“坍塌损失”。图4表明,DT-CFAR充分利用了目标散射点密度和回波幅度信息,其性能明显优于SDD-BI检测器和积累检测器。图5则表明,对于中等密集散射点目标,DT-CFAR检测器等价为积累检测器而明显优于SDD-BI检测器,这是因为当α=0.5时DT-CFAR检测器的最优第 1门限l0=0,这时DT-CFAR检测器等价为积累检测器,而 SDD-BI检测器由于没有利用幅度信息,存在一定的信杂比损失。
图3~图5的分析结果表明,当目标散射点密度已知时DT-CFAR检测器具有最佳的检测性能,而对于实际的目标,即使已知目标1维距离像(图6为某实测目标的归一化1维距离像),也很难得到目标散射点密度α的确切值,这时需要先对目标散射点密度进行估计,然后根据估计值确定检测器的最佳参数。
图7给出了DT-CFAR检测器、SDD-BI检测器的检测概率与目标散射点密度估计值的关系曲线并与积累检测器进行比较,其中,目标散射点密度估计值取1/64到0.5(间隔为1/64),信杂比SCR设定为15 dB。从图7可以看出,当目标散射点密度估计值取3/32时DT-CFAR检测器具有最大的平均检测概率,而当散射点密度估计失配时,当α<0.3时DT-CFAR检测器仍优于SDD-BI检测器和积累检测器,而当α≥0.3时DT-CFAR检测器等价为积累检测器,这是因为当α≥0.3, DT-CFAR检测器的最优第1门限l0=0,退化为积累检测器。分析图7可知,在散射点密度估计失配时SDD-BI检测器检测性能明显下降,而DT-CFAR检测器仍保持较好的检测性能,这说明DT-CFAR检测器对目标散射点密度估计失配具有较好的鲁棒性。
图6 某目标的归一化1维距离像
图7 Pd与散射点密度估计值的关系曲线
SDD-BI检测器对于稀疏散射点目标具有较优的检测性能,而对于密集散射点目标存在一定的量化损失;积累检测器对于密集散射点目标具有较好的检测性能,而对于稀疏散射点目标存在严重的“坍塌损失”。综合这两种检测器的优点,本文充分利用目标散射点的密度和幅度信息,提出了 TD-CFAR检测器,推导了其虚警概率和检测概率的解析表达式,并给出了最优第1门限的确定方法,为工程应用提供了必要的理论指导。通过比较TD-CFAR检测器、积累检测器和SDD-BI检测器对不同散射点密度目标的检测性能得出:对于不同散射点密集程度目标,与SDD-BI检测器和积累检测器相比,TDCFAR检测器都有最佳的检测性能。最后,对某实测目标1维距离像的检测结果表明,TD-CFAR检测器在散射点密度估计失配时具有较好的鲁棒性。
[1] 何友, 关键, 孟祥伟, 等. 雷达目标检测与恒虚警处理[M]. 第2版, 北京: 清华大学出版社, 2011: 1-4, 115.
He You, Guan Jian, Meng Xiang-wei,et al.. Radar Target Detection and CFAR Processing[M]. 2nd Ed, Beijing:Tsinghua University Press, 2011: 1-4, 115.
[2] Wehner D R. High-Resolution Radar[M]. 2nd Ed . Boston:MA, Artech House, 1995: 220-221.
[3] Gerlach K and Steiner M J. Adaptive detection of range distributed targets [J].IEEE Transactions on Signal Processing, 1999, 47(7): 1844-1851.
[4] Hughes P K. A high-resolution radar detection strategy [J].IEEE Transactions on Aerospace Electron System, 1983,19(5): 663-667.
[5] Bandiera F, Orlando D, and Ricci G. CFAR detection strategies for distributed targets under conic constraints [J].IEEE Transactions on Signal Processing, 2009, 57(9):3305-3316.
[6] Dai F Z, Liu H W, and Wu S J. Generalized adaptive subspace detector for range-Doppler spread target with high resolution radar[J].Science China Information Sciences, 2011,54(1): 172-181.
[7] Jian T, He Y, Su F,et al.. High resolution radar target adaptive detector and performance assessment [J].Journal of Systems Engineering and Electronics, 2011, 22(2): 212-218.
[8] Guan J and Zhang X L. Subspace detection for range and Doppler distributed targets with Rao and Wald tests [J].Signal Processing, 2011, 91(1): 51-60.
[9] 孟祥伟, 曲东才, 何友. 高斯背景下距离扩展目标的恒虚警检测[J]. 系统工程与电子技术, 2005, 27(6): 1012-1015.
Meng Xiang-wei, Qu Dong-cai, and He You. CFAR detection for range-extended target in Gaussian background[J].Systems Engineering and Electronics, 2005, 27(6): 1012-1015.
[10] 陈远征, 周剑雄, 付强. 扩展目标的最优二进制检测[J]. 系统工程与电子技术, 2011, 33(1): 26-29.
Chen Yuan-zheng, Zhou Jian-xiong, and Fu Qiang. Optimal binary detection for distributed targets [J].Systems Engineering and Electronics, 2011, 33(1): 26-29.
[11] He Y, Gu X F, Jian T,et al.. A m out of n detector based on scattering density [C]. IET International Radar Conference 2009 (IET 2009), Guilin, April 20-22, 2009: 1-4.
[12] Goldman H. Performance of excision CFAR detector in the presence of interferes [J].IEE Proceedings-F Radar and Signal Processing, 1990, 137(3): 163-171.
Double Threshold CFAR Detector of Range-spread Target and Its Performance Analysis
Gu Xin-feng Jian Tao He You
(Research Institute of Information Fusion,Naval Aeronautical and Astronautical University,Yantai264001,China)
For the problem of normal integrator exists serious “collapsing loss” when detecting spare scatterer targets, and Scatterer Density Depend-Binary Integrator (SDD-BI) exists signal to clutter ratio loss when detecting dense scatterer targets, a Double Threshold-Constant False Alarm Ratio (DT-CFAR) detector of range-spread target is proposed by using the density and amplitude of a target scatterer. The formula relating false alarm probability and detection probability to detection threshold is deduced in the Gaussian clutter background,and the optimal first threshold for different target parameters is given. Experiment results show that DT-CFAR has better detection performance for diverse scatterer density targets and has robustness when the scatter density of a target is estimated mismatching.
Target detection; Range-spread target; Double threshold; Constant False Alarm Rate (CFAR)
TN957.51
A文章编号:1009-5896(2012)06-1318-06
10.3724/SP.J.1146.2011.01094
2011-10-21收到,2011-12-23改回
国家自然科学基金(61032001,61102166)资助课题
*通信作者:何友 ncif2009@126.com
顾新锋: 男,1983年生,博士生,研究方向为高分辨雷达目标检测.
简 涛: 男,1980年生,博士生,研究方向为雷达信号检测与信号处理.
何 友: 男,1956年生,教授,研究方向为雷达自适应检测方法、多传感器信息融合、模式识别等.