基于HRRP序列的空间进动目标参数估计方法

查 林，陈大庆，吴 鹏

(1．中国电子科技集团公司第三十八研究所，安徽合肥 230088；2．太原卫星发射中心技术部，山西太原 036301；3.中国人民解放军32086部队，北京 100000)

0 引言

中段防御是导弹防御的关键阶段，因为它具有防御面积大、拦截时间长、弹道相对固定[1]的优点。然而当导弹在中段释放诱饵后，由于弹头、诱饵、残骸、假目标等组成包含真假目标的威胁目标群，并且以大致相同的速度和轨迹在大气层外做惯性飞行，很难用轨迹和速度上的差别来区分各目标；此外，其他各种先进突防手段的采用，也使得中段拦截导弹在目标识别上面临极大的挑战，这已经成为弹道导弹防御系统发展的一个瓶颈。但是由于弹头目标的特殊结构在特定的受力作用下会产生微动，所谓微动是指目标或目标的组成部分除了质心平动以外的振动、转动、翻滚和加速运动等微小运动。因为此运动状态常常是独一无二的，反映了目标的精细特征，所以可以用来作为目标识别的重要依据。因此，很多学者提出利用中段目标的微动特征来区分真假弹头[2-13]，微动参数估计的主要途径[2]有基于时频图[3]、ISAR(Inverse Synthetic Aperture Radar)图像[4-5]、RCS(Radar Cross Section)[6-7]、高分辨距离像[8-14]等。

雷达高分辨距离像包含了目标径向的空间结构信息[15]，在雷达光学区内能够较好地描述目标的散射特性，在弹道中段，因受弹头进动的影响[16-17]，由多幅连续高分辨距离像构成的一维距离像序列上的散射中心位置将会发生周期性走动，其包含着目标的进动调制信息，因此可利用目标的一维距离像序列来估计目标的进动参数和形状参数。文献[8-13]以锥体目标为研究对象，其中文献[8]在利用目标一维距离像序列变化估计进动参数过程中需要已知目标的一些结构参数，而这些参数实际中是很难获取的；文献[9-10]利用散射中心在雷达视线上投影距离之间的相关性实现目标进动特征的提取；文献[11]对目标散射中心的一维投影距离进行建模，采用非线性拟合的方法估计进动和结构参数；文献[12-13]通过雷达视角的变化来实现目标进动特征的提取，其中文献[12]从时间-距离像数据实现锥体弹头散射点三维空间位置的重构，并在此基础上进行目标进动特征的提取，文献[13]利用多视角观测一维距离像序列联合获得弹道目标各个散射中心的绝对坐标；文献[14]以旋转对称锥柱体目标为研究对象，基于静态电磁散射数据，利用雷达观测视角内锥柱体各散射中心的一维距离像序列中各散射中心间的相对位置变化的极值与目标参数之间的关系来进行目标进动和结构参数的估计，不过仅从区域Ⅱ(α<β<π/2)进行考虑，且在参数估计过程中未考虑参数耦合问题。

本文从宽带雷达体制出发，以锥体目标为研究对象，针对缺少目标结构参数等先验信息及估计过程中参数耦合问题，基于高分辨距离测量提出了一种能同时估计出空间锥体目标进动参数和形状参数的方法，仿真表明该方法可行、有效。文中第二部分对锥体目标的进动模型进行分析，得出目标姿态角随时间变化的关系式，进而得出目标径向长度序列与姿态角关系式；第三部分对如何从回波信号中提取径向长度序列作了详细介绍；第四部分给出了本方法的具体实现步骤及推导过程；第五部分通过仿真计算验证了该方法的可行性和有效性。

1 锥体目标进动模型

弹道导弹的弹头外形一般较常见的有平底锥弹头、球底锥弹头、平底锥柱弹头及球底锥柱弹头4种，本文主要对平底锥弹头进行分析，其结构如图1所示，其中α为目标半锥角，b为目标斜边长度，R为底面半径，1、2、3表示弹头目标的3个散射中心。

图1 锥体目标结构图

在弹道导弹飞行过程中，为了保证飞行的平稳性和命中的精度，弹头目标会采用自旋的方式定向，而在诱饵释放、弹箭分离时，弹头目标不可避免地会受到横向扰动，进而产生进动，称为微动。由于诱饵、假目标、残骸等其他空间目标无姿态控制系统，故不会具有规律的微动特性。弹头目标的规律微动特性导致姿态角，即雷达视线与目标自旋轴的夹角，呈现周期变化。姿态角的变化规律同目标的进动参数有关，但是姿态角一般不易直接测量得到。另一方面，如果不考虑目标的平动，姿态角的周期变化会导致空间进动锥体目标的径向长度序列呈现规律变化，其变化规律同锥体目标的形状参数和进动参数有关，其中，某次观测中锥体目标的径向长度序列，是指通过该次观测中的各次回波计算出的锥顶散射中心和锥底散射中心的径向距离的差值，并由这些差值组成的序列。因此可以利用弹头目标特有的微动特性来提取弹头目标的进动参数和形状参数。

根据以上分析，可建立弹道中段锥体目标进动模型[9]如图2所示，Oxcyczc坐标系为平动坐标系，Oxsyszs坐标系为随体坐标系，弹头绕其对称轴Ozs以角速度Ω作自旋运动，同时轴Ozs绕轴Ozc以角速度ω进动，进动角为θ，故zs为自旋轴，zc为进动轴，γ为雷达视线与进动轴夹角，即俯仰角，ω为进动角频率，β为雷达视线与自旋轴夹角，即姿态角，LOS是指雷达视线方向。根据图2所示的锥体目标的进动模型，可以得到第k次观测锥体目标的姿态角随时间t变化的表达式为

βk(t)=arccos[cosθkcosγk+

sinθksinγkcos(ωt+φk)]

(1)

式中，γk为第k次观测的俯仰角，θk为第k次观测的进动角，φk为第k次观测的进动初相角，k=1,2。

图2 锥体目标进动模型

实际雷达探测目标时一般迎着锥顶照射锥体目标，且锥形弹头和诱饵的半锥角较小，一般取 5°～15°，此时可认为锥体目标进动时姿态角β< π-α，在这一角度范围内，钝头锥对应散射中心1与2均一直可见[9]。根据上述姿态角随时间变化的表达式，可以求出第k次观测中锥体目标的径向长度序列的表达式为

Lk(n)=bcos(βk(tkn)+α)

(2)

式中，tkn表示第k次观测第n个回波的时刻，n=1,2,…,N，N为每次观测的回波数目。

2 径向长度序列提取

当目标运动速度不高时，可以认为在一个脉冲持续时间内，目标到雷达站的距离不变，即所谓的“停-跳”假设；但对于高速机动目标，如导弹目标，上述假设不再成立，此时应当首先估计每次观测时目标的轨道运动信息。令第k次观测中第i个回波时，目标到雷达站距离的表达式为

Rki(t)=Rki0+fki(t)

(3)

式中，Rki0为参考时刻目标距雷达站的距离，fki(t)为第k次观测第i个回波中目标的轨道运动表达式。下面对所获得的宽带回波信号进行平动补偿，具体方法如下：

令第k次观测中第i个回波对应的雷达发射信号的表达式为

ski(t)=pki(t)ej2πf0t=|pki(t)|ejφki(t)ej2πf0t

(4)

式中，k=1,2，i=1,2,…,N，f0为载频，c为电磁波的传播速度，pki(t)为发射信号的复包络，|pki(t)|为pki(t)的幅度，φki(t)为pki(t)的相位，|·|表示求模运算。则该次回波信号解调后的表达式为

srki(t)=Apki(t-τki(t))e-j2πf0τki(t)=

A|pki(t-τki(t))|ejφki(t-τki(t))e-j2πf0τki(t)

(5)

式中，τki(t)=2Rki(t)/c，A为该次回波信号的幅度。将目标到雷达站的距离表达式代入式(5)，则该次回波的表达式可以写作：

srki(t)=Apki(t-τki(t))e-j2πf0τki(t)=

A|pki(t-τki(t))|ejφki(t-τki(t))·

e-j2πf0τki(t)=

A|pki(t-τki(t))|ejφki1(t)·

(6)

(7)

至此，完成了对宽带回波信号的平动补偿，下面依次对各信号用MUSIC算法[16]作参数估计，提取强散射中心的距离，并形成散射中心距离历程图(一维距离像序列)，其中，散射中心距离历程图是指横坐标为各次回波的时刻，纵坐标为各次回波中强散射中心的距离的图像。然后运用图像边缘检测方法[17]，从散射中心距离历程图提取各次回波时刻锥顶散射中心和锥底散射中心的径向距离的差值，这些差值组成的序列即为锥体目标的径向长度序列，并将第一次观测的径向长度序列记为L1(n)，第二次观测的径向长度序列记为L2(n)，其中，n=1,2,…,N。

3 参数估计

l1max=bcos[(γ1+α)-θ]=

bcos(γ1+α)cosθ+bsin(γ1+α)sinθ

(8)

l1min=bcos[(γ1+α)+θ]=

bcos(γ1+α)cosθ-bsin(γ1+α)sinθ

(9)

l2max=bcos[(γ2+α)-θ]=

bcos(γ2+α)cosθ+bsin(γ2+α)sinθ

(10)

l2min=bcos[(γ2+α)+θ]=

bcos(γ2+α)cosθ-bsin(γ2+α)sinθ

(11)

进一步有以下关系成立：

l1max+l1min=2bcos(γ1+α)cosθ

(12)

l1max-l1min=2bsin(γ1+α)sinθ

(13)

l2max+l2min=2bcos(γ2+α)cosθ

(14)

l2max-l2min=2bsin(γ2+α)sinθ

(15)

为进一步推导方便，令a11=l1max+l1min，a12=l1max-l1min，a21=l2max+l2min，a22=l2max-l2min。所以有

(a11sinθ)2+(a12cosθ)2=4b2sin2θcos2θ

(16)

(a21sinθ)2+(a22cosθ)2=4b2sin2θcos2θ

(17)

因此可以解得

(18)

(19)

(20)

l1max=bcos[(θ+α)-γ1]=

bcos(θ+α)cosγ1+bsin(θ+α)sinγ1

(21)

l1min=bcos[(γ1+α)+θ]=

bcos(θ+α)cosγ1-bsin(θ+α)sinγ1

(22)

l2max=bcos[(θ+α)-γ2]=

bcos(θ+α)cosγ2+bsin(θ+α)sinγ2

(23)

l2min=bcos[(θ+α)+γ2]=

bcos(θ+α)cosγ2-bsin(θ+α)sinγ2

(24)

进一步有以下关系成立：

l1max+l1min=2bcos(θ+α)cosγ1

(25)

l1max-l1min=2bsin(θ+α)sinγ1

(26)

l2max+l2min=2bcos(θ+α)cosγ2

(27)

l2max-l2min=2bsin(θ+α)sinγ2

(28)

为进一步推导方便，令a11=l1max+l1min，a12=l1max-l1min，a21=l2max+l2min，a22=l2max-l2min。所以有

[a11sin(θ+α)]2+[a12cos(θ+α)]2=

4b2sin2(θ+α)cos2(θ+α)

(29)

[a21sin(θ+α)]2+[a22cos(θ+α)]2=

4b2sin2(θ+α)cos2(θ+α)

(30)

因此可以解得

(31)

(32)

(33)

(34)

(35)

(36)

(37)

令x=cosωt，则有

(38)

以下来分析h的变化规律：

1) 当x=±1，即ωt=kπ时，h取最小值0。

2) 下面分析h何时取得极大值。根据求极值原理可知当h对x的一阶偏导为0时h可以取得极大值。

(39)

(40)

至此，我们完成了对目标进动角、半锥角及斜边长度的估计，其计算流程如图3所示。

4 仿真与分析

实验所观测的锥体目标如图1所示，其形状参数为：目标半锥角α=14°，目标斜边长度b=1 m，底面半径R=0.25 m；进动参数为：进动角θ=10°，进动角频率ω=31.4 rad/s。第一次观测的俯仰角γ1=37°，第二次观测的俯仰角γ2=37.2°。实验所用数据为空间进动锥体目标电磁仿真数据，测量频带为9～11 GHz，频点间隔为0.02 GHz。假定目标回波中的观测噪声为高斯白噪声，观测信噪比SNR的定义为

(41)

图3 方法流程图

首先在观测信噪比20 dB条件下进行仿真实验，得如下结果：图4(a)为第一次观测所得的目标锥顶、锥底散射中心的一维距离像序列，图4(b)为第二次观测所得的目标锥顶、锥底散射中心的一维距离像序列，图4(c)、(d)分别为第一、二次观测所得的锥体目标径向长度序列。从图中可以看出，散射点对应的一维距离像序列与径向长度序列均呈正弦规律变化，且频率一致，这与前述分析一致，故可通过这种规律性变化来进行目标形状和进动参数估计；另外，从图中还可以看出此时所得的径向长度序列曲线有毛刺，并且信噪比越低，毛刺越多，即抖动越大，这是因为信噪比较低时，距离估计精度比较差。图5以第一次观测为例，给出了不同信噪比条件下目标径向长度序列曲线。

(a) 第一次观测得到的散射中心一维距离像序列

(b) 第二次观测得到的散射中心一维距离像序列

(c) 第一次观测得到的径向长度序列

(d) 第二次观测得到的径向长度序列图4 锥体目标一维距离像及径向长度序列

(a) SNR=15 dB

(b) SNR=20 dB

(c) SNR=25 dB图5 不同信噪比下的径向长度序列曲线

(a) 第一次观测得到的平滑滤波后的径向长度距离曲线

(b) 第二次观测得到的平滑滤波后的径向长度距离曲线图6 平滑滤波后的径向长度序列曲线

(a) 锥体目标径向长度微分曲线

(b) 平滑滤波后的锥体目标径向长度微分曲线图7 锥体目标径向长度微分曲线

图9 半锥角估计误差

图10 斜边长度的估计误差

5 结束语

本文从宽带雷达体制出发，以空间锥体目标为研究对象，针对缺少目标结构参数等先验信息及估计过程中参数耦合问题，通过对其进动模型的分析，根据锥顶和锥底散射中心径向距离变化曲线，利用点散射中心间相对位置的极值变化信息，提出一种基于高分辨距离测量的中段弹道导弹目标进动参数和形状参数估计方法，并给出了详细推导步骤，仿真结果表明此方法可以同时估计出空间进动锥体目标的形状参数和进动参数；另外，该方法充分利用了空间进动锥体目标的径向长度序列、进动参数和形状参数的关系，且不使用非线性曲线拟合的方法估计目标的进动参数，从而提高了参数估计方法的可实施性，改善了参数估计过程的稳定性。