Parameters of real 2-FSK and 4-FSK patterns

## 1. Introduction

Automatic recognition of modulation is rapidly evolving area of signal analysis. In recent years, much interest by academic and military research institutes has focused around the research and development of recognition algorithms modulation. There are two mains reasons to know the correct modulation type of a signal: to preserve the signal information content and to decide the suitable counter action such as jamming (Nandi & Azzouz, 1998), (Grimaldi et al, 2007), (Park & Dae, 2006).

From this viewpoint, considerable attention is being paid to the research and development of algorithms for the recognition of modulated signals. The need of practice made it necessary to solve the questions of automatic classification of samples of received signals with use of computers and available software.

In this chapter, a new original configuration of subsystems for the automatic modulation recognition of digital signals is described. The signal recognizer being developed consists of five subsystems:

adaptive antenna arrays,

pre-processing of signals,

key features extraction,

modulation recognizer and

output stage.

This chapter describes the use of Walsh–Hadamard transform (WHT) and Karhunen-Loeve transform (KLT) for the modulation recognition in high frequency (HF) and very high frequency (VHF) bands. The input real signal is pre-processed and converted to the “phase image”. The WHT and KLT is applied and the dimensionality reduction is implemented and the classifier recognized the signal. The clustering analysis method was chosen by acclamation for 2-class and 3-class recognition of 2-FSK, 4-FSK and PSK signals. The 2-class and 3-class minimum-distance modulation classifier was created in the MATLAB programme. The tests of designed algorithm were implemented on real signal patterns.

## 2. Orthogonal transforms used for modulation recognition

The utilization of orthogonal transforms for the recognition of various types of modulated signals is described in a number of reference sources. Fourier transform (Ahmed & Rao, 1975), (Jondral, 1991), Haar transform (Ahmed & Rao, 1975), discrete cosine transform (Ahmed & Rao, 1975), (Jondral, 1991), Walsh–Hadamard transform (WHT) (Ahmed & Rao, 1975), (Richterova, 1997, 2001) and Karhunen–Loeve transform (KLT) (Hua & Liu, 1998), (Richterova, 2001), (Richterova & Juracek, 2006) belong to the most frequently exploited and recommended orthogonal transforms. In this chapter, the use of WHT and KLT for the recognition of the frequency shift keying (2–FSK and 4–FSK) signals and the phase shift keying (2–PSK and 4–PSK) signals will be described.

### 2.1. Walsh-Hadamard transform

The Walsh–Hadamard transform (WHT) is perhaps the most well–known of the nonsinusoidal orthogonal transforms. The WHT has gained prominence in various digital signal processing applications, since it can essentially be computed using additions and subtractions only. WHT is used for the Walsh representation of the data sequences. Their basis functions are sampled Walsh functions which can be expressed in terms of the Hadamard matrix. The WHT is defined by relation (Ahmed & Rao, 1975),

where :

B(N) - coefficients of WHT,

N - order of the WHT,

H(N) - N–order Hadamard matrix,

X(N) - signal vector.

An algorithm for the WHT was realized in the MATLAB programme.

### 2.2. Karhunen-Loeve transform

The Karhunen-Loeve transform (Hua & Liu, 1998) (named after Kari Karhunen and Michel Loeve) is a representation of a stochastic process as an infinite linear combination of orthogonal functions, analogous to a Fourier series representation of a function on a bounded interval.

In contrast to a Fourier series, where the coefficients are real numbers and the expansion basis consists of sinusoidal functions (that is, sine and cosine functions), the coefficients in the Karhunen-Loeve transform are random variables and the expansion basis depends on the process. In fact, the orthogonal basis functions used in this representation are determined by the covariance function of the process. The KLT is a key element of many signal processing and communication tasks.

The Karhunen-Loeve Transform (KLT), also known as Hotelling Transform and Eigenvector Transform, is closely related to the Principal Component Analysis (PCA) and widely used in many fields of data analysis.

Let

or in matrix form:

As the covariance matrix

and we can construct an

satisfying

The

or in matrix form:

Here

Now, given a signal vector

where the ith component

Left multiplying

By this transform we see that the signal vector

## 3. Principle of the recognition of FSK and PSK signals

The common fundamental diagram for recognition of 2-FSK, 4-FSK and PSK signals is introduced in Fig. 1 (Richterova, 1999, 2001). General principle of this system for recognition will be described in next text.

The inquiry analog signal

The property of Karhunen-Loeve transform will be used for the recognition of 2-FSK, 4-FSK and PSK signals. All samples of signal pattern are not needed to the proper recognition; it is possible to use the dimensional reduction of the matrix. The proper classification of signal and his enlistment into corresponding group of signals follow up the block of orthogonal transform.

The minimum distance classifier will be used for the solution of the problem of the recognition of 2-FSK, 4-FSK and PSK signals. The principle of minimum distance classifier will be described in the next section.

### 3.1. Phase image

The input signal is given by sequence of the samples corresponding to the digital form of recognition signal. The input vector has the length of 2048 samples. The "phase image" of modulated signal is composed so, that they are generated of points about "the coordinates" - the value of sample and the difference between samples.

These points are mapping into the rectangular net about proportions 32 x 32 so, that a relevant point of net is allocated the number one. If more points fall through into the identical node, then is adding the number one next. These output values are standardized and quantized (Richterova, 1997, 1999, 2001), (Richterova & Juracek, 2006). The “phase images” of 2-FSK and 4-FSK signals are presented on Fig. 2.

Lower frequency of FSK signal corresponds to the ellipse, which lies near to centre of image. Higher frequency of FSK signal corresponds to the ellipse, which is on the margin of image. The “phase image” of PSK is one ellipse.

### 3.2. The 3-class minimum-distance classifier

The minimum-distance classifier is designed to operate on the following decision rule (Ahmed & Rao, 1975), (Richterova, 2001), (Richterova & Juracek, 2006):

A given pattern

Let

Simplification of

Clearly,

The classifier thus computes three numbers

### 3.3. The 2–class minimum–distance classifier

The process of the recognition of 2–FSK and 4–FSK signals by means of the 2–class minimum distance classifier is shown in Fig.4.

Now we briefly describe the block scheme from Fig. 4. A real signal inputs from a receiver via an A/D converter into a computer, where it is processed and stored in the data format WAV.

Then block is pre–processed and used for the filtering and scaling of the input vector of real pattern of 2–FSK or 4–FSK signal. The scaling of real pattern signal is performed in this case and the spectral power density signal is calculated by means of the function PSD, which is implemented in the MATLAB programme. By testing of PSD waveform is proccesed the classification in to FSK or PSK signals.

The phase image (Fig. 2) of the processed pattern of real signal is created by means of algorithms described in (Richterova, 1997, 2001).

The feature vector is a result of the pre–processing of the real signal pattern. The feature vector is classified via a minimum–distance classifier.

Two 2–class minimum–distance classifiers for 2–FSK and 4–FSK signals and two 2–class minimum–distance classifiers for 2–PSK and 4–PSK signals were designed and realized via the learning process. The learning process and the working principle of these minimum–distance classifiers are described in (Ahmed & Rao, 1975), (Richterova, 1997, 2001). The practical application of 2–class minimum–distance classifiers is presented in (Richterova, 1997, 2001).

## 4. Experimental results for 2-class classifier

The performance evaluations of the proposed 2–class minimum–distance classifier (see Fig.4) are introduced for 2–FSK and 4–FSK signals and for 2–PSK and 4–PSK. We have classified 42 real patterns of FSK signals. Parameters of real FSK signal patterns, which were used to experimental process, are presented in Table 1.

Modulation type | Carrier frequency [Hz] | Sampling frequency [Hz] | Modulation velocity [Bd] |

2-FSK | 2400 | 44100 | 150 |

4-FSK | 4800 | 44100 | 150 |

2-FSK | 4800 | 44100 | 100 |

4-FSK | 2400 | 44100 | 100 |

2-FSK | 1200 | 44100 | 50 |

4-FSK | 4800 | 44100 | 200 |

We have classified 32 real patterns of PSK signals. Parameters of real PSK signal patterns, which were used to experimental process, are presented in Table 2.

Modulation type | Carrier frequency [Hz] | Sampling frequency [Hz] | Modulation velocity [Bd] |

2-PSK | 2400 | 44100 | 150 |

4-PSK | 4800 | 44100 | 150 |

2-PSK | 4800 | 44100 | 100 |

4-PSK | 2400 | 44100 | 100 |

2-PSK | 4800 | 44100 | 200 |

4-PSK | 4800 | 44100 | 200 |

The results of the performance are summarized in Table 3 for the classifier based on the WHT and FSK signals.

Real pattern signals | Classified modulation type [%] | |

2–FSK | 4–FSK | |

2–FSK | 87,0 | 13,0 |

4–FSK | 15,0 | 85,0 |

The results of the performance are presented in Table 4 for the classifier based on the WHT and PSK signals.

Real pattern signals | Classified modulation type [%] | |

2–PSK | 4–PSK | |

2–FSK | 77,0 | 23,0 |

4–FSK | 25,0 | 75,0 |

The results of the performance are summarized in Table 5 for the classifier based on the KLT and FSK signals.

Real pattern signals | Classified modulation type [%] | |

2–FSK | 4–FSK | |

2–FSK | 87,0 | 13,0 |

4–FSK | 15,0 | 85,0 |

The results of the performance are presented in Table 6 for the classifier based on the KLT and PSK signals.

Real pattern signals | Classified modulation type [%] | |

2–PSK | 4–PSK | |

2–FSK | 77,0 | 23,0 |

4–FSK | 25,0 | 75,0 |

## 5. Experimental results for a 3-class classifier

The learning process of the 3-class minimum-distance classifier was effected for 30 realizations of the simulation patterns of 2-FSK, 4-FSK, 2-PSK and 4-PSK signals with the modulation velocity 50, 100, 150 and 200 Bd and SNR 15, 20, 30 [dB].

The experimental tests for the 3-class minimum-distance classifier were implemented on 200 realizations of real patterns of 2-FSK, 4-FSK, 2-PSK and 4-PSK signals. Parameters of real signal patterns, which were used to experimental process, are presented in Table 7.

The results of classification of modulation type for 200 realizations of 2-FSK, 4-FSK, 2-PSK and 4-PSK signals are introduced in Table 8.

Modulation type | Carrier frequency [Hz] | Sampling frequency [Hz] | Modulation velocity [Bd] |

2-FSK | 2400 | 44100 | 150 |

4-FSK | 4800 | 44100 | 150 |

2-PSK | 2400 | 44100 | 100 |

4-PSK | 2400 | 44100 | 100 |

Modulation type | Correct classification in [%] for KLT | Correct classification in [%] for WHT |

2-FSK | 75 | 78 |

4-FSK | 72 | 76 |

PSK | 64 | 62 |

The results of classification of PSK signal show, that the classifier based on orthogonal transform have not optimal solution for problems of classification of real PSK signals. We are explored other solution of classification of PSK signals. We want to introduce simple classifier of PSK signal based on cyclostationary feature detection.

This next section introduced the basic cyclostationary descriptors and on simply experiments demonstrates the effectiveness of this approach and its resistance against additive white Gaussian noise (AWGN). The high resistance level results from correlation principle of cyclostationary method because mean value of correlation of AWGN signal is zero.

There are many important aspects in cyclostationary measurement and application; one of them is high sensitivity of obtained results on input parameters that is to say the input sequence length and required cyclic and frequency resolution. In some applications, the input signal pre-processing is necessary. The future research will be focused on automatic modulation recognition.

## 6. Cyclostationary signals descriptors

A random process x(t) is said to be Nth order cyclostationary in the strict sense if its Nth order distribution function exhibits periodicity in time with period T (Gardner et al., 2006)

In practice it is often sufficient to use only second order statistics which leads to the definition of second-order cyclostationarity in wide sense. The key second order statistical characteristic is instantaneous autocorrelation function Rxx(t,τ). So, the process x(t) is said to be cyclostationary in a wide sense if its autocorrelation function is periodic in time with period T

where the instantaneous autocorrelation function is defined as (Semmlow, 2004)

where τ is the time lag and * represents the complex conjugate of the signal x(t). Because the instantaneous autocorrelation function is periodic in time (for all τ) it can be expanded as Fourier series

where

The cyclic autocorrelation function presents the key descriptor of cyclostationary signals in two-dimensional time domain. When we reformulate the equation (20) by substitution (18) we obtained (Gardner, 1991)

Now it is possible to bring out two others, but equivalent, definitions of cyclostationary signals (Gardner, 1991). Firstly, the signal x(t) appears the second order periodicity (x(t) is cyclostationary signal) if and only if the power spectral density of the delay-product signal (18) for some delays τ contains spectral lines at some nonzero frequencies

Cyclic autocorrelation function is the descriptor in time domain. For many applications, it is more useful and convenient to apply the descriptor in frequency domain. The well-known is the Wiener-Khinchin theorem that determines the power spectral density by applying Fourier transform on autocorrelation function. When we apply the Fourier transform on the cyclic autocorrelation function, we obtain the so-called spectral correlation density function

It is convenient to use the form (23) for the real calculation of SCDF

where

is spectral component of x(t) at frequency f with bandwidth 1/T.

In the next sections, there will be shown simple application of cyclostationary descriptors on the tasks of symbol timing recovery and modulation recognition. In both cases the resistibility of this approach against additive white Gaussian noise will be investigated. The experiments were simulated in the MATLAB by partial usage the MATLAB code published in (Costa, 1996).

## 7. Symbol timing recovery experiment

The symbol timing recovery plays the critical role for reliable data detection in the digital communication receiver. In most cases the timing signal is derived from the demodulated signal received x(t) (Barry et all, 2004). One of the methods used is based on spectral-line generation by applying nonlinear transformation. Because the timing signal usually has phase jitter, the phase lock loop (PLL) is used often to reduce it to the level desired. The general form of symbol timing recovery system based on spectral-line generation is shown in Fig. 6.

Some subsystems of symbol timing recovery exploit squaring of signal as nonlinear transformation. But squaring transformation works well only for certain types of data coding. For example, if signal x(t) represents binary data with levels ±1, the square transform generates constant output signal

where

Next problem is to determine the optimal delay τ or by considering CAF optimal cyclic frequency α. It depends on particular modulation format as well as the pulse shape and data correlation. The mathematical derivations for some modulation types are described in (Gardner, 1986), (Lopez-Salcedo & Vazquez, 2003).

An experiment in the SIMULINK was simulated to show the timing recovery of BPSK modulated signal based on multiplication of demodulated signal and its delayed version (the delay equals to half of symbol period

Appropriate timing behaviour of signals are displayed in Fig. 8.

Further, the simulation of timing recovery performance depending on AWGN in transmission channel was done. To the solution this problem, the ideal modulation data generated by the pulse generator (instead of random generator) were be used. The level of AWGN (

## 8. Modulation recognition experiment

The spectral correlation density function provides convenient classification criteria to be applied for modulation recognition. In contrast to traditional power spectral density function, which is often identical for different modulation type (example for BPSK, QPSK), SCDF can lead to different graphs (Gardner et all, 1987). Many works dealing with this problem were published (Gardner et all, 1987), (Qi et all, 2009). So, we try to simulate the SCDF of modulated signals BPSK and QPSK in MATLAB. Both signals were created with following parameters:

sampling frequency 8192 Hz,

carrier frequency 2048 Hz,

symbol rate 512 b/s,

data length 1024.

Obtained SCDF functions are shown in Fig. 10. Although PSK and QPSK signals have the same spectral density function, the graphic difference of both results (spectral correlation density functions of PSK and QPSK signal) is evident. We can split the graph to standard four quadrants and investigate the magnitude of maximum peak in each quadrant separately. While PSK signal has the equivalent level of magnitude for all peaks (quadrants), QPSK signal reaches in two opposed quadrants only about 50% of maximum. Exact value is dependant on statistical character of just utilized data and oscillates for each simulation.

Now we select very simply decision criteria for modulation recognition. If we make a contour of normalized SCDF on level 0,75 we obtain different results for each type modulation type. For PSK signal it has four products while for QPSK signal it has only two products. The situation illustrate Fig. 11.

What can be more surprising is that this simply decision criteria works well under low level of signal to noise ratio (SNR). If we add to the signals additive white Gaussian noise it will work well. The results of simulation for SNR = - 3 dB are shown on Fig. 12 and Fig. 13.

## 9. Conclusion

The operation of the 2–class minimum–distance classifier of 2–FSK and 4–FSK signals has been verified by means of practical programming realized in the MATLAB programme. Results of the classification show that the Walsh–Hadamard transform has better properties for the recognition of FSK compared to Karhunen–Loeve transform (Richterova & Juracek, 2006). The classification efficiency of the 2–class minimum–distance classifier is superior to the linear classifier and quadratic classifier as presented in (Richterova, 2001). Experimental results show that the principal block scheme for the recognition of real pattern of 2–FSK and 4–FSK signals can be used for the special signal analysis in the areas where it is necessary to know or recognize the signal modulation type.

The experimental results indicate, that the 3-class minimum-distance classifier based on Karhunen-Loeve transform and the 3-class minimum-distance classifier based on Walsh-Hadamard transform are able recognize of three types of digital modulated signals. The number of correct classification for real patterns of FSK signals reachs approximately 70 %. The real patterns of PSK signals are classified correctly about 60 %. The designed classifiers are not able recognize modulation type 2-PSK and 4-PSK because the characteristic features of modulation type 2-PSK and 4-PSK are very similar. If we want recognize the groups of 2-PSK and 4-PSK, we must propose other type of modulation classifier.

The practical programming solution and performed experiments were verified a working 3-class minimum-distance classifier for the recognition of 2-FSK, 4-FSK and PSK signals. Results of the classification show that the Walsh–Hadamard transform has better properties for the recognition of FSK compared to Karhunen–Loeve transform. The results show, that this principle can be used for the technical analysis of signals in the branch, where necessary is obtained the information about the modulation type in the automatic system.

Although the theory of cyclostationary signals is known for several decades the practical research in this field still grows. New applications in cognitive radio (Sebesta, 2010), in military field (for example cyclostationary can be used to detect and identify ships, submarines or torpedoes in the ocean (Costa, 1996)) or in other scientific areas are investigated. This paper introduced the basic cyclostationary descriptors and on simply experiments demonstrates the effectiveness of this approach and its resistance against AWGN noise. The high resistance level results from correlation principle of cyclostationary method (useful signal) because mean value of AWGN correlation is zero.

There are many important aspects in cyclostationary measurement and application that were not mentioned here because of bounded length of paper. One of them is high sensitivity of obtained results on input parameters that is to say the input sequence length and required cyclic and frequency resolution. In some applications this leads on necessity of input signal pre-processing. The future research will be focus on automatic modulation recognition.