一种鲁棒的压缩感知窄带自旋目标成像方法

摘要

本发明提供了一种鲁棒的压缩感知窄带自旋目标成像方法，对回波数据沿雷达视线方向进行傅里叶变换，使用时频谱自相关函数估计目标的旋转角速度和自旋角速度，将二维成像区域离散化，利用正交匹配追踪算法求解，通过非线性最小二乘方法提得到最后的目标像。由于采用将网格误差加入优化模型中，利用正交匹配追踪算法在不增加计算复杂度的情况下进行成像模型的优化重建，并通过非线性最小二乘方法提高散射点重建的准确性，相比传统压缩感知成像方法，有效解决了因网格误差造成的成像性能下降的问题，改善了目标散射点提取的准确性和鲁棒性。

主权项

一种鲁棒的压缩感知窄带自旋目标成像方法，其特征在于包括下述步骤：步骤1：对不同方位时刻下获取的回波数据，沿雷达视线方向进行傅里叶变换实现距离压缩；步骤2：在对目标进行包络对齐后，由于雷达发射的窄带信号，目标会被压缩到一个距离单元，因而信号能量最大的距离单元即为目标所在的距离单元，并将该距离单元回波数据记为s，其为N×1列向量；步骤3：使用时频谱自相关函数估计目标的旋转角速度w_s：由于自旋目标多普勒分量周期性分布的特点，对回波数据s进行时频变换，利用所得时频分布幅度的自相关函数的局部峰值估计目标的自旋角速度w_s；步骤4：将二维成像区域离散化，其横向x轴和纵向y轴对应的网格数分别记为N_x和N_y，网格单位尺寸为δr，并将该场景按行或列用一维列向量表示为其中σ_k表示第k个散射点的后向散射系数；步骤5：构建目标字典F和网格误差字典G、H，维度为N×N_xN_y，并定义辅助目标向量I′和I″，其维度同I一致，则目标信号s写成下面形式：s＝FI+GI′+HI″+E其中字典F、G和H的每一元素为：$F_{\mathrm{ni}}=-j\frac{4\pi }{\lambda }[x_i\cos {\theta }_n+y_i\sin {\theta }_n]$$G_{\mathrm{ni}}=F_{\mathrm{ni}}\cos {\theta }_n{\alpha }_1^{-1}$$H_{\mathrm{ni}}=F_{\mathrm{ni}}\sin {\theta }_n{\alpha }_2^{-1}$I′_i＝‑j4πα₁Δx_iI_i/λI″_i＝‑j4πα₂Δy_iI_i/λ${\alpha }_1={\left(\frac{1}{N}\underset{n}{\Sigma }{\cos }^2{\theta }_n\right)}^{1/2}$${\alpha }_2={\left(\frac{1}{N}\underset{n}{\Sigma }{\sin }^2{\theta }_n\right)}^{1/2},$(n＝1,2,…,N；i＝1,2,…,N_xN_y)其中F_ni为字典F的第n行第i列元素，G_ni和H_ni的定义同F_ni一致，I_i、I′_i和I″_i分别为I、I′和I″的第i个元素，E为信号噪声，其维度同s一致，λ为雷达波长，(x_i,y_i)为离散化场景中第i个网格位置，(Δx_i,Δy_i)为网格误差，θ_n＝w_st_n为目标自旋的第n个角度，t_n为方位慢时间的第n个时刻，α₁和α₂为辅助系数，代表虚数的单位，N为目标信号s的采样数；步骤6：利用改进的正交匹配追踪(OMP)算法求解步骤5中的方程，具体实施如下：步骤6‑1：初始化目标向量I＝I′＝I″＝0，剩余信号r＝s，支撑区域步骤6‑2：用剩余信号r分别向字典F、G和H中的每一列作投影，计算内积，并找到最大内积对应的向量位置i_max：$i_{\max }=\arg {\max }_i\left[\right|F_i^{*}r|+|G_i^{*}r|{(4{\mathrm{\pi \alpha }}_1\mathrm{\delta r}/\lambda )}^{-1}+|H_i^{*}r\left|{(4{\mathrm{\pi \alpha }}_2\mathrm{\delta r}/\lambda )}^{-1}\right]$其中F_i、G_i和H_i代表矩阵中的第i列，而和为F_i、G_i和H_i的共轭转置，α₁和α₂为步骤5中定义的辅助系数，δr为离散化场景网格大小，λ为雷达波长，并更新目标支撑区域Λ，即$\Lambda =\Lambda \cup \left\{i_{\max }\right\};$步骤6‑3：分别将最大内积对应的向量和记录在临时矩阵F′、G′和H′中，同时将对应向量从字典F、G和H中剔除；步骤6‑4：记录投影系数(I,I′,I″)＝arg min||F′I+G′I′+H′I″‑s||₂，并更新剩余信号为r＝s‑F′I‑G′I′‑H′I″；步骤6‑5：重复步骤6‑2～6‑4，直到剩余信号r的能量低于预设门限为止；步骤7：根据重构的I、I′和I″，可以恢复网格误差步骤8：利用改进OMP算法的重建结果，通过非线性最小二乘方法提高散射点重建的准确性，即在目标支撑区域Λ中求解下面的优化问题：$\arg \min \underset{n}{\Sigma }{|s_n-\underset{\Lambda }{\Sigma }{\sigma }_{\Lambda }\exp (-j\frac{4\pi }{\lambda }[(x_{\Lambda }+{\mathrm{\Delta x}}_{\Lambda })\cos {\theta }_n+(y_{\Lambda }+{\mathrm{\Delta y}}_{\Lambda })\sin {\theta }_n\left]\right)|}^{2}n=1,...,N$其中s_n为目标信号s的第n个采样值，x_Λ代表目标支撑区域Λ中的横向网格位置，Δx_Λ为其对应位置的横向网格误差，y_Λ代表目标支撑区域Λ中的纵向网格位置，Δy_Λ为其对应位置的纵向网格误差，步骤8中的优化问题可以采用现成的高斯牛顿方法、梯度方法等进行求解，输出结果σ_Λ即为最后的目标像。