Automatic Classification of Superimposed Modulations for 5G MIMO Two-Way Cognitive Relay Networks

Haithem Ben Chikha and Ahmad Almadhor

Jouf University,College of Computer and Information Sciences,Computer Engineering and Networks Department,Sakaka,72388,Kingdom of Saudi Arabia

Abstract: To promote reliable and secure communications in the cognitive radio network,the automatic modulation classification algorithms have been mainly proposed to estimate a single modulation.In this paper,we address the classification of superimposed modulations dedicated to 5G multipleinput multiple-output (MIMO) two-way cognitive relay network in realistic channels modeled with Nakagami-m distribution.Our purpose consists of classifying pairs of users modulations from superimposed signals.To achieve this goal,we apply the higher-order statistics in conjunction with the Multi-BoostAB classifier.We use several efficiency metrics including the true positive(TP) rate,false positive (FP) rate,precision,recall,F-Measure and receiver operating characteristic (ROC) area in order to evaluate the performance of the proposed algorithm in terms of correct superimposed modulations classification.Computer simulations prove that our proposal allows obtaining a good probability of classification for ten superimposed modulations at a low signal-to-noise ratio,including the worst case (i.e.,m=0.5),where the fading distribution follows a one-sided Gaussian distribution.We also carry out a comparative study between our proposal using MultiBoostAB classifier with the decision tree (J48) classifier.Simulation results show that the performance of MultiBoostAB on the superimposed modulations classifications outperforms the one of J48 classifier.In addition,we study the impact of the symbols number,path loss exponent and relay position on the performance of the proposed automatic classification superimposed modulations in terms of probability of correct classification.

Keywords: Automatic classification;MIMO two-way cognitive relay network;Nakagami-m channels;superimposed modulations;5G

1 Introduction

Recently,a lot of attention has been paid to the two-way relaying (TWR) scheme,which consists of the exchange information between two users via a commonly shared relay in the absence of a direct link between them [1-6].The transmission process under a TWR channel(TWRC) is performed in two-time slots.In the first time slot,the two user nodes send signals to the relay node.In the second time slot,the relay node broadcasts the received signals to the users.In this context,the physical-layer network coding (PNC) introduced in [7],is proposed to allow the relay node to decode a linear function of the received signals,and thereafter to allow each user to decode the incoming message from the other user based on the self-message.The PNC can double the throughput of a TWRC compared to the conventional one-way relay channel by decreasing the time slots for the exchange of one packet from four to two [7,8].It can acheive 1/2 bit of the capacity using a single-input single-output (SISO) Gaussian TWRC and it is assymptotically optimal in the presence of high signal-to-noise ratio (SNR) levels [9].Note that concurrent transmissions based on PNC provide a high spectrum efficiency compared to the network coding and time-division solutions [10].Consequently,the TWR finds applications in a wide range of applications envisaged for 5th generation new radio (5G NR) wireless networks and beyond,including the streaming 4K video,on-line cloud sharing and machine-to-machine communications [11,12].For the purpose of obtaining a good performance when using PNC for these applications,various studies have been carried out with a specific focus on TWRC issues,including the design of symbol mapping [13,14],channel estimation [15,16] and phase synchronization or effect of time [17-20].For satellite communication applications,a satellite can serve as a relay to enable the simultaneous information exchange between two ground stations.In this context,it is necessary to design heterogeneous modulation PNC to allow the exchange unequal amount of data between the two ground stations and/or to exchange data in asymmetric channel conditions of the two ground station-satellite links.Hence,each station may vary the used modulation type and order.Thus,each station needs the knowledge of the modulation scheme employed by the other station in order to properly demodulate the incoming signal.To guarantee an accurate data reception at the stations,a correct detection of the stations modulations from the superimposed symbols is demanded.It is well-known that multiple-input multiple-output (MIMO)system can offer a considerable gain compared to SISO system,especially in environments presenting rich-scattering.In fact,it is a popular technique for increasing the spectral efficiency and the reliability of cellular networks.For that reason,deploying many antenna elements in MIMO systems is one of the most promising technologies in 5G NR systems that can enable beamforming and spatial multiplexing [21].In a MIMO receiver node,the space-time decoder or the spatial demultiplexer and the demodulator are used to recover the transmitted binary information.In fact,the receiver node is the entity that converts the received waves into a binary stream.Consequently,the estimation of the transmitted binary information necessitates a prior knowledge of the communication parameters,such as the number of the source antennas,coding,noise variance,channel matrix,and modulation.To design efficient cognitive radio,several algorithms dedicated to the estimation of the communication parameters have been proposed in the literature.The authors in [22-24],have proposed algorithms for estimating the number of the source antennas.On the other hand,many approaches have been proposed for the detection of the coding [25,26].Other algorithms dedicated to the channel matrix estimation are available in [27,28].The most widely existing modulation classification algorithms in the literature are proposed to estimate a single M-ary modulation [29-34].However,the superposition of two M-ary modulations leads to a significant augmentation of the resulting constellation size in addition to an unusual spatial arrangement [35,36].In fact,the superposition of two modulations with orders M1and M2leads to a modulation with an order upper bounded by M1×M2.Furthermore,the superimposed modulated useful information with the noise results in the dispersion of the constellation points from their appropriate positions.In [35,36],the authors have addressed the problem of modulation classification of superimposed modulations in two-way relaying MIMO systems with PNC under Rayleigh channels.In addition,the zero-forcing (ZF) precoding technique is applied at each source node before transmission.To the best of the authors’knowledge,there is no previous work,which focused on the problem of superimposed modulations classification for MIMO two-way cognitive relay (TWCR) network under realistic channel modeled by Nakagami-m.

In this paper,we propose an algorithm dedicated to the classification of the superimposed users modulations for MIMO TWCR network under Nakagami-m channels.Here,we use a ZF processing at the relay node.The purpose is to classify pairs of users modulations for given superimposed constellations in the case of the presence of Nakagami-m fading.At the relay node,we extract the higher-order statistics (HOSs) of the equalized superimposed symbols as an input to the MultiBoostAB classifier.To evaluate the performance of this latter classifier in terms of modulation classification,several efficiency metrics,such as the true positive rate,falsepositive rate,precision,recall,F-Measure and receiver operating characteristic (ROC) area are used.Simulation results illustrate that the proposed algorithm can achieve a good classification probability of modulations pair even in the worst case of Nakagami-m (i.e.,m=0.5) at a low signal-to-noise ratio (SNR).

The rest of this paper is organized as follows: In Section 2,we model the considered MIMO TWCR network.In Section 3,we describe the proposed modulation classification algorithm.Our main results are illustrated in Section 4.Section 5 gives the conclusion of the paper.

Mathematical Notations:E[.] stands for expectation.(·)-1and(·)Hdenote the inverse and the conjugate transpose operations,respectively.CL×Crepresents the set ofL×Cmatrices over complex field.

2 Considered MIMO TWCR Network

A MIMO TWCR network is considered as shown in Fig.1,where two users denoted as U1and U2and equipped with NU1and NU2antennas,respectively,exchange information through a common relay (R) equipped with NRantennas.

For simplicity,an equal NUnumber of antennas is assumed at all nodes (i.e.,NU=NU1=NU2=NR).Here,we consider that there is no direct link between U1and U2due to the presence of a heavy shadowing.In our transmission model,the message exchange takes place in two time slots.In the first time slot,a multiple access (MAC) phase,in which both users U1and U2,simultaneously send their signals to R,is performed.Here,we consider that U1and U2use two modulations with orders M1and M2.Their superposition leads to a modulation with an order upper bounded by M1×M2.We consider ten combinations of modulation pairs reported in Tab.1.

Table 1:Considered combinations of modulation pairs

For example,in Fig.2,we show the constellation of the superposition of a 4PSK and a 16QAM that contains 64 different points.

AtR,the superimposed modulated useful information is affected by the noise and results a dispersion of the constellation points from their appropriate positions.Hence,at time instantt,we assume that the received signalis given by

wherePU1RandPU2Rare the signal powers atRfrom the usersU1andU2,respectively.HUiR∈CNR×NUirepresents the channel matrix betweenUiandR.It is given by

HUiRare modeled with Nakagami-mifading,i=1,2.Here,the factormirepresents the Nakagami-mseverity parameter of theUi-Rchannel.As the value ofmiincreases,the fading severity decreases.The case ofmi=1 corresponds to the Rayleigh fading.x1,tand x2,tare the modulated data vectors atU1andU2,respectively.Finally,nR,tis a circularly complex Gaussian noise of varianceσ2n.The SNR average of therth relay antenna from thejth antenna of userUiis expressed as

wherei∈{1,2},r∈{1,...,NR}and

In the second time slot,a linear processing is performed on the received signal at the relay nodeR.In fact,a ZF processing is applied to the superimposed signals to eliminate the interpair and the inter-user interferences.The linear processing matrix,denoted by PL∈C2NU×NR,is given by

Here,we suppose that the relay nodeRhas a perfect knowledge of HU1Rand HU2Rsince the estimation of the backward channels can be performed based on pilot signaling [37].Therefore,the transformed signal is written as

In this work,we consider the distances between nodes.Given the fact thatdSDdenotes the Euclidean distance between a source,denoted bySand a destination,denoted byD,the path loss betweenSandDis defined in [38] as

where?represents a constant depending on both the environment and the carrier wavelength,andηdenotes the path-loss exponent and usually varies between two and six.Here,we consider that the two usersU1andU2have the same transmission powerPU=PU1=PU2.To ensure a fair comparison with one-hop transmission,consider thatPUU′=PU1U2=PU2U1to be the power of the received signal at the end-node of theU1U2link (i.e.,direct link).Hence,the power of the received signal atRfrom the userUi,denoted byPUiR,can be written as

wheredU1U2denotes the distance between the user nodesU1andU2,whiledUiRrepresents the distance between the userUi,i=1,2 and the relayrepresents theUi-Rpower gain with respect toPUU′.This latter quotient is denoted byGUiR,i=1,2.The average SNRs of theUi-Rlink is written based on the normalized fading coefficients as

The quotientPUU′/σn2represents the reference SNR and we denote it by.

In the second time slot,Rbroadcasts the superimposed signal to the users with an additional overhead containing the estimated users modulations.A mandatory condition to obtain the correct information at bothU1andU2is the appropriate detection of the modulations used byU1andU2.Thus,Rshould perfectly detect the users modulations from the equalized signal zR.

In this context,we propose an algorithm for classifying superimposed modulations.It is mainly composed of two subsystems.The first subsystem allows to extract the higher-order statistics (HOSs) features from the equalized signal zR,while the second subsystem allows detecting the users modulations pair based on the extracted features and the MultiBoostAB classifier.

In the following,we describe the proposed modulation classification algorithm.

3 Proposed Superimposed Modulations Classification Algorithm

The proposed superimposed modulations classification algorithm is divided into two main steps.The first one consists of extracting a set of appropriate features,while the second one concerns the classification based on supervised machine learning techniques.In the following,we explain these two steps.

3.1 Extraction of Discriminating Features

The higher-order statistics (HOSs) composed by the higher-order moments (HOMs) and the higher-order cumulants (HOCs) have shown in several recent existing works in literature their ability to classify modulations for MIMO systems [39].In fact,each modulation scheme can be characterized by a set of HOMs and HOCs.The use of HOSs up to order eight allow the correct classification of various modulation types [40].

Thejth-order HOM of the equalized sequence at theathantennais given by [41]

An estimation of the HOMs can be expressed as

Thejth-order HOC of the z(Ra)signal can be expressed as

Thejth-order HOC may be written as a function of lower and equal ordered HOMs as follows

where Ψ runs through the list of all partitions of {1,...,j},φruns through the list of all blocks of the partition Ψ andδis the elements number of the partition Ψ.We raise each HOC to the power 2/jsince the magnitude of HOCs increases with their order [42].

The process of classification of modulations pair for the received signal yR,tis illustrated in Fig.3.A training phase is firstly launched to build a model using the MultiBoosting through the use of a learning database (LDB).Then,the test phase is done to classify the modulations pair of yR,tbased on the model that is already built with the LDB.In the following,we describe the MultiBoosting classifier.

3.2 MultiBoostAB

In this work,we use the MultiBoosting (MultiBoostAB) classifier,which is a combination of the Boosting and the Wagging techniques [43].We present in the Algorithm 1 the pseudocode of MultiBoostAB classifier.

The idea is to harness the benefits provided by both techniques.In fact,this classifier takes advantage of Wagging’s superior variance reduction in addition to the AdaBoost’s high bias and variance reduction.Here,we employ the C4.5 (J48) [44] as a base learning algorithm since with this latter MultiBoost classifier provides a good prediction comparing to the AdaBoost classifier.To prove the effectiveness of MultiBoostAB operating with J48 classifier in superimposed modulation classification,we carry out a comparative study with the J48 classifier alone that outperforms the performance of multilayer perceptron classifier trained with resilient backpropagation training algorithm [45].

Algorithm 1: MultiBoostAB algorithmimages/BZ_1765_265_498_1091_2310.png

3.3 Metrics Used for Performance Evaluation of Classifiers

In this study,we compare between classifiers using true positive (TP) rate,false positive(FP) rate,precision,recall and F-Measure metrics.The precision,recall and F-measure are given respectively as

4 Simulation Results

Simulation experiments are conducted to demonstrate the advantages of the proposed automatic classification modulation algorithm.Here,we apply our proposal to classify ten combinations of modulation pairs,i.e.,M={(16QAM,16QAM),(16QAM,64QAM),(16QAM,2PSK),(16QAM,4PSK),(64QAM,64QAM),(64QAM,2PSK),(64QAM,4PSK),(2PSK,2PSK),(2PSK,4PSK),(4PSK,4PSK)}.

For each pair inM,we construct a training set with 200 superimposed signals,where the user messages x1,tand x2,tofU1andU2,respectively,are created in a random manner.For each signal,the number of symbolsNis fixed to 10000 symbols.Here,the all channels are subject to Nakagamimfading.1000 Monte Carlo trials are generated as a test set for each superimposed pair inM.Thus,the test set contains 10000 Monte Carlo trials results (i.e.,Numtotal=10×1000).We consider a free-space path loss model withη=2 and the number of antennas at each node is fixed toNU=4.In the training and test phases,the set of features is building using the HOSs(HOCs and HOMs) of z(a)R.Finally,we use MultiBoostAB with J48 as a base classifier,where the number of training subsetNTSis equal to 10 [46].This latter setting shows that the MultiBoostAB classifier provides a good compromise between the probability of modulation classification and the speed of the training phase.

In this work,the probability of the correct classification is computed by

whereNψmis the number of trials for which the pair modulationsψm∈Mis perfectly classified.For each test trial,a collaboration between allNRantennas is made in order to take the decision.Indeed,the pair of modulations having the majority of votes represents the estimated modulations pair.

4.1 Accuracy of the MultiBoostAB Classifier

We firstly evaluate the performance of the MultiBoostAB classifier using a 10-fold crossvalidation [47] on the training set described above.Tab.2 displays the detailed accuracy by superimposed modulations.By analyzing the average of the TP rate,FP rate,precision,recall,F-Measure and receiver operating characteristic (ROC) area,it is clearly shown that the Multi-BoostAB offers a good classification performance.In fact,the values of TP rate,precision,recall,F-Measure and ROC area are very close to 1 and the value of the FP rate is very close to 0.Therefore,the MultiBoostAB classifier is efficient to automatically classify superimposed modulations.

4.2 Impact of the Nakagami-m Fading Parameter

Fig.4 shows the impact of the channel fading severity on theP,forManddU1R=dU2R=0.5.It is apparent that the probability of the classification decreases when decreasing fading parametermsince the fading becomes more and more severe.We also observe that the proposed algorithm has the ability to classify superimposed modulations at low SNR even in the case where the fading severity parameter is set to 0.5 (i.e.,worst-case).By comparing the performance of the MultiBoostAB with J48 classifier alone,it is clearly shown that the MultiBoostAB classifier offers a gain compared to J48 classifier in terms of superimposed modulations classification.For example,at 95% of superimposed modulations classification,the MultiBoostAB classifier provides a SNR gain of about 0.5 dB compared to J48 for the case wherem=0.5.Consequently,the MultiBoostAB classifier is more appropriate in superimposed modulations classification.

Table 2:Detailed accuracy by superimposed modulations for the MultiBoostAB classifeir with NU=4,N=10000,η=2,dU1R=dU2R=0.5,m=0.5 and =5 dB

Table 2:Detailed accuracy by superimposed modulations for the MultiBoostAB classifeir with NU=4,N=10000,η=2,dU1R=dU2R=0.5,m=0.5 and =5 dB

Mi TP rate FP rate Precision Recall F-measure ROC area(16QAM,16QAM) 0.976 0.001 0.994 0.976 0.985 0.999(16QAM,64QAM) 0.979 0.003 0.97 0.979 0.974 0.999(16QAM,2PSK) 0.999 0.001 0.994 0.999 0.996 1(16QAM,4PSK) 0.995 0 0.996 0.995 0.996 0.999(64QAM,64QAM) 0.991 0.002 0.984 0.991 0.988 1(64QAM,2PSK) 0.994 0 0.997 0.994 0.996 1(64QAM,4PSK) 0.996 0 0.996 0.996 0.996 0.998(2PSK 2PSK) 1 0 0.999 1 0.999 1(2PSK,4PSK) 0.998 0 0.999 0.998 0.998 0.999(4PSK,4PSK) 1 0 0.999 1 0.999 1 Average 0.993 0.001 0.993 0.993 0.993 0.999

4.3 Impact of the Symbols Number

Fig.5 presents the average probability of correct classification,P,as a function of,for many values of symbols numberNusing MultiBoostAB classifier.The increase ofNleads to an improvement inPsince the HOMs estimation accuracy calculated in (8) increases when increasingN.

4.4 Impact of the Relay Position

As seen in Eq.(6),we can incorporate different relay positions,where all distances involved in the calculation of the power gain are relative to the distance between the two usersU1andU2.In fact,we assume thatdU1R+dU2R=1.

Fig.6 illustrates the average probability of correct classification of the proposed algorithm,P,as a function ofat different relay positions forM.It is clearly shown that,when the relay is located exactly in the middle (i.e.,dU1R=dU1R=0.5),the best performance is obtained.For example,a good performance (P～=100%) is reached atdU1R=0.5,while only 97% is achieved atdU1R=0.1 for=5 dB.It can be also shown from this Figure that the proposed algorithm achieves an excellent performance (P～=100%) at all relay positions and for a low SNR (=8 dB).

4.5 Impact of the Path Loss Exponent

Fig.7 presents the average probability of correct classification of the proposed algorithm,P,as a function of,for various values ofη,usingManddU1R=dU2R=0.5.The increase ofηleads to a considerable improvements inP.Here,the gain is obtained thanks to the increase of theatRwhen increasingηfor the same referenceas seen in Eq.(6).

4.6 Impact of the Antenna Number

Fig.8 shows the average probability of correct classification of the proposed algorithm,P,as a function of,forNU=NU1=NU2=NR={4,8},usingManddU1R=dU2R=0.5.One can see thatPis significantly increased when settingNUto 8.In fact,a good performance is achieved at=1 dB forNU=8,whereas the same performance is obtained at=5 dB forNU=4.

5 Conclusion and Future Work

We have proposed an automatic classification algorithm of superimposed modulations designed for MIMO TWCR network over Nakagami-mchannels.At the relay node,we have extracted a set of HOSs from the superimposed received symbols as features extraction.Then,we have employed the MultiBoostAB classifier.Simulations were performed to show the performance of the proposed classification modulation algorithm based on several metrics.We have demonstrated that our proposal has the ability to provide good performances at a low SNR in the case where the fading severity parameter is set to 0.5 (i.e.,worst-case).We have also carried out a comparative study between our proposal using MultiBoostAB classifier with J48 classifier.Through simulation results,we have clearly showed that the performance of MultiBoostAB on the superimposed modulations classifications outperforms the one of J48 classifier.Finally,we have studied the impact of the symbols number,path loss exponent and relay position on the performance of the proposed automatic classification superimposed modulations in terms of probability of correct classification.

In the future work,we will investigate the use of deep learning based neural networks in order to further improve the probability of the correct classification of superimposed modulations at low SNR values.

Funding Statement: This work was supported by Jouf University.

Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the present.