广义逆高斯纹理海杂波背景下的自适应失配检测器

范一飞; 陈铎; 粟嘉; 郭子薰; 陶明亮; 王伶

doi:10.11999/JEIT231440

广义逆高斯纹理海杂波背景下的自适应失配检测器

doi: 10.11999/JEIT231440

范一飞^1, ,,
陈铎^{1, 2},
粟嘉¹,
郭子薰¹,
陶明亮¹,
王伶¹

1.
西北工业大学电子信息学院西安 710072
2.
中国航空工业集团公司洛阳电光设备研究所洛阳 471000

基金项目: 国家自然科学基金(62171379,62301435)，中国博士后科学基金(2023M732870)，博士后创新人才支持计划(BX20230497)，上海航天科技创新基金(SAST2023-044)

详细信息

作者简介:
范一飞：男，副研究员，博士，研究方向为雷达信号处理、海面微弱目标检测

陈铎：男，硕士，研究方向为海杂波背景下的微弱目标检测

粟嘉：男，副教授，博士，研究方向为雷达信号处理、雷达干扰抑制

郭子薰：女，博士后，博士，研究方向为雷达信号处理、海面微弱目标检测

陶明亮：男，副教授，博士，研究方向为雷达信号处理、雷达干扰抑制

王伶：男，教授，博士，研究方向为阵列信号处理、抗干扰处理

通讯作者:
范一飞　fanyifei888@163.com

中图分类号: TN95
计量
- 文章访问数: 283
- HTML全文浏览量: 166
- PDF下载量: 34
- 被引次数: 0
出版历程
- 收稿日期: 2024-01-02
- 修回日期: 2024-04-24
- 网络出版日期: 2024-05-15
- 刊出日期: 2024-09-26

Adaptive Detectors for Mismatched Signal under Sea Clutter Background with Generalized Inverse Gaussian Texture

FAN Yifei^{1
, ,},
CHEN Duo^{1, 2},
SU Jia¹,
GUO Zixun¹,
TAO Mingliang¹,
WANG Ling¹

1.
School of Electronics and Information, Northwestern Polytechnical University, Xi’an 710072, China
2.
Luoyang Institute of Electro-Optical Equipment, Aviation Industry of China, Luoyang 471000, China

Funds: The National Natural Science Foundation of China (62171379, 62301435), China Postdoctoral Science Foundation (2023M732870), The Postdoctoral Innovation Talents Support Program (BX20230497), ShangHai Aerospace Science and Technology Innovation Fund (SAST2023-044)

摘要

摘要: 针对雷达对海探测过程中理论导向矢量与实际导向矢量之间不匹配导致的虚警概率升高的问题，该文在复合高斯模型(CGM)下设计自适应失配检测器。为了抑制失配信号，在零假设中引入与理论导向矢量正交的虚拟信号，从而给出存在失配信号的目标检测模型。将CGM的纹理分量建模为广义逆高斯分布，分别基于两步广义似然比(GLRT)和最大后验GLRT(MAP GLRT)准则发展类似于自适应波束形成器正交抑制检测(ABORT)的自适应失配检测器，并通过理论证明所提失配检测器对散斑协方差矩阵和目标多普勒导向矢量具有恒虚警(CFAR)特性。仿真和实测数据实验结果表明，所提失配检测器在导向矢量匹配情况下的检测性能和失配情况下的抗失配性能之间具有良好的折衷。
- 海杂波 /
- 广义逆高斯纹理复合高斯模型 /
- 失配检测 /
- 自适应波束形成器正交抑制检测
Abstract: Considering mismatched problem between theoretical steering vector and actual steering vector causes false-alarm-rate increase in the process of maritime radar detection, the adaptive mismatched detectors are studied under Compound Gaussian Model (CGM). In order to reject mismatched signal, the fictitious signal orthogonal to theoretical steering vector is introduced in the null hypothesis, and a target detection with mismatched signal is given. The texture component of CGM is represented by generalized inverse distribution, and the Adaptive Beamformer Orthogonal Rejection Test (ABORT) is developed based on two-step Generalized Likelihood Ratio Test (GLRT) and Maximum A Posteriori GLRT (MAP GLRT) criterions respectively. Both the proposed detectors are testified to have Constant False AlaRm (CFAR) characteristics for speckle covariance matrix and target doppler steering vector. Experimental results based on simulated and real measured sea clutter data indicate that the proposed mismatched detectors show preferable target detection performance under the matched steering vector condition and anti-mismatch capability under the mismatched steering vector condition.

HTML全文

1. 引言

随着雷达分辨率的不断提升，海杂波呈现出明显的重拖尾和非高斯特征，不再满足中心极限定理的假设，导致高斯杂波背景下设计的检测器虚警和漏检概率升高^[1]。复合高斯模型(Compound Gaussian Model, CGM)将高分辨海杂波建模为慢变纹理分量调制快变散斑分量的过程，能够较好地拟合高分辨海杂波的非高斯特性。利用不同分布的随机变量建模CGM的纹理分量即可得到不同的杂波模型。当纹理分量分别服从伽马分布、逆伽马分布和逆高斯(Inverse Gaussian, IG)分布时，海杂波幅度序列的概率密度函数(Probability Density Function, PDF)分别服从K分布^[2]，Pareto分布^[3]和逆高斯纹理复合高斯模型(Compound Gaussian model with Inverse Gaussian texture, CG-IG)^[4]。然而，受海况、擦地角、雷达极化方式和载频等因素的影响，上述3种杂波模型的拟合适用范围仍然存在限制，且存在模型选择问题。针对此问题，文献[5]提出使用广义逆高斯(Generalized Inverse Gaussian, GIG)分布建模CGM的纹理分量，得到具有高泛化能力和高拟合优度的广义逆高斯纹理复合高斯模型(Compound-Gaussian model with Generalized Inverse Gaussian texture, CG-GIG)。由于伽马分布、逆伽马分布和IG分布分别是GIG在不同参数下的特例，所以在GIG纹理背景下设计的相干检测器能够最大程度避免因拟合误差造成的性能损失。

为了解决非高斯、时变杂波背景下的目标检测问题，文献[6]首次在CGM下提出纹理分量未知的归一化自适应匹配滤波(Adaptive Normalized Matched Filter, ANMF)检测器。ANMF检测器虽然适用于所有纹理结构的复合高斯杂波，但其对任意一种纹理结构的杂波均不具有最优检测性能。针对伽马纹理结构下的K分布杂波，文献[7]基于GLRT准则，提出了最优K检测器(Optimum K Detector, OKD)。针对逆伽马纹理结构下的广义Pareto分布，文献[8]提出了广义似然比线性门限检测器(Generalized Likelihood Ratio Test with Linear-Threshold Detector, GLRT-LTD)。值得注意的是，GLRT-LTD为广义Pareto分布下的最优相干检测器，且具有解析的门限结构。文献[9]在CG-IG杂波背景下，基于GLRT准则设计具有逆高斯纹理结构的广义似然比(Generalized Likelihood Ratio Test with Inverse Gaussian texture, GLRT-IG)检测器。文献[10]给出了CG-GIG杂波背景下的最优相干检测器，命名为GIG纹理下的广义似然比检测(Generalized Likelihood Ratio Test with Generalized Inverse Gaussian texture, GIG-GLRT)器。为了进一步提升检测器的性能，海杂波的频谱分布、纹理分量空间相关性和散斑协方差矩阵斜对称特征等先验信息被引入检测器的结构设计。文献[11]结合ANMF检测器，在CG-GIG杂波背景下提出近最优的alpha自适应匹配滤波(α Adaptive Matched Filter, α-AMF)检测器，实现了检测性能和计算复杂度的折衷。文献[12]将海杂波纹理分量的空间相关性引入GLRT检测器的设计，在均匀和部分均匀杂波背景下提升了目标检测性能。文献[13]将协方差矩阵的斜对称特性引入检测器设计中，一定程度上降低了检测器对参考单元杂波数据的依赖。然而，上述检测器均假定目标速度已知，并未考虑导向矢量失配对检测器的影响。在实际雷达对海探测过程中，受波束指向误差等不确定性因素的影响，实际导向矢量和理论导向矢量可能存在不匹配现象，如图1所示。导向矢量的失配可能会检测到非期望的副瓣目标，从而引起不必要的虚警。为了平衡检测器的目标检测性能以及抗失配能力，文献[14]在存在失配信号的条件下推导得到自适应波束形成器正交抑制检测(Adaptive Beamformer Orthogonal Rejection Test, ABORT)算法。文献[15]在CG-IG分布杂波背景下，提出类似于ABORT检测器结构的目标检测算法。针对距离扩展目标检测问题，文献[16]在K分布杂波背景下设计类似于ABORT的检测器。

图 1 失配示意图

下载: 全尺寸图片幻灯片

本文针对GIG纹理背景下的失配检测问题，分别基于GLRT和最大后验估计GLRT(Maximum A Posteriori GLRT, MAP GLRT)准则发展类似于ABORT的失配检测器。首先，通过在零假设中引入与理论导向矢量正交的虚拟信号来抑制失配信号。然后基于失配情况下的目标检测模型设计两个失配检测器，并证明其对散斑协方差矩阵和目标多普勒导向矢量均具有CFAR特性。最后仿真和实测海杂波数据实验结果表明，提出的失配检测器以牺牲匹配情况下较小的检测性能为代价，换取了较高的抗失配性能。相较于CG-GIG分布下的GIG-GLRT, MAP-GLRT^[17], α-AMF检测器以及不依赖纹理结构分布的ANMF和ABORT检测器，提出的失配检测器同时兼备较好的目标检测性能和抗失配性能。

2. 失配信号的目标检测模型

失配情况下的实际导向矢量 ${{\boldsymbol p}_{\mathrm{m}}}$ 通常由平行和正交于理论导向矢量 ${\boldsymbol{p}}$ 的 ${\boldsymbol p}_{\mathrm{m}}^\parallel$ 和 ${\boldsymbol p}_{\mathrm{m}}^ \bot$ 构成，即 ${{\boldsymbol p}_{\mathrm{m}}} = {\boldsymbol p}_{\mathrm{m}}^\parallel + {\boldsymbol p}_{\mathrm{m}}^ \bot$ 。其中， ${\boldsymbol p}_{\mathrm{m}}^\parallel = \alpha {\boldsymbol{p}}$ 是希望检测到的目标回波，而 ${\boldsymbol p}_{\mathrm{m}}^ \bot = \beta {\boldsymbol{q}}$ 是不希望检测到的干扰信号。为了提升检测器的抗失配性能，与ABORT检测器类似，引入虚拟信号 $\beta {\boldsymbol{q}}$ 建模零假设 ${H_0}$ 下的待检测(Cell Under Test, CUT)回波向量。因此在GIG纹理背景下导向矢量失配的目标检测模型定义为

$\left. {\begin{array}{*{20}{c}} {{H_0}:\left\{ {\begin{array}{*{20}{l}} {{\boldsymbol{z}} = \beta {\boldsymbol{q}} + {\boldsymbol{c}}} \\ {{{\boldsymbol{z}}_k} = {{\boldsymbol{c}}_k},k = 1,2, \cdots ,L} \end{array}} \right.} \\ {{H_1}:\left\{ {\begin{array}{*{20}{l}} {{\boldsymbol{z}} = \alpha {\boldsymbol{p}} + {\boldsymbol{c}}} \\ {{{\boldsymbol{z}}_k} = {{\boldsymbol{c}}_k},k = 1,2, \cdots ,L} \end{array}} \right.} \end{array}} \right\}$

(1)

其中，零假设 ${H_0}$ 表示CUT中不存在目标，而备择假设 ${H_1}$ 表示CUT中存在目标。向量 ${\boldsymbol{z}} = {\left[ {{z_1},{z_2}, \cdots ,{z_N}} \right]^{\mathrm{T}}}$ 表示CUT回波，其中包含 $N$ 个相参脉冲； ${{\boldsymbol{z}}_k},k = 1,2, \cdots ,L$ 表示CUT周围存在的 $L$ 个参考单元； $\alpha$ 表示目标回波幅度，通常被建模为Swerling 0型； ${\boldsymbol{p}} = \left[ {1,{{\mathrm{e}}^{{\mathrm{j}}2\pi {f_{\mathrm{d}}}}}, \cdots ,{{\mathrm{e}}^{{\mathrm{j}}2\pi \left( {N - 1} \right){f_{\mathrm{d}}}}}} \right]$ 表示目标回波的导向矢量，其归一化多普勒频率 ${f_{\mathrm{d}}}$ 在区间 $\left[ { - 0.5,0.5} \right]$ 内服从均匀分布； $\tau$ 和 ${{u}}$ 分别表示杂波的纹理分量和散斑分量； $\beta {\boldsymbol{q}}$ 表示人为引入 ${H_0}$ 的虚拟信号， $\beta$ 和 ${\boldsymbol{q}}$ 分别表示虚拟信号的复幅度和导向矢量，虚拟导向矢量 ${\boldsymbol{q}}$ 在白化空间中与理论导向矢量 ${\boldsymbol{p}}$ 正交，即 $\left\langle {{{\boldsymbol{M}}^{ - 1/2}}{\boldsymbol{p}}} \right\rangle = {\left\langle {{{\boldsymbol{M}}^{ - 1/2}}{\boldsymbol{q}}} \right\rangle ^ \bot }$ 。需要注意的是，在ABORT检测器中，所引入的虚拟信号 ${\boldsymbol{q}}$ 与理论导向矢量 ${\boldsymbol{p}}$ 在准白化空间中正交，即 $\left\langle {{{\boldsymbol{S}}^{ - 1/2}}{\boldsymbol{p}}} \right\rangle = {\left\langle {{{\boldsymbol{S}}^{ - 1/2}}{\boldsymbol{q}}} \right\rangle ^ \bot }$ ，其中 ${\boldsymbol{S}}$ 表示理论散斑协方差矩阵 ${\boldsymbol{M}}$ 的估计值。根据球不变随机向量模型，CUT回波向量在 ${H_0}$ 和 ${H_1}$ 假设下的条件PDF定义为

$\begin{split} \,&{f_0}\left( {{\boldsymbol{z}}|\beta ,{\boldsymbol{M}},\tau ;{H_0}} \right) = \frac{1}{{{{\left( {\pi \tau } \right)}^N}\left| {\boldsymbol{M}} \right|}}\\ & \qquad \cdot \exp \left( { - \frac{{{{\left( {{\boldsymbol{z}} - \beta {\boldsymbol{q}}} \right)}^{\mathrm{H}}}{{\boldsymbol{M}}^{ - 1}}\left( {{\boldsymbol{z}} - \beta {\boldsymbol{q}}} \right)}}{\tau }} \right) \end{split}$

(2)

$\begin{split} & {f_1}\left( {{\boldsymbol{z}}|\alpha ,{\boldsymbol{M}},\tau ;{H_1}} \right) = \frac{1}{{{{\left( {\pi \tau } \right)}^N}\left| {\boldsymbol{M }}\right|}}\\ & \qquad \cdot\exp \left( { - \frac{{{{\left( {{\boldsymbol{z}} - \alpha {\boldsymbol{p}}} \right)}^{\mathrm{H}}}{{\boldsymbol{M}}^{ - 1}}\left( {{\boldsymbol{z}} - \alpha {\boldsymbol{p}}} \right)}}{\tau }} \right) \end{split}$

(3)

其中， $\tau$ 表示GIG分布的随机数，其PDF表达式 $f\left( \tau \right)$ 为

$f\left( \tau \right) = \frac{{{{\left( {a/b} \right)}^{p/2}}{\tau ^{p - 1}}}}{{2{K_p}\left( {\sqrt {ab} } \right)}}{{\mathrm{e}}^{ - \frac{{\left( {a\tau + b/\tau } \right)}}{2}}},a \gt 0,b \gt 0$

(4)

其中， $a$ 和 $b$ 分别表示GIG分布的形状参数和尺度参数， ${K_p}\left( \cdot \right)$ 表示 $p$ 阶第2类修正Bessel函数。伽马和逆伽马分布分别对应形状参数 $a = 0$ 和尺度参数 $b = 0$ 时的GIG分布；当 $p = - {1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}$ , $b = \lambda u$ 和 $a = {\lambda \mathord{\left/ {\vphantom {\lambda u}} \right. } u}$ 时，GIG分布退化为形状参数为 $\lambda$ 尺度参数为 $\mu$ 的GIG分布。

3. 基于GIG纹理结构的失配检测器

本节考虑导向矢量失配的影响，基于两步GLRT和MAP-GLRT准则在广义逆高斯杂波背景下设计失配检测器，并通过理论证明所提检测器对散斑协方差矩阵 ${\boldsymbol{M}}$ 以及目标多普勒导向矢量 ${\boldsymbol{p}}$ 具有CFAR特性。

3.1 GLRT准则下的失配检测器

本节基于两步GLRT准则设计失配检测器。首先假设杂波散斑分量 ${{u}}$ 的协方差矩阵 ${\boldsymbol{M}}$ 已知，那么导向矢量失配情况下的GLRT准则定义为

$\frac{{\mathop {\max }\limits_\alpha \displaystyle\int {{f_1}\left( {{\boldsymbol{z}}|\alpha ,\tau ;{H_1}} \right)f\left( \tau \right){\mathrm{d}}\tau } }}{{\mathop {\max }\limits_\beta \displaystyle\int {{f_0}\left( {{\boldsymbol{z}}|\beta ,\tau ;{H_0}} \right)f\left( \tau \right){\mathrm{d}}\tau } }}\mathop \gtrless \limits_{{H_0}}^{{H_1}} {\xi _{\mathrm{g}}}$

(5)

其中， ${\xi _{\mathrm{g}}}$ 为GLRT检测器的检测门限，式(5)中参数 $\alpha$ 和 $\beta$ 的最大似然(Maximum Likelihood, ML)估计为 $\hat \alpha$ 与 $\hat \beta$ ，定义分别为

$\hat \alpha = \frac{{{{\boldsymbol{p}}^{\mathrm{H}}}{{\boldsymbol{M}}^{ - 1}}{\boldsymbol{z}}}}{{{{\boldsymbol{p}}^{\mathrm{H}}}{{\boldsymbol{M}}^{ - 1}}{\boldsymbol{p}}}}$

(6)

$\hat \beta = \frac{{{{\boldsymbol{q}}^{\mathrm{H}}}{{\boldsymbol{M}}^{ - 1}}{\boldsymbol{z}}}}{{{{\boldsymbol{q}}^{\mathrm{H}}}{{\boldsymbol{M}}^{ - 1}}{\boldsymbol{q}}}}$

(7)

将式(6)和式(7)分别代入式(2)和式(3)的条件PDF中，可以得到

$\begin{split} \mathop {\max }\limits_\beta {f_0}\left( {{\boldsymbol{z}}|\beta ,\tau ;{H_0}} \right) = \,& {f_0}\left( {{\boldsymbol{z}}|\hat \beta ,\tau ;{H_0}} \right) = \frac{1}{{{{\left( {\pi \tau } \right)}^N}\left| {\boldsymbol{M}} \right|}}\\ & \cdot\exp \left( { - \frac{{{l_0}}}{\tau }} \right)\\[-1pt] \end{split}$

(8)

$\begin{split} \mathop {\max }\limits_\alpha {f_1}\left( {{\boldsymbol{z}}|\alpha ,\tau ;{H_1}} \right) =\,& {f_0}\left( {{\boldsymbol{z}}|\hat \alpha ,\tau ;{H_1}} \right) = \frac{1}{{{{\left( {\pi \tau } \right)}^N}\left| {\boldsymbol{M}} \right|}}\\ & \cdot \exp \left( { - \frac{{{l_1}}}{\tau }} \right) \\[-1pt] \end{split}$

(9)

其中 ${l_0}$ 和 ${l_1}$ 分别定义为

${l_0} = {{\boldsymbol{z}}^{\mathrm{H}}}{{\boldsymbol{M}}^{ - 1}}{\boldsymbol{z}} - \frac{{{{\left| {{{\boldsymbol{q}}^{\mathrm{H}}}{{\boldsymbol{M}}^{ - 1}}{\boldsymbol{z}}} \right|}^2}}}{{{{\boldsymbol{q}}^{\mathrm{H}}}{{\boldsymbol{M}}^{ - 1}}{\boldsymbol{q}}}}$

(10)

${l_1} = {{\boldsymbol{z}}^{\mathrm{H}}}{{\boldsymbol{M}}^{ - 1}}{\boldsymbol{z}} - \frac{{{{\left| {{{\boldsymbol{p}}^{\mathrm{H}}}{{\boldsymbol{M}}^{ - 1}}{\boldsymbol{z}}} \right|}^2}}}{{{{\boldsymbol{p}}^{\mathrm{H}}}{{\boldsymbol{M}}^{ - 1}}{\boldsymbol{p}}}}$

(11)

设矩阵 ${\boldsymbol{A}}$ 的投影矩阵为 ${{\boldsymbol{P}}_{\boldsymbol{A}}} = {\boldsymbol{A}}{\left( {{{\boldsymbol{A}}^{\mathrm{H}}}{\boldsymbol{A}}} \right)^{ - 1}}{{\boldsymbol{A}}^{\mathrm{H}}}$ ，那么投影矩阵 ${{\boldsymbol{P}}_{\boldsymbol{A}}}$ 的正交补矩阵定义为 ${\boldsymbol{P}}_{\boldsymbol{A}}^ \bot = {\boldsymbol{I}} - {{\boldsymbol{P}}_{\boldsymbol{A}}} = {\boldsymbol{I}} - {\boldsymbol{A}}{\left( {{{\boldsymbol{A}}^{\mathrm{H}}}{\boldsymbol{A}}} \right)^{ - 1}}{{\boldsymbol{A}}^{\mathrm{H}}}$ ，其中 ${\boldsymbol{I}}$ 表示与矩阵 ${\boldsymbol{A}}$ 行维数相同的单位矩阵。利用虚拟导向矢量 ${\boldsymbol{q}}$ 与理论导向矢量 ${\boldsymbol{p}}$ 在白化空间中相互正交的性质， $\left\langle {{{\boldsymbol{M}}^{ - 1/2}}{\boldsymbol{p}}} \right\rangle = {\left\langle {{{\boldsymbol{M}}^{ - 1/2}}{\boldsymbol{q}}} \right\rangle ^ \bot }$ 可以改写为

${{\boldsymbol{I}}_N} - \frac{{{{\boldsymbol{M}}^{ - 1/2}}{\boldsymbol{q}}{{\boldsymbol{q}}^{\mathrm{H}}}{{\boldsymbol{M}}^{ - 1/2}}}}{{{{\boldsymbol{q}}^{\mathrm{H}}}{{\boldsymbol{M}}^{ - 1}}{\boldsymbol{q}}}} = \frac{{{{\boldsymbol{M}}^{ - 1/2}}{\boldsymbol{p}}{{\boldsymbol{p}}^{\mathrm{H}}}{{\boldsymbol{M}}^{ - 1/2}}}}{{{{\boldsymbol{p}}^{\mathrm{H}}}{{\boldsymbol{M}}^{ - 1}}{\boldsymbol{p}}}}$

(12)

将式(12)代入式(10)，化简 ${l_0}$ 为

$\begin{split} {l_0} \,& = {{\boldsymbol{z}}^{\mathrm{H}}}{{\boldsymbol{M}}^{ - 1/2}}\left( {{{\boldsymbol{I}}_N} - \frac{{{{\boldsymbol{M}}^{ - 1/2}}{\boldsymbol{q}}{{\boldsymbol{q}}^{\mathrm{H}}}{{\boldsymbol{M}}^{ - 1/2}}}}{{{{\boldsymbol{q}}^{\mathrm{H}}}{{\boldsymbol{M}}^{ - 1}}{\boldsymbol{q}}}}} \right){{\boldsymbol{M}}^{ - 1/2}}{\boldsymbol{z}} \\ & = \frac{{{{\left| {{{\boldsymbol{p}}^{\mathrm{H}}}{{\boldsymbol{M}}^{ - 1}}{\boldsymbol{z}}} \right|}^2}}}{{{{\boldsymbol{p}}^{\mathrm{H}}}{{\boldsymbol{M}}^{ - 1}}{\boldsymbol{p}}}}\\[-1pt] \end{split}$

(13)

将式(4)和式(9)代入式(5)，即可得到式(5)分子部分关于纹理分量 $\tau$ 的积分结果，如式(14)所示

$\begin{split} \int {{f_1}\left( {{\boldsymbol{z}}|\hat \alpha ,\tau ;{H_1}} \right)f\left( \tau \right){\mathrm{d}}\tau } =\,& \frac{{{a^{{N / 2}}}{b^{ - {p / 2}}}{{\left( {2{l_1} + b} \right)}^{\frac{{p - N}}{2}}}}}{{{\pi ^N}\left| {\boldsymbol{M}} \right| {K_p}\left( {\sqrt {ab} } \right)}}\\ & \cdot{K_{N - p}}\left( {\sqrt {a\left( {2{l_1} + b} \right)} } \right) \end{split}$

(14)

同理，式(5)分母部分关于纹理分量 $\tau$ 的积分结果如式(15)所示

$\begin{split} \int {{f_0}\left( {{\boldsymbol{z}}|\hat \beta ,\tau ;{H_0}} \right)f\left( \tau \right){\mathrm{d}}\tau } =\,& \frac{{{a^{{N / 2}}}{b^{ - {p / 2}}}{{\left( {2{l_0} + b} \right)}^{\frac{{p - N}}{2}}}}}{{{\pi ^N}\left| {\boldsymbol{M}} \right|{K_p}\left( {\sqrt {ab} } \right)}}\\ & \cdot {K_{N - p}}\left( {\sqrt {a\left( {2{l_0} + b} \right)} } \right) \end{split}$

(15)

将式(14)和式(15)代入式(5)，即可得到广义逆高斯纹理复合高斯杂波背景下类似于AOBRT(ABORT-like Generalized Compound Gaussian with Inverse Gaussian, A-CGGIG)检测器的检验统计量，如式(16)所示

$\frac{{{{\left( {2{l_1} + b} \right)}^{\frac{{P - N}}{2}}}{K_{N - p}}\left( {\sqrt {a\left( {2{l_1} + b} \right)} } \right)}}{{{{\left( {2{l_0} + b} \right)}^{\frac{{P - N}}{2}}}{K_{N - p}}\left( {\sqrt {a\left( {2{l_0} + b} \right)} } \right)}}\mathop \gtrless \limits_{{H_0}}^{{H_1}} {\xi _{\mathrm{g}}}$

(16)

第2步，利用参考单元数据 ${{\boldsymbol{z}}_k},k = 1,2, \cdots ,L$ 估计散斑分量的协方差矩阵实现自适应相干检测。本文采用约束渐进最大似然估计器(Constrained Approximate Maximum Likelihood Estimator, CAMLE)估计 ${\boldsymbol{M}}$ ，降低海杂波纹理分量 $\tau$ 对散斑协方差矩阵估计的影响，CAMLE的定义为

$\left. \begin{aligned} & {\hat {\boldsymbol{M}}\left( m \right) = \frac{N}{L}\sum\limits_{k = 1}^L {\frac{{{{\boldsymbol{z}}_k}{\boldsymbol{z}}_k^{\mathrm{H}}}}{{{\boldsymbol{z}}_k^{\mathrm{H}}{{\hat {\boldsymbol{M}}}^{ - 1}}\left( m \right){{\boldsymbol{z}}_k}}}} } \\ & {\hat {\boldsymbol{M}}\left( m \right) = \frac{N}{{{\mathrm{Tr}}\left( {\hat {\boldsymbol{M}}\left( m \right)} \right)}}{{\hat {\boldsymbol{M}}}_{{\mathrm{CAMLE}}}}\left( {{m}} \right)} \end{aligned} \right\}$

(17)

其中， $m$ 表示迭代次数，运算符 ${\mathrm{Tr}}\left( \cdot \right)$ 表示求矩阵的迹。CAMLE的初值通常采用归一化样本协方差矩阵(Normalized Sample Covariance Matrix, NSCM)，通过1～3次迭代即可得到较为准确的估计结果。将散斑协方差矩阵的估计结果 $\hat {\boldsymbol{M}}$ 带入式(16)，得到自适应的A-CGGIG检测器

$\frac{{{{\left( {2{{\hat l}_1} + b} \right)}^{\frac{{P - N}}{2}}}{K_{N - p}}\left( {\sqrt {a\left( {2{{\hat l}_1} + b} \right)} } \right)}}{{{{\left( {2{{\hat l}_0} + b} \right)}^{\frac{{P - N}}{2}}}{K_{N - p}}\left( {\sqrt {a\left( {2{{\hat l}_0} + b} \right)} } \right)}}\mathop \gtrless \limits_{{H_0}}^{{H_1}} {\xi _{\mathrm{g}}}$

(18)

其中， ${\hat l_1}$ 和 ${\hat l_0}$ 分别是 ${l_1}$ 和 ${l_0}$ 的估计值，定义为

$\left. \begin{aligned} & {{{\hat l}_1} = {{\boldsymbol{z}}^{\mathrm{H}}}{{\hat {\boldsymbol{M}}}^{ - 1}}{\boldsymbol{z}} - {{\left| {{{\boldsymbol{p}}^{\mathrm{H}}}{{\hat {\boldsymbol{M}}}^{ - 1}}{\boldsymbol{z}}} \right|}^2}/\left( {{{\boldsymbol{p}}^{\mathrm{H}}}{{\hat {\boldsymbol{M}}}^{ - 1}}{\boldsymbol{p}}} \right)} \\ & {{{\hat l}_0} = {{\left| {{{\boldsymbol{p}}^{\mathrm{H}}}{{\hat {\boldsymbol{M}}}^{ - 1}}{\boldsymbol{z}}} \right|}^2}/\left( {{{\boldsymbol{p}}^{\mathrm{H}}}{{\hat {\boldsymbol{M}}}^{ - 1}}{\boldsymbol{p}}} \right)} \end{aligned} \right\}$

(19)

3.2 MAP GLRT准则下的失配检测器

本节基于两步MAP GLRT准则设计失配检测器。首先假设杂波散斑分量 ${{u}}$ 的协方差矩阵 ${\boldsymbol{M}}$ 已知，那么导向矢量失配情况下的MAP GLRT准则定义为

$\frac{{\mathop {\max }\limits_{{\tau _1},\alpha } {f_1}\left( {{\boldsymbol{z}}|\alpha ,{\tau _1};{H_1}} \right)f\left( \tau \right)}}{{\mathop {\max }\limits_{{\tau _0},\beta } {f_0}\left( {{\boldsymbol{z}}|\beta ,{\tau _0};{H_0}} \right)f\left( \tau \right)}}\mathop \gtrless \limits_{{H_0}}^{{H_1}} {\xi _{\mathrm{h}}}$

(20)

其中， ${\xi _{\mathrm{h}}}$ 表示MAP GLRT检测器的检测门限。参数 $\alpha$ 和 $\beta$ 的ML估计 $\hat \alpha$ 和 $\hat \beta$ 同式(6)和式(7)。将 $\hat \alpha$ 和 $\hat \beta$ 代入式(2)和式(3)，并乘以式(4)所示的广义逆高斯分布，即可得到式(21)和式(22)

$\begin{split} {f_0}\left( {{\boldsymbol{z}}|\hat \beta ,{\tau _0};{H_0}} \right)f\left( {{\tau _0}} \right) =\,& \frac{{{{\left( {a/b} \right)}^{p/2}}}}{{2{K_p}\left( {\sqrt {ab} } \right){\pi ^N}\left| {\boldsymbol{M}} \right|}}\tau _0^{p - N - 1}\\ & \cdot \exp \left[ { - \left( {\frac{a}{2}{\tau _0} + \frac{{{l_0} + {b \mathord{\left/ {\vphantom {b 2}} \right. } 2}}}{{{\tau _0}}}} \right)} \right] \end{split}$

(21)

$\begin{split} {f_1}\left( {{\boldsymbol{z}}|\hat \alpha ,{\tau _1};{H_1}} \right)f\left( {{\tau _1}} \right) =\,& \frac{{{{\left( {a/b} \right)}^{p/2}}}}{{2{K_p}\left( {\sqrt {ab} } \right){\pi ^N}\left| {\boldsymbol{M}} \right|}}\tau _1^{p - N - 1}\\ & \cdot\exp \left[ { - \left( {\frac{a}{2}{\tau _1} + \frac{{{l_1} + {b \mathord{\left/ {\vphantom {b 2}} \right. } 2}}}{{{\tau _1}}}} \right)} \right] \end{split}$

(22)

分别计算式(21)和式(22)的对数表达式，并对纹理分量 $\tau$ 求偏导，将其结果置0，即可以分别得到零假设 ${H_0}$ 和备择假设 ${H_1}$ 下纹理分量的MAP估计 ${\hat \tau _0}$ 和 ${\hat \tau _1}$ ，如式(23)和式(24)所示

$\begin{split} \quad\;\; {\hat \tau _0} = \,& \left( - \left( {N + 1 - p} \right) + {\mathrm{sqrt}}\Bigr( {{\left( {N + 1 - p} \right)}^2} \right.\\ & \left.+ a\left( {2{l_0} + b} \right) \Bigr) \right)/a \end{split}$

(23)

$\begin{split} \quad \;\; {\hat \tau _1} =\,& \left( - \left( {N + 1 - p} \right) + {\mathrm{sqrt}}\Bigr( {{\left( {N + 1 - p} \right)}^2} \right.\\ & \left. + a\left( {2{l_1} + b} \right) \Bigr) \right)/a \end{split}$

(24)

将式(11)、式(13)、式(23)和式(24)代入式(20)，即可得到基于MAP准则的广义逆高斯纹理复合高斯杂波背景下类似于AOBRT(ABORT-like Generalized Compound Gaussian with Inverse Gaussian, AM-CGGIG)检测器的检验统计量，如式(25)所示

$\frac{{{{\hat \tau }_1}^{p - N - 1}\exp \left[ { - \left( {\dfrac{a}{2}{{\hat \tau }_1} + \dfrac{{{l_1} + {b \mathord{\left/ {\vphantom {b 2}} \right. } 2}}}{{{{\hat \tau }_1}}}} \right)} \right]}}{{{{\hat \tau }_0}^{p - N - 1}\exp \left[ { - \left( {\dfrac{a}{2}{{\hat \tau }_0} + \dfrac{{{l_0} + {b \mathord{\left/ {\vphantom {b 2}} \right. } 2}}}{{{{\hat \tau }_0}}}} \right)} \right]}}\mathop \gtrless \limits_{{H_2}}^{{H_1}} {\xi _{\mathrm{h}}}$

(25)

第2步，将式(17)中的CAMLE的估计结果代入式(25)，即可得到AM-CGGIG检测器的自适应形式，定义为

$\frac{{{{\hat \tau }_1}^{p - N - 1}\exp \left[ { - \left( {\dfrac{a}{2}{{\hat \tau }_1} + \dfrac{{{{\hat l}_1} + {b \mathord{\left/ {\vphantom {b 2}} \right. } 2}}}{{{{\hat \tau }_1}}}} \right)} \right]}}{{{{\hat \tau }_0}^{p - N - 1}\exp \left[ { - \left( {\dfrac{a}{2}{{\hat \tau }_0} + \dfrac{{{{\hat l}_0} + {b \mathord{\left/ {\vphantom {b 2}} \right. } 2}}}{{{{\hat \tau }_0}}}} \right)} \right]}}\mathop \gtrless \limits_{{H_2}}^{{H_1}} {\xi _{\mathrm{h}}}$

(26)

对比式(18)和式(26)的检验统计量可以发现，次优的AM-CGGIG检测器由于利用MAP算法估计海杂波的纹理结构，其检验统计量中不存在第2类修正Bessel函数，因此降低了计算复杂度。

3.3 失配检测器的CFAR特性分析

本节分析提出的A-CGGIG和AM-CGGIG检测器的CFAR特性。可以发现式(18)和式(26)均为变量 ${\hat l_0}$ 和 ${\hat l_1}$ 的函数，因此仅需分析 ${\hat l_0}$ 和 ${\hat l_1}$ 与失配检测器的关系即可，首先将这两个变量分别重写为

${\hat l_0} = \frac{{{{\left| {{{\left( {{{\boldsymbol{M}}^{ - {1 / 2}}}{\boldsymbol{p}}} \right)}^{\mathrm{H}}}{{\boldsymbol{M}}^{{1 / 2}}}{{\hat M}^{ - 1}}{{\boldsymbol{M}}^{{1 / 2}}}\left( {{{\boldsymbol{M}}^{ - {1 / 2}}}{\boldsymbol{z}}} \right)} \right|}^2}}}{{{{\left( {{{\boldsymbol{M}}^{ - {1/ 2}}}{\boldsymbol{p}}} \right)}^{\mathrm{H}}}{{\boldsymbol{M}}^{{1/ 2}}}{{\hat {\boldsymbol{M}}}^{ - 1}}{{\boldsymbol{M}}^{{1 / 2}}}\left( {{{\boldsymbol{M}}^{ - {1 / 2}}}{\boldsymbol{p}}} \right)}}$

(27)

${\hat l_1} = {\left( {{{\boldsymbol{M}}^{ - {1 / 2}}}{\boldsymbol{z}}} \right)^{\mathrm{H}}}{{\boldsymbol{M}}^{{1 / 2}}}{\hat {\boldsymbol{M}}^{ - 1}}{{\boldsymbol{M}}^{{1 / 2}}}\left( {{{\boldsymbol{M}}^{ - {1 / 2}}}{\boldsymbol{z}}} \right) - {\hat l_0}$

(28)

对于重塑后的导向矢量 ${{\boldsymbol{M}}^{{{ - 1} \mathord{\left/ {\vphantom {{ - 1} 2}} \right. } 2}}}{\boldsymbol{p}}$ 而言，总存在一个Householder矩阵 ${\boldsymbol{P}}$ ，能够将其转化为 ${\boldsymbol{P}}{{\boldsymbol{M}}^{{{ - 1} / 2}}}{\boldsymbol{p}} = \left| {{{\boldsymbol{M}}^{{{ - 1} / 2}}}{\boldsymbol{p}}} \right|{\boldsymbol{v}}$ ，其中 ${\boldsymbol{v}} = {\left( {1,0,0, \cdots ,0} \right)^{\mathrm{T}}}$ 。基于球不变随机向量模型，分别化简 ${\hat l_0}$ 和 ${\hat l_1}$ 为

$\begin{split} & {\hat l_0} \\ & = \frac{{{{\left| {{{\left( {{\boldsymbol{P}}{{\boldsymbol{M}}^{ - {1 / 2}}}{\boldsymbol{p}}} \right)}^{\mathrm{H}}} {\boldsymbol{P}}{{\boldsymbol{M}}^{{1 / 2}}} {{\hat {\boldsymbol{M}}}^{ - 1}} {{\boldsymbol{M}}^{{1 / 2}}}{{\boldsymbol{P}}^{\mathrm{H}}} \left( {{\boldsymbol{P}}{{\boldsymbol{M}}^{ - {1 / 2}}}{\boldsymbol{z}}} \right)} \right|}^2}}}{{{{\left( {{\boldsymbol{P}}{{\boldsymbol{M}}^{ - {1 /2}}}{\boldsymbol{p}}} \right)}^{\mathrm{H}}}{\boldsymbol{P}}{{\boldsymbol{M}}^{{1 / 2}}}{{\hat{\boldsymbol{M}}}^{ - 1}} {{\boldsymbol{M}}^{{1 / 2}}}{{\boldsymbol{P}}^{\mathrm{H}}} \left( {{\boldsymbol{P}}{{\boldsymbol{M}}^{ - {1 / 2}}}{\boldsymbol{p}}} \right)}} \\ & = \tau \frac{{{{\left| {{{\boldsymbol{v}}^{\mathrm{H}}}{{\bar {\boldsymbol{M}}}^{ - 1}}{\boldsymbol{n}}} \right|}^2}}}{{{{\boldsymbol{v}}^{\mathrm{H}}}{{\bar {\boldsymbol{M}}}^{ - 1}}{\boldsymbol{v}}}}\\[-1pt] \end{split}$

(29)

$\begin{split} {\hat l_1} = \,&{\left( {{\boldsymbol{P}}{{\boldsymbol{M}}^{ - {1 / 2}}}{\boldsymbol{z}}} \right)^{\mathrm{H}}}{\boldsymbol{P}}{{\boldsymbol{M}}^{{1 / 2}}}{\hat{\boldsymbol{ M}}^{ - 1}}{{\boldsymbol{M}}^{{1 / 2}}}{{\boldsymbol{P}}^{\mathrm{H}}}\\ & \cdot\left( {{\boldsymbol{P}}{{\boldsymbol{M}}^{ - {1 / 2}}}{\boldsymbol{z}}} \right) - {\hat l_0} = \tau \left( {{{\boldsymbol{n}}^{\mathrm{H}}}{{\bar {\boldsymbol{M}}}^{ - 1}}{\boldsymbol{n}}} \right) - {\hat l_0} \end{split}$

(30)

其中， ${\boldsymbol{n}} = {\boldsymbol{P}}{{\boldsymbol{M}}^{ - {1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}}}{{u}}$ 和 $\bar {\boldsymbol{M}} = {\boldsymbol{P}}{{\boldsymbol{M}}^{ - {1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}}}\hat {\boldsymbol{M}}{{\boldsymbol{M}}^{ - {1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}}} {{\boldsymbol{P}}^{\mathrm{H}}}$ 。杂波向量 ${\boldsymbol{n}}$ 服从 ${\boldsymbol{n}}{\text{～}}{\mathrm{CN}}\left( {0,{\boldsymbol{I}}} \right)$ 。此外，由于采用CAMLE，矩阵 $\bar {\boldsymbol{M}}$ 与理论散斑协方差矩阵 ${\boldsymbol{M}}$ 之间相互独立^[18]。将式(29)和式(30)代入式(18)和式(26)中，可以发现检测器A-CGGIG和AM-CGGIG的检验统计量均独立于 ${\boldsymbol{M}}$ 和 ${\boldsymbol{p}}$ ，这表明A-CGGIG和AM-CGGIG检测器对 ${\boldsymbol{M}}$ 和 ${\boldsymbol{p}}$ 具有CFAR特性。

4. 失配检测器性能评估

本节在导向矢量匹配和出现失配两种情况下，利用仿真和实测海杂波数据验证提出的A-CGGIG和AM-CGGIG检测器性能及其CFAR特性。理论导向矢量 ${\boldsymbol{p}}$ 和实际导向矢量 ${{\boldsymbol{p}}_{\mathrm{m}}}$ 失配角 $\theta$ 余弦值的平方定义为

${\cos ^2}\theta = \frac{{{{\left| {{{\boldsymbol{p}}^{\mathrm{H}}}{{\boldsymbol{M}}^{ - 1}}{{\boldsymbol{p}}_{\mathrm{m}}}} \right|}^2}}}{{\left( {{{\boldsymbol{p}}^{\mathrm{H}}}{{\boldsymbol{M}}^{ - 1}}{\boldsymbol{p}}} \right)\left( {{\boldsymbol{p}}_{\mathrm{m}}^{\mathrm{H}}{{\boldsymbol{M}}^{ - 1}}{{\boldsymbol{p}}_{\mathrm{m}}}} \right)}}$

(31)

当 ${\cos ^2}\theta = 1$ 时，表示实际导向矢量 ${{\boldsymbol{p}}_{\mathrm{m}}}$ 和理论导向矢量 ${\boldsymbol{p}}$ 相匹配；当 $0 \lt {\cos ^2}\theta \lt 1$ 时，表示理论导向矢量 ${\boldsymbol{p}}$ 失配于实际导向矢量 ${{\boldsymbol{p}}_{\mathrm{m}}}$ 。

4.1 仿真数据实验结果及性能分析

在仿真杂波数据实验中，散斑协方差矩阵 ${\boldsymbol{M}}$ 被建模为1阶滞后衰减系数为 $\rho$ 的指数相关型协方差矩阵，即 ${{\boldsymbol{M}}_{i,j}} = {\rho ^{\left| {i - j} \right|}},1 \le i,j \le N$ ，对于海杂波而言 $\rho \in \left[ {0.9,0.99} \right]$ ，本节设定 $\rho = 0.9$ 。此外，设定脉冲积累数目 $N = 6$ ，参考单元数目 $L = 20$ ，仿真杂波参数为 $a = 1,b = 1,p = 2$ ；目标归一化多普勒频率为 ${f_{\mathrm{d}}} = 0.3$ ，目标检测中的虚警概率为 ${P_{{\mathrm{fa}}}} = {10^{ - 4}}$ ，蒙特卡洛实验次数为10⁵，信杂比(Signal Clutter Ratio, SCR)为

${\mathrm{SCR}} = 10\lg \left( {\frac{{{{\left| \alpha \right|}^2}}}{{{S_\rho }\left( {{f_{\mathrm{d}}}} \right)}}} \right)$

(32)

其中， ${S_\rho }\left( {{f_{\mathrm{d}}}} \right)$ 表示海杂波的多普勒功率谱密度，定义为

${S_\rho }\left( {{f_{\mathrm{d}}}} \right) = 1 + 2\sum\limits_{n = 1}^{N - 1} {\left( {1 - \frac{n}{N}} \right){\rho ^n}} \cos \left( {2\pi {f_{\mathrm{d}}}n} \right)$

(33)

展示了提出的A-CGGIG和AM-CGGIG检测器在导向矢量匹配时，即 ${{\boldsymbol{p}}_{\mathrm{m}}} = {\boldsymbol{p}}$ ，与ANMF, ABORT, GIG-GLRT, MAP-GLRT以及α-AMF检测器的检测性能曲线；当信杂比 ${\mathrm{SCR}} = 6\;{\mathrm{dB}}$ 时，展示了提出的检测器和对比检测器的接收机工作特性(Receiver Operating Characteristic, ROC)曲线。列举了A-CGGIG和AM-CGGIG检测器与对比检测器目标检测概率达到 $90\%$ 时所需要的SCR。

图 2 在匹配情况下提出的检测器与对比检测器的目标检测性能

下载: 全尺寸图片幻灯片

表 1 导向矢量匹配时上述检测器检测概率达到

$90\%$ 时所需要的SCR(dB)

检测器类型	ANMF	ABORT	GIG-GLRT	MAP-GLRT	α-AMF	A-CGGIG	AM-GIGIG
信杂比	8.40	7.61	5.70	5.71	5.82	5.80	5.93

下载: 导出CSV

| 显示表格

综合分析图2和表1，可以发现提出的A-CGGIG和AM-CGGIG检测器在导向矢量匹配时的检测性能与GIG-GLRT和MAP-GLRT检测器的差距小于0.5 dB，且优于α-AMF, ANMF以及ABORT检测器。

在导向矢量出现失配时，所检测到的目标并非是期望方向上的目标，此时希望检测器的性能越差越好。提出的A-CGGIG和AM-CGGIG检测器与对比检测器的抗失配性能对比如图3所示，图3(a)表征检测器在不同SCR和失配角下抗失配性能的台面图，图3中0.9, 0.5, 0.1表示检测概率。从图3中可以发现提出的A-CGGIG和AM-CGGIG的抗失配性能优于GIG-GLRT, MAP-GLRT以及α-AMF检测器。

图 3 在失配情况下提出的检测器与对比检测器的目标检测性能

下载: 全尺寸图片幻灯片

沿着纵坐标取点，即可获取检测器在相应失配角下的检测性能曲线。和分别展示了当 ${\cos ^2}\theta = 0.8$ 和 ${\cos ^2}\theta = 0.6$ 时，上述检测器的性能曲线。从图3中可以发现，所有检测器的性能均会随着失配角度的增加而下降。一般而言，检测器的抗失配性能提升是以牺牲检测性能为代价的，而提出的检测器在导向矢量匹配时的检测性能与最优检测器差距小于0.5 dB，且对导向矢量的失配较为敏感，属于同时兼顾目标检测性能和抗失配能力的检测器。

分析了所提A-CGGIG和AM-CGGIG检测器关于目标多普勒频率 ${f_{\mathrm{d}}}$ 和1阶迟滞衰减系数 $\rho$ 的CFAR特性。从和中可以发现，自适应A-CGGIG和AM-CGGIG检测器的 ${P_{{\mathrm{fa}}}}$ 在不同 ${f_{\mathrm{d}}}$ 和 $\rho$ 保持为10^–4不变，与预先设置的 ${P_{{\mathrm{fa}}}}$ 一致。这表明 ${f_{\mathrm{d}}}$ 和 $M$ 不影响自适应A-CGGIG和AM-CGGIG检测器的门限大小，表明提出的检测器对 ${f_{\mathrm{d}}}$ 和 $M$ 具有CFAR特性。

图 4 检测器虚警概率曲线

下载: 全尺寸图片幻灯片

4.2 实测数据实验结果及性能分析

在实测数据实验中，IPIX雷达1998年19980205_171437_ANTSTEP文件VV极化数据被用于验证提出的检测器性能。19980205_171437_ANTSTEP数据共具有28个距离单元，为了方便实验，剔除第6～9受目标回波影响的距离单元，在剩余的纯海杂波数据中添加相干脉冲数目 $N = 6$ 的仿真目标信号，通过蒙特-卡洛实验获得检测概率曲线。展示了CG-GIG分布对实测海杂波数据的拟合结果，可以发现相较于K分布、广义Pareto分布以及CG-IG分布，CG-GIG分布具有较好的拟合效果，其参数估计结果为 $\hat a = 0.15,\hat b = 1.35, \hat p = - 0.7$ 。图5(b)给出了所提A-CGGIG和AM-CGGIG检测器与对比检测器在导向矢量匹配情况时的目标检测性能。

图 5 实测海杂波数据背景下检测器性能分析

下载: 全尺寸图片幻灯片

列举了A-CGGIG和AM-CGGIG检测器与对比检测器目标检测概率达到 $90\%$ 时所需要的SCR。从图5(b)和表2可以发现，A-CGGIG和AM-CGGIG检测器与GIG-GLRT检测器在导向矢量匹配时性能差距小于0.5 dB。

表 2 导向矢量匹配时上述检测器检测概率达到

$90\%$ 时所需要的SCR(dB)

检测器类型	ANMF	ABORT	GIG-GLRT	MAP-GLRT	α-AMF		A-CGGIG	AM-GIGIG
信杂比	12.88	12.16	11.34	11.50	11.43		11.37	11.53

下载: 导出CSV

| 显示表格

图6分析了A-CGGIG和AM-CGGIG检测器与对比检测器在导向矢量失配情况下的检测性能。从图6中可以发现，提出的检测器的抗失配性优于GIG-GLRT和MAP-GLRT检测器。综上所述，提出的检测器仅以牺牲匹配情况下较小的检测性能为代价，提升了出现失配信号时的抗失配性能。

图 6 检测器抗失配性能对比

下载: 全尺寸图片幻灯片

5. 结论

本文针对导向矢量失配情况下的目标检测问题，为了抑制失配信号，在 ${H_0}$ 假设中人为引入虚拟信号 $\beta {\boldsymbol{q}}$ ，并基于GLRT和MAP GLRT准则在GIG纹理背景下发展对散斑协方差矩阵 ${\boldsymbol{M}}$ 和导向矢量 ${\boldsymbol{p}}$ 具有CFAR特性的A-CGGIG和AM-CGGIG检测器。仿真和实测数据实验结果表明，提出的检测器兼顾匹配情况下的目标检测性能和失配情况下的抗失配性能。

图 1 失配示意图

下载: 全尺寸图片幻灯片

图 2 在匹配情况下提出的检测器与对比检测器的目标检测性能

下载: 全尺寸图片幻灯片

图 3 在失配情况下提出的检测器与对比检测器的目标检测性能

下载: 全尺寸图片幻灯片

图 4 检测器虚警概率曲线

下载: 全尺寸图片幻灯片

图 5 实测海杂波数据背景下检测器性能分析

下载: 全尺寸图片幻灯片

图 6 检测器抗失配性能对比

下载: 全尺寸图片幻灯片

导向矢量匹配时上述检测器检测概率达到 $90\%$ 时所需要的SCR(dB)

检测器类型 ANMF ABORT GIG-GLRT MAP-GLRT α-AMF A-CGGIG AM-GIGIG

信杂比 8.40 7.61 5.70 5.71 5.82 5.80 5.93

下载: 导出CSV

导向矢量匹配时上述检测器检测概率达到 $90\%$ 时所需要的SCR(dB)

检测器类型 ANMF ABORT GIG-GLRT MAP-GLRT α-AMF A-CGGIG AM-GIGIG

信杂比 12.88 12.16 11.34 11.50 11.43 11.37 11.53

下载: 导出CSV

参考文献(18)

[1]	许述文, 白晓惠, 郭子薰, 等. 海杂波背景下雷达目标特征检测方法的现状与展望[J]. 雷达学报, 2020, 9(4): 684–714. doi: 10.12000/JR20084. XU Shuwen, BAI Xiaohui, GUO Zixun, et al. Status and prospects of feature-based detection methods for floating targets on the sea surface[J]. Journal of Radars, 2020, 9(4): 684–714. doi: 10.12000/JR20084.
[2]	XIN Zhihui, LIAO Guisheng, YANG Zhiwei, et al. Analysis of distribution using graphical goodness of fit for airborne SAR sea-clutter data[J]. IEEE Transactions on Geoscience and Remote Sensing, 2017, 55(10): 5719–5728. doi: 10.1109/TGRS.2017.2712700.
[3]	LIANG Xiang, YU Han, ZOU Pengjia, et al. Multiscan recursive Bayesian parameter estimation of large-scene spatial-temporally varying generalized Pareto distribution model of sea clutter[J]. IEEE Transactions on Geoscience and Remote Sensing, 2022, 60: 5115416. doi: 10.1109/TGRS.2022.3191467.
[4]	水鹏朗, 田超, 封天. 逆高斯纹理复合高斯杂波对异常样本稳健的三分位点估计方法[J]. 电子与信息学报, 2023, 45(2): 542–549. doi: 10.11999/JEIT211483. SHUI Penglang, TIAN Chao, and FENG Tian. Outlier-robust tri-percentile parameter estimation method of compound-Gaussian clutter with inverse Gaussian textures[J]. Journal of Electronics & Information Technology, 2023, 45(2): 542–549. doi: 10.11999/JEIT211483.
[5]	XUE Jian, XU Shuwen, LIU Jun, et al. Model for non-Gaussian sea clutter amplitudes using generalized inverse Gaussian texture[J]. IEEE Geoscience and Remote Sensing Letters, 2019, 16(6): 892–896. doi: 10.1109/LGRS.2018.2886782.
[6]	CONTE E, LOPS M, and RICCI G. Asymptotically optimum radar detection in compound-Gaussian clutter[J]. IEEE Transactions on Aerospace and Electronic Systems, 1995, 31(2): 617–625. doi: 10.1109/7.381910.
[7]	DONG Yunhan. Optimal coherent radar detection in a K-distributed clutter environment[J]. IET Radar, Sonar & Navigation, 2012, 6(5): 283–292. doi: 10.1049/iet-rsn.2011.0273.
[8]	SANGSTON K J, GINI F, and GRECO M S. Coherent radar target detection in heavy-tailed compound-Gaussian clutter[J]. IEEE Transactions on Aerospace and Electronic Systems, 2012, 48(1): 64–77. doi: 10.1109/TAES.2012.6129621.
[9]	CHEN Sijia, KONG Lingjiang, and YANG Jianyu. Adaptive detection in compound-Gaussian clutter with inverse Gaussian texture[J]. Progress in Electromagnetics Research M, 2013, 28: 157–167. doi: 10.2528/PIERM12121209.
[10]	薛健. 复合高斯海杂波背景雷达目标检测算法[D]. [博士论文], 西安电子科技大学, 2020: 80–88. doi: 10.27389/d.cnki.gxadu.2020.000039. XUE Jian. Radar target detection methods in compound-Gaussian sea clutter[D]. [Ph. D. dissertation], Xidian University, 2020: 80–88. doi: 10.27389/d.cnki.gxadu.2020.000039.
[11]	XU Shuwen, WANG Zhexiang, BAI Xiaohui, et al. Optimum and near-optimum coherent CFAR detection of radar targets in compound-Gaussian clutter with generalized inverse Gaussian texture[J]. IEEE Transactions on Aerospace and Electronic Systems, 2022, 58(3): 1692–1706. doi: 10.1109/TAES.2021.3120045.
[12]	施赛楠, 水鹏朗, 刘明. 基于复合高斯杂波纹理结构的相干检测[J]. 电子与信息学报, 2016, 38(8): 1969–1976. doi: 10.11999/JEIT151194. SHI Sainan, SHUI Penglang, and LIU Ming. Coherent detection based on texture structure in compound-Gaussian clutter[J]. Journal of Electronics & Information Technology, 2016, 38(8): 1969–1976. doi: 10.11999/JEIT151194.
[13]	WANG Zhihang, HE Zishu, HE Qin, et al. Persymmetric range-spread targets detection in compound Gaussian sea clutter with inverse Gaussian texture[J]. IEEE Geoscience and Remote Sensing Letters, 2022, 19: 8018305. doi: 10.1109/LGRS.2021.3101369.
[14]	PULSONE N B and RADER C M. Adaptive beamformer orthogonal rejection test[J]. IEEE Transactions on Signal Processing, 2001, 49(3): 521–529. doi: 10.1109/78.905870.
[15]	WANG Zhihang, HE Zishu, HE Qin, et al. Adaptive CFAR detectors for mismatched signal in compound Gaussian sea clutter with inverse Gaussian texture[J]. IEEE Geoscience and Remote Sensing Letters, 2022, 19: 3502705. doi: 10.1109/LGRS.2020.3047390.
[16]	WANG Zeyu, GANG Li, and LI Ming. Adaptive detection of distributed target in the presence of signal mismatch in compound Gaussian clutter[J]. Digital Signal Processing, 2020, 102: 102755. doi: 10.1016/j.dsp.2020.102755.
[17]	王喆祥. 海杂波背景下雷达目标自适应检测算法研究[D]. [硕士论文], 西安电子科技大学, 2021: 12, 16. doi: 10.27389/d.cnki.gxadu.2021.001059. WANG Zhexiang. Research on adaptive detection algorithm of radar targets in sea clutter[D]. [Master dissertation], Xidian University, 2021: 12, 16. doi: 10.27389/d.cnki.gxadu.2021.001059.
[18]	PASCAL F, FORSTER P, OVARLEZ J P, et al. Performance analysis of covariance matrix estimates in impulsive noise[J]. IEEE Transactions on Signal Processing, 2008, 56(6): 2206–2217. doi: 10.1109/TSP.2007.914311.

施引文献

资源附件(0)

访问统计

图(6) / 表(2)

计量

文章访问数: 283
HTML全文浏览量: 166
PDF下载量: 34
被引次数: 0

1. 引言
2. 失配信号的目标检测模型
3. 基于GIG纹理结构的失配检测器
3.1 GLRT准则下的失配检测器
3.2 MAP GLRT准则下的失配检测器
3.3 失配检测器的CFAR特性分析
4. 失配检测器性能评估
4.1 仿真数据实验结果及性能分析
4.2 实测数据实验结果及性能分析
5. 结论

1. 引言
2. 失配信号的目标检测模型
3. 基于GIG纹理结构的失配检测器
3.1 GLRT准则下的失配检测器
3.2 MAP GLRT准则下的失配检测器
3.3 失配检测器的CFAR特性分析
4. 失配检测器性能评估
4.1 仿真数据实验结果及性能分析
4.2 实测数据实验结果及性能分析
5. 结论

参考文献(18)

施引文献

资源附件(0)

访问统计

姓名
邮箱
手机号码
标题
留言内容
验证码

留言板

广义逆高斯纹理海杂波背景下的自适应失配检测器

doi: 10.11999/JEIT231440

通讯作者:
范一飞　fanyifei888@163.com

计量

Adaptive Detectors for Mismatched Signal under Sea Clutter Background with Generalized Inverse Gaussian Texture

1. 引言

2. 失配信号的目标检测模型

3. 基于GIG纹理结构的失配检测器

3.1 GLRT准则下的失配检测器

3.2 MAP GLRT准则下的失配检测器

3.3 失配检测器的CFAR特性分析

4. 失配检测器性能评估

4.1 仿真数据实验结果及性能分析

4.2 实测数据实验结果及性能分析

5. 结论

计量

目录

1. 引言

2. 失配信号的目标检测模型

3. 基于GIG纹理结构的失配检测器

3.1 GLRT准则下的失配检测器

3.2 MAP GLRT准则下的失配检测器

3.3 失配检测器的CFAR特性分析

4. 失配检测器性能评估

4.1 仿真数据实验结果及性能分析

4.2 实测数据实验结果及性能分析

5. 结论

留言板

广义逆高斯纹理海杂波背景下的自适应失配检测器

doi: 10.11999/JEIT231440

通讯作者: 范一飞 fanyifei888@163.com

计量

出版历程

Adaptive Detectors for Mismatched Signal under Sea Clutter Background with Generalized Inverse Gaussian Texture

1. 引言

2. 失配信号的目标检测模型

3. 基于GIG纹理结构的失配检测器

3.1 GLRT准则下的失配检测器

3.2 MAP GLRT准则下的失配检测器

3.3 失配检测器的CFAR特性分析

4. 失配检测器性能评估

4.1 仿真数据实验结果及性能分析

4.2 实测数据实验结果及性能分析

5. 结论

计量

出版历程

目录

1. 引言

2. 失配信号的目标检测模型

3. 基于GIG纹理结构的失配检测器

3.1 GLRT准则下的失配检测器

3.2 MAP GLRT准则下的失配检测器

3.3 失配检测器的CFAR特性分析

4. 失配检测器性能评估

4.1 仿真数据实验结果及性能分析

4.2 实测数据实验结果及性能分析

5. 结论

通讯作者:
范一飞　fanyifei888@163.com