Next Article in Journal
Anti-Disturbance Control for Quadrotor UAV Manipulator Attitude System Based on Fuzzy Adaptive Saturation Super-Twisting Sliding Mode Observer
Previous Article in Journal
Soil Is Still an Unknown Biological System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Performance of the Multiple Active Antenna Spatial Modulation with 3-Dimensional Constellation

1
Department of EEC Engineering, KoreaTech, Cheonan 31253, Korea
2
Department of CS, Northeastern Illinois University, 5500 N St Louis Avenue, Chicago, IL 60625, USA
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2020, 10(11), 3718; https://doi.org/10.3390/app10113718
Submission received: 1 May 2020 / Revised: 16 May 2020 / Accepted: 25 May 2020 / Published: 27 May 2020
(This article belongs to the Section Electrical, Electronics and Communications Engineering)

Abstract

:
In spatial modulation (SM), a single signal symbol is transmitted from a given physical antenna, where both the signal symbol and the antenna index carry information. SM with multiple active antennas (MA-SM) transmits several signal symbols from a combination of antennas at each channel use, thereby increasing the spectral efficiency. MA-SM is proposed in combination with a new 3-dimensional constellation, where signal symbols transmitted from a given antenna combination are rotated before transmission. In this paper, we derived an upper-bound on the error probability of the MA-SM as a function of the rotation angles. The search for the optimal rotation angles is modeled as a multi-objective optimization problem. We concluded based on both analytical and simulation results that the 3-dimensional constellation with the optimal angles achieved negligible improvement. Therefore, we do not recommend using the 3-dimensional constellation with the MA-SM system.

1. Introduction

Index modulation (IM) is a relatively new technique in which information is conveyed to the receiver using the conventional signal symbols and/or the index(es) of additional system resources [1]. The index represents a physical antenna in spatial modulation (SM) and space shift keying (SSK) [2,3], spreading codes [4], polarities [5], sub-carriers [6], combinations of angles [7], virtual channels [8], among others. SM attracted an increasing attention because it achieves a multiple-input multiple-output (MIMO) gain in terms of spectral efficiency while keeping the transmitter as simple as that of the single-input systems.
The generalized SM (GSM) reduces the number of transmit antennas required to achieve a given spectral efficiency as compared to SM. GSM transmits a single signal symbol from a combination of antennas, instead of a single antenna as in SM [9]. The parallel SM (PSM) divides the available antenna set into groups, and the conventional SM is independently performed in each group using the same signal symbol. Therefore, the PSM reduces the number of transmit antennas at no cost in terms of the spectral efficiency [10,11]. Recently, a new constellation that reduces the asymptotic error performance was designed in [12].
Quadrature SM (QSM) is an extension of SM, where the real and imaginary part of a single signal symbol are transmitted from an in-phase and quadrature spatial dimension, respectively, leading to an improvement in the spectral efficiency [13]. To further increase the spectral efficiency, an improved QSM (IQSM) is presented to transmits two signal symbols at each channel use, where the real parts and imaginary parts are transmitted from antenna combinations of length of 2 on the orthogonal cosine and sine carriers, respectively [14]. A high-performance parallel implementation of the improved QSM was recently proposed in [15].
A complex QSM (CQSM) was proposed in [16,17], where two signal symbols are transmitted from two physical antennas whose indices also carry information. The two signal symbols are drawn from two distinct constellations, where the second is a rotated version of the first. The optimization and the convergence tendencies of the constellation are investigated in [17]. It was shown that the convergence of rotation angle is dominated by the probability of the event that the two symbols are transmitted from the same antenna. An improved CQSM is proposed in [18], where the transmitter is equipped with an additional antenna that is used only when the two signal symbols are supposed to be transmitted in the CQSM. The first signal symbol is transmitted from its designated antennas; the second signal symbol from the additional antenna. The design of the modulation set of the improved CQSM is addressed in [19]. In this category of algorithms, the optimization of the rotation angles plays a very crucial role in determining the error performance of the system. As such, the constellation should be carefully optimized, as mentioned in the above works.
A generalized multiple active antennas SM (MA-SM) was proposed in [20], where each activated antenna transmits an independent signal symbol, leading to further improvement in the spectral efficiency. In addition to the idea of MA-SM, authors in [20] proposed two additional algorithms for the MA-SM: 1) A simple detection algorithm, and 2) 3-dimensional constellation. The proposed simple alternative to the optimal maximum-likelihood (ML) detector comes at a cost in terms of the bit-error rate (BER) and the diversity order. Building upon the pioneering work in [20], Several works have addressed the detection process of the MA-SM and proposed more efficient detectors with low computational complexity (c.f. [21,22]). On the other hand and to the best of our knowledge, the merits of the 3-dimensional constellation was not verified in the original work nor by any other researchers.
Contributions: We believe that the verification of the proposed technique in the literature should be as important as introducing new techniques. This work is a complete effort to verify the merits of the 3-dimensional constellation proposed in [20] for the MA-SM system. To this end, we derive the analytical formulae of the upper bound on the pair-wise error probability (PEP) of the MA-SM for the system parameters provided in Table I of [20]. The derived formulae are functions of the rotation angles introduced through the 3-dimensional constellation. We model the search of the rotation angles in [20] as a multi-objective optimization that is solved to obtain the optimal rotation angles that reduce the error probability. Based on the provided analytical analysis and simulation results, we conclude that the 3-dimensional constellation obtained from the optimal angles derived in this paper slightly improve the BER performances of the MA-SM with signal-to-noise ratio (SNR) gains of less than 0.5 dB, as compared to the performance of MA-SM when the rotation angles that were given in [20] are used. Overall, we recommend that the 3-dimensional constellation should not be used with the MA-SM due to the marginal improvement in the BER performance.
In his paper, we used exactly the same transmitter structure described in [20]. Specifically, we used the same number of transmit antennas, the same values of the angles, and the same number of active antennas at each channel use. The work in [20] is pioneering and of great importance to the development of spatial modulation techniques. Our work should be regarded as a continuity rather than a replacement.
The rest of this paper is organized as follows. In Section 2 and Section 3, we introduce the related work and system model of MA-SM, respectively. In Section 4, we derived the upper bound on the PEP of the MA-SM in terms of the rotation angles for two system configurations and support the analytical performance with numerical results. Simulation results were presented in Section 5 and we draw the paper’s conclusions in Section 6.

2. Related Work

2.1. Spatial Modulation

We consider a single base station with N T transmit antennas that communicates with a single user equipped with n r receive antennas. In the conventional SM [2], a single signal symbol, s k is transmitted from a single antenna, where both the signal symbol’s index and the antenna index carry information. As such, the system equation is given by:
y = h i s k + n
where h i is the i-th column of the channel matrix H of size n r × N T , n is the additive white Gaussian noise vector of length n r and y is the received vector of length n r . The elements of H are independent and follow a centered circularly symmetric Gaussian distribution with variance of 1. The signal symbol s k is drawn from a conventional phase shift keying (PSK) or quadrature amplitude modulation (QAM) constellation Ω , whose size is L, and q = l o g 2 ( L ) is the number of bits conveyed by each signal symbol.

2.2. Generalized Spatial Modulation

SM is generalized in [9] to address the limitations of the log-2 number of transmit antennas. In GSM, two or more transmit antennas are activated to transmit a single signal symbol. The transmit antennas are selected such that l o g 2 C N a N T is the number of possible antenna combinations, where N a is the number of active antennas and · is the floor operator. The spectral efficiency of GSM is accordingly given by:
S E g s m = q + l o g 2 C N a N T
and the system equation is given by:
y = s k N a i = 1 N a h i + n

3. System Model of MA-SM

In order to increase the spectral efficiency, several independent signal symbols are transmitted from a combination of activated antennas in the MA-SM system at a cost of additional radio frequency (RF) chains. Therefore, the MA-SM combines the advantage of low complexity of SM system and the high spectral efficiency of spatial multiplexing (SMux). Let the number of active antennas at each channel use be denoted by N P . The spectral efficiency of MA-SM achieved by the spatial symbols, i.e., the indexes of the available combinations, is given by:
S E s p a = l o g 2 C N P N T
where C N P N T is the set of combinations of length N P out of the available N T transmit antennas. The set of all antenna combinations that can be used for transmission is denoted by Γ = { γ 1 , , γ N } . Following the notation in [20], we define ξ j as a vector of length N T with binary entries of 0 or 1. The values associated with these vectors are explained in light of the following example. Let N T = 4 and N P = 2 , then:
γ 1 = ( 1 , 2 ) ξ 1 = ( 1 , 1 , 0 , 0 ) γ 2 = ( 1 , 3 ) ξ 2 = ( 1 , 0 , 1 , 0 ) γ 3 = ( 1 , 4 ) ξ 3 = ( 1 , 0 , 0 , 1 ) γ 4 = ( 2 , 3 ) ξ 4 = ( 0 , 1 , 1 , 0 ) γ 5 = ( 2 , 4 ) ξ 5 = ( 0 , 1 , 0 , 1 ) γ 6 = ( 3 , 4 ) ξ 6 = ( 0 , 0 , 1 , 1 )
Because N should be a power of two, only N = 4 combinations are used for transmission. The system equation is then given by:
y = 1 N P H s + n ,
where the vector s contains exactly N P non-zero elements whose locations are determined by γ .
The selection of the antenna combinations used for transmission out of the available ones is based on maximizing the Hamming distance among the ξ vectors. In the sequel, the same set of combinations selected in [20] will be used to conduct a fair comparison and analysis. As a contribution in [20], a rotation angle is associated with each of the used combinations, where all symbols transmitted from a given combination of antennas are rotated with the same angle. Table 1 summarizes these parameters for three configurations considered in [20].
Based on the above system description, the total spectral efficiency of the MA-SM system is given by:
S E m a s m = S E s i g + S E s p a = N P × q + l o g 2 C N P N T

4. Performance Analysis of the MA-SM

4.1. Performance Analysis of MA-SM with Maximum-Likelihood Receiver

In the following, we derive the closed-form formula for the error probability of the MA-SM system. For the sake of simplicity, the upper bound on the error probability is derived assuming N P = 2 . The generalization is straightforward. Let
g m 1 = ( h i 1 s k 1 + h i 2 s k 2 ) e j θ m 1 and g m 2 = ( h i ^ 1 s k ^ 1 + h i ^ 2 s k ^ 2 ) e j θ m 2
be defined as two received codewords, where θ m 1 and θ m 2 are the rotation angles associated with the antenna combination ( i 1 , i 2 ) and ( i ^ 1 , i ^ 2 ) , respectively. The PEP is defined as:
P r g m 1 g m 2 = Q g m 1 g m 2 2 2 σ n 2 ,
where Q ( · ) is the Gaussian tail function, or simply the Q-function. The unconditional pairwise error probability (UPEP) is obtained by taking the expectation of (8) over the channel matrix H . A closed-form formula of (8) is given as follows [2]:
P r g m 1 g m 2 = μ m 1 , m 2 n R l = 0 n R 1 n R 1 + l l 1 μ m 1 , m 2 l ,
where:
μ m 1 , m 2 = 1 2 1 ρ · d m 1 , m 2 ( t x ) 2 4 + ρ · d m 1 , m 2 ( t x ) 2 ,
and the square Euclidean distance between the two vector symbols is given by:
d m 1 , m 2 ( t x ) 2 = s m 1 s m 2 2 .
In (10), ρ is the signal-to-noise ratio (SNR). At high SNR, the asymptotic pairwise error probability is approximated as follows:
P r g m 1 g m 2 2 n R Γ ( n R + 0.5 ) π ( n R ) ! 1 ρ d m 1 , m 2 ( t x ) 2 n R = 2 n R 1 n R ρ n R d m 1 , m 2 ( t x ) 2 n R
Finally, the union bound on the UPEP is obtained by averaging over all the hypotheses of the transmit vectors as follows:
P r e 1 2 M m 1 = 1 2 M m 2 = 1 2 M P r g m 1 g m 2 = 2 n R 1 n R ρ n R 2 M m 1 = 1 2 M m 2 = 1 2 M d m 1 , m 2 ( t x ) 2 n R = 2 n R 1 n R ρ n R 2 M i 1 = 1 N p e r m s i ^ 1 = 1 N p e r m s k 1 L k ^ 1 L k 2 L k ^ 2 L d m 1 , m 2 ( t x ) 2 n R = 2 n R 1 n R ρ n R L 2 k = 1 B f k Ω k
where M is the spectral efficiency of the MA-SM system, N p e r m s is the number of antenna combinations that can be used for transmission, and f k is the frequency associated with the term Ω k . The derivation of the Ω terms depends on the distance,
δ ( ξ i , ξ j ) = X O R ( ξ i , ξ j )
where X O R ( · ) is the logical exclusive or operation.

4.2. Analysis of the MA-SM Performance for Configuration I and II

In this subsection, the performance of configuration I and II from Table 1 are analyzed because they have the same N P and N T . The distance, δ ( ξ i , ξ j ) has three distinct values: 0, 2 and 4. These values are analyzed in the following.
Case 1: δ ( ξ i , ξ j ) = 0 .
This case occurs when the combinations γ i = ( i 1 , i 2 ) and γ j = ( i ^ 1 , i ^ 2 ) are equal, which implies that there is no error in the spatial symbol. In this case, the two symbols in the transmitted vector and those in the estimated vector are located at the same index. In other words, the signal symbols in both vectors are transmitted from the same antennas. Accordingly, ( i 1 = i ^ 1 , i 2 = i ^ 2 ) and θ i = θ j . Based on the first and second columns of Table 1, the frequency of this term is equal to 4. The Ω term is given by:
Ω 1 = k 1 L k 2 L k ^ 1 L k ^ 2 L | s k 1 s k ^ 1 | 2 + | s k 2 s k ^ 2 | 2 n R
Therefore, the error probability will depend on the Euclidean distance among the signal symbols transmitted from the same antenna.
Case 2: δ ( ξ i , ξ j ) = 2 .
In this case, one signal symbol from the transmitted vector and one in the estimated vector are located at the same index. The second symbols in the transmitted and estimated vectors are located at different indexes. This corresponds to one of the following equiprobable four cases:
  • i 1 = i ^ 1 , i 2 i ^ 2 , γ 1 = ( 1 , 3 ) γ 3 = ( 1 , 4 ) and f 2 = 2 . The corresponding Ω term is given by:
    Ω 2 = k 1 L k 2 L k ^ 1 L k ^ 2 L | s k 1 e j θ 1 s k ^ 1 e j θ 3 | 2 + | s k 2 | 2 + | s k ^ 2 | 2 n R
  • i 1 = i ^ 1 , i 2 i ^ 2 , γ 2 = ( 2 , 4 ) γ 4 = ( 2 , 3 ) and f 3 = 2 . The corresponding Ω term is given by:
    Ω 3 = k 1 L k 2 L k ^ 1 L k ^ 2 L | s k 1 e j θ 2 s k ^ 1 e j θ 4 | 2 + | s k 2 | 2 + | s k ^ 2 | 2 n R
  • i 1 i ^ 1 , i 2 = i ^ 2 , γ 1 = ( 1 , 3 ) γ 4 = ( 2 , 3 ) and f 4 = 2 , leading to
    Ω 4 = k 1 L k 2 L k ^ 1 L k ^ 2 L | s k 2 e j θ 1 s k ^ 2 e j θ 4 | 2 + | s k 1 | 2 + | s k ^ 1 | 2 n R
  • i 1 i ^ 1 , i 2 = i ^ 2 , γ 2 = ( 2 , 4 ) γ 3 = ( 1 , 4 ) and f 5 = 2 . The Ω term is given by:
    Ω 5 = k 1 L k 2 L k ^ 1 L k ^ 2 L | s k 2 e j θ 2 s k ^ 2 j θ 3 | 2 + | s k 1 | 2 + | s k ^ 1 | 2 n R
Case 3: δ ( ξ i , ξ j ) = 4 .
In this case, the two signal symbols in the transmitted vector and those in the estimated vector are located at different indexes. The two possibilities leading to δ ( ξ i , ξ j ) = 4 are: i 1 i ^ 1 , i 2 i ^ 2 , γ 1 = ( 1 , 3 ) γ 2 = ( 2 , 4 ) , γ 3 = ( 1 , 4 ) γ 4 = ( 2 , 3 ) and f 6 = 4 . The Ω term is given by:
Ω 6 = k 1 L k 2 L k ^ 1 L k ^ 2 L | s k 1 | 2 + | s k ^ 1 | 2 + | s k 2 | 2 + | s k ^ 2 | 2 n R
Note that the two signal symbols are drawn from the same PSK/QAM constellation. Based on (15)–(20), Ω 1 and Ω 6 are excluded from the cost function of the optimization problem because they do not depend on the rotation angles. Also, the weights associated with each of the terms Ω 2 to Ω 5 are equal. Therefore, they can be removed from the multi-objective optimization without affecting the end result. However, they should be used in the evaluation of error performance as they might dominate the sum in (13). Finding the set of angles that reduces the error probability is formulated as a multi-objective optimization problem as follows:
θ * = arg min 0 θ i π / 2 , i = 1 , , 4 { Ω 2 + Ω 3 + Ω 4 + Ω 5 }
where θ is a vector of length of 4. The optimal angles are obtained using the Matlab optimization tool. First, we analyze the obtained angles in the case of PSK modulation. The general form of the optimal angles obtained using (21) is given by:
| θ 1 θ 3 | = | θ 2 θ 4 | = π L , with θ 1 = θ 2 , θ 3 = θ 4
To verify the formula in (22), a three-dimensional graph is generated for four different PSK modulation schemes and illustrated in Figure 1. The value of the cost function depends on the difference between the angles rather than on individual values. As presented in Figure 1, we can conclude that the difference in the value of the cost function in (21) using the optimal angles and non-optimal angles vanishes as the modulation order increases. Table 2 depicts the minimum value of the cost function when the optimal angles are used and the maximum value of the cost function assuming θ i = 0 , i for several modulation schemes. The cost values are obtained using (21) by substituting the corresponding values of the rotation angles. For modulation orders 16 or higher, the value of the cost function becomes independent of the rotation angles. As shown in Figure 1, the change in the function’s value (y-axis value) is marginal for the case of 16-PSK.
The error performance of the MA-SM, assuming the whole set of Ω terms, is considered and analyzed using PSK modulation. Ω 6 is simplified as follows:
Ω 6 = L 4 4 n R
Figure 2 depicts the weighted values of Ω 2 Ω 5 , weighted by their corresponding frequency f, on the left y-axis and the value of Ω 1 weighted by f 1 on the right y-axis. The x-axis represents the iteration number of the optimization process given in (21). The 8-PSK modulation is used for the generation of both graphs. The number of receive antennas is equal to 3 and 5 in the left and right sub-figure, respectively. Based on (13), the upper bound on the pairwise error probability depends on the weighted sum of the six Ω terms. When the optimization algorithm converges to the optimal solution for n R = 3 , f 1 Ω 1 is about 6094 ( f 1 =4, Ω 1 = 1523.5 , from (15)) and f 6 Ω 6 = 256 ( f 6 = 4 and Ω 6 = 64 , from (23)), and the sum of the terms that depend on the rotation angles is about 1180. Note that the values of Ω are obtained by substituting the values of the optimal rotation angles in (15)–(20) for any modulation set. Therefore, the error probability will be dominated by the terms that do not depend on the rotation angles. Further dominance is remarked as n R increases, where
f 1 Ω 1 + f 6 Ω 6 i = 2 5 f i Ω i 75.4 ,
From the previous analysis and for n R = 5 , we conclude that:
  • For a given modulation order, the effect of the rotation angle on the upper bound on the pairwise error probability vanishes as n R increases.
  • For a given n R , the effect of the rotation angle on the upper bound on the pairwise error probability vanishes as the modulation order increases.
It is therefore expected that the change in the rotation angles would have minimal impact, if any, on the error performance of the MA-SM system. The error performance is verified in Section 5 using simulations.
The optimization of the rotation angles is analyzed for the case when both signal symbols are drawn from the same L-QAM constellation with L > 4 . Based on the result of the optimization process, the optimal rotation angle is given as:
| θ 1 θ 3 | = | θ 2 θ 4 | = π 4 , with θ 1 = θ 2 , θ 3 = θ 4
Equation (22) is satisfied only for n R = 2 or 3. For n R = 4 and 6, we remarked that | θ 1 θ 3 | = | θ 2 θ 4 | 65 leads to a better performance than values in (22). Several other sets of angles that are not symmetric achieve similar performances.
Figure 3 depicts the weighted Ω terms that depend on the rotation angles on the left y-axis, and the weighted Ω 1 on the right y-axis, versus the iteration number of the multi-objective optimization process. The value of Ω 6 is neglected as it is much smaller than Ω 1 . For instance, Ω 1 / Ω 6 approximately equals 72 for 64-QAM. Based on the depicted graphs, we can draw the same conclusion of PSK modulation that the error probability will be dominated by the terms that do not depend on the rotation angles. Therefore, the rotation angles would have negligible effect on the error performance. Table 3 presents a comparison between the value of the cost function in (21) using the optimal rotation angles suggested in [20] and the conventional system without rotation for 16 and 64-QAM schemes. For both scenarios, the optimal rotation angles minimize the cost function.

4.3. Analysis of the MA-SM Performance for Configuration III

Based on the last configuration in Table 1, the size of the set Γ is equal to 8. Therefore, there are 64 codeword error events. The event γ i γ i , where the spatial symbol is received without error, is repeated eight times. This leads to a single Ω term, denoted by Ω 1 in Table 4. Note that Ω 1 is the only term that does not depend on the rotation angles. Also, since the probability of the error event γ i γ j is equal to γ j γ i , there are 28 other Ω terms, each is repeated twice. Out of the 28 terms, 19 have a similar form to Ω 2 , where there are two terms of the squared Euclidean distance between two rotated signal symbols and two power terms of a signal symbol. The remaining nine terms include a single term of the squared Euclidean distance between two rotated signal symbols and four power terms, as in Ω 8 . Table 4 gives a partial list of the Ω expressions.
In Table 4, the term with more occurrences of the rotation angles has more impact on the end result of the optimization process. For instance, Ω 8 is a function of only two rotation angles each with a single frequency. On the other hand, Ω 2 is a function of θ 1 and θ 2 , where each angle occurred twice. We can therefore conclude that Ω 2 has more impact on the obtained optimal rotation angle that does Ω 8 .
Table 5 shows the obtained optimal rotation angles for QPSK and 16-QAM and several values of receive antennas. The values of the cost function (the sum term in (13) with B = 29 , f 1 = 8 , and f i = 2 , for i = 2,⋯, 29) is also included using the optimal angles, the angles given in [20] and the all-zero angles. In all scenarios, the optimal rotation angles derived in this paper reduced the values of the cost function. We remarked that several sets of rotation angles produce almost the same value of the multi-objective cost function. Table 5 lists a single example. Unlike in the case of configuration I and II, it is difficult to obtain a closed-form expression for the rotation angles.

5. Simulation Results and Discussion

The effect of the rotation angles on the BER performance of the MA-SM system is evaluated for several system configurations. We adapted in this paper the same transmitter design in [20] in terms of the number of transmit antennas, the number of active transmit antennas at each channel use, the modulation schemes used for bit-to-symbol mapping, and the channel and noise models. The elements of the channel matrix and the noise vector are follow a cyclically symmetric Gaussian distribution with zero mean. The channel state information (CSI) is assumed to be perfectly known for the study of this paper. This is a general assumption in the literature when the channel estimation is not addressed in the provided work. This paper and [20] developed a baseband implementation of the communication system. We also assume that all signal symbols are drawn from the same modulation set.
Figure 4 depicts the performance of MA-SM with N T = 4 , N P = 2 and several values of n R . QPSK and 16-QAM are used to generate the results depicted in the left and right sub-figure, respectively. For both cases, we remarked that the performances of MA-SM using the optimal angles, the angles proposed in [20], and with the angles set to zero (without rotation) are identical for n R = 3 and 4. In the case of n R = 6 , the BER performance of MA-SM using the optimal angles slightly outperforms the performance of the other two angle sets, especially in the case of QPSK. However, this improvement is marginal and barely noticeable.
Figure 5 depicts the performance of MA-SM with N T = 5 , N P = 3 and several values of n R . The performance of MA-SM using the three angle sets coincides for n R = 3 . In the case of n R = 4 , a slight improvement in the BER performance is remarked in the case of QPSK. At n R = 6 , the performance of MA-SM with QPSK is improved due to the optimal rotation angle, and a slight improvement is also observed when 16-QAM is used. Again, the improvements in the BER using the 3-dimensional constellation are marginal and the MA-SM system can be used without applying any rotation to of the combinations while keeping the BER performance identical to that when the optimal rotation angles are applied. This paper developed a baseband implementation of the communication system. As a future work, we would like to implement the system using software-defined radio. This will permit the evaluation of the system under more realistic environment and taking into consideration system impairments.

6. Conclusions

The goal of this work is to validate the application of the 3-dimensional constellation that was introduced in the literature for the generalized SM with multiple active antennas (MA-SM). In this paper, we derived closed-form upper bound formulae of the pair-wise error probability of the MA-SM system as a function of the rotation angles used in the 3-dimensional constellation. To this end, the search problem is modeled as a multi-objective optimization. The obtained optimal set of angles are then evaluated against those which were proposed in the literature and the MA-SM system without the 3-dimensional constellation. We finally showed, using analytical and simulation results, that only a slight improvement of less than 0.5 dB in the error performance is achieved using the optimal angles. Therefore, we do not recommend the use of 3-dimensional constellation with the MA-SM system.

Author Contributions

Conceptualization, S.M., M.M., K.-S.C.; methodology, S.M., T.H., N.M., M.M. and K.-S.C.; software, S.M., T.H. and N.M.; validation, S.M., T.H., N.M., M.M. and K.-S.C.; formal analysis, M.M. and S.M.; investigation, S.M., T.H., N.M., M.M. and K.-S.C.; resources, K.-S.C. and M.M.; data curation, S.M., T.H. and N.M.; writing–original draft preparation, S.M., M.M. and K.-S.C.; writing–review and editing, S.M., T.H., N.M., M.M. and K.-S.C.; visualization, S.M., T.H. and N.M.; supervision, K.-S.C. and M.M.; project administration, M.M.; funding acquisition, K.-S.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the BK21 plus program through the National Research Foundation (NRF) funded by the Ministry of Education of Korea.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
BERBit-error rate
CQSMComplex quadrature spatial modulation
CSIChannel state information
GSMGeneralized spatial modulation
IMIndex modulation
IQSMIndex quadrature spatial modulation
MA-SMMultiple active antennas spatial modulation
MIMOMultiple input multiple output
MLMaximum-likelihood
PEPPair-wise error probability
PSKPhase shift keying
QAMQuadrature amplitude modulation
QPSKQuadrature phase shift keying
QSMQuadrature spatial modulation
RFRadio frequency
SMSpatial modulation
SMuxSpatial multiplexing
SSKSpace shift keying
SNRSignal-to-noise ratio
UPEPUnconditional pairwise error probability

References

  1. Ishikawa, N.; Sugiura, S.; Hanzo, L. 50 years of permutation, spatial and index modulation: From classic RF to visible light communications and data storage. IEEE Commun. Surv. Tutor. 2018, 20, 142–149. [Google Scholar] [CrossRef] [Green Version]
  2. Mesleh, R.; Haas, H.; Sinanovic, S.; Ahn, C.; Yun, S. Spatial modulation. IEEE Trans. Veh. Technol. 2008, 57, 2228–2241. [Google Scholar] [CrossRef]
  3. Jeganathan, J.; Ghrayeb, A.; Szczecinski, L.; Ceron, A. Space Shift Keying Modulation for MIMO Channels. IEEE Trans. Wirel. Commun. 2009, 8, 3692–3703. [Google Scholar]
  4. Kaddoum, G.; Ahmed, M.F.A.; Nijsure, Y. Code Index Modulation: A High Data Rate and Energy Efficient Communication System. IEEE Commun. Lett. 2015, 19, 175–178. [Google Scholar]
  5. Zhang, J.; Wang, Y.; Zhang, J.; Ding, L. Polarization Shift Keying (PolarSK): System Scheme and Performance Analysis. IEEE Trans. Veh. Technol. 2017, 66, 10139–10155. [Google Scholar] [CrossRef]
  6. Basar, E.; Aygolu, U.; Panayirci, E.; Poor, H.V. Orthogonal Frequency Division Multiplexing With Index Modulation. IEEE Trans. Signal Process. 2013, 61, 5536–5549. [Google Scholar] [CrossRef]
  7. Zhu, X.; Wang, Z.; Wang, Q.; Haas, H. Virtual Spatial Modulation for MIMO Systems. In 2016 IEEE Global Communications Conference (GLOBECOM); IEEE: Washington, DC, USA, 2016; pp. 1–6. [Google Scholar]
  8. Li, J.; Wen, M.; Zhang, M.; Cheng, X. Virtual Spatial Modulation. IEEE Access 2016, 4, 6929–6938. [Google Scholar] [CrossRef]
  9. Younis, A.; Mesleh, R.; Di Renzo, M.; Haas, H. Generalised spatial modulation for large-scale MIMO. In Proceedings of the 2014 22nd European Signal Processing Conference (EUSIPCO), Lisbon, Portugal, 1–5 September 2014; pp. 346–350. [Google Scholar]
  10. Ju, P.; Zhang, M.; Cheng, X.; Wang, C.X.; Yang, L. Generalized spatial modulation with transmit antenna grouping for correlated channels. In Proceedings of the 2016 IEEE International Conference on Communications (ICC), Kuala Lamput, Malaysia, 22–27 May 2016. [Google Scholar]
  11. AbuTayeh, S.; Alsalahat, M.; Kaddumi, I.; Alqannas, Y.; Althunibat, S.; Mesleh, R. A half-full transmit-diversity spatial modulation scheme. In Proceedings of the International Conference on Broadband Communications. Networks and Systems (BROADNETS), Faro, Portugal, 19–20 September 2018; pp. 257–266. [Google Scholar]
  12. Mohaisen, M. Constellation design and performance analysis of the parallel spatial modulation. Int. J. Commun. Syst. 2019, 32, e4165. [Google Scholar]
  13. Mesleh, R.; Ikki, S.S.; Aggoune, H.M. Quadrature Spatial Modulation. IEEE Trans. Veh. Technol. 2015, 64, 2738–2742. [Google Scholar]
  14. Vo, B.; Nguyen, H.H. Improved Quadrature Spatial Modulation. In 2017 IEEE 86th Veh. Technol. Conference (VTC-Fall); IEEE: Toronto, ON, USA, 2017; pp. 1–5. [Google Scholar]
  15. Holoubi, T.; Murtala, S.; Muchena, N.; Mohaisen, M. On the performance of improved quadrature spatial modulation. ETRI J. 2020. [Google Scholar] [CrossRef]
  16. Yigit, Z.; Basar, E. Double Spatial Modulation: A High-Rate Index Modulation Scheme for MIMO Systems. In 2016 International Symposium on Wireless Communication Systems (ISWCS); IEEE: Poznan, Poland, 2016; pp. 347–351. [Google Scholar]
  17. Mohaisen, M.; Lee, S. Complex Quadrature Spatial Modulation. ETRI J. 2017, 39, 514–524. [Google Scholar] [CrossRef]
  18. Mohaisen, M. Increasing the Minimum Euclidean Distance of the Complex Quadrature Spatial Modulation. IET Commun. 2018, 12, 854–860. [Google Scholar] [CrossRef]
  19. Iqbal, A.; Mohaisen, M.; Kwak, K.S. Modulation Set Optimization for the Improved Complex Quadrature SM. Wirel. Commun. Mobile Comput. 2018, 1–12. [Google Scholar] [CrossRef]
  20. Wang, J.; Jia, S.; Song, J. Generalised Spatial Modulation System with Multiple Active Transmit Antennas and Low Complexity Detection Scheme. IEEE Trans. Wirel. Commun. 2012, 11, 1605–1615. [Google Scholar] [CrossRef]
  21. Altın, G.; Çelebi, M.E. A Simple Low-Complexity Algorithm for Generalized Spatial Modulation. AEU Int. J. Electron. Commun. 2018, 97, 63–67. [Google Scholar] [CrossRef]
  22. Zheng, B.; Wen, M.; Chen, F.; Huang, N.; Ji, F.; Yu, H. The K-Best Sphere Decoding for Soft Detection of Generalized Spatial Modulation. IEEE Trans. Commun. 2017, 65, 4803–4816. [Google Scholar] [CrossRef]
Figure 1. Three-dimensional graphs of the cost function in (21) as a function of θ 1 and θ 3 for BPSK, QPSK, 8-PSK, and 16-PSK modulation schemes and n R = 3 .
Figure 1. Three-dimensional graphs of the cost function in (21) as a function of θ 1 and θ 3 for BPSK, QPSK, 8-PSK, and 16-PSK modulation schemes and n R = 3 .
Applsci 10 03718 g001
Figure 2. The value of the weighted Ω terms that are affected by the rotation angle (left y-axis) and the value of the rotation-independent Ω 1 (right y-axis) for n R = 3 (left sub-figure) and n R = 5 (right sub-figure), assuming 8-PSK.
Figure 2. The value of the weighted Ω terms that are affected by the rotation angle (left y-axis) and the value of the rotation-independent Ω 1 (right y-axis) for n R = 3 (left sub-figure) and n R = 5 (right sub-figure), assuming 8-PSK.
Applsci 10 03718 g002
Figure 3. The value of the weighted Ω terms that depend on the rotation angle (left y-axis) and the value of the rotation-independent Ω 1 (right y-axis) for n R = 3 and 16-QAM (left sub-figure) and 64-QAM (right sub-figure).
Figure 3. The value of the weighted Ω terms that depend on the rotation angle (left y-axis) and the value of the rotation-independent Ω 1 (right y-axis) for n R = 3 and 16-QAM (left sub-figure) and 64-QAM (right sub-figure).
Applsci 10 03718 g003
Figure 4. BER of MA-SM using the optimal rotation angles, the angles derived in [20], and the all-zero angle set, for N T = 4 , N P = 2 , n R = 3, 4, and 6, using QPSK (left sub-figure) and 16-QAM (right-sub-figure).
Figure 4. BER of MA-SM using the optimal rotation angles, the angles derived in [20], and the all-zero angle set, for N T = 4 , N P = 2 , n R = 3, 4, and 6, using QPSK (left sub-figure) and 16-QAM (right-sub-figure).
Applsci 10 03718 g004
Figure 5. BER of MA-SM using the optimal rotation angles, the angles derived in [20], and the all-zero angle set, for N T = 5 , N P = 3 , n R = 3, 4, and 6, using QPSK (left sub-figure) and 16-QAM (right-sub-figure).
Figure 5. BER of MA-SM using the optimal rotation angles, the angles derived in [20], and the all-zero angle set, for N T = 5 , N P = 3 , n R = 3, 4, and 6, using QPSK (left sub-figure) and 16-QAM (right-sub-figure).
Applsci 10 03718 g005
Table 1. MA-SM system parameters.
Table 1. MA-SM system parameters.
Configuration IConfiguration IIConfiguration III
N T 445
N P 223
ModulationQPSK16-QAM16-QAM
Γ (1,3),(1,3),(1,2,3), (1,2,4),
(2,4),(2,4),(1,2,5), (1,3,4),
(1,4),(1,4),(1,3,5), (1,4,5),
(2,3)(2,3)(2,3,5), (2,4,5)
Rotation angle θ 1 = 0 , θ 2 = π / 8 , θ 1 = 0 , θ 2 = π / 8 , θ 1 = 0 , θ 2 = π / 16 ,
θ 3 = π / 4 , θ 4 = 3 π / 8 θ 3 = π / 4 , θ 4 = 3 π / 8 ⋯, θ 8 = 7 π / 16
Table 2. The value of the cost function in (21) for optimal angles and θ i = 0 for several PSK modulation schemes and n R = 3 .
Table 2. The value of the cost function in (21) for optimal angles and θ i = 0 for several PSK modulation schemes and n R = 3 .
ModulationQPSK8-PSK16-PSK32-PSK
θ optimal θ i = 0 optimal θ i = 0 optimal θ i = 0 optimal θ i = 0
Cost32.841.2590.2592.29459.39459.31.5135 × 10 5 1.5135 × 10 5
Table 3. A comparison between the minimum value, maximum value, and the one corresponding to the angles in [20] of the cost function in (21) for 16 and 64-QAM with n R = 3 .
Table 3. A comparison between the minimum value, maximum value, and the one corresponding to the angles in [20] of the cost function in (21) for 16 and 64-QAM with n R = 3 .
Modulation16-QAM64-QAM
θ optimal θ i = 0  [20]optimal θ i = 0  [20]
Cost3119535883313131.4027 × 10 7 1.5975 × 10 7 1.4225 × 10 7
Table 4. Selected partial list of the Ω terms obtained for configuration III.
Table 4. Selected partial list of the Ω terms obtained for configuration III.
Error EventExpression
γ i γ i Ω 1 = | s k 1 s k ^ 1 | 2 + | s k 2 s k ^ 2 | 2 + | s k 3 s k ^ 3 | 2 n R
γ 1 γ 2 Ω 2 = | s k 1 e j θ 1 s k ^ 1 e j θ 2 | 2 + | s k 2 e j θ 1 s k ^ 2 e j θ 2 | 2 + | s k 3 | 2 + | s k ^ 3 | 2 n R
γ 1 γ 8 Ω 8 = | s k 2 e j θ 1 s k ^ 1 e j θ 8 | 2 + | s k 1 | 2 + | s k 3 | 2 + | s k ^ 2 | 2 + s k ^ 3 | 2 n R
γ 2 γ 4 Ω 10 = | s k 1 e j θ 2 s k ^ 1 e j θ 4 | 2 + | s k 3 e j θ 2 s k ^ 3 e j θ 4 | 2 + | s k 2 | 2 + | s k ^ 2 | 2 n R
γ 2 γ 8 Ω 12 = | s k 2 e j θ 2 s k ^ 1 e j θ 8 | 2 + | s k 3 e j θ 2 s k ^ 2 e j θ 8 | 2 + | s k 1 | 2 + | s k ^ 3 | 2 n R
γ 3 γ 6 Ω 16 = | s k 1 e j θ 3 s k ^ 1 e j θ 6 | 2 + | s k 3 e j θ 3 s k ^ 3 e j θ 6 | 2 + | s k 2 | 2 + | s k ^ 2 | 2 n R
γ 3 γ 8 Ω 19 = | s k 2 e j θ 3 s k ^ 1 e j θ 8 | 2 + | s k 3 e j θ 3 s k ^ 3 e j θ 8 | 2 + | s k 1 | 2 + | s k ^ 2 | 2 n R
γ 4 γ 7 Ω 21 = | s k 2 e j θ 4 s k ^ 2 e j θ 7 | 2 + | s k 1 | 2 + | s k 3 | 2 + | s k ^ 1 | 2 + s k ^ 3 | 2 n R
γ 5 γ 6 Ω 24 = | s k 1 e j θ 5 s k ^ 1 e j θ 6 | 2 + | s k 3 e j θ 5 s k ^ 3 e j θ 6 | 2 + | s k 2 | 2 + | s k ^ 2 | 2 n R
γ 6 γ 8 Ω 27 = | s k 2 e j θ 6 s k ^ 2 e j θ 8 | 2 + | s k 3 e j θ 6 s k ^ 3 e j θ 8 | 2 + | s k 1 | 2 + | s k ^ 1 | 2 n R
γ 7 γ 8 Ω 29 = | s k 1 e j θ 7 s k ^ 1 e j θ 8 | 2 + | s k 3 e j θ 7 s k ^ 3 e j θ 8 | 2 + | s k 2 | 2 + | s k ^ 2 | 2 n R
Table 5. Examples of the obtained optimal rotation angles (in radians) for QPSK and 16-QAM and several n R values. The values cost (opt), cost ([20]), and cost (zeros) are the value of the sum in (13) using the obtained optimal rotation angles, the angles given in [20], and the all-zero angles, respectively.
Table 5. Examples of the obtained optimal rotation angles (in radians) for QPSK and 16-QAM and several n R values. The values cost (opt), cost ([20]), and cost (zeros) are the value of the sum in (13) using the obtained optimal rotation angles, the angles given in [20], and the all-zero angles, respectively.
n R = 3 n R = 4 n R = 6
QPSK16-QAMQPSK16-QAMQPSK16-QAM
θ 1 1.17711.46670.43310.90550.63230.4048
θ 2 0.63870.81410.97590.25740.08171.1770
θ 3 0.03710.35230.01.36681.07120.0069
θ 4 0.03600.30561.57081.36781.06700.0002
θ 5 0.46030.72711.14480.15051.49850.7996
θ 6 1.00351.22250.59661.01500.48380.3971
θ 7 0.79851.08420.80760.44010.26641.1897
θ 8 1.55190.06880.05910.72781.00310.7830
cost (opt)2246.41.7022 × 10 7 628.61.0412 × 10 7 65.91.3906 × 10 7
cost ([20]2344.21.7218 × 10 7 695.51.0720 × 10 7 88.91.6092 × 10 7
cost (zeros)2674.61.8555 × 10 7 920.71.4343 × 10 7 168.34.3682 × 10 7

Share and Cite

MDPI and ACS Style

Murtala, S.; Holoubi, T.; Muchena, N.; Mohaisen, M.; Choi, K.-S. On the Performance of the Multiple Active Antenna Spatial Modulation with 3-Dimensional Constellation. Appl. Sci. 2020, 10, 3718. https://doi.org/10.3390/app10113718

AMA Style

Murtala S, Holoubi T, Muchena N, Mohaisen M, Choi K-S. On the Performance of the Multiple Active Antenna Spatial Modulation with 3-Dimensional Constellation. Applied Sciences. 2020; 10(11):3718. https://doi.org/10.3390/app10113718

Chicago/Turabian Style

Murtala, Sheriff, Tasnim Holoubi, Nishal Muchena, Manar Mohaisen, and Kang-Sun Choi. 2020. "On the Performance of the Multiple Active Antenna Spatial Modulation with 3-Dimensional Constellation" Applied Sciences 10, no. 11: 3718. https://doi.org/10.3390/app10113718

APA Style

Murtala, S., Holoubi, T., Muchena, N., Mohaisen, M., & Choi, K. -S. (2020). On the Performance of the Multiple Active Antenna Spatial Modulation with 3-Dimensional Constellation. Applied Sciences, 10(11), 3718. https://doi.org/10.3390/app10113718

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