Next Article in Journal
Transformerless High Step-Up DC-DC Converters with Switched-Capacitor Network
Previous Article in Journal
Optimized Design of Wind Turbine Blade Receptors Based on Electrostatic Field Theory
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Emitter Signal Waveform Classification Based on Autocorrelation and Time-Frequency Analysis

1
College of Physical Science and Technology, Central China Normal University, No. 152 Luoyu Road, Wuhan 430079, China
2
Department of Electronic Technology, Naval University of Engineering, Wuhan 430033, China
*
Author to whom correspondence should be addressed.
Electronics 2019, 8(12), 1419; https://doi.org/10.3390/electronics8121419
Submission received: 30 October 2019 / Revised: 21 November 2019 / Accepted: 25 November 2019 / Published: 28 November 2019
(This article belongs to the Section Circuit and Signal Processing)

Abstract

:
Emitter signal waveform recognition and classification are necessary survival techniques in electronic warfare systems. The emitters use various techniques for power management and complex intra-pulse modulations, which can create what looks like a noisy signal to an intercept receiver, so emitter signal waveform recognition at a low signal-to-noise ratio (SNR) has gained increased attention. In this study, we propose an autocorrelation feature image construction technique (ACFICT) combined with a convolutional neural network (CNN) to maintain the unique feature of each signal, and a structure optimization for CNN input layer called hybrid model is designed to achieve image enhancement of the signal autocorrelation, which is different from using a single image combined with CNN to complete classification. We demonstrate the performance of ACFICT by comparing feature images generated by different signal pre-processing algorithms, and the evaluation indicators are signal recognition rate, image stability degree, and image restoration degree. This paper simulates six types of the signals by combining ACFICT with three types of hybrid model, the simulation results compared with the literature show that the proposed methods not only has a high universality, but also better adapts to waveform recognition at low SNR environment. When the SNR is –6 dB, the overall recognition rate of the method reaches 88%.

1. Introduction

Electronic warfare (EW) is a military action whose objective is the control of the electromagnetic spectrum (EMS). This objective is achieved through offensive electronic attack (EA), defensive electronic protection (EP), intelligence gathering, and threat recognition electronic warfare support (ES) actions. Electronic intelligence (ELINT) receiver via prolonged and accurate measurement of all the characteristics of a radar emitter (waveform and antenna patterns) in order to provide the necessary data for its analysis and modeling of the associated weapon system, as well as for its identification to be logged in the emitter libraries [1]. In practice, the automatic emitter waveform recognition technique is a core survival technique for an intercept receiver performing threat recognition and radar emitter identification. As the premise and basis of recognition, the emitter signal waveform classification is an important link to ELINT.
In the literature, there have been a few signal waveform recognition techniques that utilize feature extraction techniques and classification techniques to extract features from the intercepted signal, and to classify the intercepted signal based on the extracted features, respectively. To accurately identify emitter signals requires the extraction of features that can represent signals. Time-frequency analysis (TFA) [2] is widely used in emitter signal waveform recognition and classification. Combined with deep learning in the field of computer vision [3], and models of neural network structures [4,5], researchers have obtained better recognition results from the time-frequency feature of signals [6,7,8]. However, some problems restrict such developments in environment monitoring, where signals generated by different types of electromagnetic sources are in many situations noisy, misshaped, or changing in relation to the weather condition, task, and application [9,10]. Time-frequency images (TFI) will be polluted at low SNRs, making signals more difficult to classify. To reduce the influence of noise on TFI requires the signal to be processed by an effective filtering algorithm. Some signal-filtering algorithms have focused on one or two dimensions of signal processing. In one-dimensional space, combined with some traditional theory based on the characteristics of random noise with a zero mean [11,12], such as sampling integrals and digital averages, a digital average method to realize low-SNR signal recognition was proposed [8]. With the recent development of the theories of adaptive filtering, autocorrelation detection, and wavelet transforms, more new noise-reduction algorithms have been proposed. Wavelet transforms and autocorrelation detection were combined to improve the detection ability of weak LFM signals [13]. Time-frequency peak filtering (TFPF) was proposed to suppress random noise [14]. Empirical mode decomposition (EMD) and discrete wavelet transforms were combined to reduce noise [15]. In two-dimensional space, noise-reduction algorithms mainly deal with time-frequency matrices. A stacked convolution denoising automatic encoder (SCDAE) was proposed, in which time-frequency data were reconstructed to enhance the signal component in the TFI [7]. A method of image morphology (IM) erosion and expansion to process time-frequency data was proposed [16]. Singular value decomposition was used to decompose the time-frequency matrices of a signal into noise and signal subspaces to reduce the influence of noise [17]. The threshold value of a time-frequency matrix was defined to remove noise as much as possible through threshold filtering [18]. IM and threshold filtering were combined to reduce noise in TFI [19].
However, most noise-reduction algorithms are difficult to apply at a low-SNR environment, such as digital average denoising [11], it did not give further explanation on how to judge the trigger condition of the digital average. In fact, the literature pointed out that finding digital average trigger conditions is also faced with a complex environment with low SNR [8]. For adaptive filtering [12], the noise component completely submerges the signal component in the time-frequency feature image of the signal at low SNRs. In addition, for the two-dimensional noise reduction algorithm, the overall recognition rate of the de-noised signal in [7] is less than 65% in the SNR of –9 dB. The way of reducing time-frequency domain noise by image processing can reduce the noise of TFI [16,18], [19], however, the processing of IM and threshold filtering tends to lose the signal components contained in TFI. The noise-reduction effect for singular value decomposition was simulated at SNR > 0 dB [17], while simulations were not conducted for SNR < 0 dB. After the observation of some denoising algorithms, it can be found that some were merely applied with a certain style of signal, such as LFM [13], seismic signals [14], and electrocardiogram (ECG) signals [15]. However, considering the application of noise-reduction algorithms to specific signals in unknown space, the premise is still to understand the modulation styles of specific signals.
In view of the high requirement of signal de-noising at low SNR and the practical requirement of emitter signal recognition task, our previous studies used autocorrelation images and multiple parallel CNNs to improve signal recognition rates, but lacked theoretical analysis of the causes of autocorrelation images [18]. Since the classification task are different from signal parameter measurements, the techniques provided in this manuscript are intended to meet the requirements for signal classification at low SNR, but the measurement of signal parameters is not considered. Therefore, we propose a signal autocorrelation feature image construction technique (ACFICT) which a signal’s autocorrelation sequence is acquired, and the sequences are converted to the two-dimensional space after TFA. Unlike the TFI, which is obtained directly from the TFA, a TFI combined with the ACFICT can maintain the unique feature of each signal. This manuscript further explains the construction techniques of autocorrelation images on a theoretical level. The feature images formed by different algorithms are comprehensively evaluated according to the three image evaluation indicators in signal recognition rate, image stability degree, and image restoration degree. Then, unlike taking the single CNN to classify, we combine two groups of CNN and one group of bi-directional long short-term memory (BiLSTM) as classifier. For the convenience of description, the structure is called hybrid model in this manuscript. The three groups of models are not connected in the feature extraction stage. Additionally, the input-layer structure optimization is further analyzed by the signal recognition rate and time consumption. We verified the universality of ACFICT by changing the network structure of CNN in the hybrid model. In comparison with the competition literatures, the simulation results further prove that ACFICT combined with the hybrid model is a more effective signal modulation type classification method.
The remainder of this manuscript is organized as follows. Section 2 introduces the basic framework of emitter signal recognition based on deep learning. Section 3 proposes a pre-processing method for signals and the evaluation indices defined by the signal pre-processing and emitter-signal recognition tasks. Section 4 presents the results of a simulation experiment, and the dataset-generation mode of the experiment is given. Section 5 relates our conclusions.

2. Basic Recognition Framework

In this section, we present the basic framework of signal classification based on deep learning. The fundamental task of signal classification is to classify various signals captured from the space. The methods based on the deep learning gain a better performance currently [11,20,21]. The specific process to signal is shown in Figure 1, the space signal is captured into the signal receiver via the antenna. The intermediate frequency (IF) sampling signals y ( k ) will be pre-processed which can decrease the effect of spatial noise on the recognition task. Additionally, traditional pre-processing may be performed in one-dimensional space, or it may be combined with signal transformation to multi-dimensional space and then pre-processed. The authors of [7,19,21] denoise the signal by analysing TFI in the time-frequency transform domain. In [11], the two pre-processing modes are combined. First, digital averaging is performed on y ( k ) , and then the TFI of the signal is processed in the time-frequency transform domain. The ACFICT proposed in this manuscript also combine the two pre-processing procedures together, where y ( k ) dealt with autocorrelation is converted into TFI that has positive robustness to noise, then the TFI is processed by down sampling to reduce image size and the memory consumption on the graphics processing unit (GPU) and central processing unit (CPU). Finally, the pre-processed signal is input into the pre-designed and pre-trained classifier to classify the signal modulation mode. To stress that this manuscript is similar to [19,22,23,24], and assume that we sample a complete data in one pulse of the signal.

3. Propose Signal Pre-Processing Method

In this section, an image-construction algorithm of the signal is proposed to construct a feature image that can describe a signal with high robustness at a low SNR. Combined with autocorrelation detection and time-frequency transform theory, the algorithm can produce a feature image and achieve image enhancement. To evaluate the performance between the proposed algorithm and other algorithms, considering different signal classification methods based on images, they all have similar requirements for signal feature images, this manuscript proposes three evaluation indices: Stability, recoverability, and recognition rate.

3.1. Signal Pre-Processing

3.1.1. Emitter Signals

The signal from the receiving module can be generally expressed as follows:
y ( k ) = x ( k ) + n ( k ) = A ( k ) e j θ ( k ) + n ( k )
where, x ( k ) is the ideal discrete signal coming after if sampling, n ( k ) is additive white Gaussian noise, k is the index value increasing with sampling interval sequentially, A ( k ) is the instantaneous envelope of the ideal sampling signal, θ ( k ) is the instantaneous phase of the ideal sampling signal. Additionally, instantaneous phase θ ( k ) can be computational expressed further by instantaneous frequency f ( k ) and instantaneous phase offset ϕ ( k ) as follows:
θ ( k ) = 2 π f ( k ) ( k T s ) + ϕ ( k )
where, T s represents the signal sampling interval. In practice, various signals are usually realized by changing the frequency (frequency modulation) and phase (phase modulation) of the signal. In this manuscript, six types of signals are simulated to verify the algorithm which are conventional phase (CP), linear frequency modulation (LFM), nonlinear cosine phase modulation (NCPM), binary phase shift keying (BPSK), binary frequency shift keying (BFSK), and quadrature phase shift keying (QFSK). Some detailed information such as instantaneous envelope, instantaneous frequency, and instantaneous phase offset of the six types above are shown in Table 1.
As shown in Table 1, the carrier frequencies f c , f 1 , f 2 , f 3 , and f 4 are fixed values. In the modulation of LFM, B is the signal bandwidth, τ p w is the signal pulse width, k T s is the discrete sampling time. In the modulation of NCPM, f k is the modulation frequency of the instantaneous phase.

3.1.2. Time-Frequency Transformation

The transformation used in this manuscript is the Choi–Williams distribution (CWD). When applied to signal x ( t ) , the time-frequency transformation with the bilinear form is expressed as
C x ( t , Ω ) = π σ τ 2 x ( u + τ / 2 ) x ( u τ / 2 )   e π 2 σ ( u t ) 2 / 4 τ 2 j Ω τ d u d τ
where C x ( t , Ω ) is the result of TFA, t and Ω are the time and frequency axis, respectively, and the scale factor σ is used to control the distribution of CWD cross-terms. When σ inhibits the cross-terms in CWD, the frequency resolution decreases. The value σ = 1 is used to balance the relationship between CWD cross-term suppression and the frequency resolution.
Six types of signals are simulated with an SNR of 9 dB in this manuscript. The TFIs of the six signals that can be obtained by TFA are shown in Figure 2.

3.1.3. Signal Feature Images Construction

Given a signal x ( t ) and zero-mean Gaussian white noise n ( t ) , the observable signal is expressed as y ( t ) = x ( t ) + n ( t ) , So, the autocorrelation is simplified to
R y ( τ ) = E [ y ( t ) y ( t τ ) ] = E { [ x ( t ) + n ( t ) ] [ x ( t τ ) + n ( t τ ) ] } = E [ x ( t ) x ( t τ ) ] + E [ n ( t ) n ( t τ ) ] + E [ x ( t ) n ( t τ ) ] + E [ n ( t ) x ( t τ ) ] = R x ( τ ) + R n ( τ ) + R x n ( τ ) + R n x ( τ )
Since the noise n ( t ) is not related to the signal x ( t ) , we can conclude that R x n ( τ ) = R n x ( τ ) = 0 , and then R y ( τ ) = R x ( τ ) + R n ( τ ) . For zero-mean Gaussian white noise n ( t ) with wide bandwidth, its autocorrelation function R n ( τ ) mainly affects the nearby position τ = 0 , and when the value of τ is large, then R x ( τ ) can be characterized approximately by R y ( τ ) , that is, R x ( τ ) R y ( τ ) . For convenience, we will write R y ( τ ) as R y ( t ) . Then, combined with the TFA result of the autocorrelation function, R y ( t ) is calculated after being taken into Equation (3):
C R ( t , Ω ) = π σ τ 2 R y ( u + τ / 2 ) R y ( u τ / 2 )   e π 2 σ ( u t ) 2 / 4 τ 2 j Ω τ d u d τ
According to Equation (5), the result obtained by the Choi–Williams transformation reflects the change of the autocorrelation function R y ( t ) . As R y ( t ) is determined only by the signal x ( t ) , the obtained time-frequency transform C R ( t , Ω ) is a response of x ( t ) in the time-frequency plane after the autocorrelation domain, so it can represent the signal x ( t ) uniquely.
In the actual signal reception, the IF sampling of y ( t ) is always the discrete value noted as y ( k ) , and the autocorrelation function of the discrete values is
R ^ y ( k ) = n = 0 N 1 y ( n ) y ( n k )
where N is the number of sampling points, and R ^ y ( k ) is the estimated value of R y ( τ ) . When n < 0 , then y ( n ) = 0 . So, the discrete transformation R ^ y ( k ) has the relation
C R ( i , l ) = m k π σ m 2 R ^ y ( k + m ) R ^ y ( k m )   e π 2 σ ( k i ) 2 / 16 m 2 j 4 π m l / N
Then, C R ( i , l ) will be inverted into the pixel range considering the need to convert C R ( i , l ) to a picture, and the pixel value range of the image is 0 ~ 255.
For convenience of expression, C i j represents the data in the two-dimensional matrix C R ( i , l ) of row i and column j , C R represents C R ( i , l ) , and the equation can be written as
P i c t u r e = 255 C i j min ( C R ) max ( C R ) min ( C R )
where max ( ) and min ( ) respectively represent the maximum and minimum points of the two-dimensional matrix, and P i c t u r e is the feature image of various signals constructed in this manuscript that can represent the image formed by mapping the matrix C R ( i , l ) to the pixel interval. Finally, the feature image is compressed to size 64 × 64 by extracting pixels at equal intervals to reduce the computational load of the CPU and GPU.
Combined with signal autocorrelation and TFA, the autocorrelation images of the six types of signals are obtained as shown in Figure 3 at an SNR of 9 dB. It can be seen that the noise appearing in Figure 2 disappears significantly after the signals have been autocorrelation processed.
According to Equation (6), the autocorrelation value R ^ y ( 0 ) reaches the maximum when k = 0 and decreases as k increases. The case found from the observation of Figure 3 is that the autocorrelation images of the signals have the largest pixel value at the left edge, while the pixel gradually weakens and disappears along the right side, conforming to the distribution rule of autocorrelation values in Equation (6). However, due to the effect of autocorrelation calculation, the problem is that we can’t extract more features from every pixel in the image, while only few areas of the image have useful information. To solve this problem, the value range of k during the construction of signal feature images is analyzed in Section 3.1.4 based on the distribution law of signal autocorrelation values.

3.1.4. Signal Autocorrelation Deviation Analysis

In reality, R ^ y ( k ) is used as an estimate of R y ( τ ) because the signal after sampling is composed of one-dimensional vectors of finite length.
We further discuss the estimation deviation, which occurs when calculating R ^ y ( k ) , as this will affect the pixel distribution of the signal autocorrelation image. We first take the continuous signal y ( t ) as an example to evaluate the estimation deviation R ^ y ( k ) . The calculation is usually taken within the limited time T , so the autocorrelation function can be expressed as
R ^ y ( τ ) = 1 T 0 T y ( t ) y ( t τ ) d t
where R ^ y ( τ ) is the estimated value of R y ( τ ) within the finite time, and R ^ y ( k ) can be obtained after sampling within a time range of integration.
The mean square deviation of R ^ y ( τ ) can be expressed as
var [ R ^ y ( τ ) ] = E [ ( R ^ y ( τ ) R y ( τ ) ) 2 ]
When the zero-mean white noise n ( t ) in the received signal y ( t ) has a Gaussian distribution with bandwidth of B , then the variance of R ^ y ( τ ) can be calculated as [22].
var [ R ^ y ( τ ) ] 1 2 B T [ R ^ y 2 ( 0 ) + R ^ y 2 ( τ ) ]
When R ^ y ( τ ) 0 , then the normalized root mean square deviation of R ^ y ( τ ) can be expressed as
ε ( τ ) = var [ R ^ y ( τ ) ] R ^ y ( τ ) = 1 2 B T 1 + ρ y 2 ( τ ) ρ y ( τ )
After analyzing the deviation between the estimated value R ^ y ( τ ) and the theoretical value R y ( τ ) based on Equation (12), it can be found that the normalized deviation ε ( τ ) is jointly determined by signal bandwidth B , integral time T , and signal autocorrelation coefficient ρ y ( τ ) . For the fact that R ^ y ( k ) can be obtained by sampling and autocorrelation within the integral time, so the deviation of R ^ y ( k ) should have the expression similar to Equation (12). In fact, according to the Nyquist sampling theorem, the sampling rate of the signal should be greater than twice its maximum angular frequency, that is, f s 2 B , where f s is the signal sampling rate. The number of sampling points N of the integration time T must obey the following expression to maintain information integrity in the continuous-to-discrete transformation, where T s is the sampling interval time:
N = T T s = f s T 2 B T
When f s obeys the Nyquist Sampling theorem, the deviation of R ^ y ( k ) (i.e., ε ( k ) ) is calculated as
ε ( k ) = 1 N 1 + ρ y 2 ( k ) ρ y ( k )
where ρ y ( k ) = R ^ y ( k ) / R ^ y ( 0 ) . The deviation of R ^ y ( k ) can be estimated reasonably using Equation (14). Different deviations ε ( k ) can be obtained based on the autocorrelation coefficient ρ y ( k ) and number of sampling points N . The value of ε ( k ) changes with k . Six types of signals are simulated under parameters satisfying Equation (13), with an SNR of 9 dB, to analyze the value and change trend of ε ( k ) . The change curves of ε ( k ) are shown in Figure 4.
The number of signal sampling points in Figure 4 is set to 256. It can be found that there are relatively mild changes in ε ( k ) at k 125 , while there will be dramatic deviations in R ^ y ( k ) , which fail to be used to construct the feature images at k > 125 . In practice, when k > 125 , the estimation deviation for R ^ y ( k ) is large because it is calculated from few discrete values, which is the reason why the pixels representing the signal in the autocorrelation images of Figure 3 gradually weaken until they disappear.
In order to obtain more abundant autocorrelation values which can form autocorrelation feature images, we should explore the most appropriate range of k in R ^ y ( k ) . So, the calculation method for signal’s autocorrelation result is changed to extend the range of R ^ y ( k ) . Additionally, H ^ y ( k ) is calculated as
H ^ y ( k ) = { n = 0 N 1 y ( n ) y ( n k ) k 0 H ^ y ( k ) k < 0
R ^ y ( k ) is calculated from H ^ y ( k ) as
R ^ y ( k ) = H ^ y ( k N ) k = 1 , 2 , , 2 N 1
The deviation curves for each signal are updated based on Equation (16) as shown in Figure 5.
The number of sampling points in Figure 5 is 511, and it can be found that the deviation value changes slightly when 150 k 350 , but beyond that range, the estimated values R ^ y ( k ) are calculated from fewer discrete values according to Equations (15) and (16), and thus, the dramatic change will occur again.
After comparing Figure 4 with Figure 5, it can be found that the range 150 k 350 in Figure 5 has more reasonable autocorrelation values, so it is more suitable for the construction of the signal autocorrelation feature image. We can find that when R ^ y ( k ) appears to be an unstable change, a large deviation value will appear at the position far from k = 0 , and the autocorrelation deviation is relatively small near k = 0 , so this manuscript can further determine the reasonable range of k based on these factors.
For a given sampling signal y ( k ) with length N , an autocorrelation sequence R ^ y ( k ) with length N is formed from the range ( N / 2 , 3 N / 2 ) of autocorrelation sequence R ^ y ( k ) with length 2 N 1 , acquired from Equation (16). For convenience, the function c l i p ( ) is used to represent the generation method for the above autocorrelation sequence.
The deviation ε ( k ) is constrained in a smaller range and, in turn, the autocorrelation image reconstructed by R ^ y ( k ) can more effectively characterize the signal.
We modify the autocorrelation results based on the analysis in this section, and R ^ y ( k ) is taken into Equation (6) to construct an autocorrelation feature image with high robustness toward the noise effect. The autocorrelation feature image formed by R ^ y ( k ) at an SNR of 9 dB is shown in Figure 6.
Comparing Figure 3 and Figure 6, it can be seen that extending the range of k is efficient for the images in Figure 6, the images have better modification in resolution than those in Figure 3.
In this study, the autocorrelation sequence is formed by intercepting the autocorrelation value of the interval ( N / 2 , 3 N / 2 ) , which avoids the part of the autocorrelation sequence with large deviation values, so that autocorrelation images can fully reflect the overall change of the signal after the autocorrelation domain, and at the same time, remove the invalid part of the image caused by the larger deviation.

3.1.5. Signal Feature Image Enhancement

In this manuscript, the feature image based on autocorrelation theory can effectively represent six types of signals, but its image is still affected by low SNR. Specifically, the pixel strength of the image begins to weaken at low SNR. Figure 7 shows the changes of autocorrelation images of NCPM signals in different SNR environments. It can be found that with the decrease of SNR, the image of NCPM tends to blur. Therefore, in order to obtain feature images that can adapt to low SNR environments, this section is devoted to image enhancement of autocorrelation images of signals.
In order to achieve image enhancement, the pixels representing the signal should have larger pixel values. The pixel values in this manuscript are the reflection of the signal in time-frequency plane after passing through the autocorrelation domain. If we want to enhance the pixels representing the signal, we should consider improving the autocorrelation values which can represent the signal. Considering that the iteration calculation of autocorrelation can improve the autocorrelation value of signals, this manuscript proposes to implement image enhancement by multiple autocorrelation of signals. For convenience, the result of m times autocorrelation is expressed as H ^ y m ( k ) and R ^ y m ( k ) , the length of R ^ y m ( k ) is noted as N m , and the primitive signal y ( k ) is noted as y 0 ( k ) , the length of y 0 ( k ) is noted as N 0 . If R ^ y m ( k ) is considered as a new signal y m ( k ) , the autocorrelation result of the sequence y m ( k ) is given below.
H ^ y m m + 1 ( k ) = { n = 0 N m 1 y m ( n ) y m ( n k ) k 0 H ^ y m m + 1 ( k ) k < 0
R ^ y m m + 1 ( k ) = H ^ y m m + 1 ( k N m ) k = 1 , 2 , , 2 N m 1
y m + 1 ( k ) = R ^ y m m + 1 ( k )
Through the iterative calculation of Equations (17)~(19), the autocorrelation result R ^ y m ( k ) of the signal after m times autocorrelation can be obtained. However, considering the change of the length of the autocorrelation sequence by Equation (18), iterative calculation will make the sequence length rapid grow, so this section also uses the c l i p ( ) function proposed in Section 3.1.4 to process multiple autocorrelation sequences, so that the length of multiple autocorrelation is also N . Based on the image-enhancement theory in this study, the two times autocorrelation image of the signal is shown in Figure 8. Comparing this with Figure 7, it can be found that the feature image of the signal obtained through two times autocorrelation can extract more features from every pixel in the image.

3.2. Model Input Layer Optimization

According to Section 3.1.5, multiple autocorrelation can achieve image enhancement of the signal autocorrelation image, and each autocorrelation can obtain an autocorrelation image that characterizes the signal. In fact, three autocorrelation images characterizing the signal can be obtained by TFA after three times autocorrelation of the signal. More primitive information can be provided if the underlying information of autocorrelation image is used wisely, so we propose an optimized usage for the multiple images as an input of the CNN classification.
As shown in Figure 9, the current TFA-based classification network only takes a single channel TFI as input [7,11,19,21,23], because TFI acquired by TFA is usually a single channel image. However, this manuscript obtains three feature images which can represent signals by iterating the autocorrelation value of signals, which is different from using a single image combined with CNN to complete classification. It means that this method provides more initial information, but the network structure also needs to be further improved. As a result, a deep learning network is designed as shown in Figure 9. In this network, we combine two groups of CNN and one group of BiLSTM to realize the feature extraction of three feature images, and finally complete the signal classification. The three groups of models are not connected in the feature extraction stage. Only after extracting the high-dimensional features, we use the full connection layer to synthesize all the neuron information of each model for classification. For ease of explanation, the network structure is called the hybrid model.
It should be emphasized that considering the different hyper parameters of CNN and the depth of the network layer have different effects on the classification results, this manuscript discusses the CNN model used in the classification of signal modulation methods [11,19]. At the same time, a deeper CNN structure is proposed to analyze the impact of different networks on the classification of signals, which will be further explained in Section 4.2.2. For the BiLSTM structure. Zachary [25] points out that if the forgetting gate bias is set to one during model initialization, the output explosion or disappearance problem caused by gradient in the process of model transmission at the beginning can be effectively avoided, so this manuscript sets it to one.

3.3. Performance Evaluation of Algorithm

The performance of the algorithm in this manuscript is mainly evaluated according to the signal recognition rate, image stability degree, and image restoration degree. The image stability degree and image restoration degree are indices to evaluate different pre-processing algorithms, which are defined in this section. Aditionally, the signal recognition rate is the final evaluation index of the signal classification task. This manuscript compares it with similar competitive literature.
For the image stability degree, y ( k ) is the signal to be processed, and feature images A and B are formed from the effects of noise interference A and B , respectively. When feature images A and B have high similarity, the given pre-processing algorithm can be proven to be slightly influenced by different noises, i.e., it has high stability.
For the image restoration degree, feature images A and C in Figure 10 are respectively formed with and without noise interference. If A and C have high similarity, it shows that the signal pre-processing algorithm can effectively suppress the influence of noise and restore the image formed by the signal without noise.
We use (peak signal-to-noise-ratio) PSNR, which is the logarithmic expression of mean square error, to evaluate the similarity of two images.
For convenience, in feature images P and Q of size m × n , the pixel points in row i and column j are noted as P ( i , j ) and Q ( i , j ) , respectively, and the mean square error is calculated as
M S E = 1 m n i = 1 m j = 1 n | | P ( i , j ) Q ( i , j ) | | 2
PSNR is calculated as
PSNR = 10 log 10 ( 255 2 M S E )
According to Equations (20) and (21), the higher the value of PSNR the higher the similarity of the two images. Therefore, a higher value of PSNR indicates a greater image stability degree and image restoration degree.

4. Evaluation and Analysis of ACFICT

In this section, the performance of ACFICT is evaluated. We first explain the generation method of the dataset in the simulation experiment. Then, we introduce common pre-processing methods and classification networks. Finally, the proposed algorithm is contrasted with the given pre-processing methods and the classification networks. The performance of ACFICT is evaluated in terms of image stability degree, image restoration degree, signal recognition rate, and time consumption. Additionally, the input-layer structure optimization is further analyzed by the signal recognition rate and time consumption.

4.1. Simulation Environment and Data Set Generation

In a real environment, the carrier frequency and signal parameters of various types of signals change all the time. Hence, the values of the parameters of the simulated signals are randomly selected within the specified range in order to keep the simulation consistent with the actual electromagnetic environment.
For a given range ( a , b ) , U ( a , b ) is used to indicate that the corresponding parameter values are uniformly distributed in the interval from a to b. The parameter value ranges of the six types of emitter signal signals are shown in Table 2 by using the note U ( a , b ) and the parameters in Table 1.
In Table 2, f c , f 1 , f 2 , f 3 , and f 4 are different carrier frequencies; τ p w is the pulse width; f k is the modulation frequency of the NCPM signal; L c is the code length of the BPSK; N is the number of sampling points in the signal; and L is the code length. Within the encoding cycle of a BFSK, two-frequency conversion of a BFSK signal depends on the parameter N _ r a n f , where g ( r a n d ( 1 , 50 ) ) represents that 50 random numbers are generated between 0 and 1 and divided by the boundary of 0.5. Similarly, in the encoding cycle of a QFSK, four-frequency conversion of a QFSK signal depends on the parameter N _ r a n f , h ( r a n d ( 1 , 50 ) ) represents that 50 random numbers are generated between 0 and 1 and divided by the boundary of 0.25, 0.5, and 0.75.
The training and test sets are generated separately, based on the parameter configuration of the six types of signals. The data of the training set are generated from SNRs of –9 to 9 dB, in steps of 3 dB. Each SNR environment has 500 training samples for each type of signal, for a total of 21,000 training samples. Similarly, the data of the test set are generated from SNRs of –9 to 6 dB in 3 dB steps. Every SNR has 300 testing samples for each type of signal, for a total of 12,600 testing samples. The model was built in the Python Tensorflow framework with the Nvidia 1050Ti GPU.

4.2. Signal Pre-Processing Algorithm and Classification Network Model

4.2.1. Signal Pre-Processing Algorithm

To reasonably compare the performance differences between various algorithms, the signal pre-processing algorithms commonly used in signal classification are presented. Since the CNN processes two-dimensional data, the CWD transformation mentioned in Section 3.1.2 was applied to various signal pre-processing algorithms to convert one-dimensional data to two-dimensional data to meet the input data requirements of the CNN.
Common signal pre-processing algorithms can be categorized as having one or two processing dimensions. Commonly used one-dimensional algorithms are the adaptive filter [8,26] and wavelet transform (WT) [27]. The literature about [28] in comparison of LMS adaptive filtering algorithm, points out that the comprehensive performance of the normalized least mean square (NLMS) algorithm is the best. The ACFICT proposed in this manuscript is mainly processed at the one-dimensional level. After the signal is processed by the one-dimensional signal processing algorithm, the data are mapped to the pixel interval to form a feature image that can represent the signal. On the two-dimensional level, the original signal is directly transformed into two-dimensional data by CWD, and the two-dimensional signal processing algorithm is processed on two-dimensional data. The commonly used signal processing algorithm is image threshold denoising (ITD) [19]. Literature [20] proposed combining image morphology with ITD, and the two-dimensional data obtained by signal processing algorithm are mapped to the pixel interval to become the feature image of the signal.
To sum up, we present the signal pre-processing algorithms in Table 3, including LMS, NLMS, wavelet threshold denoising, and the proposed ACFICT at the one-dimensional level. The two-dimensional level includes image threshold denoising and combining ITD with IM. In addition, since CWD can be used as a signal pre-processing algorithm, we consider its performance when used alone.

4.2.2. Network Model

Among the different deep-learning classifiers, CNN has been widely used in [7,20,21,23]. Since the performance of the proposed signal pre-processing algorithm on each network model needs to be evaluated reasonably, in order to express the structure of the network model easily, the structure of CNN is modelled in this manuscript. It performs feature extraction and dimensionality reduction of images through alternate convolution and pooling operations, and its classical network structure usually consists of five parts: The input layer, the convolution layer, (i.e., down sampling layer), the fully connected layer, and the output layer.
As Ian G has already given the specific calculation methods of convolution, pooling, and full connection [28], this manuscript will not go into details about their calculation.
The performance of the algorithm in this manuscript was evaluated with three signal-classification networks. In the network model of signal classification, literature [20] and literature [11] both designed the corresponding network structure and achieved a good recognition rate through super-parameter tuning. The network structure is shown in Figure 11 and Figure 12. Considering the development of the current deep network model, we designed the network structure shown in Figure 13. Among the three networks, the network structure in reference [11] (Figure 12) has more convolution kernels and feature maps than Figure 11). The network presented in this study (Figure 13) has deeper network layers than the others. By designing three network structures, we comprehensively evaluated the performance of the algorithm on different network classifiers. It should be emphasized that in the three networks, the unified input image size of this manuscript is 64 × 64 , the convolution fill is uniformly zero padded, the block size is set to 256, and the learning rate is set to 0.01. For convenience, the network structure of Figure 11 is called CNN1, the network structure of Figure 12 is called CNN2, and the network structure of Figure 13 is called CNN3.
Let S i j n be a symbol to describe i -th pooling layer of the j -th feature image in size n × n ; Let C i k n be a symbol to describe i -th convolution layer of the j -th feature image in size n × n ; Let W i k j l be a symbol to describe the j -th convolution kernel of the k -th group of convolution kernel in i -th convolution layer, the convolution kernel size is l × l , accordingly. Since the described image in the RGB mode, each input image has three channels. So, when i = 0 , then the value of j can be 1, 2, and 3.
Convolution layers and pooling layers are operated alternatively when classifying the image. In the pooling operation, the k -th feature image of the ( i + 1 ) -th pooling layer can be expressed as follows:
C ( i + 1 ) k n = j ( S i j n W ( i + 1 ) k j l ) + b k
where, “ ” represents the convolution operation. b k is a bias of the k -th convolution kernel.
Then, the ( i + 1 ) -th convolution layer goes into the pooling operation:
S ( i + 2 ) k n / 2 = P o o l i n g ( C ( i + 1 ) k n )
where, P o o l i n g ( ) is a pooling function of CNN which can do pooling operation with the n × n feature image in the pooling layer.
In this manuscript, the size of the pooled size is defined as 2 × 2 , the pooling mode is the maximum pooling. After maximum pooling, the number of feature maps remains the same, but the size shrinks into n / 2 × n / 2 .
After a successive of convolution and pooling operations, a series of feature maps which is converted into a unidimensional vector will be put into a fully connection network to form classification results as follows:
[ y 1 , y 2 , , y m 1 , y m ] = F u l l c o n n e c t ( S i j n )
where, F u l l c o n n e c t ( ) represents a fully connected function.
In CNN1, it includes a two-layer convolutional layer and a two-layer pooling layer. According to the parameters of CNN1, it can be modeled as
{ C 1 q 1 64 = q 0 = 1 3 ( S 0 q 0 64 W 1 q 1 q 0 5 ) + b q 1 S 2 q 1 32 = P o o l i n g ( C 1 q 1 64 ) q 1 = 1 , 2 , , 6 C 3 q 2 32 = q 1 = 1 6 ( S 2 q 1 32 W 3 q 2 q 1 2 ) + b q 2 S 4 q 2 16 = P o o l i n g ( C 3 q 2 32 ) q 2 = 1 , 2 , , 16 [ y 1 , y 2 , , y m 1 , y m ] = F u l l c o n n e c t ( S 4 q 2 16 )
The model of CNN2, with two convolution layers and two pooling layers, is expressed as
{ C 1 q 1 64 = q 0 = 1 3 ( S 0 q 0 64 W 1 q 1 q 0 7 ) + b q 1 S 2 q 1 32 = P o o l i n g ( C 1 q 1 64 ) q 1 = 1 , 2 , , 30 C 3 q 2 32 = q 1 = 1 30 ( S 2 q 1 32 W 3 q 2 q 1 3 ) + b q 2 S 4 q 2 16 = P o o l i n g ( C 3 q 2 32 ) q 2 = 1 , 2 , , 60 [ y 1 , y 2 , , y m 1 , y m ] = F u l l c o n n e c t (   S 4 q 2 16 )
The model of CNN3, with three convolution layers and three pooling layers, is expressed as
{ C 1 q 1 64 = q 0 = 1 3 ( S 0 q 0 64 W 1 q 1 q 0 3 ) + b q 1 S 2 q 1 32 = P o o l i n g ( C 1 q 1 64 ) q 1 = 1 , 2 , , 32 C 3 q 2 32 = q 1 = 1 32 ( S 2 q 1 32 W 3 q 2 q 1 3 ) + b q 2 S 4 q 2 16 = P o o l i n g ( C 3 q 2 32 ) q 2 = 1 , 2 , , 64 C 5 q 3 16 = q 2 = 1 64 ( S 4 q 2 16 W 5 q 3 q 2 3 ) + b q 3 S 6 q 3 8 = P o o l i n g ( C 5 q 3 16 ) q 3 = 1 , 2 , , 64 C 7 q 4 8 = q 3 = 1 64 ( S 6 q 3 8 W 7 q 4 q 3 3 ) + b q 4 S 8 q 4 4 = P o o l i n g ( C 7 q 4 8 ) q 4 = 1 , 2 , , 64 [ y 1 , y 2 , , y m 1 , y m ] = F u l l c o n n e c t ( S 8 q 4 4 )
We must emphasize that the number q 0 of channels of the input layer image varies with the algorithm, and in the input layer of the algorithm proposed in this manuscript, q 0 = 1 or q 0 = 3 , while the other algorithms use q 0 = 1 .

4.3. Comparison of Algorithm Performance

4.3.1. Image Restoration Degree and Image Stability Degree

According to the evaluation formula defined in Section 3.3, this manuscript first simulates the signal processing algorithms given in Table 3, and evaluates the image stability degree and image restoration degree of each algorithm. Figure 14 shows the performance of each algorithm on two indices.
From Figure 14, it can be seen that with the increase of the SNR, the PSNR values of each indice are improved.
The PSNR in Figure 14a was used to analyze the image restoration degree of each algorithm. Based on the definition in Section 3.3, the image restoration degree describes the difference between the feature image formed by the pre-processing algorithm with and without noise interference, so the indice can be used to evaluate the noise-suppression effect of each algorithm in different noise environments. A larger PSNR indicates better noise suppression.
In Figure 14a, we can see that WT and ACFICT have more advantages in image restoration degree than the other five algorithms. For WT and ACFICT, when SNR = –9 dB, PSNR = –21 dB in WT, and PSNR = –17 dB in ACFICT. When SNR = 9 dB, PSNR = –8 dB in WT, and PSNR = 9 dB in ACFICT. For the other five algorithms, including LMS, NLMS, CWD, and two algorithms based on ITD, when SNR rises from –9 to 9 dB, their PSNR rises from –25~–20 to –15~–10 dB. It can be found that when the SNR is –9 dB, ACFICT still outperforms the PSNR value of 4 dB than WT algorithm. Therefore, by synthesizing the above analysis, ACFICT achieves the best result in image restoration degree compared with the other six algorithms.
As shown in Figure 14b, the image stability degree was analyzed by PSNR. According to the definition of Section 3.3, image stability degree represents the difference between the two images formed by signal pre-processing algorithm under two independent noise disturbances. Therefore, this indice can be used to evaluate the stability of image features formed by the algorithm in different noise environments. The more stable each pixel characteristic of the signal is, the stronger the anti-noise ability of the algorithm is, and the bigger the PSNR value is.
Adaptive filtering (NLMS, LMS), CWD, and the method combining CWD and ITD have poor stability and share trends like that in Figure 14b. When SNR = –9 dB, then the range of their PSNR is –23~–21 dB; when SNR = 0 dB, their PSNR is –19~–17 dB; and when SNR = 9 dB, their PSNR is –14~–11 dB. WT, IM, and ACFICT have greater stability than the four algorithms above.
It can be concluded that ACFICT has greater image stability than the other six algorithms because the values of PSNR stay higher than 0 dB in the range –9~9 dB. When SNR = –9 dB, then PSNR = 1 dB, and when SNR = 9 dB, then PSNR = 11 dB. However, observing the ACFICT curve, it can be found that its PSNR value is relatively low and the change is small at a SNR of –9~–3 dB. In fact, under the SNR of –9~–3 dB, the noise almost drowns the signal, and the characteristic image composed of ACFICT gradually loses the pixel characteristics of the signal. At this time, it focuses more on reflecting the change after noise autocorrelation, so when SNR < –3 dB, its PSNR value hardly changes. When SNR > –3 dB, the feature image composed of ACFICT starts to focus more on reflecting the change after signal autocorrelation, which makes the PSNR value of ACFICT start to rise.
In summary, from the comprehensive consideration of the image restoration degree and image stability degree, it can be concluded that ACFICT is a more effective signal feature image construction algorithm.

4.3.2. Signal Recognition Rate

Test sets generated according to the simulation environment in Section 4.1 were used in the three networks to verify the seven pre-processing algorithms. The simulation results of the signal recognition rate are shown in Figure 15.
From Figure 15, it can be noted that when SNR = –9 dB, the recognition rates of ACFICT all surpass 60%, whereas ACFICT reaches 60% in CNN1 and CNN2 and 70% in CNN3, while the other six algorithms fall below 50%. In the six algorithms, CWD and the method combining CWD and ITD perform better [7,15,19], all reaching 40% in CNN1 and CNN2, and almost reaching 60% in CNN3, while the other four algorithms fall below 40%.
With the increase of SNR, there are some improvements to all the algorithms. When SNR > 3 dB, except for WT, the recognition rates of the other algorithms are higher, but the recognition rate of each algorithm has a little different in different models. ACFICT has the highest recognition rate in the simple model CNN1. ACFICT, CWD and the method combined with CWD and ITD reach approximately 100% in the CNN2 which has larger convolution kernel and more feature images. ACFICT, CWD, CWD+ITD, and CWD+ITD+IM reach approximately 100% in the CNN3 which has deeper layer.
It can be found that although the recognition rate of each algorithm is different in different models, the trend of change is the same as a whole. At the same time, under the three network models, the ACFICT proposed in this manuscript has better recognition performance, which further illustrates that the better recognition performance of ACFICT is not limited to specific network models. When we use CNN to identify signal modulation mode, the ACFICT is a better signal pre-processing algorithm.

4.3.3. Competitive Literature Comparison

Test sets generated according to the simulation environment in Section 4.1 were used in the comparison with the same type of competition literature. It is noted that reference [7] also studies the classification of modulation types of signals, which explores the recognition rates of random projections and sparse classification (RPSC) [29] and SCDAEs [7], and its simulation signals and research ideas are similar to this manuscript. Therefore, the recognition rate of this manuscript is further compared with that of the technology in [7]. After setting the same simulation conditions, the simulation results of the signal recognition rate are shown in Figure 16.
It can be found in Figure 16 that the recognition rate of the hybrid model in CP, LFM, NCPM, QFSK, and BPSK is better than that of RPSC and SCDAEs, regardless of which CNN structure is adopted. So, from the overall perspective, the hybrid model gets a higher recognition rate.
It is noted that when the SNR is –6 dB, the recognition rate obtained by constructing the hybrid model with CNN3 is higher than that of SCDAEs and RPSC. Overall, when SNR = –6 dB, the recognition rate of CNN3 is 88%, and SCDAEs is 74%. With the increase of SNR, the recognition rate of hybrid model and SCDAEs gradually increases. However, the recognition rate of SCDAEs is comparable to that of the hybrid model only when the SNR is > 6 dB. This further illustrates that the hybrid model proposed in this manuscript has better recognition performance at low SNR on the recognition rate indicator.
It should not be overlooked that different CNN structures will affect the recognition result of the signal modulation mode, which is particularly evident in Figure 16. Among the three hybrid models, CNN3 with deeper layer has the best result, while CNN2 with large convolution kernel has the lowest recognition rate. This shows that although a better network structure has been obtained through hyper parameter adjustment in literature [11] and [20], it is still only applicable to the corresponding signal pre-processing method. In the pre-processing method proposed in this manuscript, a better network parameter configuration is shown in Figure 13.

4.4. Optimization Analysis of Proposing Method

4.4.1. Influence of Changing Input Layer on Signal Recognition Rate

In Section 3.2, based on ACFICT, an optimization scheme for the image of the CNN’s input layer was proposed. Specifically, ACFICT was used to obtain three feature images that can represent signals, and then three images are input into three parallel deep learning models to complete the classification task. Figure 17 shows the recognition results based on a set of and three sets of parallel deep learning models. For a more intuitive comparison, we also draw the recognition results of literature [16] into them.
It can be seen from Figure 17 that based on ACFICT, the recognition result obtained by using a single set of CNN is higher than that proposed by the literature [7]. In the comparison of CNN, no matter which structure is adopted, the hybrid model can obtain better recognition rate improvement than the original CNN structure. This shows that the improvement of the recognition rate of the hybrid model is not due to a specific deep learning model, which can better integrate multiple input images to obtain better recognition results. Specifically, when CNN1 is adopted, the hybrid model increases the recognition result of the original CNN1 by 7%, and for CNN2 and CNN3, it increases by 2% and 3%, respectively. As the SNR increases, the recognition rates of CNN and hybrid models tend to be consistent at SNR = 3 dB.

4.4.2. Recognition Speed Analysis

Table 4 and Table 5 show the processing speeds obtained when changing the number of input feature images and the classification model. The GPU for processing is GTX 1050Ti and the CPU is i7-8565U. It should be emphasized that the simulation results are the average values obtained after 100 Monte Carlo simulations.
It can be found that the time consumption of the method proposed in this manuscript is mainly in the signal pre-processing stage. Although more input information is introduced into multiple autocorrelation feature maps, it inevitably leads to the increase of calculation time. For each additional autocorrelation feature map, the pre-processing speed of ACFICT will increase by 0.02~0.03 s. For the recognition of the classifier, its processing speed is much faster than that of the signal pre-processing stage. However, it cannot be ignored that the high recognition results obtained by using multiple groups of deep learning models in this manuscript also lead to a reduction of one order of magnitude in the recognition speed of the classifier. In terms of the recognition speed, no matter which CNN structure is used, the hybrid model is 10 times lower than the single group CNN. Therefore, according to the method proposed in this manuscript, a conclusion reflecting the simulation results of this manuscript is that ACFICT proposed in this manuscript is an effective signal pre-processing algorithm, which provides a new feature image for signal classification under low SNR, and further improves the signal recognition rate through the hybrid model. However, the average time obtained in signal pre-processing and classifier recognition reveals that ACFICT and hybrid model bear more computing load. Increasing the number of autocorrelation feature images and parallel deep learning models in hybrid model will reduce the recognition speed of the proposed method.

5. Conclusions

ACFICT was proposed to construct signal feature images that are robust enough to represent signals at low SNR. According to the results of deviation analysis, we propose to combine the autocorrelation sequences with TFA in a specified range to form TFI. As for the decrease of pixel intensity of feature images, which may be caused by SNR reduction, we further present a multiple-autocorrelation method to overcome the influence of low-SNR signal feature images and proposed the optimization scheme of an input-layer image based on multiple-autocorrelation theory. In the testing stage, ACFICT is compared with six types of signal pre-processing algorithms on three indices, which are image restoration degree, image stability degree, and signal recognition rate. At the same time, the input layer optimization scheme based on ACFICT can further improve the signal recognition rate at low SNR. Three images that generated by ACFICT are input into three parallel deep learning models called hybrid model. In the comparison of CNN, no matter which structure is adopted, the hybrid model can obtain better recognition rate improvement than the original CNN structure. This shows that the improvement of the recognition rate of the hybrid model is not due to a specific deep learning model, which can better integrate multiple input images to obtain better recognition results. In comparison with the competition literatures, the simulation results further prove that ACFICT combined with the hybrid model is a more effective signal modulation type classification method. When the SNR is –6 dB, the overall recognition rate of the method reaches 88%. So, the hybrid model combined with the ACFICT can maintain the unique feature of each signal and is not susceptible to noise interference at the low-SNR environment, which meets the requirements of signal classification and is easier to classify signals.

Author Contributions

Guidance of theoretical analysis and writing, Z.M. and Z.H.; Operation of the experiments, analysis, and writing of the manuscript, Z.M.; Guidance of theoretical analysis, G.H.; Guidance and optimizing of experiments, A.L.

Funding

This work was supported in part by the National Natural Science Foundation of China under Grant 61501484.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. De Martino, A. Introduction to EW Systems, 2nd ed.; Artech House: Norwood, MA, USA, 2012. [Google Scholar]
  2. Zhu, J.D.; Zhao, Y.J.; Tang, J. Automatic recognition of radar signals based on time-frequency image character. In Proceedings of the IET International Radar Conference, 14–16 April 2013; pp. 1–6. [Google Scholar]
  3. Krizhevsky, A.; Sutskever, I.; Hinton, G.E. ImageNet Classification with Deep Convolutional Neural Networks. In Proceedings of the Neural Information Processing Systems (NIPS), Lake Tahoe, NV, USA, 3–6 December 2012; pp. 1097–1105. [Google Scholar]
  4. Matuszewski, J. Radar signal identification using a neural network and pattern recognition methods. In Proceedings of the 14th International Conference on Advanced Trends in Radioelectronics, Slavske, Ukraine, 20–24 February 2018; pp. 79–83. [Google Scholar]
  5. Matuszewski, J.; Huang, D.P. Recognition of electromagnetic sources with the use of deep neural networks. In Proceedings of the 12th Conference on Reconnaissance and Electronic Warfare Systems (CREWS), Oltarzew, Poland, 19–21 December 2018; pp. 98653–98667. [Google Scholar]
  6. Zhou, Z.; Huang, G.; Gao, J. Radar Emitter Recognition algorithm based on deep learning. J. Xidian Univ. 2017, 44, 77–82. [Google Scholar]
  7. Zhou, Z.; Huang, G.; Chen, H. Automatic Radar waveform recognition Based on Deep Convolutional Denoising Auto-encoders. Circuits Syst. Signal Process. 2018, 37, 4034–4048. [Google Scholar] [CrossRef]
  8. Kong, S.; Kim, M.; Hoang, L.M. Automatic LPI Radar waveform recognition Using CNN. IEEE Access 2018, 6, 4207–4219. [Google Scholar] [CrossRef]
  9. Matuszewski, J.; Sikorska-Łukasiewicz, K. Neural network application for emitter identification. In Proceedings of the International Radar Symposium, Prague, Czech Republic, 28–30 June 2017; pp. 28–30. [Google Scholar]
  10. Janusz Dudczyk, A.K. Adaptive Forming of the Beam Pattern of Microstrip Antenna with the Use of an Artificial Neural Network. Int. J. Antennas Propag. 2012, 2012, 935073. [Google Scholar]
  11. Gao, J.Z. Detection of Weak Signals; Tsinghua University Press: Beijing, China, 2004. [Google Scholar]
  12. Bismor, D.; Czyz, K.; Ogonowski, Z. Review and Comparison of Variable Step-Size LMS Algorithms. Int. J. Acoust. Vib. 2016, 21, 24–39. [Google Scholar] [CrossRef]
  13. Le, B.; Liu, Z.; Gu, T. Weak LFM signal detection based on wavelet transform modulus maxima denoising and other techniques. Int. J. Wavel. Multiresolution Inf. Process. 2010, 8, 313–326. [Google Scholar] [CrossRef]
  14. Li, J.; Meng, K.; Li, Y. Adaptive linear TFPF for seismic random noise attenuation. J. Pet. Explor. Prod. Technol. 2018, 8, 1443–1453. [Google Scholar] [CrossRef]
  15. Kabir, M.A.; Shahnaz, C. Denoising of ECG signals based on noise reduction algorithms in EMD and wavelet domains. Biomed. Signal Process. Control. 2012, 7, 481–489. [Google Scholar] [CrossRef]
  16. Zilberman, E.R.; Pace, P.E. Autonomous Time-Frequency Morphological Feature Extraction Algorithm for LPI Radar Modulation Classification. In Proceedings of the International Conference on Image Processing, Atlanta, GA, USA, 8–11 October 2006; pp. 2321–2324. [Google Scholar]
  17. Hassanpour, H. A Time-Frequency Approach for Noise Reduction; Academic Press: New York, NY, USA, 2008. [Google Scholar]
  18. Huang, Z.; Ma, Z.Y.; Huang, G.M. Radar Waveform Recognition Based on Multiple Autocorrelation Images. IEEE Access 2019, 7, 98653–98667. [Google Scholar] [CrossRef]
  19. Liu, L.; Wang, S.; Zhao, Z. Radar waveform recognition Based on Time-Frequency Analysis and Artificial Bee Colony-Support Vector Machine. Electronics 2018, 7, 59. [Google Scholar] [CrossRef]
  20. Zhang, M.; Diao, M.; Guo, L. Convolutional Neural Networks for Automatic Cognitive Radio waveform recognition. IEEE Access 2017, 5, 11074–11082. [Google Scholar] [CrossRef]
  21. Wang, Q.; Du, P.; Yang, J. Transferred deep learning based on waveform recognition for cognitive passive radar. Signal Process. 2019, 155, 259–267. [Google Scholar] [CrossRef]
  22. Kishore, T.R.; Rao, K.D. Automatic Intrapulse Modulation Classification of Advanced LPI Radar Waveforms’. IEEE Trans. Aerosp. Electron. Syst. 2017, 53, 901–914. [Google Scholar] [CrossRef]
  23. Zhang, M.; Liu, L.; Diao, M. LPI Radar waveform recognition Based on Time-Frequency Distribution. Sensors 2016, 16, 1682. [Google Scholar] [CrossRef] [PubMed]
  24. Lunden, J.; Koivunen, V. Automatic Radar waveform recognition. IEEE J. Sel. Top. Signal Process. 2007, 1, 124–136. [Google Scholar] [CrossRef]
  25. Lipton, Z.C. A Critical Review of Recurrent Neural Networks for Sequence Learning. arXiv 2015, arXiv:1506.00019. [Google Scholar]
  26. Bendat, J.S.; Piersol, A.G. Engineering Application of Correlation and Spectrum Analysis; John Wiley & Sons: New York, NY, USA, 1980. [Google Scholar]
  27. Sarkar, T.K.; Salazar-Palma, M.; Wicks, M.C. Wavelet Applications in Engineering Electromagnetics; Artech House: Norwood, MA, USA, 2002. [Google Scholar]
  28. Ian, G.; Yoshua, B.; Aaron, C. Deep Learning; MIT Press: Cambridge, MA, USA, 2016. [Google Scholar]
  29. Huang, G.M.; Zuo, W.; Ma, J. Robust radar waveform recognition algorithm based on random projections and sparse classification. IET Radar Sonar Navig. 2014, 8, 290–296. [Google Scholar]
Figure 1. The framework of emitter signal waveform classification.
Figure 1. The framework of emitter signal waveform classification.
Electronics 08 01419 g001
Figure 2. Time-frequency images of six types of signals: (a) Conventional phase (CP); (b) linear frequency modulation (LFM); (c) nonlinear cosine phase modulation (NCPM); (d) binary phase shift keying (BPSK); (e) binary frequency shift keying (BFSK); (f) quadrature phase shift keying (QFSK).
Figure 2. Time-frequency images of six types of signals: (a) Conventional phase (CP); (b) linear frequency modulation (LFM); (c) nonlinear cosine phase modulation (NCPM); (d) binary phase shift keying (BPSK); (e) binary frequency shift keying (BFSK); (f) quadrature phase shift keying (QFSK).
Electronics 08 01419 g002
Figure 3. Autocorrelation images of signals: (a) CP; (b) LFM; (c) NCPM; (d) BPSK; (e) BFSK; (f) QFSK.
Figure 3. Autocorrelation images of signals: (a) CP; (b) LFM; (c) NCPM; (d) BPSK; (e) BFSK; (f) QFSK.
Electronics 08 01419 g003
Figure 4. The signal ε ( k ) curves obtained from Equations (6) and (14): (a) CP; (b) LFM; (c) NCPM; (d) BPSK; (e) BFSK; (f) QFSK.
Figure 4. The signal ε ( k ) curves obtained from Equations (6) and (14): (a) CP; (b) LFM; (c) NCPM; (d) BPSK; (e) BFSK; (f) QFSK.
Electronics 08 01419 g004
Figure 5. The signal ε ( k ) curves obtained from Equations (6), (15), and (16): (a) CP; (b) LFM; (c) NCPM; (d) BPSK; (e) BFSK; (f) QFSK.
Figure 5. The signal ε ( k ) curves obtained from Equations (6), (15), and (16): (a) CP; (b) LFM; (c) NCPM; (d) BPSK; (e) BFSK; (f) QFSK.
Electronics 08 01419 g005aElectronics 08 01419 g005b
Figure 6. Autocorrelation images of signals formed by c l i p ( ) : (a) CP; (b) LFM; (c) NCPM; (d) BPSK; (e) BFSK; (f) QFSK.
Figure 6. Autocorrelation images of signals formed by c l i p ( ) : (a) CP; (b) LFM; (c) NCPM; (d) BPSK; (e) BFSK; (f) QFSK.
Electronics 08 01419 g006
Figure 7. The change of NCPM autocorrelation image at different SNR.
Figure 7. The change of NCPM autocorrelation image at different SNR.
Electronics 08 01419 g007
Figure 8. The change of NCPM two times autocorrelation image at different SNR.
Figure 8. The change of NCPM two times autocorrelation image at different SNR.
Electronics 08 01419 g008
Figure 9. Structure optimization for CNN input layer.
Figure 9. Structure optimization for CNN input layer.
Electronics 08 01419 g009
Figure 10. Image optimization for CNN input layer.
Figure 10. Image optimization for CNN input layer.
Electronics 08 01419 g010
Figure 11. Network structure diagram of CNN1.
Figure 11. Network structure diagram of CNN1.
Electronics 08 01419 g011
Figure 12. Network structure diagram of CNN2.
Figure 12. Network structure diagram of CNN2.
Electronics 08 01419 g012
Figure 13. Network structure diagram of CNN3.
Figure 13. Network structure diagram of CNN3.
Electronics 08 01419 g013
Figure 14. Evaluation method of the signal pre-processing algorithms: (a) Image restoration degree; (b) image stability degree.
Figure 14. Evaluation method of the signal pre-processing algorithms: (a) Image restoration degree; (b) image stability degree.
Electronics 08 01419 g014
Figure 15. Signal recognition results of algorithms under three network models: (a) CNN1; (b) CNN2; (c) CNN3.
Figure 15. Signal recognition results of algorithms under three network models: (a) CNN1; (b) CNN2; (c) CNN3.
Electronics 08 01419 g015
Figure 16. Comparison of the recognition rate between the proposed hybrid model and competitive literature.
Figure 16. Comparison of the recognition rate between the proposed hybrid model and competitive literature.
Electronics 08 01419 g016
Figure 17. Analysis of the influence of changing input layer on signal recognition rate. (a), (b) and (c) are the recognition results obtained by constructing the hybrid model with CNN1, CNN2, and CNN3, respectively.
Figure 17. Analysis of the influence of changing input layer on signal recognition rate. (a), (b) and (c) are the recognition results obtained by constructing the hybrid model with CNN1, CNN2, and CNN3, respectively.
Electronics 08 01419 g017
Table 1. Parameters of signal.
Table 1. Parameters of signal.
Modulation Mode f ( k ) [Hz] ϕ ( k ) [rad] A ( k )
CP f c constant1
LFM f c + B ( k T s ) / τ p w constant1
NCPM f c f c f k sin 2 π f k ( k T s ) 1
BPSK f c 0 or π 1
QFSK f 1 , f 2 , f 3 , f 4 constant1
BFSK f 1 , f 2 constant1
Table 2. Range of signal parameter values.
Table 2. Range of signal parameter values.
Modulation TypeParameterRange
CP f c U ( f s / 6 , f s / 4 )
τ p w U ( 5 e 6 , 20 e 6 )
LFM f c U ( f s / 6 , f s / 4 )
B U ( f s / 20 , f s / 10 )
τ p w U ( 5 e 6 , 7.5 e 6 )
NCPM f c U ( f s / 6 , f s / 4 )
τ p w U ( 5 e 6 , 20 e 6 )
f k U ( f s / 100 , f s / 50 )
BPSK f c U ( f s / 6 , f s / 4 )
L c { 7 , 11 , 13 }
BFSK f 1 U ( f s / 6 , f s / 3 )
f 2 U ( f s / 3 , f s / 2 )
N U ( 256 , 512 )
L U ( 30 , 40 )
N _ r a n f g ( r a n d ( 1 , 50 ) )
QFSK f 1 U ( f s / 6 , f s / 4 )
f 2 U ( f s / 4 , f s / 3 )
f 3 U ( f s / 3 , 5 f s / 12 )
f 4 U ( 5 f s / 12 , f s / 2 )
N U ( 256 , 512 )
L U ( 30 , 40 )
N _ r a n f h ( r a n d ( 1 , 50 ) )
Table 3. Signal pre-processing algorithms.
Table 3. Signal pre-processing algorithms.
No
ProcedureOne DimensionTwo Dimensions
Dimension
1LMSNLMSWTACFICTCWD
2CWDITDITDNone
3NoneNoneIMNone
Table 4. Signal pre-processing speed analysis.
Table 4. Signal pre-processing speed analysis.
Number of Feature ImagesProcessing Speed
10.0356 s
20.0535 s
30.0845 s
Table 5. Speed analysis of classifier recognition.
Table 5. Speed analysis of classifier recognition.
ClassifierProcessing Speed
CNN10.00027 s
CNN20.00081 s
CNN30.00089 s
2CNN1+BiLSTM0.00253 s
2CNN2+BiLSTM0.00341 s
2CNN3+BiLSTM0.00346 s

Share and Cite

MDPI and ACS Style

Ma, Z.; Huang, Z.; Lin, A.; Huang, G. Emitter Signal Waveform Classification Based on Autocorrelation and Time-Frequency Analysis. Electronics 2019, 8, 1419. https://doi.org/10.3390/electronics8121419

AMA Style

Ma Z, Huang Z, Lin A, Huang G. Emitter Signal Waveform Classification Based on Autocorrelation and Time-Frequency Analysis. Electronics. 2019; 8(12):1419. https://doi.org/10.3390/electronics8121419

Chicago/Turabian Style

Ma, Zhiyuan, Zhi Huang, Anni Lin, and Guangming Huang. 2019. "Emitter Signal Waveform Classification Based on Autocorrelation and Time-Frequency Analysis" Electronics 8, no. 12: 1419. https://doi.org/10.3390/electronics8121419

APA Style

Ma, Z., Huang, Z., Lin, A., & Huang, G. (2019). Emitter Signal Waveform Classification Based on Autocorrelation and Time-Frequency Analysis. Electronics, 8(12), 1419. https://doi.org/10.3390/electronics8121419

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop