2算法介绍与模型描述 本文假设接收噪声是高斯白噪声,平方律检波输出服从 指数分布,目标类型是SwedingII型,各参考单元统计独立同 分布(IDD)。有目标时检测单元的概率密度函数(PDF)为 1 V(i)=(n—r—i+1)( ( )一 (f一1)),(i=1,…,n)(5) 根据式(5)可得 J , 南数(MGF)为 ,( ,… )(6) (7) 将P用 n}来表示,经不复杂的推导可得r的矩产生函 ,( ) 一 , >0 ( ) 无目标时各参考单元所服从的概率密度函数(PDF)和 概率分布函数(CDF)为: _厂( )= e~: ̄/xi, ( ):E[e-Up]=[1+ >0 (2) 可以看出,UMVE方法与基于 —r个独立样本的 方 法性能相同,相当于基于 —r个独立样本的最大似然估计。 F(x)=1一e-x/ix, >0 (3) 对式(7)取反拉氏变换,可得P的概率密度函数为 其中 是噪声功率水平,A是信号与噪声的平均功率之 , …, 一 _=£ 比。 ,(p)= “ ( /.t)… 南1 n—r, (8) 2.1 无偏最小方差估计算法 2.2 MUMCA—CFAR检测器模型 为了解决CM方法没有充分利用剩余样本的问题,文献 MUMCA—CFAR检测器的原理图如图1所示。图中的停 [5]提出了无偏最小方差估计算法,它是将参考单元采样值 止移位控制器 和参考单元移位寄存器一起实现自动筛选 }按幅值排序,得到排序后的序列{ ( )≤…≤ ( )≤… 功能,即将检测为目标的样本剔除掉,使之不进入后沿参考 ≤ ㈨},删除最大的r个样值,对剩余样值采用加权平均方 滑窗。D为检测单元,M和Ⅳ分别为前、后沿参考滑窗长度。x 法得到背景功率水平估计,所采用的统计量为 和Y分别是对前后沿滑窗数据采用UMVE方法和cA方法得 1 n一1 p (∑%)+(r+1) ) (4) 到的两个局部估计,Z:X+l,是总的背景噪声估计, 是根 据相应的恒虚警算法和设定的虚警概率求出的门限系数。自 显然,经排序后的样值之间不再满足统计独立,为了求 H1 出P的肋F,引入下面的线性变换,产生统计独立的随机变 适应判决准则为D享TZ,其中 ,表示目标存在,Ⅳ0表示目 量 标不存。检测器的虚警概率为 图1 MUMCA—CFAR检测器原理框图 :P( {no) 』) ( ) 一 dxdz MUMCA—CFAR检测器使用如下两个局部估计: l e- (z)dz=E ]I :T/p= (u)1 : (9) ( ,+(r+ ) c ),y 其中 ( )是背景功率水平估计z的概率密度函数, (11) 其中X㈤是将前沿参考滑窗采样值X 按幅值由小到大 (M)是Z的磐产生函数。当 T 时,式(9)变转化 排序后得到的序列。显然,X和Y是独立的,由式(7),采用 成均匀背景下检测概率的表达式 UMVE方法得到的局部估计X的矩产生函数为 =P(H l H。)=M (u)l :— (10) (u)=【l+M-r] (12) 可见求得背景功率估计z的矩产生函数 (1L)是分析 由CA方法得到的局部估计y的矩产生函数为 检测器性能的关键。 Mr(u (13) 3 均匀背景下的MUMCA—CFAR检测器 由于两个独立随机变量和的矩母函数等于各随机变量 ・--——322---—— 维普资讯 http://www.cqvip.com 矩母函数的乘积 ,所以z= +Y的矩产生函数为 (M)=Mx( )×M (“) =表2 MUMCA和OS在多目标环境下的CFAR损失(dB) ,,J,/R 0,0 1,0 2,0 3,0 4,0 5,0 6,0 [ + r × ㈤ 吾=[ + r × (14) MUMCA 1.2923 1.3445 1.6947 1.8521 2.3101 2.6800 3.7440 0S 1.397l 1.6232 2.1072 2.3406 3.0774 3.7026 4.6983 均匀背景下,MUMCA—CFAR检测器的虚警概率和检测 概率为 一= 表中 MUMCA:m=n=16,r=6 OS:R=32。k=26 从表2可以看出,在多目标干扰情况下MUMCA的 (15) CFAR损失比OS要小,而且随着干扰目标数的增多,MUMCA 的CFAR损失增长缓慢,相对于OS性能优势越明显,如IL= 6,IR=0,时,它比OS改善了近ldB。 一=Mz(u)J : ! (1+A)(M—r) 5结论 本文基于无偏最小方差估计(UMVE)和单元平均(CA) 对任意设定的 值,利用数值方法解式(15),可以得到 门限系数 。 己,, 提出的一种新的恒虚警检测器(MUMCA),它更为有效地利 用了参考单元所提供的信息,具有更少的恒虚警损失。在均 匀背景和多目标环境下,MOSUM的性能比OS和MOSCA获 得了改善。同时MUMCA也是一种更容易实现的CFAR方 在均匀背景中,表1给出了参考单元长度M=N=16, 检测概率P =0.5,不同虚警概率时MUMCA,MOSCA和0S3 种检测器的CFAR损失。对于MUMCA和MOSCA—CFAR,取 r=6,它们具有容纳6个干扰目标的能力,对于OS—CFAR, 为了与MOSUM—CFAR具有相同的干扰目标容许限度,取k :案,它的样本排序时间只有Os和CM的四分之一。 参考文献: 『1] H M Finn and R S Johnson.Adaptive detection mode with thresh— old control as a function of spafiMly sampled clutter——level esti・ -26。从表1的结果可以看出,均匀背景下,MUMCA的性能 比MOSCA和OS获得了改善。 表1 3种检测器在均匀背景中的CFAR损失(dB) mates[J].RCA Review,1968,29(9):414—464. [2] J T Richard and G M Dillard.Adaptive detection algorithms for multiple—Target situations(J].IEEE Trans.On AES,1977,13 (4):338—343. [3]何友,关键.一种基于排序和平均的新恒虚警检测器[J].现 代雷达,1995,17(4):32—36. [4] 马健,等.一种新的基于删除平均和单元平均的恒虚警率检测 器[J].空军工程大学学报,2003,4(1):26—29. 表中 MUMCA:m=n=16,r=6 MOSCA:m=n=16,r:6 OS: =32.k:26 [5] D T Nagle,J Saniie.Performance analysis of linearly combined or- der statisitc CFAR detectors[J].IEEE Trans.On AES,1995,32 (2):522—533. [6]He You.Performance of some generliased modiifed order statistics 4 多目标环境下MUMCA—CFAR检测器的性能 在多目标环境下,为了分析方便,本文仅讨论强干扰目 标的影响。所谓强干扰是指干扰目标回波与背景杂波功率比 INR无限大,从而使干扰目标回波样本总是占据最高的有序 采样,一般来说,这是一种最差的情况,对于有限的INR,检 CFAR detectors with automatic censoring technique in multiple tar- get situations[C].IEE Pmc.一F,Radar and Signal Processing, I994,141(4):205—212. [7]何友等.雷达自动检测与恒虚警处理(第二章)[M].北京:清 华大学出版社.1999—5. 测损失将会变小。在采用自动筛选技术后,强干扰目标进入 后沿滑窗的概率很小,所以,仅考虑干扰目标位于前沿滑窗 的情况(/R=0)。当进入前沿滑窗的干扰目标数/L<r时,可 [作者简介] 曲 超(1977一),男(汉族),黑龙江人,博士生,主 要研究方向:信号检测与估计; 以用M一儿,r一儿分别代替式(16)中的M和r来获得 MUMCA在多目标环境中的检测性能,当儿=r时,前沿参考 郝程鹏(1975一),男(汉族),黑龙江宝清人,博士, 副研究员,主要研究方向:信号检测与估计; 滑窗剩余样本统计独立,对剩余样本采用CA方法即可实现 无偏最小方差估计。表2给出了多目标环境下MUMCA和OS 的CFAR检测损失,其中 =10硼。,P :0.5。 与识别。 杨树元(1942一),男(汉族),河北人,主任研究员、 博士生导师,主要研究方向:数字信号处理、阵列信号处理,图像处理 323.--—.—
因篇幅问题不能全部显示,请点此查看更多更全内容