Automatic Jamming Modulation Classification Exploiting Convolutional Neural Network for Cognitive Radar
Sensing the external complex electromagnetic environment is an important function for cognitive radar, and the concept of cognition has attracted wide attention in the field of radar since it was proposed. In this paper, a novel method based on an idea of multidimensional feature map and convolutional neural network (CNN) is proposed to realize the automatic modulation classification of jamming entering the cognitive radar system. The multidimensional feature map consists of two envelope maps before and after the pulse compression processing and a time-frequency map of the receiving beam signal. Drawing the one-dimensional envelope in a 2-dimensional plane and quantizing the time-frequency data to a 2-dimensional plane, we treat the combination of the three planes (multidimensional feature map) as one picture. A CNN-based algorithm with linear kernel sensing the three planes simultaneously is selected to accomplish jamming classification. The classification of jamming, such as noise frequency modulation jamming, noise amplitude modulation jamming, slice jamming, and dense repeat jamming, is validated by computer simulation. A performance comparison study on convolutional kernels in different size demonstrates the advantage of selecting the linear kernel.
Cognitive radar proposed in  senses the external electromagnetic environment by recognizing the external emitters and accomplishes antijamming, clutter suppression, and target detection based on the environment information and prior knowledge. Hence, it has significant potential for the optimization and improvement of the performance of the radar. Radar for military applications usually faces a complex electromagnetic environment, especially jamming emitted from electronic countermeasure systems, which must be rejected before target detection. Jamming environment perception is an essential part of the antijamming technology [2–4] for cognitive radar, and it is the premise and foundation on how to employ the various antijamming measures. For example, frequency agility  requiring jamming power measurement on each working frequency point can be regarded as a type of simple environment-cognition waveform transmission technology. The perception of jamming refers to the acquirement of the following parameters and status of the interferences: jamming number, azimuth, modulation type, power, frequency bandwidth, and so on. The direction of arrival and number estimation of the interferences have been discussed widely in the field of array signal processing [6–8]. For instance, multiple signal classification (MUSIC) algorithm  is usually used to measure the direction of the interferences, and the Gerschgorin disk algorithm  is employed to estimate the number of the interferences. Measurement of the jamming power and working frequency bandwidth is widely in use in modern radar systems. However, few papers are devoted to classifying the modulation classes of the radar jamming.
Different antijamming techniques are suitable for different jamming modulation mode. Therefore, in order to make the selection of the antijamming measures more intelligent for cognitive radar, it is essential to recognize the modulation class of different types of interferences. We give an explanation of the antijamming process in a cognitive radar in Figure 1. The antijamming techniques including sidelobe cancellation (SLC), sidelobe blanking (SLB), constant false alarm rate (CFAR) processing, and antijamming transmitting waveform are controlled by the results of jamming sensing. The perception of jamming is implemented after digital beamforming (DBF) or pulse compression of the array data. In this example, the sensing of jamming is composed of four parts: (1) jamming power measurement using DBF; (2) jamming bandwidth measurement using fast Fourier transform (FFT); (3) direction of arrival estimation using MUSIC; (4) automatic modulation classification of jamming (AMCOJ). We focus on the automatic modulation classification of jamming in this paper.
The classes of jamming modulation can be divided into noise frequency modulation jamming , repeater jamming , pulse jamming , frequency shift jamming , and so on. Since there are few studies on AMCOJ in the field of radar, the automatic modulation classification (AMC) of communication and radar signals in the related areas can be used as references.
In the related fields of communication and electronic warfare, there are a great number of studies on automatic modulation classification of communication signals and radar signals. AMC is regarded as an issue of pattern recognition with two processing stages. Firstly, the features of the signals are extracted by various processing algorithms, such as cumulants , radar ambiguity function , and time-frequency analysis . Secondly, the features are used as the input of certain classification algorithms, such as support vector machine  and fuzzy clustering . Recently, some papers propose approaches based on deep learning accomplishing feature extraction and modulation classification simultaneously. The study in  adopts the method of heterogeneous deep model fusion, combining convolutional neural network (CNN), and long short-term memory and realizes AMC of multiple single carrier-modulated signals. An algorithm based on deep learning self-coding network is introduced in , which achieves AMC of five classes of signals: binary phase shift keying (BPSK), 8-phase shift keying (8PSK), 4-quadrature amplitude modulation (4QAM), 16-quadrature amplitude modulation (16QAM), and 64-quadrature amplitude modulation (64QAM). Compared with the traditional sparse autoencoder, the nonnegative constrained autoencoder improves the sparsity and classification accuracy. In , a k-sparse autoencoder-based classifier with characteristics of disordered input data called disordered deep neural network is proposed, reducing the complexity of the actual AMC system and making the classification accuracy close to optimum. The study in  gives a deep learning model based on one-dimensional convolution neural network to classify large-scale coded high-dimensional sequences, thus realizing radar emitter classification. A blind modulation classification method of combining time-frequency distribution and CNN is employed in . By processing time-frequency diagram as an image, a series of frequency-modulated signals, including linear frequency-modulated signal, sinusoidal frequency signal, and polynomial frequency signal, are identified and classified. The study in  accomplishes feature extraction using short-time Fourier transform and AMC based on support vector machine of various types signals, such as single frequency, liner frequency modulation, binary frequency shift keying, 4-frequency shift keying, BPSK, and 16QAM. The superiority of the proposed algorithm lies in that it can classify different modulation types, not just for one specific type. A method based on CNN is carried out to discriminate seven types of pulse repetitive interval modulation in , which is robust to the environment with lost and spurious pulses.
From the above literature, it can be seen that using deep learning to achieve AMC of signals is a hot research topic. As a new research direction of machine learning, deep learning [27, 28] has made breakthroughs in many fields in recent years. Deep learning achieves complex nonlinear function approximation, which has stronger generalization ability, and depicts more abundant essential information of the data samples compared with the conventional shallow model. Due to its excellent performance, deep learning is widely applied in computer vision , speech recognition , image recognition , radar target recognition , classification of radar and communication signals in electronic warfare systems , and so on. Accordingly, on the basis of the complex jamming environment faced by cognitive radar, this paper utilizes the deep learning technology based on CNN to carry out the classification of jamming modulation automatically, realizing jamming environment perception for cognitive radar.
Besides cognitive radar, modulation classification has been considered in cognitive radio systems, for military application in communication systems, and more frequently in optical communication systems [34, 35]. The study of this paper can be extended to modulation classification in these areas .
The rest of this paper is organized as follows. Section 2 introduces the principle of AMCOJ based on CNN for cognitive radar. Section 3 details the multidimensional feature map construction method for seven types of jamming. We design the CNN network structure in Section 4. Section 5 verifies the performance of the proposed classification algorithm through computer simulation. The conclusion is drawn in Section 6.
2. Principle of AMCOJ
The AMCOJ for a cognitive radar consists of multidimensional feature map construction and CNN design. This paper mainly aims at the AMC of seven types of common jamming, namely, noise frequency modulation jamming (NFMJ), noise amplitude modulation jamming (NAMJ), monofrequency pulse jamming (MFPJ), noise frequency-modulated pulse jamming (NFMPJ), dense repeater jamming (DRJ), sparse repeater jamming (SRJ), and slice jamming (SJ, also known as interrupted-sampling and repeater jamming (ISRJ)). The characteristics of the above interferences are as listed in Table 1.
Figure 2 shows the diagram of the jamming sensing network in a cognitive radar. The jamming to a cognitive radar can be roughly divided into two categories, coherent or noncoherent. Since it is coherent to the radar, the waveforms of the repeater jamming before and after pulse compression change significantly. Hence, the data before and after PC processing is selected as the input of the jamming sensing network as shown in Figure 2.
This paper mainly focuses on the radar environment perception processing consisting of multidimensional feature map construction and CNN model design. The differences of different jamming types can be observed in time and frequency domain, so that maps of short-time Fourier transform (STFT) and envelopes before and after PC can be used as features. This paper proposes the method of constructing multidimensional feature map to ensure the correct recognition rate, since single feature map is not effective in discriminating some jamming types. For instance, it is difficult to discriminate noise frequency modulation jamming from noise amplitude modulation jamming in the time-frequency domain. However, there is an obvious difference between them in the time domain before PC. As shown in Figure 2, the feature of time frequency and envelopes can be obtained by STFT analysis and envelope calculation, respectively. The one-dimensional envelope before and after PC can be graphically processed in the time domain to form a two-dimensional image suitable for the use of CNN. The STFT result can be transformed to be a two-dimensional feature map by threshold quantization. The three feature maps are merged to form a multidimensional feature map. The CNN adopting the classical 2-dimensional convolution network structure treats the two-dimensional data of the multidimensional feature map as its input and employs the convolutional kernel to synthetically perceive the input. Through training, the network is continuously updated to its optimum and finally, the seven classes of jamming are recognized.
3. Construction of the Multidimensional Feature Map
Deep learning based on CNN is now widely employed in the field of image recognition. If we transform the jamming features into images in a two-dimensional map, CNN will be possible to be applied to extract their differences as its recognition process of an image. In order to improve the classification accuracy of the jamming using CNN approach, this paper proposes the idea of extracting multiple jamming features from a multidimensional feature map. Envelope features and time-frequency spectrum features are selected to form the feature map. The one-dimensional envelope features before and after pulse compression are plotted in the time domain as planar maps, and the time-frequency features are quantified in amplitude and turn into a two-dimensional planar map, either. Then the above planar maps are combined to form a multidimensional feature map used as the input of CNN. The construction process of each feature map is discussed as follows.
3.1. Time-Frequency Feature Map
Obvious differences can be observed in the time-frequency map of the interferences. The short-time Fourier transform (STFT) can be used to obtain the time-frequency spectrum and can be written aswhere , , , , and represent the discrete-time instant, the discrete frequency, the sample number, the total frequency number, and time delay, respectively, and is the rectangular window function. Formula (1) indicates that is 2-dimensional data and can be transformed into a time-frequency feature map after quantization.
The preprocessing of ( denotes the absolute value) mainly includes image binarization and geometric transformation. Image binarization refers to setting the gray value of the pixels on the image to 0 or 1 by comparing the values of the map with a threshold. Geometric transformation means that the scaling of an image is transformed into a desired size. Taking slice jamming as an example, the time-frequency spectrum gained by STFT is shown in Figure 3(a). The absolute value of STFT of slice jamming is a two-dimensional matrix data. By setting threshold (which is set to the half of the peak value), the value larger than in the matrix is quantized to 1 and vice versa to 0. The binary feature map of slice jamming is shown in Figure 3(b). The analysis of the other types of radar jamming is similar to that of slice jamming, and the binary feature maps are shown in Figures 4(a)–4(f), respectively. Similar to that of noise frequency modulation jamming in Figure 3(a), the binary feature map of noise amplitude modulation jamming is disordered as shown in Figure 3(b). One frequency can be observed in the monofrequency pulse jamming map in Figure 4(c). Since noise frequency modulation is used to modulate the pulse, the frequency variation of noise frequency-modulated pulse jamming is irregular in Figure 4(d). Figures 4(e) and 4(f) show that the distribution of binary feature maps of dense repeater jamming and sparse repeater jamming is still similar to the transmitted linear frequency modulation signal. Therefore, we can conclude that it is difficult to distinguish between noise frequency modulation jamming, noise amplitude modulation jamming, and noise frequency-modulated pulse jamming according to the time-frequency feature alone.
3.2. Envelope Feature Map before PC
According to the characteristics of the seven classes of jamming in the time domain, the envelope features before PC can be extracted. The data before PC (after DBF) is a one-dimensional vector; it can be depicted as a binary map by drawing the envelope in the time domain. The envelope feature map before PC is termed as . We establish a 2-dimensional decision matrix latticed by rows and columns according to the required map size and judge whether there is a plot falling within the corresponding grids. If there is a point falling in the grids, the pixel in the envelope map is recorded as 1. Otherwise, it is recorded as 0. The schematic diagram of the decision matrix is shown in Figure 5. Taking slice jamming as an example, the envelope of slice jamming before PC (after DBF) process is shown in Figure 6(a). By means of the graphical processing, the envelope feature map is illustrated in Figure 6(b). Similarly, the creation of the remaining jamming maps is similar to slice jamming and can be observed in Figure 7. Due to the similarity of monofrequency pulse jamming, noise frequency-modulated pulse jamming, and sparse repeater jamming in the envelope characteristics, we have trouble to distinguish between these three types of jamming barely according to this feature.
3.3. Envelope Feature Map after PC
We further extract the envelope feature after PC. Since the data processed by PC is a one-dimensional vector, it is essential to form a binary map by the above graphical processing. Still taking slice jamming as an example, the waveform and binary feature map after PC are shown in Figures 8(a) and 8(b), respectively. Similarly, the map construction of the remaining six types of jamming is the same and can be observed in Figure 7. However, the similarity of noise frequency modulation jamming and noise amplitude modulation jamming in the envelope feature makes it difficult to discriminate between the two.
3.4. Multidimensional Feature Map
Viewed from the time-frequency feature map, noise frequency modulation jamming, noise amplitude modulation jamming, and noise frequency-modulated pulse jamming are not easy to be distinguished. From the perspective of the envelope before PC, monofrequency pulse jamming, noise frequency-modulated pulse jamming, and sparse repeater jamming are prone to be confused and indistinguishable. Through the envelope feature map after PC, it is difficult to differentiate between noise frequency modulation jamming and noise amplitude modulation jamming. The schematic diagram of multidimensional feature map , as shown in Figure 9, is constructed by combining the above three single feature maps, which mark the upper part as “time frequency,” middle part as “envelope before PC,” and bottom part as “envelope after PC.” The multidimensional feature map of the seven classes of jamming is observed in Figure 7, it can be seen that noise frequency modulation jamming, noise amplitude modulation jamming, and noise frequency-modulated pulse jamming can be distinguished according to the envelope feature map before PC, and monofrequency pulse jamming, noise frequency-modulated pulse jamming, and sparse repeater jamming are apt to discriminate through the time-frequency feature map. Therefore, seven classes of jamming can be classified simultaneously using the multidimensional feature map.
4. Design and Analysis of the CNN Model
The architecture of a typical CNN is generally composed of five stages: input layer, convolution layer, sampling layer, full connection layer, and output layer . We propose a 2-dimensional CNN model according to the complexity of our data. Without using the conventional local field perception in many CNN applications, we design convolution kernels to synthetically perceive the features in the three maps simultaneously. As shown in Figure 9, the network consists of one input layer, one convolution layer, one max-pooling layer, one fully connected layer, and one output layer. The first layer is the input layer, and its inputs are the three feature maps with dimensions of 64 × 512. The second layer is the convolution layer C1 with six convolution kernels of size 192 × 1. Considering the dimension of input feature data, the linear convolution kernel extracts the multiple feature information synthetically from the three feature maps. After the process in this layer, the 2-dimensional feature map is transformed into 6 1-dimensional feature maps with the size of 1 × 512. The third layer S1 is the sampling layer, compressing the feature map into 1 × 64 by max-pooling with the size of 1 × 8 and reducing the complexity of the model. The fourth layer F1 is the full connection layer. Since the feature map generated is 1-dimensional, the feature maps can be connected successively to form a 384 × 1 feature vector. Finally, the output layer maps the feature data to 0 and 1 by a sigmoid function and outputs the corresponding maximum probability value as the jamming class.
The training process of CNN mainly includes two stages: forward propagation and backward feedback. Firstly, the training samples carry out forward propagation network according to the flow of Figure 10. Secondly, comparing the errors between the output and the expectation values, the parameter sensitivity is calculated to update the weights and biases of each layer, that is, to enter the backward feedback network. By constantly updating weights and biases, the network errors decrease as the cost function converges, and the training of the network model is completed. In the testing phase, the weights and biases trained in the training phase are used to classify the test data.
5. Numerical Simulations
5.1. Parameter Setting
In this section, we present a variety of simulation experiments to demonstrate the performance of our proposed AMCOJ algorithm using CNN based on multidimensional feature map (MD-CNN). Seven different modulation candidates are considered here, namely, noise frequency modulation jamming, noise amplitude modulation jamming, monofrequency pulse jamming, noise frequency-modulated pulse jamming, dense repeater jamming, sparse repeater jamming, and slice jamming. Since the robustness of the MD-CNN algorithm is deeply affected by the training data samples, training samples in different conditions should be included in the sample database. Hence, we generate the training samples via changing the jamming principal parameters such as time delay, slice period, and repeat interval. The jamming parameter settings are shown in Table 2. Linear frequency modulation signal is adopted as the radar transmitting signal, the bandwidth of which varies from 4 kHz to 1 MHz and the pulse width is from 40 μs to 120 μs. In addition, the pulse repetition period of the transmitted signal is set as 600 μs and sampling frequency is 2 MHz. The jamming-to-noise ratio (JNR) ranges from 3 to 15 dB subject to the additive white Gaussian noise.
In the network training stage, the feature maps are gained by the processing of the envelope and time-frequency spectrum. We merge three feature maps with the size of 64 × 512 separately into one 192 × 512 multidimensional feature map and use it as the sample input for CNN. The iteration number of CNN is set to 50. The network model is trained by 350 samples and tested by another 175 samples.
5.2. Performance Analysis
The correct recognition rates of the seven jamming classes using MD-CNN algorithm under different JNRs are illustrated in Table 3. Errors are prone to occur under low JNR circumstances due to the high similarity between noise frequency-modulated pulse jamming and noise amplitude modulation jamming in the time domain and time-frequency domain. It can be seen from the data in Table 3 that the numbers of mispredicted samples of each jamming degrade with the increase of JNR.
We compare the performance of the proposed method with the single feature map method using STFT feature map and envelope feature map before and after PC. The total correct recognition rates based on the single feature maps and the multidimensional feature map are plotted in Figure 11. Under the same JNR conditions, the classification accuracy of our proposed algorithm is higher than that of the single feature map algorithm. The result demonstrates that the correct recognition rate of the jamming is effectively improved by constructing the multidimensional feature map. In addition, under the scenario of JNR = 8 dB and above, the correct recognition rate of our proposed method reaches over 95%. The results indicate that the proposed algorithm is effective in discriminating the modulated jamming candidates.
5.3. Comparison Study with Different Convolutional Kernels and Feature Extraction Approach
5.3.1. Comparison of CNN Model with Different Convolutional Kernels
The selection of the CNN kernel is of great significance to the CNN performance. As shown in Figure 12, the CNN algorithms with several types of convolution kernel are simulated as a comparison to our proposed CNN with linear kernel. Under the scenario of JNR = 15 dB, the parameter setting and correct recognition rates of the CNN models with different convolution kernels are shown in Table 4. As can be seen from Table 4, CNN1, CNN2, and CNN3 with convolution kernel of 1-channel 192 × n1 dimension have a better performance than those of CNN4 and CNN5, which means that multifeature comprehensive perception of the three maps simultaneously is superior to local perception of them in succession.
The same conclusion can be drawn that the performance of CNN6, CNN7, and CNN8 with three channels sensing the three feature maps degrades slightly in comparison with CNNs1–3. We take the CNN6 model with three channels as an example to analyze the numbers of mispredicted samples of each jamming class in detail, as shown in Table 5. Compared with the data in Table 3, we can see that the performance of CNN1 with one channel is superior to that of CNN6 with three channels for most of the jamming classes.
5.3.2. Comparison Study of the Feature Extraction Approach Using Radar Ambiguity Function
As a conventional approach, radar ambiguity function has been utilized to extract the modulation type of the jamming feature. Here, as a comparison study, we use the combination of radar ambiguity function and a type of classical CNN called LeNet5 (termed as RA + LeNet5) to classify different jamming signals.
LeNet5 here used is composed of 7 layers: (1) the first convolution layer with 5 × 5 kernels; (2) the first pooling layer with 6 2 × 2 kernels; (3) the second convolution layer with 16 5 × 5 kernels; (4) the second pooling layer with 16 5 × 5 kernels; (5) the third convolution layer with 120 5 × 5 kernels; (6) full connection layer with 120 output; (7) output layer with 7 types of output.
We test the performance of RA + LeNet5 and our proposed method of CNN1 in Table 4 with 350 samples after a training process of 700 samples. The numbers of mispredicted samples of each jamming class are illustrated in Table 6. Fewer mistakes can be observed with CNN1 of this paper at JNR of 6 dB and 8 dB, respectively. As we have pointed out, the superior performance of our proposed CNN1 lies mainly in that it can perceive three features in three maps simultaneously.
The algorithm complexity of RA + LeNet5 and the method of this paper can be analyzed in two aspects, the feature extraction phase and the classification phase.
Algorithm complexity in feature extraction map construction phase is as follows:(1)The method of this paper: there are 64-point FFT of 30 times to calculate the short Fourier transform and 2N multiplications for amplitude calculation without rooting operation, where N refers to the data length(2)RA + LeNet5: to calculate the radar ambiguity function, there is a 64-point FFT of 30 times in the time delay domain and 64-point FFT of 30 times in the Doppler domain
The number of convolutions is used to compare the algorithm complexity in the classification phase, since convolutions can represent the main complexity for these phases:(1)Classification method of this paper: convolution kernel of 192 × 1(2)RA + LeNet5: 6 5 × 5 convolution kernels + 16 5 × 5 convolution kernels + 120 5 × 5 convolution kernels
We can conclude from this comparison that, in the feature extraction map construction phase, the complexity of the two methods is the same, while in the classification phase, the multikernel and multiconvolution layer structure makes RA + LeNet5 more complex than that of our method from the viewpoint of the convolution numbers.
We employ the CNN approach to solve the recognition of jamming for cognitive radar. Considering the differences between different types of jamming in envelope and time-frequency spectrum, the multidimensional feature map is constructed by combining three features of the received signal including amplitude before and after PC and STFT analysis. After the training process, the CNN observes the local features on the combined map and gives the classification result automatically. Computer simulation reveals that the sensing and classification of the jamming classes based on the multidimensional feature map achieve superior performance in comparison with the single feature map method. In addition, multifeature comprehensive perception of the three pictures simultaneously with a linear kernel CNN is superior to CNN with other kernel types.
The data used to support the findings of this study are from the simulation conditions of open literature, and the authors made computer simulations to generate all the simulation data according to these conditions.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
This work was partially supported by the Foundation of Science and Technology on Electronic Information Control Laboratory.
G. Huang and L. Yang, “A radar anti-jamming technology based on blind source separation,” Journal of Circuits & Systems, vol. 3, no. 6, pp. 2021–2024, 2004.View at: Google Scholar
M. A. Richards, J. A. Scheer, and W. A. Holm, Principles of Modern Radar, SciTech Publishing, Raleigh, NC, USA, 2010.
L. Zhao, M. Amin, A. Lindsey et al., “Subspace array processing for the suppression of FM jammers in GPS receivers,” IEEE Transactions on Aerospace & Electronics Systems, vol. 40, no. 1, pp. 80–92, 2004.View at: Google Scholar
H. T. Wu, J. F. Yang, and F. K. Chen, “Source number estimators using transformed Gerschgorin radii,” IEEE Transactions on Signal Processing, vol. 43, no. 6, pp. 1325–1333, 1995.View at: Google Scholar
C. Jiang, M. Gao, Y. Shi et al., “A new method for anti-noise FM interference,” Wireless Sensor Network, vol. 1, no. 4, pp. 294–299, 2009.View at: Google Scholar
S. W. Chen, “The 2-D cosinusoidal phase modulation repeater jamming of SAR,” Journal of Electronics & Information Technology, vol. 31, no. 8, pp. 1862–1866, 2009.View at: Google Scholar
Y. H. Lee and S. J. Kim, “Median-prefiltering-based robust acquisition of direct-sequence spread-spectrum signals in wide-band pulse jamming,” IEEE Transactions on Vehicular Technology, vol. 51, no. 1, pp. 171–179, 2002.View at: Google Scholar
H. X. Huang, Z. T. Huang, and Y. Y. Zhou, “A study on the shift-frequency jamming to SAR,” Journal of Astronautics, vol. 27, no. 3, pp. 463–468, 2006.View at: Google Scholar
G. Wannberg, A. Pellinen-Wannberg, and A. Westman, “An ambiguity-function-based method for analysis of Doppler decompressed radar signals applied to EISCAT measurements of oblique UHF-VHF meteor echoes,” Radio Science, vol. 31, no. 3, pp. 497–518, 2016.View at: Google Scholar
Z. Duona, D. Wenrui, Z. Baochang et al., “Automatic modulation classification based on deep learning for unmanned aerial vehicles,” Sensors, vol. 18, no. 3, p. 924, 2018.View at: Google Scholar
S. Jun, X. Guangluan, R. Wenjuan et al., “Radar emitter classification based on unidimensional convolutional neural network,” IET Radar, Sonar & Navigation, vol. 12, no. 8, pp. 862–867, 2018.View at: Google Scholar
L. Xueqiong, H. Zhitao, W. Fenghua et al., “Towards convolutional neural networks on pulse repetition interval modulation recognition,” IEEE Communications Letters, vol. 22, no. 11, pp. 2286–2289, 2018.View at: Google Scholar