一种估计空间弹头进动周期的目标微动特征提取方法

文档序号:9908441阅读:731来源:国知局
一种估计空间弹头进动周期的目标微动特征提取方法
【专利说明】一种估计空间弹头进动周期的目标微动特征提取方法 【技术领域】
[0001 ]本发明属于雷达目标探测识别领域,具体涉及雷达目标微动特征分析和提取技 术,特别是一种估计空间弹头进动周期的目标微动特征提取方法。 【【背景技术】】
[0002] 目标微动特征反映了目标的电磁散射特性、几何结构特性以及运动特性等。弹头 在空间飞行过程中要保持稳定,一般均需要采取姿态控制技术,导致弹头目标常常在宏观 运动的同时伴有小角度的进动现象,而伴随弹头飞行的其他空间目标如碎片、诱饵等,由于 没有姿态控制,将会出现随机翻滚、摆动等运动现象。弹头、碎片、诱饵等空间目标在运动类 型、运动状态以及运动参数等上的差异是雷达识别真实弹头的重要依据。
[0003] 弹头进动对雷达回波产生周期性规律调制,表现在频域就是微多普勒频率,在时 频二维成像平面上则对应着相应的周期性曲线。空间弹头目标进动特征提取方法的一般有 两类方法:一类是对时频面上的周期性曲线进行准确建模并求解曲线参数得到弹头进动特 征;另一类是利用图像域一参数域变换处理算法,对变换后的参数域进行穷举搜索提取弹 头进动特征。这些方法主要存在如下问题:第一,空间弹头目标散射点在雷达观测时间内极 易发生越距离单元走动,导致微动曲线在时频成像时出现断续现象,传统基于连续曲线参 数模型的目标特征提取方法精度受到影响;第二,空间进动弹头目标圆环处散射点是典型 滑动型散射中心,进动时频曲线是多条非正弦周期曲线相叠加形式,基于理想正弦曲线模 型的目标特征提取方法误差可能较大;第三,弹头目标雷达回波属于多分量信号,二次型时 频分析方法虽然分辨率高但存在交叉耦合项,传统单分量信号假设模型下的目标特征提取 方法可能难以实用,显然有待改进。 【
【发明内容】

[0004] 本发明要解决的技术问题是:在弹头目标圆环处散射点表现出的非正弦时频曲线 特性情况下,利用循环幅度差矩阵和平均循环幅度差系数完成对弹头目标进动周期的精确 估计。
[0005] 本发明采取的技术方案如下:
[0006] 第一步,弹头进动回波时频二维矩阵获取
[0007] 利用线性时频分析方法(如短时傅里叶分析)对空间弹头目标回波进行时频变换, 以频率序列为行(记作m,me [0,Μ-1],M是总行数),以时间序列为列(记作n,ne [0,Ν-1],N 是总列数)排列成弹头进动回波时频二维矩阵f(m,n)。
[0008] 第二步,计算循环幅度差矩阵
[0009] 从时频二维矩阵f(m,n)中第一列开始,按列计算所选取的列与矩阵f(m,n)中所有 列的循环幅度差,然后取各自最小值,存成一行;依次进行上述计算,得到时频二维矩阵f (m,η)对应的循环幅度差矩阵MA,矩阵大小为N X N。
[0010]第三步,计算平均循环幅度差系数
[0011] 从Ma中第一个对角线开始,依次取出相应对角线全部元素并计算均值,得到包含N 个元素的平均循环幅度差系数4^?),?6[0,^1]是11\中对角线序号。
[0012] 第四步,弹头进动周期估计
[0013] 对六"?)的波谷位置进行搜索,然后对所有搜索到的波谷位置取平均值可得进度 周期估计值。
[0014] 本发明的有益效果主要包括:
[0015] 第一,适应非正弦以及时频曲线断续情况。弹头圆环处散射为滑动型散射,其进动 时频曲线为非正弦曲线,并且由于遮挡效应等因素影响,时频曲线呈现断续现象,而本发明 方法对这些不理想因素具有较好的鲁棒性。
[0016] 第二,适应多曲线叠加耦合。弹头目标常常是多分量信号,各散射点回波对应的时 频曲线在时频面叠加耦合情况比较严重,本发明方法能有效适应多分量信号带来的曲线叠 加耦合效应。
[0017] 第三,目标特征提取结果精确。循环幅度差算法本身具有较好的抗虚假极值点抗 干扰等有利特性,从而进一步保证了弹头进动周期估计结果准确可靠。 【【附图说明】】
[0018] 图1是空间弹头模型示意图。
[0019] 图2是暗室测量所用弹头进动示意图。
[0020] 图3是时频二维矩阵图。
[0021 ]图4是循环幅度差矩阵结果图。
[0022]图5是平均循环幅度差系数结果图。 【【具体实施方式】】
[0023]下面结合附图1-5对本发明作进一步的说明。步骤如下:
[0024] 第一步,对弹头目标回波进行短时傅里叶变换。进动弹头目标圆环处散射表现为 滑动型散射,如图1所示,在坐标系O-Xyz (〇是弹头质心亦为坐标系中心,X是横轴,y是纵轴, z是竖轴)中,弹头以0为支点围绕z轴进动,进动周期为^ γ,fd是进动频率,进动角大小 Id 为Θ。由雷达目标特性可知弹头回波主要由A、B、C、D和E五个强散射点回波构成。依据离散傅 里叶变换形式,计算弹头目标进动回波时频二维矩阵为
[0025]
公式一)其中,m代表频率,η代表时间,1 是时间变量,s[ ·]是弹头目标回波序列,h[ ·]是窗函数,e是指数函数,j是虚数单位。 [0026]图2是实测暗室测量中所用弹头进动测量示意图,进动周期Td ? 3.846s。如图3所 示,空间进动弹头圆环处散射点回波时频曲线是非理想正弦曲线形式,fB是散射点B对应时 频曲线理论值,fc是散射点C对应时频曲线理论值。弹头目标回波是多分量信号,B、C两点时 频曲线相互耦合叠加并出现断续现象,可能导致传统目标特征提取方法出现误差。
[0027] 第二步,从f (m,n)的η列中任意取两列m和Π2,记为f (m,m)和f (m,n2),计算两者的 循环幅度差矩阵。循环幅度差公式为
[0028]
(公式二)其中,ke
[0,M-I ]是频率变量,modM(m+k)表示m+k对M取模后的余数,f (modM(m+k),Π2)表示f (m,Π2)中 第modM(m+k)个元素。
[0029] 取循环幅度差最小值并进行归一化
[0030]
(公式三)其中 mink(A(k;m, n2))表示按k取A (k; m,n2)中所有元素的最小值。
[0031] 按公式二和公式三顺序对f(m,n)中所有列两两计算循环幅度差,可得循环幅度差 矩阵为
[0032] (公式四)其 中,Am(0,0)表示ηι = 0,Π 2 = 0时对应的循环幅度差最小值,Am(0,1)表
当前第1页1 2 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1