一种针对高阶机动目标雷达回波的检测方法与流程

本发明属于雷达测量技术领域，具体涉及一种针对高阶机动目标雷达回波的检测方法。

背景技术：

因在通信、声纳和电子情报侦察等领域的广泛应用，高速机动目标的长时间积累和参数估计问题越来越受到重视。但是在长时间积累过程中，高速机动目标出现的距离徙动和多普勒扩散问题会引入积累损失，进而影响目标参数估计精度。在1990年《proceedingsoftheieee1990nationalaerospaceandelectronicsconference》第210页至214页，wux等人发表的“simultaneousimagingofmultipletargetsinaninversesyntheticapertureradar”一文中提出了对多目标采用最大似然估计方法进行参数估计和运动补偿，虽然能提供较高的多目标分辨能力，但需要进行多维参数搜索，计算量大。在2012年《ieeetransactionsonsignalprocessing》第60卷第12期第6190页至6201页，xuj等人发表的“radon-fouriertransformforradartargetdetection(i):generalizeddopplerfilterbank”一文中提出了一种新的积累算法—检测前聚焦算法，该算法可同时解决运动目标的距离徙动和多普勒扩散问题，实现相参积累。但是，该算法盲速旁瓣问题突出，严重影响目标的能量积累，并且由于需要对多维参数进行联合搜索，计算量巨大。

技术实现要素：

有鉴于此，本发明的目的是提供一种针对高阶机动目标雷达回波的检测方法，能够对目标的距离徙动和多普勒扩散进行补偿，同时能够解决盲速旁瓣问题并获得精确的参数估计结果。

一种针对高阶机动目标雷达回波的检测方法，包括如下步骤：

步骤1、对雷达回波数据进行脉压处理，具体方法如下：

假设观测场景中存在k个点目标，针对第k个目标，雷达发射线性调频信号表达式为：

其中，k＝1,2,...,k；tp为脉冲宽度，fc为载波频率，γ为调频斜率，τ为快时间，即距离时间，t为慢时间，且t＝nt，n＝0,1,....n-1，n为相参积累时间内发射的脉冲数，t为脉冲重复周期；

雷达接收信号经下变频和距离脉冲压缩处理后变为：

其中，σ0k为目标的反射系数，g为距离压缩增益，b为线性调频信号带宽，c为光速，sinc(x)＝sin(x)/x，λ＝c/fc为波长。

忽略目标的三阶以上运动参数引起的距离徙动和多普勒频谱展宽，雷达到第k个点目标的瞬时斜距为：

其中目标在t0时刻的距离、速度、加速度和加加速度的3阶运动参数为

将式(3)代入式(2)得：

步骤2、首先将式(2)变换到距离频率域，并将其分为两个子带信号，使得其中一个子信号的频域范围为另一个子信号的频域范围为然后分别将这两个信号进行距离向频移δfr/2，得到两个子带信号经过距离向频移后的时域信号分别为：

其中，δfr表示两个子带信号中心频率之间的频率差；g1和g2分别为两子带信号的距离向脉压处理增益，且满足

对这两个子带信号进行共轭相乘，得到合成信号：

步骤3、将代入公式(7)得：

其中，是目标运动轨迹搜索方程，r0、v、a和分别表示目标的初始距离、速度、加速及加加速度；

再把公式(3)代入(8)式得：

则k个目标的脉间多目标信号表达式为：

其中：

为各目标的后向散射系数，为待估计的目标参数，满足ta＝t；

s'r(ta)的参数化瞬时自相关函数定义为：

其中，*表示取共轭；l为时延变量，为与变标操作相关的固定时延；表示瞬时自相关函数交叉相关项；表示瞬时自相关函数的自相关项，表达式为：

对ta进行广义keystone变换，变换表达式为：

其中t'a为变标后的时间变量，η为变标因子；

将式(14)代入式(12)得到：

对公式(15)进行离散傅里叶变换得：

将代入(16)式得：

对式(16-1)采用非均匀傅里叶变换来完成信号能量积累，即：

其中，nufft(·)表示非均匀傅里叶变换算子；

步骤4、针对搜索参数r0、v、a和在各自取值范围内任意取一个值，得到一个参数组合将搜索参数r0、v、a和代入到目标运动轨迹搜索方程中，计算得到r(t)；然后返回步骤3，将r(t)的值代入到(8)式中，执行步骤3后最终得到式(17)的值；通过计算不同的搜索参数组合得到r(t)并返回步骤3，得到各参数组合对应的式(17)的值，则得到的最大值对应的一组参数即为目标的参数值，实现目标检测。

本发明具有如下有益效果：

本发明针对脉冲多普勒体制下变加速运动目标模型，提出了一种针对高阶机动目标雷达回波的检测方法，该方法首先采用子带双频共轭处理降低合成信号的等效频率，然后通过对距离、速度、加速度和加加速度的联合搜索实现对合成信号的数据抽取，再沿抽取的运动轨迹采用调频率-二次调频率变换进行相参积累，该方法能够同时实现距离徙动校正和多普勒扩散补偿，从而获得精确的参数估计结果。

附图说明

图1为子带信号构造示意图；

图2为本发明一种针对高阶机动目标雷达回波的检测方法流程图。

具体实施方式

下面结合附图并举实施例，对本发明进行详细描述。

如图2所示，本发明一种针对高阶机动目标雷达回波的检测方法，具体包括如下步骤：

步骤1、对回波数据进行脉压处理，具体方法如下:

假设观测场景中存在k个点目标。不失一般性，对第k个目标的回波模型进行分析。假设雷达发射线性调频信号为

其中tp为脉冲宽度，fc为载波频率，γ为调频斜率，τ为快时间，即距离时间，t为慢时间，且t＝nt,(n＝0,1,....n-1)，n为相参积累时间内发射的脉冲数，t为脉冲重复周期。

雷达接收信号经下变频和距离脉冲压缩处理后变为

其中σ0k为目标的反射系数，g为距离压缩增益，b为线性调频信号带宽，c为光速，sinc(x)＝sin(x)/x，λ＝c/fc为波长。

由于目标的三阶以上运动参数引起的距离徙动和多普勒频谱展宽均不会跨越一个相应的分辨单元，因此可以忽略。为雷达到第k个点目标的瞬时斜距:

其中目标在t0时刻的距离、速度、加速度和加加速度的3阶运动参数为

将式(3)带入式(2)可得

步骤2:构造两子带信号，并对两子带信号进行子带双频共轭处理，具体方法如下:

首先将式(4)变换到距离频率域并将其分为两个子带信号，如图1所示，然后分别将这两个信号进行距离向频移δfr/2，其中δfr＝|fc1-fc2|表示两个子带信号中心频率之间的频率差。这两个子带信号经过距离向频移后的时域信号分别为

其中g1和g2分别为两子带信号的距离向脉压处理增益，且满足

对这两个子带信号进行共轭相乘，得到合成信号的自项表达式为：

可见，合成信号scomp(t,τ)的等效频率变为δfr，且满足δfr＜＜fc。

步骤3、将步骤2处理后的信号进行调频率-二次调频率变换处理，具体方法如下:

是目标运动轨迹搜索方程，将带入公式(7)可得：

再把公式(3)代入上式可得：

上式可看成多项式相位信号。

由于k个目标同时存在，脉间多目标信号可视为多个多项式相位信号之和，表达式如下

其中：

为各目标的后向散射系数，ta＝t，为待估计的目标参数，满足由于满足δfr＜＜fc，因此可以增大参数估计范围。

s'r(ta)的参数化瞬时自相关函数定义为：

其中，*表示取共轭；l为时延变量，为与变标操作相关的固定时延；

表示瞬时自相关函数交叉相关项；表示瞬时自相关函数的自相关项，表达式为：

ta和l之间存在的耦合，因此能量不能完全积累。为消除ta和l之间存在的耦合，进行广义keystone变换，即其中t'a为变标后的时间变量，η为变标因子。

将式(14)带入式(12)可得：

从式(15)可以看出，原时间变量ta和时延变量l之间的耦合被消除。对公式(15)进行离散傅里叶变换得：

将代入上式得：

从式(16)可以看出，此时的信号能量沿直线分布，沿l'轴即可进行能量的积累处理。由于l'轴数据为非均匀采样数据，为此需采用非均匀傅里叶变换来完成信号能量积累，即：

其中，nufft(·)表示非均匀傅里叶变换算子。

步骤4、针对搜索参数r0、v、a和在各自取值范围内任意取一个值，得到一个参数组合将搜索参数r0、v、a和代入到目标运动轨迹搜索方程中，计算得到r(t)；然后返回步骤3，将r(t)的值代入到(8)式中，执行步骤3后最终得到式(17)的值；通过计算不同的搜索参数组合得到r(t)并返回步骤3，得到各参数组合对应的式(17)的值，则得到的最大值对应的一组参数即为目标的参数值，实现目标检测。

实施例：

本实例中，雷达系统仿真参数如下:载频fc＝10ghz,脉宽tp＝100μs，线性调频带宽b＝150mhz,采样频率fs＝375mhz,脉冲重复频率prf＝1000hz，相参积累脉冲数n＝2048。仿真目标位置与运动参数如下:初始径向距离r01＝80km，r02＝80km，径向速度v1＝600m/s，v2＝550m/s，径向加速度a1＝50m/s²，a2＝60m/s²，径向加加速度雷达接收到的两目标回波的原始信噪比为snr1＝-22db，snr2＝-25db。

本实施例的步骤如下:

步骤1、对回波数据进行脉压处理，得到：

步骤2、首先将式(2)变换到距离频率域，并将其分为两个子带信号，使得其中一个子信号的频域范围为另一个子信号的频域范围为然后分别将这两个信号进行距离向频移δfr/2，得到两个子带信号经过距离向频移后的时域信号分别为：

步骤3、将代入公式(7)得：

再把公式(3)代入(8)式得：

则k个目标的脉间多目标信号表达式为：

其中：

为各目标的后向散射系数，为待估计的目标参数，满足ta＝t；

s'r(ta)的参数化瞬时自相关函数定义为：

表示瞬时自相关函数交叉相关项；表示瞬时自相关函数的自相关项，表达式为：

对ta进行广义keystone变换，变换表达式为：

其中t'a为变标后的时间变量，η为变标因子；

将式(14)代入式(12)得到：

对公式(15)进行离散傅里叶变换得：

将代入(16)式得：

对式(16-1)采用非均匀傅里叶变换来完成信号能量积累，即：

其中，nufft(·)表示非均匀傅里叶变换算子；

步骤4、针对搜索参数r0、v、a和在各自取值范围内任意取一个值，得到一个参数组合将搜索参数r0、v、a和代入到目标运动轨迹搜索方程中，计算得到r(t)；然后返回步骤3，将r(t)的值代入到(8)式中，执行步骤3后最终得到式(17)的值；通过计算不同的搜索参数组合得到r(t)并返回步骤3，得到各参数组合对应的式(17)的值，则得到的最大值对应的一组参数即为目标的参数值，实现目标检测。

本实施例中，经过不断的重复执行步骤3，最后得到2个目标的搜索参数：径向速度径向加速度径向加加速度

综上所述，以上仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。