Multi-radiation source combination pulse sequencing method
Technical Field
The invention belongs to the field of radar target signal generation and simulation, and mainly relates to a sequencing method for output pulse signals of a multi-radiation source combined circuit.
Background
When the function of the airborne radar warning equipment is tested on the ground, an intensive electromagnetic environment needs to be constructed firstly. With the rapid development of computer simulation technology in recent years, more researches tend to fuse multiple radar signals into a high-density pulse description word sequence by using software programming, and the sequence can be used for simulating the multi-source signals by only one radar signal simulator. However, in the above process, there is a high possibility that a plurality of radar signals may generate pulse overlapping in time, and too much overlapping may cause a part of pulses to be lost or parameters to be changed, so that the warning device cannot rapidly and effectively sort and identify electromagnetic signals. Therefore, how to effectively handle the overlapping phenomenon in the pulse sequencing is an urgent problem to be researched in the process of constructing the intensive electromagnetic environment by utilizing computer simulation.
Disclosure of Invention
In view of the above, on the basis of analyzing the existing multi-radar signal sorting method, the invention designs a new pulse sorting algorithm to reduce the pulse overlapping probability as much as possible, preserve the pulse integrity, make the constructed electromagnetic environment denser, make the equipment more quickly and accurately complete the sorting and identification of the radiation source, and finally improve the reliability of the radar warning equipment detection.
In order to achieve the above object, the present invention provides a method for sequencing combined pulses from multiple radiation sources, including:
step S1, according to the parameters of the multiple radiation source signals, the pulse arrival time of each path of signals is changed in turn to generate a combined signal sequence, the optimal output effect is sought, and the pulse delay data of each path is recorded;
and step S2, sequencing according to each path of pulse delay data corresponding to the optimal output effect, and generating a high-density pulse description word sequence.
2. The method as claimed in claim 1, wherein the signal parameters in step S1 include pulse width, pulse repetition interval, pulse amplitude, and carrier frequency of each radar radiation signal, and the parameter values of each radar radiation signal are constant during the sequencing process.
3. The method of claim 1, wherein the generating of the combined signal sequence of step S1 includes:
step S11, first, selecting one path as reference sequence according to the multi-path signals provided by the scene, wherein the sequence will not be delayed; then, the signal sequence with the largest PRI is selected from the multi-path signals, and the PRI is pmTruncating n from time 0mFor a pulse repetition period, the sequence length L can be expressed as:
L=pm*nm
the length of each sequence is L, and the number of pulses n contained in the signal sequence with the sequence number of k is nkComprises the following steps:
nk=INT(L/pk)
wherein p iskIndicates the PRI value of the signal sequence with the sequence number k and satisfies pk≤pm. INT denotes a rounding operation.
Step S12, constructing a multi-source signal sequence group, wherein the array dimension is 3 and is represented by dlist [ k, t, j ], and the meaning of each dimension is as follows:
k: the radar number is more than or equal to 0 and less than rnum, and rnum is the number of radars;
t: the signal delay amount is 1 μ s. Let reference sequence PRI be p0By using the periodicity of the pulse signal, the delay range of other sequences only needs to generate the delay length of a repetition period relative to the reference signal, and the optimal combining effect can be achieved, namely t is more than or equal to 0 and less than p0;
j: discrete time points in the signal sequence have the unit of 1 mu s, and j is more than or equal to 0 and less than L.
Step S13, implementing triple traversal using k, t, j as variables, to obtain a temporary combination signal sequence plist [ j ]:
in the above formula, the summation symbol represents the corresponding summation of array elements, so plist [ j ] length is also L. Because of the pulse overlapping phenomenon, plist [ j ] does not represent the final combined output result and is only used as an important criterion for pulse optimization selection and deletion.
And step S14, taking the radar of the current emission signal as a main line, and performing pulse optimization selection and rejection to obtain an effective combination sequence alist.
In step S15, a reference value maxt is set, with an initial value of 0. And calculating the number of combined output pulses or the duty ratio n when the current delay amount of each path is tau [ k ], if n is more than maxt, the maxt is n, and recording the delay amount tau [ k ] of each path at the moment.
4. The method according to claim 3, wherein the pulse optimization in step S14 is selected from the group consisting of:
(1) the radar serial number lsa of the transmitted signal in the initial state is 0;
(2) setting i as a circulation variable, satisfying (for i ═ 0; i < L; i + +;
(3) when plist [ i ] is equal to 0 while plist [ i +1 ]! When the pulse number is equal to 0, a new pulse transmitting condition is met, the radar serial number k at the moment is recorded, and lsa is equal to k;
(4) when (3) the new pulse emission condition is not satisfied, if there is plist [ i ]! 0 and plist [ i +1 ]! 0 and plist [ i +1 ]! The combination sequence value is not zero and changes as plist [ i ] and dlist [ lsa ] [ t ] [ i +1] ═ 0; and dlist [ lsa ] [ t ] [ i +1] ═ 0, that is, if the change occurs at the trailing edge of the current reserved pulse sequence, the pulse is reserved, and the radar serial number k corresponding to the pulse is recorded, and lsa ═ k;
5. the method of claim 1, wherein the generating of the combined signal sequence in step S1 is performed by: when a program is designed, the core time-consuming operation is changed into multi-thread parallel operation, each thread occupies an independent CPU core, the operation is not interfered with each other, and an operator can still operate a software main interface to perform other work in the operation process.
6. The method of claim 1, wherein the pulse description words in step S2 include pulse arrival time, pulse width, pulse repetition interval, pulse amplitude, and carrier frequency of each pulse; combining the message head and the message tail, the pulse description word can be used as a message for the FPGA to solve in real time; all parameters are coded in a 16-system mode, wherein the content of a message header and the content of a message tail are fixed, the rest data are generated according to the principle that a low data bit is in front of a high data bit and a high data bit is behind the low data bit, and based on the rules, the combined output pulse sequence can be converted into a PDW sequence file for hardware to read.
Drawings
Fig. 1 is a general flowchart of a method for sequencing combined pulses from multiple radiation sources according to the present invention.
Fig. 2 is a schematic diagram of a multi-radiation source combined pulse sequencing method according to the present invention, showing the selection scheme when pulses are overlapped.
Fig. 3 is a waveform diagram of a combined output after being sequenced by the method for sequencing combined pulses from multiple radiation sources according to the present invention.
Fig. 4 is a pulse description word sequence diagram generated by a multi-radiation source combined pulse sequencing method according to the invention.
Detailed Description
In order to make the implementation objects, technical solutions and advantages of the present application clearer, the technical solutions in the embodiments of the present application will be described in more detail below with reference to the drawings in the embodiments of the present application. In the drawings, the same or similar reference numerals denote the same or similar elements or elements having the same or similar functions throughout. The described embodiments are some, but not all embodiments of the present application. The embodiments described below with reference to the drawings are exemplary and intended to be used for explaining the present application, and should not be construed as limiting the present application. All other embodiments obtained by a person of ordinary skill in the art based on the embodiments of the present application without any creative effort belong to the protection scope of the present application. Embodiments of the present application will be described in detail below with reference to the drawings.
With reference to fig. 1, the technical solution of the present invention is to design a method for sequencing combined pulses of multiple radiation sources, which includes the following steps:
step S1, according to the parameters of the multiple radiation source signals, the pulse arrival time of each path of signals is delayed in turn to generate a combined signal, the optimal output effect is sought, and the pulse delay data of each path is recorded;
and step S2, sequencing according to each path of pulse delay data corresponding to the optimal output effect, and generating a high-density pulse description word sequence.
Further, the signal parameters described in step S1 should include pulse width, pulse repetition interval, pulse amplitude, and carrier frequency, and the present invention focuses on the research of pulse sequencing algorithm, so that the conventional pulse signal with a simpler form is selected.
Further, the generating of the combined signal in step S1 includes the following:
step S11, first, the rootAccording to the multipath signals provided by the scene, one path is selected as a reference sequence, and the sequence is not subjected to delay processing. Then, the signal sequence with the largest PRI is selected from the multi-path signals, and the PRI is pmTruncating n from time 0mFor a pulse repetition period, the sequence length L can be expressed as:
L=pm*nm
the length of each sequence is L, and the number of pulses n contained in the signal sequence with the sequence number of k is nkComprises the following steps:
nk=INT(L/pk)
wherein p iskIndicates the PRI value of the signal sequence with the sequence number k and satisfies pk≤pm. INT denotes a rounding operation.
Step S12, constructing a signal sequence group by a plurality of radar transmission signal parameters in scene setting, including information such as carrier frequency, pulse width, pulse repetition interval, amplitude and the like, wherein the dimension of the group is 3, and is represented by dlist [ k, t, j ], and the meaning of each dimension is as follows:
k: the radar number is more than or equal to 0 and less than rnum, and rnum is the number of radars;
t: the signal delay amount is 1 μ s. Let reference sequence PRI be p0By using the periodicity of the pulse signal, the delay range of other sequences only needs to generate the delay length of a repetition period relative to the reference signal, and the optimal combining effect can be achieved, namely t is more than or equal to 0 and less than p0。
j: discrete time points in the signal sequence have the unit of 1 mu s, and j is more than or equal to 0 and less than L.
The pulse delay array is now complete and contains all the delay possibilities of the signals relative to the reference signal.
Step S13, implementing triple traversal using k, t, j as variables, to obtain a temporary combination signal sequence plist [ j ]:
in the above formula, the summation symbol represents the corresponding summation of array elements, so plist [ j ] length is also L. Because of the pulse overlapping phenomenon, plist [ j ] does not represent the final combined output result and is only used as an important criterion for pulse optimization selection and deletion.
And step S14, taking the radar of the current emission signal as a main line, and performing pulse optimization selection and rejection to obtain an effective combination sequence alist.
Further, the pulse optimization trade-off in step S14 adopts the trade-off scheme shown in fig. 2, and the algorithm steps are as follows:
(1) the radar serial number lsa of the transmitted signal in the initial state is 0;
(2) setting i as a cyclic variable, satisfying (for i ═ 0; i < L; i + +).
(3) When plist [ i ] is equal to 0 while plist [ i +1 ]! When the pulse number is equal to 0, a new pulse transmitting condition is met, the radar serial number k at the moment is recorded, and lsa is equal to k;
(4) when (3) the new pulse emission condition is not satisfied, if there is plist [ i ]! 0 and plist [ i +1 ]! 0 and plist [ i +1 ]! The combination sequence value is not zero and changes as plist [ i ] and dlist [ lsa ] [ t ] [ i +1] ═ 0; and dlist [ lsa ] [ t ] [ i +1] ═ 0, that is, if the change occurs at the trailing edge of the current reserved pulse sequence, the pulse is reserved, and the radar serial number k corresponding to the pulse is recorded, and lsa ═ k;
in step S15, a reference value maxt is set, with an initial value of 0. And calculating the pulse number or duty ratio n of the alist output by combining when the delay amount of each current path is tau [ k ], if n is larger than maxt, the maxt is n, and recording the delay amount tau [ k ] of each path at the moment. Based on the above method, the value τ [ k ] when the total pulse number (duty ratio) in the combined pulse sequence is maximum is found, and the combined pulse sequence generated based on the delay value is the optimal combined output sequence, as shown in fig. 3 (duty ratio is maximum).
Further, the generating of the combining signal in step S1 has a large calculation amount due to the operations involving flows and factors such as multidimensional array traversal, loop nesting, superposition and rejection, complex structure, and the like; in addition, the hardware conditions such as the internal memory of the computer, the core number of the CPU and the like are limited, so that the problem of too long operation time due to too large operation amount occurs, and the efficiency verification of the alarm device cannot be completed quickly. In order to solve the problemMultithreading is to be used. When a program is designed, the core time-consuming operation is changed into multi-thread parallel operation, each thread occupies an independent CPU core, the operation is not interfered with each other, and an operator can still operate a software main interface to perform other work in the operation process. For example, a computer CPU is configured as an 8-core, the program follows a traversal method with dlist k, t, j in each thread]The traversal range of the medium variable t is only provided with the reference signal p 01/8 of (1); this greatly shortens the overall time consumption 7/8.
Further, the pulse descriptor in step S2 should include pulse arrival time, pulse width, pulse repetition interval, pulse amplitude, and carrier frequency of each pulse; combining the message head and the message tail, the pulse description word can be used as a message for the FPGA to solve in real time; all parameters are coded in a 16-system mode, wherein the content of a message header and the content of a message tail are fixed, the rest data are generated according to the principle that a low data bit is in front of a high data bit and a high data bit is behind the low data bit, and based on the rules, the combined output pulse sequence can be converted into a PDW sequence file, for example, as shown in the attached figure 4, and the PDW sequence file can be read by hardware.