CN114578293B - Electric scanning radar signal identification method by using intercepted signal amplitude value - Google Patents
Electric scanning radar signal identification method by using intercepted signal amplitude value Download PDFInfo
- Publication number
- CN114578293B CN114578293B CN202210183193.3A CN202210183193A CN114578293B CN 114578293 B CN114578293 B CN 114578293B CN 202210183193 A CN202210183193 A CN 202210183193A CN 114578293 B CN114578293 B CN 114578293B
- Authority
- CN
- China
- Prior art keywords
- amplitude
- value
- signal
- sequence
- main beam
- 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
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/28—Details of pulse systems
- G01S7/285—Receivers
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Abstract
The invention provides a method for identifying an electric scanning radar signal by using an intercepted signal amplitude value, which comprises the following steps: acquiring a pulse signal amplitude-arrival time sequence of radar target radiation intercepted by an electronic reconnaissance receiver, and sequentially performing unit conversion and normalization processing on each pulse signal amplitude; obtaining a maximum main lobe sequence; obtaining a self-correlation function value sequence; determining a first minimum value point on the left side and a first minimum value point on the right side of the autocorrelation function value of 0 shift digit in the autocorrelation function value sequence and calculating a first time length and a second time length; determining a first main beam signal and a second main beam signal; selecting the maximum main beam signal as a final main beam signal; and acquiring the flatness ratio of the final main beam signal and identifying the electric scanning radar signal. The method can realize the rapid identification of the radar target of an enemy in an electronic scanning mode and the rapid early warning of the radar target with high threat degree in the electronic reconnaissance process.
Description
Technical Field
The invention belongs to the technical field of radar signal identification, and particularly relates to an electric scanning radar signal identification method by utilizing an intercepted signal amplitude value.
Background
In military operations, enemy radars adopting an electric scanning antenna scanning mode are key targets of electronic countermeasures and reconnaissance, and phased array radars adopting an electric scanning antenna scanning mode are often adopted by enemy new system phased array radars. In the signal data intercepted by the reconnaissance receiver, the electric scanning radar signals are identified, and further the next signal parameter analysis is carried out, so that precious time can be won for battlefield preparation and electronic battle actions.
The radar signal adopting the electric scanning antenna scanning mode has the characteristics different from mechanical scanning radar signals, is mainly reflected in the amplitude value change characteristics of the intercepted signals, judges whether the intercepted radar signals are radiated by the electric scanning radar or not by extracting the amplitude value of the radar signal intercepted by the reconnaissance receiver and analyzing and processing the change characteristics of the amplitude value, can quickly judge whether the radar is a high-threat target or not, and formulates corresponding electronic warfare measures. The method utilizes the amplitude value to identify the electric scanning radar signal, and has the advantages of high identification speed, high identification precision, insensitivity to noise and the like compared with other parameters of the intercepted radar signal.
The method is characterized in that characteristic parameters corresponding to amplitude values of intercepted signals are extracted to judge the scanning types of the radar antenna, the prior art can extract characteristic parameters such as scanning period, kurtosis, the number of main lobes, amplitude variation of main lobe peak values of a pulse group, main lobe interval variation, main lobe and side lobe gain ratio and the like by utilizing pulse amplitude value data of the intercepted radar signals, and the judgment of the antenna scanning types corresponding to the mechanical scanning radar signals and the electric scanning radar signals is realized by utilizing the difference of different antenna scanning types on the characteristic parameters and adopting a decision tree algorithm of a support vector machine. The prior art implementation scheme has the following steps:
the first step is as follows: and inputting a data processing process.
Firstly, normalizing an intercepted signal amplitude sequence; then, extracting a maximum main lobe sequence by using the normalized amplitude sequence; and extracting the maximum value of the first-order difference absolute value aiming at the amplitude sequence corresponding to the maximum main lobe, calculating the ratio of the pulse number of the same wave bit, and taking the maximum value and the ratio as characteristic parameters for distinguishing a subsequent electronic scanning mode and a subsequent mechanical scanning mode.
The second step is that: and (4) identification processing of an electronic scanning mode.
Firstly, for an amplitude sequence in which the weighted sum of the pulse number ratio of the same wave bit and the maximum value of the first-order difference absolute value is greater than a set threshold value, a corresponding signal is judged as an electronic scanning radar signal. Then, extracting the mean square error of the amplitude sequence, taking the mean square error as a characteristic parameter to distinguish one-dimensional electronic scanning from two-dimensional electronic scanning, and judging as a two-dimensional electronic scanning mode when the mean square error is smaller than a set threshold value; otherwise, the method is determined as a one-dimensional electronic scanning mode.
The third step: and (4) identifying a mechanical scanning mode.
Firstly, for an amplitude sequence in which the weighted sum of the pulse number ratio of the same wave bit and the maximum value of the first-order difference absolute value is smaller than a set threshold value, a corresponding signal is judged as a mechanical scanning radar signal. Then, carrying out uniform sampling rate, antenna scanning period estimation and characteristic parameter extraction processing on the amplitude sequence to obtain kurtosis, the number of main lobes, the maximum difference value of the main lobe amplitudes and the maximum ratio of the main lobe interval of the amplitude sequence, taking the parameters as input data of a decision tree classifier of a support vector machine, distinguishing mechanical scanning modes of different change types, and identifying the mechanical scanning modes of antennas such as circular scanning, sector scanning, spiral scanning, raster scanning, conical scanning and the like.
In the prior art, the ratio of a first-order difference absolute value to the pulse number of the same wave position is adopted to judge whether the signal is an electric scanning radar signal, when the number of captured signal pulses is large, the calculated amount of two characteristic parameters is large, and the timeliness of the prior art is poor. In the prior art, processing processes such as electronic scanning mode identification processing, mechanical scanning mode identification processing and the like need to be carried out, the time consumption is long, the requirements on software and hardware cost are high, the cost performance of the prior art is low for the condition that only electric scanning radar signal warning needs to be realized, and a large amount of time cost and technical cost are applied to identification of mechanical scanning types.
Disclosure of Invention
The invention aims to provide a method for identifying a radar signal by using an electric scanning radar signal amplitude value, which can realize the quick identification of an enemy radar target by adopting an electronic scanning mode and realize the quick early warning of the radar target with high threat degree in the electronic reconnaissance process.
In order to achieve the purpose, the invention adopts the following technical scheme:
a method for identifying electric scanning radar signals by using intercepted signal amplitude values comprises the following steps:
acquiring a pulse signal amplitude-arrival time sequence of radar target radiation intercepted by an electronic reconnaissance receiver, and sequentially performing unit conversion and normalization processing on each pulse signal amplitude;
step two, acquiring a maximum main lobe sequence according to the pulse signal amplitude-arrival time sequence after normalization processing;
step three, carrying out autocorrelation processing on the amplitude of the pulse signal in the maximum main lobe sequence to obtain an autocorrelation function value sequence;
step four, determining a first minimum value point on the left side and a first minimum value point on the right side of the autocorrelation function value of 0 shift digit in the autocorrelation function value sequence, and calculating a first time length and a second time length;
step five, respectively determining a first main beam signal corresponding to the first time length and a second main beam signal corresponding to the second time length;
selecting the maximum main beam signal from the first main beam signal and the second main beam signal as a final main beam signal;
and seventhly, acquiring the flatness ratio of the final main beam signal and identifying the electric scanning radar signal.
Further, in the first step, the unit conversion formula is:
b i ’＝10 bi/20 ；
wherein, b i ' the amplitude of the ith pulse signal after unit conversion; b i I =1,2,3, …, N is the number of pulse signal amplitudes in the pulse signal amplitude sequence before unit conversion.
Further, in the first step, the formula of the normalization process is:
a i ＝b i ’/b max ；
b max ＝max[1,2,…,b i ’,…,b N ’]；
wherein, a i The amplitude of the ith pulse signal after normalization processing is obtained; b is a mixture of i ' the amplitude of the ith pulse signal after unit conversion; b max Maximum pulse signal in pulse signal amplitude sequence after unit conversionMagnitude of the sign.
Further, the specific implementation process of the step two is as follows:
step 21, obtaining the maximum pulse signal amplitude from the pulse signal amplitude-arrival time sequence after normalization processing;
step 22, reducing the left side and the right side of the maximum pulse signal amplitude to a preset pulse signal amplitude, and then obtaining corresponding 2 nearest time points;
and step 23, taking the normalized pulse signal amplitude-arrival time sequence between the 2 nearest time points as a maximum main lobe sequence.
Further, in step three, the autocorrelation function value is:
wherein R < R >]The autocorrelation function value with the shift digit r; "+" is the conjugate operation; a is j J =1,2, … and Q-r, wherein Q is the number of pulse signal amplitudes in the maximum main lobe sequence, and r is the number of shift bits of two autocorrelation pulse signal amplitudes.
Further, in step four, the autocorrelation function values of the left first minimum value point and the right first minimum value point respectively satisfy the following conditions:
R[r left ]<R[0]，R[r left ]<R[r left -1]；
R[r right ]<R[0]，R[r right ]<R[r right +1]；
wherein R < R > left ]And R < R > right ]The autocorrelation function values of the left first minimum value point and the right first minimum value point are respectively; r0]、R[r left -1]And R < R > right +1]Are 0 shift digit, r left -1 shift bit number, r right +1 shift the autocorrelation function value corresponding to the number of bits; r is left And r right And the shift digits are corresponding to the first minimum value point on the left side and the first minimum value point on the right side respectively.
In the fourth step, the first time length and the second time length are respectively:
wherein, t-legth 1 And t-legth 2 Respectively a first time length and a second time length;
and
are respectively the r-th right +1 and | r left The arrival time value corresponding to | +1 pulses, "|" is an absolute value; t is t 1 The time of arrival value corresponding to the 1 st pulse.
Further, in step five, the determining process of the first main beam signal is:
step 511, let j =0;
step 512, calculate
Step 513, letting j = j +1;
step 514, determine if j is less than Q-r right If so, then calculate
Step 515 is entered; if not, let m = j, go to step 516;
step 515, judge P new (j) Whether or not greater than P density1 If yes, let P density1 ＝P new (j) Returning to step 513; if not, return to step 513;
step 516, will
As the first main beam signal.
Further, in step five, the determining process of the second main beam signal is as follows:
step 521, let j =0;
step 522, calculate
Step 523, let j = j +1;
step 524, judging whether j is less than Q-r right If so, then calculate
Entering step 525; if not, let m = j, go to step 526;
step 525, judge P new (j) Whether or not greater than P density2 If yes, let P density2 ＝P new2 (j) Returning to step 523; if not, return to step 523;
step 526, will
As the second main beam signal.
Further, the specific implementation process of the step seven is as follows:
step 710, calculating the maximum amplitude value a of the final main beam signal max Minimum amplitude value a min Average amplitude value and average amplitude deviation value sigma;
711, according to the maximum amplitude value a max Minimum amplitude value a min And average amplitude deviation value sigma to determine the initial value a of amplitude value collection begin ＝a min σ, end value a end ＝a max + σ, statistical range a width ＝a end -a begin And counting the number of windows
Is to round up upwards;
step 712, setting the initial value of the amplitude value collection condition statistical vector sequence number z to be 1; setting a z-th amplitude value collection condition statistical vector a sta [z]Is 0; setting initial value a of amplitude parameter comparison window window For the initial value a in the case of a collection of amplitude values begin ；
Step 713, sorting the final main beam signal amplitude value sequence from large to small according to amplitude values;
714, setting an initial value of a sequence number e in the sequenced final main beam signal amplitude value sequence as 1;
step 715, judging whether z is larger than the number of the statistical windows, if so, entering step 717; if not, go to step 716;
step 716, when e<L and a are window <a e <a window When + sigma, let a sta [z]＝a sta [z]+1,e = e +1, return to step 715; when a is more than or equal to L, entering step 717; when a is e ≥a window When + sigma, let e = e +1,a window ＝a window + σ, z = z +1, return to step 715;
wherein, L is the number of amplitude values of the sequenced final main beam signal amplitude value sequence;
step 717, counting statistics vector a of collection of amplitude values sta [z]The medium maximum amplitude value and the final flatness of the main beam signal are determined;
step 718, judging whether the flatness of the final main beam signal is greater than or equal to a threshold value, if so, determining that the radar signal is an electric scanning radar signal; if not, the radar signal is not an electric scanning radar signal.
The invention has the beneficial effects that:
the method comprises the steps of performing unit conversion and normalization processing on the amplitude of each pulse signal in a pulse signal amplitude-arrival time sequence of radar target radiation intercepted by an electronic reconnaissance receiver to obtain a maximum main lobe sequence; carrying out autocorrelation processing on the amplitude of the pulse signal in the maximum main lobe sequence to obtain an autocorrelation function value sequence; determining a first minimum value point on the left side and a first minimum value point on the right side of the autocorrelation function value of 0 shift digit in the autocorrelation function value sequence and calculating a first time length and a second time length; respectively determining a first main beam signal corresponding to a first time length and a second main beam signal corresponding to a second time length; selecting the largest main beam signal from the first main beam signal and the second main beam signal as a final main beam signal; according to the pulse signal amplitude sequence of the final main beam signal, the flatness ratio of the main beam signal is calculated, electric scanning radar signal identification is carried out, rapid identification of radar targets of enemies in an electronic scanning mode is achieved, and rapid early warning of the radar targets with high threat degree in an electronic reconnaissance process is achieved.
Drawings
Fig. 1 is a schematic flow chart of an electric scanning radar signal identification method using an intercepted signal amplitude value according to the present invention.
Detailed Description
The following detailed description of the embodiments of the invention refers to the accompanying drawings.
In this embodiment, a method for identifying an electrical scanning radar signal by using an amplitude value of an intercepted signal is provided, and with reference to fig. 1, the method for identifying an electrical scanning radar signal includes the following steps:
s1, obtaining a pulse signal amplitude-arrival time sequence of radar target radiation intercepted by an electronic reconnaissance receiver, and sequentially performing unit conversion and normalization processing on each pulse signal amplitude.
Pulse signal amplitude value-arrival time sequence P = [ (b) of certain radar target radiation intercepted by electronic reconnaissance receiver of the embodiment 1 ,t 1 ),(b 2 ,t 2 ),…,(b i ,t i ),…,(b N ,t N )]I =1,2, …, N is the number of pulses, b i Is the i-th pulse signal amplitude, t, of the radar signal i Is the i-th pulse arrival time value of the radar signal, i.e. the time value of the radar signal arriving at the receiver.
The present embodiment performs unit conversion according to the following formula:
b i ’＝10 bi/20 ；
wherein, b i ' the amplitude of the ith pulse signal after unit conversion; b i I =1,2,3, …, N is the number of pulse signal amplitudes in the pulse signal amplitude sequence before unit conversion.
Converting the amplitude value from dBW unit to V unit to obtain pulse signal amplitude value-arrival time sequence P' = [ (b) expressed by absolute value 1 ’,t 1 ),(b 2 ’,t 2 ),…,(b i ’,t i ),…,(b N ’,t N )]。
And converting the amplitude value from a relative value to an absolute value, removing the pulse corresponding to the invalid pulse amplitude data, and normalizing the data. The pulse amplitude data is normalized according to the following formula:
a i ＝b i ’/b max ；
b max ＝max[1,2,…,b i ’,…,b N ’]；
wherein, a i The amplitude of the ith pulse signal after normalization processing is obtained; b i ' the amplitude of the ith pulse signal after unit conversion; b max The maximum pulse signal amplitude in the pulse signal amplitude sequence after unit conversion is adopted.
And S2, acquiring a maximum main lobe sequence according to the pulse signal amplitude-arrival time sequence after normalization processing.
And extracting pulse amplitude data corresponding to the main beam signal from the pulse signal amplitude-arrival time sequence after normalization processing. The specific implementation process is as follows:
and step 21, obtaining the maximum pulse signal amplitude from the pulse signal amplitude-arrival time sequence after the normalization processing.
From the normalized pulse signal amplitude-arrival time series p new The maximum pulse signal amplitude is extracted.
Step 22, reducing the left side and the right side of the maximum pulse signal amplitude to a preset pulse signal amplitude, and then acquiring corresponding 2 nearest time points;
and searching two corresponding nearest time points when the amplitude of the left and right side pulses of the maximum pulse signal amplitude is reduced to 0.09 times of the maximum value.
And step 23, taking the normalized pulse signal amplitude-arrival time sequence between the 2 nearest time points as a maximum main lobe sequence.
Taking the amplitude of Q pulse signals between two nearest time points as a maximum main lobe sequence p mainew ：
p mainew ＝{p[x]，t[x]}；
p[x]＝[a 1 ,a 2 ,…,a q ,…,a Q ]；
t[x]＝[t 1 ,t 2 ,…,t q ,…,t Q ]；
Wherein, a q The amplitude of the q pulse signal in the maximum main lobe sequence is obtained; t is t q The q pulse arrival time value in the maximum main lobe sequence; q =1,2,3, …, Q is the number of pulse signal amplitudes in the maximum main lobe sequence.
And S3, carrying out autocorrelation processing on the amplitude of the pulse signal in the maximum main lobe sequence to obtain an autocorrelation function value sequence.
In this embodiment, the autocorrelation function value is:
wherein R < R >]The autocorrelation function value with the shift digit r; "x" is the conjugate operation; a is j J =1,2, … and Q-r, wherein Q is the number of pulse signal amplitudes in the maximum main lobe sequence, and r is the number of shift bits of two autocorrelation pulse signal amplitudes.
And S4, determining a first minimum value point on the left side and a first minimum value point on the right side of the autocorrelation function value of 0 shift digit in the autocorrelation function value sequence, and calculating a first time length and a second time length.
In this embodiment, the autocorrelation function values of the left first minimum value point and the right first minimum value point respectively satisfy the following conditions:
R[r left ]<R[0]，R[r left ]<R[r left -1]；
R[r right ]<R[0]，R[r right ]<R[r right +1]；
wherein R < R > left ]And R < R > right ]The autocorrelation function values of the left first minimum value point and the right first minimum value point are respectively; r0]、R[r left -1]And R < R > right +1]Are 0 shift digit, r left -1 shift bit number, r right +1 shift the autocorrelation function value corresponding to the number of bits; r is left And r right The shift digits corresponding to the first minimum value point on the left side and the first minimum value point on the right side are respectively.
In this embodiment, the first time length and the second time length are respectively:
wherein, t-legth 1 And t-legth 2 Respectively a first time length and a second time length;
and
are respectively the r-th right +1 and | r left The arrival time values corresponding to | +1 pulses; "| x |" is an absolute value; t is t 1 The time of arrival value corresponding to the 1 st pulse.
And S5, respectively determining a first main beam signal corresponding to the first time length and a second main beam signal corresponding to the second time length.
The determination process of the first main beam signal in this embodiment is as follows:
step 511, let j =0;
step 512, calculate
Step 513, let j = j +1;
step 514, determine if j is less than Q-r right If so, then calculate
Step 515 is entered; if not, let m = j, go to step 516;
step 515, judge P new (j) Whether or not it is greater than P density1 If yes, let P density1 ＝P new (j) Returning to step 513; if not, return to step 513;
step 516, will
As the first main beam signal.
The determination process of the second main beam signal in this embodiment is as follows:
step 521, letting j =0;
step 522, calculate
Step 523, let j = j +1;
step 524, judge whether j is less than Q-r right If so, then calculate
Entering step 525; if not, let m = j, go to step 526;
step 525, judge P new (j) Whether or not greater than P density2 If yes, let P density2 ＝P new2 (j) Returning to step 523; if not, returning to the step523；
Step 526, will
As the second main beam signal.
And S6, selecting the maximum main beam signal from the first main beam signal and the second main beam signal as a final main beam signal.
When P is present density1 Greater than or equal to P density2 Then the first main beam signal P signal1 Is the final main beam signal; otherwise, the second main beam signal P signal2 Is the final main beam signal.
And S7, acquiring the flatness ratio of the final main beam signal and identifying the electric scanning radar signal.
And calculating the flatness ratio of the final main beam signal according to the amplitude value convergence condition of the final main beam signal, and judging whether the electric scanning radar signal exists according to the flatness ratio result. The specific implementation process is as follows:
step 710, calculating the maximum amplitude value a of the final main beam signal max Minimum amplitude value a min The average amplitude value and the average amplitude deviation value sigma;
the final main beam signal average amplitude value is:
wherein, the first and the second end of the pipe are connected with each other,
the average amplitude value of the final main beam signal; a is l L =1,2, …, L is the number of pulses of the final main beam signal.
The average amplitude deviation value σ is:
711, according to the maximum amplitude value a max Minimum amplitude value a min And average amplitude deviation value sigma, determining initial value a of amplitude value collection begin ＝a min σ, end value a end ＝a max + σ, statistical range a width ＝a end -a begin And counting the number of windows
Is to round up upwards;
step 712, setting the initial value of the amplitude value collection condition statistical vector sequence number z to be 1; setting a z-th amplitude value collection condition statistical vector a sta [z]Is 0; setting initial value a of amplitude parameter comparison window window For the initial value a in the case of a collection of amplitude values begin ；
Step 713, sorting the final main beam signal amplitude value sequence from large to small according to amplitude values;
714, setting an initial value of a sequence number e in the sequenced final main beam signal amplitude value sequence as 1;
step 715, judging whether z is larger than the number of the statistical windows, if so, entering step 717; if not, go to step 716;
step 716, when e<L, and a window <a e <a window When + sigma, let a sta [z]＝a sta [z]+1,e = e +1, return to step 715; when a is larger than or equal to L, entering step 717; when a is e ≥a window When + sigma, let e = e +1,a window ＝a window + σ, z = z +1, return to step 715;
wherein, L is the number of amplitude values of the sequenced final main beam signal amplitude value sequence;
step 717, counting statistics vector a of collection of amplitude values sta [z]The medium maximum amplitude value and the final flatness of the main beam signal are determined;
calculate the final principal according to the following formulaFlatness ratio F of beam signal signal ：
F signal ＝max(a sta [z])*100/L；
Wherein, F signal The final main beam signal flatness ratio is obtained; max (a) sta [z]) Counting a vector a for the collection of statistical amplitude values sta [z]A medium maximum amplitude value; and L is the number of amplitude values of the sequenced final main beam signal amplitude value sequence.
Step 718, judging whether the flatness of the final main beam signal is greater than or equal to a threshold value, if so, determining that the radar signal is an electric scanning radar signal; if not, the radar signal is not an electric scanning radar signal.
If F signal If the radar signal is larger than or equal to the threshold (such as 70), the radar signal is judged to be an electric scanning radar signal; otherwise, the radar signal is judged not to be the electric scanning radar signal.
In the embodiment, a maximum main lobe sequence is obtained by performing unit conversion and normalization processing on the amplitude of each pulse signal in a pulse signal amplitude-arrival time sequence of radar target radiation intercepted by an electronic reconnaissance receiver; carrying out autocorrelation processing on the amplitude of the pulse signal in the maximum main lobe sequence to obtain an autocorrelation function value sequence; determining a first minimum value point on the left side and a first minimum value point on the right side of the autocorrelation function value of 0 shift digit in the autocorrelation function value sequence and calculating a first time length and a second time length; respectively determining a first main beam signal corresponding to a first time length and a second main beam signal corresponding to a second time length; selecting the largest main beam signal from the first main beam signal and the second main beam signal as a final main beam signal; according to the pulse signal amplitude sequence of the final main beam signal, the flatness ratio of the main beam signal is calculated and the electric scanning radar signal is identified, so that the radar target of an enemy adopting an electronic scanning mode is quickly identified, and the quick early warning of the radar target with high threat degree in the electronic reconnaissance process is realized.
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 (7)
1. A method for identifying electric scanning radar signals by using intercepted signal amplitude values is characterized by comprising the following steps:
acquiring a pulse signal amplitude-arrival time sequence of radar target radiation intercepted by an electronic reconnaissance receiver, and sequentially performing unit conversion and normalization processing on each pulse signal amplitude;
step two, acquiring a maximum main lobe sequence according to the pulse signal amplitude-arrival time sequence after normalization processing;
step three, carrying out autocorrelation processing on the amplitude of the pulse signal in the maximum main lobe sequence to obtain an autocorrelation function value sequence;
step four, determining a first minimum value point on the left side and a first minimum value point on the right side of the autocorrelation function value of 0 shift digit in the autocorrelation function value sequence, and calculating a first time length and a second time length;
step five, respectively determining a first main beam signal corresponding to the first time length and a second main beam signal corresponding to the second time length;
selecting the maximum main beam signal from the first main beam signal and the second main beam signal as a final main beam signal;
seventhly, acquiring the flatness ratio of the final main beam signal and identifying an electric scanning radar signal;
the concrete implementation process of the seventh step is as follows:
step 710, calculating the maximum amplitude value a of the final main beam signal max Minimum amplitude value a min The average amplitude value and the average amplitude deviation value sigma;
711, according to the maximum amplitude value a max Minimum amplitude value a min And average amplitude deviation value sigma to determine the initial value a of amplitude value collection begin ＝a min σ, end value a end ＝a max + σ, statistical range a width ＝a end -a begin And counting the number of windows
Rounding up;
step 712, setting an initial value of the statistical vector sequence number z of the amplitude value aggregation condition to be 1; setting a z-th amplitude value collection condition statistical vector a sta [z]Is 0; setting initial value a of amplitude parameter comparison window window For the initial value a in the case of a collection of amplitude values begin ；
Step 713, sorting the final main beam signal amplitude value sequence from large to small according to amplitude values;
714, setting the initial value of the sequence number e in the sequenced final main beam signal amplitude value sequence as 1;
step 715, judging whether z is larger than the number of the statistical windows, if so, entering step 717; if not, go to step 716;
step 716, when e<L, and a window <a e <a window When + sigma, let a sta [z]＝a sta [z]+1,e = e +1, return to step 715; when a is more than or equal to L, entering step 717; when a is e ≥a window When + sigma, let e = e +1,a window ＝a window + σ, z = z +1, return to step 715;
wherein, L is the number of amplitude values of the sequenced final main beam signal amplitude value sequence;
step 717, counting statistics vector a of collection of amplitude values sta [z]The medium maximum amplitude value and the final flatness of the main beam signal are determined;
step 718, judging whether the flatness of the final main beam signal is greater than or equal to a threshold value, if so, determining that the radar signal is an electric scanning radar signal; if not, the radar signal is not an electric scanning radar signal.
2. An electric scanning radar signal identification method according to claim 1, characterized in that in step one, the unit conversion formula:
b i ’＝10 bi/20 ；
wherein, b i ' the amplitude of the ith pulse signal after unit conversion; bi is the ith pulse signal amplitude before unit conversion, i =1,2,3, …, N, N is the number of pulse signal amplitudes in the pulse signal amplitude sequence before unit conversion.
3. An electric scanning radar signal identification method according to claim 2, wherein in step one, the formula of the normalization process is:
a i ＝b i ’/b max ；
b max ＝max[1,2,…,b i ’,…,b N ’]；
wherein, a i The amplitude of the ith pulse signal after normalization processing is obtained; b i ' the amplitude of the ith pulse signal after unit conversion; b is a mixture of max The maximum pulse signal amplitude in the pulse signal amplitude sequence after unit conversion is adopted.
4. The method for identifying electric scanning radar signals according to any one of claims 1 to 3, wherein the concrete implementation process of the second step is as follows:
step 21, obtaining the maximum pulse signal amplitude from the pulse signal amplitude-arrival time sequence after normalization processing;
step 22, reducing the left side and the right side of the maximum pulse signal amplitude to a preset pulse signal amplitude, and then acquiring corresponding 2 nearest time points;
and step 23, taking the pulse signal amplitude-arrival time sequence after the normalization processing between the 2 nearest time points as a maximum main lobe sequence.
5. An electric scanning radar signal identification method according to any one of claims 1 to 3, wherein in step three, the autocorrelation function values are:
wherein R < R >]The autocorrelation function value with the shift digit r; "x" is the conjugate operation; a is a j J =1,2, … and Q-r, wherein Q is the number of pulse signal amplitudes in the maximum main lobe sequence, and r is the number of shift bits of two autocorrelation pulse signal amplitudes.
6. The method for identifying an electric scanning radar signal according to any one of claims 1 to 3, wherein in step four, the autocorrelation function values of the left first minimum value point and the right first minimum value point respectively satisfy the following conditions:
R[r left ]<R[0]，R[r left ]<R[r left -1]；
R[r right ]<R[0]，R[r right ]<R[r right +1]；
wherein R < R > left ]And R < R > right ]The autocorrelation function values are respectively the first minimum value point on the left side and the first minimum value point on the right side; r0]、R[r left -1]And R < R > right +1]Are 0 shift digit, r left -1 shift bit number, r right +1 shift the autocorrelation function value corresponding to the number of bits; r is left And r right And the shift digits are corresponding to the first minimum value point on the left side and the first minimum value point on the right side respectively.
7. An electric scanning radar signal identification method according to claim 6, wherein in step four, the first time length and the second time length are respectively:
wherein, t-legth 1 And t-legth 2 Respectively a first time length and a second time length;
and
are respectively the r-th right +1 and | r left The arrival time values corresponding to | +1 pulses; "| x |" is an absolute value; t is t 1 The time of arrival value corresponding to the 1 st pulse.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210183193.3A CN114578293B (en) | 2022-02-28 | 2022-02-28 | Electric scanning radar signal identification method by using intercepted signal amplitude value |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210183193.3A CN114578293B (en) | 2022-02-28 | 2022-02-28 | Electric scanning radar signal identification method by using intercepted signal amplitude value |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114578293A CN114578293A (en) | 2022-06-03 |
| CN114578293B true CN114578293B (en) | 2023-01-13 |
Family
ID=81771995
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210183193.3A Active CN114578293B (en) | 2022-02-28 | 2022-02-28 | Electric scanning radar signal identification method by using intercepted signal amplitude value |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114578293B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114814734A (en) * | 2022-07-01 | 2022-07-29 | 中国人民解放军国防科技大学 | Radar antenna scanning type identification method and system under low signal-to-noise ratio condition |
| CN117075061B (en) * | 2023-08-15 | 2024-01-30 | 扬州宇安电子科技有限公司 | Radar antenna measurement management system and method based on cloud computing |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107144835B (en) * | 2017-04-28 | 2020-01-21 | 安徽四创电子股份有限公司 | Low-altitude target monitoring method |
| CN108572353B (en) * | 2018-04-28 | 2023-03-28 | 西安电子科技大学 | Pulse time sequence planning method for low-interception radar |
| CN110426680B (en) * | 2019-07-20 | 2022-05-24 | 中国船舶重工集团公司第七二四研究所 | Multidimensional feature joint classification method based on interference signal time frequency and correlation coefficient |
| CN111707995A (en) * | 2020-06-19 | 2020-09-25 | 雷震烁 | Radar antenna scanning mode identification method based on limited-penetration visual effect |
| CN113109782B (en) * | 2021-04-15 | 2023-08-15 | 中国人民解放军空军航空大学 | Classification method directly applied to radar radiation source amplitude sequence |
| CN113239524B (en) * | 2021-04-26 | 2022-07-01 | 北京理工大学 | Method for calculating phase center of electric scanning array antenna based on moving rotating reference system simulation |
-
2022
- 2022-02-28 CN CN202210183193.3A patent/CN114578293B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN114578293A (en) | 2022-06-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114578293B (en) | Electric scanning radar signal identification method by using intercepted signal amplitude value | |
| CN109444820B (en) | Method for detecting target after interference suppression of multi-channel radar when clutter and interference coexist | |
| CN109444869B (en) | Radar extension target parameter adjustable detector for signal mismatch | |
| CN111965632B (en) | Radar target detection method based on Riemann manifold dimensionality reduction | |
| CN107015205B (en) | False target elimination method for distributed MIMO radar detection | |
| CN110232371B (en) | High-precision HRRP radar multi-target identification method based on small samples | |
| CN113253251B (en) | FDA-MIMO radar detection method and system when target speed is unknown | |
| CN112526456B (en) | Radar signal identification method | |
| CN113221631B (en) | Sequence pulse anti-interference target detection method based on convolutional neural network | |
| CN108872961B (en) | Radar weak target detection method based on low threshold | |
| CN110907908B (en) | Navigation radar signal sorting method based on envelope analysis | |
| CN109932727B (en) | Method for improving long-distance measurement precision in laser distance measurement system | |
| CN108535704B (en) | Signal pre-sorting method based on self-adaptive two-dimensional clustering | |
| CN112698267B (en) | Interference source testing and positioning method | |
| CN111796266B (en) | Tracking method before plane detection of uniformly accelerated moving target RD | |
| CN110988856B (en) | Target detection trace agglomeration algorithm based on density clustering | |
| CN107589419A (en) | A kind of method of energy peak joint wideband Range Profile Objective extraction | |
| CN116148794A (en) | Method and system for identifying airborne fire control radar scanning system | |
| CN114814734A (en) | Radar antenna scanning type identification method and system under low signal-to-noise ratio condition | |
| CN115356692A (en) | Radar signal sorting and batch increasing processing method based on non-overlapping time slices on PRI interval tree | |
| CN111123274B (en) | Target detection method of underwater sonar imaging system | |
| CN106371095A (en) | Pulse compression technique-based range imaging method and range imaging system | |
| CN111103588B (en) | Triangular wave multi-target identification method by utilizing signal energy | |
| CN111537976B (en) | Moving target detection method and device of airborne radar | |
| CN113341371B (en) | DOA estimation method based on L array and two-dimensional ESPRIT algorithm |
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 |