CN111624575B - Method and system for rapidly extracting unknown radar target pulse sequence - Google Patents

Method and system for rapidly extracting unknown radar target pulse sequence Download PDF

Info

Publication number
CN111624575B
CN111624575B CN202010747915.4A CN202010747915A CN111624575B CN 111624575 B CN111624575 B CN 111624575B CN 202010747915 A CN202010747915 A CN 202010747915A CN 111624575 B CN111624575 B CN 111624575B
Authority
CN
China
Prior art keywords
pulse
value
doa
peak
block
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.)
Active
Application number
CN202010747915.4A
Other languages
Chinese (zh)
Other versions
CN111624575A (en
Inventor
孟祥豪
陈晶
安永旺
毛毅
王红军
段永胜
李媛丽
张奎
王昊
张坤峰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN202010747915.4A priority Critical patent/CN111624575B/en
Publication of CN111624575A publication Critical patent/CN111624575A/en
Application granted granted Critical
Publication of CN111624575B publication Critical patent/CN111624575B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • 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/28—Details of pulse systems
    • G01S7/2813—Means providing a modification of the radiation pattern for cancelling noise, clutter or interfering signals, e.g. side lobe suppression, side lobe blanking, null-steering arrays

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention discloses a method for quickly extracting an unknown radar target pulse sequence, which comprises the following steps: acquiring a staggered full pulse sequence; obtaining a DOA parameter peak value of each pulse data in the staggered full pulse sequence and carrying out DOA blocking processing; obtaining a PW parameter peak value of each pulse data in each DOA block and carrying out PW block processing; acquiring an RF parameter peak value of each pulse data in each PW block and performing RF block; judging whether the absolute value of the difference between the RF parameter value of the unknown radar target pulse and the RF parameter peak value is smaller than a threshold value, if so, sequentially performing RF-PRI projection sample image extraction and unknown radar target pulse matching processing on the first pulse data block; and if not, sequentially carrying out PRI projection sample image extraction and unknown radar target pulse matching processing on the second pulse data block. The invention also discloses a storage medium and a system. The pulse sample graph extracted by the method has high accuracy, small calculated amount and high operation efficiency.

Description

Method and system for rapidly extracting unknown radar target pulse sequence
Technical Field
The invention belongs to the technical field of radar pulse stream processing, and particularly relates to a method and a system for quickly extracting an unknown radar target pulse sequence.
Background
In the dense radar pulse stream, pulse signals of a radar in a certain working mode are quickly screened out, signal parameters which are as comprehensive and accurate as possible are obtained, precious time can be won, the target can be interfered in a targeted and quick mode, meanwhile, the information such as parameters, positions, activity conditions, startup and shutdown rules and the like of key targets are updated in real time, the threat level and the electromagnetic environment situation of the radar target can be conveniently and quickly identified, and quick threat warning is realized.
Aiming at the situation that the type of radar signals contained in full-pulse data is unknown, a correlation method for quickly extracting radar signals from the existing radar signals, such as a radar signal quick extraction method based on a self-extraction pulse sample diagram, is characterized in that the method is based on the pulse sample diagram, utilizes the matching correlation among radar signal pulses of the same kind, and screens radar radiation source signals which can be described by the pulse sample diagram while realizing the self-extraction of the pulse sample diagram through the integral translation of the pulses. A pulse sample diagram is extracted by using a cyclic correlation method, and four-dimensional characteristic parameters of a single pulse of a radar signal, namely RF, PW, PRI and MOP, are extracted. The matching sorting of the signals is carried out by using the characteristic parameters, but the method has the following defects:
1) when the pulse sample image is extracted, the similarity measurement of two pulses needs to calculate all characteristic parameters, and for a radar signal (such as an RF agile type) with a certain characteristic parameter changing randomly, the change of the RF parameter has no rule, so that the accuracy of the similarity measurement is influenced, and the extraction accuracy is reduced.
2) When the matching of the two pulses is carried out, all the feature parameter similarities are calculated every time, the calculation amount is large, and the consumed time is long.
Disclosure of Invention
One of the objectives of the present invention is to provide a method for rapidly extracting an unknown radar target pulse sequence, where an extracted pulse sample graph has high accuracy, strong pertinence to different variation types of unknown radar target signal parameters, small calculation amount, and high operation efficiency
It is a second object of the present invention to provide a storage medium.
The third object of the present invention is to provide a system for rapidly extracting an unknown radar target pulse sequence.
In order to achieve one of the purposes, the invention adopts the following technical scheme:
a rapid extraction method of an unknown radar target pulse sequence comprises the following steps:
step one, acquiring a staggered full pulse sequence;
step two, acquiring a signal arrival azimuth DOA parameter peak value of each pulse data in the staggered full pulse sequence and performing DOA blocking processing;
step three, acquiring a pulse width PW parameter peak value of each pulse data in each DOA block and carrying out PW block processing; step four, obtaining a carrier frequency RF parameter peak value of each pulse data in each PW block and carrying out RF block processing to obtain a first pulse data block corresponding to the RF parameter peak value and a second pulse data block corresponding to a non-RF parameter peak value;
step five, judging whether the absolute value of the difference between the RF parameter value of the unknown radar target pulse and the RF parameter peak value is smaller than a threshold value, if so, sequentially performing RF-PRI projection sample image extraction and unknown radar target pulse matching processing on the first pulse data block; and if not, sequentially carrying out PRI projection sample image extraction and unknown radar target pulse matching processing on the second pulse data block.
Further, the specific implementation process of the step two is as follows:
step 21, reading the pulse data in the interleaved full pulse sequence and sequencing according to the DOA parameter value to obtain the maximum value DOA of the DOA parameter valuemaxAnd minimum value DOAmin；
Step 22, according to the preset DOA value resolution alphagateAnd maximum value of DOA parameter value DOAmaxAnd minimum value DOAminCalculating the DOA peak value statistical initial value DOAbeginDOA peak value statistic end value DOAendAnd DOA peak statistical range DOAwidthAnd DOA Peak statistical Window number n _ DOAwidth；
Step 23, setting the initial value of the DOA peak value statistical vector serial number i as 1, and setting the ith DOA peak value statistical vector peak _ DOA [ i [ ]]Is 0; DOA parameter comparison window alphawindowIs initialized to DOAbeginThe initial value of the DOA peak value breakpoint pulse sequence number k is 1;
step 24, setting the initial value of the pulse serial number j selected from the staggered full pulse sequence as 1;
step 25, judging whether i is larger than n _ DOAwidthIf yes, go to step 27; if not, go to step 26;
step 26, when j is<N and alphawindow<x1(j)≤αwindow+αgateThen peak _ DOA [ i ]]Adding 1 and then assigning peak _ DOA [ i]J is added with 1 and then is given to j, and the step 25 is returned; when j is>N, then go to step 27; when x is1(j)>αwindow+αgateAdding 1 to k, adding 1 to j, and adding 1 to j, alphawindowPlus alphagateAfter giving alphawindowAdding 1 to i, assigning i, and returning to the step 25;
wherein x is1(j) The DOA parameter value corresponding to the jth pulse; n is the number of pulses in the staggered full pulse sequence;
step 27, counting the number M of pulses with the DOA parameter value being the wave crest in the DOA peak value statistical vector sequence; and the interleaved full pulse sequence is divided into M DOA blocks.
Further, the specific implementation process of the step three is as follows:
step 31, obtaining the number of pulses N _ DOA in each DOA blockpid(ii) a Sequencing the pulse data in each DOA block according to the pulse width PW parameter value to obtain the maximum PW of the pulse PW parameter values in the corresponding DOA blockmaxAnd minimum value PWmin(ii) a Step 32, according to the preset PW value resolution taugateAnd the maximum value PW of the pulse PW parameter values in each DOA blockmaxAnd minimum value PWminCalculating the PW peak value statistic initial value PW in the corresponding DOA blockbeginPW peak value statistical end value PWendAnd PW peak statistical range PWwidthAnd the number of PW peak statistical windows n _ PWwidth；
Step 33, setting the initial value of the PW peak value statistical vector serial number i ' in each DOA block to be 1, and setting the ith ' PW peak value statistical vector peak _ PW [ i ']Is 0; PW parameter comparison window τwindowHas an initial value of PWbeginThe initial value of PW peak breakpoint pulse number k' is 1;
step 34, setting the initial value of the pulse serial number j' selected from the staggered full pulse sequence in the corresponding DOA block as 1;
step 35, judging whether i' is larger than n _ PWwidthIf yes, go to step 37; if not, go to step 36;
step 36, when j'<N_DOApidAnd τwindow<x2(j')≤τwindow+τgateThen peak _ PW [ i']Adding 1 and then giving peak _ PW [ i']Adding 1 to j ', and then assigning j', and returning to the step 35; when j'>N_DOApidThen go to step 37; when x is2(j')>τwindow+ τgateAdding 1 to k ', j ', and then adding 1 to j ', tauwindowPlus taugatePost-assigning to tauwindowAdding 1 to i ', assigning i', and returning to the step 35;
wherein x is2(j') is a PW parameter value corresponding to the jth pulse in the corresponding DOA block;
step 37, counting the number L of pulses with wave crest as the PW parameter value in the PW peak value statistical vector sequence in each DOA block; and dividing the corresponding DOA block-wise intra-interleaved full pulse sequence into L PW blocks.
Further, the specific implementation process of the step four is as follows:
step 41, obtaining the number N _ PW of pulses in each PW blockpid(ii) a And sequencing the pulse data in each PW block according to the pulse width RF parameter valueObtaining the maximum value RF of the pulse RF parameter value in the corresponding PW blockmaxAnd minimum value RFmin；
Step 42, resolution f according to the preset RF valuegateAnd maximum value RF of the pulse RF parameter value within each PW blockmaxAnd minimum value RFminCalculating the statistical initial value RF of the RF peak value in the corresponding PW blockbeginStatistical end of RF peak value RFendAnd RF peak statistical range RFwidthAnd number of RF peak statistics windows n _ RFwidth；
Step 43, setting the initial value of the RF peak value statistic vector sequence number i '' in each PW block to be 1, and setting the ith '' RF peak value statistic vector peak _ PW [ i '']Is 0; RF parameter comparison window fwindowIs RF as an initial valuebeginThe initial value of the RF peak breakpoint pulse number k' is 1;
step 44, setting the initial value of the pulse serial number j' selected from the staggered full pulse sequence in the corresponding PW block as 1;
step 45, determine if i' is greater than n _ RFwidthIf yes, go to step 47; if not, go to step 46;
step 46, when j'<N_PWpidAnd f iswindow<x3(j'')≤fwindow+fgateThen peak _ RF [ i']Adding 1 and then giving peak _ RF [ i']Adding 1 to j 'and assigning j', and returning to step 45; when j'>N_PWpidThen go to step 37; when x is3(j'')>fwindow+fgateAdding 1 to k ', adding 1 to j', and fwindowPlus fgateRear given fwindowAdding 1 to i 'and assigning it to i', and returning to step 45;
wherein x is3(j ') is the value of the RF parameter corresponding to the jth' pulse within the corresponding PW block;
and step 47, classifying the pulse data with the peak value of the RF parameter in the corresponding PW sub-block into a first pulse data block, and classifying the residual pulse data into a second pulse data block.
Further, in the fifth step, the specific process of sequentially performing RF-PRI projection sample map extraction and unknown radar target pulse matching processing on the first pulse data block includes:
step 511, setting the initial value of the full pulse shift digit i ' ' ' _ shift to 1, and setting the initial value of the pulse matching flag word corresponding to each pulse to zero;
step 512, setting the initial value of the pulse sequence k ' ″ in the first pulse data block to be 1, setting the initial value of the pulse sequence j ' ″ in the shifted first pulse data block to be i ' _ shift +1, and setting the initial value of the Number _ matched of the first matching pulses corresponding to the shift to be 0;
step 513, judge | x3(k''')-x3(j ' ' ') if the resolution f of the RF value is less than or equal togateIf yes, go to step 514; if not, add 1 to j '' 'and give it to j' '', go to step 515;
wherein x is3(k ' ' ') and x3(j '' ') the values of the RF parameters corresponding to the k' '' th pulse in the first pulse data block, respectively; x is the number of3(j '' ') is the RF parameter value corresponding to the jth' '' pulse in the shifted first pulse data block;
step 514, determine | [ x |)4(k''')-x4(1)]-[x4(j''')-x4(1+i'''_shift)]| whether or not the PRI value resolution T is less than or equal togateIf yes, then x is4(k ' ' ') the corresponding pulse matching flag is set to 1, and Number _ matched is added with 1 and then given to Number _ matched, k ' ' ' is added with 1 and then given to k ' ' ', j ' ' ' is added with 1 and then given to j ' ' ', and step 515 is entered; if not, add 1 to j '' 'and give it to j' '', go to step 515;
wherein x is4(k ' ' ') and x4(1) PRI parameter values corresponding to the k ' ' ' and 1 pulse, respectively, in the first pulse data block; x is the number of4(j ' ' ') and x4(1+ i '' '_ shift) is a PRI parameter value corresponding to the jth' '' 'and 1+ i' '' _ shift pulses in the shifted first pulse data block respectively; step 515, determine if j' ″ is less than N _ number, if yes, return to step 513; if not, go to step 516; wherein, N _ number is the total number of pulses in the first pulse data block;
step 516, judging whether the Number _ mapped is larger than the Number _ Threshold, if so, then the corresponding RF-PRI projection sample at the momentThe skeleton period of the picture is T _ frame ═ x4(1+i'''_shift)-x4(1) Step 517 is entered; if not, adding 1 to i '' '_ shift, and then giving the i' '' _ shift, and returning to the step 513;
wherein, Number _ Threshold is a first matching pulse Number Threshold value;
517, extracting the pulse with the pulse matching flag word set as 1, and determining the pulse with the minimum arrival time value as a first pulse z (1); step 518, a pulse matching search is performed with the first pulse z (1) as the start pulse and the corresponding RF-PRI projection sample map skeleton period T _ frame as the time delay.
Further, in the fifth step, the specific process of sequentially performing PRI projection sample image extraction and unknown radar target pulse matching processing on the second pulse data block is as follows:
step 521, setting an initial value of a full pulse shift digit i ' ' ' _ shift to be 1, and setting an initial value of a pulse matching flag word corresponding to each pulse to be zero;
step 522, setting the initial value of the pulse sequence k '' '' in the second pulse data block to 1, setting the initial value of the pulse sequence j '' '' in the shifted second pulse data block to i '' '_ shift +1, and setting the initial value of the Number _ matched' of the second matching pulses corresponding to the shift to 0;
step 523, judge | [ x4'(k'''')-x4'(1)]-[x4'(j'''')-x4'(1+i''''_shift)]| whether or not the PRI value resolution T is less than or equal togateIf yes, then x is4'(k' '') has a corresponding pulse matching reference character set to 1, and gives the Number _ matched 'after adding 1 to the Number _ matched', the k '' '' after adding 1 to the Number _ matched ', and the j' '' 'after adding 1 to the Number _ matched', the j '' '' goes to step 524; if not, add 1 to j '' ', and then give j' '', go to step 524;
wherein x is4'(k' '') and x4' (1) PRI parameter values corresponding to kth ' ' ' ' and 1 st pulses, respectively, in the second pulse data block; x is the number of4'(j' '') and x4' (1+ i ' ' ' _ shift) are the PRI parameter values corresponding to the jth ' ' ' ' and 1+ i ' ' ' _ shift pulses in the shifted second pulse data block, respectively;
step 524, judging whether j ' ' ' ' is smaller than N _ number ', if yes, returning to step 523; if not, go to step 525; wherein, N' _ number is the total number of pulses in the first pulse data block;
step 525, determine whether Number _ matched 'is greater than Number _ Threshold', if so, the corresponding skeleton period of the PRI projection sample map is T _ frame ═ x at this time4'(1+i''''_shift)-x4' (1), proceed to step 526; if not, add 1 to i ' ' ' ' _ shift and give it to i ' ' ' _ shift, return to step 523;
wherein, Number _ Threshold' is a second matching pulse Number Threshold value;
step 526, extracting the pulse with the pulse matching flag word set to be 1, and determining the pulse with the minimum arrival time value as the first pulse n (1);
and step 527, taking the first pulse n (1) as an initial pulse, and taking the corresponding PRI projection sample image skeleton period T _ frame' as a time delay to perform unknown radar target pulse matching search.
Further, in the third step, PW partitioning is performed on multiple DOA partitions by using multiple processes.
Further, in step four, multiple processes are adopted to simultaneously perform RF blocking processing on multiple PW blocks.
In order to achieve the second purpose, the invention adopts the following technical scheme:
a storage medium storing program instructions; the program instructions, when executed, implement the fast extraction method described above.
In order to achieve the third purpose, the invention adopts the following technical scheme:
a rapid extraction system of unknown radar target pulse sequences comprises the storage medium.
The invention has the beneficial effects that:
according to the invention, the staggered full-pulse sequence DOA block division and the PW block division are realized by the signal arrival azimuth DOA parameter peak value and the pulse width PW parameter peak value, so that outlier points (including DOA outlier points and PW block outlier points) and noise interference are eliminated; through carrier frequency RF parameter peak values, radar target pulses (RF fixed types, RF hopping types and RF pulse group agility types, namely first pulse data blocks) with regularly changed RF parameters in data in the blocks are separated from radar target pulses (RF inter-pulse agility types and second pulse data blocks) with randomly changed RF parameters, so that on one hand, the phenomenon that alarm leakage is caused due to the fact that target signals of the RF inter-pulse agility types are lost because a full pulse sequence with the peak values which are not counted is removed is avoided; dividing all the pulses corresponding to the RF parameter peak values into one block, and then carrying out subsequent processing, so that the target signals of RF hopping and RF pulse group agility types are separated into a plurality of data blocks, and the target data batching phenomenon is caused; and the first pulse data block is subjected to RF-PRI projection sample image extraction and unknown radar target pulse matching processing in sequence, and the second pulse data block is subjected to PRI projection sample image extraction and unknown radar target pulse matching processing in sequence, so that the operation time of irregularly-changed characteristic parameters during matching is reduced, the calculated amount is reduced, the processing time is saved, the processing speed is improved, the influence of randomly-changed characteristic parameters on the accuracy of a matching value during similarity measurement is avoided, and the extraction accuracy is improved.
Drawings
FIG. 1 is a schematic flow chart of a method for rapidly extracting an unknown radar target pulse sequence according to the present invention;
fig. 2 is a schematic flow chart of the process of sequentially performing RF-PRI projection sample map extraction and unknown radar target pulse matching on a first pulse data block.
Detailed Description
The following detailed description of embodiments of the invention refers to the accompanying drawings.
The embodiment provides a method for rapidly extracting an unknown radar target pulse sequence, and referring to fig. 1, the method for rapidly extracting the unknown radar target pulse sequence includes the following steps:
step one, acquiring a staggered full pulse sequence.
The radar pulse description word in this embodiment includes a signal azimuth of arrival (DOA) parameter, a carrier frequency (RF) parameter, a Pulse Width (PW) parameter, and a time of arrival (TOA) parameter. The full pulse sequence mathematical model Y is represented as:
Figure GDA0002676661470000061
wherein Y is a full pulse sequence mathematical model; y isnIs a characteristic parameter vector of the nth pulse, yn＝[x1(n),x2(n),x3(n),x4(n)]， x1(n)～x4(n) respectively representing the DOA parameter value, the PW parameter value, the RF parameter value and the TOA parameter value corresponding to the nth pulse; n is 1,2,3, and N is the number of pulses in the interleaved full pulse train.
And step two, acquiring a signal arrival azimuth DOA parameter peak value of each pulse data in the staggered full pulse sequence and performing DOA blocking processing.
For pulse blocking corresponding to the DOA parameter peak value, assuming that the number of the peak values is M, dividing the full pulse data into M blocks, judging the pulse which is not distributed to any data block as a outlier or an interference pulse and eliminating the outlier and a disordered pulse, reducing the operation amount, improving the operation efficiency, and specifically realizing the DOA blocking processing course as follows:
step 21, reading the pulse data in the interleaved full pulse sequence and sequencing according to the DOA parameter value to obtain the maximum value DOA of the DOA parameter valuemaxAnd minimum value DOAmin。
In this embodiment, the full-pulse data is sorted according to the order of the DOA parameter values from small to large, and the maximum value and the minimum value of the obtained DOA parameter values are respectively DOAmaxAnd DOAmin。
When the DOA parameter value is [ -alpha ]limit，αlimit]When the data is internal, the corresponding full pulse data is valid; otherwise, the corresponding full pulse data, alpha, is rejectedlimitThe maximum DOA parameter threshold is typically 30 deg..
Step 22, according to the preset DOA value resolution alphagateAnd maximum value of DOA parameter value DOAmaxAnd minimum value DOAminCalculating the DOA peak value statistical initial value DOAbeginDOA peak value statistic end value DOAendAnd DOA peak statisticsExtent of DOAwidthAnd DOA Peak statistical Window number n _ DOAwid。
In this step, the DOA peak value statistics initial value DOAbeginAnd a statistical end value DOAendRespectively as follows:
DOAbegin＝DOAmin-αgate-1(°)；
DOAend＝DOAmax+αgate+1(°)；
the statistical range of the DOA peak value is DOAwidth＝DOAend-DOAbegin；
DOA peak statistical window number of
Figure GDA0002676661470000071
Wherein, DOA value resolution agate＝0.5°，
Figure GDA0002676661470000072
Indicating rounding up.
Step 23, setting the initial value of the DOA peak value statistical vector serial number i as 1, and setting the ith DOA peak value statistical vector peak _ DOA [ i [ ]]Is 0; DOA parameter comparison window alphawindowIs initialized to DOAbeginThe initial value of the DOA peak value breakpoint pulse sequence number k is 1;
step 24, setting the initial value of the pulse serial number j selected from the staggered full pulse sequence as 1;
step 25, judging whether i is larger than n _ DOAwidthIf yes, go to step 27; if not, go to step 26;
step 26, when j is<N and alphawindow<x1(j)≤αwindow+αgateThen peak _ DOA [ i ]]Adding 1 and then assigning peak _ DOA [ i]J is added with 1 and then is given to j, and the step 25 is returned; when j is>N, then go to step 27; when x is1(j)>αwindow+αgateAdding 1 to k, adding 1 to j, and adding 1 to j, alphawindowPlus alphagateAfter giving alphawindowAdding 1 to i, assigning i, and returning to the step 25;
wherein x is1(j) The DOA parameter value corresponding to the jth pulse; n is the number of pulses in the interleaved full pulse sequence.
Step 27, counting the number M of pulses with the DOA parameter value being the wave crest in the DOA peak value statistical vector sequence; and the interleaved full pulse sequence is divided into M DOA blocks.
The peak conditions of this embodiment are: peak _ DOA [ i ] is greater than peak _ DOA [ i-1], and peak _ DOA [ i ] is greater than or equal to peak _ DOA [ i +1 ]; or when peak _ DOA [ i ] is equal to peak _ DOA [ i-1] and peak _ DOA [ i ] is greater than peak _ DOA [ i +1 ].
The conditions of the wave trough are: peak _ DOA [ i ] is less than peak _ DOA [ i-1], and peak _ DOA [ i ] is less than or equal to peak _ DOA [ i +1 ]; or peak _ DOA [ i ] is equal to peak _ DOA [ i-1], and peak _ DOA [ i ] is less than peak _ DOA [ i +1 ].
And step three, acquiring the pulse width PW parameter peak value of each pulse data in each DOA block and carrying out PW block processing.
Pulse blocking corresponding to a PW parameter peak value, wherein the number of the peak values is assumed to be L, full pulse data is divided into M blocks, pulse which is not distributed to any data block is judged to be outlier points or interference pulses and is eliminated, PW outlier points and chaotic pulses are eliminated, the operation amount is reduced, the operation efficiency is improved, and the specific implementation process of PW blocking processing is as follows:
step 31, obtaining the number of pulses N _ DOA in each DOA blockpid(ii) a Sequencing the pulse data in each DOA block according to the pulse width PW parameter value to obtain the maximum PW of the pulse PW parameter values in the corresponding DOA blockmaxAnd minimum value PWmin；
Step 32, according to the preset PW value resolution taugateAnd the maximum value PW of the pulse PW parameter values in each DOA blockmaxAnd minimum value PWminCalculating the PW peak value statistic initial value PW in the corresponding DOA blockbeginPW peak value statistical end value PWendAnd PW peak statistical range PWwidthAnd the number of PW peak statistical windows n _ PWwidth。
PWbegin＝PWmin-τgate；
PWend＝PWmax+τgate；
PWwidth＝PWend-PWbegin；
Figure GDA0002676661470000082
Wherein the PW value resolution τgate＝0.1μs，
Figure GDA0002676661470000081
Indicating rounding up.
Step 33, setting the initial value of the PW peak value statistical vector serial number i ' in each DOA block to be 1, and setting the ith ' PW peak value statistical vector peak _ PW [ i ']Is 0; PW parameter comparison window τwindowHas an initial value of PWbeginThe initial value of PW peak breakpoint pulse number k' is 1;
step 34, setting the initial value of the pulse sequence number j' selected from the staggered full pulse sequence in the corresponding DOA block to be 1;
step 35, judging whether i' is larger than n _ PWwidthIf yes, go to step 37; if not, go to step 36;
step 36, when j'<N_DOApidAnd τwindow<x2(j')≤τwindow+τgateThen peak _ PW [ i']Adding 1 and then giving peak _ PW [ i']Adding 1 to j ', and then assigning j', and returning to the step 35; when j'>N_DOApidThen go to step 37; when x is2(j')>τwindow+ τgateAdding 1 to k ', j ', and then adding 1 to j ', tauwindowPlus taugatePost-assigning to tauwindowAdding 1 to i ', assigning i', and returning to the step 35; wherein x is2(j') is a PW parameter value corresponding to the jth pulse in the corresponding DOA block;
step 37, counting the number L of pulses with wave crest as the PW parameter value in the PW peak value statistical vector sequence in each DOA block; and dividing the corresponding DOA block-wise intra-interleaved full pulse sequence into L PW blocks.
The peak conditions of this embodiment are: peak _ PW [ i '] is greater than peak _ PW [ i' -1], and peak _ PW [ i '] is greater than or equal to peak _ PW [ i' +1 ]; or when peak _ PW [ i '] is equal to peak _ PW [ i' -1], and peak _ PW [ i '] is greater than peak _ PW [ i' +1 ].
The conditions of the wave trough are: peak _ PW [ i '] is less than peak _ PW [ i' -1], and peak _ PW [ i '] is less than or equal to peak _ PW [ i' +1 ]; or peak _ PW [ i '] is equal to peak _ PW [ i' -1], and peak _ PW [ i '] is less than peak _ PW [ i' +1 ].
And step four, obtaining a carrier frequency RF parameter peak value of each pulse data in each PW block and carrying out RF block processing to obtain a first pulse data block corresponding to the RF parameter peak value and a second pulse data block corresponding to the non-RF parameter peak value.
According to the method, through RF parameter peak value statistics, radar target pulses (RF fixed types, RF jumping types and RF pulse group agility types) with regular RF parameter changes in data in a block are separated from radar target pulses (RF interpulse agility types) with random RF parameter changes, so that firstly, a full pulse sequence without the counted peak value can be prevented from being removed, and therefore target signals of the RF interpulse agility types are lost, and the phenomenon of alarm missing is caused; and secondly, the phenomenon that target data are batched because target signals of RF hopping and RF pulse group agility types are separated into a plurality of data blocks due to the fact that each pulse corresponding to the peak value is divided into one block and then subsequent processing is carried out is avoided. The specific implementation process of the RF partitioning process is as follows:
step 41, obtaining the number N _ PW of pulses in each PW blockpid(ii) a Sequencing the pulse data in each PW block according to the pulse width RF parameter value to obtain the maximum value RF of the pulse RF parameter value in the corresponding PW blockmaxAnd minimum value RFmin；
Step 42, resolution f according to the preset RF valuegateAnd maximum value RF of the pulse RF parameter value within each PW blockmaxAnd minimum value RFminCalculating the statistical initial value RF of the RF peak value in the corresponding PW blockbeginStatistical end of RF peak value RFendAnd RF peak statistical range RFwidthAnd number of RF peak statistics windows n _ RFwidth；
Resolution f of RF value of the present embodimentgate3MHz, the resolution of RF parameter peak statistics, if the difference between the RF values of two pulses is less than fgateThen, it is determined that the RF values are the same。
Step 43, setting the initial value of the RF peak value statistic vector sequence number i '' in each PW block to be 1, and setting the ith '' RF peak value statistic vector peak _ PW [ i '']Is 0; RF parameter comparison window fwindowIs RF as an initial valuebeginThe initial value of the RF peak breakpoint pulse number k' is 1;
step 44, setting the initial value of the pulse sequence number j' selected from the interleaved full pulse sequence in the corresponding PW block to 1;
step 45, determine if i' is greater than n _ RFwidthIf yes, go to step 47; if not, go to step 46;
step 46, when j'<N_PWpidAnd f iswindow<x3(j'')≤fwindow+fgateThen peak _ RF [ i']Adding 1 and then giving peak _ RF [ i']Adding 1 to j 'and assigning j', and returning to step 45; when j'>N_PWpidThen go to step 37; when x is3(j'')>fwindow+fgateAdding 1 to k ', adding 1 to j', and fwindowPlus fgateRear given fwindowAdding 1 to i 'and assigning it to i', and returning to step 45; wherein x is3(j ') is the value of the RF parameter corresponding to the jth' pulse within the corresponding PW block;
and step 47, classifying the pulse data with the peak value of the RF parameter in the corresponding PW sub-block into a first pulse data block, and classifying the residual pulse data into a second pulse data block.
In order to save time and cost and increase processing speed, in this embodiment, the number of sub-processes is L according to the PW block number L, which is respectively numbered as sub-process 1, sub-process 2, and sub-process …, and a plurality of sub-processes behind the sub-process L perform RF block processing synchronously.
Step five, judging whether the absolute value of the difference between the RF parameter value of the unknown radar target pulse and the RF parameter peak value is smaller than a threshold value, if so, sequentially performing RF-PRI projection sample image extraction and unknown radar target pulse matching processing on the first pulse data block; and if not, sequentially carrying out PRI projection sample image extraction and unknown radar target pulse matching processing on the second pulse data block.
Referring to fig. 2, the specific process of sequentially performing RF-PRI projection sample map extraction and unknown radar target pulse matching processing on the first pulse data block in this embodiment is as follows:
step 511, setting the initial value of the full pulse shift digit i ' ' ' _ shift to 1, and setting the initial value of the pulse matching flag word corresponding to each pulse to zero;
step 512, setting the initial value of the pulse sequence k ' ″ in the first pulse data block to be 1, setting the initial value of the pulse sequence j ' ″ in the shifted first pulse data block to be i ' _ shift +1, and setting the initial value of the Number _ matched of the first matching pulses corresponding to the shift to be 0;
step 513, judge | x3(k''')-x3(j ' ' ') if the resolution f of the RF value is less than or equal togateIf yes, go to step 514; if not, add 1 to j '' 'and give it to j' '', go to step 515;
wherein x is3(k ' ' ') and x3(j '' ') the values of the RF parameters corresponding to the k' '' th pulse in the first pulse data block, respectively; x is the number of3(j '' ') is the RF parameter value corresponding to the jth' '' pulse in the shifted first pulse data block;
step 514, determine | [ x |)4(k''')-x4(1)]-[x4(j''')-x4(1+i'''_shift)]| whether or not the PRI value resolution T is less than or equal togateIf yes, then x is4(k ' ' ') the corresponding pulse matching flag is set to 1, and Number _ matched is added with 1 and then given to Number _ matched, k ' ' ' is added with 1 and then given to k ' ' ', j ' ' ' is added with 1 and then given to j ' ' ', and step 515 is entered; if not, add 1 to j '' 'and give it to j' '', go to step 515;
wherein x is4(k ' ' ') and x4(1) PRI parameter values corresponding to the k ' ' ' and 1 pulse, respectively, in the first pulse data block; x is the number of4(j ' ' ') and x4(1+ i '' '_ shift) is a PRI parameter value corresponding to the jth' '' 'and 1+ i' '' _ shift pulses in the shifted first pulse data block respectively; PRI value resolution T of the present embodiment gate1 μ s. Wherein, the radar target maximum possible and minimum possible PRI values are respectively Tlimit_Max20000 μ s and tlimit_Min＝5μs。
Step 515, determine if j' ″ is less than N _ number, if yes, return to step 513; if not, go to step 516;
wherein, N _ number is the total number of pulses in the first pulse data block;
step 516, determining whether Number _ mapped is greater than Number _ Threshold, if yes, the corresponding RF-PRI projection sample graph skeleton period is T _ frame ═ x at this time4(1+i'''_shift)-x4(1) Step 517 is entered; if not, adding 1 to i '' '_ shift, and then giving the i' '' _ shift, and returning to the step 513;
wherein, Number _ Threshold is a first matching pulse Number Threshold value;
517, extracting the pulse with the pulse matching flag word set as 1, and determining the pulse with the minimum arrival time value as a first pulse z (1);
step 518, a pulse matching search is performed with the first pulse z (1) as the start pulse and the corresponding RF-PRI projection sample map skeleton period T _ frame as the time delay.
In this embodiment, the specific process of sequentially performing PRI projection sample image extraction and unknown radar target pulse matching processing on the second pulse data block is as follows:
step 521, setting an initial value of a full pulse shift digit i ' ' ' _ shift to be 1, and setting an initial value of a pulse matching flag word corresponding to each pulse to be zero;
step 522, setting the initial value of the pulse sequence k '' '' in the second pulse data block to 1, setting the initial value of the pulse sequence j '' '' in the shifted second pulse data block to i '' '_ shift +1, and setting the initial value of the Number _ matched' of the second matching pulses corresponding to the shift to 0;
step 523, judge | [ x4'(k'''')-x4'(1)]-[x4'(j'''')-x4'(1+i''''_shift)]| whether or not the PRI value resolution T is less than or equal togateIf yes, then x is4'(k' '') has a corresponding pulse matching reference character set to 1, and gives the Number _ matched 'after adding 1 to the Number _ matched', the k '' '' after adding 1 to the Number _ matched ', and the j' '' 'after adding 1 to the Number _ matched', the j '' '' goes to step 524; if not, add 1 to j '' ', and then give j' '', go to step 524;
wherein x is4'(k' '') and x4' (1) PRI parameter values corresponding to kth ' ' ' ' and 1 pulse, respectively, in the second pulse data block; x is the number of4'(j' '') and x4' (1+ i ' ' ' _ shift) is a PRI parameter value corresponding to jth ' ' ' ' and 1+ i ' ' ' _ shift pulses in the shifted second pulse data block, respectively;
step 524, judging whether j ' ' ' ' is smaller than N _ number ', if yes, returning to step 523; if not, go to step 525; wherein, N' _ number is the total number of pulses in the first pulse data block;
step 525, determine whether Number _ matched 'is greater than Number _ Threshold', if so, the corresponding skeleton period of the PRI projection sample map is T _ frame ═ x at this time4'(1+i''''_shift)-x4' (1), proceed to step 526; if not, add 1 to i ' ' ' ' _ shift and give it to i ' ' ' _ shift, return to step 523;
wherein, Number _ Threshold' is a second matching pulse Number Threshold value;
step 526, extracting the pulse with the pulse matching flag word set to be 1, and determining the pulse with the minimum arrival time value as the first pulse n (1);
and step 527, taking the first pulse n (1) as an initial pulse, and taking the corresponding PRI projection sample image skeleton period T _ frame' as a time delay to perform unknown radar target pulse matching search.
In the embodiment, the staggered full-pulse sequence DOA blocking and PW blocking are realized by the signal arrival azimuth DOA parameter peak value and the pulse width PW parameter peak value, so that outlier points (including DOA outlier points and PW blocked outlier points) and noise interference are eliminated; through carrier frequency RF parameter peak values, radar target pulses (RF fixed types, RF hopping types and RF pulse group agility types, namely first pulse data blocks) with regularly changed RF parameters in data in the blocks are separated from radar target pulses (RF inter-pulse agility types and second pulse data blocks) with randomly changed RF parameters, so that on one hand, the phenomenon that alarm leakage is caused due to the fact that target signals of the RF inter-pulse agility types are lost because a full pulse sequence with the peak values which are not counted is removed is avoided; dividing all the pulses corresponding to the RF parameter peak values into one block, and then carrying out subsequent processing, so that the target signals of RF hopping and RF pulse group agility types are separated into a plurality of data blocks, and the target data batching phenomenon is caused; and the first pulse data block is subjected to RF-PRI projection sample image extraction and unknown radar target pulse matching processing in sequence, and the second pulse data block is subjected to PRI projection sample image extraction and unknown radar target pulse matching processing in sequence, so that the operation time of irregularly-changed characteristic parameters during matching is reduced, the calculated amount is reduced, the processing time is saved, the processing speed is improved, the influence of randomly-changed characteristic parameters on the accuracy of a matching value during similarity measurement is avoided, and the extraction accuracy is improved.
Another embodiment provides a storage medium storing program instructions; the program instructions, when executed, implement the fast extraction method provided by the above embodiments.
Yet another embodiment provides a fast extraction system of unknown radar target pulse sequences, the fast extraction system comprising the storage medium of the above embodiment.
Although the embodiments of the present invention have been described in detail with reference to the preferred embodiments, it should be understood by those skilled in the art that various changes may be made and equivalents may be substituted without departing from the spirit and scope of the embodiments of the present invention.

Claims (10)

1. A rapid extraction method of an unknown radar target pulse sequence is characterized by comprising the following steps:
step one, acquiring a staggered full pulse sequence;
step two, acquiring a signal arrival azimuth DOA parameter peak value of each pulse data in the staggered full pulse sequence and performing DOA blocking processing;
step three, acquiring a pulse width PW parameter peak value of each pulse data in each DOA block and carrying out PW block processing;
step four, obtaining a carrier frequency RF parameter peak value of each pulse data in each PW block and carrying out RF block processing to obtain a first pulse data block corresponding to the RF parameter peak value and a second pulse data block corresponding to a non-RF parameter peak value;
step five, judging whether the absolute value of the difference between the RF parameter value of the unknown radar target pulse and the RF parameter peak value is smaller than a threshold value, if so, sequentially performing RF-PRI projection sample image extraction and unknown radar target pulse matching processing on the first pulse data block; and if not, sequentially carrying out PRI projection sample image extraction and unknown radar target pulse matching processing on the second pulse data block.
2. The rapid extraction method according to claim 1, wherein the second step is implemented by the following steps:
step 21, reading the pulse data in the interleaved full pulse sequence and sequencing according to the DOA parameter value to obtain the maximum value DOA of the DOA parameter valuemaxAnd minimum value DOAmin；
Step 22, according to the preset DOA value resolution alphagateAnd maximum value of DOA parameter value DOAmaxAnd minimum value DOAminCalculating the DOA peak value statistical initial value DOAbeginDOA peak value statistic end value DOAendAnd DOA peak statistical range DOAwidthAnd DOA Peak statistical Window number n _ DOAwidth；
Step 23, setting the initial value of the DOA peak value statistical vector serial number i as 1, and setting the ith DOA peak value statistical vector peak _ DOA [ i [ ]]Is 0; DOA parameter comparison window alphawindowIs initialized to DOAbeginThe initial value of the DOA peak value breakpoint pulse sequence number k is 1;
step 24, setting the initial value of the pulse serial number j selected from the staggered full pulse sequence as 1;
step 25, judging whether i is larger than n _ DOAwidthIf yes, go to step 27; if not, go to step 26;
step 26, when j is<N and alphawindow<x1(j)≤αwindow+αgateThen peak _ DOA [ i ]]Adding 1 and then assigning peak _ DOA [ i]J is added with 1 and then is given to j, and the step 25 is returned; when j is>N, thenEntering step 27; when x is1(j)>αwindow+αgateAdding 1 to k, adding 1 to j, and adding 1 to j, alphawindowPlus alphagateAfter giving alphawindowAdding 1 to i, assigning i, and returning to the step 25;
wherein x is1(j) The DOA parameter value corresponding to the jth pulse; n is the number of pulses in the staggered full pulse sequence;
step 27, counting the number M of pulses with the DOA parameter value being the wave crest in the DOA peak value statistical vector sequence; and the interleaved full pulse sequence is divided into M DOA blocks.
3. The rapid extraction method according to claim 2, wherein the third step is realized by the following steps:
step 31, obtaining the number of pulses N _ DOA in each DOA blockpid(ii) a Sequencing the pulse data in each DOA block according to the pulse width PW parameter value to obtain the maximum PW of the pulse PW parameter values in the corresponding DOA blockmaxAnd minimum value PWmin；
Step 32, according to the preset PW value resolution taugateAnd the maximum value PW of the pulse PW parameter values in each DOA blockmaxAnd minimum value PWminCalculating the PW peak value statistic initial value PW in the corresponding DOA blockbeginPW peak value statistical end value PWendAnd PW peak statistical range PWwidthAnd the number of PW peak statistical windows n _ PWwidth；
Step 33, setting the initial value of the PW peak value statistical vector serial number i ' in each DOA block to be 1, and setting the ith ' PW peak value statistical vector peak _ PW [ i ']Is 0; PW parameter comparison window τwindowHas an initial value of PWbeginThe initial value of PW peak breakpoint pulse number k' is 1;
step 34, setting the initial value of the pulse serial number j' selected from the staggered full pulse sequence in the corresponding DOA block as 1;
step 35, judging whether i' is larger than n _ PWwidthIf yes, go to step 37; if not, go to step 36;
step 36, when j'<N_DOApidAnd τwindow<x2(j')≤τwindow+τgateThen peak _ PW [ i']Adding 1 and then giving peak _ PW [ i']Adding 1 to j ', and then assigning j', and returning to the step 35; when j'>N_DOApidThen go to step 37; when x is2(j')>τwindow+τgateAdding 1 to k ', j ', and then adding 1 to j ', tauwindowPlus taugatePost-assigning to tauwindowAdding 1 to i ', assigning i', and returning to the step 35;
wherein x is2(j') is a PW parameter value corresponding to the jth pulse in the corresponding DOA block;
step 37, counting the number L of pulses with wave crest as the PW parameter value in the PW peak value statistical vector sequence in each DOA block; and dividing the corresponding DOA block-wise intra-interleaved full pulse sequence into L PW blocks.
4. The rapid extraction method according to claim 1, wherein the specific implementation process of the step four is as follows:
step 41, obtaining the number N _ PW of pulses in each PW blockpid(ii) a Sequencing the pulse data in each PW block according to the pulse width RF parameter value to obtain the maximum value RF of the pulse RF parameter value in the corresponding PW blockmaxAnd minimum value RFmin；
Step 42, resolution f according to the preset RF valuegateAnd maximum value RF of the pulse RF parameter value within each PW blockmaxAnd minimum value RFminCalculating the statistical initial value RF of the RF peak value in the corresponding PW blockbeginStatistical end of RF peak value RFendAnd RF peak statistical range RFwidthAnd number of RF peak statistics windows n _ RFwidth；
Step 43, setting the initial value of the RF peak value statistic vector sequence number i '' in each PW block to be 1, and setting the ith '' RF peak value statistic vector peak _ PW [ i '']Is 0; RF parameter comparison window fwindowIs RF as an initial valuebeginThe initial value of the RF peak breakpoint pulse number k' is 1;
step 44, setting the initial value of the pulse serial number j' selected from the staggered full pulse sequence in the corresponding PW block as 1;
step 45, determine if i' is greater than n _ RFwidthIf yes, go to step 47; if not, go to step 46;
step 46, when j'<N_PWpidAnd f iswindow<x3(j'')≤fwindow+fgateThen peak _ RF [ i']Adding 1 and then giving peak _ RF [ i']Adding 1 to j 'and assigning j', and returning to step 45; when j'>N_PWpidThen go to step 37; when x is3(j'')>fwindow+fgateAdding 1 to k ', adding 1 to j', and fwindowPlus fgateRear given fwindowAdding 1 to i 'and assigning it to i', and returning to step 45;
wherein x is3(j ') is the value of the RF parameter corresponding to the jth' pulse within the corresponding PW block;
and step 47, classifying the pulse data with the peak value of the RF parameter in the corresponding PW sub-block into a first pulse data block, and classifying the residual pulse data into a second pulse data block.
5. The fast extraction method according to any one of claims 1 to 4, wherein in the fifth step, the specific process of sequentially performing RF-PRI projection sample map extraction and unknown radar target pulse matching processing on the first pulse data block is as follows:
step 511, setting the initial value of the full pulse shift digit i ' ' ' _ shift to 1, and setting the initial value of the pulse matching flag word corresponding to each pulse to zero;
step 512, setting the initial value of the pulse sequence k ' ″ in the first pulse data block to be 1, setting the initial value of the pulse sequence j ' ″ in the shifted first pulse data block to be i ' _ shift +1, and setting the initial value of the Number _ matched of the first matching pulses corresponding to the shift to be 0;
step 513, judge | x3(k''')-x3(j ' ' ') if the resolution f of the RF value is less than or equal togateIf yes, go to step 514; if not, add 1 to j '' 'and give it to j' '', go to step 515;
wherein x is3(k ' ' ') is the kth pulse in the first pulse data block'' RF parameter values corresponding to the pulses; x is the number of3(j '' ') is the RF parameter value corresponding to the jth' '' pulse in the shifted first pulse data block;
step 514, determine | [ x |)4(k''')-x4(1)]-[x4(j''')-x4(1+i'''_shift)]| whether or not the PRI value resolution T is less than or equal togateIf yes, then x is4(k ' ' ') the corresponding pulse matching flag is set to 1, and Number _ matched is added with 1 and then given to Number _ matched, k ' ' ' is added with 1 and then given to k ' ' ', j ' ' ' is added with 1 and then given to j ' ' ', and step 515 is entered; if not, add 1 to j '' 'and give it to j' '', go to step 515;
wherein x is4(k ' ' ') and x4(1) PRI parameter values corresponding to the k ' ' ' and 1 pulse, respectively, in the first pulse data block; x is the number of4(j ' ' ') and x4(1+ i '' '_ shift) is a PRI parameter value corresponding to the jth' '' 'and 1+ i' '' _ shift pulses in the shifted first pulse data block respectively;
step 515, determine if j' ″ is less than N _ number, if yes, return to step 513; if not, go to step 516;
wherein, N _ number is the total number of pulses in the first pulse data block;
step 516, determining whether Number _ mapped is greater than Number _ Threshold, if yes, the corresponding RF-PRI projection sample graph skeleton period is T _ frame ═ x at this time4(1+i'''_shift)-x4(1) Step 517 is entered; if not, adding 1 to i '' '_ shift, and then giving the i' '' _ shift, and returning to the step 513;
wherein, Number _ Threshold is a first matching pulse Number Threshold value;
517, extracting the pulse with the pulse matching flag word set as 1, and determining the pulse with the minimum arrival time value as a first pulse z (1);
step 518, a pulse matching search is performed with the first pulse z (1) as the start pulse and the corresponding RF-PRI projection sample map skeleton period T _ frame as the time delay.
6. The fast extraction method according to any one of claims 1 to 4, wherein in the fifth step, the specific process of sequentially performing PRI projection sample map extraction and unknown radar target pulse matching processing on the second pulse data block is as follows:
step 521, setting an initial value of a full pulse shift digit i ' ' ' _ shift to be 1, and setting an initial value of a pulse matching flag word corresponding to each pulse to be zero;
step 522, setting the initial value of the pulse sequence k '' '' in the second pulse data block to 1, setting the initial value of the pulse sequence j '' '' in the shifted second pulse data block to i '' '_ shift +1, and setting the initial value of the Number _ matched' of the second matching pulses corresponding to the shift to 0;
step 523, judge | [ x4'(k'''')-x4'(1)]-[x4'(j'''')-x4'(1+i''''_shift)]| whether or not the PRI value resolution T is less than or equal togateIf yes, then x is4'(k' '') has a corresponding pulse matching reference character set to 1, and gives the Number _ matched 'after adding 1 to the Number _ matched', the k '' '' after adding 1 to the Number _ matched ', and the j' '' 'after adding 1 to the Number _ matched', the j '' '' goes to step 524; if not, add 1 to j '' ', and then give j' '', go to step 524;
wherein x is4'(k' '') and x4' (1) PRI parameter values corresponding to kth ' ' ' ' and 1 st pulses, respectively, in the second pulse data block; x is the number of4'(j' '') and x4' (1+ i ' ' ' _ shift) are the PRI parameter values corresponding to the jth ' ' ' ' and 1+ i ' ' ' _ shift pulses in the shifted second pulse data block, respectively;
step 524, judging whether j ' ' ' ' is less than N ' _ number, if yes, returning to step 523; if not, go to step 525; wherein, N' _ number is the total number of pulses in the second pulse data block;
step 525, determine whether Number _ matched 'is greater than Number _ Threshold', if so, the corresponding skeleton period of the PRI projection sample map is T _ frame ═ x at this time4'(1+i''''_shift)-x4' (1), proceed to step 526; if not, add 1 to i ' ' ' ' _ shift and give it to i ' ' ' _ shift, return to step 523;
wherein, Number _ Threshold' is a second matching pulse Number Threshold value;
step 526, extracting the pulse with the pulse matching flag word set to be 1, and determining the pulse with the minimum arrival time value as the first pulse n (1);
and step 527, taking the first pulse n (1) as an initial pulse, and taking the corresponding PRI projection sample image skeleton period T _ frame' as a time delay to perform unknown radar target pulse matching search.
7. The fast extraction method according to any one of claims 1 to 4, characterized in that in step three, PW block processing is performed on a plurality of DOA blocks simultaneously by using multiple processes.
8. The fast extraction method according to any one of claims 1 to 4, wherein in step four, multiple processes are used to simultaneously perform RF blocking on multiple PW blocks.
9. A storage medium, wherein the storage medium stores program instructions; the program instructions, when executed, implement the fast extraction method of any one of claims 1 to 8.
10. A fast extraction system of an unknown radar target pulse sequence, characterized in that it comprises a storage medium according to claim 9.
CN202010747915.4A 2020-07-30 2020-07-30 Method and system for rapidly extracting unknown radar target pulse sequence Active CN111624575B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010747915.4A CN111624575B (en) 2020-07-30 2020-07-30 Method and system for rapidly extracting unknown radar target pulse sequence

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010747915.4A CN111624575B (en) 2020-07-30 2020-07-30 Method and system for rapidly extracting unknown radar target pulse sequence

Publications (2)

Publication Number Publication Date
CN111624575A CN111624575A (en) 2020-09-04
CN111624575B true CN111624575B (en) 2020-10-23

Family

ID=72273126

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010747915.4A Active CN111624575B (en) 2020-07-30 2020-07-30 Method and system for rapidly extracting unknown radar target pulse sequence

Country Status (1)

Country Link
CN (1) CN111624575B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112363118B (en) * 2021-01-12 2021-03-26 中国人民解放军国防科技大学 Rapid high-precision batching method and system for satellite reconnaissance radar signals

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109270497A (en) * 2018-10-28 2019-01-25 西南电子技术研究所（中国电子科技集团公司第十研究所） The multi-Dimensional parameters Pre-sorting method of radar pulse signal
CN109613486A (en) * 2018-12-03 2019-04-12 中国人民解放军空军工程大学 A kind of Radar Signal Sorting Method based on core cluster support vector clustering
CN109683143A (en) * 2019-03-07 2019-04-26 西安电子科技大学 Radar Signal Sorting Method, device, computer equipment and can storage medium
US10345432B1 (en) * 2016-03-29 2019-07-09 Lockheed Martin Corporation Flexible wideband radar receiver architecture with enhanced sensitivity, dynamic range, instantaneous bandwidth, and range window capability
CN110376553A (en) * 2019-01-17 2019-10-25 西安电子科技大学 Low SNR signal parameter extracting method based on pulse cutting
CN110806563A (en) * 2019-11-19 2020-02-18 西南交通大学 Radiation source signal clustering and sorting method based on radar pulse aliasing degree judgment

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10345432B1 (en) * 2016-03-29 2019-07-09 Lockheed Martin Corporation Flexible wideband radar receiver architecture with enhanced sensitivity, dynamic range, instantaneous bandwidth, and range window capability
CN109270497A (en) * 2018-10-28 2019-01-25 西南电子技术研究所（中国电子科技集团公司第十研究所） The multi-Dimensional parameters Pre-sorting method of radar pulse signal
CN109613486A (en) * 2018-12-03 2019-04-12 中国人民解放军空军工程大学 A kind of Radar Signal Sorting Method based on core cluster support vector clustering
CN110376553A (en) * 2019-01-17 2019-10-25 西安电子科技大学 Low SNR signal parameter extracting method based on pulse cutting
CN109683143A (en) * 2019-03-07 2019-04-26 西安电子科技大学 Radar Signal Sorting Method, device, computer equipment and can storage medium
CN110806563A (en) * 2019-11-19 2020-02-18 西南交通大学 Radiation source signal clustering and sorting method based on radar pulse aliasing degree judgment

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于动态合批处理的重点雷达信号快速识别;朱永成 等;《电子信息对抗技术》;20180131;第33卷(第1期);第23-28页 *
基于合成特征参数和脉冲样本图的雷达信号分选;孟祥豪 等;《火力与指挥控制》;20120930;第37卷(第9期);第50-54页 *

Also Published As

Publication number Publication date
CN111624575A (en) 2020-09-04

Similar Documents

Publication Publication Date Title
CN113721219B (en) Radar signal sorting method and system based on multi-parameter clustering
Cain et al. Convolutional neural networks for radar emitter classification
CN112986928A (en) Signal sorting multi-source fusion processing method in complex electromagnetic environment
CN111624575B (en) Method and system for rapidly extracting unknown radar target pulse sequence
CN109613483B (en) Multi-target track initiation method based on Hough transformation
CN112906737B (en) Method for clustering and identifying based on density features based on multiple radiation sources
Xi et al. An improved SDIF algorithm for anti-radiation radar using dynamic sequence search
CN111487592B (en) Radar signal sorting method based on GPU and CPU platform drift mean clustering
CN108535704B (en) Signal pre-sorting method based on self-adaptive two-dimensional clustering
Shapero et al. Identifying agile waveforms with neural networks
CN113608193A (en) Radar multi-target distance and speed estimation method based on UNet
CN110531362A (en) A kind of object detection method of high-resolution moving sonar Knowledge-based
CN115629365A (en) Method and device for sorting few sample signals in high-density background signals
CN116299195A (en) Radar signal processing method based on TOA sequence relativity
CN115980689A (en) Point cloud detection-based radiation source signal sorting method, device, equipment and medium
CN111257839B (en) Radar signal sorting method
CN112363118B (en) Rapid high-precision batching method and system for satellite reconnaissance radar signals
CN113075620A (en) Signal sorting method based on multi-station time difference grid clustering
Behar et al. Hough detector with adaptive non-coherent integration for target detection in pulse jamming
CN111811515A (en) Multi-target track extraction method based on Gaussian mixture probability hypothesis density filter
CN110880012A (en) Correlation method for frequency information of agile radar radiation source between multiple scout platforms
CN110780275B (en) Method for removing batch increase by signal sorting
Davey et al. A markov model for initiating tracks with the probabilistic multi-hypothesis tracker
CN113702919B (en) Method and device for estimating PRI value and extracting pulse sequence
CN117310636B (en) Fixed pulse repetition interval measurement method, device and medium

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