基于GTD参数提取的跳频RCS重构方法与流程

基于gtd参数提取的跳频rcs重构方法
技术领域
1.本技术涉及雷达技术领域，特别是涉及一种基于gtd参数提取的跳频rcs重构方法。


背景技术：

2.雷达散射截面积(radar cross section,rcs)为雷达科学研究和应用实践提供了数据基础，对雷达目标特性反演、雷达目标成像和雷达目标智能识别等至关重要。随着大数据时代的到来和在现代雷达超分辨技术的发展，人们迫切需要越来越丰富和完备的雷达目标rcs数据。
3.当前，雷达目标rcs数据获取的主要方式包括电磁仿真计算、暗室缩比测量和外场实际测量。与其它两种方式相比，暗室缩比测量方式虽具有真实性好、保密性强等优势，但传统的暗室步进频rcs测量采用扫频模式，需要对每个目标姿态的所有步进频点进行依次扫频，具有测量周期长和数据存储空间大等问题。rcs重构可以在一定程度上解决上述问题，其主要是通过提取散射中心参数来完成的。基于现代谱估计的gtd(geometric theory of diffraction,gtd)散射参数提取已经取得了较好的成果。然而，在实现本发明过程中，发明人发现前述传统的rcs重构技术，存在着无法有效处理非均匀不完备rcs数据的技术问题。


技术实现要素：

4.基于此，有必要针对上述技术问题，提供一种基于gtd参数提取的跳频rcs重构方法、一种基于gtd参数提取的跳频rcs重构装置、一种计算机设备以及一种计算机可读存储介质，能够有效处理非均匀不完备rcs数据，实现rcs重构。
5.为了实现上述目的，本发明实施例采用以下技术方案：
6.一方面，本发明实施例提供一种基于gtd参数提取的跳频rcs重构方法，包括步骤：
7.对跳频模式下暗室步进频rcs测量数据进行稀疏建模，得到雷达目标的后向散射场gtd散射模型的稀疏表征；
8.采用irls算法对gtd散射模型的矩阵形式进行优化求解，得到稀疏最优解；
9.对稀疏最优解进行门限截断处理后，进行gtd散射参数提取，得到等效散射中心的散射参数；散射参数包括散射位置、散射强度和散射类型；
10.将散射参数代入gtd散射模型，得到暗室rcs测量跳频模式下的rcs重构结果。
11.另一方面，还提供一种基于gtd参数提取的跳频rcs重构装置，包括：
12.稀疏处理模块，用于对跳频模式下暗室步进频rcs测量数据进行稀疏建模，得到雷达目标的后向散射场gtd散射模型的稀疏表征；
13.优化求解模块，用于采用irls算法对gtd散射模型的矩阵形式进行优化求解，得到稀疏最优解；
14.参数提取模块，用于对稀疏最优解进行门限截断处理后，进行gtd散射参数提取，
得到等效散射中心的散射参数；散射参数包括散射位置、散射强度和散射类型；
15.重构处理模块，用于将散射参数代入gtd散射模型，得到暗室rcs测量跳频模式下的rcs重构结果。
16.又一方面，还提供一种计算机设备，包括存储器和处理器，存储器存储有计算机程序，处理器执行计算机程序时实现任一项的上述基于gtd参数提取的跳频rcs重构方法的步骤。
17.再一方面，还提供一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现任一项的上述基于gtd参数提取的跳频rcs重构方法的步骤。
18.上述技术方案中的一个技术方案具有如下优点和有益效果：
19.上述基于gtd参数提取的跳频rcs重构方法，将稀疏重构理论与gtd散射模型相结合对雷达目标rcs数据进行建模，利用rcs数据在空域上是稀疏的这一事实，基于gtd散射模型给出了跳频模式下雷达目标rcs数据的稀疏表征，并且将gtd散射参数提取问题转化为l
p
范数优化问题。进一步采用以l2范数为准则的迭代加权最小二乘(iteratively
‑
reweighed
‑
least
‑
squares，irls)算法对上述l
p
范数进行求解，实现目标不同姿态下的gtd散射参数提取。最后，将提取的散射参数代入gtd散射模型即可实现rcs重构。能够在rcs数据非均匀不完备的条件下提取gtd散射参数，进而实现了跳频模式下基于gtd散射参数提取的rcs重构的目的，对分析数据受限条件下的目标散射特性、对缩减暗室步进频rcs测量的周期和和扩增雷达目标rcs数据有较高的工程应用价值。
附图说明
20.图1为一个实施例中基于gtd参数提取的跳频rcs重构方法的流程示意图；
21.图2为一个实施例中rcs数据稀疏建模处理的流程示意图；
22.图3为一个实施例中采用irls算法进行优化求解的流程示意图；
23.图4为一个实施例中基于gtd参数提取的跳频rcs重构方法的仿真实施流程示意图；
24.图5为一个实施例中信噪比在25db条件下的仿真散射参数的平均结果示意图；其中，(a)为稀疏度ρ＝100.0％时100次mont carlo仿真散射参数的平均结果，(b)为稀疏度ρ＝50.0％时100次mont carlo仿真散射参数的平均结果；
25.图6为一个实施例中信噪比在25db条件下的仿真rcs重构的平均结果示意图；其中，(a)为稀疏度ρ＝100.0％时100次mont carlo仿真rcs重构的平均结果，(b)为稀疏度ρ＝50.0％时100次mont carlo仿真rcs重构的平均结果；
26.图7为一个实施例中两段锥组合体目标的cad模型示意图；
27.图8为一个实施例中雷达方位角θ＝1.05
°
时对两段锥组合体目标的参数提取与重构结果示意图；其中，(a)为散射参数的提取结果，(b)为rcs序列的重构结果；
28.图9为一个实施例中两段锥组合体目标的rcs序列示意图；其中，(a)为原始rcs序列，(b)为重构后的rcs序列(稀疏度ρ＝50％)；
29.图10为一个实施例中基于gtd参数提取的跳频rcs重构装置的模块结构示意图。
具体实施方式
30.为了使本技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本技术进行进一步详细说明。应当理解，此处描述的具体实施例仅仅用以解释本技术，并不用于限定本技术。
31.除非另有定义，本文所使用的所有的技术和科学术语与属于本技术的技术领域的技术人员通常理解的含义相同。本文中在本技术的说明书中所使用的术语只是为了描述具体的实施例的目的，不是旨在于限制本技术。
32.另外，本发明各个实施例之间的技术方案可以相互结合，但是必须是以本领域普通技术人员能够实现为基础，当技术方案的结合出现相互矛盾或无法实现时，应当认为这种技术方案的结合不存在，也不在本发明要求的保护范围之内。
33.传统的rcs重构技术无法用于处理非均匀不完备rcs数据，受到奈奎斯特采样定理的制约，且其频率依赖因子在级数展开时会引入了高阶近似误差。此外，通过测量获取完备的rcs数据也需要消耗巨大的时空开销。因此，研究如何有效地缩减暗室rcs测量的时空开销、扩增雷达rcs数据具有较大的应用前景。
34.本发明针对传统的rcs重构技术存在着的无法有效处理非均匀不完备rcs数据的技术问题，提出了一种基于gtd参数提取的跳频rcs重构方法，总体设计构思是将稀疏重构理论与gtd散射模型相结合对雷达目标rcs数据进行建模，利用rcs数据在空域上是稀疏的这一事实，基于gtd散射模型给出了跳频模式下雷达目标rcs数据的稀疏表征，并将gtd散射参数提取问题转化为l
p
范数优化问题。进一步采用以l2范数为准则的irls算法对上述l
p
范数优化问题进行求解，以实现目标不同姿态下的gtd散射参数提取。最后，将提取的散射参数代入gtd散射模型即可实现rcs重构。
35.请参阅图1，一方面，本发明提供一种基于gtd参数提取的跳频rcs重构方法，包括如下步骤s12至s18：
36.s12，对跳频模式下暗室步进频rcs测量数据进行稀疏建模，得到雷达目标的后向散射场gtd散射模型的稀疏表征。
37.可以理解，在光学区，雷达目标电磁散射可等效为多个局部强散射中心的相干合成。
38.请参阅图2，在一些实施方式中，步骤s12可以包括如下处理步骤s121至s124：
39.s121，基于gtd散射模型对暗室rcs测量数据进行表征。
40.可以理解，一般的gtd散射模型为：
[0041][0042]
其中，y表示雷达目标的后向散射回波，i表示等效散射中心的个数，f0表示发射信号的初始频率，c表示电磁波的传播速度，a
i
表示第i个散射中心的散射强度，α
i
＝{
‑
1,
‑
0.5,0,0.5,1}表示第i个散射中心的散射类型，r
i
表示第i个散射中心在雷达视线上的相位距离。
[0043]
具体的，步进频暗室rcs测量系统通过依次发射相干步进频信号工作。假设相干步进频信号的频率间隔为δf，则第n个步进频信号的频率为f
n
＝f0 nδf，(n＝0,1,
…
,n
‑
1)。因此，步进频暗室rcs测量等效于对式(1)作均匀离散化处理，故，暗室步进频rcs测量对应
的散射数据采用如下gtd散射模型的稀疏表征：
[0044][0045]
其中，y
n
表示雷达目标在频率f
n
处的rcs数值，δf表示相干步进频信号的频率间隔，则第n个步进频信号的频率为f
n
＝f0 nδf，(n＝0,1,
…
,n
‑
1)。测试信号的带宽为b＝nδf。
[0046]
s122，对雷达最大无模糊距离进行均匀网格化处理。
[0047]
可以理解，常规雷达的径向距离分辨率为δ
r
＝c/2b。为了获取更为精确的散射位置分布，对最大无模糊距离r
u
＝c/2δf作均匀网格化处理。
[0048]
s123，根据均匀网格化处理后的散射位置的候选网格对均匀离散化处理后的gtd散射模型进行推导，得到网格化后的gtd散射模型。
[0049]
具体的，假设l是距离分辨率的提高倍数，则均匀网格划分后的分辨率为散射位置的候选网格为其中当网格化后的距离分辨率足够高时，所有散射中心均可以近似分布在离散距离网格上，即暗室步进频rcs测量数据采用如下离散网格化后的gtd散射模型表征，因此式(2)可以进一步推导为：
[0050][0051]
因此，网格化后的gtd散射模型即为式(3)所示。y
n
表示雷达目标在频率f
n
处的rcs数值，r
m
表示第m个距离网格对应的散射位置、表示第m个距离网格对应的散射强度，α
m
表示第m个距离网格对应的散射类型，ln表示离散距离网格的数量。当时，可以认为该距离单元上不存在强散射分布。
[0052]
s124，用网格化后的gtd散射模型对跳频rcs测量模式下的rcs数据进行稀疏表征并转换为矩阵形式。
[0053]
可以理解，步进频暗室rcs测量的扫频模式需要对目标的每个姿态进行全频带扫频测量。在待测带宽为b＝nδf的假设下，每个目标姿态均需要进行n次测量。而在跳频rcs测量模式下，只需要从上述n个待测步进频点中随机抽取q个频点进行rcs测量，之后通过rcs重构方法进而实现全频带rcs获取的目的。
[0054]
具体的，假设v为跳频rcs测量模式下实际被测频率的索引集合，则v共包含q个元素且为了便于描述，定义稀疏度ρ：
[0055][0056]
由于同一个距离网格r
m
上可选择5种不同的散射类型考虑系
统加性噪声的影响，将式(3)写成矩阵形式：
[0057]
y＝sφx w
ꢀꢀ
(3)
[0058]
其中，表示随机跳频rcs测量模式下由被测rcs构成的列向量，简称跳频rcs序列，q表示从n个待测步进频点中随机抽取的频点数量。
[0059]
表示由ln个离散距离网格内所有候选散射中心构成的列向量，包含了散射位置、散射强度和散射类型等信息。
[0060]
表示跳频rcs测量模式对应的下采样矩阵，且v
q
表示跳频rcs测量模式下实际被测频点q的索引。
[0061]
表示加性复高斯白噪声列向量。表示的是一个冗余字典矩阵：
[0062][0063]
其中，表示的是一个子矩阵，共包含了5种候选的散射类型。因此，可进一步写成：
[0064][0065]
式中是冗余字典矩阵φ的基向量。根据式(3),可被下式(8)确定：
[0066][0067]
为了便于分析，记：
[0068]
ψ＝sφ
ꢀꢀ
(9)
[0069]
式中ψ∈q
×
5(ln
‑
1)表示随机跳频字典矩阵。
[0070]
s14，采用irls算法对gtd散射模型的矩阵形式进行优化求解，得到稀疏最优解。
[0071]
可以理解，列向量x表征了雷达目标rcs序列的稀疏分量，通过跳频rcs序列对x进行优化求解即可实现gtd散射参数提取。式(5)是一个线性欠定逆问题，包含了无穷多组解。为了避免直接对l0范数优化求解，irls算法将l0范数优化问题转化为l
p
范数优化问题：
[0072][0073]
式中在0≤p≤1的前提下，式(10)又可以通过l2范数近似逼近，进一步转化为如下优化问题：
[0074][0075]
假设ω
k(t)
是x
k
第t次迭代更新的加权系数，取为：
[0076][0077]
式(12)两边同时乘以ω
k(t)
并求和可得：
[0078][0079]
式(13)表明irls算法可以通过优化最小l2范数快速近似求解最小l
p
范数优化问题。因此，通过对目标函数(11)进行交替更新即可求得gtd散射参数对应的稀疏最优解
[0080]
请参阅图3，在一些实施方式中，上述的步骤s14具体可以包括如下处理步骤s141至s144：
[0081]
s141，对矩阵形式的迭代初始次数、正则化参数和列向量进行初始化。
[0082]
具体的，采用irls算法对解向量x进行求解时，首先需要对迭代初始次数t、正则化参数ε和解向量x进行初始化：
[0083][0084]
s142，更新irls算法的迭代加权系数矩阵。
[0085]
可以理解，理论上根据式(12)即可确定对应的权重系数ω
k(t)
，但为了避免因ω
k(t)
→
0时式(17)中矩阵运算(ψ
h
w
(t)
ψ)
‑1出现病态，通常会引入正则化参数ε∈[0,1]。具体的，权重系数ω
k(t)
可以通过下式确定：
[0086][0087]
相应地，权重系数矩阵w
(t)
可由构成：
[0088][0089]
式中ε
(t
‑
1)
表示第(t
‑
1)次迭代时的正则化参数，表示第(t
‑
1)次迭代时列向量x
(t
‑
1)
第k个元素的大小，ln表示离散距离网格的数量。
[0090]
s143，更新矩阵形式的列向量。
[0091]
具体的，正则化参数的更新方式，将式(16)中的权重系数矩阵代入式(17)即可实现解向量x的第t次更新x
(t)
：
[0092]
x
(t)
＝w
(t)
ψ
h
(ψ
h
w(
t
)ψ)
‑1y
ꢀꢀ
(17)
[0093]
式中x
(t)
表示第t次迭代的列向量，ψ表示随机跳频字典矩阵，(
·
)
h
表示共轭转置运算，(
·
)
‑1表示矩阵求逆运算。
[0094]
s144，更新正则化参数，若相邻解向量差值的二范数与解向量二范数的比值不超过设定的求解精度或者迭代次数超过最大迭代次数，则终止迭代并将最终更新的列向量作
为稀疏最优解。
[0095]
具体的，正则化参数的更新方式为：
[0096][0097]
其中，ε
(t)
表示第t次迭代的正则化参数。为了加快x
(t)
的收敛速度，根据式(18)更新正则化参数ε
(t)
。
[0098]
依次迭代更新式(15)、(16)、(17)和(18)并执行t＝t 1操作，直至相邻解向量差值的二范数与解向量二范数的比值，也即||x
(t)
‑
x
(t
‑
1)
||2/||x
(t)
||2＜δ(δ为求解精度，例如可以但不限于设定为10
‑6)，或者t≥t(t为最大迭代次数，例如可以但不限于设为300)时终止迭代，则最终更新的x
(end)
可视为稀疏最优解
[0099]
s16，对稀疏最优解进行门限截断处理后，进行gtd散射参数提取，得到等效散射中心的散射参数；散射参数包括散射位置、散射强度和散射类型。
[0100]
可以理解，为了确保gtd散射参数提取的有效性，对稀疏最优解进行门限截断，以去因除旁瓣或噪声干扰造成的虚拟散射中心。
[0101]
具体的，对稀疏最优解进行门限截断处理的方式为：
[0102][0103]
其中，表示门限截断后的解向量的第k个元素，η表示截断门限，表示所述稀疏最优解，表示所述稀疏最优解的第k个元素。因此，中的任一非零元素均表征一个散射中心。
[0104]
s18，将散射参数代入gtd散射模型，得到暗室rcs测量跳频模式下的rcs重构结果。
[0105]
可以理解，假定x中的第k个元素表征第i个散射中心，则该散射中心所在的距离网格m为：
[0106][0107]
式中floor(
·
)表示向下取整函数，mod(
·
)表示取余函数。
[0108]
相应地，其散射位置可以确定为：
[0109][0110]
其散射强度为：
[0111][0112]
根据式(6)、式(7)和式(8)，其散射类型可以判断为：
[0113]
[0114]
等效散射中心的位置强度和类型均可以根据式(21)、(22)和(23)分别确定。将上述提取的散射参数代入式(2)即可实现暗室rcs测量跳频模式下的rcs重构：
[0115][0116]
上述基于gtd参数提取的跳频rcs重构方法，将稀疏重构理论与gtd散射模型相结合对雷达目标rcs数据进行建模，利用rcs数据在空域上是稀疏的这一事实，基于gtd散射模型给出了跳频模式下雷达目标rcs数据的稀疏表征，并且将gtd散射参数提取问题转化为l
p
范数优化问题。
[0117]
进一步采用以l2范数为准则的迭代加权最小二乘算法对上述l
p
范数进行求解，实现目标不同姿态下的gtd散射参数提取。最后，将提取的散射参数代入gtd散射模型即可实现rcs重构。能够在rcs数据非均匀不完备的条件下提取gtd散射参数，进而实现了跳频模式下基于gtd散射参数提取的rcs重构的目的，对分析数据受限条件下的目标散射特性、对缩减暗室步进频rcs测量的周期和和扩增雷达目标rcs数据有较高的工程应用价值。
[0118]
在一个实施例中，如图4所示，为了更直观且全面地说明上述基于gtd参数提取的跳频rcs重构方法，下面给出了其中一个应用上述基于gtd参数提取的跳频rcs重构方法进行仿真实验的示例。
[0119]
需要说明的是，本说明书中给出的实施示例仅为示意性的，并非为本发明具体实施方式的唯一限定，本领域技术人员可以在本发明提供的实施示例的示意下，同理采用上述提供的基于gtd参数提取的跳频rcs重构方法，实现对不同应用场景下的实验与应用分析。
[0120]
首先，以仿真数据进行实验测试。
[0121]
表1列出了仿真目标的散射参数，该仿真目标共包含4个散射中心。假定雷达工作在x波段，初始频率f0＝8ghz，雷达带宽为b＝2ghz，步进频率间隔为δf＝20mhz，高分辨率提升倍数设置为l＝4。
[0122]
表1
[0123][0124]
表2列出了信噪比在25db条件下100次mont carlo仿真散射参数仿真结果的均值。可见，本技术可以在稀疏度为ρ＝100.0％和ρ＝50.0％的条件下实现散射参数提取，且两种条件下散射参数的提取结果相近且与真实值相逼近。
[0125]
表2
[0126][0127]
图5为信噪比在25db条件下100次mont carlo仿真散射参数的平均结果。
[0128]
图5中(a)和(b)分别为信噪比为25db时在稀疏度ρ＝100.0％和ρ＝50.0％的条件下100次mont carlo仿真散射参数的平均结果。可见，两种稀疏条件下s1、s2、s3、s4和s5的散射位置、散射强度和散射类型均被正确估计，仅散射强度存在一定偏差，验证了本技术可以在rcs数据非均匀不完备的条件下提取gtd散射参数。
[0129]
图6中(a)和(b)分别为信噪比为25db时在稀疏度ρ＝100.0％和ρ＝50.0％的条件下100次mont carlo仿真rcs重构的平均结果。可见，两种稀疏条件下重构的rcs序列和和真实的rcs序列高度吻合，仅在峰值和峰谷处存在一定偏差，验证了本技术可以在rcs数据非均匀不完备的条件下重构rcs序列。
[0130]
进一步的，以电磁计算数据进行试验验证。
[0131]
图7为两段锥组合体电磁计算目标的cad模型。雷达初始视线平行于z轴并指向电磁计算目标的鼻锥方向。雷达仿真参数设如下：初始频率f0＝8ghz，雷达带宽为b＝2ghz，步进频率间隔为δf＝20mhz，方位角间隔为δθ＝0.05
°
，方位累积角为0～6
°
[0132]
图8中(a)为雷达方位角θ＝1.05
°
时在疏度ρ＝50.0％的条件下gtd散射参数的提取结果。可见，散射回波主要由5个等效散射中心构成，其散射位置和散射强度散与雷达一维距离像(range profile,rp)相吻合。此外，散射类型为镜面反射和边缘绕射的散射分量占据主导，与理论分析相吻合。图8中(b)为雷达方位角θ＝1.05
°
时在疏度ρ＝50.0％的条件下rcs序列的重构结果。可见，重构的rcs序列与原始rcs序列高度吻合，仅在峰值和峰谷处存在一定偏差，进一步验证了本技术的效性。
[0133]
图9中(a)和(b)分别为两段锥组合体目标的原始rcs序列和在稀疏度ρ＝50％条件下重构后的rcs序列。可见，重构后的rcs序列与原始rcs序列高度吻合，其峰值信噪比高达(peak
‑
signal
‑
to
‑
noise ratio,psnr)39.0999db。主观的视觉效果和客观的psnr指标均验证了本技术的有效性。
[0134]
综上，本技术可以在非均匀不完备rcs数据的条件下提取gtd散射参数，进而实现跳频模式下rcs重构。本技术对分析数据受限条件下的目标散射特性、对缩减暗室步进频rcs测量的周期和和扩增雷达目标rcs数据具有较高的工程应用价值。
[0135]
应该理解的是，虽然图1至图3流程图中的各个步骤按照箭头的指示依次显示，但
是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且图1至图3的至少一部分步骤可以包括多个子步骤或者多个阶段，这些子步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些子步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
[0136]
请参阅图10，在一个实施例中，还提供了一种基于gtd参数提取的跳频rcs重构装置100，包括稀疏处理模块11、优化求解模块13、参数提取模块15和重构处理模块17。其中，稀疏处理模块11用于对跳频模式下暗室步进频rcs测量数据进行稀疏建模，得到雷达目标的后向散射场gtd散射模型的稀疏表征。优化求解模块13用于采用irls算法对gtd散射模型的矩阵形式进行优化求解，得到稀疏最优解。参数提取模块15用于对稀疏最优解进行门限截断处理后，进行gtd散射参数提取，得到等效散射中心的散射参数；散射参数包括散射位置、散射强度和散射类型。重构处理模块17用于将散射参数代入gtd散射模型，得到暗室rcs测量跳频模式下的rcs重构结果。
[0137]
上述基于gtd参数提取的跳频rcs重构装置100，通过各模块的协作，将稀疏重构理论与gtd散射模型相结合对雷达目标rcs数据进行建模，利用rcs数据在空域上是稀疏的这一事实，基于gtd散射模型给出了跳频模式下雷达目标rcs数据的稀疏表征，并且将gtd散射参数提取问题转化为l
p
范数优化问题。进一步采用以l2范数为准则的迭代加权最小二乘irls算法对上述l
p
范数进行求解，实现目标不同姿态下的gtd散射参数提取。最后，将提取的散射参数代入gtd散射模型即可实现rcs重构。能够在rcs数据非均匀不完备的条件下提取gtd散射参数，进而实现了跳频模式下基于gtd散射参数提取的rcs重构的目的，对分析数据受限条件下的目标散射特性、对缩减暗室步进频rcs测量的周期和和扩增雷达目标rcs数据有较高的工程应用价值。
[0138]
在一个实施例中，上述稀疏处理模块11的各模块，还可以用于实现上述基于gtd参数提取的跳频rcs重构方法各实施例中的其他相应子步骤。
[0139]
在一个实施例中，上述优化求解模块13的各模块，还可以用于实现上述基于gtd参数提取的跳频rcs重构方法各实施例中的其他相应子步骤。
[0140]
关于基于gtd参数提取的跳频rcs重构装置100的具体限定，可以参见上文中基于gtd参数提取的跳频rcs重构方法的相应限定，在此不再赘述。上述基于gtd参数提取的跳频rcs重构装置100中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于具体数据处理功能的设备中，也可以软件形式存储于前述设备的存储器中，以便于处理器调用执行以上各个模块对应的操作，前述设备可以是但不限于本领域已有的各型雷达数据计算分析设备。
[0141]
又一方面，还提供一种计算机设备，包括存储器和处理器，存储器存储有计算机程序，处理器执行计算机程序时可以实现以下步骤：对跳频模式下暗室步进频rcs测量数据进行稀疏建模，得到雷达目标的后向散射场gtd散射模型的稀疏表征；采用irls算法对gtd散射模型的矩阵形式进行优化求解，得到稀疏最优解；对稀疏最优解进行门限截断处理后，进行gtd散射参数提取，得到等效散射中心的散射参数；散射参数包括散射位置、散射强度和散射类型；将散射参数代入gtd散射模型，得到暗室rcs测量跳频模式下的rcs重构结果。
[0142]
在一个实施例中，处理器执行计算机程序时还可以实现上述基于gtd参数提取的跳频rcs重构方法各实施例中增加的步骤或者子步骤。
[0143]
再一方面，还提供一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：对跳频模式下暗室步进频rcs测量数据进行稀疏建模，得到雷达目标的后向散射场gtd散射模型的稀疏表征；采用irls算法对gtd散射模型的矩阵形式进行优化求解，得到稀疏最优解；对稀疏最优解进行门限截断处理后，进行gtd散射参数提取，得到等效散射中心的散射参数；散射参数包括散射位置、散射强度和散射类型；将散射参数代入gtd散射模型，得到暗室rcs测量跳频模式下的rcs重构结果。
[0144]
在一个实施例中，计算机程序被处理器执行时，还可以实现上述基于gtd参数提取的跳频rcs重构方法各实施例中增加的步骤或者子步骤。
[0145]
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成的，计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本技术所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(rom)、可编程rom(prom)、电可编程rom(eprom)、电可擦除可编程rom(eeprom)或闪存。易失性存储器可包括随机存取存储器(ram)或者外部高速缓冲存储器。作为说明而非局限，ram以多种形式可得，诸如静态ram(sram)、动态ram(dram)、同步dram(sdram)、双数据率sdram(ddrsdram)、增强型sdram(esdram)、同步链路(synchlink)dram(sldram)、存储器总线式动态随机存储器(rambus dram，简称rdram)以及接口动态随机存储器(drdram)等。
[0146]
以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
[0147]
以上实施例仅表达了本技术的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本技术构思的前提下，还可做出若干变形和改进，都属于本技术保护范围。因此本技术专利的保护范围应以所附权利要求为准。