Abstract

This paper concentrates on the designing of a robust Tomlinson-Harashima Precoder (THP) over multiple-input multiple-output (MIMO) channels in wireless communication systems with assumption of imperfect channel state information (CSI) at the transmitter side. With the assumption that the covariance matrix of channel estimation error is available at the transmitter side, we design a THP that presents robustness against channel uncertainties. In the proposed robust THP, the transmit power is further minimized by using the Tilted constellation concept. This power minimization reduces the interchannel Interference (ICI) between subchannels and, furthermore, recovers some part of the THP's power loss. The bit error rate (BER) of the proposed system is further improved by using a power loading technique. Finally, the simulation results compare the performance of our proposed robust THP with a conventional MIMO-THP.

1. Introduction

Recently, due to the impressive capacity of multiple-input multiple-output (MIMO) channel, many companies and researches have been attracted toward designing transceivers for MIMO systems and especially robust transceivers in the case of channel uncertainties. It is well known that inter channel interference (ICI) is one of the important problems of the MIMO systems. It is possible to preeliminate the ICI at the transmitter side by using precoding techniques and consequently reduce transmission power. Due to channel estimation errors, channel variations during time, and channel feedback errors, the assumption of perfect CSI at the transmitter side is not always true and hence the CSI at the transmitter side has uncertainties [1]. Due to these uncertainties, the ICI will be increased and consequently leads to BER growth. Linear precoder in [2] is designed for imperfect CSI, but it could only slightly reduce uncertainties effects by using statistical information of the channel. Nonlinear Tomlinson-Harashima Precoder (THP) in [3] has been considered for imperfect CSI assumption, and we combine it with the tilted constellation method in our work.

In this paper, we propose a design of a robust MIMO-THP that can reduce uncertainties effects more than the proposed methods in [2, 3] and subsequently can further reduce BER. Similar to [3], we utilize second-order statistics of the uncertainties of the MIMO channels at the transmitter side. Moreover, by using tilted constellation concept in [4], we can significantly reduce transmit power. It should be noted that transmit power reduction can decrease ICI between subchannels and in addition recover some parts of the THP’s power loss. This ICI reduction between subchannels is one of the main challenges in the designing of our precoder.

Tilted constellation method is introduced in [4] for SISO channel. Based on this idea, we have extended it for the MIMO channel. Here, the power reduction is done by tilting or rotating the ordinary constellation through a proper angle. The desired angle is one that minimizes the transmit power and is selected from a set of possible angles. The results show that the power reduction will be greater at high SNRs values.

To do the tilted constellation method, we use different angles for tilting the constellation in each antenna. Then, we find a set of proper angles, which are selected in a manner to minimize the transmit power in all antennas. The transmit power reduction can reduce the transmission block power, lower than the power of modulated symbols in each block. Finally, we can send proper angles to the receiver within the main data frame. We will show that the proposed method outperforms the conventional MIMO-THP in most SNRs and also in all symbol block lengths.

We can adjust transmit power of each antenna when the CSI (perfect/imperfect) is available at the transmitter side which is known as power loading [5]. It should be noted that in the MIMO systems some of the corresponding parallel subchannels (corresponding to each antenna) might have very low channel condition or might be useless for sometimes [5]. In this situation, the transmitter can adjust transmit power of each antenna by redistributing the available transmit power to get a better average error rate [5]. Moreover, by utilizing power loading, we improve the BER performance of the tilted MIMO-THP. We also study the complexity of our proposed method with respect to some previous methods.

This paper is organized as follows. In Section 2, the system model of the Tilted MIMO-THP in imperfect CSI scenario is introduced. In Section 3, the proposed robust tilted MIMO-THP is described for imperfect CSI. Section 4 presents the power loading technique to improve the performance of the robust precoder. In Section 5, the simulation results are given to show the performance of the proposed method and its comparison with the conventional MIMO-THP. Finally, conclusions are drawn in Section 6.

2. System Model

In this paper, an MIMO communication system with transmit antenna at the base station and users at the receivers side, each equipped with one receive antenna, are considered. The block diagram of the considered system is shown in Figure 1.

The input data symbols and the transmitted signals are denoted by vector and vector , respectively. Both vectors are of size . The channel matrix is modeled as , where is the estimated channel matrix whose elements are independent and identically distributed (i.i.d.) zero mean complex Gaussian random variables with variance , and is an channel uncertainties matrix whose elements are i.i.d. zero mean complex Gaussian random variables with variance . Thus, can be modeled as . The matrix is the proper rotation angles matrix which is named as tilting matrix, to minimize the transmission power. The received signal can be denoted as follows: where is the noise vector whose elements are assumed to be i.i.d. complex Gaussian random variables with zero mean and variance (i.e. .

As we see in Figure 1, the dimensional vector corresponds to input data, and dimensional vector represents transmitted signals. The feedback matrix is added to preeliminate the interference from previous users, and the tilting matrix is added to reduce the transmit power. Finally, in order to keep the power of transmitted symbols within the original constellation boundary, we employ modulo operator at the forward loop. The boundary region of the constellation is related to the order of modulation constellation, which for rectangular -ary QAM modulation is . Using tilted constellation, the signal constellation is further rotated by a set of appropriate angles, represented by the diagonal elements of the Tilting matrix , to reduce the ICI and power loss of THP.

The resultant signal is then passed through a unitary feed-forward filter , to eliminate the residual interference. Finally, the precoded signal is sent through the MIMO channel. All processes to eliminate the interference are performed at the transmitter side; hence, receivers at the user side are left with some simple operations including power scaling (i.e., elements of the diagonal matrix ), de-modulo operation, constellation tilting with reverse angle, and a single user detection.

By using the second-order statistics of the uncertainties of the MIMO channels at the transmitter side, we design a THP that presents robustness against channel uncertainties. We can describe the difference between impressive vector and input vector to the decision block as follows [3]:

MMSE solution should minimize the signal error as: where is the total transmitted power. Since, we assume the matrix is unitary thus, the power constraint is guaranteed.

Instead of analysis of (3) we can use orthogonality property as follow [3]:

Hence from (4) we have the following:

We can calculate and using (1) as follows:

By substituting (6) into (5), and after some manipulations, we have the following:

Since is unitary, we have the following: where . Matrix can be obtained by doing the Cholesky factorization of (8). We can calculate matrixes , , and as follows:

Substituting , , and into (9) the error covariance matrix is equal to the following: where the scalar represents the .

It can be seen in (10) that, while the channel has uncertainty, the error covariance matrix has an additional term of   with respect to its counterpart in perfect CSI case.

In this paper, we assume the previous relations and matrixes as the conventional robust MIMO-THP, that is, without using tilted constellations.

3. Robust Tilted MIMO-THP

As mentioned before, in the general model of Figure 1, each diagonal element of the Tilting matrix determines the rotation angle of its corresponding antenna signal.

In conventional THP, the main constellation is only extended by the modulo operator on its boundary regions. However, in the tilted-THP, the signal constellation is further rotated by a set of appropriate angles. To do this, first a block of symbols with length is divided to groups each of length . Then, for every group of symbols which is transmitted from each antenna, the transmitter chooses an appropriate angle from a set of possible angles in such a way that the transmitted power is minimized based on some predefined criteria. From [4], the ordinary constellation is tilted by possible angles as then the optimal angle is selected by the following equation: where is the group of symbols in each transmit antenna, and illustrates the interference sequence at each antenna due to the previous antennas and is calculated similar to the conventional THP. The difference between the conventional (i.e. un-tilted) and tilted THP is in rotating of the transmitted signal of the previous antennas by their optimal angles. Transmitted power and resulting optimal angle of each antenna can be calculated as follows [4]:

The vector , which is transmitted from the first transmit antenna, is not affected by any interference signal, that is why we transmit without any interference. The optimal tilted angle is resulted by minimizing the transmitted power of the first antenna; that is:

For , which is transmitted from the second antenna, the optimal transmit angle is achieved by computing the minimum transmitted power similar to the , but, here, there is an interference due to the , that is, , and hence we have the following:

Similarly, for the antenna, we have the following: where is the optimal tilted angle for the antenna. At the receiver side, the received signal for the antenna can be shown as follows: where denotes the row of the matrix . By using (15), (16) can be written as

Since , after modulo operator we obtain the following:

As shown in [3], in the tilted constellation case, the matrix is similar to the untilted or regular MIMO-THP case. Its difference is only in rotation of the phase of the noise. Since the noise has circularly symmetric complex Gaussian distribution, for the conventional case, all equations and results for the decomposition of matrices in the tilted constellation case are the same.

If we use the tilted constellation for detecting the original symbols directly, the performance will be degraded due to rotation and reduction of minimum distance between symbols. Hence, we should send the optimal angles to the receiver to detect the desired symbols. Upon receipt, first, the symbols are reversely rotated at each antenna and then are detected in a way similar to the conventional THP.

4. Power Loading

The SER of MIMO-THP in imperfect CSI can be approximated as follows [6]: where is the transmitted power of transmit antenna and is noise power [3]:

Here, , , and . With the assumption of small error, that is, [3], we can approximate the value as follows: where is the total tilted power. At the above equation, we assumed the worst-case (i.e., ). We distribute the power to minimize SER in imperfect CSI by using the Lagrange method, so that it can be written as follows:

Hence, as in [3], we have the following:

Unfortunately, we did not find any closed form solution for (23), so we should use some suboptimal solution as in [3].

By Solving for , we have the following: where , , , and is real Lambert function [3].

Thus, the power vector can be obtained by using the following unique repeated solution [3].(I)A small positive value is chosen so that (II)The value is calculated as follows: (III)If is not close enough to , multiply and be returned to the stage II.(IV)Based on (24), vector is calculated.

This algorithm is present as pseudocode in Algorithm 1.

   ; % initial value for
   = nT; % Total transmit power
  Calculate ,
Loop: Calculate
  If then end
  Else if , Break;
  Else , go to loop
  end

Now, we discuss the complexity of the proposed method. We realize that the proposed method adds some parts to the conventional THP method. To compare the complexities of the proposed method and the conventional THP, we only consider those parts that are changed in our method. For example, for the case of , the complexity of the conventional THP is at the transmitter side and is at the receiver side while for the proposed Tilted-THP method this complexity is at the transmitter side and is at the receiver side, where is the length of the input symbol to each antenna. As it can be seen, the increased complexity is suitable. On the other hand, as we will show in the next section in simulation results, the performance gain that we attain by the proposed method is more considerable with respect to the above complexity increasing.

5. Simulation Results

To evaluate the performance of our proposed method, we use some simulations. To do this we assume an MIMO-THP system with 4 transmit and 4 receive antennas with 4-QAM modulation. We assume that the channel estimation error covariance and 0.05. In our figures, the labels “conv,” “Rob,” and “P.L.” are corresponding to conventional, robust, and power loading methods, respectively.

Figure 2 shows the transmitted power of the MMSE-THP system versus the number of tilted angles. In this figure, the SNR is assumed to be 10 dB. It can be seen that by increasing the number of available tilted angles, the transmitted power reduces until it converges to a constant floor. Also, it can be seen that by decreasing the block length, the power reduction will be more and the power converges faster to its floor. The reason of this floor is that we assumed that the matrix is unitary. As can be seen, for example, the transmitted power for the conventional MIMO-MMSE-THP for is 4.91 but, for the Tilted constellation and with the same parameters, the transmitted power decreases to 4 and 4.24 when the block length increases to 5 and 10, respectively.

Figure 3 shows BER of the MMSE and robust-tilted-MMSE precoders versus the Eb/No for different block lengths and for with and without the power loading. As it is seen in this figure, with the Robust-Tilted THP, the BER decreases compared to the conventional one. This BER reduction is more when we use the power loading. Moreover, because of more power reduction in smaller block lengths, the BER reduction will be more in smaller block length. For example, for  dB and in the conventional MIMO-MMSE-THP, we have while in the robust Tilted constellation MIMO-MMSE-THP for with and without power loading we have and , respectively.

Figure 4 shows BER of the MMSE and Robust-Tilted-MMSE precoders versus the Eb/No for different block lengths and for with and without the power loading. As it is seen in this figure, with the Robust-Tilted constellation, BER decreased compared to the conventional MIMO-THP which is similar to Figure 3. Again, BER reduction is more when we use the power loading. Due to more power reduction of Tilted constellation for shorter block length, the resultant BER has reduced more for shorter block lengths. More channel uncertainties with respect to Figure 3 lead to BER incensement. For example, at in the conventional MIMO-MMSE-THP, BER = 0.011 while in the robust Tilted constellation MIMO-MMSE-THP for with and without power loading and , respectively.

6. Conclusion

In this paper, we proposed a robust-tilted-MIMO-THP scheme for MIMO channels based on MMSE criterion with the assumption that the imperfect CSI is available at the transmitter side. We used second-order terms of the uncertainties of the MIMO channels at the transmitter to design a robust precoder. In the proposed robust method, we minimized the transmit power by using the Tilted constellation. The power minimization will result in ICI reduction between subchannels and in addition recovering some of the THP’s power loss. Consequently, this ICI reduction and power recovering result led to BER improvement, especially at high SNRs. The transmission power is more reduced by utilizing more Tilted angles and smaller symbol block length. Also, we demonstrated that the achieved performance gain by using the proposed method is more considerable with respect to its higher complexity. Finally, we showed that we can improve our precoder performance further, by utilizing power loading.

Supplementary Materials

The proposed Tilted MIMO-THP in decentralized scenario.

The proposed power loading algorithm

Transmitted power vs. the number of Tilted angles in different block length N=5,10,20. (SNR=10dB, = nR =4, 4QAM). BER vs. SNR and several block length in MIMO-MMSE-THP with and without Power Loading (ρ=0.05). BER vs. SNR and several block length in MIMO-MMSE-THP with and without Power Loading (ρ=0.1).

  1. Supplementary Material