一种基于近场稀疏成像外推的目标RCS测量方法与流程

本发明涉及微波测量相关技术领域，具体的说，是涉及一种基于近场稀疏成像外推的目标rcs测量方法。

背景技术

雷达散射截面(rcs)测量对于设计和评估武器装备的隐身性能非常重要。常规的rcs测量方法主要包括室外rcs测量、紧缩场测量和室内微波暗室测量，均已发展到相当成熟的阶段。大型室外场和紧缩场均可以满足目标rcs远场测量条件，但是室外场受天气和环境的影响较大，获得高精度数据代价巨大，且保密性不好；紧缩场要求很高的加工精度，使得系统的代价非常昂贵。室内微波暗室测量具有很强的抗干扰性，保密性好，测试效率高，但室内微波暗室的空间有限，对于电大尺寸目标往往难以满足远场测量条件。

针对室内微波暗室rcs测量不满足远场条件的问题，国内外开展了许多利用近场散射数据来获得目标rcs的研究工作。目前的测量方法在提取目标散射中心时，需要通过clean算法迭代提取强散射中心信息，由此引起的积累误差将降低散射中心位置和强度信息的提取精度，此外，现有方法并未考虑天线方向图和距离衰减对回波功率造成的影响，rcs外推结果存在一定的误差。

技术实现要素：

本发明为了解决上述问题，提出了一种基于近场稀疏成像外推的目标rcs测量方法，采用稀疏重建算法进行目标高分辨成像，生成图像中非零像素即为目标的散射中心，避免了clean算法迭代过程带来的积累误差问题，并且在稀疏字典构造过程中考虑天线方向图和距离衰减因素，提高了rcs测量的精度。

为了实现上述目的，本发明采用如下技术方案：

一种基于近场稀疏成像外推的目标rcs测量方法，包括如下步骤：

步骤1：设定近场测量参数：包括近场测量的测试频率、测试角度、天线到转台中心的距离r0，确定测试天线方向图参数g²(θ,f)。

步骤2：根据设定的近场测量的测试频率、测试角度、天线到转台中心距离r0和天线方向图参数构建反映目标实际散射特性的稀疏字典。

步骤3：在近场测试条件下，根据设定的测量参数检测不同频率、不同方位角度下的空暗室和目标近场散射信号数据，所述近场散射信号数据包括幅度和相位信息。

步骤4：在与步骤3同样的测试条件下，利用转台获得不同频率、不同方位角度下的空暗室和定标体散射数据，所述近场散射信号数据包括幅度和相位信息。

步骤5：将步骤3和4测量的近场散射信号数据进行预处理。

步骤6：将处理后的数据利用稀疏重建算法分别生成目标高分辨图像和定标体高分辨图像。

步骤7：根据步骤6生成的目标高分辨图像和定标体高分辨图像分别计算目标的远场散射场和定标体的重建散射场。

步骤8：通过定标体的rcs根据计算的目标体的远场散射数据和定标体的重建散射场计算目标体的rcs。

进一步的，所述步骤2构建反映目标实际散射特性的稀疏字典的方法具体为：将二维近场散射信号数据矩阵和二维成像场景散射系数矩阵按列重排为一维列向量，并将近场散射信号数据采用矩阵表示：

enear＝φa

enear表示近场散射信号数据列向量，a表示成像场景散射系数γ列向量，φ表示稀疏字典，将第i个二维成像场景位置对应散射中心的单位幅度近场散射信号数据列向量作为稀疏字典的第i个列向量构建稀疏字典φ。

进一步的，所述步骤4预处理为滤除噪声和干扰，所述预处理的方法为背景对消和软件时域门方法。

进一步的，所述步骤6利用稀疏重建算法生成目标高分辨图像和定标体高分辨图像的方法为：

步骤601：构建目标矩阵φ1，目标矩阵φ1是稀疏字典的同型矩阵，初始化目标矩阵φ1为空，初始化剩余能量为近场散射信号的初始能量，设定剩余信号能量阈值为近场散射信号的初始能量的0.05；

步骤602：求解剩余信号与稀疏字典φ中列向量最大相关系数的位置索引，并将所述位置索引在稀疏字典φ中对应的向量更新到目标矩阵φ1相应的位置，将所述稀疏字典φ中对应的向量置零；

步骤603：根据目标矩阵φ1计算成像场景散射系数a，然后根据成像场景散射系数a计算剩余信号能量，将剩余信号能量与设定的剩余信号能量阈值比较，如果剩余信号能量高于阈值，执行步骤602，否则执行步骤604；

步骤604：稀疏重建过程结束，计算得到的a即为生成的高分辨图像。

进一步的，所述步骤7具体为：由步骤6生成的目标体高分辨图像和定标体高分辨图像中的非零像素对应提取出来的目标强散射中心，利用这些散射中心分别重构目标的远场散射场和定标体的重建散射场。

一种基于近场稀疏成像外推的目标rcs的测量系统，包括：

用于设定近场测量参数的装置。

用于根据设定的近场测量的测试频率、测试角度、天线到转台中心距离r0和天线方向图参数构建反映目标实际散射特性的稀疏字典的装置。

用于测量近场散射信号数据的测量装置。

用于将近场散射信号数据进行预处理的装置。

用于将处理后的数据利用稀疏重建算法分别生成目标高分辨图像和定标体高分辨图像的装置。

用于将目标高分辨图像和定标体高分辨图像分别计算目标的远场散射场和定标体的重建散射场的装置。

用于通过定标体的rcs根据计算的目标体的远场散射数据和定标体的重建散射场计算目标体的rcs的装置。

与现有技术相比，本发明的有益效果为：

(1)相比于传统基于转台成像的rcs外推方法，本发明避免了clean算法迭代过程中因误差积累带来的rcs外推精度降低问题，同时考虑天线方向图和距离衰减因素，进一步提高了rcs测量精度。

(2)本发明在采样数据不满足奈奎斯特采样定律或存在数据缺损情况下，依然能够保证rcs外推精度。

附图说明

构成本申请的一部分的说明书附图用来提供对本申请的进一步理解，本申请的示意性实施例及其说明用于解释本申请，并不构成对本申请的限定。

图1是本发明的近场测量系统模型图；

图2是本发明的基于近场稀疏成像外推的目标rcs测量方法的流程图。

具体实施方式：

下面结合附图与实施例对本发明作进一步说明。

应该指出，以下详细说明都是示例性的，旨在对本申请提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。

需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本申请的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。

下述实施例为本申请的一种典型的实施方式，如图2所示，一种基于近场稀疏成像外推的目标rcs测量方法，包括如下步骤：

步骤1：确定近场关键参数：包括天线到转台中心的距离r0和测试天线方向图参数。测试天线方向图参数可通过查找说明或进行测试获取，天线到转台中心的距离r0可通过一维距离像直接获得。

步骤2：根据近场测量的测试频率、测试角度、天线到转台中心距离r0和天线方向图参数构建反映目标实际散射特性的稀疏字典。

步骤3：在不满足远场测试条件下，利用转台获得不同频率、不同方位角度下的空暗室和目标散射数据，包括幅度和相位信息。

步骤4：在同样的测试条件下，利用转台获得不同频率、不同方位角度下的空暗室和定标体散射数据，包括幅度和相位信息。

步骤5：将步骤3和4测量的散射数据进行预处理。

步骤6：利用稀疏重建算法生成目标体高分辨图像和定标体高分辨图像。

步骤7：根据步骤6生成的目标体高分辨图像和定标体高分辨图像分别计算目标体的远场散射数据和定标体的重建散射场。

步骤8：通过定标体的rcs根据计算的目标体的远场散射数据和定标体的重建散射场计算目标体的rcs。

所述步骤2构建反映目标实际散射特性的稀疏字典的方法具体为：

如图1所示，由于收发天线放置的很近，因此双站角很小，可以近似为单站系统。单站系统通过接收天线获得的目标近场散射信号为：

如图1所示θ是转台相对于初始位置的转角，目标上的散射中心相对于转台中心的位置坐标为(x,y)，目标上的散射中心对应的目标散射系数为γ(x,y)，f表示发射频率，r表示目标散射点(x,y)到天线的距离，c表示光速，g²(θ,f)为测试天线方向图参数，d表示待成像区域。

将二维近场散射信号数据矩阵和二维成像场景散射系数矩阵按列重排为一维列向量，将公式(1)表示成矩阵形式：

enear＝φa

enear表示近场散射信号数据列向量，a表示成像场景散射系数γ列向量，φ表示稀疏字典，将第i个二维成像场景位置对应散射中心的单位幅度近场散射信号数据列向量作为稀疏字典的第i个列向量构建稀疏字典φ。稀疏字典φ中的第i个列向量表示为：

其中，vec[]表示列向量化。

所述步骤5预处理为滤除噪声和干扰，所述预处理的方法为背景对消和软件时域门方法。

所述步骤6利用稀疏重建算法生成目标高分辨图像和定标体高分辨图像的方法为：

步骤601：构建目标矩阵φ1，目标矩阵φ1是稀疏字典的同型矩阵，初始化目标矩阵φ1为空，初始化剩余能量为近场散射信号的初始能量，设定剩余信号能量阈值为近场散射信号的初始能量的0.05；

步骤602：求解剩余信号与稀疏字典φ中列向量最大相关系数的位置索引，并将所述位置索引在稀疏字典φ中对应的向量更新到目标矩阵φ1相应的位置，将所述稀疏字典φ中对应的向量置零；

步骤603：根据目标矩阵φ1计算成像场景散射系数a，然后根据成像场景散射系数a计算剩余信号能量，将剩余信号能量与设定的剩余信号能量阈值比较，如果剩余信号能量高于阈值，执行步骤602，否则执行步骤604；

步骤604：稀疏重建过程结束，计算得到的a即为生成的高分辨图像。

在高频区，目标散射信号可以表示为若干强散射中心散射信号的和，这体现了目标信号的空域稀疏特性，因而可以利用稀疏重建算法实现目标图像的高分辨重构，本实施例步骤6的具体过程为：

(1)初始化剩余信号r0＝enear，迭代次数n＝1，目标矩阵φ1为空；

(2)求解最大相关系数的位置索引，即其中为稀疏字典φ的第i个列向量；

(3)将稀疏字典φ中ind对应的向量更新到目标矩阵φ1中，并将字典φ中的置零；

(4)求解投影系数：其中h表示共轭转置，并更新剩余信号rn＝enear-φ1a；

(5)n＝n+1，重复过程(2)至过程(4)，直到剩余信号能量小于初始信号能量的0.05，得到目标矩阵φ1。由目标矩阵φ1计算得到的a即为生成的高分辨图像。

进一步的，所述步骤7具体为：

目标rcs远场数据是目标的各散射中心随频率和方位的变化，步骤6生成的目标体高分辨图像和定标体高分辨图像中的非零像素对应提取出来的目标强散射中心，利用这些散射中心分别重构目标的远场散射场和定标体的重建散射场，通过如下公式计算：

其中，fmin和fmax分别为发射频率的最小值和最大值，θ表示转角θ的范围，an表示第n个散射中心的强度，(xn,yn)为第n个散射中心的坐标位置，n为生成高分辨图像中非零像素的个数。

所述步骤8，通过定标体的rcs根据计算的目标体的远场散射数据和定标体的重建散射场计算目标体的rcs。通过公式(4)计算的目标的远场散射场为efar(f,θ)，定标体的重建散射场为e0(f,θ)，通过如下公式计算目标rcs：

其中，σ为被测目标的rcs，σ0为定标体的rcs。

以上所述仅为本申请的优选实施例而已，并不用于限制本申请，对于本领域的技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本申请的保护范围之内。

上述虽然结合附图对本发明的具体实施方式进行了描述，但并非对本发明保护范围的限制，所属领域技术人员应该明白，在本发明的技术方案的基础上，本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。