法律状态公告日
法律状态信息
法律状态
2013-11-20
专利权的转移 IPC(主分类):G01S13/04 变更前: 变更后: 登记生效日:20131101 申请日:20080730
专利申请权、专利权的转移
2010-12-29
授权
授权
2009-02-18
实质审查的生效
实质审查的生效
2008-12-24
公开
公开
技术领域
本发明属于雷达目标检测技术领域,涉及一种雷达目标恒虚警检测方法。
背景技术
目标检测是雷达的首要任务,恒虚警检测是一种在雷达自动检测系统中给检测策略提供检测阈值并且使杂波和干扰对系统虚警概率影响最小的信号处理方法。经过30余年的发展,恒虚警检测已逐渐成为现代雷达的一项标准技术,并在各种类型的雷达中获得广泛应用。
经典的恒虚警检测方法是先利用待检单元前后滑窗中的参考单元产生一个自适应门限,然后,比较待检单元信号强度与自适应门限的大小,如果待检单元信号强度超过门限,则判断该单元中有目标存在。这种经典的恒虚警检测方法的框图如图1所示。常用的单元半均类恒虚警(CA-CFAR)检测器以及顺序统计类恒虚警(OS-CFAR)检测器都采用这种方案。CA-CFAR检测器根据参考单元信号强度的均值生成自适应门限,在均匀高斯背景下,单元信号强度的均值是杂波强度的最优估计,因此,它具有最佳的检测性能。然而,在非高斯背景下,单元信号强度的均值不再是杂波强度的最优估计,其检测性能有所下降;在多目标干扰情况下,干扰信号的存在会抬高杂波强度的估计,导致检测性能迅速下降。
OS-CFAR检测器先将参考单元的信号强度排序,选择排序后的第k个最小值作为杂波强度水平估计。该方法屏蔽了干扰信号对杂波强度估计的影响,因此,在多目标干扰环境下,能保持较好的检测性能。然而,在均匀环境下,由于利用排序后的第k个最小值估计杂波水平不是最优的,其检测性能会有损失。因此,就需要一种适应性更好的检测方法,能够在多目标干扰环境下和非高斯环境下,都保持较优的检测性能。
发明内容
本发明提供一种基于拟合优度检验的雷达目标恒虚警检测方法。本发明利用的是目标回波与背景杂波分布特性的差异性,和传统的基于自适应门限的检测方法相比,它受背景分布特性和干扰目标的影响很小,对非高斯环境和多目标干扰环境,具有很强的适应性。另外,该方法还通过对参考单元样本做排序截尾处理,能够进一步抑制干扰目标的影响。
本发明利用了目标回波与背景杂波不同的统计分布特性,通过检验待检区的回波样本是否服从背景分布,从而判断是否有目标存在。也就是说,如果待检区的回波样本服从背景分布,则有理由相信,待检区回波源于背景杂波,从而判断没有目标存在;否则,将判断有目标存在。假如背景分布函数为F(x),雷达接收机通过脉冲积累,得到N个来自于某待检单元的回波样本Z={z1,z2,...,zN},目标检测可以转化为以下拟合优度检验问题:
为了解决(1)式所定义的基于拟合优度检验的目标检测问题,并保持目标检测的恒虚警性,本发明技术方案如下:
一种基于拟合优度检验的雷达目标恒虚警检测方法,如图2所示,包括以下步骤:
步骤1通过雷达接收系统接收N个脉冲,形成快慢时间域上的N行R列的数据矩阵,其中R表示雷达探测空间总的距离单元数。
步骤2将步骤1所得的数据矩阵进行对数放大,将威布尔(Weibull)型的背景分布转化为极大值分布型的位置-尺度(Location-Scale)类背景分布,得到对数放大后的N行R列的数据矩阵。引入对数放大器的目的是利用位置-尺度(Location-Scale)型背景分布的一些特性以保证本发明的恒虚警性。
步骤3根据步骤2所得的对数放大后的N行R列的数据矩阵,对于任一待检距离单元,都对应有N个待检样本Z={z1,z2,...,zN};选择前后相邻的M个距离单元作为参考单元,得到L=M×N个背景样本。
步骤4利用步骤3所得的背景样本估计背景分布的位置参数和尺度参数,具体包括以下步骤:
步骤4-1对L=M×N个背景样本从小到大(或从大到小)排序并删除后面(或前面)的r个样本以屏蔽干扰目标影响,得到有序背景样本序列Y=(y(1),...,y(L-r))T。
步骤4-2采用最优线性无偏估计器(BLUE),按下式计算得到背景分布位置参数估计值和尺度参数估计值
其中,D=(IL-r u0)是(L-r)×2的辅助矩阵,而IL-r是一个(L-r)维的单位向量;u0和C0分别是(L-r)维标准有序向量Y0的均值和协方差矩阵。
步骤5对步骤3所得的任一待检距离单元的N个待检样本Z={z1,z2,...,zN}按下式进行归一化处理,得到归一化后的待检样本Z′={z′1,z′2,...,z′N};
由于最优线性无偏估计和是等变化估计,归一化后的待检样本Z′={z′1,z′2,...,z′N}的统计特性与背景真实参数无关。
步骤5对归一化后的待检距离单元样本Z′={z′1,z′2,...,z′N}采用Anderson-Darling(AD)拟合优度检验:
其中,F(·)表示极大值分布函数。
如果检验统计量A2大于所设门限η,H1假设被接受,即样本Z′={z′1,z′2,...,z′N}不服从极大值分布型背景分布,从而判断出该待检距离单元有目标存在;如果检验统计量A2小于所设门限η,H0假设被接受,即样本Z′={z′1,z′2,...,z′N}服从极大值分布型背景分布,从而判断出该待检距离单元没有目标存在。
需要说明的是:
1、因归一化后的待检样本Z′={z′1,z′2,...,z′N}的统计特性与背景真实参数无关,步骤5中的检验统计量A2的统计特性也与背景真实参数无关,因此,基于(4)式判决时的虚警概率(PFA)也与真实的背景参数θS、θL无关,这就证明了本发明的恒虚警性。
2、本发明技术方案中,脉冲数
3、步骤5中所述门限η是在预先设定的最终虚警概率PFA(通常PFA=10-5或10-6)下,通过蒙特卡罗(Monte Carlo)仿真求得的。即在H0假设下,做Np(Np=102/PFA)次蒙特卡罗(Monte Carol)试验,然后搜索门限η,使得
Num(A2>η|H0)/Np=PFA (5)
其中,Num(A2>η|H0)为H0假设下,检验统计量A2超过门限的次数。
4、本发明在对应的预设虚警概率PFA下求得门限η后,可按下式计算检测概率PD
PD=Num(A2>η|H1)/Np (6)
本发明的有益效果是:
本发明将目标检测问题转化为拟合优度检验问题,通过检验待检单元样本与背景分布特性的差异性来检测目标,它受非高斯背景和干扰目标的影响较小;同时,对背景单元样本从小到大(或从大到小)排序并截尾(或去头)处理方法进一步抑制干扰对背景分布参数估计的影响,因此,本发现提供的方法具有更强的适应性。
附图说明
图1为经典的基于滑窗自适应门限的恒虚警检测方案示意图。
图2为本发明的流程图。
图3针对本发明具体实施方式给出的仿真例子,比较了Weibull背景下,虚警率PFA=10-5时,基于拟合优度检验的恒虚警检测(GoF CFAR)和经典的基于滑窗自适应门限的最优线性无偏估计恒虚警检测(BLUE CFAR)的性能。
图4针对本发明具体实施方式给出的仿真例子,比较了Weibull背景下,虚警率PFA=10-5时,存在两个干扰目标时,GoF CFAR和BLUE CFAR在几种不同截尾(或去头)深度下的性能。
具体实施方式
下面结合一个仿真例子给出本发明的具体实施方式,并比较其与经典检测方法的性能。
设某雷达波束宽度θB=2度,脉冲重复频率fp=940Hz,天线转速wr=46度/秒根据
可知驻留时间内积累的脉冲数目N=40。选择待测距离单元前后M=8个距离单元做为参考单元。背景为较为尖锐的山地杂波,服从Weibull分布,形状参数c=0.63。经对数放大器后,Weibull分布被转化为位置-尺度型的Gumbel分布。设目标辐度波动服从瑞利分布。
图4比较了本发明提出的基于拟合优度检验的恒虚警检测方法(GoF CFAR)和经典的基于滑窗自适应门限的恒虚警检测方法的性能,为了让比较更加客观,在此比较中,选择最优线性无偏估计恒虚警检测方法(BLUE CFAR)做为精典检测方法的代表,其采用的自适应门限为:
其中,和分别为背景分布的形状参数和尺度参数的最优线性无偏估计。
由图3可知,在较为尖锐的Weibull杂波环境下,在虚警概率PFA=10-5,检测概率PD=0.8以上时,GoF CFAR比BLUE CFAR信杂比改善3dB。
图4比较了存在2个干扰目标影响下,GoF CFAR和BLUE CFAR在r=2N、r=4N、r=0这几种不同截尾(或去头)深度时的检测性能。由图4可以发现,在恰当的截尾(或去头)深度下(r=2N),GoF CFAR和BLUE CFAR的检测性能和均匀环境下的检测性能相比,只有微小下降,但是,可以成功抑制两个干扰目标。当截尾深度变深时(r=4N),一半的参考单元被删除,此时尽管可以抑制4个干扰目标,但GoF CFAR和BLUE CFAR的检测性能都有相当明显下降。在没有截尾筛选的情况下(r=0),用污染的参考单元估计形状参数和尺度参数时,两种检测方法都不能再保证是虚警性的,检测性能下降很快。
机译: 基于空间的合成孔径雷达目标检测方法
机译: 一种基于动态光强度的车辆检测方法,夜间车辆检测方法及系统
机译: 一种基于探针的PCR检测方法来测量循环的脱甲基化的贝塔细胞衍生DNA的水平作为糖尿病中贝塔细胞损失的一种测量方法