CN117554903A

CN117554903A - Sequence splicing and sorting method, system and storage medium based on pulse amplitude continuity

Info

Publication number: CN117554903A
Application number: CN202311520782.7A
Authority: CN
Inventors: 姬利海; 尹俊平; 沈卫超
Original assignee: Shanghai Zhangjiang Institute Of Mathematics
Current assignee: Shanghai Zhangjiang Institute Of Mathematics
Priority date: 2023-11-14
Filing date: 2023-11-14
Publication date: 2024-02-13

Abstract

The invention provides a sequence splicing and sorting method, a system and a storage medium based on pulse width continuity, wherein the method comprises the following steps: acquiring a plurality of initial radar patterns, and carrying out similarity comparison aiming at carrier frequency typical values, pulse width typical values, repetition interval typical values and arrival azimuth angle typical values of pulse sequences of different initial radar patterns; and combining the amplitude extension characteristic and the arrival time extension characteristic of the pulse sequence, combining partial initial radar patterns according to a preset rule for the initial radar patterns with fixed repetition intervals and pulse group types, and reducing the batch rate for subsequent sorting. The invention also discloses a sequence splicing and sorting system based on pulse amplitude continuity and a computer readable storage medium. The invention can reduce the sorting and batch increasing rate, can obtain more accurate carrier frequency, pulse width, repetition interval and other pattern parameters, has higher pulse loss resistance and pulse jitter resistance, and has better engineering application prospect.

Description

Sequence splicing and sorting method, system and storage medium based on pulse amplitude continuity

Technical Field

The invention relates to the technical field of counterreconnaissance, in particular to a sequence splicing and sorting method, a system and a storage medium based on pulse width continuity.

Background

The radar signal sorting is to separate out the pulses belonging to different radar radiation sources from the intercepted dense radar pulse streams (namely full pulses, each pulse is expressed by pulse description words), which is the basis of radar electronic countermeasure information analysis, and the characteristic parameters of the radar can be accurately measured and analyzed only after the random overlapped pulse streams are sorted into individual pulse sequences of each radar, so that the functional purpose, platform type, threat level and other attributes of the radars are judged, and the accurate interference is carried out on the enemy threat radiation sources.

Most of the radar signal sorting methods in the prior art are developed around five large parameters (arrival time TOA, carrier frequency RF, pulse width PW, pulse amplitude PA and arrival azimuth DOA). At present, parameters commonly used in a radar signal sorting algorithm are mainly carrier frequency, pulse width, pulse repetition interval and azimuth, and the method mainly comprises the steps of firstly carrying out pulse clustering by using the carrier frequency, the pulse width and the azimuth, and then determining a sorting pattern and a corresponding pulse sequence by using the repetition interval as a main sorting. The reason why the pulse width parameter is not easily applied is that: 1) The pulse width curve of each pattern is morphologically irregular (in comparison, RF, repetition intervals PRI, PW have some known pattern type); 2) The pulse amplitude value is closely related to the limit conditions such as the measurement range, the precision and the like of the reconnaissance equipment, the pulse amplitude values of different radar signal patterns are not greatly different, and the characteristics of a large number of flat tops and more wild points are often present. In a complex electromagnetic environment, pulse amplitude data are mixed together, disordered and difficult to be used as the basis of signal sorting. However, while the pulse width data is difficult to use directly in a clustering algorithm or a main sorting algorithm, it has a significant advantage: for a section of PA data with high enough confidence, it can be effectively determined whether the initially sorted intersecting pulse sequence belongs to the same pattern, because some PAs have continuous and clear envelopes, pulses obviously belonging to the same envelope can be determined to belong to the same pattern.

The known PRI sorting algorithm comprises sequence searching, traditional statistical histogram, CDIF, SDIF, PRI transformation and the like, and under the conditions of complex signal environments such as pulse loss, pulse jitter and the like, a pulse sequence belonging to the same pattern is often split into pulse sequences corresponding to a plurality of patterns, and the detection performance is required to be improved.

Therefore, the invention considers that the pulse amplitude data is used for further verification of the initial sorting pattern, and when the pulse amplitude continuity condition is met, the initial pattern can be combined into a more accurate new pattern, so that a more accurate RF, PRI, PW pattern parameter is obtained.

Disclosure of Invention

Aiming at the problems, the invention aims to provide a sequence splicing and sorting method, a system and a storage medium based on pulse amplitude continuity, which are used for carrying out similarity comparison on carrier frequencies, pulse widths, repetition intervals of different initial radar patterns and arrival azimuth angles of pulse sequences, and simultaneously considering amplitude extension characteristics of the pulse sequences, particularly PRI fixed and pulse group type patterns, and combining some initial patterns into more accurate new patterns according to relevant rules, so that sorting and batch increasing rate is reduced, and more accurate pattern parameters such as carrier frequencies, pulse widths, repetition intervals and the like are obtained.

In order to solve the technical problems, the invention provides the following technical scheme:

in one aspect, a sequence stitching and sorting method based on pulse width continuity is provided, which comprises the following steps:

acquiring a plurality of initial radar patterns, and carrying out similarity comparison aiming at carrier frequency RF typical values, pulse width PW typical values, repetition interval PRI typical values and arrival azimuth DOA typical values of pulse sequences of different initial radar patterns;

combining the amplitude PA extension characteristic and the arrival time TOA extension characteristic of the pulse sequence, combining partial initial radar patterns according to a preset rule for the PRI fixed and pulse group type initial radar patterns, and reducing the batch rate for subsequent sorting.

Preferably, the acquiring a plurality of initial radar patterns specifically includes:

s1, acquiring vectors formed by a plurality of initial radar patterns, and marking the vectors as initial radar pattern vectors A; each initial radar pattern includes: pulse number sequence, repetition interval PRI type and typical value, carrier frequency RF type and typical value, pulse width PW type and typical value, pulse arrival azimuth DOA typical value;

acquiring a vector formed by pulse descriptors to be sorted, and marking the vector as a pulse descriptor vector; each pulse descriptor includes: pulse number, pulse arrival time TOA, carrier frequency RF, pulse width PW, pulse amplitude PA, pulse arrival azimuth angle DOA.

Preferably, the performing similarity comparison for the carrier frequency RF typical values, the pulse width PW typical values, the repetition interval PRI typical values, and the arrival azimuth angle DOA typical values of the pulse sequence for different initial radar patterns specifically includes:

s2, the length of the initial radar pattern vector A is n, the initialized vector mark is a zero vector with the length of n, i is initialized to be 0, and the integer mark2 is 0;

step S3, when mark [ i ] is non-zero, making i=i+1, turning to step S3 for cyclic execution, otherwise, initializing jList as empty vector, self-increasing mark2 by 1, initializing j=i+1;

step S4, when mark [ j ] is non-zero, let j=j+1, go to step S4 to be executed circularly, otherwise, extract the repeat interval parameter types of A [ i ], A [ j ], if both are fixed type or pulse group type, execute step S5, otherwise, let j=j+1, go to step S4 to be executed circularly;

step S5, carrier frequency similarity comparison method similar_rf_PA and pulse width similarity comparison method similar_pw are applied to compare carrier frequency similarity and pulse width similarity of A [ i ] and A [ j ], when the carrier frequency similarity and the pulse width similarity are all satisfied, step S6 is executed, otherwise, j=j+1 is caused, and step S4 is executed;

s6, extracting the pulse number sequence corresponding to the known patterns A [ i ], A [ j ] respectively, and recording as t1, t2; taking the pulse number sequence as an index, obtaining a patterned azimuth vector from the pulse description word vector, comparing azimuth similarity, and executing step S7 if azimuth measurement values do not conflict;

and S7, analyzing the typical values of the repetition intervals of A [ i ], A [ j ], finding out the maximum value mxpri and the minimum value mnpri, when the difference between mxpri and mnpri meets the limiting condition, making j=j+1, turning to the step S4, otherwise, executing the subsequent steps.

Preferably, combining the amplitude PA extension characteristic and the arrival time TOA extension characteristic of the pulse sequence, and combining part of the initial radar patterns according to a predetermined rule for the initial radar patterns of the PRI fixed and pulse group types, wherein the method specifically comprises the following steps:

step S8, judging whether pulses with the same pulse number exist in the sequences t1 and t2, if so, enabling the variable SameID to be true, otherwise, enabling the variable SameID to be false;

step S9, combining t1 and t2, arranging according to ascending order of all pulse numbers to form a new sequence, recording by using st, analyzing the positions of pulses belonging to t1 and t2 in st to obtain a segmented splicing relationship of t1 and t2 in st, and recording pulse numbers of all splicing positions to form a temp vector; initializing the whole to 0 when the length of temp meets the preset condition and is smaller than the preset threshold value, wherein the whole is used for recording the times that the adjacent pulses of the sequences t1 and t2 meet the pulse arrival time and the pulse amplitude at the splicing position are consistent; initializing k to 0, and circularly executing step S9 until k is increased to the length of temp minus 1; when the length of temp does not meet the above condition, let j=j+1, go to step S4;

step S10, if SameID is false and all is 0, let j=j+1, go to step S4, otherwise, execute step S11;

step S11, jList vector increment element j, the value of mark2 is given to mark [ j ], let j=j+1, go to step S4;

step S12, circularly executing the steps S4 to S11 until j equals n-1 and jumps out of the cycle of the steps S4 to S11, combining the patterns with the A [ i ] pattern according to the pattern sequence number to be combined recorded by jList, updating the pulse number sequence of the A [ i ], updating the carrier frequency RF, the repetition interval PRI, the parameter type and the typical value of the pulse width PW, and updating the pulse to reach the azimuth DOA typical value;

step S13, circularly executing the steps S3 to S12 until i is equal to n-1, and deleting a pattern with mark marks being not 0 from the initial radar pattern vector A after the circle from the step S3 to the step S12 is skipped.

Preferably, the step S6 specifically includes the following steps:

step S61, considering the limitation of the precision of the angle measuring equipment, setting the azimuth values of a plurality of pulses in the pulse descriptor vector as invalid values, and setting the azimuth values as the same invalid azimuth value;

step S62, judging DOA typical value vector of the pattern: setting azimuth tolerance, wherein azimuth values within a tolerance allowable range are considered to be the same azimuth value; respectively analyzing azimuth vectors corresponding to the sequences t1 and t2, wherein invalid values are also included, and selecting azimuth typical value vectors with most 1-5 azimuth values to form respective patterns;

step S63, comparing the azimuth similarity: if the azimuth representative value vectors of A [ i ], A [ j ] each contain an effective value of the angle, but the effective values do not have the same value within the allowable range of azimuth tolerance, j=j+1 is caused to go to step S4, otherwise, the azimuth measurement values of A [ i ], A [ j ] are considered to be not in conflict in other cases, and step S7 is executed.

Preferably, in the step S7, when the difference between mxpri and mnpri satisfies the constraint condition means that: the difference between mxpri and mnpri is not allowed to be too large, in particular, if mnpri >100 μsec, mxpri/mnpri >2 is considered to be too large; if 10< mnpri < = 100 microseconds, the difference is considered too large when mxpri/mnpri > = 3.

Preferably, the step S9 specifically includes the following steps:

step S91, for temp [ k ], extracting 10 pulses which are positioned before temp [ k ] and contain temp [ k ] per se in st, and analyzing whether they belong to a single pattern in A [ i ] or A [ j ] at the same time; setting continuous 1 as true when meeting, otherwise setting false;

extracting 10 pulses located after temp k in st, analyzing whether they belong to a single pattern in A i or A j at the same time; setting continuous 2 as true when meeting, otherwise, setting false;

when temp, continue, continuous 2 meet the condition, go to step S92; otherwise, let k=k+1, go to step S91 to execute circularly;

preferably, in the step S91, the conditions to be satisfied are: if temp is 4 or less, it is required that both of the continue1 and the continue2 are true, if temp is greater than 4, only one of the continue1 and continue2 is true.

Step S92, setting a time continuity mark a as false, comparing an arrival time difference dDOA and a pulse amplitude difference dPA of pulses with the numbers of temp [ k ] and temp [ k ] +1, and if dDOA < = 0.02 seconds, dPA < = 0.6dB, setting a as true; let a be true if 0.02 seconds < dTOA < = 0.03 seconds, dPA < = 1.0 dB; when a is true, step S93 is executed, otherwise step S94 is executed.

Step S93, comparing the maximum time difference dTOA1 and the maximum pulse amplitude difference dPA1 of 10 pulses before temp [ k ] in the st sequence; comparing the maximum time difference dTOA2 and the maximum pulse amplitude difference dPA2 of 10 pulses after temp [ k ] in the st sequence; if the conditions are satisfied, the whole is self-increased by 1, and the step S10 is performed; if not, let k=k+1, go to step S91 to execute circularly;

preferably, in the step S93, the conditions to be satisfied are: meanwhile, the difference between dTOA1< = 0.04 seconds, dTOA2< = 0.04 seconds and dPA1 and dPA2 is less than or equal to 2dB.

Step S94, comparing the time difference between temp [ k ] and temp [ k-1] in the st sequence and the time difference between temp [ k+1] and temp [ k ], if the condition is satisfied, making k=k+1, and turning to step S91 for cyclic execution; otherwise, executing step S95;

preferably, in the step S94, the conditions to be satisfied are: at least one of the two time differences is greater than 0.4 seconds.

Step S95, calculating dTOA1, dTOA2, dPA1 and dPA, and simultaneously calculating intervals of 10 pulse amplitudes before and after temp [ k ] and 20 pulse amplitudes before and after temp [ k ]; if the conditions are met, enabling the alloys to be increased by 1, and turning to the step S10; otherwise, let k=k+1, go to step S91;

preferably, in the step S95, the conditions to be satisfied are: if dPA <6dB, dPA2<6dB, dPA1 and dPA have a difference of less than or equal to 2dB; or if dPA <6dB, dPA2<6dB, PA1 and PA2 are completely contained one by the other; or if dPA <6dB, dPA2<6dB, dPA1 and dPA2 are equal to or greater than 0.5dB.

In another aspect, a sequence stitching and sorting system based on pulse width continuity is provided, and the system comprises a processor module for executing the sequence stitching and sorting method based on pulse width continuity.

In another aspect, a computer readable storage medium having stored therein at least one instruction loaded and executed by a processor to implement the pulse width continuation based sequence stitching and sorting method is provided.

Compared with the prior art, the technical scheme provided by the invention has the following beneficial effects:

1. according to the invention, the PRI at the repetition interval and the pulse group type can be fixed, and some initial patterns can be combined according to the relevant rules, so that on one hand, the sorting and batch increasing rate is reduced, and on the other hand, more accurate pattern parameters are obtained.

2. The invention is based on the similarity comparison of carrier frequency RF typical values, pulse width PW typical values, repetition interval PRI typical values and azimuth DOA typical values of pulse sequences in different modes, and also considers the amplitude extension characteristics of the pulse sequences.

3. According to the method, pulse amplitude PA data are mainly considered to be used for further verification of the initial sorting patterns, when pulse amplitude continuity conditions are met, the initial patterns can be combined into more accurate new patterns, more accurate RF, PRI, PW pattern parameters are obtained, and radar targets possibly contained in the subsequent sorted patterns can be identified.

4. The invention applies the PA rule from the angle of optimizing and sorting, and compared with the existing sorting method utilizing the PA rule, the invention has higher pulse loss resistance and pulse jitter resistance, thereby having better engineering application prospect.

5. The invention does not need complex parameter selection, thereby greatly reducing the workload of operators and professional requirements.

Drawings

In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings required for the description of the embodiments will be briefly described below, and it is apparent that the drawings in the following description are only some embodiments of the present invention, and other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.

FIG. 1 is a diagram of radar signal pulse amplitude signatures analyzed in accordance with an embodiment of the present invention;

fig. 2 is a flowchart of a sequence stitching and sorting method based on pulse width continuity according to an embodiment of the present invention.

Detailed Description

In order to make the objects, technical solutions and advantages of the embodiments of the present invention more clear, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings of the embodiments of the present invention. It will be apparent that the described embodiments are some, but not all, embodiments of the invention. All other embodiments, which can be made by a person skilled in the art without creative efforts, based on the described embodiments of the present invention fall within the protection scope of the present invention.

The pulse width in the scout data is complex, and the following 2 strips can be used as the principle: 1) In the conventional mechanical scanning radar work, pulse amplitude in the reconnaissance data can show pulse amplitude characteristics caused by antenna scanning, and the pulse amplitude form is a quadratic function envelope shape formed by single points, as shown in (a) of fig. 1; 2) In one-dimensional phased array radar operation, the pulse amplitude can be embodied in a multi-point, flat-top, multi-layer, partial signal shape with quadratic function envelope as a whole, as shown in fig. 1 (b).

In the two-dimensional phased array radar, the pulse width morphology is not significantly regular, but pulses belonging to the same pattern also tend to be similar in pulse width morphology on the time axis, as shown in fig. 1 (c).

The PA of a single pattern has some characteristics, but each pattern in actual scout data has a complex background signal, and multiple conditions such as signal overlapping, high pulse density, pulse loss and the like commonly exist in the scout data. For example, the pulse density is too high, the signal is submerged by a dense signal set, and even if there is an envelope, the signal is almost covered all the time, as shown in (d) of fig. 1, and for example, the signal is limited to be lost due to shielding during reconnaissance, and a large number of pulses are lost, but the signal exists on intermittent time segments, and the envelope is incomplete, as shown in (e) of fig. 1.

The application of the CDIF, SDIF, PRI transformation equal sorting algorithm searches for pulse sequences and does not effectively sort a pattern through. For the conventional signal, due to the reasons of complex signal environment, missing pulse and the like, when the selected signal is blocked and hidden as shown in (d) and (e) in fig. 1, such a single sample search pulse is easy to split at certain time points, especially when the repetition interval extracted by two split pulse sequences and the carrier frequency typical values are not completely consistent, the two pulse sequences can be used as two patterns with different typical value parameters, so that sorting and batching are caused.

For phased array signals, the parameters of the patterns themselves are more variable, such as the cases of (b) and (c) in fig. 1, and the common cases of overlapping of a plurality of signals, high pulse density, pulse leakage and the like are also faced, as shown in (f) in fig. 1, so that the possibility of splitting the phased array patterns into different patterns is higher.

The invention is a method researched based on the problems faced by the practice, focuses on further optimizing the patterns which are sorted by other main sorting algorithms, combines and sorts certain signals according to the general rule of radar signals on the pulse amplitude, obtains more accurate sorting patterns and corresponding combined pulse sequences, and effectively reduces the sorting batch increasing rate.

The embodiment of the invention provides a sequence splicing and sorting method based on pulse amplitude continuity, which comprises the following steps of:

The method combines similarity comparison and TOA and PA continuity analysis to classify the signals of the initial radar patterns, so that the sorting and batch increasing rate can be reduced on one hand, and more accurate pattern parameters such as carrier frequency, pulse width, repetition interval and the like can be obtained on the other hand.

In an alternative embodiment, as shown in fig. 2, the acquiring a plurality of initial radar patterns specifically includes:

Further, the performing similarity comparison on the carrier frequency RF typical values, the pulse width PW typical values, the repetition interval PRI typical values, and the arrival azimuth angle DOA typical values of the pulse sequence for different initial radar patterns specifically includes:

step S2, the length of the initial radar pattern vector A is n, the initialized vector mark is a zero vector with the length of n, the initialized setting i is 0, and the integer mark2 is 0 (step S3 to step S12 are circularly executed until i is equal to n-1, and step S13 is shifted).

And S3, when mark [ i ] is non-zero, making i=i+1, turning to step S3 for cyclic execution, otherwise, initializing jList as a null vector, and self-increasing mark2 by 1, and initializing j=i+1 (cyclically executing steps S4 to S11 until j is equal to n-1, turning to step S12).

And S4, when mark [ j ] is non-zero, making j=j+1, turning to step S4 for cyclic execution, otherwise, extracting the repetition interval parameter types of A [ i ] and A [ j ], if both are of a 'fixed' or 'pulse group' type, executing step S5, otherwise, making j=j+1, turning to step S4 for cyclic execution.

And S5, comparing the carrier frequency similarity and the pulse width similarity of A [ i ], A [ j ] by using a carrier frequency similarity comparison method similar_rf_PA and a pulse width similarity comparison method similar_pw, and executing the step S6 when the carrier frequency similarity and the pulse width similarity are similar, otherwise, making j=j+1 and turning to the step S4.

S6, extracting the pulse number sequence corresponding to the known patterns A [ i ], A [ j ] respectively, and recording as t1, t2; and (7) taking the pulse number sequence as an index, obtaining the azimuth vector of the pattern from the pulse description word vector, comparing the azimuth similarity, and executing the step S7 if the azimuth measurement values do not conflict.

Further, the step S6 specifically includes the following steps:

in step S61, in consideration of the limitation of the accuracy of the angle measuring device, the azimuth values of a plurality of (sometimes almost half) pulses existing in the pulse descriptor vector are set as invalid values, and these azimuth values are set as the same invalid azimuth value.

In this embodiment, an invalid azimuth value is indicated by 380 °.

Step S62, judging DOA typical value vector of the pattern: setting azimuth tolerance, wherein azimuth values within a tolerance allowable range are considered to be the same azimuth value; and respectively analyzing azimuth vectors corresponding to the sequences t1 and t2, wherein an invalid value of 380 degrees is also included, and selecting azimuth typical value vectors with most 1-5 azimuth values to form respective patterns.

The orientation tolerance in this embodiment is chosen to be 8.

Step S63, comparing the azimuth similarity: if the azimuth representative value vectors of A [ i ], A [ j ] each contain valid values of angles other than 380 degrees, but the valid values do not have the same values within the allowable range of azimuth tolerance, j=j+1 is caused to go to step S4, otherwise, the azimuth measurement values of A [ i ], A [ j ] are considered not to conflict, and step S7 is executed.

And S7, analyzing the typical values of the repetition intervals of A [ i ], A [ j ], finding out the maximum value mxpri and the minimum value mnpri, when the difference between mxpri and mnpri meets the limiting condition, turning to the step S4, otherwise, executing the subsequent step S8.

The limiting conditions of this embodiment are: the difference between mxpri and mnpri is not allowed to be too large, in particular, if mnpri >100 μsec, mxpri/mnpri >2 is considered to be too large; when 10< mnpri < = 100, the difference is considered too large when mxpri/mnpri > = 3.

Further, combining the amplitude PA extension characteristic and the arrival time TOA extension characteristic of the pulse sequence, and combining part of the initial radar patterns according to a predetermined rule for the initial radar patterns of the PRI fixed and pulse group types, wherein the method specifically comprises the following steps:

step S8, judging whether pulses with the same pulse number exist in the sequence of t1 and t2 (pri formed by the pulses nearby and pulse repetition interval change rules of two patterns are met simultaneously), if yes, setting a variable SameID to true, otherwise, setting false.

Step S9, combining t1 and t2, arranging according to ascending order of all pulse numbers to form a new sequence, recording by st, analyzing the positions of pulses belonging to t1 and t2 in st to obtain a segmented splicing (mutual interleaving) relation of t1 and t2 in st, and recording pulse numbers at each splicing position to form a temp vector; initializing the whole to 0 when the length of temp meets the preset condition and is smaller than the preset threshold value, wherein the whole is used for recording the times that the adjacent pulses of the sequences t1 and t2 meet the pulse arrival time and the pulse amplitude at the splicing position are consistent; initializing k to 0, and circularly executing step S91 until k is increased to the length of temp minus 1; when the length of temp does not satisfy the above condition, let j=j+1, go to step S4.

In this embodiment, the length of temp needs to satisfy the following conditions: less than half of the length of t1 and less than half of the length of t2, the threshold conditions to be satisfied are: 500-1000 a.k.a. 800 a..

Further, the step S9 specifically includes the following steps:

when temp, continue, continuous 2 meet the condition, go to step S92; otherwise, let k=k+1, go to step S91 for loop execution.

In the step, the conditions to be satisfied are: if temp is 4 or less, it is required that both of the continue1 and the continue2 are true, if temp is greater than 4, only one of the continue1 and continue2 is true.

Step S93, comparing the maximum time difference dTOA1 and the maximum pulse amplitude difference dPA1 of 10 pulses before temp [ k ] in the st sequence; comparing the maximum time difference dTOA2 and the maximum pulse amplitude difference dPA2 of 10 pulses after temp [ k ] in the st sequence; if the conditions are satisfied, the whole is self-increased by 1, and the step S10 is performed; if not, let k=k+1, go to step S91 for loop execution.

In the step, the conditions to be satisfied are: meanwhile, the difference between dTOA1< = 0.04 seconds, dTOA2< = 0.04 seconds and dPA1 and dPA2 is less than or equal to 2dB.

Step S94, comparing the time difference between temp [ k ] and temp [ k-1] in the st sequence and the time difference between temp [ k+1] and temp [ k ], if the condition is satisfied, making k=k+1, and turning to step S91 for cyclic execution; otherwise, step S95 is performed.

In the step, the conditions to be satisfied are: at least one of the two time differences is greater than 0.4 seconds.

Step S95, calculating dTOA1, dTOA2, dPA1 and dPA, and simultaneously calculating intervals of 10 pulse amplitudes before and after temp [ k ] and 20 pulse amplitudes before and after temp [ k ]; if the conditions are met, enabling the alloys to be increased by 1, and turning to the step S10; otherwise, let k=k+1, go to step 91.

In the step, the conditions to be satisfied are: if dPA <6dB, dPA2<6dB, dPA1 and dPA have a difference of less than or equal to 2dB, enabling the alloys to increase by 1; or if dPA <6dB, dPA2<6dB, PA1 and PA2 are 2dB contained completely by one another, make the aliles self-increment by 1; or if dPA <6dB, dPA2<6dB, dPA1 and dPA have a difference of 0.5dB or more, the whole is self-increased by 1.

Step S10, if SameID is false and all is 0, let j=j+1, go to step S4, otherwise, go to step S11.

Step S11, jList vector increment element j, assign mark2 value to mark [ j ], let j=j+1, go to step S4.

And step S12, circularly executing the steps S4 to S11 until j is equal to n-1, skipping the cycle of the steps S4 to S11, combining the patterns with the A [ i ] pattern according to the pattern sequence number to be combined recorded by jList, updating the pulse number sequence of the A [ i ], updating the carrier frequency RF, the repetition interval PRI, the parameter type and the typical value of the pulse width PW, and updating the pulse to reach the azimuth DOA typical value.

In this embodiment, the updated azimuth typical value still selects the 1 to 5 most-occurring azimuth values.

Further, the specific process of the carrier frequency similarity comparison method similar_rf_pa in this embodiment includes:

step one, setting rfa as a vector formed by the carrier frequency type and the typical value of the pattern a, and setting rfb as a vector formed by the carrier frequency type and the typical value of the pattern b;

step two, when rfa and rfb are both fixed, the difference value of the two is within the tolerance range, and the two are judged to be similar, the new typical value is the typical value of rfa, otherwise, the two are judged to be dissimilar;

step three, rfa is fixed, rfb is pulse grouping or carrier frequency conversion (change among a plurality of typical values), if the typical value of rfa is compared with a certain typical value of rfb within a tolerance range, the typical value of rfa is judged to be similar to the typical value of rfb, the new typical value vector is a typical value vector of rfb, otherwise, the typical value vector and the typical value vector are judged to be dissimilar;

step four, when the rfa and the rfb are pulse groups or carrier frequency conversion, if the interval where all typical values of the rfa are located and the interval where all typical values of the rfb are located do not intersect, judging that the two are dissimilar; otherwise, judging that the two are similar, extracting the union set of all typical values of the two, and de-duplicating in the tolerance range to be used as a new typical value vector;

and fifthly, if one or both of rfa and rfb are jittering or agile, returning a result of dissimilar two, and reducing merging errors of a subsequent algorithm.

Further, the specific process of the pulse width similarity_pw method in this embodiment is similar to the carrier frequency similarity determination, but when there are multiple typical values of pulse width, the condition of similarity is relaxed, and only the typical values of pulse width are close, and even if there is no intersection, the result of "both are similar" is returned.

The method combines similarity comparison and TOA and PA continuity analysis to classify the signals of the initial radar patterns, so that the sorting and batch increasing rate is reduced on one hand, and more accurate pattern parameters such as carrier frequency, pulse width, repetition interval and the like are obtained on the other hand.

Correspondingly, the embodiment of the invention also provides a sequence splicing and sorting system based on the pulse width continuity, which comprises a processor module, wherein the processor module is used for executing the sequence splicing and sorting method based on the pulse width continuity. The system of the present embodiment may be used to implement the technical solution of the method embodiment shown in fig. 1, and its implementation principle and technical effects are similar, and are not described here again.

Embodiments of the present invention also provide a computer readable storage medium, such as a memory comprising instructions executable by a processor in a terminal to perform the above-described pulse width continuation based sequence stitching and sorting method. For example, the computer readable storage medium may be ROM, random Access Memory (RAM), CD-ROM, magnetic tape, floppy disk, optical data storage device, etc.

It should be noted that, in this document, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or terminal that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or terminal. Without further limitation, an element defined by the phrase "comprising one … …" does not exclude the presence of other like elements in a process, method, article or terminal device comprising the element.

References in the specification to "one embodiment," "an example embodiment," "some embodiments," etc., indicate that the embodiment described may include a particular feature, structure, or characteristic, but every embodiment may not necessarily include the particular feature, structure, or characteristic. Further, when a particular feature, structure, or characteristic is described in connection with an embodiment, it is submitted that it is within the knowledge of one skilled in the relevant art to effect such feature, structure, or characteristic in connection with other embodiments whether or not explicitly described.

The invention is intended to cover any alternatives, modifications, equivalents, and variations that fall within the spirit and scope of the invention. In the following description of preferred embodiments of the invention, specific details are set forth in order to provide a thorough understanding of the invention, and the invention will be fully understood to those skilled in the art without such details. In other instances, well-known methods, procedures, flows, components, circuits, and the like have not been described in detail so as not to unnecessarily obscure aspects of the present invention.

Those of ordinary skill in the art will appreciate that all or a portion of the steps in implementing the methods of the embodiments described above may be implemented by a program that instructs associated hardware, and the program may be stored on a computer readable storage medium, such as: ROM/RAM, magnetic disks, optical disks, etc.

The foregoing description of the preferred embodiments of the invention is not intended to limit the invention to the precise form disclosed, and any such modifications, equivalents, and alternatives falling within the spirit and scope of the invention are intended to be included within the scope of the invention.

Claims

1. The sequence splicing and sorting method based on pulse amplitude continuity is characterized by comprising the following steps of:

2. The sequence stitching and sorting method based on pulse width continuity according to claim 1, wherein the acquiring a plurality of initial radar patterns specifically comprises:

3. The pulse width continuity-based sequence stitching and sorting method according to claim 2, wherein the performing similarity comparison between the carrier frequency RF typical values, the pulse width PW typical values, the repetition interval PRI typical values, and the arrival azimuth angle DOA typical values of the pulse sequence for different initial radar patterns specifically comprises:

s5, comparing carrier frequency similarity and pulse width similarity of A [ i ], A [ j ] by carrier frequency similarity comparison method similar_rf_PA and pulse width similarity comparison method similar_pw, when the similarity is satisfied, executing

Step S6, otherwise, let j=j+1, go to step S4;

4. The pulse amplitude continuity-based sequence stitching and sorting method according to claim 3, wherein combining the amplitude PA extension characteristic and the arrival time TOA extension characteristic of the pulse sequence, for initial radar patterns of the PRI fixed and pulse group type, combines part of the initial radar patterns according to a predetermined rule, specifically comprising:

5. The sequence stitching and sorting method based on pulse width continuity according to claim 3, wherein the step S6 specifically comprises the following steps:

6. The sequence stitching and sorting method based on pulse width continuity according to claim 3, wherein in the step S7, when the difference between mxpri and mnpri satisfies a constraint condition means that: the difference between mxpri and mnpri is not allowed to be too large, in particular, if mnpri >100 μsec, mxpri/mnpri >2 is considered to be too large; if 10< mnpri < = 100 microseconds, the difference is considered too large when mxpri/mnpri > = 3.

7. The sequence stitching and sorting method based on pulse width continuity according to claim 4, wherein the step S9 specifically comprises the steps of:

when temp, continue, continuous 2 meet the condition, go to step S92; otherwise, let k=k+1, go to step S91 to execute circularly; wherein, the conditions to be satisfied are: if temp is 4 or less, it is required that both of the continue1 and the continue2 are true, if temp length is greater than 4, only one of the continue1 and continue2 is true;

step S92, setting a time continuity mark a as false, comparing an arrival time difference dDOA and a pulse amplitude difference dPA of pulses with the numbers of temp [ k ] and temp [ k ] +1, and if dDOA < = 0.02 seconds, dPA < = 0.6dB, setting a as true; let a be true if 0.02 seconds < dTOA < = 0.03 seconds, dPA < = 1.0 dB; when a is true, executing step S93, otherwise executing step S94;

step S95, calculating dTOA1, dTOA2, dPA1 and dPA, and simultaneously calculating intervals of 10 pulse amplitudes before and after temp [ k ] and 20 pulse amplitudes before and after temp [ k ]; if the conditions are met, enabling the alloys to be increased by 1, and turning to the step S10; otherwise, let k=k+1, go to step S91.

8. The sequence stitching and sorting method based on pulse width continuity according to claim 7, wherein in the step S93, the condition to be satisfied is: meanwhile, the difference between dTOA1< = 0.04 seconds, dTOA2< = 0.04 seconds and dPA1 and dPA2 is less than or equal to 2dB;

in the step S94, the conditions to be satisfied are: at least one of the two time differences is greater than 0.4 seconds;

in the step S95, the conditions to be satisfied are: if dPA <6dB, dPA2<6dB, dPA1 and dPA have a difference of less than or equal to 2dB; or if dPA <6dB, dPA2<6dB, PA1 and PA2 are completely contained one by the other; or if dPA <6dB, dPA2<6dB, dPA1 and dPA2 are equal to or greater than 0.5dB.

9. A sequence stitching and sorting system based on pulse width continuity, characterized by comprising a processor module for executing the sequence stitching and sorting method based on pulse width continuity according to any one of claims 1 to 8.

10. A computer readable storage medium having stored therein at least one instruction that is loaded and executed by a processor to implement the pulse width continuation based sequence stitching sorting method of any one of claims 1-8.