本发明隶属于雷达信号处理领域，具体涉及一种基于线性融合的雷达目标自适应检测方法。
背景技术：
雷达目标自适应检测主要涉及两个方面：未知杂波谱特性的准确估计和自适应检测器结构的合理构建；未知杂波协方差矩阵估计是杂波谱特性估计中的重要内容，其估计结果的好坏直接影响最终的检测性能，实际中可利用与待检测距离单元邻近的不包含目标的参考数据，对未知的杂波谱特性进行估计。
对低分辨率雷达场景，杂波分布满足中心极限定理，不同距离单元杂波统计特性可用高斯分布进行建模。但随着杂波环境的复杂化、雷达分辨率的提高或在低掠地角条件下，背景杂波呈现出较强的非高斯特性，一般可采用复合高斯分布进行描述和建模；根据中心极限定理的局部有效性，复合高斯分布可表示为时空慢变的非负纹理分量与时空快变的复高斯散斑分量之乘积。针对经典的复高斯背景下的雷达点目标检测问题，广义似然比检验检测器通过对所有未知参数进行联合估计，充分利用待检测距离单元数据和参考数据中的杂波信息，在高斯环境中能获得最优的检测性能，且对高斯杂波具有恒虚警率(CFAR)特性；复高斯背景下未知杂波协方差矩阵的最大似然估计可表示为采样协方差矩阵形式，与广义似然比检验检测器相比，基于采样协方差矩阵的自适应匹配滤波器能以较小的运算量在高斯环境中获得次优的检测性能，且对高斯杂波具有CFAR特性。以上两种点目标检测器在高斯背景下各有所长，广义似然比检验检测器能获得最优检测性能，但运算量略高，而自适应匹配滤波器在降低运算量的同时仍能保持次优的检测性能，但二者的检测性能和CFAR特性均随着杂波非高斯特性的增强而变差。在杂波服从复合高斯分布的非高斯背景下，未知杂波协方差矩阵的最大似然估计没有闭型表达式，且涉及到超越方程的求解，归一化采样协方差矩阵等均为次优估计方法；以两步法广义似然比检验为基础，基于归一化采样协方差矩阵的自适应归一化匹配滤波器在非高斯环境中具有较好的检测性能，且对复合高斯杂波表现出一定的CFAR鲁棒性，但其检测性能和CFAR特性在高斯杂波背景下有所退化。虽然针对高斯和复合高斯杂波背景存在最优或次优的杂波协方差矩阵估计方法和点目标检测器结构，但实际环境中，杂波的非高斯特性往往在时间和空间上随环境的改变而变化，上述特定杂波背景下的最优或次优杂波协方差矩阵估计方法和点目标检测器难以适应杂波环境的快速变化，导致相应的检测性能和CFAR特性恶化。
针对实际环境中杂波非高斯程度的时空易变特点，如何综合利用特定杂波背景下最优或次优杂波协方差矩阵估计方法和点目标检测器，合理利用杂波非高斯特征信息，通过简洁有效的融合规则设计合适的点目标自适应检测器结构，在增强检测器CFAR鲁棒性的前提下，进一步提高未知杂波环境下的雷达目标检测性能，是提升复杂电磁环境下雷达探测能力的关键，也是目前急需解决的难题。
技术实现要素：
针对实际环境中杂波非高斯程度的时空易变特点，如何综合利用特定杂波背景下最优或次优杂波协方差矩阵估计方法，充分发挥自适应匹配滤波器和广义似然比检验检测器对高斯杂波环境的优越性能，以及自适应归一化匹配滤波器对非高斯杂波环境的鲁棒性，合理利用杂波非高斯特征信息，通过简洁有效的融合规则，设计合适的点目标自适应检测器结构，兼顾实际检测环境对检测性能和算法运算量的需求，增强检测器的CFAR鲁棒性和对杂波环境的自适应性，提高未知杂波环境下的雷达目标检测性能。
本发明所述基于线性融合的雷达目标自适应检测方法包括以下技术措施：
步骤1获取单个待检测距离单元雷达回波向量x0作为待检测数据向量，以单个待检测距离单元为中心，在其前后分别连续取一定数量距离单元雷达回波观测数据，构成R个参考数据向量xm,m＝1,2,…,R，其中，x0和xm,m＝1,2,…,R均为N×1维的向量，N表示雷达接收阵元数与相干处理脉冲数的乘积；
步骤2利用参考数据向量xm,m＝1,2,…,R，分别计算未知协方差矩阵的采样协方差矩阵估计值和近似最大似然估计值
其中，未知协方差矩阵的采样协方差矩阵估计值表示为
上式中，上标H表示矩阵的共轭转置运算；
未知协方差矩阵的近似最大似然估计值通过如下的迭代过程实现
其中，上标-1表示矩阵的求逆运算，k＝0,1,...,K-1，K表示最终的迭代次数，表示第k次迭代的估计矩阵；为增强估计过程的普适性，式(2)的初始化矩阵采用单位阵，而最终的迭代结果即为近似最大似然估计值
步骤3根据实际检测环境对检测性能和算法运算量的需求，基于采样协方差矩阵估计值计算回波向量x0的均匀背景检测器检测统计量λh，基于近似最大似然估计值计算回波向量x0的自适应归一化匹配滤波器检测统计量λANMF；根据实际杂波的非高斯程度设定比例因子α，基于均匀背景检测器和自适应归一化匹配滤波器的线性融合，构建回波向量x0的自适应检测统计量λα；
其中，基于采样协方差矩阵估计值回波向量x0的均匀背景检测器的形式采用自适应匹配滤波器或广义似然比检验检测器；
自适应匹配滤波器检测统计量λAMF表示为
上式中，p表示已知的空时导向矢量，是一个N×1维的单位向量，根据雷达工作参数确定；
广义似然比检验检测器检测统计量λGLRT表示为
基于近似最大似然估计值回波向量x0的自适应归一化匹配滤波器检测统计量λANMF表示为
基于λh和λANMF的线性融合，构建回波向量x0的自适应检测统计量λα，即
λα＝αλh+(1-α)λANMF (6)
上式中，比例因子α满足0≤α≤1，且根据实际杂波的非高斯程度设定，使得检测器具有很强的泛化能力，其中杂波的非高斯程度越大，则α的取值越小；λh根据实际需要取λGLRT或λAMF，若检测器对检测性能要求较高则取λh＝λGLRT，若检测器对降低运算量要求较高则取λh＝λAMF；
步骤4为保持检测方法的恒虚警率特性，根据预设的虚警概率设置检测门限T；将回波向量x0的自适应检测统计量λα与门限T进行比较，若λα≥T，则判定待检测距离单元存在目标，回波向量x0不作为其他距离单元的参考数据；反之若λα<T，则判定待检测距离单元不存在目标，回波向量x0作为后续其他距离单元的参考数据。
与背景技术相比，本发明的有益效果是：1)合理利用杂波非高斯特征信息，能根据实际环境中杂波非高斯程度自适应设定线性融合中的比例因子；2)通过综合利用特定杂波背景下最优或次优杂波协方差矩阵估计方法，能将高斯背景下自适应匹配滤波器或广义似然比检验检测器与复合高斯背景下自适应归一化匹配滤波器进行简洁高效的线性融合，充分发挥自适应匹配滤波器和广义似然比检验检测器对高斯杂波环境的优越性能，以及自适应归一化匹配滤波器对非高斯杂波环境的鲁棒性；3)所设计的点目标检测器结构简洁，具有很强的泛化能力，能同时兼顾检测性能、算法运算量等多方面需求，能增强检测器的CFAR鲁棒性和对杂波环境的自适应性，提高未知杂波环境下的雷达目标检测的综合性能。
附图说明
图1是本发明所提出的基于线性融合的雷达目标自适应检测方法的功能模块图。图1中，1.采样协方差矩阵计算模块，2.近似最大似然估计矩阵计算模块，3.均匀背景检测器构建模块，4.自适应归一化匹配滤波器构建模块，5.线性融合计算模块，6.检测判决模块。
具体实施方式
下面结合附图对本发明作进一步描述。本发明实施例用来解释说明本发明，而不是对本发明进行限制，在本发明的精神和权利要求的保护范围内，对本发明做出的任何修改和改变，都落入本发明的保护范围。
参照说明书附图1，本发明的具体实施方式分为以下几个步骤：
步骤1待检测距离单元为中心，在其前后分别连续取一定数量距离单元雷达回波观测数据，构成R个参考数据向量xm,m＝1,2,…,R，其中，xm,m＝1,2,…,R均为N×1维的向量，N表示雷达接收阵元数与相干处理脉冲数的乘积；并将参考数据向量送至采样协方差矩阵计算模块(1)和近似最大似然估计矩阵计算模块(2)；
步骤2利用参考数据向量xm,m＝1,2,…,R，在采样协方差矩阵计算模块(1)中，根据式(1)计算未知协方差矩阵的采样协方差矩阵估计值并将送至均匀背景检测器构建模块(3)；在近似最大似然估计矩阵计算模块(2)中，根据式(2)计算未知协方差矩阵的近似最大似然估计值并将送至自适应归一化匹配滤波器构建模块(4)；
值得注意的是，当K＝1，根据式(2)计算的近似最大似然估计值表示为
由式(7)可知，由于式(2)的初始化矩阵采用单位阵，经过一次迭代后获得的近似最大似然估计值恰好为归一化采样协方差矩阵，说明归一化采样协方差矩阵为近似最大似然估计值在K＝1时的特例，显示了式(2)所示的近似最大似然估计过程具有很强的普适性；对于非高斯杂波环境，当K≥3时，近似最大似然估计能获得较好的估计精度，综合考虑算法运算量和检测器设计对协方差矩阵估计精度的要求，取K＝3的迭代结果作为近似最大似然估计值；
步骤3获取单个待检测距离单元的N×1维雷达回波向量x0作为待检测数据向量，并将x0送至均匀背景检测器构建模块(3)和自适应归一化匹配滤波器构建模块(4)；在均匀背景检测器构建模块(3)中，根据实际检测环境对检测性能和算法运算量的需求，基于计算回波向量x0的均匀背景检测器检测统计量λh，若检测器对检测性能要求较高则利用式(4)取λh＝λGLRT，若检测器对降低运算量要求较高则利用式(3)取λh＝λAMF；在自适应归一化匹配滤波器构建模块(4)中，根据式(5)利用计算回波向量x0的自适应归一化匹配滤波器检测统计量λANMF，并将λh和λANMF送至线性融合计算模块(5)；在线性融合计算模块(5)中，根据实际杂波的非高斯程度设定比例因子α，杂波的非高斯程度越大，则α的取值越小，基于均匀背景检测器和自适应归一化匹配滤波器的线性融合，根据式(6)构建回波向量x0的自适应检测统计量λα，将λα送至检测判决模块(6)；
值得注意的是，当背景杂波分布的高斯性较强时，α的取值较大，在极端的高斯背景下，取α＝1，此时由式(6)可知λα＝λh，说明基于采样协方差矩阵的自适应匹配滤波器或广义似然比检验检测器为本发明所提出的检测方法在α＝1时的特例，且本发明所提出的检测方法可根据实际检测环境对检测性能和算法运算量的需求自动选择两种检测器中的一种；而当背景杂波分布的非高斯性较强时，α的取值较小，在极端的复合高斯分布的非高斯背景下，取α＝0，此时由式(6)可知λα＝λANMF，说明自适应归一化匹配滤波器为本发明所提出的检测方法在α＝0时的特例；综上所述，本发明所提出的检测方法能合理利用杂波非高斯特征信息，根据实际杂波分布的非高斯程度自适应调整参数，使检测器结构能适应于实际杂波环境的变化，具有很强的泛化能力，通过将高斯和复合高斯背景下最优或次优杂波协方差矩阵估计方法和检测器进行匹配利用，充分发挥广义似然比检验检测器在高斯杂波背景下的最优检测性能、自适应匹配滤波器兼顾运算量和检测性能需求的优异综合性能、以及自适应归一化匹配滤波器在非高斯杂波背景下的CFAR鲁棒性和检测性能优势，且能涵盖现有三种检测器结构；
步骤4在检测判决模块(6)中进行检测判决并输出检测结果，为保持检测方法的恒虚警率特性，根据预设的虚警概率设置检测门限T；将单个待检测距离单元对应的回波向量x0的自适应检测统计量λα与门限T进行比较，若λα≥T，则判定待检测距离单元存在点目标，回波向量x0不能作为其他距离单元的参考数据；反之若λα<T，则判定待检测距离单元不存在点目标，回波向量x0可作为后续其他距离单元的参考数据。