其中，首先建立进动目 标丰吴型，之后，
[0173] 第一雷达在视角1下利用最优路径法从目标一维距离像获取目标径向长度信息， 并通过GHT (广义Hough变换)得到其目标径向长度曲线。
[0174] 第二雷达在视角2下利用最优路径法从目标一维距离像获取目标径向长度信息， 并通过GHT (广义Hough变换)得到其目标径向长度曲线。
[0175] 最后，经过公式推导及步骤S3、S4得到目标进动角及进动轴视线角差值，最终得到 目标长度、目标底面半径。
[0176] 本发明提供的基于多雷达的进动目标结构参数提取方法，能够以较小的误差精确 提取合作/非合作目标的结构参数及进动角，同时大幅减小提取进动特征所需的观测时间。
[0177] 本领域普通技术人员可以理解实现上述实施例方法中的全部或部分步骤是可以 通过程序来指令相关的硬件来完成，该程序可以存储于一计算机可读取存储介质中，如： R0M/RAM、磁碟、光盘等。
[0178] 以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人 员来说，在不脱离本发明原理的前提下，还可以作出若干改进和润饰，这些改进和润饰也应 视为本发明的保护范围。
【主权项】
1. 一种进动目标结构参数提取方法，其特征在于，包括：51. 通过第一雷达的目标一维距离像获取第一雷达目标径向长度信息；通过第二雷达 的目标一维距离像获取第二雷达目标径向长度信息；52. 根据第一雷达目标径向长度信息计算第一雷达目标径向长度曲线;根据第二雷达 目标径向长度信息计算第二雷达目标径向长度曲线；53. 利用第一雷达目标径向长度曲线、及第二雷达目标径向长度曲线计算目标进动角、 及第一雷达与第二雷达的进动轴视线角差值；54. 基于目标进动角、进动轴视线角差值、第一雷达目标径向长度曲线及第二雷达目标 径向长度曲线计算目标长度及目标底面半径； 所述进动轴视线角为雷达视线与目标进动轴的夹角。2. 如权利要求1所述的方法，其特征在于，步骤S3具体包括：531. 利用第一雷达目标径向长度曲线、第二雷达目标径向长度曲线及公式1计算目标 进动角；532. 利用第一雷达目标径向长度曲线、第二雷达目标径向长度曲线、目标进动角及公 式2计算第一雷达与第二雷达的进动轴视线角差值；其中，Θ为目标进动角;Xi为第一雷达在第一时间段1^内，目标径向长度的最大值与最小 值的和;Yi为第一雷达在第一时间段!^内，目标径向长度的最大值与最小值的差;X2为第二 雷达在第二时间段Τ2内，目标径向长度的最大值与最小值的和;Υ2为第二雷达在第二时间段 Τ2内，目标径向长度的最大值与最小值的差;为第一雷达的进动轴视线角，且~从扩到90° 按照第一步进值设置为与A对应的第二雷达的进动轴视线角；i>l、且ieN。3. 如权利要求2所述的方法，其特征在于，步骤S4具体包括： S41.基于目标进动角、进动轴视线角差值及公式3、公式4计算对应于每个第一雷达进 动轴视线角<^的目标样本长度U及目标样本底面半径心；S42.根据目标样本长度U及目标样本底面半径Ri提取目标长度及目标底面半径； 其中，Τη…仏为!^内的N个时间点，Τ2γ··Τ2μ为T2内的Μ个时间点;M>1，N>1，MEN，NEN; Len(Tii)…Len(TiN)为与T11…Tin对应的目标径向长度，Len(T2i)…Len(T2M)为与T21…T2M对 应的目标径向长度；ω为目标进动角速度；φβι、φ β2分别为第一雷达、第二雷达视角下的目 标进动初相角；_[11'"1:川为与1'11'"1'1〃对应的第一雷达视线角，Τ21···τ2Μ为与Τ21···Τ2Μ对应的第 二雷达视线角;且4. 如权利要求3所述的方法，其特征在于，步骤S42具体包括： S421.根据公式5计算每一组目标样本长度心及目标样本底面半径U的误差值e(ai);S422 .比较所有误差值的大小，并提取最小的误差值对应的目标样本长度作为目标长 度，提取最小的误差值对应的目标样本底面半径作为目标底面半径。5. 如权利要求4所述的方法，其特征在于，步骤S2具体包括：521. 根据第一雷达目标径向长度信息，确定第一雷达目标径向长度曲线的参数ln、l12、 和的取值范围；根据第二雷达目标径向长度信息，确定第二雷达目标径向长度曲线的参数 121、 l22、q>2的取值范围；522. 基于111、112、91的取值范围及精度要求，设置分别与111、1 12、中1对应的取值间隔， 对1 η、112、φ>:进行离散取值，生成1 η、112、％参数数组;并 基于121、122、φ2的取值范围及精度要求，设置分别与121、1 22、φ2对应的取值间隔，对121、 122、 和进行离散取值，生成121、122、Φ_2参数数组；523. 利用第一雷达获取的目标径向长度信息、及公式6检验In Jmq)!参数数组;根据 检测结果确定第一雷达目标径向长度曲线的参数;并 利用第二雷达获取的目标径向长度信息、及公式7检验121、122、卟参数数组;根据检测 结果确定第二雷达目标径向长度曲线的参数；所述第一雷达获取的目标径向长度信息包括时间点tn…t1P及与该时间点对应的目标 径向长度Len(tii)…Len(tip); 所述第二雷达获取的目标径向长度信息包括时间点加…加及与该时间点对应的目标 径向长度 Len(t2i)."Len(t2Q);P>l，Q>l，PEN，QeN〇6. 如权利要求5所述的方法，其特征在于，步骤S23具体包括： 5231. 选取In、112、((?参数数组中任一数组，将tn···t1P及该任一数组代入公式6，得到P 个运算结果;分别计算P个运算结果与Len(tn)…Len(t1P)的差值，得到111、1 12、中1参数数组 中该任一数组的P个差值;并 选取l21、l22、q>2参数数组中任一数组，将t2r-t 2Q&该任一数组代入公式7,得到Q个运算 结果;分别计算Q个运算结果与Len(t21)…Len(t2Q)的差值，得到121、1 22、(内参数数组中该任 一数组的Q个差值； 5232. 比较1 η、112、恥参数数组中任一数组的Ρ个差值与第一阈值的大小，统计Ρ个差值 中小于第一阈值的差值数量，并将该差值数量作为 1η、1?2、φι参数数组中该任一数组的合 格数;并 比较121、122、φ?参数数组中任一数组的Q个差值与第二阈值的大小，统计Q个差值中小 于第二阈值的差值数量，并将该差值数量作为121、122、恥参数数组中该任一数组的合格数； 5233. 选取111、112、(卩1参数数组中合格数最大的数组，将该数组的参数作为第一雷达目 标径向长度曲线的参数;并 选取121、122、92参数数组中合格数最大的数组，将该数组的参数作为第二雷达目标径 向长度曲线的参数； 所述第一阈值根据第一雷达精度要求设置； 所述第二阈值根据第二雷达精度要求设置。7. 如权利要求1-6任一所述的方法，其特征在于，步骤S1具体为： 通过第一雷达的目标一维距离像，利用最优路径法获取第一雷达目标径向长度信息； 通过第二雷达的目标一维距离像，利用最优路径法获取第二雷达目标径向长度信息。8. 如权利要求7所述的方法，其特征在于，在步骤S1之前，所述方法还包括： SO.从U部雷达中任取两部合为一组，形成V组雷达;根据所述V组雷达中任一组中的两 部雷达提取目标长度k及目标底面半径Rj; 在步骤S4之后，所述方法还包括： S5.根据目标长度k及目标底面半径&计算最佳目标长度及最佳目标底面半径;其中，U >2,且UEN; ; 1 < j < V，且jEN。9. 如权利要求8所述的方法，其特征在于，步骤S5具体为： 根据公式8计算k的平均值L#，将k的平均值1^#作为最佳目标长度； 根据公式9计算吣的平均值1?#，将心的平均值浐作为最佳目标底面半径；10.如权利要求8所述的方法，其特征在于，步骤S5具体为： 根据公式10计算k的加权平均值ΙΛ将k的加权平均值f作为最佳目标长度； 根据公式11计算Rj的加权平均值If，将Rj的加权平均值浐作为最佳目标底面半径；其中，W」为V组雷达中第j组雷达的权值，W」由第j组中两部雷达中的每一雷达的性能参 数确定。
【专利摘要】本发明公开了一种进动目标结构参数提取方法，其特征在于，包括：分别通过两部雷达的目标一维距离像获取两部雷达的目标径向长度信息；根据获取的目标径向长度信息计算两部雷达的目标径向长度曲线；利用两部雷达的目标径向长度曲线计算目标进动角、及两部雷达的进动轴视线角差值；基于目标进动角、两部雷达的进动轴视线角差值、及两部雷达的目标径向长度曲线计算目标长度及目标底面半径。与现有技术相比，本发明能够以较小的误差精确提取合作/非合作目标的结构参数及进动角，同时大幅减小提取进动特征所需的观测时间。
【IPC分类】G01S7/41
【公开号】CN105676200
【申请号】CN201610009296
【发明人】唐传子, 任红梅, 盛晶
【申请人】北京环境特性研究所
【公开日】2016年6月15日
【申请日】2016年1月7日