[发明专利]一种基于奇异谱分解的海杂波压制及目标检测方法

说明书

一种基于奇异谱分解(Singular Spectrum Analysis,SSA)的海杂波压制及目标检测方法，首先计算回波数据的自相关函数，以自相关函数的第一过零点位置作为重构雷达回波数据轨迹矩阵的维数，构建回波数据的Toeplitz矩阵、求得主分量(Principal Component Analysis,PCA)数量，依据PCA数量判断回波中是否包含目标；而后对雷达回波进行奇异谱分解，重构主分量和次分量(Minor Component Analysis,MCA)对应的时间序列，当回波中包含目标时，取最大主分量对应的时间序列、舍去其他分量对应的时间序列；反之，则选择最小次分量对应的时间序列而舍去其他分量对应的时间序列，这一处理过程可被视为不依赖回波统计学模型的自适应滤波，经处理后的雷达回波数据中海杂波强度被有效压制，有利于提高后续弱小目标的检测概率。

技术领域

本发明属于目标检测技术领域，涉及海杂波背景下的目标检测，为一种基于奇异谱分解的海杂波压制及目标检测方法。

背景技术

海杂波是指对海观测雷达波束照射到海面时所产生的回波，海杂波的特性与雷达波束覆盖的区域面积、海浪高度(即海况)、雷达工作主频及带宽等因素密切相关，通常海杂波反映了由重力引起的排浪的运动特性和海面张力引起的管状细碎波的微动特性，即海杂波包含2个主分量，海杂波的幅度、带宽等参数随海况的变化呈现出非平稳、非线性特性，如附图1中(a)和(b)所示，为利用加拿大McMaster大学IPIX雷达测量的海杂波数据(http://soma.ece.mcmaster.ca/ipix/dartmouth/datasets.html)的时域幅度图，存在明显的尖峰，且尖峰的强度与数量随海况的增高而增加，附图2为上述海杂波数据的频谱图，如其中(a)和(b)所示，可见高海况时由于海面张力波引起的微动特性加剧而导致其带宽显著展宽。对海杂波背景下弱小目标检测方法的研究一直是雷达信号处理领域的热点，在民用领域内被用于检测海面漂浮物(浮冰、小型船舶和飞行器残骸等)，为民用船舶导航和海上搜救等提供技术支持。在军事领域中，该方法被用于探测海上具有隐身特性的舰船、潜艇潜望镜和贴近海面飞行的飞行器等，传统的方法是将海杂波视为一种复杂的平稳随机过程，通过大量的观测数据建立统计学意义上的概率模型，如Weibull分布，log-normal分布、复合K分布模型和Pareto分布等，进而利用成熟的检测方法实现目标检测。然而，对海杂波平稳性的假设是建立在短观测时间基础上的，为提高弱小目标的检测概率，雷达不得不增加对目标区域的观测时间，以提高经过相干积累后目标回波能量，但是随观测时间的增长海杂波将不再是平稳随机过程，如附图3中(a)和(b)所示，为南非CSIR TFC15_023.01海杂波数据集的Pareto分布参数，其形状参数k和尺度参数sigma均随距离单元的不同而随机波动，若采用固定模型参数必然会造成较大的测量误差。

综上，对海观测雷达为探测海上弱小目标，通常以增加对海面目标区域的观测时间来增强目标的回波能量，然而在大相干处理时间内所接收的海面回波(海杂波)强度也相应增强，且呈现出非线性、非高斯和非平稳特性，传统的海杂波背景下的弱小目标检测方法均是将海杂波视为服从某种概率密度分布，以有无目标时的概率密度分布特征之间的差异检测目标，传统方法忽略了海杂波的非平稳特性，在高海况时虚警概率过高，不适用于大相干处理时间的目标检测。

为解决上述问题，雷达工程领域内的专家学者进行了大量的研究工作，已有一定数量的研究成果。经典的解决方案之一是对雷达回波的统计学特性进行动态跟踪，按特性相近的原则对其进行分段处理，在工程层面上解决了部分问题，但随着海况的提高，雷达回波数据的分段越来越短，最终与短时观测的情形无本质的差异，致使对弱小目标探测时仍存在较高的虚警概率、较低的发现概率，因此，基于统计学特性的信号处理方法在高海况条件下已不满足工程的需要，必须寻找不依赖统计学模型的海杂波抑制方法，进而解决这一难题。

发明内容