阅读说明：本技术 一种小范围重频抖动信号的谐波抑制方法 (Harmonic suppression method for small-range repetition frequency jitter signal ) 是由 田田 韩向清 匡华星 于 2020-06-12 设计创作，主要内容包括：本发明属于雷达数据处理及信号分选技术领域,特别涉及一种小范围重频抖动信号的谐波抑制方法,主要解决小范围重频抖动信号的谐波被误分选为重频固定信号的问题；其实现步骤是：1.遍历输入脉冲缓存并进行到达时间差直方图统计,记录所有过门限的直方图小箱；2.对于每个过门限的小箱k,PRI值标识为p<Sub>k</Sub>,寻找对应PRI值为<Image he="133" wi="125" file="DDA0002536242340000011.GIF" imgContent="drawing" imgFormat="GIF" orientation="portrait" inline="no"></Image>的直方图小箱r,获取小箱r相邻m个小箱范围内的直方图统计值分布并求和；3.依次计算N不同取值的情形,若任一情形下的直方图统计值分布及求和同时满足重频抖动谐波判断条件,则判定小箱k对应的PRI值p<Sub>k</Sub>为重频抖动信号的谐波；本发明有利于提高信号分选准确率,可用于雷达侦测系统。(The invention belongs to the technical field of radar data processing and signal sorting, and particularly relates to a harmonic suppression method for a small-range repetition frequency jitter signal, which mainly solves the problem that the harmonic of the small-range repetition frequency jitter signal is mistakenly sorted into a repetition frequency fixed signal; the method comprises the following implementation steps: 1. traversing the input pulse buffer, carrying out arrival time difference histogram statistics, and recording all histogram bins which pass through a threshold; 2. for each bin k that crosses the threshold, the PRI value is identified as p k Find a corresponding PRI value of The histogram bin r, a histogram system for obtaining the range of m adjacent bins of the bin rEvaluating distribution and summing; 3. sequentially calculating the conditions of different values of N, and if the distribution and summation of the histogram statistical values in any condition simultaneously meet the judgment condition of the repetition frequency jitter harmonic wave, judging the PRI value p corresponding to the small box k k Is a harmonic of the re-frequency jittered signal; the invention is beneficial to improving the signal sorting accuracy and can be used for a radar detection system.)

一种小范围重频抖动信号的谐波抑制方法

技术领域

本发明属于雷达数据处理及信号分选技术领域。

背景技术

重频抖动信号是辐射源进行复杂时域调制的一种重要表现形式，其时域特征可表示如下：t_n＝t_n-1+p(1±ε_n)n＝1,2,...,N-1，其中，p是抖动中心，ε是抖动半径。当ε取值较小时，例如小于5％，则称之为小范围抖动信号。重频抖动信号的信号分选是一个难题。实际环境中，由于信号类型混杂，且重频固定信号的存在会降低抖动信号分选准确率，因此通常采用先分选重频固定信号再分选重频抖动信号的分选流程，即对于输入的混合信号PDW数据，先进行重频固定信号的分选处理，剩余PDW数据再进行重频抖动信号分选处理。

但是，对于小范围重频抖动信号，由于其脉冲重复间隔的变化范围较小，例如只有几十微秒甚至十几微妙，导致其二次或者更高次谐波在一段较短时间内有一定概率与重频固定信号具有类似的时域特征，从而被误分选为重频固定信号，降低了抖动信号分选正确率。现有谐波抑制方法，通常只能应对原始信号及其谐波属于相同信号类型的情况。对于小范围重频抖动信号，由于其谐波的信号类型可能发生改变，因此并不能很好地抑制。如何抑制小范围重频抖动信号的谐波及其导致的误分选是一个亟待解决的问题。

发明内容

本发明的目的在于提供一种小范围重频抖动信号的谐波抑制方法，对通过直方图统计处理生成的重频固定信号进行谐波判断，剔除可能是小范围重频抖动信号谐波的PRI值，提高信号分选准确率。

为实现上述技术目的，本发明采用如下步骤予以实现：

步骤1：遍历输入脉冲缓存，对任意两个脉冲计算到达时间差并进行直方图统计，对于每个过门限的直方图小箱k，记录其直方图统计值v_k以及对应的PRI值p_k；

步骤2：对于每个过门限的小箱k，寻找直方图小箱r，使其对应PRI值为

其中N为正整数且N_min≤N≤N_max，获取以小箱r为中心[r-m,r+m]范围内所有小箱的直方图统计值，m为正整数，并计算统计值的和v_sum；

步骤3：对于过门限小箱k，依次计算N不同取值的情形，获取对应的直方图小箱r，若任一情形下同时满足以下两个条件：(1)v_sum*β＞v_k，其中0＜β≤1；(2)小箱r为中心[r-m,r+m]范围内的所有小箱中存在q个小箱的直方图统计值大于等于阈值v_thr，其中q≥2，v_thr≥1；则判定过门限小箱k对应的PRI值p_k为重频抖动信号的谐波。

本发明的有益效果为：针对小范围重频抖动信号的谐波在信号分选时可能被误分选为重频固定信号的问题，对通过直方图统计处理生成的重频固定信号进行谐波判断，剔除可能是小范围重频抖动信号谐波的PRI值，提高信号分选准确率。

附图说明

图1为本发明实施例流程图。

具体实施方式

图1为本发明实施例的流程图。结合附图，本发明实施例的步骤为：

步骤1：遍历输入脉冲缓存，对任意两个脉冲计算到达时间差并进行直方图统计，对于每个过门限的直方图小箱k，记录其直方图统计值v_k以及对应的PRI值p_k；

(1.1)对于任意两个脉冲，计算其到达时间差dtoa以及直方图统计时所属的小箱b，

其中τ_min是信号分选算法关注的PRI范围的最小值，x是小箱宽度；将小箱b的直方图统计值加一；

(1.2)遍历输入脉冲缓存，重复上述步骤直到遍历结束；

(1.3)对于小箱b，计算门限值

其中α是权系数，0＜α＜1，T是观察时间，其值等于输入脉冲缓存中任意两个脉冲到达时间差的最大值，τ是小箱b对应的PRI值，

(1.4)记录所有直方图统计值大于门限值的直方图小箱，对于每个过门限的直方图小箱k，直方图统计值用v_k表示，对应的PRI值p_k表示，

步骤2：对于每个过门限的小箱k，寻找直方图小箱r，使其对应PRI值为其中N为正整数且N_min≤N≤N_max，获取以小箱r为中心[r-m,r+m]范围内所有小箱的直方图统计值，m为正整数，并计算统计值的和v_sum；

(2.1)对于每个过门限的直方图小箱k，设其对应的PRI值p_k，寻找对应PRI值为的直方图小箱r，

其中τ_min是信号分选算法关注的PRI范围的最小值，x是小箱宽度；

(2.2)查看以小箱r为中心[r-m,r+m]范围内的所有直方图小箱，将上述范围内各个小箱的直方图统计值用v_i表示，计算统计值的和

步骤3：对于过门限小箱k，依次计算N不同取值的情形，获取对应的直方图小箱r，若任一情形下同时满足以下两个条件，则判定小箱k对应的PRI值p_k为重频抖动信号的谐波：

(1)v_sum*β＞v_k，其中0＜β≤1；

(2)小箱r为中心[r-m,r+m]范围内的所有小箱中存在q个小箱的直方图统计值大于等于阈值v_thr，其中q≥2，v_thr≥1，则判定过门限小箱k对应的PRI值p_k为重频抖动信号的谐波。