阅读说明：本技术 一种单脉冲信号的实时配对分选方法 (Real-time pairing sorting method for single pulse signals ) 是由 赵闯 姜宏志 刘成城 赵拥军 杨静 胡德秀 于 2020-12-09 设计创作，主要内容包括：本发明涉及一种单脉冲信号的实时配对分选方法,属于雷达脉冲信号分选技术领域。本发明首先利用主、副站接收信号的时差范围,构造时差窗,实现对主、副站接收的脉冲之间的粗匹配,缩小了脉冲匹配的范围,降低了分选的计算量；然后联合多参数信息对粗脉冲对序列进行精确配对,改善了以往分选高重频脉冲信号的模糊问题和分选超低重频脉冲信号困难的问题；利用接收同一辐射源同一脉冲的时差不变性对精脉冲对进行分选,保证了分选的正确率。此外本发明仅对主站实时接收的单个脉冲进行分选,克服了以往多站时差分选需要进行时差累积构建时差直方图而导致无法进行实时分选的缺点,为信号实时分选提供了有效的途径。(The invention relates to a real-time pairing sorting method for monopulse signals, and belongs to the technical field of radar pulse signal sorting. The invention firstly utilizes the time difference range of the signals received by the primary station and the secondary station to construct a time difference window, realizes the rough matching between the pulses received by the primary station and the secondary station, reduces the range of pulse matching and reduces the calculation amount of sorting; then, the coarse pulse pair sequences are accurately paired by combining multi-parameter information, so that the problems of the fuzzy problem of sorting high-repetition-frequency pulse signals and the difficulty of sorting ultra-low-repetition-frequency pulse signals in the prior art are solved; the accurate pulse pair is sorted by utilizing the invariance of the time difference of receiving the same pulse of the same radiation source, so that the sorting accuracy is ensured. In addition, the invention only sorts the single pulse received by the main station in real time, overcomes the defect that the multi-station time difference sorting in the prior art cannot carry out real-time sorting because time difference accumulation is needed to construct a time difference histogram, and provides an effective way for the real-time sorting of signals.)

1. A real-time pairing sorting method of single pulse signals is characterized by comprising the following steps:

1) determining a time difference range between the main station and each secondary station, and taking the time difference range as a time difference window;

2) respectively carrying out coarse matching on the single pulse signal received by the main station in real time and the pulse signal stream received by each secondary station in the time difference window by using the determined time difference window to obtain a coarse pulse pair sequence of the main station single pulse matched with each secondary station pulse stream;

3) based on the characteristic that the intra-pulse parameters of the same pulse signal of the same radiation source received by the main station and the secondary station have similarity, the coarse pulse pair sequences between the main station and the secondary station are screened by combining multi-parameter information to obtain the fine pulse pairs of the main station and the secondary stations;

4) and calculating the time difference between the fine pulse pairs of the main station and different secondary stations, taking the time difference between the fine pulse pairs of the main station and the different secondary stations as the time difference pair, and sorting the fine pulse pairs by using the time difference within the error tolerance according to the characteristic that the time difference of the same pulse signal of the same radiation source reaching each station has invariance.

2. The real-time paired sorting method for monopulse signals according to claim 1, wherein the time difference between the primary station and the secondary station in the step 1) is in the range of:

Σ(n_A) N received for the primary station A_AThe range of the time difference corresponding to each pulse,for the nth pulse in the stream received by the primary station A_AThe arrival time of each pulse, σ, represents the time required for the electromagnetic wave to propagate from the secondary station to the primary station A, and μ represents the error caused by the relative position change of the motion of the primary station A and the secondary station in the duration of the detection signal.

3. The method for real-time paired sorting of monopulse signals according to claim 1 or 2, wherein the fine pulse pair screening in step 3) is performed as follows:

A. calculating the parameter distance of each coarse pulse pair in the coarse pulse sequence;

B. calculating multi-parameter similarity matching factors between each coarse pulse pair of the main station and the secondary station according to the reference distance and the parameter distance of each coarse pulse pair;

C. and selecting a coarse pulse pair with a multi-parameter similarity matching factor larger than a set threshold value as a fine pulse pair.

4. The method for real-time paired sorting of monopulse signals according to claim 3, wherein the formula for calculating the parameter distance of the coarse pulse pair in step A is as follows:

wherein the content of the first and second substances,is the nth of the main station A_APulse and sub-station nth_BParameter distance between pulsesFrom the master station A n_APulse and sub-station nth_BThe pulse and a coarse pulse pair, | | | | | represents solving 2-norm, W is a weighting matrix, a symmetric matrix is adopted,is the nth of the main station A_AThe intra-pulse parameters of the individual pulses,is the n-th of the secondary station_BIntra-pulse parameters of individual pulses.

5. The method for real-time paired sorting of monopulse signals according to claim 3, wherein said step B is performed by using a multi-parameter similarity matching factorThe calculation formula of (2) is as follows:

wherein r is₀Is a reference distance representing the minimum multi-parameter distance between the primary station a and secondary station B pulses when uncorrelated; ,is the nth of the main station A_APulse and sub-station nth_BThe parameter distance between pulses, the nth of the master station A_APulse and sub-station nth_BOne pulse pair with one coarse pulse pair.

6. The method for real-time paired sorting of monopulse signals according to claim 3, wherein said set threshold is greater than 0.5 and less than 1.

7. The method for real-time paired sorting of monopulse signals according to claim 1 or 2, wherein the sorting process of step 4) comprises the steps of:

a. forming a time difference pair by the time difference of the fine pulse pair of the main station and two different secondary stations;

b. respectively calculating the distance between the formed time difference pair and the time difference pair corresponding to each sorted pulse sequence, selecting the minimum distance, and if the minimum distance is smaller than the error tolerance, classifying the pulse pair corresponding to the time difference pair into the sorted pulse sequence corresponding to the minimum distance; and if the minimum distance is not less than the error tolerance, a pulse sequence is newly built, and the pulse pair corresponding to the time difference pair is classified into the newly built pulse sequence.

8. The method for real-time paired sorting of monopulse signals according to claim 7, wherein the distance in step b is calculated by the formula:

wherein d is_iFor the distance between the time difference pair corresponding to the pulse sequence under the selected ith pulse source,pulse n received for the master station A_ATime of (t)_CncPulse n received for secondary station C_cAt the time of the day,pulse n received for secondary station B_BTime (DTOA)_iAB,DTOA_iAC) And i is 1,2, …, m is the arithmetic mean value of the corresponding time difference pairs of each pulse pair under each selected radiation source, and m is the number of the radiation sources.

Technical Field

The invention relates to a real-time pairing sorting method for monopulse signals, and belongs to the technical field of radar pulse signal sorting.

Background

The modern battlefield electromagnetic environment is highly complex, on one hand, pulse signals are overlapped in different degrees in time domain, space domain and frequency domain parameters, and the phenomena of batch missing and batch increasing of the signals are more and more serious; on the other hand, the statistical rule of the pulse repetition intervals presents random non-stationarity, which causes great difficulty for the electronic reconnaissance third party to receive; in addition, the signal environment is highly intensive, which increases the amount of calculation for sorting processing, making it difficult to perform sorting in real time. How to correctly separate the pulse signals of each radar from the randomly overlapped pulse streams in real time is a necessary condition for realizing the rapid identification and passive positioning of radar radiation sources, and is an important ring for electronic reconnaissance.

Aiming at the problem of radar signal sorting, the traditional single-station pulse sorting method is based on algorithms of TOA and PRI, and is difficult to effectively and real-timely sort complex pulse signals modulated between pulses; the multi-station time difference histogram sorting method is high in false alarm rate and false alarm rate, the time difference and multi-parameter combined sorting and positioning method relies on positioning calculation, and sorting operation can be performed only by acquiring observation signals with enough duration through all the algorithms, so that instantaneity is poor. Therefore, a new sorting and pairing method is proposed, and a paper named as a sorting and pairing method of the expanded time difference histogram pulses under the constraint criterion (the author: Liu Zhi Xin, Zhao Jun, published in the academic newspaper of the university of electronic science and technology of Western Ann in 2019) discloses a sorting and pairing method of single pulses. Although the accuracy of the sorting by the method of the spread time difference histogram is improved by introducing the multi-parameter constraint criterion, the complexity of the algorithm is increased, and the real-time performance of the algorithm is poor.

Therefore, for the sorting of the complex modulation radar pulse signals, the existing method has the problems of poor sorting precision, poor real-time performance and difficulty in sorting the complex modulation pulse signals.

Disclosure of Invention

The invention aims to provide a real-time pairing sorting method for single pulse signals, which aims to solve the problems of complex algorithm and poor real-time performance of the existing sorting method.

The invention provides a real-time pairing and sorting method of a single pulse signal for solving the technical problems, which comprises the following steps:

1) determining a time difference range between the main station and each secondary station, and taking the time difference range as a time difference window;

2) respectively carrying out coarse matching on the single pulse signal received by the main station in real time and the pulse signal stream received by each secondary station in the time difference window by using the determined time difference window to obtain a coarse pulse pair sequence of the main station single pulse matched with each secondary station pulse stream;

3) based on the characteristic that the intra-pulse parameters of the same pulse signal of the same radiation source received by the main station and the secondary station have similarity, the coarse pulse pair sequences between the main station and the secondary station are screened by combining multi-parameter information to obtain the fine pulse pairs of the main station and the secondary stations;

4) and calculating the time difference between the fine pulse pairs of the main station and different secondary stations, taking the time difference between the fine pulse pairs of the main station and the different secondary stations as the time difference pair, and sorting the fine pulse pairs by using the time difference within the error tolerance according to the characteristic that the time difference of the same pulse signal of the same radiation source reaching each station has invariance.

The invention constructs the time difference window by utilizing the time difference range of the signals received by the primary station and the secondary station, realizes the rough matching between the pulses received by the primary station and the secondary station, reduces the range of pulse matching and reduces the calculation amount of sorting; the coarse pulse pairs are accurately paired by combining multi-parameter information, so that the problems of the fuzzy problem of sorting high repetition frequency pulse signals and the difficulty of sorting ultra-low repetition frequency pulse signals in the prior art are solved; the accurate pulse pair is sorted by utilizing the invariance of the time difference of receiving the same pulse of the same radiation source, so that the sorting accuracy is ensured. The whole sorting process is simple and high in instantaneity.

Further, in order to improve the accuracy of determining the time difference range, the time difference range between the primary station and the secondary station in step 1) is:

Σ(n_A) N received for the primary station A_AThe range of the time difference corresponding to each pulse,for the nth pulse in the stream received by the primary station A_AThe arrival time of each pulse, σ, represents the time required for the electromagnetic wave to propagate from the secondary station to the primary station A, and μ represents the error caused by the relative position change of the motion of the primary station A and the secondary station in the duration of the detection signal.

Further, to ensure accurate registration of the coarse pulse pair, the fine pulse pair screening in step 3) is performed as follows:

A. calculating the parameter distance of each coarse pulse pair in the coarse pulse sequence;

B. calculating multi-parameter similarity matching factors between each coarse pulse pair of the main station and the secondary station according to the reference distance and the parameter distance of each coarse pulse pair;

C. and selecting a coarse pulse pair with a multi-parameter similarity matching factor larger than a set threshold value as a fine pulse pair.

Further, the formula for calculating the parameter distance of the coarse pulse pair in step a is as follows:

wherein the content of the first and second substances,is the nth of the main station A_APulse and sub-station nth_BThe parameter distance between pulses, the nth of the master station A_APulse and sub-station nth_BThe pulse and a coarse pulse pair, | | | | | represents solving 2-norm, W is a weighting matrix, a symmetric matrix is adopted,is the nth of the main station A_AThe intra-pulse parameters of the individual pulses,is the n-th of the secondary station_BIntra-pulse parameters of individual pulses.

Further, the multi-parameter similarity matching factor in the step BThe calculation formula of (2) is as follows:

wherein r is₀Is a reference distance representing the minimum multi-parameter distance between the primary station a and secondary station B pulses when uncorrelated; ,is the nth of the main station A_APulse and sub-station nth_BThe parameter distance between pulses, the nth of the master station A_APulse and sub-station nth_BOne pulse pair with one coarse pulse pair.

Further, the set threshold is greater than 0.5 and less than 1.

Further, the sorting process of the step 4) comprises the following steps:

a. forming a time difference pair by the time difference of the fine pulse pair of the main station and two different secondary stations;

b. respectively calculating the distance between the formed time difference pair and the time difference pair corresponding to each sorted pulse sequence, selecting the minimum distance, and if the minimum distance is smaller than the error tolerance, classifying the pulse pair corresponding to the time difference pair into the sorted pulse sequence corresponding to the minimum distance; and if the minimum distance is not less than the error tolerance, a pulse sequence is newly built, and the pulse pair corresponding to the time difference pair is classified into the newly built pulse sequence.

Further, in order to accurately realize sorting, the distance in step b is calculated by the following formula:

wherein d is_iFor the distance between the time difference pair corresponding to the pulse sequence under the selected ith pulse source,pulse n received for the master station A_ATime of (t)_CncPulse n received for secondary station C_cAt the time of the day,pulse n received for secondary station B_BTime (DTOA)_iAB,DTOA_iAC) And i is 1,2, …, m is the arithmetic mean value of the corresponding time difference pairs of each pulse pair under each selected radiation source, and m is the number of the radiation sources.

Drawings

FIG. 1 is a flow chart of a method for real-time paired sorting of monopulse signals according to the present invention;

FIG. 2 is a schematic diagram of a description of radar pulses received by stations in an embodiment of the present invention;

FIG. 3 is a schematic diagram of a radar pulse stream received by a receiver of each station in an embodiment of the present invention;

FIG. 4 is a schematic diagram of the positions of receivers of each station and a target radar in the embodiment of the present invention;

FIG. 5 shows the nth station A in an embodiment of the present invention_ASchematic time difference window of each pulse;

FIG. 6 is a schematic illustration of sorting based on pairs of time differences in an embodiment of the present invention;

fig. 7 is a schematic diagram of the distribution of pulse-to-time differences for successful matching in an embodiment of the present invention;

FIG. 8-a is a schematic representation of the sorting results of radiation source E1 in an embodiment of the present invention;

FIG. 8-b is a schematic representation of the sorting results of radiation source E2 in an embodiment of the present invention;

FIG. 8-c is a graphical representation of the sorting results of radiation source E3 in an embodiment of the present invention;

FIG. 8-d is a schematic representation of the sorting results of radiation source E4 in an embodiment of the present invention;

FIG. 8-E is a schematic representation of the sorting results of radiation source E5 in an embodiment of the present invention;

FIG. 8-f is a schematic representation of the sorting results of radiation source E6 in an embodiment of the present invention;

FIG. 8-g is a schematic representation of the sorting results of radiation source E7 in an embodiment of the present invention;

fig. 8-h are graphical illustrations of sorting results for radiation source E8 in an embodiment of the present invention.

Detailed Description

The following further describes embodiments of the present invention with reference to the drawings.

The invention aims at signal sorting between a main station and a plurality of secondary stations, firstly, determining a time difference range between the main station and each secondary station, taking the time difference range as a time difference window, and respectively carrying out coarse pairing on a single pulse signal received by the main station in real time and a pulse signal received by each secondary station in the time difference window based on the time difference window to obtain two pulse pair sequences which can be correctly paired with the single pulse of the main station; then, the characteristics that pulse parameters of the same pulse signal of the same radiation source received by the main station and the secondary station have similarity are utilized, pulse pairs in the pulse pair sequence are respectively subjected to pulse accurate matching, two 'pulse pairs' which are truly paired with the single pulse of the main station are obtained, time differences corresponding to the two 'pulse pairs' are calculated, and the two time differences form a 'time difference pair'; the characteristic that the time difference of the same pulse signal of the same radiation source reaching each reconnaissance station has invariance is utilized, and the real-time sorting of the single pulse signals is realized by utilizing the time difference pair within the error tolerance. The implementation flow of the method is shown in fig. 1, and the specific implementation process is as follows.

1. The range of time differences between the primary and each secondary station is determined and used as a window of time differences.

As shown in fig. 4, the multi-station radar receiving system of the present embodiment includes a primary station a, a secondary station B, and a secondary station C, and a target radar D transmits radar pulses, and the pulses received by the stations are as shown in fig. 3. The invention only takes two secondary stations as an example, so the method is also suitable for the situation that more than two secondary stations are arranged in one main station.

Each pulse is represented by a descriptor, each pulse descriptor is shown in figure 2, the master station A, the secondary station B, the secondary station C and the target radar are all in motion, the single reconnaissance time length is T, and the motion speed vectors of the master station and the secondary station are v respectively_A,v_B,v_CThe motion velocity vector of the target radiation source is v_DThe distances from the radiation source to the main station and the auxiliary station are DA, DB and DC respectively. The radar pulse propagates at the speed of light c to reach the primary and secondary stations, t_A,t_B,t_CThe time length required by the same pulse signal to reach the primary station and the secondary station is respectively as shown in the formula (1):

t_A＝DA/c,t_B＝DB/c,t_C＝DC/c (1)

as can be seen from the position of the three-station electronic reconnaissance system in FIG. 4, a triangle is formed between the positions of the primary station and the secondary station and the position of the target radiation source, and the trilateral relationship of the triangle can be obtained

Combining the above formula (1) and formula (2), it is possible to obtain:

where σ denotes the time required for the electromagnetic wave to propagate from the secondary station B to the primary station a, and ρ denotes the time required for the electromagnetic wave to propagate from the secondary station C to the primary station a. Because the main station and the secondary station may move, the position information may have a certain error, the size of the error depends on the size of the relative position change of the main station and the secondary station in the duration of the signal detection and reception, and the upper limit of the size of the error can be obtained by equation (4):

wherein v is₁＝||v_A-v_BI or v₂＝||v_A-v_CAnd | μ represents an error. Thus, equation (3) can be adjusted to:

the time difference range between the main station and the two secondary stations can be respectively determined by the formula (5), and when the synchronous signal is not adopted, the time difference range can be obtained through the position information when the detection of the main station and the secondary station is finished and the movement speed of the two stations, so that the time range of pulse pairing between the main station and the secondary station can be reduced, and the sorting instantaneity is improved.

For this embodiment, the nth station A receives the signal in the stream of pulses_AThe arrival time of each pulse is recorded asPossible in the pulse stream of the signal received by the secondary station BThe time range of arrival of the same pulse signal from the same radiation source for receiving the pulse by the master station A is as follows:

Σ(n_A) The determined time range is that the master station A isThe time difference window for the received pulse at a time instant is shown in fig. 5.

2. And carrying out coarse pairing on the single pulse signal received by the main station in real time and the pulse signal received by each secondary station in the time difference window by utilizing the determined time difference window.

For the present embodiment, the time-difference window Σ (n)_A) The pulse of the inner secondary station B is detected, and the detected pulses are all possible to be n_AThe same pulse from the same radiation source, thereby realizing that the main station A is atThe time of receipt of a single pulse is roughly matched to the time of receipt of the pulse stream by secondary station B.

Suppose that in the equation of time window ∑ (n)_A) The inner secondary stations B are respectively atIf (L +1) pulses exist at the moment, the (L +1) pulse pairs can be obtained through rough matchingThe constructed sequence is simply called a coarse pulse pair sequence.

For the same reason time difference window sigma (n)_A) The pulse of the inner secondary station C is detected, and the main station A is realizedCoarse matching between the received monopulse and the stream of pulses received by the secondary station C to obtain (K +1) pulse pairs "Called as a coarse pulse pair sequence. FIG. 5 shows the Master station A atTime difference window Σ (n) corresponding to a single pulse received at a time_A) And using the moveout window Σ (n)_A) Respectively realize the master stationA is atThe single pulse received at the moment is roughly matched with the pulse streams of the secondary station B and the secondary station C, and corresponding rough pulse pairs can be obtained.

3. And (4) screening the coarse cluster pulse pairs obtained in the step (2) to obtain fine pulse pairs.

The secondary pulses obtained by coarse matching with the time difference window may be interference pulses, pulses emitted by other radiation sources and the same pulses from the same radiation source as the single pulse received by the primary station a at the time instant, and therefore, it is necessary to further precisely match the coarse pulse pair. The following describes the exact matching process in detail by taking the primary station a and the secondary station B as an example.

The main station A is atThe pulse received at a time is in the same position as that of the secondary station BThe multi-parameter distance of the pulse received at the moment except the TOA is as follows:

wherein, | | · | | represents solving 2-norm, W is a weighting matrix, generally a symmetric matrix. If in the pulse description wordAndand if the parameters are independent, W is a diagonal matrix, and the value of the diagonal element is inversely proportional to the variance of the measurement error. Due to the fact thatAndincluding parameters in the pulse such as carrier frequency, pulse amplitude, pulse width and the like, the size of the W diagonal element is the reciprocal of the measurement error variance under the corresponding dimension of the relevant parameter.

To indicate the degree of correlation between primary A and secondary B bursts, a multi-parameter similarity matching factor is introduced

Wherein r is₀Is the reference distance, representing the minimum multi-parameter distance at which the primary a and secondary B pulses are uncorrelated. Setting the multi-parameter matching threshold as epsilon, epsilon is epsilon to [0,1 ∈]For the purpose of determining whether two pulses are from the same pulse emitted by the same radiation source, since the difference between the same parameters is small when the same pulse from the same radiation source is received by two stations, epsilon is generally set to be greater than 0.5, for example, epsilon is 0.7. If it isThe parameters are successfully matched, and the 'pulse pair' successfully matched is recorded asOtherwise, the parameter matching is unsuccessful. The coarse pulse pairs can be further screened through the process, and the screened pulse pairs are called fine pulse pairs.

Similarly, for the master station AThe pulse received at a time is in the same direction as that of the secondary station CReceived at a time in the time difference window Σ (n)_A) The pulses in the interior are accurately matched, if the multiple parameters are successfully matched, the 'pulse pair' successfully matched is recorded asThus obtaining the fine pulse pair of the A station and the C station.

4. And calculating the time difference corresponding to each of the two fine pulse pairs, taking the two time differences as a time difference pair, and realizing the real-time sorting of the single pulse signals by using the time difference pair.

For thePulse n received by the master station a at a moment_AObtaining a fine pulse pair through rough pairing and precise pairingAndthe time differences of two pulse pairs are respectively

According to the formula (9), two fine pulse pairs can be obtainedAndcorresponding to "time difference pairs" (DTOA)_AB,DTOA_Ac). In the process of sorting the single pulse signals in real time, the pulse sequences corresponding to the sorted different radiation sources are respectively recorded as psi₁,Ψ₂,Ψ₃,…,Ψ_mThe corresponding time difference pair is (DTOA)_1AB,DTOA_1AC),(DTOA_2AB,DTOA_2AC),(DTOA_3AB,DTOA_3AC),……,(DTOA_mAB,DTOA_mAC). Wherein (DTOA)_iAB,DTOA_iAC) And i is 1,2, …, and m is the arithmetic mean value of each pulse pair corresponding to the time difference pair of each selected radiation source. Time difference DTOA of pulse arrival at primary station A and secondary station B_ABTime difference DTOA of the arrival of the pulse at the primary station A and the secondary station C_ACA two-dimensional rectangular plane coordinate system is established for the longitudinal axis, and the distribution of pulse sequences corresponding to the different selected radiation sources can be represented as shown in fig. 6.

Respectively calculating the distance between the time difference pair and the time difference pair corresponding to each sorted pulse sequence, selecting the minimum distance, and if the minimum distance is smaller than the error tolerance, classifying the pulse pair corresponding to the time difference pair into the sorted pulse sequence corresponding to the minimum distance; and if the minimum distance is not less than the error tolerance, a pulse sequence is newly built, and the pulse pair corresponding to the time difference pair is classified into the newly built pulse sequence. The distance calculation formula adopted is as follows:

let Δ ═ d_iminIf delta is less than or equal to delta and delta is the error tolerance of time difference pair of the same radiation source, the pulse pair is divided intoStoring the sequence Ψ_iAnd updated by averaging the time difference pairs corresponding to all the fine pulse pairs in the new sequence (DTOA)_iAB,DTOA_iAC) (ii) a Otherwise, newly creating a pulse sequence psi_m+1Will pulse pairStoring into Ψ_m+1. Thereby realizing the master stationAnd carrying out real-time pairing sorting on the single pulse received at the moment.

By the time difference DTOA of the arrival of the pulses at the primary A and secondary B stations_ABTime difference DTOA of the arrival of the pulse at the primary station A and the secondary station C for the X axis_ACAnd establishing a three-dimensional coordinate system for the Y axis and the Z axis as the time of the pulse reaching the main station A, wherein the time difference distribution of the pulse pairs subjected to pairing sorting in the coordinate system is shown in figure 7.

Through the process, the method can realize the sorting of all radar radiation sources, the embodiment sorts 8 radiation sources, and the pulse signal parameters of the 8 radiation sources are shown in table 1.

TABLE 1

The sorting results are shown in fig. 8-a, fig. 8-b, fig. 8-c, fig. 8-d, fig. 8-e, fig. 8-f, fig. 8-g and fig. 8-h, and it can be seen that the present invention can not only sort high repetition frequency and conventional repetition frequency pulse signals, but also sort ultra-low repetition frequency pulses, and has no false radiation source and very high sorting accuracy.