进动目标结构参数提取方法

摘要

本发明公开了一种进动目标结构参数提取方法，其特征在于，包括：分别通过两部雷达的目标一维距离像获取两部雷达的目标径向长度信息；根据获取的目标径向长度信息计算两部雷达的目标径向长度曲线；利用两部雷达的目标径向长度曲线计算目标进动角、及两部雷达的进动轴视线角差值；基于目标进动角、两部雷达的进动轴视线角差值、及两部雷达的目标径向长度曲线计算目标长度及目标底面半径。与现有技术相比，本发明能够以较小的误差精确提取合作/非合作目标的结构参数及进动角，同时大幅减小提取进动特征所需的观测时间。

权利要求

1.一种进动目标结构参数提取方法，其特征在于，包括：

S1.通过第一雷达的目标一维距离像获取第一雷达目标径向长度信息；通过第二雷达的目标一维距离像获取第二雷达目标径向长度信息；

S2.根据第一雷达目标径向长度信息计算第一雷达目标径向长度曲线；根据第二雷达目标径向长度信息计算第二雷达目标径向长度曲线；

S3.利用第一雷达目标径向长度曲线、及第二雷达目标径向长度曲线计算目标进动角、及第一雷达与第二雷达的进动轴视线角差值；

S4.基于目标进动角、进动轴视线角差值、第一雷达目标径向长度曲线及第二雷达目标径向长度曲线计算目标长度及目标底面半径；

所述进动轴视线角为雷达视线与目标进动轴的夹角。

2.如权利要求1所述的方法，其特征在于，步骤S3具体包括：S31.利用第一雷达目标径向长度曲线、第二雷达目标径向长度曲线及公式1计算目标进动角；

S32.利用第一雷达目标径向长度曲线、第二雷达目标径向长度曲线、目标进动角及公式2计算第一雷达与第二雷达的进动轴视线角差值；

其中，θ为目标进动角；X1为第一雷达在第一时间段T1内，目标径向长度的最大值与最小值的和；Y1为第一雷达在第一时间段T1内，目标径向长度的最大值与最小值的差；X2为第二雷达在第二时间段T2内，目标径向长度的最大值与最小值的和；Y2为第二雷达在第二时间段T2内，目标径向长度的最大值与最小值的差；αi为第一雷达的进动轴视线角，且αi从0°到90°按照第一步进值设置；βi为与αi对应的第二雷达的进动轴视线角；i＞1、且i∈N，N为大于1的自然数。

3.如权利要求2所述的方法，其特征在于，步骤S4具体包括：S41.基于目标进动角、进动轴视线角差值及公式3、公式4计算对应于每个第一雷达进动轴视线角αi的目标样本长度Li及目标样本底面半径Ri；

其中，

S42.根据目标样本长度Li及目标样本底面半径Ri提取目标长度及目标底面半径；

其中，T11…T1N为T1内的N个时间点，T21…T2M为T2内的M个时间点；M＞1，N＞1，M∈N，N∈N；

Len(T11)…Len(T1N)为与T11…T1N对应的目标径向长度，Len(T21)…Len(T2M)为与T21…T2M对应的目标径向长度；ω为目标进动角速度； 分别为第一雷达、第二雷达视角下的目标进动初相角；τ11…τ1N为与T11…T1N对应的第一雷达视线角，τ21…τ2M为与T21…T2M对应的第二雷达视线角；且

4.如权利要求3所述的方法，其特征在于，步骤S42具体包括：S421.根据公式5计算每一组目标样本长度Li及目标样本底面半径Ri的误差值e(αi)；

S422.比较所有误差值的大小，并提取最小的误差值对应的目标样本长度作为目标长度，提取最小的误差值对应的目标样本底面半径作为目标底面半径。

5.如权利要求1所述的方法，其特征在于，步骤S1具体为：

通过第一雷达的目标一维距离像，利用最优路径法获取第一雷达目标径向长度信息；

通过第二雷达的目标一维距离像，利用最优路径法获取第二雷达目标径向长度信息。

说明书

进动目标结构参数提取方法

技术领域

[0001] 本发明涉及目标特征提取领域，尤其涉及一种进动目标结构参数提取方法。

背景技术

[0002] 弹道导弹中段相对于助推段和再入段具有飞行时间较长、空间环境相对简单等特点，是攻防对抗的重要阶段。弹头在中段的自旋、进动等运动形式构成了弹头目标的微动特性，可以作为弹头类目标识别的重要特征。对于进动目标，其周期特性十分明显，其RCS(Radar Cross Section，雷达散射截面)起伏周期与进动周期一致，而其进动周期可通过窄带RCS获得，但是目标的进动角及结构信息无法通过目标RCS直接获得。

[0003] 目前进动目标的特征提取方法主要是通过目标RCS或目标一维距离像获得，且二者主要是通过单雷达获取的信息进行提取。前者通过对目标RCS进行多项式拟合，利用圆锥的RCS估算公式作为目标的模板信息来提取进动参数，但是需要了解目标在各种姿态下的RCS特性或弹头类目标的形状，并不适用于非合作目标，局限性较大；第二种方法是通过目标一维距离像求取径向长度或重构目标散射中心来提取进动特征，但是需要较长的观测时间且误差较大。

[0004] 因此，亟需一种适用于非合作目标、不依赖于较长观测时间且误差较小的进动目标结构参数提取方法以解决上述问题。

发明内容

[0005] 本发明提供了一种进动目标结构参数提取方法，利用多部雷达在短时间内观测到的回波数据获取目标径向长度信息，通过曲线拟合得到目标径向长度曲线，在此基础上提取目标长度及底面半径。与现有技术相比，本发明能够以较小的误差精确提取合作/非合作目标的结构参数及进动角，同时大幅减小提取进动特征所需的观测时间。

[0006] 本发明提供一种进动目标结构参数提取方法，包括：S1.通过第一雷达的目标一维距离像获取第一雷达目标径向长度信息；通过第二雷达的目标一维距离像获取第二雷达目标径向长度信息；S2.根据第一雷达目标径向长度信息计算第一雷达目标径向长度曲线；根据第二雷达目标径向长度信息计算第二雷达目标径向长度曲线；S3.利用第一雷达目标径向长度曲线、及第二雷达目标径向长度曲线计算目标进动角、及第一雷达与第二雷达的进动轴视线角差值；S4.基于目标进动角、进动轴视线角差值、第一雷达目标径向长度曲线及第二雷达目标径向长度曲线计算目标长度及目标底面半径；所述进动轴视线角为雷达视线与目标进动轴的夹角。

[0007] 优选地，步骤S3具体包括：S31.利用第一雷达目标径向长度曲线、第二雷达目标径向长度曲线及公式1计算目标进动角；

[0008]     公式1

[0009] S32.利用第一雷达目标径向长度曲线、第二雷达目标径向长度曲线、目标进动角及公式2计算第一雷达与第二雷达的进动轴视线角差值；

[0010]     公式2

[0011] 其中，θ为目标进动角；X1为第一雷达在第一时间段T1内，目标径向长度的最大值与最小值的和；Y1为第一雷达在第一时间段T1内，目标径向长度的最大值与最小值的差；X2为第二雷达在第二时间段T2内，目标径向长度的最大值与最小值的和；Y2为第二雷达在第二时间段T2内，目标径向长度的最大值与最小值的差；αi为第一雷达的进动轴视线角，且αi从0°到90°按照第一步进值设置；βi为与αi对应的第二雷达的进动轴视线角；i＞1、且i∈N。

[0012] 优选地，步骤S4具体包括：S41.基于目标进动角、进动轴视线角差值及公式3、公式4计算对应于每个第一雷达进动轴视线角αi的目标样本长度Li及目标样本底面半径Ri；

[0013]    公式3

[0014] 其中，

[0015]    公式4

[0016] S42.根据目标样本长度Li及目标样本底面半径Ri提取目标长度及目标底面半径；

[0017] 其中，T11…T1N为T1内的N个时间点，T21…T2M为T2内的M个时间点；M＞1，N＞1，M∈N，N∈N；Len(T11)…Len(T1N)为与T11…T1N对应的目标径向长度，Len(T21)…Len(T2M)为与T21…T2M对应的目标径向长度；ω为目标进动角速度； 分别为第一雷达、第二雷达视角下的目标进动初相角；τ11…τ1N为与T11…T1N对应的第一雷达视线角，τ21…τ2M为与T21…T2M对应的第二雷达视线角；且

[0018]

[0019] 优选地，步骤S42具体包括：S421.根据公式5计算每一组目标样本长度Ri及目标样本底面半径Li的误差值e(αi)；

[0020] 公式5

[0021] S422.比较所有误差值的大小，并提取最小的误差值对应的目标样本长度作为目标长度，提取最小的误差值对应的目标样本底面半径作为目标底面半径。

[0022] 优选地，步骤S2具体包括：S21.根据第一雷达目标径向长度信息，确定第一雷达目标径向长度曲线的参数l11、l12、 的取值范围；根据第二雷达目标径向长度信息，确定第二雷达目标径向长度曲线的参数l21、l22、 的取值范围；S22.基于l11、l12、 的取值范围及精度要求，设置分别与l11、l12、 对应的取值间隔，对l11、l12、 进行离散取值，生成l11、l12、参数数组；并基于l21、l22、 的取值范围及精度要求，设置分别与l21、l22、 对应的取值间隔，对l21、l22、 进行离散取值，生成l21、l22、 参数数组；S23.利用第一雷达获取的目标径向长度信息、及公式6检验l11、l12、 参数数组；根据检测结果确定第一雷达目标径向长度曲线的参数；并利用第二雷达获取的目标径向长度信息、及公式7检验l21、l22、 参数数组；根据检测结果确定第二雷达目标径向长度曲线的参数；

[0023]    公式6

[0024]    公式7

[0025] 所述第一雷达获取的目标径向长度信息包括时间点t11…t1P及与该时间点对应的目标径向长度Len(t11)…Len(t1P)；所述第二雷达获取的目标径向长度信息包括时间点t21…t2Q及与该时间点对应的目标径向长度Len(t21)…Len(t2Q)；P＞1，Q＞1，P∈N，Q∈N。

[0026] 优选地，步骤S23具体包括：S231.选取l11、l12、 参数数组中任一数组，将t11…t1P及该任一数组代入公式6，得到P个运算结果；分别计算P个运算结果与Len(t11)…Len(t1P)的差值，得到l11、l12、 参数数组中该任一数组的P个差值；并选取l21、l22、 参数数组中任一数组，将t21…t2Q及该任一数组代入公式7，得到Q个运算结果；分别计算Q个运算结果与Len(t21)…Len(t2Q)的差值，得到l21、l22、 参数数组中该任一数组的Q个差值；

[0027] S232.比较l11、l12、 参数数组中任一数组的P个差值与第一阈值的大小，统计P个差值中小于第一阈值的差值数量，并将该差值数量作为l11、l12、 参数数组中该任一数组的合格数；并比较l21、l22、 参数数组中任一数组的Q个差值与第二阈值的大小，统计Q个差值中小于第二阈值的差值数量，并将该差值数量作为l21、l22、 参数数组中该任一数组的合格数；

[0028] S233.选取l11、l12、 参数数组中合格数最大的数组，将该数组的参数作为第一雷达目标径向长度曲线的参数；并选取l21、l22、 参数数组中合格数最大的数组，将该数组的参数作为第二雷达目标径向长度曲线的参数；所述第一阈值根据第一雷达精度要求设置；所述第二阈值根据第二雷达精度要求设置。

[0029] 优选地，步骤S1具体为：通过第一雷达的目标一维距离像，利用最优路径法获取第一雷达目标径向长度信息；通过第二雷达的目标一维距离像，利用最优路径法获取第二雷达目标径向长度信息。

[0030] 优选地，在步骤S1之前，所述方法还包括:S0.从U部雷达中任取两部合为一组，形成V组雷达；根据所述V组雷达中任一组中的两部雷达提取目标长度Lj及目标底面半径Rj；在步骤S4之后，所述方法还包括：S5.根据目标长度Lj及目标底面半径Rj计算最佳目标长度及最佳目标底面半径；其中，U＞2，且U∈N； 1≤j≤V，且j∈N。

[0031] 优选地，步骤S5具体为：根据公式8计算Lj的平均值L#，将Lj的平均值L#作为最佳目标长度；根据公式9计算Rj的平均值R#，将Rj的平均值R#作为最佳目标底面半径；

[0032]    公式8

[0033]    公式9。

[0034] 优选地，步骤S5具体为：根据公式10计算Lj的加权平均值L*，将Lj的加权平均值L*作为最佳目标长度；根据公式11计算Rj的加权平均值R*，将Rj的加权平均值R*作为最佳目标底面半径；

[0035]    公式10

[0036]    公式11

[0037] 其中，Wj为V组雷达中第j组雷达的权值，Wj由第j组中两部雷达中的每一雷达的性能参数确定。

[0038] 根据本发明提供的进动目标结构参数提取方法，能够以较小的误差精确提取合作/非合作目标的结构参数及进动角，同时大幅减小提取进动特征所需的观测时间。

附图说明

[0039] 图1是本发明的进动目标结构参数提取方法的进动目标模型图；

[0040] 图2是本发明的进动目标结构参数提取方法的第一流程图；

[0041] 图3是本发明的进动目标结构参数提取方法的第二流程图。

具体实施方式

[0042] 为使本发明的目的、技术方案及优点更加清楚明白，以下参照附图并举出优选实施例，对本发明进一步详细说明。然而，需要说明的是，说明书中列出的许多细节仅仅是为了使读者对本发明的一个或多个方面有一个透彻的理解，即便没有这些特定的细节也可以实现本发明的这些方面。

[0043] 现有技术中的进动目标结构参数提取方法多用于合作目标，其观测时间较长且误差较大。

[0044] 针对上述技术问题，本发明提供了一种基于多雷达的进动目标结构参数提取方法，能够以较小的误差精确提取合作/非合作目标的结构参数及进动角，同时大幅减小提取进动特征所需的观测时间。

[0045] 图1是本发明的进动目标模型图。在本发明的优选实施例中，进动目标为圆锥目标。O2-xRyRzR构成了雷达坐标系，O1-xjyjzj构成了进动轴坐标系；在进动轴坐标系中，以目标进动中心为原点，目标进动轴为Xj轴，Zj轴在电磁波入射方向和Xj轴构成的平面内，xj轴与yj轴、zj轴构成右手螺旋；由于存在攻角，故目标进动轴和目标速度方向并不一致。

[0046] 图1中的目标长度为L，底面半径为R，则目标半锥角ε＝tan-1(R/L)，目标本体轴绕进动轴以角速率ω、进动角θ逆时针旋转，锥顶进动的初相角为 在雷达坐标系中，雷达Ra的视线与目标进动轴的夹角(进动轴视线角)为α，与本体轴的夹角(即雷达视线角)为τ。圆锥目标锥顶到进动中心距离为L1，进动中心到锥底的距离为L2，显然有L＝L1+L2。根据散射中心理论，当电磁波前向入射时，圆锥目标有三个强散射中心，分别为锥顶散射中心P1，两个锥底边缘散射中心P2和P3。

[0047] 在t时刻，进动中心与锥顶P1构成的单位矢量为：

[0048]

[0049] 雷达视线与进动中心构成的单位矢量为：

[0050] RaO1＝[cosα,0,-sinα,]

[0051] 则O1P1和RaO1的夹角即为雷达视线角，且易知：

[0052]

[0053] 目标径向长度为在一维距离像中相距最远的两个散射中心在径向上的投影距离。令进动中心为相位零点，f()表示投影运算，则圆锥目标三个散射中心在雷达径向上的投影位置分别为：

[0054] f([P1,P2,P3])＝[L1cosτ,-L2cosτ+Rsinτ,-L2cosτ-Rsinτ]

[0055] 根据电磁散射理论，当雷达视线角大于目标半锥角时，散射中心P3受到遮挡而不可见。而在一般情况下，由于受攻角影响，雷达视线角大于圆锥目标半锥角，故散射中心P3是难以被观测到的。则圆锥目标径向长度为散射中心P1到P2径向投影距离：

[0056] f(P1)-f(P2)＝L1cosτ+L2cosτ-Rsinτ

[0057] ＝Lcosτ-Rsinτ

[0058] 由于雷达视线角随时间变化，故径向长度亦是一个随时间变化的序列，令径向长度为Len(t)，则其满足以下方程组：

[0059]    公式12

[0060] 在实际应用中，R远小于L，则：

[0061]    公式13

[0062] 令l1＝Lcosθcosα，l2＝Lsinθsinα，则:

[0063]    公式14

[0064] 其中ω可以通过窄带分析获得，上述方程的未知参数为l1、l2、

[0065] 至此得到了圆锥目标的径向长度理论公式及实际应用公式，在此基础上，本发明提供了一种进动目标结构参数提取方法，如图2所示，包括：

[0066] S1.通过第一雷达的目标一维距离像获取第一雷达目标径向长度信息；通过第二雷达的目标一维距离像获取第二雷达目标径向长度信息。

[0067] 在本发明优选实施例中，步骤S1具体为：

[0068] 通过第一雷达的目标一维距离像，利用最优路径法获取第一雷达目标径向长度信息；通过第二雷达的目标一维距离像，利用最优路径法获取第二雷达目标径向长度信息。最优路径法能够从目标一维距离像中有效提取目标径向长度信息，在强噪声环境中也能保证较高的准确率。

[0069] S2.根据第一雷达目标径向长度信息计算第一雷达目标径向长度曲线；

[0070] 根据第二雷达目标径向长度信息计算第二雷达目标径向长度曲线。

[0071] 在本发明优选实施例中，步骤S2具体包括：

[0072] S21.根据第一雷达目标径向长度信息，确定第一雷达目标径向长度曲线的参数l11、l12、 的取值范围；根据第二雷达目标径向长度信息，确定第二雷达目标径向长度曲线的参数l21、l22、 的取值范围。

[0073] S22.基于l11、l12、 的取值范围及精度要求，设置分别与l11、l12、 对应的取值间隔，对l11、l12、 进行离散取值，生成l11、l12、 参数数组；上述l11、l12、 参数数组包含经过离散取值后的所有l11、l12、 的组合，其中的任一数组皆为第一雷达目标径向长度曲线的待定参数；

[0074] 并且，基于l21、l22、 的取值范围及精度要求，设置分别与l21、l22、 对应的取值间隔，对l21、l22、 进行离散取值，生成l21、l22、 参数数组；上述l21、l22、 参数数组包含经过离散取值后的所有l21、l22、 的组合，其中的任一数组皆为第二雷达目标径向长度曲线的待定参数。

[0075] S23.利用第一雷达获取的目标径向长度信息、及公式6检验l11、l12、 参数数组；根据检测结果确定第一雷达目标径向长度曲线的参数；并

[0076] 利用第二雷达获取的目标径向长度信息、及公式7检验l21、l22、 参数数组；根据检测结果确定第二雷达目标径向长度曲线的参数；

[0077]    公式6

[0078]    公式7

[0079] 公式6、公式7由前述公式14得到；

[0080] 第一雷达获取的目标径向长度信息包括时间点t11…t1P及与该时间点对应的目标径向长度Len(t11)…Len(t1P)；

[0081] 第二雷达获取的目标径向长度信息包括时间点t21…t2Q及与该时间点对应的目标径向长度Len(t21)…Len(t2Q)；P＞1，Q＞1，P∈N，Q∈N。

[0082] 在本发明优选实施例中，步骤S23具体包括：

[0083] S231.选取l11、l12、φ1参数数组中任一数组，将t11…t1P及该任一数组代入公式6，得到P个运算结果；分别计算P个运算结果与Len(t11)…Len(t1P)的差值，得到l11、l12、 参数数组中该任一数组的P个差值；

[0084] 并且，选取l21、l22、 参数数组中任一数组，将t21…t2Q及该任一数组代入公式7，得到Q个运算结果；分别计算Q个运算结果与Len(t21)…Len(t2Q)的差值，得到l21、l22、 参数数组中该任一数组的Q个差值。

[0085] S232.比较l11、l12、 参数数组中任一数组的P个差值与第一阈值的大小，统计P个差值中小于第一阈值的差值数量，并将该差值数量作为l11、l12、 参数数组中该任一数组的合格数；第一阈值根据第一雷达精度要求设置；

[0086] 并且，比较l21、l22、 参数数组中任一数组的Q个差值与第二阈值的大小，统计Q个差值中小于第二阈值的差值数量，并将该差值数量作为l21、l22、 参数数组中该任一数组的合格数；第二阈值根据第二雷达精度要求设置。

[0087] S233.选取l11、l12、 参数数组中合格数最大的数组，将该数组的参数作为第一雷达目标径向长度曲线的参数；

[0088] 并且，选取l21、l22、 参数数组中合格数最大的数组，将该数组的参数作为第二雷达目标径向长度曲线的参数。

[0089] 步骤S2利用曲线拟合的方法，通过广义Hough变换(GHT)准确地确定了第一雷达及第二雷达的目标径向长度曲线。

[0090] 在本发明优选实施例中，通过广义Hough变换的另一种形式得到目标径向长度曲线，步骤S23具体为：

[0091] S234.将l11、l12、 参数数组中的每一数组逐行排列，生成a*3的l11、l12、 参数矩阵，同时构建与该参数矩阵同型的零矩阵作为第一累加矩阵；l11、l12、 参数矩阵中的每一参数数组与第一累加矩阵中每一行对应；a为l11、l12、 参数矩阵中参数数组的数量；

[0092] 并且，将l21、l22、 参数数组中的每一数组逐行排列，生成b*3的l21、l22、 参数矩阵，同时构建与该参数矩阵同型的零矩阵作为第二累加矩阵；l21、l22、 参数矩阵中的每一参数数组与第二累加矩阵中每一行对应；b为l21、l22、 参数矩阵中参数数组的数量。

[0093] S235.将第一雷达获取的目标径向长度信息中的任一时间点、及l11、l12、 参数矩阵中的每一参数数组代入公式6，得到对应于每一参数数组的a个运算结果；计算上述a个运算结果与该时间点对应的目标径向长度的差值；将得到的对应于每一参数数组的a个差值与第一阈值进行比较，标记差值小于第一阈值的参数数组，并在该参数数组对应的第一累加矩阵的行加1；

[0094] 并且，将第二雷达获取的目标径向长度信息中的任一时间点、及l21、l22、 参数矩阵中的每一参数数组代入公式7，得到对应于每一参数数组的b个运算结果；计算上述b个运算结果与该时间点对应的目标径向长度的差值；将得到的对应于每一参数数组的b个差值与第二阈值进行比较，标记差值小于第二阈值的参数数组，并在该参数数组对应的第二累加矩阵的行加1。

[0095] S236.将第一雷达获取的目标径向长度信息中的P个时间点按照步骤S235依次执行；选取第一累加矩阵中最大值所在行对应的参数数组作为第一雷达目标径向长度曲线的参数；

[0096] 并且，将第二雷达获取的目标径向长度信息中的Q个时间点按照步骤S235依次执行；选取第二累加矩阵中最大值所在行对应的参数数组作为第二雷达目标径向长度曲线的参数。

[0097] S3.利用第一雷达目标径向长度曲线、及第二雷达目标径向长度曲线计算目标进动角、及第一雷达与第二雷达的进动轴视线角差值，所述进动轴视线角为雷达视线与目标进动轴的夹角。

[0098] 在本发明优选实施例中，步骤S3具体包括：

[0099] S31.利用第一雷达目标径向长度曲线、第二雷达目标径向长度曲线及公式1计算目标进动角。

[0100]    公式1

[0101] 其中，X1为第一雷达在第一时间段T1内，目标径向长度的最大值与最小值的和；Y1为第一雷达在第一时间段T1内，目标径向长度的最大值与最小值的差；X2为第二雷达在第二时间段T2内，目标径向长度的最大值与最小值的和；Y2为第二雷达在第二时间段T2内，目标径向长度的最大值与最小值的差。

[0102] S32.利用第一雷达目标径向长度曲线、第二雷达目标径向长度曲线、目标进动角及公式2计算第一雷达与第二雷达的进动轴视线角差值。

[0103]    公式2

[0104] 其中，αi为第一雷达的进动轴视线角，且αi从0°到90°按照第一步进值设置；βi为与αi对应的第二雷达的进动轴视线角；i＞1、且i∈N。

[0105] 在本发明的优选实施例中，若第一步进值为5°，则αi为0°、5°、10°、15°、20°、25°……90°；若第一步进值为10°，则αi为10°、20°、30°……90°。

[0106] 公式1、公式2的推导过程如下：

[0107] 假设在弹道中段某短时间段T内，雷达进动轴视线角近似不变，通过步骤S2，在时间段T1内第一雷达能够得到一组径向长度Len1，在时间段T2内第二雷达能够得到一组径向长度Len2，即便T1和T2都小于一个进动周期，都能够估计径向长度曲线，从而获得一个周期T内径向长度的最大值及最小值。对于第一雷达，根据公式12，可知径向长度的最大值及最小值出现在： 和 对应的径向长度分别为：

[0108] Len(t1)＝Lcos(α-θ)-Rsin(α-θ)

[0109] Len(t2)＝Lcos(α+θ)-Rsin(α+θ)

[0110] 则在T1时间段内的径向长度最大值及最小值为：

[0111] max(Len1)＝Lcos(α-θ)-Rsin(α-θ)

[0112] min(Len1)＝Lcos(α+θ)-Rsin(α+θ)

[0113] 令：

[0114] X1＝max(Len1)+min(Len1)

[0115] Y1＝max(Len1)-min(Len1)

[0116] 可得：

[0117] X1＝2Lcosαcosθ-2Rcosθsinα

[0118] Y1＝2Lsinαsinθ+2Rcosαsinθ   公式15

[0119] 由此：

[0120]

[0121]

[0122] 由于cos2α+sin2α＝1，可得：

[0123]

[0124] 同理，对于第二雷达：

[0125]

[0126] 将上述两式联立可得：

[0127]

[0128] 稍作变形可得：

[0129]     公式1

[0130] 同时，根据公式15可得：

[0131]

[0132]

[0133] 同理，对于第二雷达可得到：

[0134]

[0135]

[0136] 将上述两式联立可得：

[0137]

[0138] 在步骤S3中，第一雷达及第二雷达在较短的观测时间内，即可准确计算目标进动角及进动轴视线角差值，在此基础上，本发明能够实现短时间对目标结构参数的提取；与现有技术相比，大幅减小了提取进动特征所需的观测时间。同时，通过公式1，本发明实现了以较高精度对目标进动角的提取。

[0139] S4.基于目标进动角、进动轴视线角差值、第一雷达目标径向长度曲线及第二雷达目标径向长度曲线计算目标长度及目标底面半径。

[0140] 在本发明优选实施例中，步骤S4具体包括：

[0141] S41.基于目标进动角、进动轴视线角差值及公式3、公式4计算对应于每一个第一雷达进动轴视线角αi的目标样本长度Li及目标样本底面半径Ri；

[0142]    公式3

[0143] 其中，

[0144]    公式4

[0145] 其中，T11…T1N为T1内的N个时间点，T21…T2M为T2内的M个时间点；M＞1，N＞1，M∈N，N∈N；Len(T11)…Len(T1N)为与T11…T1N对应的目标径向长度，Len(T21)…Len(T2M)为与T21…T2M对应的目标径向长度；ω为目标进动角速度； 分别为第一雷达、第二雷达视角下的目标进动初相角；τ11…τ1N为与T11…T1N对应的第一雷达视线角，τ21…τ2M为与T21…T2M对应的第二雷达视线角；且

[0146]

[0147] 公式3、公式4由公式12得到。

[0148] 在本发明优选实施例中，利用最小二乘法求得对应于每一个第一雷达进动轴视线角αi的目标样本长度Li及目标样本底面半径Ri。

[0149] S42.根据目标样本长度Li及目标样本底面半径Ri提取目标长度及目标底面半径。

[0150] 在本发明优选实施例中，步骤S42具体包括：

[0151] S421.根据公式5计算每一组目标样本长度Ri及目标样本底面半径Li的误差值e(αi)；

[0152]    公式5

[0153] S422.比较所有误差值的大小，并提取最小的误差值对应的目标样本长度作为目标长度，提取最小的误差值对应的目标样本底面半径作为目标底面半径。

[0154] 在步骤S4中，不需要了解进动目标在各种姿态下的RCS特性或形状，即可以较高的精度提取其结构参数。

[0155] 在本发明优选实施例中，本发明提供的进动目标结构参数提取方法，在步骤S1之前，还包括:

[0156] S0.从U部雷达中任取两部合为一组，形成V组雷达；根据所述V组雷达中任一组中的两部雷达提取目标长度Lj及目标底面半径Rj；

[0157] 在步骤S4之后，还包括：

[0158] S5.根据目标长度Lj及目标底面半径Rj计算最佳目标长度及最佳目标底面半径；其中，U＞2，且U∈N； 1≤j≤V，且j∈N。

[0159] 在本发明优选实施例中，步骤S5具体为：

[0160] 根据公式8计算Lj的平均值L#，将Lj的平均值L#作为最佳目标长度；

[0161] 根据公式9计算Rj的平均值R#，将Rj的平均值R#作为最佳目标底面半径；

[0162]    公式8

[0163]    公式9。

[0164] 在本发明优选实施例中，步骤S5具体为：

[0165] 根据公式10计算Lj的加权平均值L*，将Lj的加权平均值L*作为最佳目标长度；

[0166] 根据公式11计算Rj的加权平均值R*，将Rj的加权平均值R*作为最佳目标底面半径；

[0167]    公式10

[0168]    公式11

[0169] 其中，Wj为V组雷达中第j组的权值，Wj根据第j组中两部雷达中的每一雷达的性能参数如测量精度、探测距离、稳定性参数、可靠性参数等确定。

[0170] 在本发明优选实施例中，基于V组雷达中每一组得到的目标进动角，计算上述目标进动角的平均值或加权平均值，作为最佳目标进动角。

[0171] 由此，实现了通过多于两部的雷达提取目标结构参数及目标进动角的技术效果，进一步提高了目标结构参数及目标进动角的提取精度。

[0172] 图3是本发明的进动目标结构参数提取方法的第二流程图；其中，首先建立进动目标模型，之后，

[0173] 第一雷达在视角1下利用最优路径法从目标一维距离像获取目标径向长度信息，并通过GHT(广义Hough变换)得到其目标径向长度曲线。

[0174] 第二雷达在视角2下利用最优路径法从目标一维距离像获取目标径向长度信息，并通过GHT(广义Hough变换)得到其目标径向长度曲线。

[0175] 最后，经过公式推导及步骤S3、S4得到目标进动角及进动轴视线角差值，最终得到目标长度、目标底面半径。

[0176] 本发明提供的基于多雷达的进动目标结构参数提取方法，能够以较小的误差精确提取合作/非合作目标的结构参数及进动角，同时大幅减小提取进动特征所需的观测时间。

[0177] 本领域普通技术人员可以理解实现上述实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成，该程序可以存储于一计算机可读取存储介质中，如：ROM/RAM、磁碟、光盘等。

[0178] 以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以作出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。