Search

CN-116859342-B - Radar signal sorting method and device

CN116859342BCN 116859342 BCN116859342 BCN 116859342BCN-116859342-B

Abstract

The invention discloses a radar signal sorting method and device, the method comprises the steps of obtaining radar signal information, wherein the radar signal information comprises radar target number and arrival time sequence information of intercepted radar pulses, establishing a radar pulse repetition frequency parameter model by utilizing the radar signal information, carrying out radar label distribution processing on the arrival time sequence information of the intercepted radar pulses by utilizing the radar target number to obtain a radar pulse path set, and screening the radar pulse path set by utilizing a spurious filtering discrimination model and the radar pulse repetition frequency parameter model to obtain radar pulse sequence target information. Aiming at the phenomena of pulse missing and stray pulse, the invention inputs non-ideal interception pulse trains, restores pulse characteristics before the pulse missing and stray pulse occur by designing specific likelihood factors, and improves the accuracy rate of radar pulse signal sorting.

Inventors

  • FANG YI
  • YANG JIAN
  • CHEN KEYU
  • XIAO DEZHENG

Assignees

  • 中国人民解放军军事科学院系统工程研究院

Dates

Publication Date
20260505
Application Date
20230628

Claims (10)

  1. 1. A radar signal sorting method, comprising: S1, acquiring radar signal information, wherein the radar signal information comprises the number of radar targets and arrival time sequence information of intercepted radar pulses; s2, establishing a radar pulse repetition frequency parameter model by utilizing the radar signal information; s3, carrying out radar tag distribution processing on the arrival time sequence information of the intercepted radar pulse by utilizing the number of radar targets to obtain a radar pulse path set; And S4, screening the radar pulse path set by using the spurious filtering discrimination model and the radar pulse repetition frequency parameter model to obtain radar pulse sequence target information, wherein the radar pulse sequence target information is used for representing radar targets corresponding to the intercepted radar pulses.
  2. 2. The method of claim 1, wherein the number of radar targets is denoted by K, and the time series of arrival information of the intercepted radar pulse is denoted by , wherein, Representing the arrival time of the ith intercepted radar pulse, Representing the number of intercepted radar pulses.
  3. 3. The radar signal sorting method of claim 2, wherein the radar pulse repetition frequency parameter model is used for representing a probability density function of pulse repetition frequency of an intercepted radar pulse, and comprises a first PRI model, a second PRI model, a third PRI model and a fourth PRI model, wherein the pulse repetition frequency is represented by PRI.
  4. 4. A radar signal sorting method as claimed in claim 3, wherein the first PRI model, which intercepts the PRI of radar pulses, has a parameter space of , , For the mean value of the PRI of the intercepted radar pulse in the first PRI model, For the variance of the PRI of the intercepted radar pulse in the first PRI model, The probability density function of PRI of the intercepted radar pulse of the first PRI model is the standard deviation of PRI of the intercepted radar pulse of the first PRI model The expression of (2) is: Wherein, the In order to accumulate the normal distribution function, PRI of the radar pulse is intercepted at the moment t; the second PRI model has the PRI parameter space of intercepting radar pulse of , , For the mean value of the PRI of the intercepted radar pulse in the second PRI model, For the variance of the PRI of the intercepted radar pulse in the second PRI model, The probability density function of PRI of the intercepted radar pulse of the second PRI model is the standard deviation of PRI of the intercepted radar pulse of the second PRI model The expression of (2) is: , Wherein, the The PRI of the radar pulse is intercepted for time t-1, A sliding step length of a second PRI model; the third PRI model has a parameter space for intercepting a radar pulse sequence of , , wherein, , Is a PRI state transition matrix that is a matrix of PRI states, Representing the probability that PRI transitions from i state to j state, M is the number of states of PRI, , Is the initial state distribution of the PRI of the intercepted radar pulse sequence in the third PRI model, PRI representing the sequence of intercepted radar pulses in a third PRI model The number of initial states is set to be, , , A set of gaussian distribution models representing the PRI states, A gaussian distribution model representing the mth state of PRI, , The mean of the gaussian distribution representing the mth state of PRI, Representing the variance of the gaussian distribution of the mth state of PRI, PRI sequence of interception radar pulse of the third PRI model takes value It is indicated that the number of the elements is, , In (a) and (b) Representing the value of the 1 st PRI of the PRI sequence, and so on, In (a) and (b) A value of a kth PRI representing the PRI sequence, k representing the number of PRIs contained in the PRI sequence, a probability distribution of the PRI sequence The calculated expression of (2) is: In the formula, Is that A probability density function of the corresponding gaussian distribution, Probability values representing the PRI sequence when the 1 st PRI is in the mth state, and so on, the A probability value representing when the (t+1) th PRI of the PRI sequence is in the (m) th state, the (d) A probability value representing when a t-th PRI of the PRI sequence is in an h-th state, Representing the matrix The h th row, m th column element, Represents the value of the t+1th PRI of the PRI sequence, ; The expression of the PRI of the fourth PRI model for intercepting radar pulse is as follows: , Wherein, the For the fourth PRI model amplitude, For the fourth PRI model frequency, For the fourth PRI model amplitude offset, For the pulse sequence number of the intercepted radar pulse, Gaussian noise for the fourth PRI model amplitude, , For the initial phase of the fourth PRI model amplitude, The PRI parameter space of the intercepted radar pulse of the fourth PRI model is , Probability Density function of PRI of intercepted radar pulse of fourth PRI model The expression of (2) is: Wherein, the The variance of gaussian noise for the fourth PRI model amplitude.
  5. 5. The radar signal sorting method according to claim 4, wherein the radar tag allocation processing is performed on the arrival time sequence information of the intercepted radar pulse by using the number of radar targets to obtain a radar pulse path set, including: s31, sequencing the arrival time sequence information of the intercepted radar pulse according to time sequence to obtain a pulse time sequence; s32, sequentially allocating a radar tag set to the first arrival time information in the pulse time sequence Each radar tag in ], obtaining K+1 initialized radar pulse paths, said radar pulse paths being a sequence of several radar tags, said radar tag set [ The first K radar tags in the section are respectively used for representing the 1 st radar to the K radar to which the intercepted radar pulse belongs, and spurious is used for representing that the intercepted radar pulse is a stray pulse; S33, constructing a radar pulse path set by using the K+1 initialized radar pulse paths, wherein the upper limit of the radar pulse path set is A radar pulse path, wherein e is a preset index value, Is a positive integer of 4 or more.
  6. 6. The radar signal sorting method according to claim 5, wherein the filtering the radar pulse path set to obtain radar pulse sequence target information by using the spurious filtering discrimination model and the radar pulse repetition frequency parameter model includes: S401, initializing an iteration sequence number t, to obtain t=2; S402, sequentially allocating a radar tag set to the t-th arrival time information in the pulse time sequence Each radar tag in the list, obtaining a path tail node set, wherein the path tail node set comprises K+1 radar tags; S403, for each radar pulse path, adding each radar tag in the path tail node set to the last of the radar pulse paths, and carrying out path merging processing on the radar pulse paths and the radar tags in the path tail node set to obtain K+1 expanded radar pulse paths; s404, judging whether the iteration sequence number t is larger than e-1, and obtaining a storage judging result; if the storage discrimination result is no, executing step S405; if the storage discrimination result is yes, executing step S409; s405, storing the extended radar pulse paths obtained in the S403 into a radar pulse path set; s406, calculating to obtain a similarity value of each radar pulse path of the radar pulse path set by utilizing the radar target similarity function; s407, storing the calculated similarity value of each radar pulse path into a similarity value set, wherein the similarity value set comprises the similarity value and the radar pulse path information corresponding to the similarity value; S408, setting an iteration sequence number t to increase by 1, and executing step S402; s409, calculating the similarity value of the extended radar pulse path obtained in the S403 by utilizing a radar target similarity function; S410, sorting the similarity values of the similarity value set and the similarity values of the extended radar pulse paths according to the value sequence from big to small, and searching to obtain the value before sorting Similarity value of names and corresponding radar pulse path information; s411, performing a cleaning process on the similarity value set, and before sorting the values The similarity value of the name is stored into a similarity value set, the radar pulse path set is cleared, and the values are sequenced before The radar pulse paths corresponding to the similarity values of the names are stored in the radar pulse path set; s412, judging whether the iteration sequence number t is larger than the number N of the intercepted radar pulses, and obtaining a stop iteration judgment result; If the iteration stopping judgment result is yes, selecting and obtaining a maximum similarity value from the similarity value set, and determining a radar pulse path corresponding to the maximum similarity value; If the iteration stopping judgment result is no, setting the iteration sequence number t to increase by 1, and executing step S402.
  7. 7. The radar signal sorting method of claim 6, wherein the radar target similarity function has the expression: In the formula, In order to calculate the value of the similarity, In order to be a radar pulse path, , For the number of radar tags contained in the radar pulse path, , For the ith radar tag of the radar pulse path, A spurious pulse indicator for an ith radar tag; for radar pulse paths The radar tag k appearing in (a) is in the radar pulse path At the ii-th occurrence in the radar pulse path Is provided with a position number in the middle, Is that The number of the elements contained in the composition, The method comprises the steps of filtering coefficients for preset stray pulses; , Is the first The time interval during which the radar pulse is intercepted, Is the first The arrival time of the intercepted radar pulse.
  8. 8. A data processing system for radar signal sorting, comprising: A first memory, a first processor and a computer program stored on the first memory and executable on the first processor, the first processor implementing the radar signal sorting method according to any one of claims 1-7 when the program is executed.
  9. 9. A data processing apparatus for radar signal sorting, the apparatus comprising: A memory storing executable program code; A processor coupled to the memory; The processor invokes the executable program code stored in the memory to perform the radar signal sorting method of any one of claims 1-7.
  10. 10. A computer-storable medium storing computer instructions that, when invoked, are adapted to perform the radar signal sorting method according to any one of claims 1-7.

Description

Radar signal sorting method and device Technical Field The invention relates to the field of radar signal sorting, in particular to a radar signal sorting method and device. Background In a non-cooperative electromagnetic environment, interleaved pulse streams intercepted by a scout device come from different radar radiation sources. The correct sorting of the pulse sequence of the corresponding radar from the obtained pulse stream is an important precondition for correct analysis and identification of the radar signal and ensuring the subsequent effective interference. In recent years, with the development of radar technology, the electromagnetic environment of a battlefield is unprecedented complex, and the task of sorting radar pulse signals faces a great challenge. On one hand, with the development of a radar technology with low interception probability and the parameter overlapping of pulse signals on a time-space frequency domain, the reconnaissance equipment is difficult to intercept all pulse signals, the pulse missing phenomenon is serious, and on the other hand, under the influence of various marine ship radiation electromagnetic signals, the reconnaissance equipment is easy to intercept a large number of stray pulses irrelevant to targets in a high-density signal environment, so that the sorting difficulty is increased. In addition, the radar working mode is updated and iterated rapidly, a complex PRI modulation mode is continuously emerging, and a traditional algorithm based on sequence search is difficult to adapt to a complex interleaving pulse sorting task, so that a new sorting method needs to be researched. In particular, for non-ideal situations of pulse missing and spurious pulses, the existing method is difficult to sort radar signals correctly. Disclosure of Invention Aiming at the problems that the radar working mode is updated and iterated rapidly and the complex PRI modulation mode is continuously emerging, so that the radar signals are difficult to accurately sort by the existing radar sorting method, the invention discloses a radar signal sorting method, which comprises the following steps: S1, acquiring radar signal information, wherein the radar signal information comprises the number of radar targets and arrival time sequence information of intercepted radar pulses; s2, establishing a radar pulse repetition frequency parameter model by utilizing the radar signal information; s3, carrying out radar tag distribution processing on the arrival time sequence information of the intercepted radar pulse by utilizing the number of radar targets to obtain a radar pulse path set; And S4, screening the radar pulse path set by using the spurious filtering discrimination model and the radar pulse repetition frequency parameter model to obtain radar pulse sequence target information, wherein the radar pulse sequence target information is used for representing radar targets corresponding to the intercepted radar pulses. The number of radar targets is represented by K, and the arrival time sequence information of the intercepted radar pulse is represented by { tau i }, i=1,.. N, wherein tau i represents the arrival time of the ith intercepted radar pulse, and N represents the number of the intercepted radar pulses. The radar pulse repetition frequency parameter model is used for representing a probability density function of pulse repetition frequency of an interception radar pulse and comprises a first PRI model, a second PRI model, a third PRI model and a fourth PRI model, wherein the pulse repetition frequency is represented by PRI. The first PRI model, which intercepts the PRI of the radar pulse, has a parameter space Θ 1,Mu 1 is the mean value of the PRI of the intercepted radar pulse in the first PRI model,For the variance of the PRI of the intercepted radar pulse in the first PRI model, σ 1 is the standard deviation of the PRI of the intercepted radar pulse in the first PRI model, and the expression of the probability density function ζ 1(pt) of the PRI of the intercepted radar pulse of the first PRI model is: wherein phi () is a cumulative normal distribution function, and p t is PRI for intercepting radar pulse at t moment; The second PRI model, which intercepts the PRI of the radar pulse, has a parameter space Θ 2, Mu 2 is the mean value of the PRI intercepting the radar pulse in the second PRI model,For the variance of the PRI of the intercepted radar pulse in the second PRI model, σ 2 is the standard deviation of the PRI of the intercepted radar pulse in the second PRI model, and the expression of the probability density function ζ 2(pt) of the PRI of the intercepted radar pulse in the second PRI model is: Wherein p t-1 is PRI of intercepting radar pulse at t-1 moment, and alpha is sliding step length of the second PRI model; The third PRI model has a parameter space Θ 3,Θ3 = (a, B, pi) for intercepting a radar pulse sequence, wherein a= [ a ij]M×M, a is a PRI state transition matrix