[发明专利]基于SAR‑DPCA微多普勒时‑空‑频三维分析方法

说明书

技术领域

本发明属于一种基于SAR-DPCA微多普勒的时-空-频三维分析的方法，属于SAR雷达微多普勒数据处理方法与应用技术领域，适用于SAR雷达微多普勒数据分析的理论方法和应用技术研究。

背景技术

DPCA即偏置相位中心天线技术是多通道SAR雷达抑制杂波的研一个典型的方法。它是用两个脉冲间通过阵列天线孔径相对于飞行方向的偏移来补偿雷达平台的移动。经典的DPCA方法采用两个相位中心，在相位中心距、脉冲重复频率和平台速度之间满足一定关系的基础上进行杂波抑制。正是由于此法工程上易用的优点，DPCA方法在多通道SAR雷达系统中得到了广泛地应用。

不同运动形态的动目标，其微多普勒调制信号具有独特性。介于此，对微多普勒效应的研究能够为我们提供关于雷达动目标更多的信息，为现有的雷达动目标识别方法提供一个新的途径，或成为现有一些成熟识别技术的有力补充。

时频分析又称时频分布，是一种描述信号的频率随时间变化的信号处理方法。它采用时间-频率联合来表征信号，通过在时频域分析信号，以实现全面观测信号特征的目的。时频分布主要分为线性时频表示和二次型时频表示。线性时频表示是从傅立叶变换转换来的，其特点是呈线性，主要有短时傅立叶变换，Gabor变换和小波分析。二次型变换主要是Wigner-Ville(WVD)方法。但无论是线性时频表示还是二次型时频表示，他们只能分析时间-频率这两个维度的变化，但在有些场景，比如当两个动目标在相同的距离维度，此时便无法用时频分析来描述微多普勒的信号特征。

基于此本发明，在时间-频率的基础上，引入另一个与距离相关的变量。使得原本在二维平面上的分析，拓展到三维平面。解决了不同运动点在相同距向离维度，不同方位向维度上的微多普勒特征分析。

发明内容

本发明的目的在于针对传统的SAR雷达微多普勒时频分析方法的不足，为SAR雷达微多普勒特征的研究提供一种基于SAR-DPCA微多普勒的时-空-频三维分析的方法。

本发明方法包括以下步骤：

步骤(1)获得SAR双通道距离向压缩数据域回波数据

将SAR双通道接收的雷达回波信号进行相干检波和距离压缩处理。处理后SAR双通道的基带信号为：和其中常数C₁和C₂是根据运动目标的后向散射系数和天线的方位向图所共同决定的。τ表示时间，s₁(τ)表示第一个通道的信号，s₂(τ)表示第二个通道的信号，λ表示雷达发射信号的波长，j表示虚部单位，R₀表示雷达发射通道相位中心与动目标的距离，R₁、R₂分别表示雷达两个子孔径与动目标的距离。

步骤(2)SAR双通道相位补偿与静杂波抑制

将获得的SAR双通道基带信号分别进行相位补偿，去除两个子孔径因相位中心的不同而产生的多普勒中心频率的偏差和载机上雷达运动所产生的多普勒调频项。用于相位补偿的参考函数分别为:其中P_c1(τ)表示补偿s₁(τ)的参考函数，P_c2(τ)表示补偿s₂(τ)的参考函数，j表示虚部单位，λ表示雷达发射信号的波长，v为载机速度，B为两个子孔径之间的距离，τ表示时间，prf为雷达脉冲重复频率，且它们之间满足补偿方式为，基带信号与相应的参考函数相乘。得到处理后的两路信号s'₁(τ)和s'₂(τ)。然后用s'₂(τ)减去s'₁(τ)得s₃(τ)。

步骤(3)时-空-频三维分析

(3.1)提取动目标区域信号切片

对复数s₃(τ)取其虚部的绝对值，得s₄(τ)。然后求得s₄(τ)的最大值max。然后提取s₄(τ)中大于max*0.05的部分，得s₅(τ)为m*n维矩阵。

(3.2)时-空-频三维分析，采用短时傅里叶变换对s₅(τ)进行时频分析具体步骤

步骤1：抽取s₅(τ)的第i行的数据，即得a_i(1,n)。