基于稳健稀疏分数阶傅立叶变换的雷达机动目标快速检测方法技术

本发明专利技术涉及基于稳健稀疏分数阶傅立叶（RSFRFT）变换的雷达机动目标快速检测方法，属于雷达信号处理和检测技术领域。其步骤包括：对脉冲压缩后的慢时间维雷达复数据进行时域Chirp乘法运算、频谱重排、窗函数滤波、降采样快速傅里叶变换；在频域降采样形成的各个频率通道中，进行过门限处理；重构待检测距离单元信号多普勒频率；构造检测器，进行目标检测，频域Chirp乘法运算，得到RSFRFT谱，并估计目标运动参数。本发明专利技术只需对少量目标多普勒频点进行搜索，降低计算复杂度，并且通过两次自适应门限，无需预设信号稀疏度，剔除杂波虚警，适合复杂杂波背景下的机动目标检测。

Fast Detection Method of Radar Maneuvering Target Based on Robust Sparse Fractional Fourier Transform

【技术实现步骤摘要】
基于稳健稀疏分数阶傅立叶变换的雷达机动目标快速检测方法
本专利技术属于雷达信号处理和检测
，更具体地，本专利技术涉及一种基于稳健稀疏分数阶傅立叶变换的雷达机动目标快速检测方法，可用于强杂波背景下机动目标的雷达检测处理。
技术介绍
低可观测机动目标的快速、有效检测已成为雷达
的世界性难题。受复杂环境背景和目标复杂运动特性的影响，机动目标回波信杂噪比(Signal-to-clutterandnoiseratio,SCNR)低，且回波多普勒呈现时变特性，增加了雷达探测的难度。目标多普勒频率与目标速度近似成正比，当目标做机动时，回波具有时变特性，例如，匀加速运动目标，回波可建模为线性调频(LinearFrequencyModulation,LFM)信号。运动状态较为复杂的目标在一段观测时间内，可用LFM信号作为其一阶近似。随着相控阵、泛探雷达、多输入多输出(Multiple-inputmultiple-output,MIMO)雷达等雷达新体制的发展，对目标的观测时间大大延长，提高杂波背景下机动目标的精细化处理能力。但这种凝视观测以及泛探工作模式，使得目标回波数据量增大，并且系统采样频率的提高进一步增加了数据量，从而对算法的运算效率和系统实时性提出了更高的要求。因此，研究能适应杂波环境的、高效可靠、适用于大数据量实时处理的机动目标信号分析方法对于提升雷达动目标探测能力具有重要意义。近年来，稀疏信号处理技术的发展为雷达探测技术提供了新的研究思路。由于动目标回波信号在某个域中具有一定的稀疏特性，可将动目标检测问题，转换为稀疏域中的稀疏求解和检测问题，则可通过少量的观测样本，通过求解最优化问题，在稀疏域中实现对该信号的高分辨率表示，比较有代表性的是稀疏傅里叶变换(SparseFourierTransform,SFT)和稀疏分数阶傅立叶变换(SparsefractionalFouriertransform,SFRFT)。然而，以下两方面因素导致SFT和SFRFT的稳健性下降：一方面，SFRFT需要对信号的稀疏度K进行预设，而在实际应用中，信号的稀疏度往往是未知的或者可能发生改变的；另一方面，SFRFT在降采样FFT后仅结合稀疏度和循环过程中频点出现的概率对大值系数进行估计，这在较低SCNR的情况下难以保证重构信号的可靠性。而雷达目标检测通常会面临杂波背景，频域聚集形成的强杂波点增加了回波数据的稀疏度，大大降低了算法对目标信号的处理性能。另外，SFRFT算法在重构目标信号多普勒的同时判决目标的有无，此检测判决是在平坦窗函数滤波降低了SCNR的情况下做出的，不利于杂波背景下的目标检测问题。因此，基于SFRFT的动目标检测方法难以满足复杂环境中的雷达机动目标检测需求。目前该领域存在的问题主要有以下三个方面：一是，经典的动目标检测(MovingTargetDetection,MTD)方法仅适用于匀速运动目标，对于非平稳时变信号，难以实现可靠的检测；二是，FRFT方法，需逐个多普勒通道进行旋转角匹配搜索，运算量大；三是，SFRFT方法需要预设稀疏度，而且在低SCNR情况下，信号重构可靠性差，检测性能严重下降。为此，本专利提出了稳健SFRFT(RobustFRFT,RSFRFT)并用于雷达目标检测。
技术实现思路
要解决的技术问题包括：(1)经典的MTD方法仅适用于匀速运动目标，对于非平稳时变信号，难以实现可靠的检测；(2)FRFT动目标检测方法，需逐个多普勒通道进行旋转角匹配搜索，运算量大；(3)SFRFT动目标检测方法，需要预设稀疏度，而且在低SCNR情况下，信号重构可靠性差，检测性能严重下降。本专利技术所述的基于稳健稀疏分数阶傅立叶变换的雷达机动目标快速检测方法，其特征在于包括以下技术措施：步骤一、对脉冲压缩后的慢时间维雷达复数据进行时域Chirp乘法运算、频谱重排、窗函数滤波、降采样快速傅里叶变换；步骤二、在频域降采样形成的各个频率通道中，进行过门限处理；步骤三、重构待检测距离单元信号多普勒频率；步骤四、构造检测器，进行目标检测，频域Chirp乘法运算，得到稳健稀疏分数阶傅立叶变换谱。对比现有技术，本技术方案所述的基于稳健稀疏分数阶傅立叶变换的雷达机动目标快速检测方法，有益效果在于：(1)相比于基于FRFT的常规检测方法(需要逐个多普勒频率进行匹配搜索进而构建检测器)，RSFRFT算法只需对少量目标多普勒频点进行搜索，因而能够较大程度的降低计算复杂度；(2)相比于SFRFT方法，RSFRFT算法预设信号稀疏度，既能满足实际工程对运算量的需求，又适合复杂杂波背景下的机动目标检测。附图说明图1是基于稳健稀疏分数阶傅立叶变换的雷达机动目标快速检测方法的实施流程图。具体实施方式以下结合说明书附图对本专利技术作进一步详细描述。参照说明书附图，本专利技术的处理流程分以下步骤：1)沿距离向对脉压后的雷达数据进行时域Chirp乘法运算、频谱重排、窗函数滤波、降采样快速傅里叶变换。时域Chirp乘法运算：对脉冲压缩后的慢时间维雷达复数据x(n)与Chirp1信号相乘，式中，为Chirp1信号，T为脉冲重复周期，N为脉冲个数，α表示旋转角度；频谱重排：为使大值频点尽可能的均匀随机分布，通过对时域信号进行操作以实现频谱重排，定义频谱重排方式Pσ，以及重排后的时域序列Pσ(n)＝r[(σn)modN],n∈[1,N]式中，σ为[1,N]中随机选取的奇数，mod为取模运算；窗函数滤波：定义平坦窗函数g(n)，其频谱G(m)满足式中，ε′和ε分别为通带截断因子和阻带截断因子，δ为震荡波纹，则窗函数滤波后的信号为式中，supp表示支撑，ω为窗函数长度；降采样快速傅里叶变换：由傅里叶变换性质可知，通过时域混叠可以实现频域的降采样，降采样快速傅里叶变换后的信号为式中，表示向下取整运算，FFT{}表示快速傅里叶变换运算，B表示频域降采样后的数据点数。2)在频域降采样形成的各个频率通道中，进行过门限处理。在降采样快速傅里叶变换后增加第一门限η1来估计信号稀疏度及频点，无需对稀疏度进行预设，将Z(m)中幅值超过门限η1的频点对应的坐标归入集合J中，J＝{m∈[1,B]|Z(m)≥η1}其中，门限η1通过检测技术确定，如针对时间平稳性较低、而空间平稳性较高的杂波背景，可根据背景的均匀程度适当选择均值类、有序统计量类以及自适应恒虚警等空域恒虚警检测技术。3)重构待检测距离单元信号多普勒频率。定义哈希函数通过哈希逆映射得到J中频点在信号r(n)的频谱序列中的对应坐标并保存到集合U中，U＝{m∈[1,N]|hσ(m)∈J}，哈希逆映射需要进行M次定位循环，若目标存在，则在M次定位循环中，目标对应的多普勒频率应该具有较高的“出现次数”，因此，在重构过程中设置“出现次数阈值”γ，将“出现次数”超过该阈值的频点对应的多普勒频率称为待检测距离单元信号多普勒频率。4)构造检测器，进行目标检测，频域Chirp乘法运算，得到稳健稀疏分数阶傅立叶变换谱。针对重构后得到的待检测距离单元信号多普勒频率，结合待检测距离单元和参考距离单元中的原始观测数据构造检测器进行检测判决，设重构得到的待检测距离单元信号多普勒频率为fα,d0，则目标多普勒矢量为观测矢量为通过结合检测单元和参考单元的原始观本文档来自技高网...

【技术保护点】
1.基于稳健稀疏分数阶傅立叶变换的雷达机动目标快速检测方法，其特征在于，包括以下步骤：步骤一、对脉冲压缩后的慢时间维雷达复数据进行时域Chirp乘法运算、频谱重排、窗函数滤波、降采样快速傅里叶变换；步骤二、在频域降采样形成的各个频率通道中，进行过门限处理；步骤三、重构待检测距离单元信号多普勒频率；步骤四、构造检测器，进行目标检测，频域Chirp乘法运算，得到稳健稀疏分数阶傅立叶变换谱。

【技术特征摘要】
1.基于稳健稀疏分数阶傅立叶变换的雷达机动目标快速检测方法，其特征在于，包括以下步骤：步骤一、对脉冲压缩后的慢时间维雷达复数据进行时域Chirp乘法运算、频谱重排、窗函数滤波、降采样快速傅里叶变换；步骤二、在频域降采样形成的各个频率通道中，进行过门限处理；步骤三、重构待检测距离单元信号多普勒频率；步骤四、构造检测器，进行目标检测，频域Chirp乘法运算，得到稳健稀疏分数阶傅立叶变换谱。2.根据权利要求1所述的基于稳健稀疏分数阶傅立叶变换的雷达机动目标快速检测方法，其特征在于，步骤一所述的时域Chirp乘法运算、频谱重排计算方法为：对脉冲压缩后的慢时间维雷达复数据x(n)与Chirp1信号相乘，式中，为Chirp1信号，T为脉冲重复周期，N为脉冲个数，α表示旋转角度；定义频谱重排方式Pσ，以及重排后的时域序列Pσ(n)＝r[(σn)modN],n∈[1,N]式中，σ为[1,N]中随机选取的奇数，mod为取模运算。3.根据权利要求2所述的基于稳健稀疏分数阶傅立叶变换的雷达机动目标快速检测方法，其特征在于，步骤一所述的窗函数滤波计算方法为：定义平坦窗函数g(n)，其频谱G(m)满足式中，ε′和ε分别为通带截断因子和阻带截断因子，δ为震荡波纹，则窗函数滤波后的信号为y(n)＝g(n)·Pσ(n),n∈[1,N],式中，supp表示支撑，ω为窗函数长度。4.根据权利要求3所述的基于稳健稀疏分数阶傅立叶变换的雷达机动目标快速检测方法，其特征在于，步骤一所述的降采样快速傅里叶变换计算方法为：降采样快速傅里叶变换后的信号为式中，表示向下取整运算，FFT{}表示快速傅里叶变换运算，B表示频域降采样后的数据点数。5.根据权利要求4所述的基于稳健稀疏分数阶傅立叶变换的雷达机动目标快速检测方法，其特征在于，步骤二所述的过门限处理方法为：J＝{m∈[1,B]|Z(m)≥η1}并将Z(m)中幅值超过门限η1的频点对...

【专利技术属性】
技术研发人员：陈小龙，于晓涵，黄勇，薛永华，关键，
申请(专利权)人：中国人民解放军海军航空大学，
类型：发明
国别省市：山东,37