一种基于物理先验知识的空间群目标探测方法及系统

1.本发明属于信号处理技术领域，特别是一种基于物理先验知识的空间群目标探测方法及系统。


背景技术：

2.卷积神经网络作为当下最热门的基于深度学习理论的图像内容学习的最佳算法，已在雷达目标特征提取领域中得到广泛应用，实现了从目标时频图快速获得目标本征参数（wang s et al.，“cone-shaped space target inertia characteristics identification by deep learning with compressed dataset.”ieee transactions on antennas and propagation, 2022, 70(7): 5217
–
5226.），但是卷积神经网络仅作为黑盒子使用，缺乏对目标物理先验知识的利用。
3.专利cn110411459a公布了一种使用图像序列的空间非合作目标自旋参数的辨识方法，解决了先验信息不足时的姿态估计问题，同时获得目标的自旋参数，但该方法需要提前知道目标尺寸长宽比，并且仅针对单个目标的情况，不能处理目标之间的遮挡交叉等情况，不适用于空间群目标的参数辨识。
4.针对一个雷达波束内有多个目标的情况，从而导致无法分辨出每一个个体目标，提出了一种基于多假设跟踪算法的航迹关联方法（r. du, l. liu, x. bai, and f. zhou,
ꢀ“
a new scatterer trajectory association method for isar image sequence utilizing multiple hypothesis tracking algorithm,
”ꢀ
ieee trans. geosci. remote sensing, 2022, 60: 1
–
13），但椭圆轨迹运动假设在一定程度上限制了该算法在测量数据中的应用，需要对所提出的算法进行进一步的优化研究。


技术实现要素：

5.本发明的目的在于提供一种基于物理先验知识的空间群目标探测方法及系统，通过对复杂环境中空间群目标进行分离重构，来实现复杂电磁环境下，对雷达目标本征特征的精确提取。
6.实现本发明目的的技术解决方案为：一种基于物理先验知识的空间群目标探测方法，包括以下步骤：
7.步骤1、构建不同尺寸的空间非合作锥体目标物理仿真模型，获得不同运动状态下混叠的空间群目标随时间变化的雷达回波，利用傅里叶变换获得空间群目标随时间变化的一维距离像；
8.步骤2、基于锥体目标散射中心模型的物理先验知识，根据多个随机时刻的一维距离像估算出空间群目标散射中心的数量，结合最小二乘法，提取出所有时刻散射中心的位置信息和幅度信息；
9.步骤3、基于卡尔曼滤波算法，将混叠的锥体目标散射中心模型中的位置信息重组，减少目标之间的遮挡和交叉干扰，实现空间群目标的分离，获得单个锥体目标散射中心
模型，并重构单个锥体目标微多普勒时频图和一维距离像；
10.步骤4、利用重构获得的单个锥体目标微多普勒时频图数据集和一维距离像数据集，结合双通道卷积神经网络进行训练和测试，得到复杂环境中空间非合作锥体目标的高度、半径、进动频率。
11.一种基于物理先验知识的空间群目标探测系统，该系统用于实现所述的基于物理先验知识的空间群目标探测方法，具体包括一维距离像获取模块、散射中心信息提取模块、单个锥体目标重构模块、空间非合作锥体目标特征提取模块，其中：
12.一维距离像获取模块，用于构建不同尺寸的空间非合作锥体目标物理仿真模型，获得不同运动状态下混叠的空间群目标随时间变化的雷达回波，利用傅里叶变换获得空间群目标随时间变化的一维距离像；
13.散射中心信息提取模块，基于锥体目标散射中心模型的物理先验知识，根据多个随机时刻的一维距离像估算出空间群目标散射中心的数量，结合最小二乘法，提取出所有时刻散射中心的位置信息和幅度信息；
14.单个锥体目标重构模块，基于卡尔曼滤波算法，将混叠的锥体目标散射中心模型中的位置信息重组，减少目标之间的遮挡和交叉干扰，实现空间群目标的分离，获得单个锥体目标散射中心模型，并重构单个锥体目标微多普勒时频图和一维距离像；
15.空间非合作锥体目标特征提取模块，利用重构获得的单个锥体目标微多普勒时频图数据集和一维距离像数据集，结合双通道卷积神经网络进行训练和测试，得到复杂环境中空间非合作锥体目标的高度、半径、进动频率。
16.一种移动终端，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现所述的基于物理先验知识的空间群目标探测方法。
17.本发明与现有技术相比，其显著优点为：（1）通过散射中心模型分离空间群目标，获得单个锥体目标散射中心模型，并重构单个锥体目标微多普勒时频图和一维距离像，能够有效减少噪声的干扰，对于提高目标运动参数和结构参数的估计精度有较大的帮助（2）利用卡尔曼滤波算法，提取空间群目标散射中心的运动轨迹，避免了目标之间的交叉和遮挡，为准确辨识目标运动参数和结构参数提供帮助（3）对多种复杂环境中的空间非合作锥体目标一维距离像均有较好的分离重构能力（4）通过深度学习来学习空间非合作锥体目标一维距离像，网络训练和测试过程简单且快速，模型的泛化能力较强。
附图说明
18.图1是本发明中典型的空间非合作锥体目标进动示意图。
19.图2是本发明中典型的三锥体目标空间运动过程示意图。
20.图3是基于一维距离像的散射中心模型提取分离流程图。
21.图4是空间群目标散射中心模型与单个锥体目标散射中心模型示意图。
22.图5是重构的单个锥体目标一维距离像和微多普勒时频图。
23.图6是双通道卷积神经网络识别目标参数流程图。
24.图7是参数识别误差变化图。
具体实施方式
25.本发明提供一种基于物理先验知识的空间群目标探测方法，包括以下步骤：
26.步骤1、构建不同尺寸的空间非合作锥体目标物理仿真模型，获得不同运动状态下混叠的空间群目标随时间变化的雷达回波，利用傅里叶变换获得空间群目标随时间变化的一维距离像；
27.步骤2、基于锥体目标散射中心模型的物理先验知识，根据多个随机时刻的一维距离像估算出空间群目标散射中心的数量，结合最小二乘法，提取出所有时刻散射中心的位置信息和幅度信息；
28.步骤3、基于卡尔曼滤波算法，将混叠的锥体目标散射中心模型中的位置信息重组，减少目标之间的遮挡和交叉干扰，实现空间群目标的分离，获得单个锥体目标散射中心模型，并重构单个锥体目标微多普勒时频图和一维距离像；
29.步骤4、利用重构获得的单个锥体目标微多普勒时频图数据集和一维距离像数据集，结合双通道卷积神经网络进行训练和测试，得到复杂环境中空间非合作锥体目标的高度、半径、进动频率。
30.作为一种具体示例，步骤1中所述的构建不同尺寸的空间非合作锥体目标物理仿真模型，获得不同运动状态下混叠的空间群目标随时间变化的雷达回波，利用傅里叶变换获得空间群目标随时间变化的一维距离像，具体如下：
31.对于旋转对称锥体目标，散射中心分布在两个位置，分别为锥体目标顶端和锥体目标底部边缘；
32.锥体目标的运动参数和结构参数为：进动频率，速度，进动角，高度，半径，质心高；定义进动角频率为，为当前时间，为雷达俯仰角，为雷达方位角；
33.在以目标质心为原点的坐标系下，两个位置处的散射中心径向位移理论表达式为：
34.锥体目标顶端散射中心a的径向位移表达式为：
35.(1)
36.锥体目标底部边缘散射中心b的径向位移表达式为：
37.(2)
38.理想散射中心模型表达式为：
39.(3)
40.其中为第i个散射中心的幅度信息，为第i个散射中心的位置信息，为雷达发射波长，为雷达回波；
41.由式(3)进一步推导出混叠的锥体目标雷达回波表达式为：
42.(4)
43.其中、为第i个锥体目标锥顶和锥底的幅度信息，、为第i个锥体目标锥顶和锥底的位置信息。
44.作为一种具体示例，步骤2所述的基于锥体目标散射中心模型的物理先验知识，根据多个随机时刻的一维距离像估算出空间群目标散射中心的数量，结合最小二乘法，提取出所有时刻散射中心的位置信息和幅度信息，具体如下：
45.由锥体目标的顶端散射中心a和底部边缘散射中心b的径向位移表达式(1)、式(2)获得锥体目标散射中心模型的物理先验知识，即噪声环境下的锥体目标随时间变化的散射中心模型为：
[0046] (5)
[0047]
其中，为t时刻的雷达回波；表示t时刻噪声，、为第i个锥体目标锥顶、锥底t时刻的幅度信息，、为第i个锥体目标锥顶、锥底t时刻的位置信息，为第i个锥体目标的速度，为第i个锥体目标的进动角，为第i个锥体目标的高度，为第i个锥体目标的半径，为第i个锥体目标的质心高，为第i个锥体目标的角速度；
[0048]
对噪声环境下的锥体目标随时间变化的散射中心模型进行变化，得到t时刻一维距离像的噪声：
[0049]
(6)
[0050]
其中为t时刻一维距离像的噪声；
[0051]
通过多个随机时刻一维距离像的时间累积，获得增强散射中心的一维距离像，提高真实散射中心的强度，减小噪声干扰；基于最小二乘法和锥体目标散射中心模型物理先验知识，提取一维距离像的散射中心数量，以及散射中心的位置信息和幅度信息，代入式(6)获得t时刻一维距离像的噪声，将多个随机时刻一维距离像的噪声最大值的均值作为最终的噪声阈值；
[0052]
根据获得的空间群目标散射中心的数量和噪声阈值，基于最小二乘法和锥体目标散射中心模型物理先验知识，对所有时刻的空间群目标一维距离像提取散射中心的位置信息和幅度信息，获得空间群目标散射中心模型。
[0053]
作为一种具体示例，空间群目标一维距离像的散射中心提取具体流程如下：
[0054]
(2.1)随机选取多段连续时刻的一维距离像，将每段连续时刻的一维距离像做时间维度上的累加，获得增强散射中心的一维距离像，提高真实散射中心的强度，减小噪声干扰；
[0055]
(2.2)将增强散射中心的一维距离像的均值作为起始噪声阈值，幅值大于噪声阈值的视为散射中心，幅值低于噪声阈值的视为噪声，基于最小二乘法和锥体目标散射中心模型物理先验知识，获得步骤(2.1)中多段一维距离像的散射中心数量，以及散射中心的位置信息和幅度信息；
[0056]
(2.3)将步骤(2.2)中获得的散射中心数量中的众数作为空间群目标散射中心的数量上限m0，并且锥体目标数量为m0/2，将步骤(2.2)中获得的散射中心的位置信息和幅度信息代入式(6)，获得t时刻一维距离像的噪声，将各个时刻一维距离像的噪声最大值的均值作为最终的噪声阈值n0；
[0057]
(2.4)选取t时刻的一维距离像l，设置初始散射中心个数n=0；
[0058]
(2.5)基于最小二乘法，提取t时刻的一维距离像l中幅度最大的散射中心的位置信息和幅度信息，将获得的散射中心的位置信息和幅度信息代入式(4)，重新生成无噪声一维距离像，与原本的一维距离像相减获得提取后的一维距离像l’，且散射中心个数n=n+1；
[0059]
(2.6)判断t时刻的散射中心个数是否大于等于m0或者一维距离像l’的最大值小于n0：如果判断为是，则表示t时刻一维距离像的散射中心提取完备，令t=t+1，并返回步骤(2.4)；如果判断为否，则表示t时刻一维距离像的散射中心未提取完备，令t时刻一维距离像l=l’，并返回步骤(2.5)。
[0060]
作为一种具体示例，步骤3所述的基于卡尔曼滤波算法，将混叠的锥体目标散射中心模型中的位置信息重组，减少目标之间的遮挡和交叉干扰，实现空间群目标的分离，获得单个锥体目标散射中心模型，并重构单个锥体目标微多普勒时频图和一维距离像，具体如下：
[0061]
通过步骤2获得的空间群目标散射中心模型，基于已获得的散射中心位置信息，结合卡尔曼滤波算法和锥体目标散射中心运动轨迹为正弦函数的性质，计算锥体目标散射中心运动轨迹相关参数，获得锥体目标散射中心的运动轨迹，将存在交叉遮挡的多个散射中
心运动轨迹分离出来，锥体目标散射中心运动轨迹表达式如下所示：
[0062]
(7)
[0063]
其中，、分别为第i个锥体目标的锥顶和锥底运动轨迹，为第i个锥体目标的运动幅值，为第i个锥体目标的运动角速度，为第i个锥体目标的初始相位，、分别为第i个锥体目标的锥顶和锥底初始位置，为第i个锥体目标的径向速度。
[0064]
作为一种具体示例，分离获得单个锥体目标散射中心模型，并重构单个锥体目标微多普勒时频图和一维距离像，具体流程如下：
[0065]
(3.1)排除步骤2中提取出的散射中心数量小于m0的时刻，即散射中心被遮挡时刻，并计算时刻，所有散射中心位置信息和幅度信息相关系数，选取相关系数最大，即不存在散射中心运动轨迹交叉的时刻作为为起始时刻t，；
[0066]
(3.2)根据t时刻的第i个锥体目标散射中心的位置和锥体目标散射中心运动轨迹公式(7)，结合卡尔曼滤波算法，预测t+1时刻的第i个锥体目标散射中心的位置；
[0067]
(3.3)计算t+1时刻的预测值与t+1时刻的第j个散射中心的位置信息的误差，判断是否存在大于阈值，该阈值为锥体散射中心运动幅值的5%：如果判断为是，则结合最小二乘法，根据已知的0~t时刻的第i个锥体目标散射中心的运动轨迹，以及锥体目标散射中心运动轨迹表达式(7)，获得第i个锥体目标的运动参数，并重复步骤(3.2)~步骤(3.3)；如果判断为否，则将预测值作为第i个锥体目标散射中心t+1时刻散射中心的位置，并且t=t+1，重复步骤(3.2)~步骤(3.3)；当遍历完所有时刻后，进入步骤(3.4)；
[0068]
(3.4)当遍历完所有时刻后，根据最终获得的第i个散射中心的运动轨迹，结合最小二乘法，获得锥体目标散射中心运动轨迹的相关参数、、、、、，通过公式(7)补全步骤(3.1)中排除的被遮挡时刻散射中心的位置，并由插值获得对应时刻的散射中心的幅度，最终获得完整的第i个锥体目标的锥顶散射中心的位置信息与幅度信息、锥底散射中心的位置信息与幅度信息，即获得单个锥体目标的散射中心模型；
[0069]
(3.5)将获得的单个锥体目标的散射中心模型，代入公式(4)获得单个锥体目标的雷达回波，通过时频分析和傅里叶变换获得单个锥体目标的微多普勒时频图和一维距离像。
[0070]
本发明还提供一种基于物理先验知识的空间群目标探测系统，该系统用于实现所述的基于物理先验知识的空间群目标探测方法，具体包括一维距离像获取模块、散射中心信息提取模块、单个锥体目标重构模块、空间非合作锥体目标特征提取模块，其中：
[0071]
一维距离像获取模块，用于构建不同尺寸的空间非合作锥体目标物理仿真模型，获得不同运动状态下混叠的空间群目标随时间变化的雷达回波，利用傅里叶变换获得空间群目标随时间变化的一维距离像；
[0072]
散射中心信息提取模块，基于锥体目标散射中心模型的物理先验知识，根据多个随机时刻的一维距离像估算出空间群目标散射中心的数量，结合最小二乘法，提取出所有时刻散射中心的位置信息和幅度信息；
[0073]
单个锥体目标重构模块，基于卡尔曼滤波算法，将混叠的锥体目标散射中心模型中的位置信息重组，减少目标之间的遮挡和交叉干扰，实现空间群目标的分离，获得单个锥体目标散射中心模型，并重构单个锥体目标微多普勒时频图和一维距离像；
[0074]
空间非合作锥体目标特征提取模块，利用重构获得的单个锥体目标微多普勒时频图数据集和一维距离像数据集，结合双通道卷积神经网络进行训练和测试，得到复杂环境中空间非合作锥体目标的高度、半径、进动频率。
[0075]
本发明还提供一种移动终端，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现所述的基于物理先验知识的空间群目标探测方法。
[0076]
下面结合附图及具体实施例对本发明作进一步详细描述。
[0077]
实施例
[0078]
结合图1和图2，本实施例一种基于物理先验知识的空间群目标探测方法，步骤为：
[0079]
步骤1、构建不同尺寸的空间非合作锥体目标物理仿真模型，获得不同运动状态下,混叠的空间群目标随时间变化的雷达回波，利用傅里叶变换获得空间群目标随时间变化的一维距离像；
[0080]
本发明以旋转对称锥体目标为例，锥体目标进动示意图如图1所示，三锥体目标空间运动过程示意图如图2所示，对于旋转对称锥体目标，散射中心分布在两个位置，分别为锥体目标顶端和锥体目标底部边缘；
[0081]
锥体目标的运动参数和结构参数为：进动频率，速度，进动角，高度，半径，质心高；定义进动角频率为，为当前时间，为雷达俯仰角，为雷达方位角；
[0082]
在以目标质心为原点的坐标系下，两个位置处的散射中心径向位移理论表达式为：
[0083]
锥体目标顶端散射中心a的径向位移表达式为：
[0084][0085]
锥体目标底部边缘散射中心b的径向位移表达式为：
[0086][0087]
理想散射中心模型表达式为：
[0088][0089]
其中为第i个散射中心的幅度信息，为第i个散射中心的位置信息，为雷达发射波长，为雷达回波；
[0090]
进一步推导出混叠的锥体目标雷达回波表达式为：
[0091][0092]
其中、为第i个锥体目标锥顶、锥底的幅度信息，、为第i个锥体目标锥顶和锥底的位置信息。
[0093]
步骤2、基于锥体目标散射中心模型的物理先验知识，根据多个随机时刻的一维距离像估算出空间群目标散射中心的数量，结合最小二乘法，提取出所有时刻散射中心的位置信息和幅度信息；
[0094]
本发明是对噪声环境下，空间群目标一维距离像的散射中心进行提取，由锥体目标的顶端散射中心a和底部边缘散射中心b的径向位移表达式获得锥体目标散射中心模型的物理先验知识，即噪声环境下的锥体目标随时间变化的散射中心模型为：
[0095]
[0096]
其中表示t时刻噪声，、为第i个锥体目标锥顶和锥底t时刻的幅度信息，、为第i个锥体目标锥顶和锥底t时刻的位置信息，为第i个锥体目标的速度，为第i个锥体目标的进动角，为第i个锥体目标的高度，为第i个锥体目标的半径，质心高为第i个锥体目标的质心高，为第i个锥体目标的角速度；
[0097]
对噪声环境下的锥体目标随时间变化的散射中心模型进行变化，得到t时刻一维距离像的噪声：
[0098][0099]
其中为t时刻一维距离像的噪声；
[0100]
通过多个随机时刻一维距离像的时间累积，获得增强散射中心的一维距离像，提高真实散射中心的强度，减小噪声干扰；基于最小二乘法和锥体目标散射中心模型物理先验知识，提取一维距离像的散射中心数量，以及散射中心的位置信息和幅度信息，代入上式获得t时刻一维距离像的噪声，将多个随机时刻一维距离像的噪声最大值的均值作为最终的噪声阈值；
[0101]
根据获得的空间群目标散射中心的数量和噪声阈值，基于最小二乘法和锥体目标散射中心模型物理先验知识，对所有时刻的空间群目标一维距离像提取散射中心的位置信息和幅度信息，获得空间群目标散射中心模型；
[0102]
空间群目标一维距离像的散射中心模型提取具体流程如下：
[0103]
(1) 随机选取多段连续时刻的一维距离像，将每段连续时刻的一维距离像做时间维度上的累加，获得增强散射中心的一维距离像，提高真实散射中心的强度，减小噪声干扰；
[0104]
(2) 将增强散射中心的一维距离像的均值作为起始噪声阈值，幅值大于噪声阈值的视为散射中心，幅值低于噪声阈值的视为噪声，基于最小二乘法和锥体目标散射中心模型物理先验知识，获得步骤(1)中多段一维距离像的散射中心数量，以及散射中心的位置信息和幅度信息；
[0105]
(3) 将步骤(2)中获得的散射中心数量中的众数作为空间群目标散射中心的数量上限m0，并且锥体目标数量为m0/2，将步骤(2)中获得的散射中心的位置信息和幅度信息代入上式，获得t时刻一维距离像的噪声，将各个时刻一维距离像的噪声最大值的均值作为最终的噪声阈值n0；
[0106]
(4) 选取t时刻的一维距离像l，设置初始散射中心个数n=0；
[0107]
(5) 基于最小二乘法，提取t时刻的一维距离像l中幅度最大的散射中心的位置信息和幅度信息，将获得的散射中心的位置信息和幅度信息代入混
叠的锥体目标雷达回波表达式，重新生成无噪声一维距离像，与原本的一维距离像相减获得提取后的一维距离像l’，且散射中心个数n=n+1；
[0108]
(6) 判断t时刻的散射中心个数是否大于等于m0或者一维距离像l’的最大值小于n0，如果判断为是，则表示t时刻一维距离像的散射中心提取完备，令t=t+1，并返回步骤(4)，如果判断为否，则表示t时刻一维距离像的散射中心未提取完备，令t时刻一维距离像l=l’，并返回步骤(5)。
[0109]
步骤3、基于卡尔曼滤波算法，将混叠的锥体目标散射中心模型中的位置信息重组，减少目标之间的遮挡和交叉干扰，实现空间群目标的分离，获得单个锥体目标散射中心模型，并重构单个锥体目标微多普勒时频图和一维距离像；
[0110]
本发明是对存在交叉和遮挡的空间群目标进行分离，获得单个锥体目标的散射中心模型，并重构单个锥体目标的微多普勒时频图和一维距离像。通过步骤2获得的空间群目标散射中心模型，基于已获得的散射中心位置信息，结合卡尔曼滤波算法和锥体目标散射中心运动轨迹为正弦函数的性质，计算锥体目标散射中心运动轨迹相关参数，获得锥体目标散射中心的运动轨迹，将存在交叉遮挡的多个散射中心运动轨迹分离出来，锥体目标散射中心运动轨迹表达式如下所示：
[0111]
(21)
[0112]
其中、分别为第i个锥体目标的锥顶和锥底运动轨迹，为第i个锥体目标的运动幅值，为第i个锥体目标的运动角速度，为第i个锥体目标的初始相位，、分别为第i个锥体目标的锥顶和锥底初始位置，为第i个锥体目标的径向速度；
[0113]
结合图3说明，分离获得单个锥体目标散射中心模型，并重构单个锥体目标微多普勒时频图和一维距离像，具体流程如下：
[0114]
(1) 排除步骤2中提取出的散射中心数量小于m0的时刻，即散射中心被遮挡时刻，并计算时刻，所有散射中心位置信息和幅度信息相关系数，选取相关系数最大，即不存在散射中心运动轨迹交叉的时刻作为为起始时刻t，；
[0115]
(2) 根据t时刻的第i个锥体目标散射中心的位置和锥体目标散射中心运动轨迹公式，结合卡尔曼滤波算法，预测t+1时刻的第i个锥体目标散射中心的位置；
[0116]
(3) 计算t+1时刻的预测值与t+1时刻的第j个散射中心的位置信息的误差，判断是否存在大于阈值，阈值为锥体散射中心运动幅值的5%，如果判断为是，则结合最小二乘法，根据已知的0~t时刻的第i个锥体目标
散射中心的运动轨迹，以及锥体目标散射中心运动轨迹表达式，获得第i个锥体目标的运动参数，并重复步骤(2)(3)，如果判断为否，则将预测值作为第i个锥体目标散射中心t+1时刻散射中心的位置，并且t=t+1，重复步骤(2)(3)；
[0117]
(4) 当遍历完所有时刻后，根据最终获得的第i个散射中心的运动轨迹，结合最小二乘法，获得锥体目标散射中心运动轨迹的相关参数、、、、，通过锥体目标散射中心运动轨迹表达式补全步骤(1)中排除的被遮挡时刻散射中心的位置，并由插值获得对应时刻的散射中心的幅度，最终获得完整的第i个锥体目标的锥顶散射中心的位置信息与幅度信息和锥底散射中心的位置信息与幅度信息，即获得单个锥体目标的散射中心模型；
[0118]
(5) 将获得的单个锥体目标的散射中心模型，代入混叠的锥体目标雷达回波表达式获得单个锥体目标的雷达回波，通过时频分析和傅里叶变换获得单个锥体目标的微多普勒时频图和一维距离像。
[0119]
最终结果如图4和图5所示。
[0120]
步骤4、利用重构获得的单个锥体目标微多普勒时频图数据集和一维距离像数据集，结合双通道卷积神经网络进行训练和测试，得到复杂环境中空间非合作锥体目标的高度、半径、进动频率；
[0121]
结合图6说明双通道卷积神经网络识别参数流程。本发明利用单个锥体目标的散射中心模型，重构获得单个锥体目标一维距离像和微多普勒时频图，对双通道卷积神经网络进行训练测试，识别单个锥体目标参数，其中一个卷积通道提取一维距离像的特征，另一个卷积通道提取微多普勒时频图的特征，并通过1*1的卷积核将两种不同的散射属性特征融合，共同识别单个锥体目标参数。
[0122]
针对复杂环境中的空间群目标进动时，由于目标之间存在的交叉和遮挡等干扰，难以有效获得目标本征参数的问题，传统的目标辨识方法有许多局限性，例如要求目标之间无交叉或者只能针对单个目标进行辨识，不具有广泛的应用性。传统的运动参数提取方法也存在着提取效率慢，提取精度不高且泛化能力不强的问题。
[0123]
本发明的空间非合作锥体目标的参数设置为进动频率，间隔0.4hz，进动角，锥体高度，间隔0.2m, 锥体高度，间隔0.02m，锥体数量为3，来自复杂环境的噪声信噪比，间隔1db，以及无噪声，假定目标质量分布均匀，质点位于四分之一高处。根据不同尺寸的空间非合作锥体目标物理仿真模型的回波数据生成一维距离像数据集，并且数据集不重复。
[0124]
噪声环境下的锥体目标随时间变化的散射中心模型为：
[0125][0126]
其中表示t时刻噪声，、为第i个锥体目标锥顶和锥底t时刻的幅度信息，、为第i个锥体目标锥顶和锥底t时刻的位置信息，为第i个锥体目标的速度，为第i个锥体目标的进动角，为第i个锥体目标的高度，为第i个锥体目标的半径，质心高为第i个锥体目标的质心高，为第i个锥体目标的角速度；
[0127]
其他仿真参数为：雷达频率，间隔0.1ghz，采样时间1s，雷达视线角36
°
。总共生成5000张图片。
[0128]
空间非合作锥体目标的散射中心运动轨迹可以表示为：
[0129][0130]
其中分别代表空间非合作锥体目标的运动幅值、中心位置、进动频率、初始相位、径向速度。将锥体目标散射中心运动轨迹轨迹与卡尔曼滤波算法相结合，实现散射中心的轨迹提取，进一步将混叠的群目标分离，基于分离获得的单个锥体目标散射中心模型，重构获得单个锥体目标一维距离像和微多普勒时频图各15000张，利用分重构后的单个锥体目标一维距离像和微多普勒时频图，对双通道卷积神经网络进行训练和测试，提取单个锥体目标的进动频率，高度，半径，其中一维距离像和微多普勒时频图各13000张作为训练集，2000张作为测试集。
[0131]
本发明使用相对均方根误差作为判别微动参数预测准确性的指标。相对均方根误差公式为
[0132]
[0133]
其中k=2000，为测试集的数量，与第个测试值与真实值。最终估计结果如图7所示。
[0134]
从以上算例仿真和结果分析来看，本发明的基于散射中心模型群目标分离重构方法可以对不同尺寸和运动状态下的空间群目标一维距离像提取散射中心模型，并分离获得单个锥体目标散射中心模型，基于单个锥体目标散射中心模型重构一维距离像和微多普勒时频图，然后利用重构的单个锥体目标一维距离像微多普勒时频图，对双通道卷积神经网络进行训练，提取目标本征参数。在噪声环境下本发明所使用的方法仍具有较高的准确度和可靠性，说明本发明所使用的具有更高的抗干扰能力和泛化能力，相对于传统群目标本征参数提取方法在精度和稳定性上均有一定的提升。

技术特征：
1.一种基于物理先验知识的空间群目标探测方法，其特征在于，包括以下步骤：步骤1、构建不同尺寸的空间非合作锥体目标物理仿真模型，获得不同运动状态下混叠的空间群目标随时间变化的雷达回波，利用傅里叶变换获得空间群目标随时间变化的一维距离像；步骤2、基于锥体目标散射中心模型的物理先验知识，根据多个随机时刻的一维距离像估算出空间群目标散射中心的数量，结合最小二乘法，提取出所有时刻散射中心的位置信息和幅度信息；步骤3、基于卡尔曼滤波算法，将混叠的锥体目标散射中心模型中的位置信息重组，减少目标之间的遮挡和交叉干扰，实现空间群目标的分离，获得单个锥体目标散射中心模型，并重构单个锥体目标微多普勒时频图和一维距离像；步骤4、利用重构获得的单个锥体目标微多普勒时频图数据集和一维距离像数据集，结合双通道卷积神经网络进行训练和测试，得到复杂环境中空间非合作锥体目标的高度、半径、进动频率。2.根据权利要求1所述的基于物理先验知识的空间群目标探测方法，其特征在于，步骤1中构建不同尺寸的空间非合作锥体目标物理仿真模型，获得不同运动状态下混叠的空间群目标随时间变化的雷达回波，利用傅里叶变换获得空间群目标随时间变化的一维距离像，具体如下：对于旋转对称锥体目标，散射中心分布在两个位置，分别为锥体目标顶端和锥体目标底部边缘；锥体目标的运动参数和结构参数为：进动频率，速度，进动角，高度，半径，质心高；定义进动角频率为，为当前时间，为雷达俯仰角，为雷达方位角；在以目标质心为原点的坐标系下，两个位置处的散射中心径向位移理论表达式为：锥体目标顶端散射中心a的径向位移表达式为：(1)锥体目标底部边缘散射中心b的径向位移表达式为：(2)理想散射中心模型表达式为：(3)其中为第i个散射中心的幅度信息，为第i个散射中心的位置信息，为雷达发射
波长，为雷达回波；由式(3)进一步推导出混叠的锥体目标雷达回波表达式为：(4)其中、为第i个锥体目标锥顶和锥底的幅度信息，、为第i个锥体目标锥顶和锥底的位置信息。3.根据权利要求2所述的基于物理先验知识的空间群目标探测方法，其特征在于，步骤2中基于锥体目标散射中心模型的物理先验知识，根据多个随机时刻的一维距离像估算出空间群目标散射中心的数量，结合最小二乘法，提取出所有时刻散射中心的位置信息和幅度信息，具体如下：由锥体目标的顶端散射中心a和底部边缘散射中心b的径向位移表达式(1)、式(2)获得锥体目标散射中心模型的物理先验知识，即噪声环境下的锥体目标随时间变化的散射中心模型为：(5)其中，为t时刻的雷达回波；表示t时刻噪声，、为第i个锥体目标锥顶、锥底t时刻的幅度信息，、为第i个锥体目标锥顶、锥底t时刻的位置信息，为第i个锥体目标的速度，为第i个锥体目标的进动角，为第i个锥体目标的高度，为第i个锥体目标的半径，为第i个锥体目标的质心高，为第i个锥体目标的角速度；对噪声环境下的锥体目标随时间变化的散射中心模型进行变化，得到t时刻一维距离像的噪声：(6)其中为t时刻一维距离像的噪声；
通过多个随机时刻一维距离像的时间累积，获得增强散射中心的一维距离像；基于最小二乘法和锥体目标散射中心模型物理先验知识，提取一维距离像的散射中心数量，以及散射中心的位置信息和幅度信息，代入式(6)获得t时刻一维距离像的噪声，将多个随机时刻一维距离像的噪声最大值的均值作为最终的噪声阈值；根据获得的空间群目标散射中心的数量和噪声阈值，基于最小二乘法和锥体目标散射中心模型物理先验知识，对所有时刻的空间群目标一维距离像提取散射中心的位置信息和幅度信息，获得空间群目标散射中心模型。4.根据权利要求3所述的基于物理先验知识的空间群目标探测方法，其特征在于，空间群目标一维距离像的散射中心提取具体流程如下：(2.1)随机选取多段连续时刻的一维距离像，将每段连续时刻的一维距离像做时间维度上的累加，获得增强散射中心的一维距离像；(2.2)将增强散射中心的一维距离像的均值作为起始噪声阈值，幅值大于噪声阈值的视为散射中心，幅值低于噪声阈值的视为噪声，基于最小二乘法和锥体目标散射中心模型物理先验知识，获得步骤(2.1)中多段一维距离像的散射中心数量，以及散射中心的位置信息和幅度信息；(2.3)将步骤(2.2)中获得的散射中心数量中的众数作为空间群目标散射中心的数量上限m0，并且锥体目标数量为m0/2，将步骤(2.2)中获得的散射中心的位置信息和幅度信息代入式(6)，获得t时刻一维距离像的噪声，将各个时刻一维距离像的噪声最大值的均值作为最终的噪声阈值n0；(2.4)选取t时刻的一维距离像l，设置初始散射中心个数n=0；(2.5)基于最小二乘法，提取t时刻的一维距离像l中幅度最大的散射中心的位置信息和幅度信息，将获得的散射中心的位置信息和幅度信息代入式(4)，重新生成无噪声一维距离像，与原本的一维距离像相减获得提取后的一维距离像l’，且散射中心个数n=n+1；(2.6)判断t时刻的散射中心个数是否大于等于m0或者一维距离像l’的最大值小于n0：如果判断为是，则表示t时刻一维距离像的散射中心提取完备，令t=t+1，并返回步骤(2.4)；如果判断为否，则表示t时刻一维距离像的散射中心未提取完备，令t时刻一维距离像l=l’，并返回步骤(2.5)。5.根据权利要求1所述的基于物理先验知识的空间群目标探测方法，其特征在于，步骤3中基于卡尔曼滤波算法，将混叠的锥体目标散射中心模型中的位置信息重组，减少目标之间的遮挡和交叉干扰，实现空间群目标的分离，获得单个锥体目标散射中心模型，并重构单个锥体目标微多普勒时频图和一维距离像，具体如下：通过步骤2获得的空间群目标散射中心模型，基于已获得的散射中心位置信息，结合卡尔曼滤波算法和锥体目标散射中心运动轨迹为正弦函数的性质，计算锥体目标散射中心运动轨迹相关参数，获得锥体目标散射中心的运动轨迹，将存在交叉遮挡的多个散射中心运动轨迹分离出来，锥体目标散射中心运动轨迹表达式如下所示：
(7)其中，、分别为第i个锥体目标的锥顶和锥底运动轨迹，为第i个锥体目标的运动幅值，为第i个锥体目标的运动角速度，为第i个锥体目标的初始相位，、分别为第i个锥体目标的锥顶和锥底初始位置，为第i个锥体目标的径向速度。6.根据权利要求5所述的基于物理先验知识的空间群目标探测方法，其特征在于，分离获得单个锥体目标散射中心模型，并重构单个锥体目标微多普勒时频图和一维距离像，具体流程如下：(3.1)排除步骤2中提取出的散射中心数量小于m0的时刻，即散射中心被遮挡时刻，并计算时刻，所有散射中心位置信息和幅度信息相关系数，选取相关系数最大，即不存在散射中心运动轨迹交叉的时刻作为为起始时刻t，；(3.2)根据t时刻的第i个锥体目标散射中心的位置和锥体目标散射中心运动轨迹公式(7)，结合卡尔曼滤波算法，预测t+1时刻的第i个锥体目标散射中心的位置；(3.3)计算t+1时刻的预测值与t+1时刻的第j个散射中心的位置信息的误差，判断是否存在大于阈值，该阈值为锥体散射中心运动幅值的5%：如果判断为是，则结合最小二乘法，根据已知的0~t时刻的第i个锥体目标散射中心的运动轨迹，以及锥体目标散射中心运动轨迹表达式(7)，获得第i个锥体目标的运动参数，并重复步骤(3.2)~步骤(3.3)；如果判断为否，则将预测值作为第i个锥体目标散射中心t+1时刻散射中心的位置，并且t=t+1，重复步骤(3.2)~步骤(3.3)；当遍历完所有时刻后，进入步骤(3.4)；(3.4)当遍历完所有时刻后，根据最终获得的第i个散射中心的运动轨迹，结合最小二乘法，获得锥体目标散射中心运动轨迹的相关参数、、、、、，通过公式(7)补全步骤(3.1)中排除的被遮挡时刻散射中心的位置，并由插值获得对应时刻的散射中心的幅度，最终获得完整的第i个锥体目标的锥顶散射中心的位置信息与幅度信息、锥底散射中心的位置信息与幅度信息，即获得单个锥体目标的散射中心模型；(3.5)将获得的单个锥体目标的散射中心模型，代入公式(4)获得单个锥体目标的雷达
回波，通过时频分析和傅里叶变换获得单个锥体目标的微多普勒时频图和一维距离像。7.一种基于物理先验知识的空间群目标探测系统，其特征在于，该系统用于实现权利要求1~6任一项所述的基于物理先验知识的空间群目标探测方法，具体包括一维距离像获取模块、散射中心信息提取模块、单个锥体目标重构模块、空间非合作锥体目标特征提取模块，其中：一维距离像获取模块，用于构建不同尺寸的空间非合作锥体目标物理仿真模型，获得不同运动状态下混叠的空间群目标随时间变化的雷达回波，利用傅里叶变换获得空间群目标随时间变化的一维距离像；散射中心信息提取模块，基于锥体目标散射中心模型的物理先验知识，根据多个随机时刻的一维距离像估算出空间群目标散射中心的数量，结合最小二乘法，提取出所有时刻散射中心的位置信息和幅度信息；单个锥体目标重构模块，基于卡尔曼滤波算法，将混叠的锥体目标散射中心模型中的位置信息重组，减少目标之间的遮挡和交叉干扰，实现空间群目标的分离，获得单个锥体目标散射中心模型，并重构单个锥体目标微多普勒时频图和一维距离像；空间非合作锥体目标特征提取模块，利用重构获得的单个锥体目标微多普勒时频图数据集和一维距离像数据集，结合双通道卷积神经网络进行训练和测试，得到复杂环境中空间非合作锥体目标的高度、半径、进动频率。8.一种移动终端，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现如权利要求1~6任一项所述的基于物理先验知识的空间群目标探测方法。

技术总结
本发明公开了一种基于物理先验知识的空间群目标探测方法及系统，具体为：构建不同尺寸的空间非合作锥体目标物理仿真模型，获得不同运动状态下混叠的空间群目标随时间变化的雷达回波，利用傅里叶变换获得空间群目标随时间变化的一维距离像；提取出所有时刻散射中心的位置信息和幅度信息，将混叠的锥体目标散射中心模型中的位置信息重组，重构单个锥体目标微多普勒时频图和一维距离像；结合双通道卷积神经网络进行训练和测试，得到复杂环境中空间非合作锥体目标的高度、半径、进动频率。本发明在复杂电磁环境下仍能对空间群目标进行有效分离重构，减少目标之间的交叉和遮挡等干扰，具有更高的提取精度和抗干扰能力。具有更高的提取精度和抗干扰能力。具有更高的提取精度和抗干扰能力。


技术研发人员：赖瑞鑫 李猛猛 丁大志 王绍然 史琪 张晋伟 曾董 曹仁猛 文明浩 胡悦 卢其然
受保护的技术使用者：南京理工大学
技术研发日：2023.06.14
技术公布日：2023/7/17