ABSTRACT

New transmission scheme based on space-time code is proposed for massive MIMO system, which consists of K users and a base station. Each user is equipped with four antennas and two of them are selected to send space-time codeword with Alamouti coded by using the method of generalized spatial modulation. The base station estimates the index of transmit antennas used by each user at first, and then uses an iterative method to separate each space- time codeword and decode modulated signals of each code-word. Compared with the existing scheme for the same scene, the proposed scheme reduces the complexity of decoding, while keeping the same transmission rate. Simulation results demonstrate that the proposed scheme outperforms the existing scheme in the same scene when the base station is equipped with large amount of antennas.

Subject Areas:

Information and Communication Theory and Algorithms

Keywords:

Massive MIMO Systems, Space-Time Code, Decoding Complexity, Generalized Spatial Modulation

1. Introduction

Transmission rate and spectrum efficiency are particularly important for multi- cell communication. Massive MIMO technology have the ability to improve spectrum efficiency, cut down transmission power, reduce intra cell interference and inter-cell interference of the system without increasing spectrum resources [1] . So, it becomes one of the breakthrough technologies in the future [2] .

In massive MIMO system, the base station is equipped with hundreds of antennas. K users transmit modulated signals simultaneously, or transmit the signals after precoding at the same time in a single cell uplink system [3] . If all of the users are equipped with $n_t$ antennas, there will be $n_t$ modulation symbols sent by each user per time slot [4] [5] . Generalized spatial modulation is studied in [6] , in which not only transmit antennas send symbols, but also different antenna combinations are mapped into several bits. Thus, the transmission efficiency is improved. However, the decoding complexity of Maximum Likelihood (ML) decoding, Message Passing (MP) decoding and Channel Hardening Exploiting Message Passing (CHEMP) in [6] is too high.

In order to reduce the decoding complexity, the transmission based on space- time coding is proposed in single-cell massive MIMO uplink system. All of the users are equipped with 4 antennas and generalized spatial modulation is used at each user. According to the information sequence, two of the antennas are selected to transmit the space-time code which has Alamouti structure. Firstly, the serial number of the transmit antenna used at each user are estimated at the base station. Then, each space-time codeword is separated using the orthogonality of Alamouti code. Finally, modulation symbols in each codeword are detected. The ML decoding complexity of the proposed scheme is proportional to the square of the modulation order. The simulation results show that the reliability of the proposed scheme is better than [6] .

2. System Model

The system model of the proposed scheme is shown in Figure 1, which contains K users and a base station. All users are equipped with 4 antennas and the base station is equipped with $N\left(N\ge 4K\right)$ antennas. Generalized spatial modulation proposed in [6] is adopted. Each user selects two antennas to transmit signals according to the information sequence. When the information sequences are 00, 01, 10 and 11, the corresponding antenna combinations are (1, 2), (1, 3), (1, 4) and (2, 3), respectively.

Rate-2 space-time code proposed in [7] is used at user k, $k=1,2,\cdots ,K$ . The codeword at user k, denoted by ${X^{\prime }}_k$ , is expressed as

${X^{\prime }}_k=\left[\begin{array}{c}{x}_1\\ x_2\end{array}\right]=\left[\begin{array}{cc}{x^{\prime }}_{k,1}& -{x^{\prime }}_{k,2}^{}{}^{\ast }\\ {x^{\prime }}_{k,2}& {x^{\prime }}_{k,1}^{}{}^{\ast }\end{array}\right]=\left[\begin{array}{cc}{x}_{k,1}+e^{j\theta }{x}_{k,2}& -x_{k,3}^{}{}^{\ast }-e^{-j\theta }{x}_{k,4}\\ x_{k,3}^k+e^{j\theta }{x}_{k,4}& x_{k,1}^{}{}^{\ast }+e^{-j\theta }{x}_{k,2}^{}{}^{\ast }\end{array}\right]$ (1)

where, $x_{k,\lambda }$ represents the modulation symbols, $\lambda =1,2,3,4$ , $j=\sqrt{-1}$ , ${(\cdot )}^{*}$ represents conjugate. The value of $\theta$ makes the elements of ${X^{\prime }}_k$ nonzero.

Assume that are the antenna numbers selected by user k, ${\mu }_k\ne {\nu }_k\in \left\{1,2,3,4\right\}$ ,. Then the antenna ${\mu }_k$ at user $k$ sends $x_1$ , the antenna ${\nu }_k$ sends $x_2$ and the other antennas send $\left[\begin{array}{cc}0& 0\end{array}\right]$ . Use $X_k$ to represent the codeword send by user $k$ . The dimension of $X_k$ is $4\times 2$ . The elements of $X_k$ in the row of ${\mu }_k$ and ${\nu }_k$ are $x_1$ and $x_2$ respectively, and the elements in other rows are zero.

The received signals at the base station, denoted by $Y$ , can be expressed as

(2)

where, dimension of $Y$ and $N$ are both $N\times 2$ . $N$ represents the complex Gauss noise matrix with zero mean and variance of ${\sigma }^2$ . $H_k$ , with dimension $N\times 4$ , represents the channel fading coefficient from user $k$ to the base station. The dimension of $x_{k,1}$ and $x_{k,2}$ are both $4\times 1$ . The ${\mu }_k$ -th and ${\nu }_k$ -th elements of $x_{k,1}$ are ${x^{\prime }}_{k,1}$ and, respectively. The ${\mu }_k$ -th element and ${\nu }_k$ -th element of $x_{k,2}$ are ${x^{\prime }}_{k,1}$ and ${x^{\prime }}_{k,2}$ , respectively. The other elements of $x_{k,1}$ and $x_{k,2}$ are zeros.

3. Detection

The detection at the base station can be divided into two steps. Firstly, the serial number of the transmit antennas at each user is estimated at the base station according to the received signals and the channel state information (CSI). Then, each codeword is separated and each element of the codeword is estimated using iterative method at the base station.

Use ${y^{\prime }}_1$ to represent the first column of $Y$ . Let $g=pinv\left(H\right){y^{\prime }}_1$ , which $pinv(\cdot )$ represents pseudo inverse. $g$ is a vector with a dimension of $4K\times 1$ . Use ${\overline{g}}_k$ to represent the vector with a dimension of $4\times 1$ , which consists of the $4\left(k-1\right)+1$ -th element to $4k$ -th emelent of $g$ , $k=1,2,\cdots ,K$ . According to the estimation method of transmitting antennas’ number in (2) and [8] , the serial number of the transmit antennas for user $k$ obtained at the base station can be expressed as

(3)

where, ${\mu }_{k0}\ne {\nu }_{k0}\in \left\{1,2,3,4\right\}$ , $g_{{{\mu }^{\prime }}_k}$ and $g_{{{\nu }^{\prime }}_k}$ represents the ${\mu }^{\prime }$ -th and ${\nu }^{\prime }$ -th element of ${\overline{g}}_k$ . $\Vert \cdot \Vert$ represents norm.

After the serial number of transmit antenna is estimated, the zero elements in (2) can be removed. Then, (2) can be equivalently expressed as

(4)

where, ${X^{\prime }}_k$ represents a $2\times 2$ matrix after removing the zero elements in $X_k$ , which is expressed as (1),. ${H^{\prime }}_k$ represents a $N\times 2$ matrix, which consists of ${\mu }_k$ -th column and ${\nu }_k$ -th column of $H_k$ .

Since ${X^{\prime }}_k$ has Alamouti structure, (4) can be expressed as

(5)

where $H_i^k=\left[\begin{array}{cc}{h}_{i,1}^k& h_{i,2}^k\\ h_{i,2}^k{}^{\ast }& -h_{i,1}^k{}^{\ast }\end{array}\right]$ ,, , $y_{i,j}$ represents the element of $Y$ in row $i$ and column $j$ ,, $j=1,2$ . ${(\cdot )}^{\text{T}}$ represents the transpose.

According to the method of separating Alamouti codewords in [9] and [10] , (5) can be decomposed into sub-equations, which $\lfloor \cdot \rfloor$ represents round down operation. Each sub-equation can be expressed as

(6)

Let

$\begin{array}{c}{G}_{\mathcal{l},j}\left(l\right)={\left[\begin{array}{cc}{I}_2& F_{\mathcal{l},j}{}^H\left(l\right)\end{array}\right]}^H\\ =\left[\begin{array}{cc}{I}_2& \underset{F_{\mathcal{l},j}^{}{}^H\left(l\right)}{\underbrace{-\frac{{\overline{H}}_{\mathcal{l},j}^{K-l+1}\left(l-1\right){\overline{H}}_{\mathcal{l},K-l+1}^{K-l+1}{}^H\left(l-1\right)}{{\Vert {\overline{H}}_{\mathcal{l},K-l+1}^{K-l+1}\left(l-1\right)\Vert }^2}}}\end{array}\right],\\ l\in \left\{1,\cdots ,K-1\right\},j\in \left\{1,\cdots ,K-l|l\right\}\end{array}$ (7)

where, ${(\cdot )}^H$ represents conjugate transpose, $I_2$ represents a unit matrix with a dimension of $2\times 2$ . From (4)-(8) in [10] , codewords can be separated gradually by $\left(K-1\right)$ times iteration. The matrix, vector and element with subscript $\mathcal{l}$ belong to the $\mathcal{l}$ -th sub-equation after decomposition, $\mathcal{l}=1,\cdots ,\Lambda$ .

From (25) in [10] , we can find that, the $\mathcal{l}$ -th sub-formula from (6) can be expressed as (8) after $\left(K-1\right)$ times iterative,

${\widehat{G}}_{\mathcal{l}}{}^H{\overline{y}}_{\mathcal{l}}={\widehat{G}}_{\mathcal{l}}{}^H\left[\begin{array}{c}{H}_{\left(\mathcal{l}-1\right)K+1}^{\text{1}}\\ ⋮\\ H_{\mathcal{l}K}^{\text{1}}\end{array}\right]s_1+{\widehat{G}}_{\mathcal{l}}{}^H{\overline{n}}_{\mathcal{l}}=\underset{{\tilde{H}}_{\mathcal{l}}}{\underbrace{{\widehat{G}}_{\mathcal{l}}{}^H\left[\begin{array}{cc}{\overline{h}}_{\mathcal{l},1}& {\overline{h}}_{\mathcal{l},2}\end{array}\right]}}{s}_1+{\widehat{G}}_{\mathcal{l}}{}^H{\overline{n}}_{\mathcal{l}}$ (8)

where, ${\widehat{G}}_{\mathcal{l}}=\left[\begin{array}{c}{A}_{\mathcal{l},1}\\ A_{\mathcal{l},2}\\ A_{\mathcal{l},3}\\ ⋮\\ A_{\mathcal{l},K}\end{array}\right]=\left[\begin{array}{c}{I}_2\\ F_{\mathcal{l},1}\left(K-1\right)\\ F_{\mathcal{l},1}\left(K-2\right)+{\phi }_{\mathcal{l}}{}^{\left(1\right)}\left(A_{\mathcal{l},2}\right)F_{\mathcal{l},1}\left(K-1\right)\\ ⋮\\ {\displaystyle \sum _{l=1}^{K-1}{\phi }_{\mathcal{l}}{}^{\left(K-l\right)}\left(A_{\mathcal{l},l}\right)F_{\mathcal{l},1}\left(l\right)}\end{array}\right]$ , ${\phi }_{\mathcal{l}}{}^{\left(n\right)}\left(A_{\mathcal{l},k}\right)$

represents a function replace all of $F_{\mathcal{l},\psi }\left(l\right)$ in $A_{\mathcal{l},k}$ to $F_{\mathcal{l},\psi +n}\left(l-n\right)$ , $n\in \left\{1,\cdots ,K-1\right\}$ , $\psi \in \left\{1,\cdots ,K-1\right\}$ . ${\phi }_{\mathcal{l}}{}^{\left(n\right)}\left(I_2\right)=I_2$ , $A_{\mathcal{l},k}$ is a matrix with dimension $2\times 2$ , ${\overline{h}}_{\mathcal{l},1}$ and ${\overline{h}}_{\mathcal{l},2}$ represent the 1-th column and 2-th column of

${\left[\begin{array}{cccc}{\left(H_{\left(\mathcal{l}-1\right)K+1}^1\right)}^H& {\left(H_{\left(\mathcal{l}-1\right)K+2}^1\right)}^H& \cdots & {\left(H_{\mathcal{l}K}^1\right)}^H\end{array}\right]}^H$ respectively.

Let ${\tilde{H}}_{\mathcal{l}}={\widehat{G}}_{\mathcal{l}}{}^H\left[\begin{array}{cc}{\overline{h}}_{\mathcal{l},1}& {\overline{h}}_{\mathcal{l},2}\end{array}\right]$ , and multiply both sides of (8) by ${\tilde{H}}_{\mathcal{l}}{}^H$ to get

${\tilde{H}}_{\mathcal{l}}{}^H{\widehat{G}}_{\mathcal{l}}{}^H{\overline{y}}_{\mathcal{l}}={\tilde{H}}_{\mathcal{l}}{}^H{\tilde{H}}_{\mathcal{l}}{s}_1+{\tilde{H}}_{\mathcal{l}}{}^H{\widehat{G}}_{\mathcal{l}}{}^H{\overline{n}}_{\mathcal{l}}$ (9)

Since the matrix with Alamouti structure is closed for addition, multiplication, and conjugate transpose operations, ${\tilde{H}}_{\mathcal{l}}$ has Alamouti structure, which result in ${\tilde{H}}_{\mathcal{l}}{}^H{\tilde{H}}_{\mathcal{l}}\text{=}\left({\overline{h}}_{1,1}{}^H{\widehat{G}}_1{\widehat{G}}_1{}^H{\overline{h}}_{1,1}+\cdots +{\overline{h}}_{\Lambda ,1}{}^H{\widehat{G}}_{\Lambda }{\widehat{G}}_{\Lambda }{}^H{\overline{h}}_{\Lambda ,1}\right)I_2$ . Use ${\tilde{r}}_1$ to represent the first element of ${\tilde{H}}_{\mathcal{l}}{}^H{\widehat{G}}_{\mathcal{l}}{}^H{\overline{y}}_{\mathcal{l}}$ , then (10) can be obtained by (6)-(9).

$\begin{array}{l}{\tilde{r}}_1=\left({\overline{h}}_{1,1}{}^H{\widehat{G}}_1{\widehat{G}}_1{}^H{\overline{h}}_{1,1}+\cdots +{\overline{h}}_{\Lambda ,1}{}^H{\widehat{G}}_{\Lambda }{\widehat{G}}_{\Lambda }{}^H{\overline{h}}_{\Lambda ,1}\right){x^{\prime }}_{1,1}\\ +\left({\overline{h}}_{1,1}{}^H{\widehat{G}}_1{\widehat{G}}_1{}^H{\overline{n}}_1+\cdots +{\overline{h}}_{\Lambda ,1}{}^H{\widehat{G}}_{\Lambda }{\widehat{G}}_{\Lambda }{}^H{\overline{n}}_{\Lambda }\right)\end{array}$ (10)

The base station can get ${x^{\prime }}_{1,1}$ from ${\tilde{r}}_1$ directly. The detection of the other modulation symbols is similar. Since ${X^{\prime }}_k$ has Alamouti structure and each codeword consists of two modulation symbols, the complexity of ML decoding is proportional to the square of the modulation order.

4. Bit Error Rate Analysis

Let $P_{e1}$ and $P_{e2}$ represents the error probability of estimating the number of transmit antennas and the symbol error rate (SER) of the second step. Then the SER of the system can be expressed as

$P_e=1-\left(1-P_{e1}\right)\left(1-P_{e2}\right)$ (11)

We deduce $P_{e2}$ firstly. Assume the mathematical expectation of the modulation signal is $P_1$ . The signal-noise ratio (SNR) during the process of decoding ${x^{\prime }}_{1,1}$ , denoted by $\rho$ , can be obtained from (10).

$\rho =\frac{P_1}{{\sigma }^2}\frac{{\left({\overline{h}}_{1,1}{}^H{\widehat{G}}_1{\widehat{G}}_1{}^H{\overline{h}}_{1,1}\right)}^2+\cdots +{\left({\overline{h}}_{\Lambda ,1}{}^H{\widehat{G}}_{\Lambda }{\widehat{G}}_{\Lambda }{}^H{\overline{h}}_{\Lambda ,1}\right)}^2}{{\overline{h}}_{1,1}^H{\widehat{G}}_1{\widehat{G}}_1{}^H{\widehat{G}}_1{\widehat{G}}_1{}^H{\overline{h}}_{1,1}+\cdots +{\overline{h}}_{\Lambda ,1}^H{\widehat{G}}_{\Lambda }{\widehat{G}}_{\Lambda }{}^H{\widehat{G}}_{\Lambda }{\widehat{G}}_{\Lambda }{}^H{\overline{h}}_{\Lambda ,1}}$ (12)

According to the estimation method of SER in [10] and [11] , the probability of estimating ${x^{\prime }}_{1,1}$ as ${{\widehat{x}}^{\prime }}_{1,1}$ , denoted by $P_{e2}$ , can be expressed as

where, $\mathbb{E}(\cdot )$ represents mathematical expectation. The value of $\left|{x^{\prime }}_{1,1}-{{\widehat{x}}^{\prime }}_{1,1}\right|$ depends on the modulation mode.

In what follows, we derive $P_{e1}$ . According to the consistent bound method proposed by (12) in [12] and [13] , $P_{e1}$ satisfies the following equation

$\begin{array}{c}{P}_{e1}\le \mathbb{E}\left[{\displaystyle \sum _{\widehat{\xi }}N\left(\xi ,\widehat{\xi }\right)P\left(x_0\to {\widehat{x}}_0\right)}\right]\\ \le {\displaystyle \sum _{\xi =1}^4{\displaystyle \sum _{q=1}^M{\displaystyle \sum _{\widehat{\xi }=1}^4\frac{N\left(\xi ,\widehat{\xi }\right)P\left(x_0\to {\widehat{x}}_0\right)}{4M}}}}\end{array}$ (14)

where, M represents the modulation order, $\xi$ represents the serial number of transmit antennas, $\widehat{\xi }$ is the number of antennas estimated by the base station. $x_0$ is the first column of $X_1$ . ${\widehat{x}}_0$ is the decoding results obtained by decoding the first column of $X_1$ . $N\left(\xi ,\widehat{\xi }\right)$ represents the number of error bits between $\xi$ and $\widehat{\xi }$ .

It is known from (18) in [12] and (5) in [14] ,

$P\left(x_0\to {\widehat{x}}_0\right)=f{\left(\alpha \right)}^N{\displaystyle \sum _{\beta =0}^{N-1}\left(\begin{array}{c}N-1+\beta \\ \beta \end{array}\right)}{\left(1-f\left(\alpha \right)\right)}^{\beta }$ (15)

where, $f\left(\alpha \right)=\frac{1}{2}\left(1-\sqrt{\frac{\alpha }{1+\alpha }}\right)$ , $\alpha =\frac{1}{4{\sigma }^2}{\displaystyle \sum _{\xi =1}^4{\left|x_{\xi }-{\tilde{x}}_{\xi }\right|}^2}$ , $x_{\xi }$ represent the $\xi$ -th element of $x_0$ , $\xi \in \left\{1,2,3,4\right\}$ . Substitute (15) into (14) to get

$P_{e1}\le {\displaystyle \sum _{\xi =1}^4{\displaystyle \sum _{q=1}^M{\displaystyle \sum _{\widehat{\xi }=1}^4\frac{N\left(\xi ,\widehat{\xi }\right)f{\left(\alpha \right)}^N{\displaystyle \sum _{\beta =0}^{N-1}\left(\begin{array}{c}N-1+\beta \\ \beta \end{array}\right)}{\left(1-f\left(\alpha \right)\right)}^{\beta }}{4M}}}}$ (16)

Take (13) and (16) into (11) to get the SER of the system, as shown in (17).

5. Simulation Results and Analysis

This chapter simulates the reliability of the proposed scheme and the reliability of [6] . Assume the channel follows the Rayleigh distribution and the noise is Gaussian white noise. There are 16 users using 4QAM modulation without channel coding.

The error probability $P_{e1}$ of estimating the transmit antenna number when the base station is equipped with 64 antennas is shown is Figure 2. Simulation results show that when the SNR is 22dB, the error probability is about ${10}^{-3}$ . Since the number of transmit antennas and the transmitted symbols are estimated simultaneously in [6] , $P_{e1}$ can’t be obtained. So the error probability by estimating the transmit antenna number in [6] is not shown in Figure 2. It can be seen that when N = 128, the $P_{e1}$ is 0 while SNR ≥ 0. The reason is that the more the number of receiving antennas, the more accurate of the judgment. Therefore, the error probability of the estimate of transmit antenna number is not plot in Figure 2 when N = 128.

The SER curves of the proposed scheme and the scheme in [6] with 64 antennas and 128 antennas is shown in Figure 3 and Figure 4 respectively. Simulation results show that when N = 64, the reliability of the proposed scheme is lower than [6] ，when N = 128, the reliability of the proposed scheme is better than [6] . That is because when N = 64 the $P_{e1}$ of the proposed scheme is higher, but when N = 128 the $P_{e1}$ of the proposed scheme is zero. From Figure 4, the gain of the proposed scheme is about 2.5 dB at the SER of ${10}^{-3}$ compared with the MP-GSM scheme, the gain is about 2.7 dB compared with the CHEMP-GSM scheme, and the gain is about 2.9 dB compared with the ML scheme.

6. Conclusion

Compared with the existing scheme, each user adopts rate-2 space-time code. The transmission efficiency is the same as the existing scheme. The proposed scheme uses the orthogonal characteristics of the equivalent channel matrix corresponding to the Alamouti codeword to gradually separate each codeword, and then decodes each symbol of the codewords. The decoding complexity is greatly reduced. However, the reliability of the proposed scheme is greatly dependent on the estimated of the transmission antenna. How to reduce the dependency of the reliability of the system on the estimated transmit antenna number requires further study.

Cite this paper

Tian, X.J., Jia, W.J. and Zhang, H.T. (2017) New Transmission Scheme Based on Space-Time Code for Massive MIMO. Open Access Library Jour- nal, 4: e3712. https://doi.org/10.4236/oalib.1103712

References