基于部分高分辨距离像的旋转微动目标参数估计方法
【专利摘要】本发明基于部分高分辨距离像的旋转微动目标参数估计方法，采取的技术方案如下：第一步，雷达高分辨距离像矩阵获取；第二步，选取部分高分辨距离像；第三步，旋转微多普勒频率估计；第四步，目标旋转直径估计。本发明的有益效果在于：第一，不需要雷达对目标进行长时间观测，提高了雷达目标识别的效率；第二，能对被部分干扰的目标进行识别，提高了雷达目标识别的可靠性；第三，不需要进行大量复杂运算，仅仅利用解曲线方程的方式进行数学计算，计算量较小且过程较为便捷。
【专利说明】基于部分高分辨距离像的旋转微动目标参数估计方法 【技术领域】
[0001 ]本发明属于雷达目标识别领域，具体涉及雷达微动目标特征提取技术，特别是一 种利用部分高分辨距离像估计旋转微动目标参数的方法。 【【背景技术】】
[0002] 微动是雷达探测领域一种常见现象，雷达微动目标回波包含大量目标结构和运动 的细节特征，微动参数估计也成为雷达目标识别研究的热点。旋转微动是一种典型的目标 微动形式，常见于直升机旋翼、涡轮发动机叶片等，其特征参数主要包括旋转微多普勒频 率、旋转直径(或半径)等。
[0003] 微动目标参数估计通常有RCS序列方法、时频分析方法和高分辨距离像方法等，这 些参数估计方法共同特点是要求观测时间内包含一个或者多个微动周期，要求雷达观测的 时间要足够长。导致的缺陷主要表现在:一是由于雷达观测实时性的要求，雷达目标回波可 能并不包含一个完整的目标微动周期;二是由于噪声等因素的影响，可能导致微动回波分 量部分被干扰，只有部分雷达目标信号可用于目标识别。 【
【发明内容】
】
[0004] 本发明要解决的技术问题是:利用雷达目标部分高分辨距离像完成对旋转微动目 标参数的可靠估计，本方法不仅实时性好而且估计精度高，结果可靠。
[0005] 本发明采取的技术方案如下：
[0006] 第一步，雷达高分辨距离像矩阵获取
[0007] 利用频率步进信号合成宽带雷达发射信号s(m，n)，me [1，M]是第m个频率步进信 号，Μ是总的信号个数，ne[l，N]是单个频率步进信号中第η个子脉冲，N是总的子脉冲数。对 目标回波信号进行傅里叶变换得到雷达目标高分辨距离像，将各次雷达回波信号对应的高 分辨距离像按行排列形成二维雷达高分辨距离像矩阵H A(m，n)。
[0008] 第二步，选取部分高分辨距离像
[0009] 选取HA(m，n)中未被干扰的部分，然后取该部分旋转目标距离像上振幅最大的两 个旋转点对应的距离像曲线L A(m)和LB(m)，计算两曲线差值绝对值对应的曲线U (m) = | La (m)-LB(m)|，| · I表示取绝对值。
[0010] 第三步，旋转微多普勒频率估计
[0011]选取曲线1^〇11)上任意三点，以旋转微多普勒频率为未知数联立方程组，解方程并 计算得到旋转微多普勒频率估计值β。
[0012] 第四步，目标旋转直径估计
[0013] 利用第三步中估计得到的S，再选取ftLiOii)上任意二点，以旋转直径为未知数联 立方程组，解方程并计算得到目标旋转直径估计值F。
[0014] 本发明的有益效果主要包括：
[0015] 第一，不需要雷达对目标进行长时间观测，提高了雷达目标识别的效率；
[0016] 第二，能对被部分干扰的目标进行识别，提高了雷达目标识别的可靠性；
[0017] 第三，不需要进行大量复杂运算，仅仅利用解曲线方程的方式进行数学计
[0018] 算，计算量较小且过程较为便捷。 【【附图说明】】
[0019] 图1是旋转微动目标示意图。
[0020] 图2是完整目标高分辨距离像示意图。
[0021 ]图3是选取的部分标高分辨距离像回波示意图。
[0022]图4是旋转微多普勒频率估计结果。
[0023]图5是目标旋转直径估计结果。 【【具体实施方式】】
[0024]下面结合附图对本发明作进一步的说明。步骤如下：
[0025]第一步，获取高分辨距离像对应的曲线。如图1所示，设旋转微动目标上振幅最大 的两个点为A和B，其各自对应的高分辨距离像曲线表达式可表示为：
[0027]其中Ro是旋转中心0距离雷达的距离，γα和rB为A和B两点对应的旋转半径，f 〇是雷 达信号s(m，n)的起始频率，Af是信号s(m，n)的步进值，TR是脉冲重复周期，ω为旋转微多 普勒频率，(Κ= ωπιΝΤ??+Φο为旋转角度，Φο是初始角度
为常 数。
[0028]第二步，计算部分距离像对应的差值绝对值曲线。从图2中可见，虽然旋转微动目 标对应的高分辨距离像包含多个微动周期，但是由于大部分被噪声干扰，常规微动参数估 计方法面临困难。选取部分距离像曲线如图3所示，然后计算两曲线差值的绝对值：
[0029] Li(m)= |LA(m)_LB(m) | =r*sin( ω*ηι+Φ〇-ε)(公式二）
[0031]第三步，估计旋转微多普勒频率。以ρ为间隔，从U (m)中取三个点：U (mo-p)，U (mo)，Li(m〇+p)，由下式计算每一个p值对应下的ω*(ρ):
[0033]其中_为m=2,3, · · ·，Μ-1 中任意一点：
.，pe [1 ,Μ-mo]。再求得：
[0035]然后对ω (p)求和再取平均，可得旋转微多普勒率估计值为：
[0037]旋转微多普勒率ω(ρ)随着p的变化如图4所示，最后求得旋转微多普勒率估计值 S =32.14 (rad/S) ?与真实值30 · 79 (rad/s)相对误差为4 · 2 %。
[0038] 第四步估计旋转目标直径。从LKm)中取二个点^(nnAm)，由下式计算每 一个Am值对应下的r( Am):
[0039]
[0040] 其中，mi为m=l，2,3, · · ·，M-1中任意一点，Am=l，2,3, · · ·，M-mi。然后对r( Am)求 和再取平均，可得旋转目标直径估计值为：
[0042]旋转目标直径r( △ m)随着△ m的变化如图5所示，最后求得旋转目标直径估计值 F = 15.05m，与真实值14.63m相对误差为2.9%。
【主权项】
1. 一种基于部分高分辨距离像的旋转微动目标参数估计方法，其特征在于，包括如下 步骤： 步骤一:雷达高分辨距离像矩阵获取 利用频率步进信号合成宽带雷达发射信号3(111，11)，111￡[1，1]是第111个频率步进信号，1 是总的信号个数，ne[l，N]是单个频率步进信号中第η个子脉冲，N是总的子脉冲数;对目标 回波信号进行傅里叶变换得到雷达目标高分辨距离像，将各次雷达回波信号对应的高分辨 距离像按行排列形成二维雷达高分辨距离像矩阵化(m，n); 步骤二:选取部分高分辨距离像 选取化(m，n)中未被干扰的部分，然后取该部分旋转目标距离像上振幅最大的两个旋转 点对应的距离像曲线LA(m)和LB(m)，计算两曲线差值绝对值对应的曲线^(m)=|LA(m)-LB (m) I，I · I表示取绝对值； 步骤Ξ:旋转微多普勒频率估计 选取曲线^如)上任意Ξ点，W旋转微多普勒频率为未知数联立方程组，解方程并计算 得到旋转微多普勒频率估计值^ ; 步骤四：目标旋转直径估计 利用第Ξ步中估计得到的0)，再选取曲^如)上任意二点，W旋转直径为未知数联立方 程组，解方程并计算得到目标旋转直径估计值f。2. 根据权利要求1所述的基于部分高分辨距离像的旋转微动目标参数估计方法，其特 征在于:在步骤二中，设旋转微动目标上振幅最大的两个点为A和B，其各自对应的高分辨距 离像曲线表达式表示为：其中Ro是旋转中屯、0距离雷达的距离，Μ和η为A和B两点对应的旋转半径，时是雷达信号 s(m，n)的起始频率，Af是信号s(m，n)的步进值，Tr是脉冲重复周期，ω为旋转微多普勒频 率，Φη= ωπιΝΤκ+Φ〇为旋转角度，Φο是初始角虔为常数。3. 根据权利要求1所述的基于部分高分辨距离像的旋转微动目标参数估计方法，其特 征在于:在步骤二中，选取部分距离像曲线，然后计算两曲线差值的绝对值：(公式二） 其中，r = rA+rB为目标旋转直径，ω*= coNTr。4. 根据权利要求1所述的基于部分高分辨距离像的旋转微动目标参数估计方法，其特 征在于：在步骤Ξ中，Wp为间隔，从b(m)中取Ξ个点：Li(m〇-p)，b(m〇)，b(m〇+p)，由下式计 算每一个P值对应下的ω^ρ): V ·* /5.根据权利要求1所述的基于部分高分辨距离像的旋转微动目标参数估计方法，其特 征在于:在步骤四中，从b(m)中取二个点:Li(mi)，b(mi+ Am)，由下式计算每一个Am值对应 下的!( Am):其中，mi为m=l,2,3,. . . ,Μ-1 中任意一点，Am=l,2,3,. . . 然后对!( Am)求和再 取平均，得到旋转目标直径估计值为：
【文档编号】G01S7/41GK106093900SQ201610373341
【公开日】2016年11月9日
【申请日】2016年5月30日
【发明人】潘小义, 王伟, 冯德军, 傅其祥, 顾赵宇, 刘永才, 徐乐涛, 张文明
【申请人】中国人民解放军国防科学技术大学