CN106199524A - 基于基追踪去噪的远场宽带rcs数据采集与压缩方法 - Google Patents

Abstract

本发明提供了一种基于基追踪去噪的远场宽带RCS数据采集与压缩方法，通过在方位向与高程向上的随机降采样，降低了RCS数据的实际测量总量，能够在保证不过多损失雷达目标电磁散射特性信息完整度的情况下，缩短雷达系统对观测目标RCS数据的测量时间，从而提高整个测量系统的工作效率；根据基追踪去噪模型，利用结合了远场宽带RCS数据处理算子的稀疏重构算法，来实现对回波域RCS数据的数据压缩，不仅可以减少数据压缩过程中雷达目标电磁散射特性信息的损失，还可以进一步降低压缩算法的计算复杂度与内存损耗，从而使上述数据压缩过程变得更为精确、快速与高效。

Description

基于基追踪去噪的远场宽带RCS数据采集与压缩方法

技术领域

本发明涉及雷达目标RCS数据获取与处理技术领域，尤其涉及一种基于基追踪去噪的远场宽带RCS数据采集与压缩方法。

背景技术

目标的雷达散射截面(Radar Cross Section，RCS)是表征雷达目标对于照射电磁波散射能力的一个物理量。早在雷达系统出现之前，利用麦克斯韦方程组，就已经求得了几种典型形状的完纯导体目标的电磁散射精确解。等到20世纪30年代雷达系统出现后，雷达目标就成为雷达收、发闭合回路中的一个重要环节，而雷达散射截面便是表征雷达目标特性的一个最基本、最重要的参数。60年代初发展的洲际导弹真、假弹头识别与反识别技术，以及80年代飞行器的隐身与反隐身技术，使得对雷达散射截面的研究出现了两次高潮。在这一时期，相关领域的研究人员对各类雷达目标进行了大量的静态与动态的测量研究与理论分析，由上述研究成果建立起的理论体系促使雷达目标成为雷达领域的一个独立分支，而对特定目标的雷达散射截面研究则成为该分支中的一个重要方向。

在实际工程应用中，为全面地了解复杂目标的电磁散射特性，通常需要借助宽带高分辨雷达，来对目标物体在各个观测角度上的散射回波频谱分布进行测量。而上述操作会使测量系统长时间地工作在对同一目标RCS数据的采集过程中，从而导致整个系统的工作效率变得很低。此外，由于复杂目标的RCS数据在回波域通常不具备可压缩性，因此，存放上述数据势必会占用大量的存储空间，致使研究人员对实验设施存储设备的维护频率变得更为频繁。所以，为了提高雷达目标RCS数据测量系统的工作效率，同时降低实验设施存储设备的维护成本，亟需一种高效的RCS数据采集与压缩方法。

发明内容

(一)要解决的技术问题

为了解决现有技术问题，本发明提供了一种基于基追踪去噪的远场宽带RCS数据采集与压缩方法。

(二)技术方案

本发明提供了一种基于基追踪去噪的远场宽带RCS数据采集与压缩方法，包括：步骤A：基于观测目标散射中心分布的稀疏度，完成对观测目标远场宽带RCS数据的采集，得到观测目标的远场宽带RCS数据矩阵；步骤B：构建远场宽带RCS数据处理算子；步骤C：将所述远场宽带RCS数据处理算子与用于求解基追踪去噪模型的稀疏重构算法相结合，将所述观测目标的远场宽带RCS数据矩阵转化为观测目标的散射中心，获得所述观测目标散射中心的后向散射系数；以及步骤D：根据远场宽带RCS数据压缩精度的要求，对所述观测目标散射中心进行筛选，将筛选出的观测目标散射中心的后向散射系数与空间位置进行存储，完成远场宽带RCS数据的数据压缩。

(三)有益效果

从上述技术方案可以看出，本发明的基于基追踪去噪的远场宽带RCS数据采集与压缩方法具有以下有益效果：

(1)通过在方位向与高程向上的随机降采样，降低了RCS数据的实际测量总量，能够在保证不过多损失雷达目标电磁散射特性信息完整度的情况下，缩短雷达系统对观测目标RCS数据的测量时间，从而提高整个测量系统的工作效率；

(2)根据基追踪去噪模型，利用结合了远场宽带RCS数据处理算子的稀疏重构算法，来实现对回波域RCS数据的数据压缩，不仅可以减少数据压缩过程中雷达目标电磁散射特性信息的损失，还可以进一步降低压缩算法的计算复杂度与内存损耗，从而使上述数据压缩过程变得更为精确、快速与高效。

附图说明

图1为本发明实施例的基于基追踪去噪的远场宽带RCS数据采集与压缩方法的流程图；

图2为远场平面波的入射方向相对于观测目标所在坐标系uvw的相对位置示意图；

图3为观测目标散射中心的恢复结果：图3(a)为改进过的SPGL1算法获得的结果，图3(b)为未改进的SPGL1算法获得的结果；

图4为由电磁仿真获得的远场宽带RCS数据与重构远场宽带RCS数据的幅值与相位比较结果：图4(a)中重构远场宽带RCS数据用到的观测目标散射中心是由改进过的SPGL1算法获得的，图4(b)中重构远场宽带RCS数据用到的观测目标散射中心是由未改进的SPGL1算法获得的。

具体实施方式

在基于地基平台的雷达散射截面测量实验里，观测目标的散射中心分布通常具有较好的稀疏性，本发明将雷达目标RCS数据的采集与压缩过程建模为基追踪去噪(BasisPursuit Denoising，BPDN)。通过在方位向与高程向上的随机降采样，来降低对RCS数据的实际测量总量，然后，利用用于求解基追踪去噪模型的稀疏重构算法，把采集获得的回波数据转化为观测目标的散射中心，以实现对RCS数据的数据压缩。为减少数据压缩过程中雷达目标电磁散射特性信息的损失，选用SPGL1算法(基于L1范数惩罚项的谱投影梯度算法)作为重建观测目标散射中心的方法。

SPGL1是一种精确、有效的稀疏信号重构算法。该算法利用基于牛顿迭代的寻根方法，更新在Lasso模型帕累托曲线上的正则化参数取值，从而使Lasso问题的解逐步地逼近基追踪去噪问题的解。而为了快速、高效地获得基追踪去噪问题的解，将结合了远场宽带RCS数据处理算子的谱投影梯度算法应用于上述Lasso问题序列的求解过程之中。改进过的SPGL1算法不仅保留了原算法求解精度高、收敛速度快等优点，还进一步降低了算法本身的计算复杂度与内存损耗。因此，该算法可以快速、高效地完成对远场宽带RCS数据的数据压缩工作。

为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本发明做进一步地详细说明。

本发明的实施例提供了一种基于基追踪去噪的远场宽带RCS数据采集与压缩方法，如图1所示，其具体内容包括：

步骤A：基于观测目标散射中心分布的稀疏度，完成对观测目标远场宽带RCS数据的采集，得到观测目标的远场宽带RCS数据矩阵。

步骤A具体包括：

子步骤A1：在雷达目标特征信息数据库中查找与观测目标相类似的已有目标模型，并基于该模型的空间几何形状与电特性参数，利用电磁计算方法估计观测目标散射中心分布的稀疏度；

子步骤A2：根据观测目标散射中心分布的稀疏度，在满足稀疏信号重构条件的情况下，确定雷达系统在方位向与高程向上对观测目标RCS数据的随机降采样率，然后，依照上述降采样率，在预先设置的满采样备选网格点上，挑选出用于实现随机降采样的空间采样位置；

子步骤A3：在挑选出的空间采样位置上，对观测目标远场宽带RCS数据进行采集，得到观测目标的远场宽带RCS数据矩阵。

在方位向与高程向上的随机降采样操作降低了RCS数据的实际测量总量，能够在保证不过多损失雷达目标电磁散射特性信息完整度的情况下，缩短雷达系统对观测目标RCS数据的测量时间，从而提高整个测量系统的工作效率。

步骤B：将向量支集搜索方法与滤波逆投影原理引入到远场宽带RCS数据的处理过程中，推导出用于降低稀疏重构算法计算复杂度与内存损耗的远场宽带RCS数据处理算子。

在用于求解基追踪去噪模型的传统稀疏重构算法中，一般存在包含了观测矩阵的矩阵-向量乘法运算，由于该运算的时间复杂度与空间复杂度均为平方阶，因此当观测矩阵维度较大时，包含上述运算的稀疏重构算法会在基追踪去噪模型求解过程中花费大量的运算时间，并占用较多的系统内存。为了能有效地解决该问题，我们用远场宽带RCS数据处理算子替代包含了观测矩阵的矩阵-向量乘法运算在上述算法中的功能。

步骤B具体包括：

子步骤B1：根据向量支集搜索方法的主要思想，构建出远场宽带RCS数据处理算子的回波数据生成项，该回波数据生成项的表达式为：

其中，表示用于搜索观测目标后向散射系数向量中所有非零元素位置的算子，即这里，表示由观测目标后向散射系数向量x中所有非零元素的序号构成的集合；表示利用向量x中元素序号属于集合的观测目标的后向散射系数生成远场宽带RCS数据的算子，即

这里，表示观测目标的远场宽带RCS数据向量，L表示每个步进频脉冲信号的频点总数，M表示在方位向与高程向上的空间采样位置总数，k_l＝4πf_l/c表示对应第l个发射脉冲载波频率f_l的波数，f_l表示第l个脉冲的载波频率，且f_l＝f₁+(l-1)Δf，f₁为第1个脉冲的载波频率，其取值由雷达系统的系统参数和回波数据的应用要求决定，Δf为频率的步进长度，c表示电磁波在自由空间的传播速度。如图2所示，远场平面波的入射方向相对于观测目标所在坐标系uvw的俯仰角与方位角分别为θ_m和且-π/2≤θ_m≤π/2、R_m，n表示观测目标的第m个散射中心位置在远场平面波入射方向上的投影距离，其中，当雷达系统对三维空间坐标为(u_n，v_n，w_n)的散射中心进行距离向、方位向与高程向上的远场宽带RCS数据测量时，投影距离R_m，n的计算式可以表示成当雷达系统对二维空间坐标为(u_n，v_n)的散射中心进行距离向与方位向上的远场宽带RCS数据测量时，投影距离R_m，n的计算式可以表示成

子步骤B2：根据滤波逆投影的基本原理，推导出远场宽带RCS数据处理算子的散射中心重建项的表达式为：

其中，表示把观测目标的远场宽带RCS数据向量重新排列成与之对应的数据矩阵的算子；IFFT_r(·)表示对矩阵J⊙Y中的每一列元素进行逆快速傅里叶变换的算子，为积分坐标变换所用Jacobi行列式构成的矩阵，其中，当雷达系统对散射中心进行距离向、方位向与高程向上的远场宽带RCS数据测量时，J中元素J_l，m等于当雷达系统对散射中心进行距离向与方位向上的远场宽带RCS数据测量时，J中元素J_l，m等于k_l；C_grid(·)表示将距离向上的时域投影插值到目标区域空间网格点处的算子；表示对在不同空间采样位置上获得的子图像进行相干累加的算子；表示把观测目标的后向散射系数矩阵重新排列成向量的算子；表示时域插值操作前对距离向上的时域投影进行相位补偿的向量，其元素取值为ψ₁(l)＝exp(-jπ(L-1)l/L)；ψ_2，m表示时域插值操作后对目标区域空间网格点处的后向散射系数进行相位补偿的矩阵，该矩阵共包含N个元素，当雷达系统对三维空间坐标为(u_n，v_n，w_n)的散射中心进行距离向、方位向与高程向上的远场宽带RCS数据测量时，在相位补偿矩阵ψ_2，m中元素ψ_2，m(u_n，v_n，w_n)的计算式可以表示成当雷达系统对二维空间坐标为(u_n，v_n)的散射中心进行距离向与方位向上的远场宽带RCS数据测量时，在相位补偿矩阵ψ_2，m中元素ψ_2，m(u_n，v_n)的计算式可以表示成这里，f_c表示步进频脉冲信号的中心频率，即f_c＝(f₁+f_L)/2，N的取值与观测目标后向散射系数向量的维度相等；⊙表示哈达马乘法运算。

步骤C：将远场宽带RCS数据处理算子与用于求解基追踪去噪模型的稀疏重构算法相结合，利用该算法将观测目标的远场宽带RCS数据矩阵转化为观测目标的散射中心，获得观测目标散射中心的后向散射系数。

为便于说明上述操作的具体流程，下面将以结合了远场宽带RCS数据处理算子的SPGL1算法作为示例，解释如何通过求解基追踪去噪模型，来完成从远场宽带RCS数据到观测目标散射中心的转化工作。

步骤C具体包括：

子步骤C1：将观测目标的远场宽带RCS数据矩阵输入到改进过的SPGL1算法中，同时对所要用到的参数与临时变量进行初始化。

子步骤C1具体包括：将观测目标的远场宽带RCS数据矩阵转换成由其M个列向量拼接成的远场宽带RCS数据向量并将远场宽带RCS数据向量加性噪声电平σ∈[0，||y||₂)、最优输出容限T_op∈[0，∞)、充分下降系数η∈(0，1)、Barzilai-Borwein步长取值的下界α_min与上界α_max(0＜α_min＜α_max)、第一、二、三层迭代运算的最大迭代次数N_iter1、N_iter2与N_iter3(N_iter1≥1、N_iter2≥1、N_iter3≥1)输入到改进过的SPGL1算法中。然后，对所要用到的参数与临时变量进行初始化，使第一临时变量x_p的初始值为x₀＝0，第二临时变量r_p的初始值为r₀＝y，正则化参数τ_p的初始值为τ₀＝0，Barzilai-Borwein步长α_q的初始值为α₀＝α_max，第一层迭代运算的迭代次数p的初始值为p＝1，设置观测目标后向散射系数输出向量为第三临时变量为第四临时变量为第五临时变量为第六临时变量为第七临时变量为第八临时变量为Δx，第九临时变量为Δg，梯度下降步长为α，第二层迭代运算的迭代次数为q，第三层迭代运算的迭代次数为h，并执行子步骤C2。

子步骤C2：判断p≤N_iter1是否成立，若成立，则执行子步骤C3；否则，将x_p-1作为最终结果进行输出，即输出结果获得观测目标的散射中心，步骤C执行完毕。

子步骤C3：判断|||r_p-1||₂-σ|/||r_p-1||₂≤T_op是否成立，若成立，则将x_p-1作为最终结果进行输出，即输出结果获得观测目标的散射中心，步骤C执行完毕；否则，更新正则化参数，使τ_p＝τ_p-1+(σ-φ(r_p-1))/φ′(r_p-1)，并执行子步骤C4。

子步骤C4：判断τ_p＜τ_p-1是否成立，若成立，则使并执行子步骤C5；否则，使并执行子步骤C5。其中，由远场宽带RCS数据处理算子的回波数据生成项(公式(1))计算得到。

子步骤C5：使第二层迭代运算的迭代次数q＝1，并执行第二层迭代运算过程。其中，由远场宽带RCS数据处理算子的散射中心重建项(公式(2))计算得到。

子步骤C6：使第一层迭代运算的迭代次数p＝p+1，返回执行子步骤C2。

所述执行第二层迭代运算过程具体包括：

子分步骤C5a：判断q≤N_iter2是否成立，若成立，则执行子分步骤C5b；否则，执行子步骤C6；

子分步骤C5b：判断是否成立，若成立，则执行子步骤C6；否则，使α＝α_q-1，第三层迭代运算的迭代次数h＝1，并执行第三层迭代运算过程。

子分步骤C5c：使并执行子分步骤C5d。其中，由远场宽带RCS数据处理算子的散射中心重建项(公式(2))计算得到。

子分步骤C5d：判断Δx^HΔg≤0是否成立，若成立，则使α_q＝α_max，第二层迭代运算的迭代次数q＝q+1，并返回执行子分步骤C5a；否则，使第二层迭代运算的迭代次数q＝q+1，并返回执行子分步骤C5a。

所述执行第三层迭代运算过程具体包括：

子分步骤C5b1：判断h≤N_iter3是否成立，若成立，则 并执行子分步骤C5b2；否则，执行子分步骤C5c。其中，由远场宽带RCS数据处理算子的回波数据生成项(公式(1))计算得到。

子分步骤C5b2：判断是否成立，若成立，则执行子分步骤C5c；否则，使α＝α/2，第三层迭代运算的迭代次数h＝h+1，并返回执行子分步骤C5b1。

其中，子步骤C3中的φ(r_p-1)利用如下函数进行计算：

φ(r)＝||r||₂ (3)

φ′(r_p-1)利用如下函数进行计算：

子分步骤C5b中的利用如下函数进行计算：

其中，由远场宽带RCS数据处理算子的散射中心重建项(公式(2))计算得到；为观测目标的远场宽带RCS数据向量；||·||₂表示取向量的2-范数；||·||_∞表示取向量的无穷范数；上标H表示取向量的共轭转置。

子步骤C4中的和子分步骤C5b1中的由以下投影算子计算获得：

公式(6)所示的投影算子用于计算输入向量在可行集合{z|||z||₁≤τ}上的投影，该投影算子的数据处理流程如下所述：

步骤P-1：将向量和常数τ＞0输入到投影算子中，并执行步骤P-2；

步骤P-2：判断||β||₁≤τ是否成立，若成立，则直接将β作为最终结果进行输出，即输出结果z＝β；否则，执行步骤P-3；

步骤P-3：初始化迭代次数和临时变量，使i＝1，γ＝0，κ＝0，v＝-τ，并执行步骤P-4；

步骤P-4：将向量β中的元素模值，按照取值大小进行降序排列，得到新的向量即并执行步骤P-5；

步骤P-5：判断i≤N是否成立，若成立，则γ＝v/i，并执行步骤P-6，其中，表示向量中的第i个元素；否则，执行步骤P-7；

步骤P-6：判断是否成立，若成立，则执行步骤P-7；否则，使κ＝γ，i＝i+1，并返回执行步骤P-5；

步骤P-7：利用阈值门限为κ的复数域软阈值函数，对向量β进行处理，然后，将处理过的向量作为最终结果进行输出，即输出结果z＝SoftThreshold(β，κ)，其中，复数域软阈值函数表示的数据处理过程为

上述利用结合了远场宽带RCS数据处理算子的SPGL1算法，实现对回波域RCS数据的数据压缩，不仅可以减少数据压缩过程中雷达目标电磁散射特性信息的损失，还可以进一步降低压缩算法的计算复杂度与内存损耗，从而使上述数据压缩过程变得更为精确、快速与高效。

需要说明的是，以上仅是以SPGL1算法为例进行了说明，但本发明不以此为限，除SPGL1算法外，本发明的远场宽带RCS数据处理算子可以与多种用于求解基追踪去噪模型的稀疏重构算法相结合，实现观测目标远场宽带RCS数据矩阵向观测目标散射中心的转化。

步骤D：根据远场宽带RCS数据压缩精度的要求，对重建的观测目标散射中心进行筛选，将筛选出的观测目标散射中心的后向散射系数与空间位置进行存储，完成对远场宽带RCS数据的数据压缩。

步骤D具体包括：

子步骤D1：初始化在筛选过程中所要用到的参数，并对观测目标散射中心的后向散射系数进行归一化处理。

子步骤D1具体包括：设置远场宽带RCS数据压缩精度ε的取值(ε∈[0，∞))，确定筛选门限的起始值T₁与下降步长ΔT的取值(T₁＜0、ΔT＜0)，设定对强散射中心进行筛选的总次数G的取值(G≥1)，对步骤C获得的观测目标散射中心的后向散射系数，即输出向量中的元素模值进行归一化处理，得到单位为dB的向量即令筛选次数g＝1，并执行子步骤D2；

子步骤D2：搜索向量中取值大于T_g的元素，将取值大于T_g元素的序号构建成集合即然后，生成远场宽带RCS数据向量的估计值并计算远场宽带RCS数据向量的压缩误差执行子步骤D3；

子步骤D3：判断ε_g≤ε是否成立，若成立，则将向量中元素序号属于集合的散射中心的后向散射系数与空间位置存储起来，并同时保存实测获得的远场宽带RCS数据矩阵的最大幅值max(|Y|)，完成对观测目标远场宽带RCS数据的数据压缩；否则，使筛选次数g＝g+1，并执行子步骤D4；

子步骤D4：判断g≤G是否成立，若成立，则执行子步骤D2；否则，将向量中元素序号属于集合的散射中心的后向散射系数与空间位置存储起来，并同时保存实测获得的远场宽带RCS数据矩阵的最大幅值max(|Y|)，完成对观测目标远场宽带RCS数据的数据压缩。

根据上述具体实施方式，下面将利用由电磁仿真获得的圆柱体远场宽带RCS数据，对本实施例所述方法的有益效果进行验证。由电磁仿真获得的圆柱体远场宽带RCS数据的主要参数如表1所示。

表1

工作方式	步进扫频	观测方式	单站观测
				信号起始频率	4GHz	起始方位角	-10°
信号终止频率	6GHz	终止方位角	0°
				信号频率步长	10MHz	方位角步长	0.1°

在对由电磁仿真获得的圆柱体远场宽带RCS数据进行数据压缩时，首先，以75％的方位向降采样率，从满采样数据中随机抽选出用于恢复观测目标散射中心的回波数据，然后，分别利用改进过的SPGL1算法(N-SPGL1)与未改进的SPGL1算法(SPGL1)处理上述回波数据，获得用于数据压缩的观测目标散射中心，最后，根据数据压缩条件，即保留后向散射系数归一化幅值大于-50dB的观测目标散射中心，来对上述两种算法获得的观测目标散射中心进行筛选，从而得到用于最终存储的散射中心的后向散射系数与空间位置，完成对远场宽带RCS数据的压缩工作。

如图3所示，改进过的SPGL1算法与未改进的SPGL1算法均能较好地恢复观测目标的散射中心，并且由上述散射中心重构的远场宽带RCS数据与由电磁仿真获得的远场宽带RCS数据非常接近(图4)。可是，根据表2给出的上述算法对由电磁仿真获得的远场宽带RCS数据的处理性能，可以看出：在RCS数据压缩率与压缩精度相近时，改进过的SPGL1算法所需的运行时间与系统内存均要远低于未改进的SPGL1算法。因此，结合了远场宽带RCS数据处理算子的SPGL1算法具有更为精确、快速与高效的对远场宽带RCS数据的压缩性能。

表2

	N-SPGL1	SPGL1
			算法程序运行时间	2.96min	39.55min
处理系统内存损耗	4101MB	5611MB
			RCS数据重构误差	0.029	0.037
RCS数据压缩率	166.40	176.53

至此，已经结合附图对本发明实施例进行了详细描述。依据以上描述，本领域技术人员应当对本发明的基于基追踪去噪的远场宽带RCS数据采集与压缩方法有了清楚的认识。

本发明的基于基追踪去噪的远场宽带RCS数据采集与压缩方法，通过在方位向与高程向上的随机降采样，降低了RCS数据的实际测量总量，能够在保证不过多损失雷达目标电磁散射特性信息完整度的情况下，缩短雷达系统对观测目标RCS数据的测量时间，从而提高整个测量系统的工作效率；根据基追踪去噪模型，利用结合了远场宽带RCS数据处理算子的稀疏重构算法，来实现对回波域RCS数据的数据压缩，不仅可以减少数据压缩过程中雷达目标电磁散射特性信息的损失，还可以进一步降低压缩算法的计算复杂度与内存损耗，从而使上述数据压缩过程变得更为精确、快速与高效。

需要说明的是，在附图或说明书正文中，未绘示或描述的实现方式，均为所属技术领域中普通技术人员所知的形式，并未进行详细说明。此外，上述对各元件的定义并不仅限于实施例中提到的各种方式，本领域普通技术人员可对其进行简单地更改或替换，例如：

(1)实施例中提到的方向用语，例如“上”、“下”、“前”、“后”、“左”、“右”等，仅是参考附图的方向，并非用来限制本发明的保护范围；

(2)上述实施例可基于设计及可靠度的考虑，彼此混合搭配使用或与其他实施例混合搭配使用，即不同实施例中的技术特征可以自由组合形成更多的实施例。

以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限制本发明，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种基于基追踪去噪的远场宽带RCS数据采集与压缩方法，其特征在于，包括：

步骤A：基于观测目标散射中心分布的稀疏度，完成对观测目标远场宽带RCS数据的采集，得到观测目标的远场宽带RCS数据矩阵；

步骤B：构建远场宽带RCS数据处理算子；

步骤C：将所述远场宽带RCS数据处理算子与用于求解基追踪去噪模型的稀疏重构算法相结合，将所述观测目标的远场宽带RCS数据矩阵转化为观测目标的散射中心，获得所述观测目标散射中心的后向散射系数；以及

步骤D：根据远场宽带RCS数据压缩精度的要求，对所述观测目标散射中心进行筛选，将筛选出的观测目标散射中心的后向散射系数与空间位置进行存储，完成远场宽带RCS数据的数据压缩。

2.如权利要求1所述的基于基追踪去噪的远场宽带RCS数据采集与压缩方法，其特征在于，所述步骤B具体包括：

子步骤B1：构建所述远场宽带RCS数据处理算子的回波数据生成项以及

子步骤B2：构建所述远场宽带RCS数据处理算子的散射中心重建项

3.如权利要求2所述的基于基追踪去噪的远场宽带RCS数据采集与压缩方法，其特征在于，所述回波数据生成项的表达式为：

其中，表示用于搜索观测目标后向散射系数向量中所有非零元素位置的算子，即这里，表示由观测目标后向散射系数向量x中所有非零元素的序号构成的集合；表示利用观测目标后向散射系数向量x中元素序号属于集合的观测目标的后向散射系数生成远场宽带RCS数据的算子，即

这里，表示观测目标的远场宽带RCS数据向量，L表示每个步进频脉冲信号的频点总数，M表示在方位向与高程向上的空间采样位置总数，k_l＝4πf_l/c表示对应第l个发射脉冲载波频率f_l的波数，f_l表示第l个脉冲的载波频率，且f_l＝f₁+(l-1)Δf，f₁为第1个脉冲的载波频率，Δf为频率的步进长度，c表示电磁波在自由空间的传播速度；R_m，n表示观测目标的第m个散射中心位置在远场平面波入射方向上的投影距离。

4.如权利要求3所述的基于基追踪去噪的远场宽带RCS数据采集与压缩方法，其特征在于，

当雷达系统对三维空间坐标为(u_n，v_n，w_n)的散射中心进行距离向、方位向与高程向上的远场宽带RCS数据测量时，所述投影距离R_m，n的计算式为

当雷达系统对二维空间坐标为(u_n，v_n)的散射中心进行距离向与方位向上的远场宽带RCS数据测量时，所述投影距离R_m，n的计算式为其中，θ_m和分别为远场平面波的入射方向相对于观测目标所在坐标系uvw的俯仰角与方位角。

5.如权利要求2所述的基于基追踪去噪的远场宽带RCS数据采集与压缩方法，其特征在于，所述散射中心重建项的表达式为：

其中，表示把观测目标的远场宽带RCS数据向量重新排列成与之对应的远场宽带RCS数据矩阵的算子；IFFT_r(·)表示对矩阵J⊙Y中的每一列元素进行逆快速傅里叶变换的算子，为积分坐标变换所用Jacobi行列式构成的矩阵；表示将距离向上的时域投影插值到目标区域空间网格点处的算子；表示对在不同空间采样位置上获得的子图像进行相干累加的算子；表示把观测目标的后向散射系数矩阵重新排列成向量的算子；表示时域插值操作前对距离向上的时域投影进行相位补偿的向量，在相位补偿向量ψ₁中的元素取值为ψ₁(l)＝exp(-jπ(L-1)l/L)；Ψ_2，m表示时域插值操作后对目标区域空间网格点处的后向散射系数进行相位补偿的矩阵；⊙表示哈达马乘法运算。

6.如权利要求5所述的基于基追踪去噪的远场宽带RCS数据采集与压缩方法，其特征在于，

当雷达系统对三维空间坐标为(u_n，v_n，w_n)的散射中心进行距离向、方位向与高程向上的远场宽带RCS数据测量时，在Jacobi行列式矩阵J中的元素J_l，m等于在相位补偿矩阵Ψ_2，m中元素ψ_2，m(u_n，v_n，w_n)的计算式为

当雷达系统对二维空间坐标为(u_n，v_n)的散射中心进行距离向与方位向上的远场宽带RCS数据测量时，在Jacobi行列式矩阵J中的元素J_l，m等于k_l，在相位补偿矩阵Ψ_2，m中元素ψ_2，m(u_n，v_n)的计算式为

其中，f_c表示步进频脉冲信号的中心频率，即f_c＝(f₁+f_L)/2；θ_m和分别为远场平面波的入射方向相对于观测目标所在坐标系uvw的俯仰角与方位角。

7.如权利要求2所述的基于基追踪去噪的远场宽带RCS数据采集与压缩方法，其特征在于，所述步骤C具体包括：

子步骤C1：将观测目标的远场宽带RCS数据矩阵转换成由其M个列向量拼接成的远场宽带RCS数据向量输入远场宽带RCS数据向量y、加性噪声电平σ∈[0，||y||₂)、最优输出容限T_op∈[0，∞)、充分下降系数η∈(0，1)、Barzilai-Borwein步长取值的下界α_min与上界α_max、第一、二、三层迭代运算的最大迭代次数N_iter1、N_iter2与N_iter3；使第一临时变量x_p的初始值为x₀＝0，第二临时变量r_p的初始值为r₀＝y，正则化参数τ_p的初始值为τ₀＝0，Barzilai-Borwein步长α_q的初始值为α₀＝α_max，第一层迭代运算的迭代次数p的初始值为p＝1；设置观测目标后向散射系数输出向量为第三临时变量为第四临时变量为第五临时变量为第六临时变量为第七临时变量为第八临时变量为Δx，第九临时变量为Δg，梯度下降步长为α，第二层迭代运算的迭代次数为q，第三层迭代运算的迭代次数为h，并执行子步骤C2；

子步骤C2：判断p≤N_iter1是否成立，若成立，则执行子步骤C3；否则，输出结果获得观测目标散射中心的后向散射系数，步骤C执行完毕；

子步骤C3：判断|||r_p-1||₂-σ|/||r_p-1||2≤T_op是否成立，若成立，输出结果获得观测目标散射中心的后向散射系数，步骤C执行完毕；否则，更新正则化参数，使τ_p＝τ_p-1+(σ-φ(r_p-1))/φ′(r_p-1)，并执行子步骤C4；

子步骤C4：判断τ_p＜τ_p-1是否成立，若成立，则使 并执行子步骤C5；否则，使并执行子步骤C5；

子步骤C5：使第二层迭代运算的迭代次数q＝1，并执行第二层迭代运算过程；以及

子步骤C6：使第一层迭代运算的迭代次数p＝p+1，返回执行子步骤C2。

8.如权利要求7所述的基于基追踪去噪的远场宽带RCS数据采集与压缩方法，其特征在于，所述执行第二层迭代运算过程具体包括：

子分步骤C5a：判断q≤N_iter2是否成立，若成立，则执行子分步骤C5b；否则，执行子步骤C6；

子分步骤C5b：判断是否成立，若成立，则执行子步骤C6；否则，使α＝α_q-1，第三层迭代运算的迭代次数h＝1，并执行第三层迭代运算过程；

子分步骤C5c：使并执行子分步骤C5d；以及

子分步骤C5d：判断Δx^HΔg≤0是否成立，若成立，则使α_q＝α_max，第二层迭代运算的迭代次数q＝q+1，并返回执行子分步骤C5a；否则，使第二层迭代运算的迭代次数q＝q+1，并返回执行子分步骤C5a。

9.如权利要求8所述的基于基追踪去噪的远场宽带RCS数据采集与压缩方法，其特征在于，所述执行第三层迭代运算过程具体包括：

子分步骤C5b1：判断h≤N_iter3是否成立，若成立，则 并执行子分步骤C5b2；否则，执行子分步骤C5c；以及

子分步骤C5b2：判断是否成立，若成立，则执行子分步骤C5c；否则，使α＝α/2，第三层迭代运算的迭代次数h＝h+1，并返回执行子分步骤C5b1。

10.如权利要求7所述的基于基追踪去噪的远场宽带RCS数据采集与压缩方法，其特征在于，所述步骤D具体包括：

子步骤D1：设置远场宽带RCS数据压缩精度ε、筛选门限的起始值T₁、下降步长ΔT、对强散射中心进行筛选的总次数G，对观测目标散射中心的后向散射系数模值进行归一化处理，得到单位为dB的向量即令筛选次数g＝1，执行子步骤D2；

子步骤D2：搜索向量中取值大于T_g的元素，将取值大于T_g元素的序号构建成集合即生成远场宽带RCS数据向量的估计值并计算远场宽带RCS数据向量的压缩误差执行子步骤D3；

子步骤D3：判断ε_g≤ε是否成立，若成立，对向量中元素序号属于集合的散射中心的后向散射系数与空间位置进行存储，并保存远场宽带RCS数据矩阵的最大幅值max(|Y|)，完成观测目标远场宽带RCS数据的数据压缩；否则，使筛选次数g＝g+1，执行子步骤D4；以及

子步骤D4：判断g≤G是否成立，若成立，则执行子步骤D2；否则，对向量中元素序号属于集合的散射中心的后向散射系数与空间位置进行存储，并保存远场宽带RCS数据矩阵的最大幅值max(|Y|)，完成观测目标远场宽带RCS数据的数据压缩。