CN111474524A - Radar interference equipment interference effect monitoring and decision support system - Google Patents
Radar interference equipment interference effect monitoring and decision support system Download PDFInfo
- Publication number
- CN111474524A CN111474524A CN202010323179.XA CN202010323179A CN111474524A CN 111474524 A CN111474524 A CN 111474524A CN 202010323179 A CN202010323179 A CN 202010323179A CN 111474524 A CN111474524 A CN 111474524A
- Authority
- CN
- China
- Prior art keywords
- interference
- signal
- radar
- frequency
- signals
- 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.)
- Pending
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/38—Jamming means, e.g. producing false echoes
-
- 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/40—Means for monitoring or calibrating
Abstract
The invention discloses an interference effect monitoring and decision support system for radar interference equipment, which comprises a radar radiation source signal sampling module, an interference signal sampling module, a signal processing module, an interference effect evaluation module and an interference decision generation module; the interference effect evaluation module carries out quantitative analysis on the interference pattern selection and the interference parameters to obtain an interference effect evaluation result; the quantitative analysis process comprises distinguishing signals from different radar radiation sources, realizing pairing of a radar radiation signal source and an interference signal and analyzing the similarity of the signals. The method and the device perform interference effect evaluation on the interference signals from the angle of signal similarity, provide relatively accurate interference patterns and interference parameters according to the interference evaluation result, provide effective technical support for equipment operators to perform interference decision, really realize 'closed loop' of attack and defense drilling, and contribute to improving the fighting capacity of troops.
Description
Technical Field
The invention relates to the technical field of radar interference assessment, in particular to a radar interference equipment interference effect monitoring and decision support system.
Background
The radar is a very important military electronic device, and can obtain information of azimuth angle, pitch angle, speed, distance and the like of a target through the radar, which is called as 'eyes of electronic warfare'. With the development of radar technology, radar jamming technology has also been rapidly advanced. Radar jamming is electronic interference that degrades or even completely loses the performance of an enemy device by disturbing or spoofing the enemy's radar device. From the perspective of battle, the enemy radar system is destroyed and disturbed, so that the enemy radar system cannot accurately monitor the radar signal of the enemy, and tactics of successfully realizing the battle arrangement of the enemy are guaranteed.
The radar interference equipment is limited by various factors, does not generally emit strong-power interference signals in combat training, so that trained personnel have little knowledge of the characteristics of the interference signals, and often cannot judge whether the adopted interference patterns and the set interference signal parameters are reasonable or not, and cannot evaluate the interference effect on a target radar. Summary active radar jamming equipment presents three main problems: 1. the rationality of the interference signals by the actual installation is difficult to monitor in real time, and whether the adopted interference pattern and the set interference signal parameters are reasonable or not cannot be judged, so that the daily training of troops and the teaching of colleges mostly stay in the theoretical level, and the interference effect cannot be evaluated; 2. for a specific radar signal, no interference decision support exists; 3. the equipment terminal has no visual signal visual display, and mostly takes parameter display as a main point.
Disclosure of Invention
Aiming at the technical defects of the existing radar interference equipment, the invention aims to provide a radar interference equipment interference effect monitoring and decision support system, which provides effective technical support for equipment operators to make interference decisions by analyzing the similarity of radar radiation source signals and interference signals, and giving interference effect evaluation of the interference signals, optimal interference patterns and interference parameters.
In order to solve the technical problem, the interference effect monitoring and decision support system for radar interference equipment disclosed by the invention comprises
The radar radiation source signal sampling module is used for sampling radar radiation source signals from a signal interface of the radar interference equipment;
the interference signal sampling module is used for sampling an interference signal from a signal interface of the radar interference equipment for sampling;
the signal processing module is used for carrying out time-frequency characteristic analysis on the sampled radar radiation source signals and the sampled interference signals;
the interference effect evaluation module is used for monitoring the interference effect of the interference signal in real time according to the signal similarity, and carrying out quantitative analysis on the rationality of the selection of the interference pattern and the accuracy of the setting of the interference parameters to obtain an interference effect evaluation result; the quantitative analysis process comprises the steps of distinguishing signals from different radar radiation sources through a signal sorting method based on Q MF analysis, realizing the matching of a radar radiation signal source and an interference signal through a signal sorting method based on rising edge waveform matching, and finally completing the signal similarity analysis through waveform similarity measurement;
and the interference decision generation module is used for giving a relatively accurate interference pattern and interference parameters on the basis of the interference effect evaluation result, and providing effective technical support for equipment operators to make interference decisions.
Further, the signal processing adopts a ZFT frequency measurement method to carry out frequency domain analysis on the signal, and the specific operation is that firstly, the acquired signal X is subjected to frequency shift to be changed into X', then, the X is subjected to low-pass filtering and re-sampling, and finally, fast Fourier transform is carried out, or the sectional FFT phase weighting frequency measurement method is adopted to carry out frequency domain analysis on the signal, the specific operation is that X (n) is divided into M groups which are not overlapped, the data length is L points, and the sampling data is expressed as si＝[x(iL),x(iL+1),…,x(iL+L)]Wherein i is 0,1, …, M-1, and performing L points FFT operation on the M groups of data to obtain siAfter the amplitude is accumulated, the signal frequency is roughly measured, and then phase weighting is carried out, and the signal frequency is accurately measured.
Further, signal processing adopts a time-frequency analysis method based on WVD to carry out time-frequency analysis on signals, and the specific operations are ① obtaining sampling signals and carrying out Hilbert transformation to obtain analysis signals of the sampling signals, ② determining window types, window lengths and window moving steps, ③ moving a window, calculating instantaneous autocorrelation of data in the window by taking a window center as a reference point, ④ carrying out fast Fourier transformation on the instantaneous autocorrelation, ⑤ moving the window by window moving steps and repeating the steps ③, ④ and ⑤ until all data are processed.
Further, the specific operation of distinguishing signals from different radar radiation sources by a signal sorting method based on QMF analysis is that the QMF analysis is adopted to realize orthogonal wavelet decomposition of the signals, the input signal energy is divided into orthogonal components in frequency by a dual-channel orthogonal mirror image filter bank, then the formed QMF pairs are introduced into a tree structure with a certain number of layers, the 'blocks' in each layer have the same area size by waveform decomposition, the output components of each filter are input into the QMF pairs in the next layer, each QMF pair decomposes the input signal waveform into two parts of high-frequency components and low-frequency components by taking pi as a boundary, and the 'blocks' refer to rectangular regions containing the energy of basis functions on a time-frequency plane. The quadrature mirror filter herein employs an improved SINC filter having a transfer function of
Wherein N is more than or equal to 2 and less than or equal to N2, (N-2) N is more than or equal to 2, C is a compression variable, S is a scale variable, N is the number of convolution points, and omega (N) is a Hamming window function for inhibiting Gibbs phenomenon.
Furthermore, the matching of the radar radiation signal source and the interference signal is realized through a signal sorting method based on the rising edge waveform matching, and the method has the specific operations of ① extracting the pulse signal envelope and carrying out curve fitting, ② extracting the rising edge of the pulse signal envelope after the curve fitting, and ③ judging the similarity of the radar radiation source signal and the rising edge waveform of the interference signal by adopting a Huasdorff distance method to complete the matching of the radar radiation signal source and the interference signal.
Further, signal similarity analysis is completed through waveform similarity measurement, and the method specifically comprises the steps of ① performing overall comparison on radar radiation source signals and interference signal waveforms by combining an included angle cosine algorithm and a sliding window algorithm to obtain waveform basic similarity, ② comparing amplitudes of the radar radiation source signal waveforms and the interference signal waveforms by adopting an average absolute difference algorithm to obtain amplitude similarity, and ③ performing average weighting on the waveform basic similarity and the amplitude similarity to obtain a final similarity value of the radar radiation source signal waveforms and the interference signal waveforms.
Furthermore, the evaluation standard of the interference effect evaluation result is that the interference effect is good when the final similarity value is more than 95%, the interference effect is general when 80% -95%, the critical interference is 65% -80%, and the interference effect is no more than 65%.
Further, the radar interference equipment is provided with a display module for displaying the aircraft track, the radar receiving signal, the interference signal, the current interference pattern and interference parameter, the interference effect evaluation result, the optimal interference pattern and the interference parameter; the radar receiving signal and the interference signal are visually displayed by adopting a signal rolling method, and the current interference pattern and the interference parameter are selected through a menu bar, a shortcut key or a window; the optimal interference pattern and the interference parameters are obtained and displayed by the interference decision generation module according to the radar type and the interference effect evaluation result.
The method and the device perform interference effect evaluation on the interference signals from the angle of signal similarity, provide relatively accurate interference patterns and interference parameters according to the interference evaluation result, provide effective technical support for equipment operators to perform interference decision, really realize 'closed loop' of attack and defense drilling, contribute to improving the fighting capacity of troops, and solve the problem that most of daily training and college teaching stay at the theoretical level.
Drawings
FIG. 1 is a block diagram of a radar jamming effect monitoring and decision support system;
FIG. 2 is a schematic flow chart of a ZFT frequency measurement method;
FIG. 3 shows the frequency k0A nearby spectrum spread;
FIG. 4 is a flow chart of a WVD process;
FIG. 5 is a diagram of a quadrature mirror filterbank tree structure;
FIG. 6 is a contour plot of an L FM signal 5-level QMF analysis;
FIG. 7 is a contour plot of a Barker signal 3-layer QMF analysis;
FIG. 8 is a contour plot of costas signal 3-layer QMF analysis;
FIG. 9 is a contour plot of a 2-layer QMF analysis of the frank signal;
FIG. 10 is a waveform of a pulse envelope of a radar radiation source signal;
FIG. 11 is a flow chart of a method for sorting signals based on leading edge waveform pairing;
FIG. 12 is a diagram of a dynamic sliding window algorithm;
fig. 13 is a flow chart of a segmented FFT phase weighted frequency measurement method.
Detailed Description
The present invention will be described in further detail with reference to the accompanying drawings and specific embodiments. The embodiments of the present invention have been presented for purposes of illustration and description, and are not intended to be exhaustive or limited to the invention in the form disclosed. Many modifications and variations will be apparent to those of ordinary skill in the art. The embodiment was chosen and described in order to best explain the principles of the invention and the practical application, and to enable others of ordinary skill in the art to understand the invention for various embodiments with various modifications as are suited to the particular use contemplated.
Example 1
A radar jamming effect monitoring and decision support system, as shown in FIG. 1, comprises
The radar radiation source signal sampling module is used for sampling radar radiation source signals from a signal interface of the radar interference equipment;
the interference signal sampling module is used for sampling an interference signal from a signal interface of the radar interference equipment for sampling;
the signal processing module is used for carrying out time-frequency characteristic analysis on the sampled radar radiation source signals and the sampled interference signals;
the interference effect evaluation module is used for monitoring the interference effect of the interference signal in real time according to the signal similarity, and carrying out quantitative analysis on the rationality of the selection of the interference pattern and the accuracy of the setting of the interference parameters to obtain an interference effect evaluation result; the quantitative analysis process comprises the steps of distinguishing signals from different radar radiation sources through a signal sorting method based on Q MF analysis, realizing the matching of a radar radiation signal source and an interference signal through a signal sorting method based on rising edge waveform matching, and finally completing the signal similarity analysis through waveform similarity measurement;
and the interference decision generation module is used for giving a relatively accurate interference pattern and interference parameters on the basis of the interference effect evaluation result, and providing effective technical support for equipment operators to make interference decisions.
In the aspect of time-frequency characteristic analysis of signals, frequency domain analysis and time-frequency analysis are performed in the embodiment.
The frequency domain analysis mainly realizes high-precision frequency measurement of signals and provides signal frequency spectrums. The frequency measurement is performed on the acquired digital signal by a zero-crossing point method and a Fast Fourier Transform (FFT) method. For pure signals which are not polluted by noise, the zero crossing point method is a fast method. However, in a real environment, there is no clean radar signal. If the frequency of the radar signal contaminated by noise is measured by the zero crossing method, the error will reach an intolerable degree. FFT is a classical spectral analysis method and is widely used for spectral measurement and analysis. However, FFT is limited in high precision frequency measurement. Because the frequency resolution of the FFT, Δ f ═ fsN, wherein fsIs the sampling frequency, and N is the FFT transform length; f. ofsThe reduction of (a) is limited by the nyquist frequency and cannot be reduced without limit; the increase of N is limited by processing time and computer memory and cannot be increased without limit.
In this embodiment, a ZFFT frequency measurement method is adopted to perform frequency domain analysis on a signal, and specifically, the operation is to shift the frequency of the acquired signal X into X', then perform low-pass filtering and resampling, and finally perform fast fourier transform, as shown in fig. 2. For example, if the estimated X bandwidth is about 1MHz, the frequency-shifted and low-pass filtered y (N) may be in the range of 0-1MHz, and if the sampling frequency of re-sampling is 5MHz, the sampling frequency is reduced by 100 times compared with the original sampling frequency of 500MHz, and when the transform length N is not changed, the frequency measurement accuracy is improved by 100 times.
The frequency of the signal appears as a maximum point on the spectrum, i.e. the frequency of the signal is measurable by the position of the maximum on the spectrum. Since the output result of FFT is discrete, the maximum frequency spectrum division of FFT outputThe quantity point has a certain deviation from the true frequency of the signal (the maximum point of the spectral peak), namely, the estimation problem of the frequency spectral peak exists. In order to further improve the frequency measurement precision and obtain the true signal frequency, the output result needs to be corrected. Referring to fig. 3, it is assumed that the first maximum value of the FFT is X0Frequency value of k0The second maximum is X1Frequency value of k1. If the FT operation adopts Hanning window, the correction formula is
By the method, the signal frequency measurement precision can reach within 5 KHz.
The time-frequency analysis mainly analyzes the characteristics of signals with frequency changing along with time, and the signals comprise frequency modulation signals, frequency coding signals and the like. The embodiment completes the identification of the modulation mode in radar signal pulse based on the WVD non-stationary signal processing technology, and completes the extraction of the modulation characteristic on the basis.
The Wigner-Willi distribution (WVD) is a non-stationary signal processing technology, the analysis of signals focuses more on local characteristics, a certain resolution is achieved in both time domain and frequency domain, especially for L FM signals, the time-frequency aggregation performance is approximate to ideal, the WVD of the signals is subjected to first-order moment analysis, the instantaneous frequency of the signals can be obtained, and therefore intra-pulse phase jump points of the phase coding radar signals can be extracted.
WVD of the real signal x (t) is defined as
Wherein s (t) is an analytic signal of the signals x (t), and
then d τ is 2d λ, and the formula is substituted to obtain
If order
Can also be expressed as
This is typically a fourier transform, called the instantaneous autocorrelation of the signal.
By definition, WVD can be interpreted as that the frequency spectrum corresponding to each time instant is centered at the time instant, all parts of the signal on the left and right sides of the time instant are folded and multiplied, and then the multiplied result is fourier transformed.
The calculation of WVD is a non-causal operation. In practice, the signal is within a limited interval, which can be overcome by introducing a delay. If the delay allowed for the analysis is small relative to the signal interval, then the time t can be determined by comparing the signal0Is defined by applying WV D, and the windowed signal is sw(t,t0)＝s(t)w(t-t0). The WVD obtained by replacing the original analytic signal with the windowed signal is called pseudo-WVD.
The pseudo WVD is beneficial to implementation and is suitable for real-time processing. Firstly, discretizing an input signal, wherein the discretization WVD is defined as
Then fixing the window, moving the data in sequence, inputting the data in sequence after the data in the window is processed, and circulating in sequence until the data is processed, wherein the method specifically comprises the following steps:
1. acquiring a sampling signal, and performing Hilbert conversion to obtain an analytic signal of an original signal;
2. determining the window type, the window length and the window moving step;
3. moving the window, and calculating instantaneous autocorrelation of data in the window by taking the center of the window as a datum point;
4. performing fast Fourier transform on the instantaneous autocorrelation;
5. and (5) moving the window step by the window movement, and repeating the steps 3-5 until all data are processed.
The operation speed is greatly dependent on the size delta step of the window moving step, and in the process of detecting L FM signals by adopting a WVD method, the key is to obtain the intra-pulse time-frequency energy distribution trend, so the influence of the size of the delta step on the detection result is not large, the intra-pulse time-frequency energy distribution trend can be ensured to be kept unchanged as long as the delta step is ensured to be smaller than the window width, if the value of the delta step is larger than the window width, a single-component signal is wrongly judged as a multi-component signal.
The WVD can be used to extract the instantaneous frequency of the signal, in addition to revealing the time-frequency energy distribution of the signal. The instantaneous frequency of signal x (t) is represented as
Wherein s (t) is the analytic signal of x (t). If the analytic signal of the signal is expressed as a vector of a complex plane and the rotating speed of the amplitude angle of the vector is expressed, the first moment of the WVD of the signal is directly proportional to the instantaneous frequency of the signal.
Assuming that the sequence length of the discrete signal is N, the instantaneous frequency of the signal is
Wherein f issFor the sampling frequency, M is the number of WVD frequency points, M1(n) is the first moment of WVD,
an instantaneous discrete frequency estimate based on the W VD can be obtained as
Chirp signals and phase-coded signals are typical representatives of PC signals and have found a large number of applications in new-system radars. The two signals are detected by adopting a WVD method, and the energy distribution of the time-frequency surface of the linear frequency modulation radar signal is mainly observed from two aspects; for the phase encoding radar signal, the instantaneous frequency of the signal is mainly extracted, and the phase jump point is extracted according to the frequency jump point.
Let x (t) be cos (2 pi f)0t+kπt2) Is a linear frequency modulation signal with constant amplitude, and the analytic signal is
Then calculate its instantaneous autocorrelation to obtain
And further obtaining a linear FM signal WVD of
It can be seen that the WVD of the single component chirp signal is along the straight line f ═ f0The impact line spectrum of the + kt distribution, the amplitude of which appears concentrated on a line representing the instantaneous frequency variation of the signal, therefore WVD has the best time-frequency convergence at the frequency modulation rate at which the chirp signal is exhibited.
The phase coding signal divides the pulse width T into N sub-pulses, the sub-pulse width is tau-TN, then the sub-pulses are coded by the phase of the carrier wave, fixed phase shift exists between each sub-pulse, and the magnitude of the phase shift is related to the coding rule. Because phase jump points exist among all the subcodes, frequency jump exists on the instantaneous frequency, and the extraction of the phase jump points of the phase encoding radar signals can be completed by extracting the instantaneous frequency.
From the perspective of signal processing, it is the most essential means to destroy radar coherence or make interference signals have coherence the same as radar radiation source echo, i.e. generate interference signals with similar structure and similar characteristics of radar radiation source echo, make them effectively enter radar receiver, obtain processing gain, make radar difficult to effectively distinguish radar radiation source echo and interference in the processes of receiving, processing, judging, parameter estimation, imaging, target monitoring and classification identification, thereby a large amount of false alarms appear in high resolution SAR image monitoring, or force the constant false alarm threshold to rise, and greatly reduce the monitoring probability of target. Based on the interference mechanism of signal similarity, the interference signal must have the similar characteristics of SAR radar target echo signals, so that the effective accumulation of energy is realized, and the effective interference is implemented on the SAR radar, so that the similarity of a radar radiation source signal and the interference signal becomes an important index for examining the interference effect of the equipment. In order to obtain the signal similarity, the signal characteristics need to be extracted and characteristic analysis needs to be carried out to distinguish the signals from different radar radiation sources. With the complication of signal waveforms, the emergence of radar signals of a new system with flexible and changeable parameter characteristics only depends on a sorting method based on extracted signal parameter characteristics, and the requirements of signal sorting in modern electronic warfare are difficult to meet.
In the embodiment, a signal sorting method based on QMF analysis is adopted, and under a proper decomposition layer number, different radar radiation source signals are distinguished by extracting image characteristics of the signals, so that the problem that the sorting effect of multi-phase and multi-frequency signals is poor due to cross item interference caused by bilinear transformation existing in WVD time-frequency analysis can be solved. QMF analysis is used to perform an orthogonal wavelet decomposition of a radar signal, with the input signal energy being split in frequency into orthogonal components by a two-channel quadrature mirror filter bank. In the time-frequency plane, the high frequency wavelet transform is more sharp in time and the low frequency wavelet transform is more sharp in frequency. That is, as frequency increases, "chunks" become shorter in time and occupy more frequency bandwidth. Since the wavelet transform is linear, there is a fundamental limit to the minimum area of these "blocks". The present embodiment introduces the formed QMF pairs into a tree structure having a certain number of layers, makes "blocks" in each layer have the same area size by waveform decomposition, and inputs the output components of each filter into the QMF pairs of the next layer, as shown in fig. 5, each QMF pair inputs the input components in pi boundariesThe signal waveform is decomposed into two parts of high-frequency components and low-frequency components, and the block refers to a rectangular area containing the energy of the basis function on a time-frequency plane. The quadrature mirror filter uses an improved SINC filter having a transfer function of
Wherein N is more than or equal to (N-2)/2 and less than or equal to N2, C is a compression variable, S is a scale variable, N is the number of convolution points, and omega (N) is a Hamming window function for inhibiting Gibbs phenomenon.
Let the normalized signal be input at a rate of one sample per second with a signal bandwidth of 0, pi]The output of each layer of the tree structure in fig. 5 constitutes a matrix, the sum of squares of the elements in the matrix approximately represents the energy contained in a block in the corresponding time-frequency diagram, the length (in time) and the width (in frequency) of the coefficient matrix of the block output of each layer are respectively 2 times and 1/2 times of the coefficient matrix of the upper layer, the resolution of the structure to the signal is different, the higher the number of layers is, the lower the time resolution is, but the higher the frequency resolution is, the characteristic is that the nth layer has 2n QMF pairs, the frequency length of the time-frequency resolution window of the previous layer is 2 times that of the next layer, the time length is 1/2 of the next layer, the requirement of multi-resolution is met, if the total number of layers is L, the sampling frequency of the input signal is fsThen for the l (l ≦ L) th layer, the frequency resolution is
With a temporal resolution of
For L PI radar signals, by selecting proper number of layers, the modulation patterns and modulation parameter characteristics of different signals can be easily identified, such as the bandwidth, carrier frequency, energy distribution in 'chip', phase modulation, pulse width, position information on a time-frequency diagram and the like can be determined, and even the number of current transmitters and the type of L PI transmitters can be identified.
Selected L FM signal parameter sampling frequency fs40MHz, pulse width τ 10 μ s, signal start frequency fc2MHz, and bandwidth B4 MHz. Without noise pollution, 5-layer processing was performed by QMF analysis to obtain the results of fig. 6.
The selected Barker code signal parameters are as follows: sampling frequency fs40MHz, 13 μ s pulse width T, 1 μ s sub-pulse width T, and signal center frequency f 010 MHz. Without noise pollution, the results of fig. 7 were obtained by performing 3-layer processing through QMF analysis.
Selected 10-bit costas signal parameters: sampling frequency fs40MHz, the coding sequence is [2,4,8,5,10,9,7,3,6,1]Duration of sub-frequency T p1 mus, signal start frequency f c2 MHz. Without noise pollution, the results of fig. 8 were obtained by performing 3-layer processing through QMF analysis.
Selected frank code signal parameters: sampling frequency fs40MHz, sub-pulse width τ 1 μ s, N3, i.e. code length 329, pulse width T9 mus, signal carrier frequency f 010 MHz. Without noise pollution, the results of fig. 9 were obtained by performing a 2-layer process by QMF analysis.
For radar signals with different modulation patterns, the radar signals can be clearly distinguished through QMF analysis of appropriate layer numbers, the modulation information of the signals can be accurately identified, such as the modulation slope of L FM signals, the phase mutation point positions of barker code and frank code signals and the like, a sufficiently reliable basis is provided for signal sorting, meanwhile, certain anti-noise performance is achieved, and the requirements of signal sorting accuracy and instantaneity are met.
After signals from different radar radiation sources are distinguished by a signal sorting method based on QMF analysis, in this embodiment, waveforms with fingerprint characteristics are extracted by transforming or mapping the signals by a signal sorting method based on leading edge waveform pairing, and the waveforms are used for pairing the radar radiation source signals with interference signals. The significance of further comparing the similarity of the signals is achieved only after the pairing relation between the radar radiation source signals and the interference signals is established.
The core of radar radiation source signal sorting is feature extraction, the separability of selected features determines the effectiveness of a sorting method, parameter features between pulses and fine features in the pulses are mainly adopted in the prior art, and the problems of batch increase and batch omission generally exist in sorting results. Therefore, the method based on the radar individual identification result can realize the most effective signal sorting. Because the electronic components used in any single radar radiation source have inevitable errors and the internal mechanical structure has extremely small difference, the envelope waveform of the pulse signal shows complex nonlinear characteristics, and the difference on the pulse envelope waveform of the radar radiation source signal can be said to exist inevitably, even if two radars of the same model are produced by the same manufacturer in the same batch at the same time. Thus, the hardware construction shows that the radar radiation source signal sorting based on the difference of the signal pulse envelope waveform has certain feasibility.
The pulse envelope waveform of a certain radar radiation source signal is shown in fig. 10. The received radar radiation source signal is influenced by noise interference, channel fading, multipath effect and the like, and the envelope waveform of the pulse signal is necessarily distorted in the process of being transmitted to a receiving antenna from an electromagnetic space, so that the shape of the envelope fluctuates and changes, wherein the multipath effect has the greatest influence on the envelope of the pulse signal. Research shows that in the parameter characteristics of the envelope of the signal, the front edge waveform of the envelope is influenced the least by the multipath effect, and the back edge waveform and the pulse width of the envelope are influenced the most by the multipath effect. Therefore, only the research on the overall envelope waveform of the signal from a macroscopic angle is considered, the influence caused by the multipath effect is avoided, and the rising edge of the signal envelope is selected as the fingerprint characteristic of the radar radiation source. Referring to fig. 11, the specific steps are as follows:
1. extraction of pulse envelope
The waveform obtained by directly extracting the signal envelope is not ideal, and great error influence is caused to the extraction of the subsequent rising edge and the waveform pairing. For this problem, the embodiment performs curve fitting on the preliminarily extracted signal envelope, and the curve fitting method may adopt a cubic spline interpolation method.
2. Signal sorting based on leading edge waveform pairing
For any two groups of rising edge waveforms (radar radiation source signals and interference signals), the higher the matching degree of the waveforms is, the more similar the waveforms of the two rising edges are, and the pulse pairing is realized. The present embodiment uses the Hausdorff distance method to calculate the similarity between rising edge waveforms.
The Huasdorff is a max-min type distance function defined between two point sets, is mainly used for calculating the similarity between the two point sets, and does not need to consider whether points in the two point sets have a certain corresponding relation or not, so that the similarity between the two point sets can be calculated on the premise of no mapping relation.
Let two sets of finite points as A ═ a1,a2,…,anB ═ B1,b2,…,bnThen, the Hausdorff distance function between point sets a and B can be defined as H (a, B) ═ max { H (a, B), H (B, a) }. h (A, B) is a one-way Hausdorff distance function from point set A to point set B, having
Conversely, h (B, A) is called a one-way Hausdorff distance function from point set B to point set A, having
Where | | is a class of distance norm, which may be L2And L∞Etc., typically L2I.e. the euclidean distance. If the distance d from any point a in the point set A to the point set B is definedB(a)Set A to BMinimum value of distance of each point, i.e.
Then h (A, B) is the distance d from all points in the point set A to the point set BB(a)Is understood to mean that the maximum distance in the set consisting of all minimum distance values is taken as h (a, B). Defining the standard template as a and the waveform to be matched as B, H (a, B) can be used to measure the matching degree between A, B, and similarly H (B, a) represents the matching degree between B and a, and then H (a, B) completely measures the similarity between any two point sets A, B, and is generally called a bidirectional Hausdorff distance function between A, B.
3. Cluster analysis
In the above-described waveform matching concept, the selection of the rising edge waveform and the number of sampling points, which refers to the number of all points of a single complete envelope, has a large influence on the sorting result. For the problem, in this embodiment, a rising edge of a pulse envelope is extracted from each radar radiation source as a template waveform, a Hausdorff distance is calculated with an interference signal waveform, and a cluster analysis is performed on an obtained Hausdorff distance value to complete signal sorting. The signal sorting method based on the rising edge waveform pairing can effectively solve the problems of batch increase and batch leakage commonly existing in the current signal sorting method.
And after the radar radiation source signal and the interference signal are matched, signal similarity analysis is completed through waveform similarity measurement. The existing similarity measurement theory is a theory for describing the similarity of shapes, images, signals or data, and has great significance in many fields. At present, the similarity measurement theory has been widely applied to the similarity measurement of signal waveforms. For example, the similarity of signals is measured by using a regional correlation method, two signals are uniformly divided into a plurality of time segments, and the optimal matching segment of the two signals is searched; the dynamic time normalization method is applied to measuring the similarity of the electrocardiogram signals, is sensitive to the amplitude change of the signals, but is not sensitive to the direct current offset of the whole signals; the similarity of the waveforms is measured by using a skeleton tree matching method, the method is insensitive to nonlinear fluctuation on a signal time axis, but when the peak-valley information of the waveforms is complex and the number of nodes of the tree is increased, the operation time is prolonged; the similarity between two signals is calculated by adopting an included angle cosine algorithm, and a good effect is obtained when the included angle cosine algorithm is applied to the process of calculating the similarity of the traditional Chinese medicine fingerprint, but the included angle cosine algorithm cannot well reflect the amplitude change of the waveform, particularly the proportional change of the amplitude and the direct current offset of the signals because the included angle cosine algorithm calculates the overall correlation of the waveform, and cannot embody the local characteristics of the waveform.
The signal similarity analysis method based on the waveform similarity measurement is characterized in that a pulse signal sequence of a radar radiation source signal and an interference signal is subjected to signal preprocessing, a weighted average method is provided by combining and applying an average absolute difference algorithm and a dynamic sliding window algorithm on the basis of an included angle cosine algorithm, the similarity between the radar radiation source signal waveform and the interference signal waveform is measured, basic similarity and amplitude similarity indexes of the waveform are respectively extracted, the measured similarity values are subjected to weighted average processing, the final result of waveform similarity is obtained, and the influence caused by non-coherent signals is avoided. The method comprises the following specific steps:
1. signal pre-processing
In order to ensure the comparability between signal data, before similarity measurement is carried out on the signals, the signals need to be subjected to data processing in a non-dimensionalization mode, and the standardized formula is that
2. Measuring basic similarity of waveforms
The current radar technology mainly uses a phase-coherent system, when measuring the basic similarity of a radiation source waveform and an interference waveform, the waveform needs to be aligned in phase, and the problem in the matching calculation process is solved by combining an included angle cosine algorithm and a dynamic sliding window algorithm.
Matching calculation is performed m times on the radar radiation source signal waveform and the interference signal waveform, as shown in fig. 12. In order to keep the lengths of the waveform data participating in matching consistent while the phases are aligned, each time matching is carried out, a comparison window of interference signal waveform data slides backwards by one bit and the length of the window is reduced by one bit, and the operation of the radar radiation source signal waveform is the same, but the sliding directions are opposite. After m times of matching calculation, the largest one of the m matching values is taken as the basic similarity result of the waveform matching, namely
The purpose of matching is to obtain the similarity under the condition that the basic similarity is closest to the ideal matching condition, so as to judge whether the basic similarity is qualified or not. And in the matching values obtained by m times of matching calculation, taking the maximum value to represent the waveform matching condition when the waveform is relatively close to the ideal state, and reflecting the relatively real matching degree between the radar radiation source signal waveform and the interference signal waveform.
If the waveform data length is N, there are N periodic waveforms. For accurate calculation, let the step size of the sliding window be 1 and the total sliding distance be m. Repeated simulation experiments show that when m is 2N/N and N is more than or equal to 3, the reliability of the obtained matching value is highest.
3. Measuring waveform amplitude similarity
Because the included angle cosine algorithm is insensitive to the amplitude change of the waveform, particularly the proportional change of the amplitude and the deviation of the signal, the basic similarity value calculated by the algorithm cannot effectively reflect the change condition of the amplitude of the compared waveform. However, the waveform amplitude variation has a great influence on the waveform similarity, so that the amplitude of the interference signal waveform must be separately processed to calculate the waveform amplitude similarity.
The embodiment adopts an average absolute difference algorithm to detect overall proportional change and deviation of the waveform amplitude of the interference signal. Suppose, the interference signal waveform A and the radar radiation source signal waveThe shape M is two waveforms after phase alignment, the data length of the waveforms is h, then di＝|Ai-MiL. Here, diThe absolute difference between the amplitude of the ith element of the interference signal waveform and the radar radiation source signal waveform after phase alignment can visually represent the change condition of the amplitudes of the two waveforms, but cannot be directly used for calculating the similarity of the amplitudes of the two waveforms, and further conversion is needed. The relative variation of the waveform amplitude is
After normalization is
Wherein DminAnd DmaxThe minimum value and the maximum value of the relative variation of the waveform amplitude are obtained. Average value of amplitude similarity of interference signal waveform and radar radiation source signal waveform
4. Weighted average processing
And the final similarity value of the radar radiation source signal waveform and the interference signal waveform is obtained by performing weighted average on the basic similarity value and the amplitude similarity value, wherein the weight of the similarity index can be set in advance by manually combining the prior condition, and is adjusted according to the actual condition and the actual requirement.
Final similarity value S ═ SAM*α-Saβ (0,1), α and β respectively represent overall similarity values SAMDegree of dissimilarity S of waveform amplitudesaThe weight occupied in the final similarity value calculation.
The evaluation standard of the interference effect evaluation result is that the interference effect is good when the final similarity value is more than 95%, the interference effect is general when 80% -95%, the critical interference is 65% -80%, and the interference effect is no more than 65%.
In addition, the radar interference equipment is provided with a display module for displaying the aircraft track, the radar receiving signal, the interference signal, the current interference pattern and interference parameter, the interference effect evaluation result, the optimal interference pattern and the interference parameter; the radar receiving signal and the interference signal are visually displayed by adopting a signal rolling method, and the current interference pattern and the interference parameter are selected through a menu bar, a shortcut key or a window; the optimal interference pattern and the interference parameters are obtained and displayed by the interference decision generation module according to the radar type and the interference effect evaluation result, and are displayed based on the most intuitive signals and data of equipment operators.
Example 2
The difference between this embodiment and embodiment 1 is that the frequency domain analysis is performed on the signal by using the segmented FFT phase weighting frequency measurement method, and the specific operation includes the following steps as shown in fig. 13:
1. dividing x (n) into M groups which do not overlap with each other, wherein the data length is L points, and the sampling data is expressed as si＝[x(iL),x(iL+1),…,x(iL+L)]Wherein i is 0,1, …, M-1, FFT operation of L points is respectively carried out on the M groups of data to obtain an average spectrum
2. After amplitude accumulation, the position corresponding to the maximum amplitude value in L frequency spectrums is searched, and the signal frequency is roughly measured
3. Calculating the precise deviation of the actual frequency from the coarsely measured frequency
Wherein, ∠ SiIs the phase corresponding to the ith maximum magnitude spectrum after FFT.
4. The frequency of the signal is measured accurately,
it is to be understood that the described embodiments are merely a few embodiments of the invention, and not all embodiments. All other embodiments, which can be derived by one of ordinary skill in the art and related arts based on the embodiments of the present invention without any creative effort, shall fall within the protection scope of the present invention.
Claims (10)
1. A radar jamming equipment jamming effect monitoring and decision support system is characterized by comprising
The radar radiation source signal sampling module is used for sampling radar radiation source signals from a signal interface of the radar interference equipment;
the interference signal sampling module is used for sampling an interference signal from a signal interface of the radar interference equipment for sampling;
the signal processing module is used for carrying out time-frequency characteristic analysis on the sampled radar radiation source signals and the sampled interference signals;
the interference effect evaluation module is used for monitoring the interference effect of the interference signal in real time according to the signal similarity, and carrying out quantitative analysis on the rationality of the selection of the interference pattern and the accuracy of the setting of the interference parameters to obtain an interference effect evaluation result; the quantitative analysis process comprises the steps of distinguishing signals from different radar radiation sources through a signal sorting method based on QMF analysis, realizing the matching of a radar radiation signal source and an interference signal through a signal sorting method based on rising edge waveform matching, and finally completing the signal similarity analysis through waveform similarity measurement;
and the interference decision generation module is used for giving a relatively accurate interference pattern and interference parameters on the basis of the interference effect evaluation result, and providing effective technical support for equipment operators to make interference decisions.
2. The system of claim 1, wherein the signal processing is performed by performing frequency domain analysis on the signal using a ZFFT frequency measurement method, and specifically, the system is configured to shift a frequency of the acquired signal X to X', perform low pass filtering and resampling, and perform fast fourier transform.
3. The jamming effect monitoring and decision support system for radar jamming equipment according to claim 1, wherein signal processing employs segmented FFThe T phase weighting frequency measurement method carries out frequency domain analysis on the signal and is specifically operated to divide x (n) into M groups which are not overlapped with each other, the data length is L points, and the sampling data is expressed as si＝[x(iL),x(iL+1),…,x(iL+L)]Wherein i is 0,1, …, M-1, and performing L points FFT operation on the M groups of data to obtain siAfter the amplitude is accumulated, the signal frequency is roughly measured, and then phase weighting is carried out, and the signal frequency is accurately measured.
4. The system for monitoring and supporting decision on interference effect of radar interference equipment according to claim 1, wherein signal processing adopts a WVD-based time-frequency analysis method to perform time-frequency analysis on signals, and the method comprises the specific operations of ① obtaining sampling signals and performing Hilbert transform to obtain analysis signals of the sampling signals, ② determining window type, window length and window moving step, ③ moving a window, calculating instantaneous autocorrelation of data in the window by taking a window center as a reference point, ④ performing fast Fourier transform on the instantaneous autocorrelation, ⑤ moving the window by window moving step, and repeating the steps ③, ④ and ⑤ until all data are processed.
5. The system of claim 1, wherein the signal sorting method based on QMF analysis is used to perform orthogonal wavelet decomposition of signals, the input signal energy is divided into orthogonal components in frequency by a two-channel orthogonal mirror filter bank, the resulting QMF pairs are introduced into a tree structure with a certain number of layers, the "blocks" in each layer have the same area size by waveform decomposition, the output components of each filter are input into the QMF pairs in the next layer, each QMF pair decomposes the input signal waveform into two parts of high frequency component and low frequency component with pi as boundary, and the "blocks" refer to rectangular regions containing energy of basis functions in time-frequency plane.
6. Radar jamming equipment interference effect monitoring and decision making according to claim 5Support system, characterized in that the quadrature mirror filter uses an improved SINC filter having a transfer function of
Wherein N is more than or equal to 2 and less than or equal to N/2, (N-2) N is more than or equal to 2, C is a compression variable, S is a scale variable, N is the number of convolution points, and omega (N) is a Hamming window function for inhibiting Gibbs phenomenon.
7. The system for monitoring the interference effect and supporting the decision for the radar interference equipment according to claim 5, wherein the matching of a radar radiation signal source and an interference signal is realized through a signal sorting method based on the rising edge waveform matching, and the system is specifically operated in such a way that ① extracts a pulse signal envelope and performs curve fitting, ② extracts the rising edge of the pulse signal envelope after the curve fitting, ③ judges the similarity of the rising edge waveforms of the radar radiation source signal and the interference signal by adopting a Huasdorff distance method, and the matching of the radar radiation signal source and the interference signal is completed.
8. The system for monitoring the interference effect and supporting the decision of the radar interference equipment according to claim 7 is characterized in that signal similarity analysis is completed through waveform similarity measurement, and the specific operations are ① for comparing the waveforms of the radar radiation source signal and the interference signal integrally by using an included angle cosine algorithm and a sliding window algorithm to obtain basic waveform similarity, ② for comparing the amplitudes of the waveforms of the radar radiation source signal and the interference signal to obtain amplitude similarity by using an average absolute difference algorithm, and ③ for weighting the basic waveform similarity and the amplitude similarity averagely to obtain the final similarity value of the waveforms of the radar radiation source signal and the interference signal.
9. The system of claim 8, wherein the evaluation criteria of the interference evaluation result is that the final similarity value is greater than 95% for good interference, 80% -95% for normal interference, 65% -80% for critical interference, and less than 65% for no interference.
10. The system of claim 1, wherein the radar jamming equipment is provided with a display module for displaying an aircraft track, a radar received signal, a jamming signal, a current jamming pattern and jamming parameter, a jamming effect evaluation result, an optimal jamming pattern and jamming parameter; the radar receiving signal and the interference signal are visually displayed by adopting a signal rolling method, and the current interference pattern and the interference parameter are selected through a menu bar, a shortcut key or a window; the optimal interference pattern and the interference parameters are obtained and displayed by the interference decision generation module according to the radar type and the interference effect evaluation result.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010323179.XA CN111474524A (en) | 2020-04-22 | 2020-04-22 | Radar interference equipment interference effect monitoring and decision support system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010323179.XA CN111474524A (en) | 2020-04-22 | 2020-04-22 | Radar interference equipment interference effect monitoring and decision support system |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN111474524A true CN111474524A (en) | 2020-07-31 |
Family
ID=71760554
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010323179.XA Pending CN111474524A (en) | 2020-04-22 | 2020-04-22 | Radar interference equipment interference effect monitoring and decision support system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111474524A (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112596034A (en) * | 2020-10-30 | 2021-04-02 | 南京微毫科技有限公司 | Electronic countermeasure electromagnetic spectrum control system based on universal radar and control method thereof |
| CN112924965A (en) * | 2021-02-01 | 2021-06-08 | 桂林电子科技大学 | Clustering coherent superposition-based frequency diversity array radar target imaging method |
| CN113884993A (en) * | 2021-09-03 | 2022-01-04 | 北京航天华腾科技有限公司 | Electronic countermeasure apparatus and electronic countermeasure method |
| CN114114227A (en) * | 2022-01-27 | 2022-03-01 | 安徽京淮健锐电子科技有限公司 | Radar signal modulation type identification method based on dual-channel heterogeneous fusion network |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2014009947A (en) * | 2012-06-27 | 2014-01-20 | Japan Radio Co Ltd | Interference compensation support device |
| CN105911351A (en) * | 2016-04-11 | 2016-08-31 | 湖南科技大学 | DSP-based real time ZFFT method |
| CN106484659A (en) * | 2016-09-29 | 2017-03-08 | 哈尔滨工程大学 | A kind of Generalized Likelihood Ratio line-spectrum detection method based on discrete Fourier transform |
| CN106778610A (en) * | 2016-12-16 | 2017-05-31 | 哈尔滨工程大学 | A kind of intra-pulse modulation recognition methods based on time-frequency image feature |
| CN109116312A (en) * | 2018-07-20 | 2019-01-01 | 电子科技大学 | Radar signal classification method based on QMFB and convolutional neural networks |
| CN109212491A (en) * | 2018-10-17 | 2019-01-15 | 中国航天科工集团八五研究所 | A kind of Emitter Recognition based on radar signal fingerprint characteristic |
| CN109375060A (en) * | 2018-11-12 | 2019-02-22 | 上海金智晟东电力科技有限公司 | A kind of distribution network failure wave-form similarity calculation method |
| CN110618405A (en) * | 2019-10-16 | 2019-12-27 | 中国人民解放军海军大连舰艇学院 | Radar active interference efficiency measuring and calculating method based on interference mechanism and decision-making capability |
-
2020
- 2020-04-22 CN CN202010323179.XA patent/CN111474524A/en active Pending
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2014009947A (en) * | 2012-06-27 | 2014-01-20 | Japan Radio Co Ltd | Interference compensation support device |
| CN105911351A (en) * | 2016-04-11 | 2016-08-31 | 湖南科技大学 | DSP-based real time ZFFT method |
| CN106484659A (en) * | 2016-09-29 | 2017-03-08 | 哈尔滨工程大学 | A kind of Generalized Likelihood Ratio line-spectrum detection method based on discrete Fourier transform |
| CN106778610A (en) * | 2016-12-16 | 2017-05-31 | 哈尔滨工程大学 | A kind of intra-pulse modulation recognition methods based on time-frequency image feature |
| CN109116312A (en) * | 2018-07-20 | 2019-01-01 | 电子科技大学 | Radar signal classification method based on QMFB and convolutional neural networks |
| CN109212491A (en) * | 2018-10-17 | 2019-01-15 | 中国航天科工集团八五研究所 | A kind of Emitter Recognition based on radar signal fingerprint characteristic |
| CN109375060A (en) * | 2018-11-12 | 2019-02-22 | 上海金智晟东电力科技有限公司 | A kind of distribution network failure wave-form similarity calculation method |
| CN110618405A (en) * | 2019-10-16 | 2019-12-27 | 中国人民解放军海军大连舰艇学院 | Radar active interference efficiency measuring and calculating method based on interference mechanism and decision-making capability |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112596034A (en) * | 2020-10-30 | 2021-04-02 | 南京微毫科技有限公司 | Electronic countermeasure electromagnetic spectrum control system based on universal radar and control method thereof |
| CN112924965A (en) * | 2021-02-01 | 2021-06-08 | 桂林电子科技大学 | Clustering coherent superposition-based frequency diversity array radar target imaging method |
| CN113884993A (en) * | 2021-09-03 | 2022-01-04 | 北京航天华腾科技有限公司 | Electronic countermeasure apparatus and electronic countermeasure method |
| CN113884993B (en) * | 2021-09-03 | 2022-08-09 | 北京航天华腾科技有限公司 | Electronic countermeasure apparatus and electronic countermeasure method |
| CN114114227A (en) * | 2022-01-27 | 2022-03-01 | 安徽京淮健锐电子科技有限公司 | Radar signal modulation type identification method based on dual-channel heterogeneous fusion network |
| CN114114227B (en) * | 2022-01-27 | 2022-05-31 | 安徽京淮健锐电子科技有限公司 | Radar signal modulation type identification method based on dual-channel heterogeneous fusion network |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111427018A (en) | Radar interference equipment interference effect evaluation method | |
| CN111474524A (en) | Radar interference equipment interference effect monitoring and decision support system | |
| Duk et al. | Target detection in sea-clutter using stationary wavelet transforms | |
| CN109061589B (en) | Target motion parameter estimation method of random frequency hopping radar | |
| JP4988332B2 (en) | Techniques for incoherent integration of targets with vague velocities. | |
| Liu et al. | Deep learning and recognition of radar jamming based on CNN | |
| Yang et al. | A novel feature extraction method for ship-radiated noise | |
| CN105783974B (en) | A kind of detection of linear FM signal, method for parameter estimation and system | |
| Dong et al. | Doppler frequency estimation by parameterized time-frequency transform and phase compensation technique | |
| CN104655929B (en) | A kind of digital time-frequency measuring method of time-domain signal and corresponding target identification method | |
| CN106534014A (en) | Accurate detection and separation method for multi-component LFM signal | |
| Han et al. | Wideband interference suppression for SAR via instantaneous frequency estimation and regularized time-frequency filtering | |
| Chen et al. | LFM signal detection and estimation based on deep convolutional neural network | |
| Dias et al. | Nonparametric estimation of mean Doppler and spectral width | |
| Cosoli et al. | A real-time and offline quality control methodology for SeaSonde high-frequency radar currents | |
| Hao et al. | Novel method for non-stationary signals via high-concentration time–frequency analysis using SSTFrFT | |
| CN105824017B (en) | A kind of radial velocity acquisition methods and device based on automobile gap-filler radar | |
| CN115345216A (en) | FMCW radar interference elimination method fusing prior information | |
| Shreyamsha Kumar et al. | Target identification using harmonic wavelet based ISAR imaging | |
| CN115712115A (en) | Radar target attribute scattering center extraction method combining calculation and measurement data | |
| Wang et al. | Automated weak signal frequency estimation method based on Duffing oscillator and particle swarm optimization | |
| Yang et al. | LPI radar signal detection based on the combination of FFT and segmented autocorrelation plus PAHT | |
| US5602751A (en) | Method and apparatus for spectrum analysis by complex cumulants | |
| Guo et al. | Pulse-to-pulse periodic signal sorting features and feature extraction in radar emitter pulse sequences | |
| Xiaodong et al. | Interception algorithm of S-cubed signal model in stealth radar equipment |
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 |