CN113296071A - Multi-spread signal sorting method and device based on correlation matching method - Google Patents
Multi-spread signal sorting method and device based on correlation matching method Download PDFInfo
- Publication number
- CN113296071A CN113296071A CN202110563419.8A CN202110563419A CN113296071A CN 113296071 A CN113296071 A CN 113296071A CN 202110563419 A CN202110563419 A CN 202110563419A CN 113296071 A CN113296071 A CN 113296071A
- Authority
- CN
- China
- Prior art keywords
- pulse
- value
- sequence
- calculating
- period
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/418—Theoretical aspects
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention provides a multi-staggered signal sorting method and a device based on a correlation matching method, wherein the method comprises the following steps: s101, calculating the projection energy of a 01 sequence converted by a multi-staggered signal pulse stream in a unit orthogonal subspace; s102, calculating the corresponding measurement length of each period value according to the projection energy; s103, calculating a detection threshold corresponding to each period value; s104, selecting the period values with the corresponding measurement lengths larger than the detection threshold value to form a to-be-selected set; s105, obtaining a first PRI value according to the divisor number of the period values in the to-be-selected set and the size of each measurement length; s106, separating a first pulse sequence from the multi-staggered signal pulse stream according to the first PRI value; and S107, repeating S101-S106 until the length of the residual pulse sequence is smaller than the sorting threshold, and mixing the first pulse sequences with the same PRI value to obtain a staggered sequence.
Description
Technical Field
The present invention relates to the field of Electronic Support systems (ESM) signal sorting technologies, and in particular, to a multi-parameter signal sorting method and apparatus based on a Correlation Matching Method (CMM).
Background
Pulse Repetition Interval (PRI) is one of the basic characteristic parameters of a pulsed radar signal. The pulse signal transmitted using the spread PRI technique is referred to as a spread pulse signal. The staggered pulse signal is composed of a plurality of pulse sequences with the same PRI value and different initial phases, the same PRI value is called a frame period of the staggered signal, and a time interval between the initial phases is called a sub-pulse interval (sub-PRI) of the staggered signal. The frame period and the sub-pulse interval of the estimated stagger signal are called stagger signal de-interleaving.
The existing methods for interleaving of staggered signals have different technical defects, for example, model mismatching is easy to occur, so that the sorting precision is poor, or the method is sensitive to pulse missing, or only one staggered signal can be sorted, and when a pulse stream has a plurality of staggered signals, the signal sorting cannot be finished. In a complex electromagnetic environment, a radar signal intercepted by an ESM receiver is an approximately random radar pulse stream consisting of a plurality of radar pulse sequences, the approximately random radar pulse stream not only comprises a plurality of radar radiation sources, but also the signal of each radiation source may have pulse loss, so that the prior art cannot well solve the requirement of signal sorting in the actual complex electromagnetic environment.
Disclosure of Invention
Technical problem to be solved
Aiming at the technical problems in the prior art, the invention provides a multi-staggered signal sorting method and a device based on a correlation matching method, which are used for at least partially solving the technical problems.
(II) technical scheme
The invention provides a multi-parameter signal sorting method based on a correlation matching method, which comprises the following steps: s101, converting the multi-parameter signal pulse stream into a 01 sequence, and calculating the projection energy of the 01 sequence in a unit orthogonal subspace; s102, calculating the measurement length corresponding to each period value in the multi-spread signal pulse stream according to the projection energy; s103, calculating a detection threshold corresponding to each period value; s104, respectively comparing the measurement length corresponding to each period value with the detection threshold value, and forming a candidate set by taking the period values with the corresponding measurement lengths larger than the detection threshold value; s105, obtaining a first PRI value according to the divisor number of the period values in the set to be selected and the size of the measurement length corresponding to each period value; s106, according to the first PRI value, separating a first pulse sequence corresponding to the first PRI value from the multi-staggered signal pulse stream; and S107, repeating S101-S106 until the length of the residual pulse sequence in the multi-stagger signal pulse stream is smaller than the sorting threshold, and mixing the first pulse sequences with the same PRI value to obtain a stagger sequence.
Optionally, obtaining the first PRI value according to the divisor number of the period values in the candidate set and the size of the corresponding metric length includes: comparing the number of divisors of each period value in the set to be selected, wherein the period value with the largest divisor is a first PRI value; or, when the divisor number is the same, comparing the product of the measurement length corresponding to each period value in the candidate set and the corresponding period value, and setting the period value corresponding to the result of the maximum product as the first PRI value.
Optionally, separating, from the stream of multi-staggered signal pulses, a first pulse sequence corresponding to the first PRI value according to the first PRI value comprises: s301, calculating to obtain a first search range according to the first PRI value, searching the multi-stagger signal pulse stream by taking a first pulse of the multi-stagger signal pulse stream as a starting point of a second pulse sequence, calculating a second pulse of the second pulse sequence according to the first PRI value, and updating the PRI value to obtain a second PRI value; s302, a second search range is obtained through calculation according to a second PRI value, a multi-staggered signal pulse stream is searched by taking a second pulse as a starting point of a second pulse sequence, a third pulse of the second pulse sequence is calculated according to the second PRI value, a third PRI value is obtained through updating the PRI value, and a third search range is obtained through calculation according to the third PRI value; s303, repeating S302, and adjusting the search range until the maximum value of the search range is larger than or equal to the pulse arrival time of the latest arrival pulse in the multi-staggered signal pulse stream; s304, judging whether the second pulse sequence meets the separation condition, if so, storing the second pulse sequence, and repeating S301-S303 by taking the second pulse of the multi-parameter signal pulse stream as the starting point of the third pulse sequence; otherwise, directly taking the second pulse of the multi-parameter signal pulse stream as the starting point of the third pulse sequence, and repeating S301-S303; s305, when a certain pulse of the multiple staggered signal pulse stream is used as a starting point to search the multiple staggered signal pulse stream, and the number of the residual pulses does not meet the separation condition, stopping searching; s306, calculating the first-stage difference intervals of all stored pulse sequences, taking the average value as a sample, and finding out the center of the sample, wherein the pulse sequence corresponding to the center of the sample is the first pulse sequence.
Optionally, calculating the second pulse of the second pulse train from the first PRI value comprises: s401, judging the number of pulses of the multi-staggered signal pulse stream falling into a first search range; s402, if the pulse number is 0, calculating a fourth search range of the second pulse sequence, and repeating S401; if the number of pulses is 1, the pulse is the second pulse of the second pulse sequence; if the number of pulses is greater than 1, the difference between the pulse satisfying the condition and the first pulse of the second pulse train is calculated, and the pulse corresponding to the difference closest to the first PRI value is the second pulse.
Optionally, the determining whether the second pulse sequence satisfies the separation condition includes: calculating to obtain a sorting threshold according to the pulse arrival time of the latest arrival pulse in the multi-staggered signal pulse stream and the first PRI value; and comparing the pulse number of the second pulse sequence with the sorting threshold value, wherein the second pulse sequence larger than the sorting threshold value meets the separation condition.
Optionally, the multi-variance signal sorting method based on the correlation matching method further includes: and calculating a primary difference histogram of the staggered sequence, and estimating the sub-pulse interval of the staggered sequence according to the primary difference histogram.
Optionally, converting the multi-spread signal pulse stream into a 01 sequence, and calculating the projection energy of the 01 sequence in the unit orthogonal subspace includes: constructing a unit orthogonal matrix corresponding to each period value in the multi-parameter pulse stream according to the prior PRI value range; converting the multi-parameter signal pulse stream into a 01 sequence by interpolation according to the pulse arrival time; and segmenting the 01 sequence according to the pulse period, and calculating the projection energy of the 01 sequence in the unit orthogonal subspace according to the unit orthogonal matrix and the segmented 01 sequence.
Optionally, calculating the measurement length corresponding to each period value in the multi-spread signal pulse stream according to the projection energy comprises: and calculating to obtain each measurement length according to the infinite norm of the projection energy.
Optionally, calculating the detection threshold corresponding to each period value includes: and segmenting the observation time according to each period value, and calculating to obtain a detection threshold value according to the segmented observation time.
Another aspect of the present invention provides a multi-parameter signal sorting apparatus based on a correlation matching method, including: the projection energy acquisition module is used for converting the multi-parameter signal pulse stream into a 01 sequence and calculating the projection energy of the 01 sequence in the unit orthogonal subspace; the measurement length acquisition module is used for calculating the measurement length corresponding to each period value in the multi-staggered signal pulse stream according to the projection energy; the candidate period set acquisition module is used for calculating the detection threshold corresponding to each period value, respectively comparing the measurement length corresponding to each period value with the detection threshold, and selecting the period values with the corresponding measurement lengths larger than the detection threshold to form a candidate set; the PRI value acquisition module is used for obtaining a first PRI value according to the divisor number of the period values in the set to be selected and the size of the measurement length corresponding to each period value; the pulse sequence separation module is used for separating a first pulse sequence corresponding to a first PRI value from the multi-staggered signal pulse stream according to the first PRI value; and the staggered sequence acquisition module is used for commanding the modules to repeatedly execute respective actions until the length of the residual pulse sequence in the multi-staggered signal pulse stream is smaller than the sorting threshold value, and mixing the first pulse sequences with the same PRI value to obtain the staggered sequence.
(III) advantageous effects
The invention provides a multi-parameter signal sorting method based on a correlation matching method, which comprises the steps of converting a multi-parameter signal pulse stream into a 01 sequence, calculating and comparing the corresponding measurement length and detection threshold value of different periods, estimating the PRI value of a pulse sequence, further fusing a sequence retrieval technology and completing the separation of the pulse sequence. The sorting method provided can identify the frame period of the multiple staggered signals under the condition that the pulse stream has pulse missing, and further estimate the sub-pulse interval of the staggered sequence.
The period identification criterion provided by the multi-parameter signal sorting method based on the correlation matching method improves the accuracy of estimating the PRI value by comparing the divisor number of the period to be selected and comprehensively comparing the product of the measurement length and the corresponding period by utilizing the projection characteristic of the pulse sequence on the subspace.
Drawings
The above and other objects, features and advantages of the present invention will become more apparent from the following description of the embodiments of the present invention with reference to the accompanying drawings, in which:
fig. 1 schematically shows a flow chart of a multi-staggered signal sorting method based on a correlation matching method according to an embodiment of the invention;
FIG. 2 schematically illustrates a flow diagram of a method of calculating projection energy in accordance with an embodiment of the invention;
FIG. 3 schematically illustrates a diagram of a staggered TOA model according to an embodiment of the present invention;
FIG. 4 schematically illustrates a cycle identification criteria method flow diagram according to an embodiment of the invention;
FIG. 5 schematically shows a flow diagram of a sequence retrieval algorithm according to an embodiment of the invention;
FIG. 6 schematically illustrates a flow chart of a pulse sequence separation method according to an embodiment of the invention;
FIG. 7 schematically illustrates a first round sorting result graph according to an embodiment of the present invention;
FIG. 8 is a diagram schematically illustrating a second round of sorting results according to an embodiment of the present invention;
FIG. 9 schematically shows a third round of sorting results according to an embodiment of the invention;
FIG. 10 schematically illustrates a fourth round of sorting results according to an embodiment of the present invention;
FIG. 11 schematically shows a fifth round sorting result diagram according to an embodiment of the invention;
FIG. 12 schematically shows a sixth round of sorting results according to an embodiment of the invention;
FIG. 13 is a graph schematically illustrating sub-PRIs estimation results of the staggered pulse train 1 according to an embodiment of the present invention;
FIG. 14 is a graph schematically illustrating sub-PRIs estimation results of the spread pulse sequence 2 according to an embodiment of the present invention;
FIG. 15 is a flow diagram schematically illustrating a multi-staggered signal sorting method based on a correlation matching method according to another embodiment of the present invention;
fig. 16 schematically shows a block diagram of a multi-parameter signal sorting apparatus based on a correlation matching method according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to specific embodiments and the accompanying drawings.
It should be noted that in the drawings or description, the same drawing reference numerals are used for similar or identical parts. Features of the embodiments illustrated in the description may be freely combined to form new embodiments without conflict, and each claim may be individually referred to as an embodiment or features of the claims may be combined to form a new embodiment, and in the drawings, the shape or thickness of the embodiment may be enlarged and simplified or conveniently indicated. Further, elements or implementations not shown or described in the drawings are of a form known to those of ordinary skill in the art. Additionally, while exemplifications of parameters including particular values may be provided herein, it is to be understood that the parameters need not be exactly equal to the respective values, but may be approximated to the respective values within acceptable error margins or design constraints.
Unless a technical obstacle or contradiction exists, the above-described various embodiments of the present invention may be freely combined to form further embodiments, which are within the scope of the present invention.
Although the present invention has been described in connection with the accompanying drawings, the embodiments disclosed in the drawings are intended to be illustrative of preferred embodiments of the present invention and should not be construed as limiting the invention. The dimensional proportions in the figures are merely schematic and are not to be understood as limiting the invention.
Although a few embodiments of the present general inventive concept have been shown and described, it would be appreciated by those skilled in the art that changes may be made in these embodiments without departing from the principles and spirit of the general inventive concept, the scope of which is defined in the claims and their equivalents.
Fig. 1 schematically shows a flow chart of a multi-staggered signal sorting method based on a correlation matching method according to an embodiment of the invention.
According to the embodiment of the present invention, as shown in fig. 1, the multi-parameter signal sorting method based on the correlation matching method includes, for example:
and S101, converting the multi-parameter signal pulse stream into a 01 sequence, and calculating the projection energy of the 01 sequence in the unit orthogonal subspace.
And S102, calculating the measurement length corresponding to each period value in the multi-staggered signal pulse stream according to the projection energy.
S103, a detection threshold corresponding to each period value is calculated.
S104, comparing the measurement length corresponding to each period value with the detection threshold value, and forming a candidate set by taking the period values with the corresponding measurement lengths larger than the detection threshold value.
And S105, obtaining a first PRI value according to the divisor number of the period values in the set to be selected and the size of the measurement length corresponding to each period value.
S106, according to the first PRI value, separating a first pulse sequence corresponding to the first PRI value from the multi-staggered signal pulse stream.
And S107, repeating S101-S106 until the length of the residual pulse sequence in the multi-stagger signal pulse stream is smaller than the sorting threshold, and mixing the first pulse sequences with the same PRI value to obtain a stagger sequence.
FIG. 2 schematically shows a flow chart of a method of calculating projection energy according to an embodiment of the invention.
According to an embodiment of the present invention, as shown in fig. 1, the correlation projection energy may be calculated, for example, through operations S201 to S203.
S201, according to the prior PRI value range, constructing a unit orthogonal matrix corresponding to each period value in the multi-parameter pulse stream.
According to an embodiment of the invention, let e.g. the PRI have a priori range Pmin，Pmax]，
q∈N*Constructing an unit orthogonal matrix:
wherein q is the period of the pulse sequence,
is the euler function. Matrix array
And
the calculation results are as follows: when q is a prime number, the number of bits,
wherein the content of the first and second substances,
n-0., (q-1); when q is equal to pmWherein p is prime number, m belongs to N, and m is more than 1,
wherein, the symbol
The expression is for the function of rounding up,
when q is q1q2，gcd(q1，q2)＝1，q1＜q2，
Wherein n ═ n (n)1q2+n2q1)mod q，h1＝0，...，(q1-1)，n2＝0，...，(q2-1)；l＝(l1-1)q2+l2，
An identity matrix:
wherein the content of the first and second substances,
is one dimension of
Of diagonal matrix of diagonal elements of
Wherein the content of the first and second substances,
and S202, converting the multi-parameter signal pulse stream into a 01 sequence by interpolation according to the pulse arrival time.
Fig. 3 schematically shows a staggered sequence TOA model diagram according to an embodiment of the present invention.
According to an embodiment of the present invention, let t be assumedkIs the pulse arrival time, K1. The arrival time of a pulse can be represented by a set:
T＝[t1，t2，...，tk]T (6)
to obtain a uniformly sampled pulse sequence, we need to convert the TOA model into a 01 sequence by interpolation, and the period of the 01 sequence is proportional to the PRI value of the pulse sequence. The TOA model may be, for example, as shown in fig. 3. For simplicity, the present invention sets the proportionality coefficient to 1, for example. Assuming that the pulse arrival times are all integers, the 01 sequence conversion process is as follows:
and S203, segmenting the 01 sequence according to the pulse period, and calculating the projection energy of the 01 sequence in the unit orthogonal subspace according to the unit orthogonal matrix and the segmented 01 sequence.
According to an embodiment of the present invention, we can estimate the period of the 01 sequence by calculating the projected energy of the 01 sequence on the orthogonal subspace. For the sake of calculation, we choose a 2-norm representation of the projection energy:
wherein the content of the first and second substances,
representing vectors
Yang (Yang)
The inner product of (d). Assuming that the 01 sequence is T in length, the 01 sequence is
Segmenting according to a period q:
wherein the content of the first and second substances,
(symbol)
representing a floor function. Computing sequences
And matrix
Inner product of ith column:
wherein the content of the first and second substances,
due to the fact that
Is a unit matrix, therefore
And matrix
The projection energy in column i is calculated as follows:
assuming that the period of x (t) of a pulse sequence is q and the observation time is TT, the projection energy of the pulse sequence is calculated:
wherein the content of the first and second substances,
represents the ith column of the unit orthogonal matrix,
represents a 01 sequence, a symbol
Indicating approximately proportional.
In accordance with an embodiment of the present invention,
and matrix
The projection energy of (a) is expressed as follows:
using the infinite norm of the vector as the measurement length, obtaining the measurement length corresponding to each period value as
The risk threshold may be calculated as follows:
wherein, δ is an adjustable parameter, and the value thereof can be, for example, 0.1-0.3.
According to an embodiment of the invention, the pulse sequence is calculatedThe measuring length of the column and the corresponding detection threshold value are compared, the magnitude of the two is compared, and the period values corresponding to the pulse sequences with the measuring length larger than the detection threshold value are combined into a set Q, wherein Q is [ Q ═ Q [ [ Q ] Q1，q2，...]. If the period of the pulse sequence x (t) is q, q1Is a divisor of q, x (t) being in a period q according to equation (11)1The projection energy on the subspace of (a) is:
since q1 is a divisor of q, the projection energy calculated by equation (15) is also large, and in order to avoid identifying the divisor of a period as a period, we set a period identification criterion as shown in fig. 4.
FIG. 4 schematically illustrates a flow diagram of a cycle recognition criteria method according to an embodiment of the invention.
According to an embodiment of the invention, the set of preselected periods is determined to be Q, Q ═ Q1，q2，...]，qi( i 1, 2.) is the candidate period value. The PRI value estimated for each round of sorting can be obtained, for example, through operations S401 to S402.
S401, comparing the number of divisors of each period value in the to-be-selected set, wherein the period value with the largest divisor is the first PRI value. For example, for each element in the set O, find the respective divisor among the elements contained in O, if the jth element qjWith the largest divisor of (1), qjIs the PRI value estimated for this round of sorting.
S402, when the divisor number is the same, comparing the product of the measurement length corresponding to each period value in the to-be-selected set and the corresponding period value, wherein the period value corresponding to the result of the maximum product is the first PRI value. In this embodiment, if the jth element q isjWith the ith element qiIf the divisors are the same, the metric lengths corresponding to the periods need to be compared. Considering that the length of the measurement is approximately proportional to the pulse sequence period, if q isjQ of corresponding metric lengthjTimes greater than qiQ of corresponding metric lengthiMultiple, then qjIs PRI value estimated for this round of sortingOtherwise, then qiIs the PRI value estimated for this round of sorting.
Fig. 5 schematically shows a flow chart of a sequence retrieval algorithm according to an embodiment of the invention.
Fig. 6 schematically shows a flow chart of a pulse sequence separation method according to an embodiment of the invention.
According to the embodiments of the present invention, as shown in fig. 5 and fig. 6, the pulse sequence corresponding to the estimated PRI value can be found and separated by, for example, operating SS601 to SS 606.
And SS601, obtaining a first search range by calculation according to the first PRI value, searching the multi-stagger signal pulse stream by taking the first pulse of the multi-stagger signal pulse stream as the starting point of the second pulse sequence, calculating the second pulse of the second pulse sequence according to the first PRI value, and updating the PRI value to obtain the second PRI value.
According to an embodiment of the present invention, in the proposed sequence retrieval algorithm, the estimated PRI value is taken as the initial PRI value. The PRI value is updated once per pulse separation until the separation of the pulse sequence is completed. The sequence of pulse stream arrival times is denoted by s, the estimated PRI value is denoted by PRI, and the relevant parameters are initialized, e.g. i1, toan s (i), PRI, kmin，kmaxSck ═ 1, tt ═ 1, tmp (tt) ═ toan, where the variable k ismin、kmaxAnd sck to determine the single search range of the pulses, tt for pulse counting, tmp to store the separated pulse sequence. Pulse search range:
wherein, TOAminAnd TOAmaxRespectively representing the minimum and maximum values of the pulse search range. And toan represents the time of arrival of the ith pulse.
Calculating pulse single search results:
cnn＝find((s≥TOAmin)&(s≤TOAmax)) (17)
where cnn denotes the position of the element in the pulse stream s that satisfies the search criteria, and the length (cnn) values may be 0, 1, and greater than 1, indicating no matching pulse, one matching pulse, and more than one matching pulse, respectively.
And if the number of the first matching pulses is zero, adjusting the first search range to obtain a second search range.
According to the embodiment of the present invention, if length (cnn) ═ 0, then sck ' ═ sck +1, sck ' is adjusted to the adjusted variable sck, that is, the pulse sequence is searched again at different pulse intervals, so that the situation that the search cannot be completed due to the loss of pulses can be overcome, and the pulse search range in formula (16) is recalculated with the variable sck '.
And if the number of the first matching pulses is one, accessing the matched pulses into the pulse sequence, calculating a primary difference sequence of the sequence, and averaging to obtain a second PRI value.
According to an embodiment of the present invention, if length (cnn) is 1, the updated PRI value is:
where PRI' is the updated PRI value.
And if the number of the first matching pulses is more than one, selecting a second pulse with a pulse interval closest to the first PRI value from the plurality of matching pulses, accessing the matched pulse into the pulse sequence, calculating a first-level difference sequence of the sequence, and averaging to obtain a second PRI value.
According to an embodiment of the present invention, if length (cnn) > 1, i.e. there are multiple pulses satisfying the condition, one pulse is selected from the multiple pulses such that the adjacent pulse interval is closest to the estimated PRI value, and then the PRI value is updated as:
where PRI' is the updated PRI value and s (tpnum) represents the arrival time of the next separated pulse, the time interval between the next separated pulse and the adjacent pulse being closest to the first PRI value.
And SS602, calculating to obtain a second search range according to the second PRI value, searching the multi-staggered signal pulse stream by using the second pulse as the starting point of the second pulse sequence, calculating a third pulse of the second pulse sequence according to the second PRI value, updating the PRI value to obtain a third PRI value, and calculating to obtain a third search range according to the third PRI value.
According to the embodiment of the invention, the updated PRI value and the corresponding search range are continuously calculated according to the formula (16) to the formula (19).
And SS603, repeating SS602, and adjusting the search range until the maximum value of the search range is greater than or equal to the pulse arrival time of the latest arrival pulse in the multi-staggered signal pulse stream.
According to the embodiment of the present invention, according to the updated PRI value, the relevant initial parameters are adjusted to tt' ═ tt +1, toan ═ s (cnn) or toan ═ s (tpnum), tmp (tt) ═ toan, and sck ═ 1, and the steps of SS602 are repeated. Until TOAmax＞max(s)*kmax. Where max(s) is the pulse arrival time of the latest arriving pulse.
SS604, judge whether the second pulse sequence meets the separation condition, if meet, store the second pulse sequence, regard second pulse of the pulse stream of the multiple-parameter signal as the starting point of the third pulse sequence, repeat SS601-SS 603; otherwise, SS601-SS603 are repeated with the second pulse of the multi-staggered signal pulse stream directly as the start of the third pulse sequence.
According to an embodiment of the present invention, the sorting threshold TE may be represented by the formula:
the sorting threshold TE is a threshold value of the pulse sequence including the number of pulses, σ is a proportionality coefficient, and the value range is, for example, 0.5-0.7.
When tt > TE, the pulse sequence tmp is stored. And completing a search, updating the sequence starting point, i is i +1, and toan is s (i).
And SS605 for stopping the search when the multi-staggered signal pulse stream is searched by using a certain pulse of the multi-staggered signal pulse stream as a starting point and the number of residual pulses does not satisfy the separation condition.
According to the embodiment of the invention, after the sequence starting point updating is completed, the initialization parameters are adjusted to tt ═ 1, sck ═ 1, tmp (tt) ═ toan, and the steps SS601-SS603 are repeated, when i > length(s) -TE, it is stated that s (i) is taken as the sequence starting point at this time, the rest pulses will not satisfy the separation condition any more, and the search is completed.
And SS606, calculating the first-stage difference intervals of all the stored pulse sequences, averaging, taking the average value as a sample, finding out the center of the sample, wherein the pulse sequence corresponding to the center of the sample is the first pulse sequence.
According to an embodiment of the present invention, tmp is used to store the searched pulse sequences, and there may be a plurality of pulse sequences that satisfy the separation condition after the search process is finished. And calculating the primary difference of each pulse sequence according to the stored pulse sequences, averaging, taking the average value as the PRI value of the pulse sequence, and finding the pulse sequence with the PRI value closest to the initial PRI value, wherein the pulse sequence is the output of the current search.
According to the embodiment of the present invention, the multi-variance signal sorting method based on the correlation matching method further includes, for example: and calculating a primary difference histogram of the staggered sequence, and estimating the sub-pulse interval of the staggered sequence according to the primary difference histogram.
According to an embodiment of the present invention, in order to verify the validity of the present invention, the following simulation was performed in the present embodiment. The radar signal simulation model may be, for example, a spread signal TOA model as shown in fig. 3, and the main simulation parameters are shown in table 1:
TABLE 1 Mixed pulse sequence 1 parameter Table
| Parameter(s) | |
|
| Frame period (μ s) | 27 | 35 |
| sub-PRIs(μs) | 5，9，13 | 6，12，17 |
| Noise error | ±3％ | ±3％ |
| Pulse loss rate | ±6％～±16％ | ±6％～±16％ |
| Total Observation time (μ s) | 5000 | 5000 |
The simulation results are shown in fig. 7 to 14. The sorting in this embodiment has a total of six rounds as shown in fig. 7 to 12. The first round of sorting has 3 period values exceeding the detection threshold, and according to the period identification criterion, the first round of PRI estimation results are represented by PRI1 ═ 27 μ s. The second round of sorting has 2 period values exceeding the detection threshold, and according to the period identification criterion, the second round of PRI estimation results are obtained when the PRI2 is 27 μ s. The third round of sorting has 2 period values exceeding the detection threshold, and according to the period identification criterion, the third round of PRI estimation results is represented by PRI3 ═ 27 μ s. The fourth round sorting has 1 period value exceeding the detection threshold, and the PRI4 is 35 μ s, which is the fourth round PRI estimation result. The fifth round of sorting has 1 period value exceeding the detection threshold, and the result of the fifth round of PRI estimation is PRI5 ═ 35 μ s. The sixth round of sorting has 1 period value exceeding the detection threshold, and the result of the sixth round of PRI estimation is PRI6 ═ 35 μ s. After the six sorting rounds are completed, the staggered sequence sub-pulse intervals are estimated, and the estimation results of the sub-pulse intervals are shown in fig. 13 and 14. From the simulation results, it can be known that under the condition that the pulse stream has pulse missing, the improved Correlation Matching Method (CMM) can identify the frame period of the multi-staggered pulse sequence and estimate the sub-pulse interval of the staggered pulse sequence.
Fig. 15 schematically shows a flow chart of a multi-staggered signal sorting method based on a correlation matching method according to another embodiment of the invention.
According to the embodiment of the present invention, as shown in fig. 15, the multi-variance signal sorting method based on the correlation matching method includes, for example:
s1501, constructing a unit Orthogonal Ramanujan Subspace (Uniform Orthogonal Ramanujan Subspace, UORS).
S1502, convert the input pulse stream into 01 sequences, and calculate the projection energies of the 01 sequences on different UORS.
S1503, calculating a detection threshold of the improved CMM, comparing the measurement lengths corresponding to different periods with the detection threshold, and forming a candidate set by period values with the measurement lengths exceeding the detection threshold.
S1504, using a period identification criterion to screen a period value from the period set to be selected obtained in S1503, where the period value is the PRI estimated in the current round.
S1505, according to the PRI estimated in S1504, completes the separation of the pulse sequence corresponding to the PRI. After the pulse sequence separation is finished, judging whether the number of the remaining pulses exceeds a sorting threshold, if so, returning to the step S1502, and starting the next round of sorting; otherwise, go to S1506.
And S1506, mixing the pulse sequences corresponding to the same PRI value, and estimating sub-PRIs of the mixed pulse sequences.
In summary, the embodiment of the present invention provides a multi-spread signal sorting method based on a correlation matching method. By converting the multi-stagger signal pulse stream into a 01 sequence, fusing a sequence retrieval technology, comparing the calculated measurement length with a detection threshold value, and enabling multi-stagger signals with different periods to correspond to different detection threshold values, the pulse sequences corresponding to different PRI values can be separated, the frame period of the multi-stagger signals can be identified under the condition that the pulse stream has pulse loss, and the sub-pulse interval of the stagger sequence can be further estimated.
Fig. 16 schematically shows a block diagram of a multi-parameter signal sorting apparatus based on a correlation matching method according to an embodiment of the present invention.
In another aspect, an embodiment of the present invention provides a multi-variance signal sorting apparatus based on a correlation matching method, as shown in fig. 16, the multi-variance signal sorting apparatus 1600 based on the correlation matching method includes:
and a projection energy obtaining module 1601, configured to convert the multi-dispersion signal pulse stream into a 01 sequence, and calculate projection energy of the 01 sequence in the unit orthogonal subspace.
A measurement length obtaining module 1602, configured to calculate the measurement length corresponding to each period value in the multiple-staggered-signal pulse stream according to the projection energy.
A candidate period set obtaining module 1603, configured to calculate a detection threshold corresponding to each period value, compare the measurement length corresponding to each period value with the size of the detection threshold, and select a period value with the measurement length greater than the detection threshold to form a candidate set.
The PRI value obtaining module 1604 is configured to obtain a first PRI value according to the divisor number of the period values in the to-be-selected set and the size of the measurement length corresponding to each period value.
A pulse sequence separation module 1605, configured to separate a first pulse sequence corresponding to the first PRI value from the multi-spread signal pulse stream according to the first PRI value.
A stagger sequence obtaining module 1606, configured to instruct the modules to repeatedly perform their respective actions until the length of the remaining pulse sequences in the multiple-stagger signal pulse stream is smaller than the sorting threshold, and mix the first pulse sequences with the same PRI value to obtain a stagger sequence.
The product embodiment is similar to the method embodiment in portions where details are not given, and please refer to the method embodiment, which is not described herein again.
It should be understood that the specific order or hierarchy of steps in the processes disclosed is an example of exemplary approaches. Based upon design preferences, it is understood that the specific order or hierarchy of steps in the processes may be rearranged without departing from the scope of the present disclosure. The accompanying method claims present elements of the various steps in a sample order, and are not intended to be limited to the specific order or hierarchy.
In the foregoing detailed description, various features are grouped together in a single embodiment for the purpose of streamlining the disclosure. This method of disclosure is not to be interpreted as reflecting an intention that the claimed embodiments of the subject matter require more features than are expressly recited in each claim. Rather, as the following claims reflect, invention lies in less than all features of a single disclosed embodiment. Thus, the following claims are hereby expressly incorporated into the detailed description, with each claim standing on its own as a separate preferred embodiment of the invention.
The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are only exemplary embodiments of the present invention, and are not intended to limit the present invention, and any modifications, equivalents, improvements and the like made within the spirit and principle of the present invention should be included in the protection scope of the present invention. Furthermore, to the extent that the term "includes" is used in either the detailed description or the claims, such term is intended to be inclusive in a manner similar to the term "comprising" as "comprising" is interpreted when employed as a transitional word in a claim. Furthermore, any use of the term "or" in the specification of the claims is intended to mean a "non-exclusive or".
Claims (10)
1. A multi-spread signal sorting method based on a correlation matching method is characterized by comprising the following steps:
s101, converting the multi-parameter signal pulse stream into a 01 sequence, and calculating the projection energy of the 01 sequence in a unit orthogonal subspace;
s102, calculating the measurement length corresponding to each period value in the multi-spread signal pulse stream according to the projection energy;
s103, calculating a detection threshold corresponding to each period value;
s104, respectively comparing the measurement length corresponding to each period value with the detection threshold value, and forming a candidate set by taking the period values with the corresponding measurement lengths larger than the detection threshold value;
s105, obtaining a first PRI value according to the divisor number of the period values in the to-be-selected set and the size of the measurement length corresponding to each period value;
s106, according to the first PRI value, separating a first pulse sequence corresponding to the first PRI value from the multi-spread signal pulse stream;
and S107, repeating S101-S106 until the length of the residual pulse sequence in the multi-stagger signal pulse stream is smaller than the sorting threshold, and mixing the first pulse sequences with the same PRI value to obtain a stagger sequence.
2. The method of claim 1, wherein the deriving the first PRI value according to the divisor number of period values in the candidate set and the size of the corresponding metric length comprises:
comparing the number of divisors of all period values in the to-be-selected set, wherein the period value with the largest divisor is the first PRI value; alternatively, the first and second electrodes may be,
and when the divisor number is the same, comparing the product of the measurement length corresponding to each period value in the to-be-selected set and the corresponding period value, wherein the period value corresponding to the result of the maximum product is the first PRI value.
3. The correlation matching method-based multi-spread signal sorting method according to claim 1, wherein said separating a first pulse sequence corresponding to the first PRI value from the multi-spread signal pulse stream according to the first PRI value comprises:
s301, calculating to obtain a first search range according to the first PRI value, searching the multiple-disparity signal pulse stream by taking a first pulse of the multiple-disparity signal pulse stream as a starting point of a second pulse sequence, calculating a second pulse of the second pulse sequence according to the first PRI value, and updating the PRI value to obtain a second PRI value;
s302, calculating to obtain a second search range according to the second PRI value, searching the multi-dispersion signal pulse stream by taking the second pulse as the starting point of the second pulse sequence, calculating a third pulse of the second pulse sequence according to the second PRI value, updating the PRI value to obtain a third PRI value, and calculating to obtain a third search range according to the third PRI value;
s303, repeating S302, and adjusting the search range until the maximum value of the search range is greater than or equal to the pulse arrival time of the latest arrival pulse in the multi-staggered signal pulse stream;
s304, judging whether the second pulse sequence meets the separation condition, if so, storing the second pulse sequence, and repeating S301-S303 by taking the second pulse of the multi-staggered signal pulse stream as the starting point of a third pulse sequence; otherwise, directly taking the second pulse of the multi-staggered signal pulse stream as the starting point of the third pulse sequence, and repeating S301-S303;
s305, when the multi-staggered signal pulse stream is searched by taking a certain pulse of the multi-staggered signal pulse stream as a starting point and the number of residual pulses does not meet the separation condition, stopping searching;
s306, calculating the first-stage difference intervals of all stored pulse sequences, taking an average value, taking the average value as a sample, and finding out a sample center, wherein the pulse sequence corresponding to the sample center is the first pulse sequence.
4. The correlation matching method-based multi-spread signal sorting method according to claim 3, wherein the calculating the second pulse of the second pulse sequence according to the first PRI value comprises:
s401, judging the number of pulses of the multi-staggered signal pulse stream falling into the first search range;
s402, if the pulse number is 0, calculating a fourth search range of the second pulse sequence, and repeating S401;
if the number of pulses is 1, the pulse is the second pulse of the second pulse sequence;
and if the number of the pulses is more than 1, calculating the difference value between the pulses meeting the condition and the first pulse of the second pulse sequence, wherein the pulse corresponding to the difference value closest to the first PRI value is the second pulse.
5. The correlation matching method-based multi-spread signal sorting method according to claim 3, wherein the determining whether the second pulse sequence satisfies a separation condition comprises:
calculating to obtain the sorting threshold according to the pulse arrival time of the latest arrival pulse in the multi-staggered signal pulse stream and the first PRI value;
comparing the number of pulses of the second pulse sequence with the magnitude of the sorting threshold, the second pulse sequence larger than the sorting threshold satisfying the separation condition.
6. The correlation matching method based multi-spread signal sorting method according to claim 1, wherein the correlation matching method based multi-spread signal sorting method further comprises:
and calculating a primary difference histogram of the staggered sequence, and estimating the sub-pulse interval of the staggered sequence according to the primary difference histogram.
7. The correlation matching method-based multi-spread signal sorting method according to claim 1, wherein the converting the multi-spread signal pulse stream into 01 sequences, and the calculating the projection energy of the 01 sequences in the unit orthogonal subspace comprises:
constructing a unit orthogonal matrix corresponding to each period value in the multi-parameter signal pulse stream according to the prior PRI value range;
converting the multi-spread signal pulse stream into the 01 sequence by interpolation according to pulse arrival time;
and segmenting the 01 sequence according to the pulse period, and calculating the projection energy of the 01 sequence in the unit orthogonal subspace according to the unit orthogonal matrix and the segmented 01 sequence.
8. The correlation matching method-based multi-spread signal sorting method according to claim 1, wherein said calculating the measurement lengths corresponding to the period values in the multi-spread signal pulse streams according to the projection energies comprises:
and calculating to obtain each measurement length according to the infinite norm of the projection energy.
9. The correlation matching method-based multi-spread signal sorting method according to claim 1, wherein the calculating the detection threshold corresponding to each of the period values comprises:
and segmenting the observation time according to each period value, and calculating to obtain the detection threshold according to the segmented observation time.
10. A multi-parameter signal sorting device based on a correlation matching method is characterized by comprising the following steps:
the projection energy acquisition module is used for converting the multi-parameter signal pulse stream into a 01 sequence and calculating the projection energy of the 01 sequence in a unit orthogonal subspace;
a measurement length obtaining module, configured to calculate, according to the projection energy, a measurement length corresponding to each period value in the multiple-spread signal pulse stream;
a candidate period set obtaining module, configured to calculate a detection threshold corresponding to each period value, compare the measurement length corresponding to each period value with the size of the detection threshold, and form a candidate set by using the period values with the measurement lengths larger than the detection threshold;
a PRI value obtaining module, configured to obtain a first PRI value according to the divisor number of the period values in the to-be-selected set and the size of the measurement length corresponding to each period value;
a pulse sequence separation module, configured to separate a first pulse sequence corresponding to the first PRI value from the multiple-spread-signal pulse stream according to the first PRI value;
and the staggered sequence acquisition module is used for commanding the modules to repeatedly execute respective actions until the length of the residual pulse sequence in the multi-staggered signal pulse stream is smaller than a sorting threshold value, and mixing the first pulse sequences with the same PRI value to obtain a staggered sequence.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110563419.8A CN113296071B (en) | 2021-05-21 | 2021-05-21 | Multi-spread signal sorting method and device based on correlation matching method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110563419.8A CN113296071B (en) | 2021-05-21 | 2021-05-21 | Multi-spread signal sorting method and device based on correlation matching method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113296071A true CN113296071A (en) | 2021-08-24 |
| CN113296071B CN113296071B (en) | 2023-01-03 |
Family
ID=77324302
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110563419.8A Active CN113296071B (en) | 2021-05-21 | 2021-05-21 | Multi-spread signal sorting method and device based on correlation matching method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113296071B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109683143A (en) * | 2019-03-07 | 2019-04-26 | 西安电子科技大学 | Radar Signal Sorting Method, device, computer equipment and can storage medium |
| CN109839619A (en) * | 2019-03-15 | 2019-06-04 | 北京应用物理与计算数学研究所 | Based on radar signal rough segmentation choosing method, system and the storage medium for adaptively dividing bucket |
| EP3534178A1 (en) * | 2018-02-26 | 2019-09-04 | MBDA Deutschland GmbH | Radar apparatus with antenna array and method for source location using a two-dimensional antenna array |
| CN110334471A (en) * | 2019-07-17 | 2019-10-15 | 吉林工程技术师范学院 | The PRI subcycle extracting method of non-uniform series and arteries and veins group non-uniform series |
| CN111427046A (en) * | 2020-03-27 | 2020-07-17 | 西安交通大学 | Terahertz pulse echo positioning method for improving detection precision |
| CN112198481A (en) * | 2020-10-09 | 2021-01-08 | 哈尔滨工程大学 | PRI jitter radar signal sorting method under pulse loss aliasing condition |
| CN112415504A (en) * | 2020-11-02 | 2021-02-26 | 中国科学院空天信息创新研究院 | Radar target tracking method and device |
-
2021
- 2021-05-21 CN CN202110563419.8A patent/CN113296071B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP3534178A1 (en) * | 2018-02-26 | 2019-09-04 | MBDA Deutschland GmbH | Radar apparatus with antenna array and method for source location using a two-dimensional antenna array |
| CN109683143A (en) * | 2019-03-07 | 2019-04-26 | 西安电子科技大学 | Radar Signal Sorting Method, device, computer equipment and can storage medium |
| CN109839619A (en) * | 2019-03-15 | 2019-06-04 | 北京应用物理与计算数学研究所 | Based on radar signal rough segmentation choosing method, system and the storage medium for adaptively dividing bucket |
| CN110334471A (en) * | 2019-07-17 | 2019-10-15 | 吉林工程技术师范学院 | The PRI subcycle extracting method of non-uniform series and arteries and veins group non-uniform series |
| CN111427046A (en) * | 2020-03-27 | 2020-07-17 | 西安交通大学 | Terahertz pulse echo positioning method for improving detection precision |
| CN112198481A (en) * | 2020-10-09 | 2021-01-08 | 哈尔滨工程大学 | PRI jitter radar signal sorting method under pulse loss aliasing condition |
| CN112415504A (en) * | 2020-11-02 | 2021-02-26 | 中国科学院空天信息创新研究院 | Radar target tracking method and device |
Non-Patent Citations (1)
| Title |
|---|
| 安然等: "基于AR模型的水流速雷达信号处理方法", 《电子测量技术》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113296071B (en) | 2023-01-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109031215B (en) | Method for estimating pulse repetition period of staggered radar under high pulse loss | |
| CN111445108B (en) | Data-driven power distribution network line variation relation diagnosis method, device and system | |
| CN112986928B (en) | Signal sorting multi-source fusion processing method in complex electromagnetic environment | |
| CN109839619B (en) | Radar signal rough sorting method and system based on self-adaptive bucket separation and storage medium | |
| CN108985373B (en) | Multi-sensor data weighting fusion method | |
| CN110488226B (en) | Underwater target positioning method and device | |
| CN107205226B (en) | Indoor positioning and tracking method and system based on channel classification | |
| CN105445699A (en) | Range measurement method and range measurement system based on non line-of-sight error elimination | |
| CN113296071B (en) | Multi-spread signal sorting method and device based on correlation matching method | |
| Franchi et al. | Statistical properties of the maximum Lyapunov exponent calculated via the divergence rate method | |
| CN111238480B (en) | NLOS (non line of sight) identification method based on distance residual error and application of NLOS identification method in indoor positioning | |
| CN108680897A (en) | A kind of indoor orientation method, device, electronic equipment and storage medium | |
| CN109444899B (en) | Track association method based on pure angle information | |
| CN106488554B (en) | Fingerprint database establishing method and system | |
| CN115980689A (en) | Point cloud detection-based radiation source signal sorting method, device, equipment and medium | |
| CN113075620B (en) | Signal sorting method based on multi-station time difference grid clustering | |
| CN104462826A (en) | Method for detecting and measuring multi-sensor evidence conflicts on basis of matrix singular value decomposition | |
| Zhao et al. | A multi-station signal sorting method based on TDOA grid clustering | |
| Timoshenko et al. | Detection and estimation of parameters of a random process set in multi-Scanning radar observation based on the" track-before-detect" methods | |
| US9514361B2 (en) | Classification of range profiles | |
| CN114779225B (en) | Iterative grid-based distributed radar target detection method and device | |
| JP3346326B2 (en) | Pulse train classification device | |
| CN117310636B (en) | Fixed pulse repetition interval measurement method, device and medium | |
| CN113702919B (en) | Method and device for estimating PRI value and extracting pulse sequence | |
| Fhima et al. | Fast change point analysis on the Hurst index of piecewise fractional Brownian motion |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |