CN111551909B - 一种脉间捷变频lfm回波信号跨距离门走动的校正方法 - Google Patents

Abstract

本申请公开了一种脉间捷变频LFM回波信号跨距离门走动的校正方法，该方法包括：步骤1，获取高速运动目标的回波信号，并计算回波信号的滤波信号；步骤2，根据预设替换变量，对滤波信号中的慢时间参数进行替换，对参数替换后的滤波信号进行傅里叶逆变换以实现回波信号跨距离门走动校正。通过本申请中的技术方案，在脉间捷变频情况下，对高速运动目标回波信号的跨距离门走动问题进行数学描述，通过慢时间参数替换实现捷变频回波信号的距离走动校正，以利于后续的信号积累和目标检测，并且这种校正处理能够满足雷达系统的实时性要求。

Description

一种脉间捷变频LFM回波信号跨距离门走动的校正方法

技术领域

本申请涉及雷达信号处理的技术领域，具体而言，涉及一种脉间捷变频LFM回波信号跨距离门走动的校正方法。

背景技术

线性调频信号(linear frequency modulation，LFM)是雷达中广泛应用的多普勒不敏感信号。实际中雷达目标通常都是运动的，尽管目标运动会降低LFM回波脉压的性能指标，使得主瓣展宽、副瓣抬高，但是通常不会降低到难以接受的程度。

对探测高速运动目标而言，帧内的多个脉冲回波信号通常会发生跨距离门走动现象，当进行后续的相参/非相参积累处理时，无法将能量聚集起来，从而降低了目标的可检测性，因此，需要在积累前进行多脉冲回波的跨距离门走动校正。

雷达中采用脉间捷变频技术是对抗敌方干扰的有效措施，脉冲载频的伪随机快速变化增加了敌方捕捉雷达信号、并对雷达信号实施干扰的难度，是一种常用的抗干扰技术。

但是在该技术中，高速运动目标的捷变频LFM脉冲回波的距离走动校正是一个值得研究的问题，但是目前尚未见到这方面的研究。

发明内容

本申请的目的在于：在脉间捷变频LFM脉冲回波中，实现高速运动目标回波信号的脉间跨距离门校正对齐，以利于后续的相参或非相参处理，从而增强了目标的可检测性、雷达的抗干扰能力。

本申请的技术方案是：提供了一种脉间捷变频LFM回波信号跨距离门走动的校正方法，该方法包括：步骤1，获取高速运动目标的回波信号，并计算回波信号的滤波信号；步骤2，根据预设替换变量，对滤波信号中的慢时间参数进行替换，对替换后的滤波信号进行傅里叶逆变换以实现回波信号跨距离门走动校正。

上述任一项技术方案中，进一步地，滤波信号的计算公式为：

t_n＝nT_r

f_dn＝2vf_n/c

f_n＝f_c+d_nΔf

τ_n＝2(R₀+vt_n)/c

式中，f为频率，t_n为慢时间，即第n个脉冲重复周期，n＝0,1，2，...，N-1，N为一帧内的脉冲数量，T_r为脉冲重复周期，a_r为回波信号的幅度，rect(·)为矩形函数，f_dn为运动目标在第n个脉冲中的多普勒频率，v为目标的速度，f_n为第n个脉冲的载频，f_c为捷变频雷达的初始载频，d_n为第n个脉冲的频率调制码，Δf为最小跳频间隔，B为LFM信号的带宽，τ_n为雷达发射第n个脉冲时目标的双程延时，R₀为高速运动目标的初始距离，c为光速，K为LFM信号的调频斜率；

上述任一项技术方案中，进一步地，步骤1中，计算回波信号的滤波信号，具体包括：

步骤11，根据回波信号，计算回波信号对应的混频信号，其中，混频信号的计算公式为：

u(t)＝exp(jπKt²)

τ＝2(R₀+vt_n+vt)/c

式中，t为快时间，

为全时间，T为脉冲宽度，u(t)为LFM信号，τ为回波延时；

步骤12，对混频信号进行快时间域傅里叶变换，记作频域变换信号，其中，快时间域傅里叶变换的计算公式为：

步骤13，计算匹配滤波器的传输函数，根据匹配滤波器原理，以及传输函数和频域变换信号，计算回波信号的滤波信号。

上述任一项技术方案中，进一步地，步骤2中，滤波信号进行傅里叶逆变换的计算公式为：

f_dm＝2vf_m/c

τ₀＝2R₀/c

式中，m为虚拟脉冲序号，t为快时间，

为虚拟慢时间，f_m为第m个虚拟脉冲的载频，f_dm为高速运动目标在第m个虚拟脉冲中的多普勒频率，τ₀为目标初始距离延时。

上述任一项技术方案中，进一步地，预设替换变量为：

式中，t_n为慢时间，n＝0，1，2，…，N-1，f为频率，f_m为第m个虚拟脉冲的载频，

为虚拟慢时间。

上述任一项技术方案中，进一步地，该方法适用于雷达对高速运动目标的检测。

本申请的有益效果是：

本申请中的技术方案，在脉间捷变频情况下，对脉冲间目标回波信号跨距离门走动问题进行数学描述，计算回波信号的滤波信号，并利用所构造的慢时间参数替换公式，对滤波信号中的慢时间参数进行替换，进而通过傅里叶逆变换进行脉冲压缩，实现了目标回波信号的脉间跨距离门校正对齐，增强目标的可检测性和雷达的抗干扰能力，并且这种校正处理能够满足雷达系统的实时性要求。

附图说明

本申请的上述和/或附加方面的优点在结合下面附图对实施例的描述中将变得明显和容易理解，其中：

图1是根据本申请的一个实施例的脉间捷变频LFM回波信号跨距离门走动的校正方法的示意流程图；

图2是根据本申请的一个实施例的载频跳变时目标在脉间跨距离走动的仿真图；

图3是根据本申请的一个实施例的载频跳变时目标在脉间跨距离走动被校正的仿真图；

图4是根据本申请的一个实施例的载频跳变时未校正回波的非相参积累的仿真图；

图5是根据本申请的一个实施例的载频跳变时已校正回波的非相参积累的仿真图。

具体实施方式

为了能够更清楚地理解本申请的上述目的、特征和优点，下面结合附图和具体实施方式对本申请进行进一步的详细描述。需要说明的是，在不冲突的情况下，本申请的实施例及实施例中的特征可以相互结合。

在下面的描述中，阐述了很多具体细节以便于充分理解本申请，但是，本申请还可以采用其他不同于在此描述的其他方式来实施，因此，本申请的保护范围并不受下面公开的具体实施例的限制。

如图1所示，本实施例提供了一种脉间捷变频LFM回波信号跨距离门走动的校正方法，适用于雷达对高速运动目标的检测，该方法包括：

步骤1，获取高速运动目标的回波信号，并计算回波信号的滤波信号；

以对捷变频信号进行跨距离门走动校正为例，对本实施例中的校正方法进行说明。设定捷变频雷达的初始载频为f_c，第n个脉冲的载频为f_n，n＝0，1，2，...，N-1，N为一帧内的脉冲数量，载频f_n的计算公式为：

f_n＝f_c+d_nΔf

式中，d_n为第n个脉冲的频率调制码，是随机整数，d_n∈{0，1，...，D-1}，Δf为最小跳频间隔，则第n个脉冲发射信号的计算公式可以表述为：

u(t)＝exp[jπKt²]

-T/2≤t≤T/2

t_n＝nT_r

rect(t)＝1，-1/2≤t≤1/2

式中，s_t(t，t_n)为发射信号，t为快时间，t_n为慢时间，即第n个脉冲重复周期，

为全时间，T_r为脉冲重复周期，a_t为发射信号的幅度，T为脉冲宽度，rect(·)为矩形函数，K为LFM信号的调频斜率，f_n为第n个脉冲的载频。

因此，可以将远处一个高速运动目标的回波信号表示为：

τ＝2(R₀+vt_n+vt)/c

式中，a_r为回波信号的幅度，τ为回波时延，R₀为目标的初始距离，c为光速，v为目标速度。

进一步的，上述步骤1中，具体包括：

步骤11，根据所述回波信号，计算所述回波信号对应的混频信号。

为了便于后续的视频处理，将回波信号混频到零中频，回波信号混频后的计算公式为：

τ_n＝2(R₀+vt_n)/c

式中，τ_n为雷达发射第n个脉冲时目标的双程延时。

步骤12，对混频信号进行快时间域傅里叶变换，记作频域变换信号，其中，快时间域傅里叶变换的计算公式为：

步骤13，计算匹配滤波器的传输函数U^*(f)，上标*表示共轭运算，根据匹配滤波器原理，以及传输函数U^*(f)和频域变换信号，计算得到回波信号的滤波信号，即匹配滤波器输出，其中，滤波信号的计算公式为：

t_n＝nT_r

f_dn＝2vf_n/c

f_n＝f_c+d_nΔf

τ_n＝2(R₀+vt_n)/c

式中，f为频率，t_n为慢时间，n＝0，1，2，...，N-1，N为一帧内的脉冲数量，T_r为脉冲重复周期，a_r为回波信号的幅度，rect(·)为矩形函数，f_dn为高速运动目标在第n个脉冲中的多普勒频率，v为目标的速度，f_n为第n个脉冲的载频，f_c为捷变频雷达的初始载频，d_n为第n个脉冲的频率调制码，Δf为最小跳频间隔，B为LFM信号的带宽，τ_n为雷达发射第n个脉冲时目标的双程延时，R₀为目标的初始距离，c为光速，K为LFM信号的调频斜率。

对上述得到的滤波信号X(f，t_n)执行傅里叶逆变换，实现脉冲压缩，得到：

通过对上述计算公式进行分析可见，脉压后目标回波的峰值位于t＝τ_n+f_dn/K处，该峰值位置包含了两项，分别为：

1、回波时延τ_n＝2(R₀+vt_n)/c，当一帧内目标运动距离满足vNT_r＞c/2B时，出现跨距离门走动现象，需要校正对齐；

2、多普勒耦合项f_dn/K＝2v(f_c+d_nΔf)/c/K，这是线性调频信号的固有缺陷，会导致峰值发生偏移，且这里偏移是变化的，有可能超出一个距离门，需要避免。

步骤2，根据预设替换变量，对滤波信号中的慢时间参数进行替换，对替换后的滤波信号进行傅里叶逆变换，其中，预设替换变量为：

式中，t_n为第n个脉冲对应的慢时间，n＝0，1，2，…，N-1，f为频率，f_m为第m个虚拟脉冲的载频，

为虚拟慢时间。

通过理论分析和模拟仿真发现，可以通过第m个虚拟脉冲的载频f_m与频率f，将第n个脉冲对应的慢时间t_n替换为虚拟慢时间

以实现距离走动校正。因此，设定预设替换变量为：

之后，利用预设替换变量替换滤波信号X(f，t_n)中的慢时间t_n，替换后的滤波信号计算公式为：

f_dm＝2vf_m/c

式中，f_dm为运动目标在第m个虚拟脉冲中的多普勒频率，τ₀为目标的初始双程时延。

对上述滤波信号

执行傅里叶逆变换，以实现脉冲压缩，脉冲压缩后的信号为：

f_dm＝2vf_m/c

τ₀＝2R₀/c

式中，m为虚拟脉冲序号，t为快时间，

为虚拟慢时间，f_m为第m个虚拟脉冲的载频，f_dm为运动目标在第m个虚拟脉冲中的多普勒频率，τ₀为目标初始距离延时。

通过对上述计算公式进行分析可以得出，目标位于t＝τ₀+f_dm/κ，不再跨距离门走动，从而实现了距离走动校正。

为了验证本实施例中校正方法的有效性，进行仿真验证，设定如下参数：初始载频f₀＝10GHz，一帧内的脉冲个数N＝64，D＝32，频率调制码d_n∈{0，1，...，31}，最小跳频间隔Δf＝30MHz，B＝50MHz，脉冲宽度T＝50μs，脉冲重复周期T_r＝300μs，目标速度v＝2000m/s，目标的初始距离R₀＝240km，单脉冲脉压前信噪比SNB＝-14dB，捷变频信号的载频跳变范围为D×Δf＝960MHz。

首先在上述参数下仿真目标的回波信号，然后利用本实施例中的校正方法对回波信号进行校正，仿真图依次如图2和图3所示。

由仿真图可见，校正前，64个脉冲的回波显示为斜线，不在一个距离刻度上，在脉间跨越了大约13个距离门，同距离刻度相加效果差。但是，通过本实施例中的方法进行校正之后，64个脉冲回波显示在一条横线上，都校正到了初始距离上，表示同距离刻度，可以完美实现所有脉冲的同距离刻度叠加。

距离校正前后的非相参积累如图4和图5所示，可见校正后的目标无明显距离扩展，信噪比显著提高，可检测性显著改善。

以上结合附图详细说明了本申请的技术方案，本申请提出了一种脉间捷变频LFM回波信号跨距离门走动的校正方法，包括：步骤1，获取高速运动目标的回波信号，并计算回波信号的滤波信号；步骤2，根据预设替换变量，对滤波信号中的慢时间参数进行替换，对替换后的滤波信号进行傅里叶逆变换以实现回波信号跨距离门走动校正。通过本申请中的技术方案，在脉间捷变频情况下，对高速运动目标回波信号的跨距离门走动问题进行数学描述，通过替换慢时间参数实现捷变频回波信号的距离走动校正，以利于后续的信号积累和目标检测，并且这种校正处理能够满足雷达系统的实时性要求。

本申请中的步骤可根据实际需求进行顺序调整、合并和删减。

本申请装置中的单元可根据实际需求进行合并、划分和删减。

尽管参考附图详地公开了本申请，但应理解的是，这些描述仅仅是示例性的，并非用来限制本申请的应用。本申请的保护范围由附加权利要求限定，并可包括在不脱离本申请保护范围和精神的情况下针对发明所作的各种变型、改型及等效方案。

Claims (5)

1.一种脉间捷变频LFM回波信号跨距离门走动的校正方法，其特征在于，该方法包括：

步骤1，获取高速运动目标的回波信号，并计算所述回波信号的滤波信号；其中，所述滤波信号的计算公式为：

t_n＝nT_r

f_dn＝2vf_n/c

f_n＝f_c+d_nΔf

τ_n＝2(R₀+vt_n)/c

式中，f为频率，t_n为慢时间，即第n个脉冲重复周期，n＝0，1，2，...，N-1，N为一帧内的脉冲数量，T_r为脉冲重复周期，a_r为回波信号的幅度，rect(·)为矩形函数，f_dn为运动目标在第n个脉冲中的多普勒频率，v为目标的速度，f_n为第n个脉冲的载频，f_c为捷变频雷达的初始载频，d_n为第n个脉冲的频率调制码，Δf为最小跳频间隔，B为LFM信号的带宽，τ_n为雷达发射第n个脉冲时目标的双程延时，R₀为目标的初始距离，c为光速，K为LFM信号的调频斜率；

步骤2，根据预设替换变量，对所述滤波信号中的慢时间参数进行替换，对参数替换后的滤波信号进行傅里叶逆变换以实现回波信号跨距离门走动校正。

2.如权利要求1所述的脉间捷变频LFM回波信号跨距离门走动的校正方法，其特征在于，步骤1中，计算所述回波信号的滤波信号，具体包括：

步骤11，根据所述回波信号，计算所述回波信号对应的混频信号，其中，所述混频信号的计算公式为：

u(t)＝exp(jπKt²)

τ＝2(R₀+vt_n+vt)/c

式中，t为快时间，

为全时间，T为脉冲宽度，u(t)为LFM信号，τ为回波延时；

步骤12，对所述混频信号进行快时间域傅里叶变换，记作频域变换信号，其中，快时间域傅里叶变换的计算公式为：

步骤13，计算匹配滤波器的传输函数，根据匹配滤波器原理，以及所述传输函数和所述频域变换信号，计算所述回波信号的滤波信号。

3.如权利要求1所述的脉间捷变频LFM回波信号跨距离门走动的校正方法，其特征在于，步骤2中，滤波信号进行傅里叶逆变换的计算公式为：

f_dm＝2vf_m/c

τ₀＝2R₀/c

式中，m为虚拟脉冲序号，t为快时间，

为虚拟慢时间，f_m为第m个虚拟脉冲的载频，f_dm为高速运动目标在第m个虚拟脉冲中的多普勒频率，τ₀为目标初始距离延时。

4.如权利要求1至3中任一项所述的脉间捷变频LFM回波信号跨距离门走动的校正方法，其特征在于，所述预设替换变量为：

式中，t_n为慢时间，n＝0，1，2，…，N-1，f为频率，f_m为第m个虚拟脉冲的载频，

为虚拟慢时间。

5.如权利要求1所述的脉间捷变频LFM回波信号跨距离门走动的校正方法，其特征在于，所述方法适用于雷达对所述高速运动目标的检测。

Images