基于部分高分辨距离像的旋转微动目标参数估计方法技术

本发明专利技术基于部分高分辨距离像的旋转微动目标参数估计方法，采取的技术方案如下：第一步，雷达高分辨距离像矩阵获取；第二步，选取部分高分辨距离像；第三步，旋转微多普勒频率估计；第四步，目标旋转直径估计。本发明专利技术的有益效果在于：第一，不需要雷达对目标进行长时间观测，提高了雷达目标识别的效率；第二，能对被部分干扰的目标进行识别，提高了雷达目标识别的可靠性；第三，不需要进行大量复杂运算，仅仅利用解曲线方程的方式进行数学计算，计算量较小且过程较为便捷。

【技术实现步骤摘要】

本专利技术属于雷达目标识别领域，具体涉及雷达微动目标特征提取技术，特别是一种利用部分高分辨距离像估计旋转微动目标参数的方法。
技术介绍
微动是雷达探测领域一种常见现象，雷达微动目标回波包含大量目标结构和运动的细节特征，微动参数估计也成为雷达目标识别研究的热点。旋转微动是一种典型的目标微动形式，常见于直升机旋翼、涡轮发动机叶片等，其特征参数主要包括旋转微多普勒频率、旋转直径(或半径)等。微动目标参数估计通常有RCS序列方法、时频分析方法和高分辨距离像方法等，这些参数估计方法共同特点是要求观测时间内包含一个或者多个微动周期，要求雷达观测的时间要足够长。导致的缺陷主要表现在：一是由于雷达观测实时性的要求，雷达目标回波可能并不包含一个完整的目标微动周期；二是由于噪声等因素的影响，可能导致微动回波分量部分被干扰，只有部分雷达目标信号可用于目标识别。
技术实现思路
本专利技术要解决的技术问题是：利用雷达目标部分高分辨距离像完成对旋转微动目标参数的可靠估计，本方法不仅实时性好而且估计精度高，结果可靠。本专利技术采取的技术方案如下：第一步，雷达高分辨距离像矩阵获取利用频率步进信号合成宽带雷达发射信号s(m,n)，m∈[1,M]是第m个频率步进信号，M是总的信号个数，n∈[1,N]是单个频率步进信号中第n个子脉冲，N是总的子脉冲数。对目标回波信号进行傅里叶变换得到雷达目标高分辨距离像，将各次雷达回波信号对应的高分辨距离像按行排列形成二维雷达高分辨距离像矩阵HA(m,n)。第二步，选取部分高分辨距离像选取HA(m,n)中未被干扰的部分，然后取该部分旋转目标距离像上振幅最大的两个旋转点对应的距离像曲线LA(m)和LB(m)，计算两曲线差值绝对值对应的曲线L1(m)＝|LA(m)-LB(m)|，|·|表示取绝对值。第三步，旋转微多普勒频率估计选取曲线L1(m)上任意三点，以旋转微多普勒频率为未知数联立方程组，解方程并计算得到旋转微多普勒频率估计值第四步，目标旋转直径估计利用第三步中估计得到的再选取曲L1(m)上任意二点，以旋转直径为未知数联立方程组，解方程并计算得到目标旋转直径估计值本专利技术的有益效果主要包括：第一，不需要雷达对目标进行长时间观测，提高了雷达目标识别的效率；第二，能对被部分干扰的目标进行识别，提高了雷达目标识别的可靠性；第三，不需要进行大量复杂运算，仅仅利用解曲线方程的方式进行数学计算，计算量较小且过程较为便捷。【附图说明】图1是旋转微动目标示意图。图2是完整目标高分辨距离像示意图。图3是选取的部分标高分辨距离像回波示意图。图4是旋转微多普勒频率估计结果。图5是目标旋转直径估计结果。【具体实施方式】下面结合附图对本专利技术作进一步的说明。步骤如下：第一步，获取高分辨距离像对应的曲线。如图1所示，设旋转微动目标上振幅最大的两个点为A和B，其各自对应的高分辨距离像曲线表达式可表示为：其中R0是旋转中心O距离雷达的距离，rA和rB为A和B两点对应的旋转半径，f0是雷达信号s(m,n)的起始频率，Δf是信号s(m,n)的步进值，TR是脉冲重复周期，ω为旋转微多普勒频率，φm＝ωmNTR+φ0为旋转角度，φ0是初始角度，为常数。第二步，计算部分距离像对应的差值绝对值曲线。从图2中可见，虽然旋转微动目标对应的高分辨距离像包含多个微动周期，但是由于大部分被噪声干扰，常规微动参数估计方法面临困难。选取部分距离像曲线如图3所示，然后计算两曲线差值的绝对值：L1(m)＝|LA(m)-LB(m)|＝r*sin(ω*m+φ0-ε)(公式二)其中，r＝rA+rB为目标旋转直径，ω*＝ωNTR。第三步，估计旋转微多普勒频率。以p为间隔，从L1(m)中取三个点：L1(m0-p)，L1(m0)，L1(m0+p)，由下式计算每一个p值对应下的ω*(p)：其中m0为m＝2,3,...,M-1中任意一点，取p∈[1,M-m0]。再求得：然后对ω(p)求和再取平均，可得旋转微多普勒率估计值为：旋转微多普勒率ω(p)随着p的变化如图4所示，最后求得旋转微多普勒率估计值与真实值30.79(rad/s)相对误差为4.2％。第四步估计旋转目标直径。从L1(m)中取二个点：L1(m1)，L1(m1+Δm)，由下式计算每一个Δm值对应下的r(Δm)：其中，m1为m＝1,2,3,...,M-1中任意一点，Δm＝1,2,3,...,M-m1。然后对r(Δm)求和再取平均，可得旋转目标直径估计值为：旋转目标直径r(Δm)随着Δm的变化如图5所示，最后求得旋转目标直径估计值与真实值14.63m相对误差为2.9％。本文档来自技高网...

【技术保护点】
一种基于部分高分辨距离像的旋转微动目标参数估计方法，其特征在于，包括如下步骤：步骤一：雷达高分辨距离像矩阵获取利用频率步进信号合成宽带雷达发射信号s(m,n)，m∈[1,M]是第m个频率步进信号，M是总的信号个数，n∈[1,N]是单个频率步进信号中第n个子脉冲，N是总的子脉冲数；对目标回波信号进行傅里叶变换得到雷达目标高分辨距离像，将各次雷达回波信号对应的高分辨距离像按行排列形成二维雷达高分辨距离像矩阵HA(m,n)；步骤二：选取部分高分辨距离像选取HA(m,n)中未被干扰的部分，然后取该部分旋转目标距离像上振幅最大的两个旋转点对应的距离像曲线LA(m)和LB(m)，计算两曲线差值绝对值对应的曲线L1(m)＝|LA(m)‑LB(m)|，|·|表示取绝对值；步骤三：旋转微多普勒频率估计选取曲线L1(m)上任意三点，以旋转微多普勒频率为未知数联立方程组，解方程并计算得到旋转微多普勒频率估计值步骤四：目标旋转直径估计利用第三步中估计得到的再选取曲L1(m)上任意二点，以旋转直径为未知数联立方程组，解方程并计算得到目标旋转直径估计值

【技术特征摘要】
1.一种基于部分高分辨距离像的旋转微动目标参数估计方法，其特征在于，包括如下步骤：步骤一：雷达高分辨距离像矩阵获取利用频率步进信号合成宽带雷达发射信号s(m,n)，m∈[1,M]是第m个频率步进信号，M是总的信号个数，n∈[1,N]是单个频率步进信号中第n个子脉冲，N是总的子脉冲数；对目标回波信号进行傅里叶变换得到雷达目标高分辨距离像，将各次雷达回波信号对应的高分辨距离像按行排列形成二维雷达高分辨距离像矩阵HA(m,n)；步骤二：选取部分高分辨距离像选取HA(m,n)中未被干扰的部分，然后取该部分旋转目标距离像上振幅最大的两个旋转点对应的距离像曲线LA(m)和LB(m)，计算两曲线差值绝对值对应的曲线L1(m)＝|LA(m)-LB(m)|，|·|表示取绝对值；步骤三：旋转微多普勒频率估计选取曲线L1(m)上任意三点，以旋转微多普勒频率为未知数联立方程组，解方程并计算得到旋转微多普勒频率估计值步骤四：目标旋转直径估计利用第三步中估计得到的再选取曲L1(m)上任意二点，以旋转直径为未知数联立方程组，解方程并计算得到目标旋转直径估计值2.根据权利要求1所述的基于部分高分辨距离像的旋转微动目标参数估计方法，其特征在于：在步骤二中，设旋转微动目标上振幅最大的两个点为A和B，其各自对应的高分辨距离像曲线表达式表示为：其中R0是旋转中心O距离雷达的距离，rA和rB为A和B...

【专利技术属性】
技术研发人员：潘小义，王伟，冯德军，傅其祥，顾赵宇，刘永才，徐乐涛，张文明，
申请(专利权)人：中国人民解放军国防科学技术大学，
类型：发明
国别省市：湖南;43