一种自旋空间目标姿态估计方法、系统、装置及存储介质

1.本发明涉及雷达技术领域，尤其是一种自旋空间目标姿态估计方法、系统、装置及存储介质。


背景技术：

2.自旋空间目标瞬时状态估计对了解自旋空间目标的在轨运行状态以及判断自旋空间目标航天发展态势具有非常重要的意义。利用逆合成孔径雷达(inverse synthetic aperture radar，isar)对空间目标测量获得的图像来精确获得空间目标重要载荷部件的绝对姿态、空间目标整体运动规律的技术可实际应用于包括空间目标故障救援、威胁度评估等民用和军用领域，是目前能够实现自旋空间目标在轨状态估计的实用技术。
3.目前，测定自旋空间目标在轨状态主要存在两种方式：一是通过激光传感器测量空间目标配备的角锥反射器距离变化进而确定目标在轨运行状态，二是设置目标3d模型状态参数生成2d雷达图像，根据其与观测得到的目标图像间近似程度确定目标在轨状态。但这些现有方案大多依赖积累的数据库或方位标定等强先验条件，很少将目标自旋对观测特征影响参数化建模。因而，这些数据驱动的算法在缺少这些观测先验条件时很难应用于非合作目标，在实际空间目标姿态测量中存在局限性。


技术实现要素：

4.本发明的目的在于至少一定程度上解决现有技术中存在的技术问题之一。
5.为此，本发明实施例的一个目的在于提供一种自旋空间目标姿态估计方法，该方法对预设时间段内的多站观测得到的雷达回波信号同步后进行距离
‑
多普勒处理，获得多帧同步的isar观测图像，然后在获得的多帧isar观测图像内提取目标在距离向、多普勒向上的投影尺寸参数，结合各isar雷达系统的观测几何参数建立两个独立的优化模型，使用粒子群优化方法依次求解自旋空间目标的瞬时姿态参数以及自旋参数，完成对自旋空间目标的姿态估计。
6.本发明实施例的另一个目的在于提供一种自旋空间目标姿态估计系统。
7.为了达到上述技术目的，本发明实施例所采取的技术方案包括：
8.第一方面，本发明实施例提供了一种自旋空间目标姿态估计方法，包括以下步骤：
9.通过多个isar雷达系统获取预设时间段内的多个第一回波信号，并对所述第一回波信号进行距离
‑
多普勒处理，得到目标区域的第一雷达图像序列，所述第一雷达图像序列包括多个isar观测图像；
10.在各所述isar观测图像中均提取出两个目标投影长度特征结构，并确定所述目标投影长度特征结构在距离维度上的目标距离维度尺寸和在多普勒维度上的目标多普勒维度尺寸；
11.获取各所述isar观测图像对应的第一俯仰角和第一方位角，并根据所述第一俯仰角和所述第一方位角确定各isar雷达系统的中心雷达视线方向向量；
12.根据所述目标距离维度尺寸和所述中心雷达视线方向向量构建第一优化模型，并利用粒子群算法对所述第一优化模型进行求解得到目标瞬时姿态参数；
13.根据所述目标多普勒维度尺寸和所述目标瞬时姿态参数构建第二优化模型，并利用粒子群算法对所述第二优化模型进行求解得到目标自旋参数，进而根据所述目标瞬时姿态参数和所述目标自旋参数完成对自旋空间目标的姿态估计。
14.进一步地，在本发明的一个实施例中，所述中心雷达视线方向向量的计算公式为：
[0015][0016]
其中，表示中心雷达视线方向向量，θ(t
m
)表示第一俯仰角，φ(t
m
)表示第一方位角，t
m
表示isar观测图像的相干积累时间的中心时刻。
[0017]
进一步地，在本发明的一个实施例中，所述第一优化模型包括第一目标函数和第一约束条件，所述第一目标函数为：
[0018][0019]
其中，n∈[1，2]，p∈[1，p]，p表示isar雷达系统的数量，α
n
和β
n
表示待求解的目标瞬时姿态参数，r
n，p
表示第n个目标投影长度特征结构在第p个isar观测图像的目标距离维度尺寸，l
n
表示第n个目标投影长度特征结构的真实尺寸，表示第p个isar观测图像对应的isar雷达系统的中心雷达视线方向向量；
[0020]
所述第一约束条件为：
[0021]
tanα1tanα2+cos(β1‑
β2)＝0
[0022]
其中，α1、α2、β1和β2表示待求解的目标瞬时姿态参数。
[0023]
进一步地，在本发明的一个实施例中，所述利用粒子群算法对所述第一优化模型进行求解得到目标瞬时姿态参数这一步骤，其具体包括：
[0024]
设置最短移动距离，并在所述第一优化模型可行域内随机生成第一粒子群，所述第一粒子群的可行解为第一个体位置；
[0025]
根据第一代价函数搜索得到当前的粒子最优位置和粒子群最优位置；
[0026]
调整所述第一粒子群内粒子位置后，更新当前的粒子最优位置和粒子群最优位置，直至满足最大迭代步数或最小变化准则；
[0027]
根据当前最优的第一个体位置确定所述目标瞬时姿态参数。
[0028]
进一步地，在本发明的一个实施例中，所述第二优化模型包括第二目标函数，所述第二目标函数为：
[0029][0030]
其中，
n∈[1，2]，p∈[1，p]，p表示isar雷达系统的数量，θ
rot
、φ
rot
和ω
rot
表示待求解的目标自旋参数，表示预设的雷达视线参数，表示第p个isar观测图像对应的isar雷达系统的中心雷达视线方向向量，d
n，p
表示第n个目标投影长度特征结构在第p个isar观测图像的目标多普勒维度尺寸。
[0031]
进一步地，在本发明的一个实施例中，所述利用粒子群算法对所述第二优化模型进行求解得到目标自旋参数这一步骤，其具体包括：
[0032]
设置最短移动距离，并在所述第二优化模型可行域内随机生成第二粒子群，所述第二粒子群的可行解为第二个体位置；
[0033]
根据第二代价函数搜索得到当前的粒子最优位置和粒子群最优位置；
[0034]
调整所述第二粒子群内粒子位置后，更新当前的粒子最优位置和粒子群最优位置，直至满足最大迭代步数或最小变化准则；
[0035]
根据当前最优的第二个体位置确定所述目标自旋参数。
[0036]
进一步地，在本发明的一个实施例中，根据以下公式完成对自旋空间目标的姿态估计：
[0037][0038][0039][0040]
其中，和表示自旋空间目标的特征姿态，α1、α2、β1和β2表示目标瞬时姿态参数，表示自旋空间目标的旋转向量，θ
rot
、φ
rot
和ω
rot
表示目标自旋参数。
[0041]
第二方面，本发明实施例提供了一种自旋空间目标姿态估计系统，包括：
[0042]
第一雷达图像序列获取模块，用于通过多个isar雷达系统获取预设时间段内的多个第一回波信号，并对所述第一回波信号进行距离
‑
多普勒处理，得到目标区域的第一雷达图像序列，所述第一雷达图像序列包括多个isar观测图像；
[0043]
目标距离维度尺寸和目标多普勒维度尺寸确定模块，用于在各所述isar观测图像中均提取出两个目标投影长度特征结构，并确定所述目标投影长度特征结构在距离维度上的目标距离维度尺寸和在多普勒维度上的目标多普勒维度尺寸；
[0044]
中心雷达视线方向向量确定模块，用于获取各所述isar观测图像对应的第一俯仰角和第一方位角，并根据所述第一俯仰角和所述第一方位角确定各isar雷达系统的中心雷达视线方向向量；
[0045]
目标瞬时姿态参数确定模块，用于根据所述目标距离维度尺寸和所述中心雷达视线方向向量构建第一优化模型，并利用粒子群算法对所述第一优化模型进行求解得到目标瞬时姿态参数；
[0046]
目标自旋参数确定模块，用于根据所述目标多普勒维度尺寸和所述目标瞬时姿态参数构建第二优化模型，并利用粒子群算法对所述第二优化模型进行求解得到目标自旋参数，进而根据所述目标瞬时姿态参数和所述目标自旋参数完成对自旋空间目标的姿态估计。
[0047]
第三方面，本发明实施例提供了一种自旋空间目标姿态估计装置，包括：
[0048]
至少一个处理器；
[0049]
至少一个存储器，用于存储至少一个程序；
[0050]
当所述至少一个程序被所述至少一个处理器执行时，使得所述至少一个处理器实现上述的一种自旋空间目标姿态估计方法。
[0051]
第四方面，本发明实施例还提供了一种计算机可读存储介质，其中存储有处理器可执行的程序，所述处理器可执行的程序在由处理器执行时用于执行上述的一种自旋空间目标姿态估计方法。
[0052]
本发明的优点和有益效果将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到：
[0053]
本发明实施例对预设时间段内的多站观测得到的雷达回波信号同步后进行距离
‑
多普勒处理，获得多帧同步的isar观测图像，然后在获得的多帧isar观测图像内提取目标在距离向、多普勒向上的投影尺寸参数，结合各isar雷达系统的观测几何参数建立两个独立的优化模型，使用粒子群优化方法依次求解自旋空间目标的瞬时姿态参数以及自旋参数，完成对自旋空间目标的姿态估计。本发明实施例可实现对自旋空间目标在轨瞬时状态的精确反演，利用多站联合观测对应的观测角度资源弥补了现有单站序列成像解译方式在求解高阶目标动态参数优化方面的不足，避免了某些观测视角下目标投影特征起伏，提升了自旋空间目标姿态估计的准确性。
附图说明
[0054]
为了更清楚地说明本发明实施例中的技术方案，下面对本发明实施例中所需要使用的附图作以下介绍，应当理解的是，下面介绍中的附图仅仅为了方便清晰表述本发明的技术方案中的部分实施例，对于本领域的技术人员来说，在无需付出创造性劳动的前提下，还可以根据这些附图获取到其他附图。
[0055]
图1为本发明实施例提供的一种自旋空间目标姿态估计方法的步骤流程图；
[0056]
图2为本发明实施例提供的仿真实验采用的空间目标结构与运动模型分析图；
[0057]
图3为本发明实施例提供的仿真实验提取的目标投影长度特征结构的示意图；
[0058]
图4为本发明实施例提供的一种自旋空间目标姿态估计系统的结构框图；
[0059]
图5为本发明实施例提供的一种自旋空间目标姿态估计装置的结构框图。
具体实施方式
[0060]
下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，仅用于解释本发明，而不能理解为对本发明的限制。对于以下实施例中的步骤编号，其仅为了便于阐述说明而设置，对步骤之间的顺序不做任何限定，实施例中的各步骤的执行顺序均可根据本领域技术人员的理解来进行适应性调整。
[0061]
在本发明的描述中，多个的含义是两个或两个以上，如果有描述到第一、第二只是用于区分技术特征为目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量或者隐含指明所指示的技术特征的先后关系。此外，除非另有定义，本文所使
用的所有的技术和科学术语与本技术领域的技术人员通常理解的含义相同。
[0062]
参照图1，本发明实施例提供了一种自旋空间目标姿态估计方法，具体包括以下步骤：
[0063]
s101、通过多个isar雷达系统获取预设时间段内的多个第一回波信号，并对第一回波信号进行距离
‑
多普勒处理，得到目标区域的第一雷达图像序列，第一雷达图像序列包括多个isar观测图像。
[0064]
具体地，本发明实施例由3站isar雷达系统接收预设时间段内的第一回波信号，然后将3站isar雷达系统接收到的第一回波信号同步后进行距离
‑
多普勒处理，得到目标区域的第一雷达图像序列，该序列包括3幅目标区域的isar观测图像。
[0065]
s102、在各isar观测图像中均提取出两个目标投影长度特征结构，并确定目标投影长度特征结构在距离维度上的目标距离维度尺寸和在多普勒维度上的目标多普勒维度尺寸。
[0066]
具体地，在这3帧isar观测图像内分别提取2个目标投影长度特征结构，并记录每个特征结构在距离、多普勒两个维度上的投影长度，即为目标距离维度尺寸和目标多普勒维度尺寸。
[0067]
s103、获取各isar观测图像对应的第一俯仰角和第一方位角，并根据第一俯仰角和第一方位角确定各isar雷达系统的中心雷达视线方向向量。
[0068]
具体地，根据各帧isar观测图像对应的雷达观测俯仰角、方位角信息，计算对应的isar雷达系统的中心雷达视线方向向量。
[0069]
进一步作为可选的实施方式，中心雷达视线方向向量的计算公式为：
[0070][0071]
其中，表示中心雷达视线方向向量，θ(t
m
)表示第一俯仰角，φ(t
m
)表示第一方位角，t
m
表示isar观测图像的相干积累时间的中心时刻。
[0072]
具体地，第一俯仰角是雷达视线与本体坐标系xoy平面的夹角，第一方位角是雷达视线在xoy平面的投影与y轴的夹角。
[0073]
s104、根据目标距离维度尺寸和中心雷达视线方向向量构建第一优化模型，并利用粒子群算法对第一优化模型进行求解得到目标瞬时姿态参数。
[0074]
进一步作为可选的实施方式，第一优化模型包括第一目标函数和第一约束条件，第一目标函数为：
[0075][0076]
其中，n∈[1，2]，p∈[1，p]，p表示isar雷达系统的数量，α
n
和β
n
表示待求解的目标瞬时姿态参数，r
n，p
表示第n个目标投影长度特征结构在第p个isar观测图像的目标距离维度尺寸，l
n
表示第n个目标投影长度特征结构的真实尺寸，表示第p个isar观测图像对应的isar雷达系统的中心雷达视线方向向量；
[0077]
第一约束条件为：
[0078]
tanα1tanα2+cos(β1‑
β2)＝0
[0079]
其中，α1、α2、β1和β2表示待求解的目标瞬时姿态参数。
[0080]
具体地，本发明实施例中p＝3，α
n
表示第n个目标投影长度特征与目标本体坐标系xoy平面的夹角，β
n
表示第n个目标投影长度特征在xoy平面的投影与y轴的夹角。可以认识到，通过第一优化模型求解得到α1、α1、β1和β2后，即可确定自旋空间目标的瞬时姿态。
[0081]
进一步作为可选的实施方式，利用粒子群算法对第一优化模型进行求解得到目标瞬时姿态参数这一步骤，其具体包括：
[0082]
设置最短移动距离，并在第一优化模型可行域内随机生成第一粒子群，第一粒子群的可行解为第一个体位置；
[0083]
根据第一代价函数搜索得到当前的粒子最优位置和粒子群最优位置；
[0084]
调整第一粒子群内粒子位置后，更新当前的粒子最优位置和粒子群最优位置，直至满足最大迭代步数或最小变化准则；
[0085]
根据当前最优的第一个体位置确定目标瞬时姿态参数。
[0086]
具体地，利用粒子群算法对第一优化模型进行求解得到目标瞬时姿态参数的步骤如下：
[0087]
s1041、设置最短移动距离，在第一优化模型可行域内随机生成第一粒子群，其可行解被定义为第一个体位置x
i
＝(α1，β1，α2，β2)
t
；
[0088]
s1042、根据下列第一代价函数搜索当前粒子最优位置pbest和粒子群最优位置gbest：
[0089][0090]
其中，置信因子a1决定先验条件在优化求解中的约束强弱，在发明实施例中a1设置为0.93：
[0091]
s1043、根据如下更新公式，调整第一粒子群内粒子位置：
[0092]
v
i
(t+1)＝a2v
i
(t)+a3rand1(pbest
‑
x
i
(t))+a4rand2(gbest
‑
x
i
(t))
[0093]
x
i
(t+1)＝x
i
(t)+v
i
(t)
[0094]
其中，v
i
(t)和x
i
(t)为第i次迭代中第t个个体的运动及位置，a2为惯性因子，a3和a4是用来平衡个体经验及群体经验的权重因子，随机参数rand1和rand2服从[0，1]均匀分布。
[0095]
s1044、根据第一代价函数重新计算个体目标函数值，更新当前粒子最优位置和粒子群最优位置。如果满足最大迭代步数或最小变化准则，停止迭代跳转步骤s1045；否则，跳转步骤s1043；最小变化准则指的在更新过程中，粒子最优位置和粒子群最优位置的变化需要大于最小移动阈值。
[0096]
s1045、停止迭代，输出当前最优的第一个体位置
[0097]
s105、根据目标多普勒维度尺寸和目标瞬时姿态参数构建第二优化模型，并利用粒子群算法对第二优化模型进行求解得到目标自旋参数，进而根据目标瞬时姿态参数和目标自旋参数完成对自旋空间目标的姿态估计。
[0098]
进一步作为可选的实施方式，第二优化模型包括第二目标函数，第二目标函数为：
[0099][0100]
其中，其中，n∈[1，2]，p∈[1，p]，p表示isar雷达系统的数量，θ
rot
、φ
rot
和ω
rot
表示待求解的目标自旋参数，表示预设的雷达视线参数，表示第p个isar观测图像对应的isar雷达系统的中心雷达视线方向向量，d
n，p
表示第n个目标投影长度特征结构在第p个isar观测图像的目标多普勒维度尺寸。
[0101]
具体地，发明实施例中p＝3，θ
rot
表示自旋空间目标的旋转轴与目标本体坐标系xoy平面的夹角，φ
rot
表示自旋空间目标的旋转轴在xoy平面的投影与y轴的夹角，ω
rot
表示自旋空间目标的平均转速。可以认识到，通过第二优化模型求解得到θ
rot
、φ
rot
和ω
rot
后，即可确定自旋空间目标的自旋姿态。
[0102]
进一步作为可选的实施方式，利用粒子群算法对第二优化模型进行求解得到目标自旋参数这一步骤，其具体包括：
[0103]
设置最短移动距离，并在第二优化模型可行域内随机生成第二粒子群，第二粒子群的可行解为第二个体位置；
[0104]
根据第二代价函数搜索得到当前的粒子最优位置和粒子群最优位置；
[0105]
调整第二粒子群内粒子位置后，更新当前的粒子最优位置和粒子群最优位置，直至满足最大迭代步数或最小变化准则；
[0106]
根据当前最优的第二个体位置确定目标自旋参数。
[0107]
具体地，利用粒子群算法对第二优化模型进行求解得到目标自旋参数的步骤如下：
[0108]
s1051、设置最短移动距离，在第二优化模型可行域内随机生成第二粒子群，其可行解被定义为第二个体位置x
i
＝(θ
rot
，φ
rot
，ω
rot
)
t
；
[0109]
s1052、根据以下第二代价函数搜索当前粒子最优位置pbest和粒子群最优位置gbest：
[0110][0111]
s1053、根据如下更新公式，调整第二粒子群内粒子位置：
[0112]
v
i
(t+1)＝a5v
i
(t)+a6rand1(pbest
‑
x
i
(t))+a7rand2(gbest
‑
x
i
(t))
[0113]
x
i
(t+1)＝x
i
(t)+v
i
(t)
[0114]
其中，v
i
(t)和x
i
(t)为第i次迭代中第t个个体的运动及位置，a5为惯性因子，a6和a7是用来平衡个体经验及群体经验的权重因子，随机参数rand1和rand2服从[0，1]均匀分布。
[0115]
s1054、根据第二代价函数重新计算个体目标函数值，更新当前粒子最优位置和粒子群最优位置，如果满足最大迭代步数或最小变化准则，停止迭代跳转步骤s1055；否则，跳
转步骤s1054；最小变化准则指的在更新过程中，粒子最优位置和粒子群最优位置的变化需要大于最小移动阈值。
[0116]
s1055、停止迭代，输出当前最优的第二个体位置x
i
＝(θ
rot
，φ
rot
，ω
rot
)
t
。
[0117]
进一步作为可选的实施方式，根据以下公式完成对自旋空间目标的姿态估计：
[0118][0119][0120][0121]
其中，和表示自旋空间目标的特征姿态，α1、α2、β1和β2表示目标瞬时姿态参数，表示自旋空间目标的旋转向量，θ
rot
、φ
rot
和ω
rot
表示目标自旋参数。
[0122]
本发明实施例对预设时间段内的多站观测得到的雷达回波信号同步后进行距离
‑
多普勒处理，获得多帧同步的isar观测图像，然后在获得的多帧isar观测图像内提取目标在距离向、多普勒向上的投影尺寸参数，结合各isar雷达系统的观测几何参数建立两个独立的优化模型，使用粒子群优化方法依次求解自旋空间目标的瞬时姿态参数以及自旋参数，完成对自旋空间目标的姿态估计。相较于现有技术，本发明实施例具有以下优点：
[0123]
1)本发明实施例利用自旋空间目标的多站isar同步图像的相对投影变化关系，结合空间投影理论，通过二维isar图像序列对空间目标的在轨瞬时状态实现精确反演，本发明具有很高的实用性。
[0124]
2)本发明实施例利用站联合观测对应的观测角度资源弥补了现有单站序列成像解译方式在求解高阶目标动态参数优化方面的不足，避免了某些观测视角下目标投影特征起伏，提升了本发明在实际合作、非合作空间目标状态估计应用中的鲁棒性。
[0125]
为进一步验证本发明实施例的准确性，下面结合仿真实验对本发明实施例的效果作进一步说明。
[0126]
本发明实施例的仿真实验采用的空间目标的结构如图2所示，其中θ表示雷达观测俯仰角，也即雷达视线与本体坐标系xoy平面的夹角，φ表示雷达观测方位角，也即雷达视线在xoy平面的投影与y轴的夹角，los
center
表示中心雷达视线方向，表示中心雷达视线方向向量，该目标的太阳能翼边界有明显的线性结构。本发明实施例的仿真实验的主要参数如下表1所示。
[0127]
[0128][0129]
表1
[0130]
仿真实验1：采用本发明的方法对图2的空间目标的3幅同步isar观测图像进行目标投影长度特征结构的提取，图3(a)、图3(b)和图3(c)示出了这3幅isar观测图像提取出的目标投影长度特征结构，分别用特征1和特征2表示两个目标投影长度特征结构。然后，根据获得的3站图像目标距离维度尺寸、目标多普勒维度尺寸以及三站中心雷达视线方向向量构建的两个优化模型，并使用粒子群算法求解目标各状态参数，结果如下表2所示。
[0131]
真实值估计值误差特征1姿态(0.2387，
‑
0.2002，
‑
0.9502)(
‑
0.2620，0.2198，0.9397)1.87度特征2姿态(0.6505，0.7595，0.0034)(0.6428，0.7660，
‑
0.0000)0.61度自旋轴指向(1.0000，0.0000，0.0000)(0.9995，0.0298，0.0094)1.81度转速0.0150rad/s0.0151rad/s0.0001rad/s
[0132]
表2
[0133]
从图3(a)、图3(b)和图3(c)中可以看出，空间目标的特征结构的稳定提取基本可以得到保证。从表2中可以看出，估计的特征结构的指向与真实的特征结构的指向基本一致，平均误差在3度之内，估计的目标自旋方向与真实的目标自旋方向基本一致，估计的目标自旋速度与真实的目标自旋速度在数值上较为接近，则空间目标的在轨状态可以确定。
[0134]
仿真实验2：采用本发明方法在目标进动状态下进行测试，即目标在本体坐标系内绕着某一固定轴做匀速转动。目标自转速度被设定为0.015rad/s，转轴转动速度被设定为0rad/s。根据60秒连续观测估计目标6个瞬时状态参数，目标状态估计结果如下表3所示。
[0135][0136]
表3
[0137]
从表3中可以看出，在目标在轨进动状态下，估计的特征结构的指向与真实的特征结构的指向基本一致，平均误差在3度之内，估计的目标自旋方向与真实的目标自旋方向基本一致，估计的目标自旋速度与真实的目标自旋速度在数值上较为接近，则空间目标的在轨状态可以确定。
[0138]
仿真实验3：采用本发明方法在目标章动状态下进行测试，即目标除基本进动外其自转轴也绕某一固定轴进行匀速转动。目标自转速度被设定为0.015rad/s，转轴转动速度被设定为0.005rad/s。根据60秒连续观测估计目标6个瞬时状态参数，目标状态估计结果如下表4所示。
[0139][0140]
表4
[0141]
从表4中可以看出，在目标在轨章动状态下，估计的特征结构的指向与真实的特征结构的指向基本一致，平均误差在3度之内，估计的目标自旋方向与真实的目标自旋方向基本一致，估计的目标自旋速度与真实的目标自旋速度在数值上较为接近，则空间目标的在轨状态可以确定。
[0142]
仿真实验4：采用本发明方法在目标翻滚状态下进行测试，即目标自转轴及自转速度不断变化。根据60秒连续观测估计目标6个瞬时状态参数，目标状态估计结果如下表5所示。
[0143][0144]
表5
[0145]
从表5中可以看出，在目标翻滚状态下，估计的特征结构的指向与真实的特征结构的指向基本一致，平均误差在3度之内，估计的目标自旋方向与真实的目标自旋方向基本一致，估计的目标自旋速度与真实的目标自旋速度在数值上较为接近，则空间目标的在轨状态可以确定。
[0146]
参照图4，本发明实施例提供了一种自旋空间目标姿态估计系统，包括：
[0147]
第一雷达图像序列获取模块，用于通过多个isar雷达系统获取预设时间段内的多个第一回波信号，并对第一回波信号进行距离
‑
多普勒处理，得到目标区域的第一雷达图像序列，第一雷达图像序列包括多个isar观测图像；
[0148]
目标距离维度尺寸和目标多普勒维度尺寸确定模块，用于在各isar观测图像中均提取出两个目标投影长度特征结构，并确定目标投影长度特征结构在距离维度上的目标距离维度尺寸和在多普勒维度上的目标多普勒维度尺寸；
[0149]
中心雷达视线方向向量确定模块，用于获取各isar观测图像对应的第一俯仰角和
第一方位角，并根据第一俯仰角和第一方位角确定各isar雷达系统的中心雷达视线方向向量；
[0150]
目标瞬时姿态参数确定模块，用于根据目标距离维度尺寸和中心雷达视线方向向量构建第一优化模型，并利用粒子群算法对第一优化模型进行求解得到目标瞬时姿态参数；
[0151]
目标自旋参数确定模块，用于根据目标多普勒维度尺寸和目标瞬时姿态参数构建第二优化模型，并利用粒子群算法对第二优化模型进行求解得到目标自旋参数，进而根据目标瞬时姿态参数和目标自旋参数完成对自旋空间目标的姿态估计。
[0152]
上述方法实施例中的内容均适用于本系统实施例中，本系统实施例所具体实现的功能与上述方法实施例相同，并且达到的有益效果与上述方法实施例所达到的有益效果也相同。
[0153]
参照图5，本发明实施例提供了一种自旋空间目标姿态估计装置，包括：
[0154]
至少一个处理器；
[0155]
至少一个存储器，用于存储至少一个程序；
[0156]
当上述至少一个程序被上述至少一个处理器执行时，使得上述至少一个处理器实现上述的一种自旋空间目标姿态估计方法。
[0157]
上述方法实施例中的内容均适用于本装置实施例中，本装置实施例所具体实现的功能与上述方法实施例相同，并且达到的有益效果与上述方法实施例所达到的有益效果也相同。
[0158]
本发明实施例还提供了一种计算机可读存储介质，其中存储有处理器可执行的程序，该处理器可执行的程序在由处理器执行时用于执行上述一种自旋空间目标姿态估计方法。
[0159]
本发明实施例的一种计算机可读存储介质，可执行本发明方法实施例所提供的一种自旋空间目标姿态估计方法，可执行方法实施例的任意组合实施步骤，具备该方法相应的功能和有益效果。
[0160]
本发明实施例还公开了一种计算机程序产品或计算机程序，该计算机程序产品或计算机程序包括计算机指令，该计算机指令存储在计算机可读存介质中。计算机设备的处理器可以从计算机可读存储介质读取该计算机指令，处理器执行该计算机指令，使得该计算机设备执行图1所示的方法。
[0161]
在一些可选择的实施例中，在方框图中提到的功能/操作可以不按照操作示图提到的顺序发生。例如，取决于所涉及的功能/操作，连续示出的两个方框实际上可以被大体上同时地执行或上述方框有时能以相反顺序被执行。此外，在本发明的流程图中所呈现和描述的实施例以示例的方式被提供，目的在于提供对技术更全面的理解。所公开的方法不限于本文所呈现的操作和逻辑流程。可选择的实施例是可预期的，其中各种操作的顺序被改变以及其中被描述为较大操作的一部分的子操作被独立地执行。
[0162]
此外，虽然在功能性模块的背景下描述了本发明，但应当理解的是，除非另有相反说明，上述的功能和/或特征中的一个或多个可以被集成在单个物理装置和/或软件模块中，或者一个或多个功能和/或特征可以在单独的物理装置或软件模块中被实现。还可以理解的是，有关每个模块的实际实现的详细讨论对于理解本发明是不必要的。更确切地说，考
虑到在本文中公开的装置中各种功能模块的属性、功能和内部关系的情况下，在工程师的常规技术内将会了解该模块的实际实现。因此，本领域技术人员运用普通技术就能够在无需过度试验的情况下实现在权利要求书中所阐明的本发明。还可以理解的是，所公开的特定概念仅仅是说明性的，并不意在限制本发明的范围，本发明的范围由所附权利要求书及其等同方案的全部范围来决定。
[0163]
上述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例上述方法的全部或部分步骤。而前述的存储介质包括：u盘、移动硬盘、只读存储器(rom，read
‑
only memory)、随机存取存储器(ram，random access memory)、磁碟或者光盘等各种可以存储程序代码的介质。
[0164]
在流程图中表示或在此以其他方式描述的逻辑和/或步骤，例如，可以被认为是用于实现逻辑功能的可执行指令的定序列表，可以具体实现在任何计算机可读介质中，以供指令执行系统、装置或设备(如基于计算机的系统、包括处理器的系统或其他可以从指令执行系统、装置或设备取指令并执行指令的系统)使用，或结合这些指令执行系统、装置或设备而使用。就本说明书而言，“计算机可读介质”可以是任何可以包含、存储、通信、传播或传输程序以供指令执行系统、装置或设备或结合这些指令执行系统、装置或设备而使用的装置。
[0165]
计算机可读介质的更具体的示例(非穷尽性列表)包括以下：具有一个或多个布线的电连接部(电子装置)，便携式计算机盘盒(磁装置)，随机存取存储器(ram)，只读存储器(rom)，可擦除可编辑只读存储器(eprom或闪速存储器)，光纤装置，以及便携式光盘只读存储器(cdrom)。另外，计算机可读介质甚至可以是可在其上打印上述程序的纸或其他合适的介质，因为可以例如通过对纸或其他介质进行光学扫描，接着进行编辑、解译或必要时以其他合适方式进行处理来以电子方式获得上述程序，然后将其存储在计算机存储器中。
[0166]
应当理解，本发明的各部分可以用硬件、软件、固件或它们的组合来实现。在上述实施方式中，多个步骤或方法可以用存储在存储器中且由合适的指令执行系统执行的软件或固件来实现。例如，如果用硬件来实现，和在另一实施方式中一样，可用本领域公知的下列技术中的任一项或他们的组合来实现：具有用于对数据信号实现逻辑功能的逻辑门电路的离散逻辑电路，具有合适的组合逻辑门电路的专用集成电路，可编程门阵列(pga)，现场可编程门阵列(fpga)等。
[0167]
在本说明书的上述描述中，参考术语“一个实施方式/实施例”、“另一实施方式/实施例”或“某些实施方式/实施例”等的描述意指结合实施方式或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施方式或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施方式或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施方式或示例中以合适的方式结合。
[0168]
尽管已经示出和描述了本发明的实施方式，本领域的普通技术人员可以理解：在不脱离本发明的原理和宗旨的情况下可以对这些实施方式进行多种变化、修改、替换和变型，本发明的范围由权利要求及其等同物限定。
[0169]
以上是对本发明的较佳实施进行了具体说明，但本发明并不限于上述实施例，熟悉本领域的技术人员在不违背本发明精神的前提下还可做作出种种的等同变形或替换，这些等同的变形或替换均包含在本技术权利要求所限定的范围内。