Next Article in Journal
Machine Learning-Based Peripheral Artery Disease Identification Using Laboratory-Based Gait Data
Next Article in Special Issue
Efficient Informative Path Planning via Normalized Utility in Unknown Environments Exploration
Previous Article in Journal
Internet of Things: Security and Solutions Survey
Previous Article in Special Issue
The Concept of Using the Decision-Robustness Function in Integrated Navigation Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Application of Continuous Wavelet Transform and Artificial Naural Network for Automatic Radar Signal Recognition

Faculty of Mechatronics, Armament and Aerospace, Military University of Technology, 00-908 Warsaw, Poland
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(19), 7434; https://doi.org/10.3390/s22197434
Submission received: 29 August 2022 / Revised: 24 September 2022 / Accepted: 28 September 2022 / Published: 30 September 2022
(This article belongs to the Collection Navigation Systems and Sensors)

Abstract

:
This article aims to propose an algorithm for the automatic recognition of selected radar signals. The algorithm can find application in areas such as Electronic Warfare (EW), where automatic recognition of the type of intra-pulse modulation or the type of emitter operation mode can aid the decision-making process. The simulations carried out included the analysis of the classification possibilities of linear frequency modulated pulsed waveform (LFMPW), stepped frequency modulated pulsed waveform (SFMPW), phase coded pulsed waveform (PCPW), rectangular pulsed waveforms (RPW), frequency modulated continuous wave (FMCW), continuous wave (CW), Stepped Frequency Continuous Wave SFCW) and Phase Coded Continuous Waveform (PCCW). The algorithm proposed in this paper is based on the use of continuous wavelet transform (CWT) coefficients and higher-order statistics (HOS) in the feature determination of selected signals. The Principal Component Analysis (PCA) method was used for dimensionality reduction. An artificial neural network was then used as a classifier. Simulation studies took into account the presence of noise interference with signal-to-noise ratio (SNR) in the range from −5 to 10 dB. Finally, the obtained classification efficiency is presented in the form of a confusion matrix. The simulation results show a high recognition test accuracy, above 99% with a signal-to-noise ratio greater than 0 dB. The article also deals with the selection of the type and parameters of the wavelet. The authors also point to the problems encountered during the research and examples of how to solve them.

1. Introduction

The modern battlefield requires both an effective threat detection system and a system enabling their correct classification. The richness of signals in the radio space means that, both in the case of ELINT and ES, which are elements of EW [1], the development of an effective and fast algorithm for recognizing enemy radars is a key element contributing to the success of the mission. Fast and appropriate recognition of the signal type and the related ability to identify the sources of emissions and the mode of their work is one of the key aspects determining the deployment of forces in military operations. When it comes to military navigation, the efficient ES allows to recognise and avoid dangers in a timely manner. The systems, which provide knowledge of the environment and the deployment of military forces and assets, allow a safe route to be determined, which is particularly important when navigating in areas of military operations. It can also help to identify navigation systems interferences. That is why radio-electronic reconnaissance is an integral part of the fire command and control systems.
In electronic intelligence (ELINT) systems, many emission source parameters, such as radio frequency (RF), direction of arrival (DOA), time of arrival (TOA), pulse width (PW), pulse repetition frequency (PRF), intra-pulse modulation, etc. are determined. In paper [2] an instantaneous frequency profile is used to measure the exotic modulations and their parameters. The paper considers linear frequency modulation (LFM), bi-phase modulation (BPM) and stepped frequency modulation (SFM).
The problem of measuring the parameters of specific radar signals in real time still remains a challenge. The emission identification process is most often carried out using a knowledge-based approach. Paper [3] discusses methods for determining specific radar signal parameters (signatures).
A sorting algorithm for pulse radar based on wavelet transform is proposed in [4]. The authors of the study point to the key role of the accuracy of TOA determination, the estimation of which affects the correctness of the determination of parameters such as PW and PRF. They propose the use of the wavelet transform and Haar wavelet as a tool for sorting radar signals.
A method of recognizing four types of radar signals with three characteristic parameters determined from: CWT eigenvalue, frequency domain moment, kurtosis coefficient, and frequency domain moment skewness coefficient are presented in [5]. The following signals were considered by authors: polynomial phase signal, pseudo code phase modulation and sinusoidal frequency modulation, product composite of pseudo code phase modulation and LFM, convolution composite of pseudo code phase modulation and LFM. For signal to noise ratio (SNR) higher than 0 dB, the probability of signal recognition is greater than 98%.
A proposal for a radar signal recognition method, this time based on the time-frequency transform and high order spectra analysis, is presented in [6]. Time-frequency image (TFI) determined using Wigner–Ville Distribution (WVD) is discussed in [7].
The interference classification and recognition based algorithm is presented in [8] and considers recognition of single-tone interference, multi-tone interference, narrow-band interference and the linear chirp, and impulse and comb spectrum interference.
Another aspect important from the point of view of the real-time recognition process is the classification method. Currently, due to the dynamic development of the field of machine learning, it is more and more often proposed to use artificial neural networks (ANN) as a classifier [9,10]. Using feedforward networks is proposed in [11] and convolutional neural networks (CNN) in [12], usually with features in the form of an TFI, are discussed in [7,13,14,15].
The signal recognition process typically involves pre-processing the radar signal and its analysis to identify signal features, that allow the emitter to be identified. Due to the increasing number of new types of radar signals and their complexity, the electromagnetic environment is becoming more difficult to analyse. Hence, the need to research and develop new algorithms with a wide spectrum of recognition, becomes a pressing problem to solve.
The aim of the article is to present an algorithm for automatic recognition of selected radar signals. The proposed solution based on the use of CWT and HOS is characterised by both high flexibility (related to the variety of possible wavelets and the selection of their parameters) and versatility. A significant advantage of the analysed radar signal classification algorithm is also the possibility of using CWT at individual stages of signal analysis, e.g., to estimate selected parameters of a radar signal. These conditions are met by the feature extraction methods using CWT, which is successfully used in algorithms for sorting radar pulses and determining parameters such as pulse repetition time (PRT), pulse width (PW) or time of arrival (TOA) [5].
A very important aspect discussed in the article is also the influence of parameters used in simulations of analytical wavelets on the properties of the proposed algorithm. Thanks to the analysis of the aforementioned parameters, the algorithm achieves additional adaptation possibilities depending on the specificity of the signals that will be taken into account during the classification. The usefulness of CWT is mainly due to the ability to observe the instantaneous amplitude, phase and frequency of signals simultaneously, without having to plot each of these parameters separately in the time domain. The appropriate selection of the wavelet allows, in turn, to obtain the required resolution in time and frequency, and to emphasise certain features of the analysed signals.
At the same time, the necessity to make online decisions requires the use of classifiers that meet specific time requirements. Therefore, it seems obvious to use artificial intelligence techniques. Although the process of training such a model is time-consuming and requires a large number of training vectors, the trained model allows for signal classification in immediate mode.
The article presents the process of recognizing selected types of radar signals using CWT, HOS and ANN, taking into account analytical wavelets and the influence of their parameters on the obtained results. This new approach to the process of recognizing radar signals allows for the flexibility of the algorithm and the possibility of its comprehensive application by changing parameters such as the signal observation time or the type and parameters of the wavelet used.
So far, the simulations necessary to test the proposed algorithm and confirm its effectiveness in the Matlab environment have been carried out. Ultimately, however, the authors plan to test the algorithm on real signals.

2. Materials and Methods

The effectiveness of algorithms using CWT and ANN is confirmed by research into algorithms for recognising modulations typical of signals used in communications such as amplitude shift keying (ASK), phase shift keying (PSK), frequency shift keying (FSK) or quadrature amplitude modulation (QAM) [16,17,18].
The algorithm proposed in this work identifies specific radar signals (Figure 1): Linear Frequency Modulated Pulsed Waveform (LFMPW), Stepped Frequency Modulated Pulsed Waveform (SFMCW), Phase coded pulsed waveform (PCPW), Rectangular pulsed waveform (RPW), Frequency Modulated Continuous Wave (FMCW), CW (Continuous Wave), SFCW (Stepped Frequency Continuous Wave) and Phase Coded Continuous Waveform (PCCW).
This section will present the methodology for determining radar signals features using CWT, HOS and ANN and the overall structure of proposed radar signal recognition process.

2.1. The Continuous Wavelet Transform (CWT) of Radar Waveforms

The continuous wavelet transform (CWT) enables extraction of transient information associated with amplitude/frequency changes and phase shifts which are characteristics of modulated signals. CWT of a signal x(t) is defined as [19]:
C W T ( τ , a ) = x ( t ) ψ a * ( t ) d t = 1 | a | x ( t ) ψ * t τ a d t
where the ψ ( t ) is called mother wavelet and a is the scaling constant. ψ a * is called baby wavelet and is the translated and scaled version of ψ ( t ) .
In case of digital implementation of CWT, the integral in Equation (1) can be replaced by summation. Assuming that t = kT = k, z = nT = n and scale is an even integer, we could write that [18]:
C W T ( n , a ) = 1 a x ( k ) ψ * k n a
The simulation studies analysed the use of the following analytical wavelets: generalised Morse, Bump and Morlet. The properties of the Bump wavelet are described in [20]. The generalised Morse wavelet, according to [20,21,22], is defined as follows:
ψ β , γ ( ω ) = U ( ω ) a β , γ ω β e ω γ
where a β , γ is a normalization constant, U ( ω ) is the unit step function, and β and γ are parameters controlling the wavelet form. Normalization constant is defined [21]:
a β , γ 2 ( e γ β ) β γ
It follows from the analysis presented in the [20] that by varying β and γ parameters, the generalised Morse wavelets can take a wide variety of forms. For example the γ = 1 family corresponds to the Cauchy (or Paul) wavelet, the γ = 2 correspond to analytic Derivative of Gaussian wavelets, and γ = 3 corresponds to the Airy wavelets family [21].
In Matlab Wavelet Toolbox parameters for a parameterised analytic Morse wavelet are defined as γ , the time-bandwidth product P 2 = β γ and they correspond to the analysis in [20,21,22].
Based on the analysis presented above, simulation studies were carried out to select a suitable wavelet to best reflect the characteristics of the recognised signals. We use the scalogram, which is the absolute value of the CWT of a signal, plotted as a function of time and frequency. Scale is assumed to be proportional to the inverse of the frequency [21].
Figure 2, Figure 3 and Figure 4 show the scalogram for LFMCW, PCCW and SFMCW using the following waveforms:
  • Morse with parameters γ = 27 , β = 27 presented in [20];
  • Bump;
  • Morse with default parameters γ = 3 , β = 20 ;
The best resolution in frequency was observed for the Morse wavelet with parameters γ = 27 , β = 27 . The use of a Morse wavelet with the default parameters results in a scalogram which is fuzzy in frequency. The scalograms for PCCW clearly show the moments at which the phase shift occurred. In the case of the Morse wavelet with the default parameters, the vertical lines corresponding to phase shifts are clearly visible over the entire frequency range.
It would seem, therefore, that the use of the Morse wavelet with default parameters values is appropriate for the identification of phase-coded signals. However, it is clear from Figure 5 that, in the case of stronger noise interference, the use of this wavelet will no longer be effective. For the Bump and Morse with γ = 27 , β = 27 moments of phase change are still visible.
For SFM waveforms the Morse wavelet with γ = 27 , β = 27 is recommended for best frequency resolution (Figure 5).
In order to reflect the characteristics of the individual types of waveforms and to enable their identification after calculating the CWT coefficients, frequency cuts for selected scales can be used, as shown in [17]. Using a median filter for smoothing purposes gives relevant features for classification [23]. The frequency cuts for frequency close to signal carrier frequency are presented in Figure 6.
In the case of signal recognition algorithms, where the CWT is used for feature extraction, an important parameter influencing the correctness of the classification is the scale a. The problem of how to choose the optimal scale for modulation recognition purposes has been addressed in papers [23,24,25]. One of the least complicated solutions to the scale selection problem consisting in summing CWT coefficients for several selected scales (frequencies) was proposed in [24].
During the simulation studies, it was assumed that the exact value of the carrier frequency and bandwidth for a given signal were unknown and needed to be estimated. In order to determine the characteristics of the signals, the CWT frequency cut | C W T ( f 0 , t ) | and the sum of the CWT coefficients for the two selected frequency (scale) ranges: Δ f 1 and Δ f 2 are used. Whereby the f o f c and f c is estimated signal carrier frequency. Examples of CWT frequency cut and sum of CWT coefficients for PCPW, SFCW nad LFMCW are presented respectively in Figure 7, Figure 8 and Figure 9.

2.2. Higher Order Statistics as Signal Features

Higher order statistics (HOS) for feature extraction purposes are considered in [17,26,27]. To determine the characteristics of the specific signals, selected statistical parameters, such as mean value, standard deviation, moments and cumulants of higher orders are proposed. These parameters calculated for CWT frequency cut and summed CWT coefficients for Δ f 1 and Δ f 2 are part of features on the basis of which the classification process is carried out in the proposed recognition algorithm. The cumulants of any order n, assuming that the real signal is analysed, can be calculated using the following expression [27]:
C n = M n m = 1 n 1 ( n 1 ) ! ( m 1 ) ! ( n m ) ! C m M n m
where the moment of n is defined as follows:
M n = 1 N k = 1 N ( x μ ) n
where μ is the mean value of the random variable x. According to the [27] for Gaussian zero-mean random variables, odd numbered cumulants are zero.
In turn, the following expressions for determining higher order mixed moments for complex variables are presented in [26]:
M p q = E [ ( x ( n ) p q x * ( n ) q ]
where * denotes conjugate, q is an integer less than p/2, p is an order of moment.
According to [28] cumulant of n variables x 1 , , x n is defined as:
C u m ( x 1 , , x n ) = π ( | π | 1 ) ! ( 1 ) | π | 1 B π E ( i B x i )
where the | π | goes through the list of all partitions { x 1 , , x n } , B goes through the list of all blocks of the partition | π | , E (   . ) —expectation.
Second, fourth and sixth order cumulants can be represented as [26,29]:
C 20 = C u m ( x , x ) = M 20
C 21 = C u m ( x , x * ) = M 21
C 40 = C u m ( x , x , x , x ) = M 40 3 M 20 2
C 41 = C u m ( x , x , x , x * ) = M 41 3 M 21 M 20
C 42 = C u m ( x , x , x * , x * ) = M 42 | M 20 | 2 2 M 21 2
C 63 = M 63 9 C 42 C 21 6 ( C 21 ) 3
According to [28] cumulants possess one ability which is not present in signal moments: if a variable follows Gaussian statistics, all its cumulants of order higher than two are equal to zero. That is why they are quite often used in automatic modulation recognition context, since they are directly applicable to Gaussian noise. In [28,30] normalised higher-order cumulants are proposed and this form of cumulants was used in our simulation studies.
When calculating the features of modulated signals with the use of cumulants, it is worth using the concept of normalised cumulants mentioned in the article [28,30]:
f e a t u r e 1 = | C 40 | | C 21 | 2
f e a t u r e 2 = | C 41 | | C 21 | 2
f e a t u r e 3 = | C 42 | | C 21 | 2
where C 21 is the signal energy.

2.3. Algorithm Principle

The proposed algorithm for recognising selected radar signals is shown in Figure 10 and consists of the following steps:
1.
signal pre-processing;
2.
signal features extraction based on CWT coefficients, HOS and additional signal features form CWT frequency profile;
3.
feature normalization and reduction of dimensions with PCA;
4.
classification using ANN.
The process of determining the features is as follows:
  • calculation of the CWT coefficients for the frequency range Δ f a l l ;
  • determination of CWT frequency cut for f 0 F ˜ c and its modulus;
  • calculating the sum of the CWT coefficients for each frequency range:
    Δ f 1 C W T ( f , τ ) , Δ f 2 C W T ( f , τ ) , Δ f a l l C W T ( f , τ )
    and the sum of their moduli:
    Δ f 1 | C W T ( f , τ ) | , Δ f 2 | C W T ( f , τ ) | , Δ f a l l | C W T ( f , τ ) | ;
  • calculating the mean values and standard deviations for the:
    Δ f 1 | C W T ( f , τ ) | , Δ f 2 | C W T ( f , τ ) | , Δ f a l l | C W T ( f , τ ) |
    both with and without the median filter applied;
  • determination of joint moments and cumulants for the:
    Δ f 1 C W T ( f , τ ) , Δ f 2 C W T ( f , τ ) , Δ f a l l C W T ( f , τ ) ;
  • additional signal feature form CWT frequency profile defined as: τ C W T ( f , τ ) .
Since CWT coefficients are complex in the calculation of higher order statistics, the Equations used were (6), (8)–(14). Formulas for real signals (4) and (5) can also be used, assuming that they are calculated for CWT coefficient modules. However, this option was not considered in the current study.
The number of determined features is 48 statistical parameters and one additional feature related to the number of profile maxima, for a total of 49. The use of the CWT frequency profile and the calculation of the number of maxima is due to the need for additional features to isolate SFM signals from LFM signals. The suitability of the CWT frequency profile is closely linked to the CWT frequency resolution and therefore to the choice of wavelet type and its parameters. If other types of wavelets are to be used, it is necessary to analyse whether the wavelet provides sufficient frequency resolution first.
The idea of using a profile is used in the analysis of images and is presented, among others, in the work of [31]. This is used for the methods proposed to protect documents against copying. By calculating a profile of a document page in the form of an image, it is possible to determine the individual lines and then to determine the position of the individual words for a given line. According to [31] a profile is the projection of a two-dimensional array into a single dimension. Example of frequency profiles for LFMCW, SFCW and PCCW respectively, are shown on Figure 11, Figure 12 and Figure 13. By analysing the shape of the frequency profile and finding the local maxima of greatest prominence, number of frequencies in multi-carrier signals, like SFM, can be established. The prominence of a local maxima reflects how the peak stands out with respect to its height and location relative to other peaks.
Once all the features have been calculated, normalisation is carried out and then dimensionality reduction is performed using PCA. The feed-forward neural network is proposed as classifier. The default neural network structure has been applied, fully connected, with one hidden layer.
In order to demonstrate the effectiveness of the classification, a set of train and test features was prepared according to the described feature extraction method. Based on these, the learning and testing process of the ANN was carried out. The results of the simulation tests carried out are presented in the next section.

3. Results

During the study, it was assumed that the signals are observed in fixed frequency range and that the carrier frequency for each analysed signal is estimated. Simulation parameters and their values are presented in Table 1.
Simulation studies included verification of the signal feature extraction process by observing the mutual distribution of the first three PCA components of each waveform type. The distribution of the PCA components is shown in Figure 14 and Figure 15 for SNR taking values from 0 to 10 dB and from −5 to 0 dB, respectively.
Comparing the distribution of the PCA components on Figure 15 and Figure 16, it is clear that for S N R 0 dB, the features of the individual signals mix more than for signals with S N R 0 dB.
At first, the features of the LFM and SFM signals overlapped the most. For this reason, during the research, it was decided to use an additional feature. It is related to the CWT frequency profile, and determined on the basis of the number of local maxima.
Figure 16 shows the confusion matrix for the proposed algorithm when considering SNR between 0 and 10 dB. A 100% classification correctness was obtained for all signals except PCCW, where 2% was classified as RCW. Test classification accuracy was equal to 0.9975.
Figure 17 shows the confusion matrix for the proposed algorithm when considering SNR between −5 and 0 dB. There is a noticeable saucer of classification accuracy and so for RCW 93% was correctly classified and 7% was recognised as PCCW, for PCCW 96% was classified correctly and 4% was identified as RCW. Further, for SFMPW 99% was correctly recognised 1% was classified as SFCW. For RPW and PCPW 98% was correctly recognised and 2% was respectively identifies as PCPW and RPW. Test classification accuracy was equal to 0.98.
The last confusion matrix on Figure 18 shows how the classification accuracy of individual signals develops without the additional feature related to the number of maxima of the CWT frequency profile. A slight decrease in classification accuracy for frequency modulated signals, in comparison to confusion matrix presented in Figure 16, is apparent. Test classification accuracy was equal to 0.99.

4. Discussion

The simulation results indicate a high classification accuracy (above 99%) of the selected radar signals in particular for SNR above 0 dB. The decrease in performance for SNR < 0 dB particularly relates to the ability to separate the features of CW signals without modulation from those with phase coding. The problem with SFM and LFM signals separation was solved by the use of additional methods, such as an analysis of the occurrence of local extrema of the CWT frequency profile. A key parameter influencing the ability to determine local extrema is the appropriate choice of wavelet type and its parameters to achieve the highest possible resolution in frequency. In presented work Morse wavelet with γ = 27 and β = 27 is considered.
Furthermore, the method of determining the local extrema of the CWT frequency profile can be used to determine the optimum scales (frequencies) to perform additional CWT frequency cuts and to determine other signal parameters such as the number of frequency levels used in multi-carrier signals. Furthermore, performing CWT frequency cuts at the appropriate scale (frequency) can also improve recognition performance for phase-coded signals. In view of the above, a further development of the algorithm is planned and further simulation tests are to be carried out. At the same time, it is planned to collect real data (signals) in order to confirm the effectiveness of the proposed algorithm in real conditions.
The type of wavelet used also requires further research. The spuriousness of phase-coded signal recognition may require the use of other types of wavelets, e.g., Bump or Paul (analytic Morse with γ = 1 ) which have good time resolution properties.
Analysing the distribution of the PCA components shown in Figure 14 and Figure 15 indicates very good separation using the proposed algorithms between CW and PW signals and within these groups between waveforms with and without frequency modulation.
It is therefore worth considering the application of the proposed algorithm under conditions where SNR > 0 while in the case of large noise interference, its application is limited to a more general classification of signals, e.g., into four groups:
  • CW with FM;
  • CW without FM;
  • PW with FM;
  • PW without FM.
Further analysis is also required to verify the applicability of the algorithm for a wider range of signal parameters such as PRF, PW or code type in PC waveforms.
The undeniable advantage of the proposed algorithm is the ability to use once calculated CWT coefficients for the needs of other algorithms (Figure 19), e.g., to determine the parameters of signals such as PRF, PW, TOA [4], as well as by algorithms using CNN [7,13,14,15], where the CWT coefficients can act as TFI to recognise intra-pulse modulation. The performance of the signal recognition algorithm using CWT, HOS and ANN with the signals used in communication (such as ASK, PSK, FSK and QAM) was confirmed in [17]. Therefore, once computed, the CWT coefficient matrix can be used in parallel by other processing methods, which can significantly affect the functionality of the entire recognition system.
During the simulation tests, an artificial neural network was used with the default network structure, fully connected, with one hidden layer. An open issue that requires further research is the appropriate selection of the ANN structure type and learning model parameters ensuring the highest classification efficiency.

5. Conclusions

In this paper recognition algorithm using CWT, HOS and ANN was considered. The obtained results confirm the effectiveness of using algorithms with CWT and ANN when SNR is above 0 dB. The percentage of correct signal type classification is above 99 % (test accuracy). The decrease in performance for SNR < 0 dB particularly relates to the ability to separate the features of CW signals without modulation from those with phase coding.
There were the PCA used for feature dimension reduction. The time of signal observation included about 20 pulses. The feed-forward network was considered for decision process.
The method of feature extraction based of CWT and HOS seems to be very versatile. The results presented in [17] shows that this is effective way to recognise signals as ASK, PSK, FSK, QAM and the results presented in this article show, that it is also effective for radar signal classification.
The main advantages of the proposed method include:
  • universality (can be used both in the SIGINT and ELINT area);
  • flexibility related to the possibility of selecting the type and parameters of the wavelet depending on the needs;
  • immediate decision (related to the use of ANN);
  • good separation between pulsed and continuous signals;
  • possibility of parallel use of the calculated matrix of CWT coefficients by other algorithms (e.g., signal/emitter parameter estimation).
In turn, the disadvantages of the proposed solution include:
  • visible decline in classification correctness for SNR close to or lower than 0 dB;
  • the need for using additional methods (CWT frequency profile or histogram) to improve the separation abilities within group of signals with frequency modulation;
  • additional research is required on the influence of the applied wavelet on the possibility of classifying individual types of signals (group of signals with FM, group of signals with phase coding, etc.).
Future work will include wider range of selected signals with specific parameters such as pulse width, pulse envelope, pulse repetition time, intra-pulse modulation, phase coding types, etc. Additional analyses related to the selection of the wavelet type and parameters should be carried out. It is also planned to conduct research with the use of real radar signals.

Author Contributions

Conceptualization, M.W.; methodology, M.W.; software, M.W.; validation, M.W.; formal analysis, M.W.; investigation, M.W.; resources, M.W.; data curation, M.W.; writing—original draft preparation, M.W.; writing—review and editing, M.W. and A.K.; visualization, M.W.; supervision, A.K.; project administration, A.K.; funding acquisition, M.W. All authors have read and agreed to the published version of the manuscript.

Funding

The methods and results presented in this paper have been obtained during research works conducted within the university research project entitled “Selected problems of detection, recognition and classification of signals as well as analysis of imaging and radar data”. This work was financed by the Military University of Technology (Warsaw, PL), in 2022, under research project UGB-778.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Adamy, D.L. EW 104: EW against a New Generation of Threats; Artech House: Boston, MA, USA; London, UK, 2015. [Google Scholar]
  2. Niranjan, R.; Rao, C.R.; Singh, A. Real-Time Identification of Exotic Modulated Radar Signals for Electronic Intelligence Systems. In Proceedings of the Emerging Trends in Industry 4.0 (ETI 4.0), Raigarh, India, 19–21 May 2021; pp. 1–4. [Google Scholar] [CrossRef]
  3. Matuszewski, J. The radar signature in recognition system database. In Proceedings of the 19th International Conference on Microwaves, Radar and Wireless Communications, Warsaw, Poland, 21–23 May 2012; pp. 617–622. [Google Scholar] [CrossRef]
  4. Shi, Q.; Zhang, J. Radar Emitter Signal Identification Based on Intra-pulse Features. In Proceedings of the IEEE 6th Information Technology and Mechatronics Engineering Conference (ITOEC), Chongqing, China, 4–6 March 2022; pp. 256–260. [Google Scholar] [CrossRef]
  5. Zheng, Z.; Qi, C.; Duan, X. Sorting algorithm for pulse radar based on wavelet transform. In Proceedings of the IEEE 2nd Information Technology, Networking, Electronic and Automation Control Conference (ITNEC), Chengdu, China, 15–17 December 2017; pp. 1166–1169. [Google Scholar] [CrossRef]
  6. Mingqiu, R.; Yiheng, T. Radar signal cognition based time-frequency transform and high order spectra analysis. In Proceedings of the IEEE International Conference on Signal Processing, Communications and Computing (ICSPCC), Xiamen, China, 22–25 October 2017; pp. 1–4. [Google Scholar] [CrossRef]
  7. Wang, C.; Wang, J.; Zhang, X. Automatic radar waveform recognition based on time-frequency analysis and convolutional neural network. In Proceedings of the IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), New Orleans, LA, USA, 5–9 March 2017; pp. 2437–2441. [Google Scholar] [CrossRef]
  8. Wang, G.; Ren, Q.; Su, Y. The interference classification and recognition based on SF-SVM algorithm. In Proceedings of the IEEE 9th International Conference on Communication Software and Networks (ICCSN), Guangzhou, China, 6–8 May 2017; pp. 835–841. [Google Scholar] [CrossRef]
  9. Granger, E.; Rubin, M.A.; Grossberg, S.; Lavoie, P. A What-and-Where Fusion Neural Network for Recognition and Tracking of Multiple Radar Emitters. Neural Netw. 2001, 14, 325–344. [Google Scholar] [CrossRef]
  10. Hassan, K.; Dayoub, I.; Hamouda, W.; Berbineau, M. Automatic modulation recognition using wavelet transform and neural network. In Proceedings of the 9th International Conference on Intelligent Transport Systems Telecommunications, (ITST), Lille, France, 20–22 October 2009; pp. 234–238. [Google Scholar] [CrossRef]
  11. Nedyalko, P.; Jordanov, I.; Roe, J. Radar Emitter Signals Recognition and Classification with Feedforward Networks. Procedia Comput. Sci. 2013, 22, 1192–1200. [Google Scholar] [CrossRef] [Green Version]
  12. Kiskin, I.; Zilli, D.; Li, Y.; Sinka, M.; Willis, K.; Roberts, S. Bioacoustic detection with wavelet-conditioned convolutional neural networks. Neural Comput. Appl. 2020, 32, 915–927. [Google Scholar] [CrossRef]
  13. Xia, Y.; Ma, Z.; Huang, Z. Over-the-Air Radar Emitter Signal Classification Based on SDR. In Proceedings of the 6th International Conference on Intelligent Computing and Signal Processing (ICSP), Xi’an, China, 9–11 April 2021; pp. 403–408. [Google Scholar] [CrossRef]
  14. Xu, T.; Darwazeh, I. Deep Learning for Over-the-Air Non-Orthogonal Signal Classification. In Proceedings of the IEEE 91st Vehicular Technology Conference (VTC2020-Spring), Virtual Event, 25–28 May 2020; pp. 1–5. [Google Scholar] [CrossRef]
  15. Wang, X.Z. Electronic radar signal recognition based on wavelet transform and convolution neural network. Alex. Eng. J. 2022, 61, 3559–3569. [Google Scholar] [CrossRef]
  16. Hassan, K.; Dayoub, I.; Hamouda, W.; Berbineau, M. Automatic Modulation Recognition Using Wavelet Transform and Neural Networks in Wireless Systems. EURASIP J. Adv. Signal Process. 2010, 2010, 532898. [Google Scholar] [CrossRef]
  17. Walenczykowska, M.; Kawalec, A. Type of modulation identification using Wavelet Transform and Neural Network. Bull. Pol. Acad. Sci. Tech. Sci. 2016, 64, 257–261. [Google Scholar] [CrossRef]
  18. Ho, K.C.; Prokopiw, W.; Chan, Y. Modulation identification of digital signals by the wavelet transform. IEEE Proc. Radar Sonar Navig. 2000, 147, 4. [Google Scholar] [CrossRef]
  19. Addison, P.S. The Illustrated Wavelet Handbook, Introuctory Theory and Applications in Science, Engeneering, Medicine and Finance; Taylor and Francis Group, LLC: Milton Park, UK, 2002. [Google Scholar]
  20. Lilly, J.M.; Olhede, S.C. Generalized Morse Wavelets as a Superfamily of Analytic Wavelets. IEEE Trans. Signal Process. 2012, 60, 6036–6041. [Google Scholar] [CrossRef]
  21. Lilly, J.M.; Olhede, S.C. Higher-Order Properties of Analytic Wavelets. IEEE Trans. Signal Process. 2009, 57, 146–160. [Google Scholar] [CrossRef]
  22. Lilly, J.M.; Olhede, S.C. On the Analytic Wavelet Transform. IEEE Trans. Inf. Theory 2010, 56, 4135–4156. [Google Scholar] [CrossRef] [Green Version]
  23. Gao, Y.; Li, M.; Huang, Z.; Lu, J. A symbol rate estimation algorithm based on Morlet wavelet transform and autocorrelation. In Proceedings of the IEEE Youth Conference on Information, Computing and Telecommunication, Beijing, China, 20–21 September 2009; pp. 239–242. [Google Scholar] [CrossRef]
  24. Chan, Y.T.; Plews, J.W.; Ho, K.C. Symbol Rate Estimation by the Wavelet Transform. In Proceedings of the IEEE International Symposium on Circuits and Systems (ISCAS), Hong Kong, China, 12 June 1997; Volume 1, pp. 177–180. [Google Scholar] [CrossRef]
  25. Lv, Y.; Liu, Y.; Liu, F.; Gao, J.; Liu, K.; Xie, G. Automatic Modulation Recognition of Digital Signals Using CWT Based on Optimal Scales. In Proceedings of the IEEE International Conference on Computer and Information Technology, Xi’an, China, 11–13 September 2014; pp. 430–434. [Google Scholar] [CrossRef]
  26. Zhao, Y.; Xu, Y.; Jiang, H.; Luo, Y.; Wang, Z. Recognition of digital modulation signals based on high-order cumulants. In Proceedings of the International Conference on Wireless Communications and Signal Processing (WCSP), Nanjing, China, 15–17 October 2015; pp. 1–5. [Google Scholar] [CrossRef]
  27. Sobolewski, S.; Adams, W.L.; Sankar, R. Universal Nonhierarchical Automatic Modulation Recognition Techniques For Distinguishing Bandpass Modulated Waveforms Based On Signal Statistics, Cumulant, Cyclostationary, Multifractal and Fourier-Wavelet Transforms Features. In Proceedings of the IEEE Military Communications Conference, Baltimore, MD, USA, 6–8 October 2014; pp. 748–753. [Google Scholar] [CrossRef]
  28. Orlic, V.D.; Dukic, M.L. Automatic modulation classification: Sixth-order cumulant features as a solution for real-world challenges. In Proceedings of the 2012 20th Telecommunications Forum (TELFOR), Belgrade, Serbia, 20–22 November 2012; pp. 392–399. [Google Scholar] [CrossRef]
  29. Mingzhu, L.; Yue, Z.; Lin, S.; Jingwei, D. Research on Recognition Algorithm of Digital Modulation by Higher Order Cumulants. In Proceedings of the Fourth International Conference on Instrumentation and Measurement, Computer, Communication and Control, Harbin, China, 18–20 September 2014; pp. 686–690. [CrossRef]
  30. Han, Y.; Wei, G.; Song, C.; Lai, L. Hierarchical digital modulation recognition based on higher-order cumulants. In Proceedings of the Second International Conference on Instrumentation, Measurement, Computer, Communication and Control, Harbin, China, 8–10 December 2012; pp. 1645–1648. [Google Scholar] [CrossRef]
  31. Brassil, J.T.; Low, S.; Maxemchuk, N.F. Copyright protection for the electronic distribution of text documents. Proc. IEEE 1999, 87, 1181–1196. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Considered type of waveforms.
Figure 1. Considered type of waveforms.
Sensors 22 07434 g001
Figure 2. Scalogram for the LFMCW signal made using wavelets: (a) Morse γ = 27 , β = 27 . (b) Bump (c) Morse (default parameters) γ = 3 , β = 20 .
Figure 2. Scalogram for the LFMCW signal made using wavelets: (a) Morse γ = 27 , β = 27 . (b) Bump (c) Morse (default parameters) γ = 3 , β = 20 .
Sensors 22 07434 g002
Figure 3. Scalogram for the PCCW signal made using wavelets: (a) Morse γ = 27 , β = 27 . (b) Bump (c) Morse (default parameters) γ = 3 , β = 20 .
Figure 3. Scalogram for the PCCW signal made using wavelets: (a) Morse γ = 27 , β = 27 . (b) Bump (c) Morse (default parameters) γ = 3 , β = 20 .
Sensors 22 07434 g003
Figure 4. Scalogram for the SFM pulsed signal made using wavelets: (a) Morse γ = 27 , β = 27 . (b) Bump (c) Morse (default parameters) γ = 3 , β = 20 .
Figure 4. Scalogram for the SFM pulsed signal made using wavelets: (a) Morse γ = 27 , β = 27 . (b) Bump (c) Morse (default parameters) γ = 3 , β = 20 .
Sensors 22 07434 g004
Figure 5. Scalogram for the one pulse phase coded (PC) with Barker code signal made using wavelets: (a) Morse γ = 27 , β = 27 . (b) Bump (c) Morse (default parameters) γ = 3 , β = 20 .
Figure 5. Scalogram for the one pulse phase coded (PC) with Barker code signal made using wavelets: (a) Morse γ = 27 , β = 27 . (b) Bump (c) Morse (default parameters) γ = 3 , β = 20 .
Sensors 22 07434 g005
Figure 6. | C W T ( f 0 , t ) | for f o f c (Morse wavelet used with γ = 27 , β = 27 ): (a) PCCW. (b) LFMCW. (c) SFMCW.
Figure 6. | C W T ( f 0 , t ) | for f o f c (Morse wavelet used with γ = 27 , β = 27 ): (a) PCCW. (b) LFMCW. (c) SFMCW.
Sensors 22 07434 g006
Figure 7. CWT frequency cut for f o f c and sum of CWT coefficients range Δ f 1 and Δ f 2 . PCPW.
Figure 7. CWT frequency cut for f o f c and sum of CWT coefficients range Δ f 1 and Δ f 2 . PCPW.
Sensors 22 07434 g007
Figure 8. CWT frequency cut for f o f c and sum of CWT coefficients range Δ f 1 and Δ f 2 . SFCW.
Figure 8. CWT frequency cut for f o f c and sum of CWT coefficients range Δ f 1 and Δ f 2 . SFCW.
Sensors 22 07434 g008
Figure 9. CWT frequency cut for f o f c and sum of CWT coefficients in range Δ f 1 and Δ f 2 . LFMCW.
Figure 9. CWT frequency cut for f o f c and sum of CWT coefficients in range Δ f 1 and Δ f 2 . LFMCW.
Sensors 22 07434 g009
Figure 10. The overall structure of radar signal recognition process.
Figure 10. The overall structure of radar signal recognition process.
Sensors 22 07434 g010
Figure 11. Frequency CWT profile for LFMCW.
Figure 11. Frequency CWT profile for LFMCW.
Sensors 22 07434 g011
Figure 12. Frequency CWT profile for SFCW.
Figure 12. Frequency CWT profile for SFCW.
Sensors 22 07434 g012
Figure 13. Frequency CWT profile for PCCW.
Figure 13. Frequency CWT profile for PCCW.
Sensors 22 07434 g013
Figure 14. First three components after feature dimension reduction using PCA for SNR in range [0:10] dB.
Figure 14. First three components after feature dimension reduction using PCA for SNR in range [0:10] dB.
Sensors 22 07434 g014
Figure 15. First three components after feature dimension reduction using PCA for SNR in range [−5:0] dB.
Figure 15. First three components after feature dimension reduction using PCA for SNR in range [−5:0] dB.
Sensors 22 07434 g015
Figure 16. Confusion matrix for SNR in range [0:10] dB.
Figure 16. Confusion matrix for SNR in range [0:10] dB.
Sensors 22 07434 g016
Figure 17. Confusion matrix for SNR in range [−5:0] dB.
Figure 17. Confusion matrix for SNR in range [−5:0] dB.
Sensors 22 07434 g017
Figure 18. Confusion matrix for SNR in range [0:10] dB for algorithm version without CWT frequency profile (without local maxima detection).
Figure 18. Confusion matrix for SNR in range [0:10] dB for algorithm version without CWT frequency profile (without local maxima detection).
Sensors 22 07434 g018
Figure 19. Proposal to use CWT coefficients with other methods used in signal recognition.
Figure 19. Proposal to use CWT coefficients with other methods used in signal recognition.
Sensors 22 07434 g019
Table 1. This is a table caption. Tables should be placed in the main text near to the first time they are cited.
Table 1. This is a table caption. Tables should be placed in the main text near to the first time they are cited.
Parameter NameParameter Value
F c 5 MHz ± 0.5 MHz
F s 30 MHz
Bandwidth (for LFM Waveforms)from 2.5 to 5 MHz
Sweep frequency (for SFM Waveforms)from 0.5 to 1 MHz
SNRfrom −5 to 0 [dB] or from 0 to 10 [dB]
Number of pulses in signal≈20
Pulse Width≈0.1 ms for pulsed waveforms
Pulse repetition frequency≈500 Hz
Number of each type of signal500 (400 for training, 100 for testing)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Walenczykowska, M.; Kawalec, A. Application of Continuous Wavelet Transform and Artificial Naural Network for Automatic Radar Signal Recognition. Sensors 2022, 22, 7434. https://doi.org/10.3390/s22197434

AMA Style

Walenczykowska M, Kawalec A. Application of Continuous Wavelet Transform and Artificial Naural Network for Automatic Radar Signal Recognition. Sensors. 2022; 22(19):7434. https://doi.org/10.3390/s22197434

Chicago/Turabian Style

Walenczykowska, Marta, and Adam Kawalec. 2022. "Application of Continuous Wavelet Transform and Artificial Naural Network for Automatic Radar Signal Recognition" Sensors 22, no. 19: 7434. https://doi.org/10.3390/s22197434

APA Style

Walenczykowska, M., & Kawalec, A. (2022). Application of Continuous Wavelet Transform and Artificial Naural Network for Automatic Radar Signal Recognition. Sensors, 22(19), 7434. https://doi.org/10.3390/s22197434

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