CN111781564B - Rapid prejudging method for jitter signal pulse sequence - Google Patents
Rapid prejudging method for jitter signal pulse sequence Download PDFInfo
- Publication number
- CN111781564B CN111781564B CN202010533434.3A CN202010533434A CN111781564B CN 111781564 B CN111781564 B CN 111781564B CN 202010533434 A CN202010533434 A CN 202010533434A CN 111781564 B CN111781564 B CN 111781564B
- Authority
- CN
- China
- Prior art keywords
- jitter
- pulse sequence
- sequence
- signal
- max
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/285—Receivers
- G01S7/292—Extracting wanted echo-signals
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention belongs to the technical field of radar data processing and signal sorting, and particularly relates to a rapid prejudging method for a jitter signal pulse sequence, which solves the problem of false signals caused by the false entering of the jitter signal sequence into a conventional sorting module by carrying out sequence type prejudging before signal sorting; the method comprises the following implementation steps: (1) Performing first-level difference histogram statistics of arrival time differences on the input pulse sequence; (2) Judging whether the grid width is reduced or not and carrying out primary difference histogram statistics again; (3) Performing a first conditional decision of the shaking pulse sequence according to the statistical result; (4) Finding out all pulses which meet the jitter search condition in the pulse sequence according to the first judgment result; 5. and executing secondary condition judgment of the jitter pulse sequence and outputting a final sequence type pre-judgment result. The invention is beneficial to reducing false signals and improving the sorting accuracy of the jitter signals, and can be used for a radar detection system.
Description
Technical Field
The invention belongs to the technical field of radar data processing and signal sorting, and particularly relates to a rapid prejudging method for a jitter signal pulse sequence.
Background
In signal sorting, for conventional signals, including PRI fixed signals and spread signals, histogram statistics are usually performed by using histogram bins with fixed widths. In contrast, the repetition frequency jitter signal needs to be counted by using a histogram grid with overlapping ranges. The two processing flows are different, and engineering implementation is generally divided into two different sorting modules. In an actual environment, because the existence of the conventional signal can reduce the sorting accuracy of the jitter signal, a sorting flow of sorting the conventional signal first and then sorting the jitter signal is generally adopted, that is, the input pulse buffer firstly enters the conventional signal sorting module and then enters the jitter signal sorting module. However, the jittered signal pulse, especially the jittered signal with a small jitter range, may have similar time domain characteristics as the regular signal in a short time, so that there is a small probability that a false regular signal will be generated when entering the regular signal sorting module. The existing signal sorting method can only reduce the occurrence probability of the situations by adopting targeted processing measures such as increasing sequence searching condition limitation and the like, and cannot fundamentally solve the problems.
Therefore, the type of the input pulse buffer sequence needs to be pre-judged before signal sorting, if the jitter pulse sequence is judged, the conventional sorting module is skipped over to directly enter the jitter sorting module, and the false signal problem caused by the fact that the jitter signal sequence enters the conventional sorting module is avoided. Moreover, due to the requirement of high real-time performance of signal sorting, the sequence type pre-judging method is not suitable for consuming more computing resources, and how to realize rapid pre-judging is also the key point.
Disclosure of Invention
The invention aims to provide a method for quickly pre-judging a jitter signal pulse sequence, which is used for pre-judging the type of an input pulse buffer sequence before signal sorting, directly introducing pulses judged to be the jitter pulse sequence into a jitter signal sorting module, reducing false signals and improving the sorting accuracy of the jitter signals.
In order to achieve the technical purpose, the invention comprises the following steps:
step 1: the arrival time difference primary difference histogram statistics is carried out on the input pulse sequence, and the grid width is set to be x 1 Sorting the statistical values from big to small, wherein the statistical value set and the subscript set of the first N lattices are respectively marked as N topN And idx topN N is a positive integer, and the statistic of the maximum statistic grid is marked as N max Corresponding PRI identification max ；
Step 2: calculating idx topN The difference value between the subscript of the middle and maximum grids and the subscript of the minimum grid is less than or equal to a set threshold n th Then, the width of the grid is reduced to x 2 ，x 2 ＜x 1 And carrying out arrival time difference primary difference histogram statistics on the input sequence again, and updating n topN 、idx topN 、pri max 、n max ；
And step 3: computing
If r is satisfied 1 > 0.5 or r 2 If any condition is more than 0.4, the sequence is judged not to be a shaking pulse sequence, wherein N total Is the total number of pulses contained in the input pulse sequence; if not, entering step 4;
and 4, step 4: traversing the pulse sequence, calculating the arrival time difference dtoa of each pulse and the previous pulse, and finding out the pulse sequence satisfying the jitter search condition of | dtoa-pri max |≤pri max *ratio j *2, the number of which is designated n jitter Wherein ratio j Is the maximum jitter radius of the supported jitter signal;
and 5: calculating n found jitter Mean arrival time difference of pulses dTAA mean Sum standard deviation dTAA std Calculating
Judging whether the condition is satisfied
If the signal is satisfied, the signal is judged to be a jitter pulse sequence, and if the signal is not satisfied, the signal is judged to be a non-jitter sequence.
The invention has the beneficial effects that: aiming at the problem that the conventional signal is possibly generated by missorting when a small-range jitter signal enters a conventional sorting module, the method for quickly prejudging the jitter signal pulse sequence is provided.
Drawings
FIG. 1 is a flow chart of an embodiment of the present invention.
Detailed Description
FIG. 1 is a flow chart of an embodiment of the present invention. With reference to the accompanying drawings, the specific steps of this embodiment are:
step 1: the arrival time difference of the input pulse sequence is subjected to one-level difference histogram statistics, namely, the first-level difference is that only the arrival time difference of adjacent pulses is calculated, and the grid width is set as x during the statistics 1 ；
After statistics is completed, the statistics values are sorted from big to small, and the statistics value set of the first N grids is marked as N topN The subscript set of the first N lattices is identified as idx topN N is a positive integer; the statistic of the maximum grid is marked as n max Corresponding PRI identification max ；
Step 2: compute set idx topN The difference value between the subscript of the middle and maximum grids and the subscript of the minimum grid is less than or equal to a set threshold n th The width of the grid is reduced to x 2 ，x 2 ＜x 1 And carrying out arrival time difference primary difference histogram statistics on the input pulse sequence again, and updating n topN 、idx topN 、pri max 、n max ；
And step 3: computing
If r is satisfied 1 > 0.5 or r 2 If any condition is more than 0.4, the sequence is judged not to be a shaking pulse sequence, wherein N total Is the total number of pulses contained in the input pulse sequence; if not, entering step 4;
(3.1) calculating r 1 If r is 1 If the pulse sequence is more than 0.5, judging that the sequence is not a shaking pulse sequence;
(3.2) calculating r 2 If r is 2 If the pulse sequence is more than 0.4, judging that the sequence is not a shaking pulse sequence;
(3.3) if neither of the conditions in 3.1 and 3.2 are met, entering step 4;
and 4, step 4: traversing the pulse sequence, calculating the arrival time difference dtoa of each pulse and the previous pulse, and finding out the pulse sequence satisfying the jitter search condition of | dtoa-pri max |≤pri max *ratio j *2, the number of which is designated n jitter Wherein ratio j Is the maximum jitter radius of the supported jitter signal;
and 5: calculating the searched n jitter Mean arrival time difference of pulses dTAA mean Sum standard deviation dTAA std Calculating
Judging whether the condition is satisfied
If the signal is satisfied, the signal is judged to be a jitter pulse sequence, and if the signal is not satisfied, the signal is judged to be a non-jitter sequence.
(5.1) calculating the searched n jitter Mean difference of arrival (dTAA) of a sequence of shaking pulses mean Sum standard deviation dTAA std ；
(5.2) calculating r 3 ；
And (5.3) judging whether the conditions shown in the formula 4 are met, if so, judging the sequence to be a jitter pulse sequence, and if not, judging the sequence to be a non-jitter sequence.
Claims (1)
1. A method for quickly prejudging a jitter signal pulse sequence is characterized by comprising the following steps:
step 1: the arrival time difference primary difference histogram statistics is carried out on the input pulse sequence, and the grid width is set as x 1 Sorting the statistical values from big to small, and respectively identifying the statistical value set and the subscript set of the first N lattices as N topN And idx topN N is a positive integer, and the statistic of the maximum statistic grid is marked as N max Corresponding PRI identification max ；
Step 2: calculating idx topN The difference value between the subscript of the middle and maximum grids and the subscript of the minimum grid is less than or equal to a set threshold n th The width of the grid is reduced to x 2 ，x 2 <x 1 And carrying out arrival time difference primary difference histogram statistics on the input pulse sequence again, and updating n topN 、idx topN 、pri max 、n max ；
And step 3: computing
If r is satisfied 1 >0.5 or r 2 >0.4, it is determined that the input pulse train is not a dither pulseBurst sequence, wherein N total Is the total number of pulses contained in the input pulse sequence; if not, entering step 4;
and 4, step 4: traversing the input pulse sequence, calculating the arrival time difference dtoa of each pulse and the previous pulse, and finding out the pulse meeting the jitter search condition | dtoa-pri max |≤pri max *ratio j *2, the number of which is denoted n jitter Wherein ratio j Is the maximum jitter radius of the supported jitter signal;
and 5: calculating n found jitter Mean arrival time difference of pulses dTAA mean Sum standard deviation dTAA std Calculating
Judging whether the condition is satisfied
If the two sequences are satisfied, the sequence is determined to be a jitter pulse sequence, and if the two sequences are not satisfied, the sequence is determined to be a non-jitter sequence.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010533434.3A CN111781564B (en) | 2020-06-12 | 2020-06-12 | Rapid prejudging method for jitter signal pulse sequence |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010533434.3A CN111781564B (en) | 2020-06-12 | 2020-06-12 | Rapid prejudging method for jitter signal pulse sequence |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111781564A CN111781564A (en) | 2020-10-16 |
| CN111781564B true CN111781564B (en) | 2022-11-18 |
Family
ID=72756240
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010533434.3A Active CN111781564B (en) | 2020-06-12 | 2020-06-12 | Rapid prejudging method for jitter signal pulse sequence |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111781564B (en) |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3559205B2 (en) * | 1999-09-06 | 2004-08-25 | 三菱電機株式会社 | Pulse train separation apparatus and pulse train separation method |
| US6832172B2 (en) * | 2001-06-15 | 2004-12-14 | Tektronix, Inc. | Apparatus and method for spectrum analysis-based serial data jitter measurement |
| CN106295500B (en) * | 2016-07-22 | 2019-06-18 | 中国船舶重工集团公司第七二四研究所 | A kind of repetition dither signal and normal signal separation method |
| CN106842149B (en) * | 2017-03-31 | 2019-08-16 | 中国科学院电子学研究所 | With the improved sequential difference histogram method for separating of Clustering |
| JP2018179548A (en) * | 2017-04-04 | 2018-11-15 | 株式会社東芝 | Jitter pulse train analyzer and jitter pulse train analysis method |
| CN108181613B (en) * | 2017-12-04 | 2020-07-14 | 山东航天电子技术研究所 | Non-uniform quantitative sorting method for sequence difference values of PRI (pulse repetition index) jittering signals |
| CN110764063B (en) * | 2019-10-15 | 2024-04-30 | 哈尔滨工程大学 | Radar signal sorting method based on combination of SDIF and PRI transformation method |
-
2020
- 2020-06-12 CN CN202010533434.3A patent/CN111781564B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN111781564A (en) | 2020-10-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110764063B (en) | Radar signal sorting method based on combination of SDIF and PRI transformation method | |
| CN109270497B (en) | Multidimensional parameter pre-sorting method for radar pulse signals | |
| Meinshausen et al. | Estimating the proportion of false null hypotheses among a large number of independently tested hypotheses | |
| US20210203605A1 (en) | Method and apparatus for detecting abnormal traffic pattern | |
| CN108226869B (en) | Detection method of radar signal with multiple frequency spread | |
| US20100283666A1 (en) | Radar signals clustering method using frequency modulation characteristics and combination characteristics of signals, and system for receiving and processing radar signals using the same | |
| CN114417943B (en) | Identification method of repetition frequency modulation type | |
| CN112799053B (en) | Real-time pairing and sorting method for monopulse signals | |
| CN108181613B (en) | Non-uniform quantitative sorting method for sequence difference values of PRI (pulse repetition index) jittering signals | |
| CN112986928A (en) | Signal sorting multi-source fusion processing method in complex electromagnetic environment | |
| CN109727428B (en) | Repeated alarm suppression method based on deep learning | |
| CN110932807A (en) | Spectrum sensing method of MIMO (multiple input multiple output) system under non-Gaussian noise | |
| CN110929842A (en) | Accurate intelligent detection method for burst time region of non-cooperative radio signal | |
| CN110596652B (en) | FPGA-based zero intermediate frequency secondary radar response pulse extraction method | |
| CN111781564B (en) | Rapid prejudging method for jitter signal pulse sequence | |
| US7916814B2 (en) | Method and apparatus for real-time pulse parameter estimator | |
| CN110632588A (en) | Zero intermediate frequency secondary radar direct current offset compensation algorithm based on FPGA | |
| CN113030932B (en) | Robust adaptive detection method and system for extended target | |
| US11846662B2 (en) | Pulse train signal cycle estimation device, pulse train signal cycle estimation method, and pulse train signal cycle estimation program | |
| CN106877901B (en) | A kind of detection method of low noise than direct sequence signal | |
| CN112198481A (en) | PRI jitter radar signal sorting method under pulse loss aliasing condition | |
| CN115980689A (en) | Point cloud detection-based radiation source signal sorting method, device, equipment and medium | |
| KR101468560B1 (en) | A method of detecting outlier pulses of radar signal based on mahalanobis distance | |
| CN110426696B (en) | Pulse defect radar signal characteristic sequence searching method | |
| Guillén et al. | Two-dimensional determination of the decision boundary for a radar detection method in the moment space |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |