- Research
- Open Access
- Published:

# Optimization of a MIMO amplify-and- forward relay system with channel state information estimation error and feedback delay

*EURASIP Journal on Advances in Signal Processing*
**volume 2016**, Article number: 117 (2016)

## Abstract

This paper addresses the robust design of a multiple-input multiple-output amplify-and-forward relay system against channel state information (CSI) mismatch due to estimation error and feedback delay. The estimation error and feedback delay are expressed using appropriate models, from which we derive the conditional mean square error (MSE) between the desired and the received signals upon the estimated CSI. The conditional MSE is then minimized to optimize the relay beamforming matrix with relay transmission power constraint. It is shown that the proposed optimization problem reduces to the conventional minimum MSE problem when CSI mismatch vanishes. By analyzing the structure of the optimal beamforming matrix, the optimization problem is simplified so that it can be directly solved using the genetic algorithm (GA). To further reduce the computational load, we develop a relaxed version of the optimization problem. It is found that the relaxation enables us to efficiently solve the problem using water filling strategy. Computer simulations show that both GA and water filling solutions are superior to conventional solutions without CSI mismatch consideration, while the water filling is 1000 times faster than the GA.

## Introduction

Relaying technique is capable of extending communication range and coverage by providing link to shadowed users via relay nodes. From the perspective of signal processing, the cooperative relays can be viewed as a virtual antenna array which provides spatial diversity to combat frequency/time fading of channels. Attracted by its obvious merits, relaying technique has received extensive study in recent years [1–3].

The advantages of multiple-input multiple-output (MIMO) systems can be exploited in relay communications by accommodating multiple antennas at the nodes [4]. Recently, many works concentrate on designing non-regenerative MIMO relay systems. With perfect channel state information (CSI) assumption of all hops, optimal designs of MIMO amplify-and-forward (AF) relay networks were proposed [5–9]. To take CSI mismatch into account, robust design of a three-node MIMO relay system against CSI mismatch for linear non-regenerative MIMO relays was established in [10], where quality of service was attained by minimizing the averaged transmission power subject to mean square error (MSE) constraints at each data stream. Assuming the CSI uncertainty lies in a norm bounded region, two performance metrics, mutual information and MSE, were adopted to design the MIMO AF relay precoders in [11]. With the goal to minimize the MSE between the transmitted and the estimated symbol, Zhang proposed a robust precoder for MIMO AF relay systems against channel estimation error [12]. In [13], MSE minimization criterion was proposed to deal with CSI mismatch, and a closed-form solution was derived. In [14], two schemes aiming to maximize the signal-to-interference noise ratio were proposed to deal with CSI mismatch. In [15], the channel quantization error was considered and robust precoding schemes were proposed based on zero forcing and minimum MSE criterion. Also, a novel precoding scheme for MIMO relay systems in the presence of imperfect CSI was introduced in [16], where a base station precoding matrix and a relay station precoding matrix were created.

All the aforementioned algorithms focus on enhancing robustness against CSI estimation error. Besides CSI estimation error, other factors can also cause mismatch. In [17], quatized CSI feedback was considered and a strategy of scaling quantization quality of both two-hop links was proposed. Besides, channel feedback delay also has a significant influence on the performance of an AF relay system [18]. Therefore, it is highly desirable to consider both CSI estimation error and feedback delay when designing the MIMO AF relay system.

This paper presents a robust design of a MIMO AF relay system against CSI feedback delay and estimation error. The conditional MSE between the desired and received signals upon the estimated CSI is derived based on models of CSI feedback delay and estimation error. The conditional MSE is minimized subject to relay power constraint to optimize the beamforming matrix. Two solutions, one global and one relaxed, are proposed. Computer simulations show that both solutions outperform the MMSE strategy in terms of bit error rate (BER). It is also shown that ignoring the feedback delay results in higher BER. Although the relaxed solution is not as good as the global solution, it requires 1000 times less CPU time than the global solution. Hence, the relaxed solution provides a fast option for practical use.

## Relay system mode

Figure 1 depicts a relay network consisting of a source node, a relay node, and a destination node, where the source and destination nodes are equipped with *M* antennas, and the relay node is equipped with *L* antennas. It is assumed that *L*≥*M*. Direct link between the source node and the destination node is supposed to be absent. Frequency flat fading channels are considered. The source to relay and relay to destination channel matrices are represented by **H**
_{1} and **H**
_{2}, respectively. It is assumed that all elements of the two matrices are independent and identically distributed complex Gaussian variables with zero mean and variances \(\sigma _{h_{1}}^{2}\) and \(\sigma _{h_{2}}^{2}\).

The output of the system, denoted by \(\widehat {\mathbf {s}}\), is given by

where **s** is the transmitted signal, **W** denotes the equalization matrix at the destination node, **Q** denotes the relay beamforming matrix and **n**=**H**
_{2}
**Q**
**n**
_{1}+**n**
_{2} is the additive noise at the destination node. **n**
_{1} and **n**
_{2} are the additive complex Gaussian noises with zero mean and covariance \({\sigma _{1}^{2}}\mathbf {I}_{N_{r}}\) and \({\sigma _{2}^{2}}\mathbf {I}_{N_{s}}\), respectively.

In the conventional MMSE strategy, **Q** and **W** are optimized by minimizing \(E\left (\parallel \mathbf {s}-\mathbf {\widehat {s}}{\parallel _{2}^{2}}\right)_{\mathbf {s},\mathbf {n}}\) subject to relay transmission power constraint [19]. This criterion is optimal if the estimated channel matrices are the same as the real ones. However, due to estimation error and feedback delay, there is usually certain mismatch between the estimated and the real channel matrices in practice. Therefore, a more realistic algorithm should take the CSI mismatch into account when designing **Q** and **W**.

## Models of CSI feedback delay and estimation error

### Source to relay CSI

Denote the (*i*,*j*)^{th} element of **H**
_{1} as \(h_{1}^{(i,j)}\). The maximum likelihood (ML) estimate of \(h_{1}^{(i,j)}\) is represented by \(\widehat {h}_{1}^{i,j}\). According to the stochastic error model [20], \(h_{1}^{(i,j)}\) and \(\widehat {h}_{1}^{i,j}\) are related by

where \(\epsilon _{1}^{(i,j)}\) denotes the estimation error which is modeled by a complex Gaussian variable with zero mean and variance \(\sigma _{\epsilon _{1}}^{2}={{\sigma _{1}^{2}}}/{(N_{1}-1)}\). *N*
_{1} is the number of training samples to obtain \(\widehat {h}_{1}^{(i,j)}\), and \({\sigma _{1}^{2}}\) is the noise level at relay node.

Based on the Bayesian theory, the conditional expectation of \(h_{1}^{(i,j)}\) upon \(\widehat {h}_{1}^{i,j}\) is given by

Conditional correlation of \(h_{1}^{(i,j)}\) and \(h_{1}^{(k,l)}\) upon \(\widehat {h}_{1}^{(i,j)}\) and \(\widehat {h}_{1}^{(k,l)}\) is given by

where \(\sigma ^{2}_{h_{1}|\widehat {h}_{1}}=\sigma _{h_{1}}^{2}-{\sigma _{h_{1}}^{4}}{/\left (\sigma _{h_{1}}^{2}+\sigma _{\epsilon _{1}}^{2}\right)}\).

### Relay to destination CSI

Since the estimated CSI needs to be fed back from the destination to relay, due to channel temporal variation, the estimated CSI is actually an early version of the real CSI. Denote the (*i*,*j*)^{th} element of **H**
_{2} as \(h_{2}^{(i,j)}\). According to the channel temporal variation model, \(h_{2}^{\tau (i,j)}\) which represents an early version of \(h_{2}^{(i,j)}\) is related to \(h_{2}^{(i,j)}\) via [21]

where *ρ* is the normalized correlation coefficient calculated as *ρ*=*J*
_{0}(2*π*
*f*
_{
d
}
*τ*). *J*
_{0} denotes the first kind Bessel function of order zero, and *f*
_{
d
} is the maximum Doppler frequency. *ζ*
^{(i,j)} in (5) is a complex Gaussian variable with zero mean and variance \(\left (1-\rho ^{2}\right)\sigma _{h_{2}}^{2}\).

Denote the estimate of \(h_{2}^{\tau (i,j)}\) by \(\widehat {h}_{2}^{\tau (i,j)}\). Based on the stochastic error model, the relationship between \(\widehat {h}_{2}^{\tau (i,j)}\) and \(h_{2}^{\tau (i,j)}\) is given by

where \(\epsilon _{2}^{(i,j)}\) is a complex Gaussian variable with zero mean and variance \(\sigma ^{2}_{\epsilon _{2}}={{\sigma _{2}^{2}}}/{(N_{2}-1)}\). *N*
_{2} and and \({\sigma _{2}^{2}}\) are the number of training samples and the noise level in estimating \(\widehat {h}_{2}^{\tau (i,j)}\), respectively.

With (5) and (6), using the Bayesian theory, the conditional expectation of \(h_{2}^{(i,j)}\) upon \(\widehat {h}_{2}^{\tau (i,j)}\) is given by [21]

The conditional correlation of \(h_{2}^{(i,j)}\) and \(h_{2}^{(k,l)}\) upon \(\widehat {h}_{2}^{(i,j)}\) and \(\widehat {h}_{2}^{(k,l)}\) is given by

where \(\sigma ^{2}_{h_{2}|\widehat {h}_{2}}=\sigma _{h_{2}}^{2}-{\rho ^{2}\sigma _{h_{2}}^{4}}/{\left (\sigma _{h_{2}}^{2}+\sigma _{\epsilon _{2}}^{2}\right)}\).

## Robust design of a MIMO AF relay system

### Minimum conditional MSE criterion

Taking CSI mismatch of **H**
_{1} and **H**
_{2} into account, the proposed objective function is to minimize the conditional MSE between the transmitted and received signals subject to relay power constraint, which is given as follows:

where *P*
_{
r
} is the upper bound of relay transmission power.

Assuming that **s** and **n** are independent, based on (1), setting the gradient of (9a) with respect to **W**
^{H} to zero yields the optimal **W**

where **W**
_{opt} denotes the optimal solution of **W**, and \(\mathbf {R}_{n}={\sigma _{s}^{2}}E\left [\mathbf {HH}^{H}|\mathbf {\widehat {H}}_{1},\mathbf {\widehat {H}}_{2}\right ] +{\sigma _{1}^{2}}E\left [\mathbf {H}_{2}\mathbf {QQ}^{H}\mathbf {H}_{2}^{H}|\mathbf {\widehat {H}}_{2}\right ] +{\sigma _{2}^{2}}\mathbf {I}_{N_{s}}\).

Define \(\mathbf {\overline {H}}=E\left [\mathbf {H}|\mathbf {\widehat {H}}_{1},\mathbf {\widehat {H}}_{2}\right ]\), \(\mathbf {\overline {H}}_{1}=E\left [\mathbf {H}_{1}|\mathbf {\widehat {H}}_{1}\right ]\) and \(\mathbf {\overline {H}}_{2}=E\left [\mathbf {H}_{2}|\mathbf {\widehat {H}}_{2}\right ]\). Substituting (1) and (10) into (9a) yields (11) (shown on the top of next page).

It is observed from (11a) that **Q** is contained in the inversion manipulation; therefore, direct optimization of (11) is difficult. To facilitate the solution of (11), the structure of optimal **Q** is analyzed and it is found that optimal **Q** has the form of

where **Q**
_{opt} is the optimal **Q**, *Φ*_{1} is an *M*×*M* diagonal matrix, and \(\mathbf {V}_{2}^{H}\) and **U**
_{1} are unitary matrices constituted by right- and left-singular vectors of \(\overline {\mathbf {H}}_{2}\) and \(\overline {\mathbf {H}}_{1}\), respectively. The proof of (12) is provided in Appendix A. Using (12), (1) is equivalently expressed as

where

*ϕ*
_{
i
} is the *i*
^{th} diagonal element of *Φ*_{1}, \(\alpha _{i}=\lambda _{1,i}^{2}\lambda _{2,i}^{2}\), \(\beta _{i}={\sigma _{s}^{2}}\left (\lambda _{1,i}^{2}+\sigma _{h_{1}|\widehat {h}_{1}}^{2}\right)\lambda _{2,i}^{2}+{\sigma _{1}^{2}}\lambda _{2,i}^{2}\), \(\gamma _{i}={\sigma _{s}^{2}}\mid \lambda _{1,i}\mid ^{2}+\)
\({\sigma _{s}^{2}}\sigma _{h_{1}|\widehat {h}_{1}}^{2}+{\sigma _{1}^{2}},\)

and *λ*
_{1,i
} and *λ*
_{2,i
} denote the *i*
^{th} singular values of \(\overline {\mathbf {H}}_{1}\) and \(\overline {\mathbf {H}}_{2}\), respectively. Details of transformation from (11) to (13) are presented in Appendix B. The optimization problem (13) is easier to solve than (11). Once the optimal *Φ*_{1} is derived, **Q**
_{opt} and **W**
_{opt} are computed using (12) and (10), respectively.

When \(\sigma _{h_{1}|\widehat {h}_{1}}^{2}=0\) and \(\sigma _{h_{2}|\widehat {h}_{2}}^{2}=0\), *b*=0, *c*=0, \(\beta _{i}={\sigma _{s}^{2}}\lambda _{1,i}^{2}\lambda _{2,i}^{2}+{\sigma _{1}^{2}}\lambda _{2,i}^{2}\), and \(\gamma _{i}={\sigma _{s}^{2}}\mid \lambda _{1,i}\mid ^{2}+{\sigma _{1}^{2}},\). (13) becomes

which are equivalent to (24) and (25) of [19]. Therefore, the minimum conditional MSE criterion reduces to the MSE criterion when the CSI mismatch vanishes.

### Global solution by genetic algorithm

It is observed from (13) that multiplying *ϕ*
_{
i
} by a constant maximizes the relay power and the value of the objective function (13a). Therefore, the optimal *ϕ*
_{
i
} can be obtained while (13b) achieves equality. Based on this observation, the values of chromosomes in genetic algorithm (GA) is optimized within 0 and 1, then multiplied with a constant *α* which is obtained when (13b) achieves equality, i.e.,

Substituting *α*
*ϕ*
_{
i
} into (13a) yields the value of the fitness function of GA.

### Relaxed solution by water filling strategy

It is observed from (13a) that the terms of *b* and *c* contain *ϕ*
_{
i
},*i*=1,…,*M*, which prevent deriving an analytical solution to (13). To avoid high computational loads of using the global searching algorithm, a relaxed version of (13) is proposed here.

From (13a), it is noted that increasing the values of *b* and *c* reduces the value of *J*(*Φ*_{1}), which means

where *J*(*Φ*_{1})_{max} is computed from (13a) using *b*
_{max} and *c*
_{max}. Here, *b*
_{max} and *c*
_{max} denote the maximum values of *b* and *c*, respectively. From the relay power constraint (13b), the possible values of *b*
_{max} and *c*
_{max} are straightforward to derive and are given below:

Substituting (20) and (21) into (13), and using the Lagrange multiplier technique, the solution of the relaxed version of (13) is given by

where

and (*x*)^{+}=max(*x*,0), *μ* is the Lagrange constant which should be chosen such that (13b) is satisfied.

## Computer simulations

To demonstrate validity and advantages of the proposed strategy, the following simulation scenarios are devised. The channel coefficients are assumed to be complex Gaussian variables with zero mean and variance \(\sigma _{h_{1}}^{2}\) and \(\sigma _{h_{2}}^{2}\). The signal-to-noise ratio (SNR) of the backward and forward channels are defined as \({\mathrm {SNR_{1}}}={\sigma _{s}^{2}}\sigma _{h_{1}}^{2}/{\sigma _{1}^{2}}\) and \({\mathrm {SNR_{2}}}=P_{r}\sigma _{h_{2}}^{2}/{\sigma _{2}^{2}}L\), respectively. The source is generated from a QPSK constellation. It is assumed that the number of data samples is 120, containing *N*
_{
tr
}=20 training samples. The number of *M* is 4. The relay transmission power *P*
_{
r
} is 0 dB. The BER is derived from 100 independent trials for all the plots. The proposed algorithm is compared with the conventional MMSE strategy and the robust method which only considers CSI estimation error.

Figure 2 plots the BER versus SN*R*
_{1} with *L*=4, SN*R*
_{2}=10 dB, and *ρ*=0.90. Because CSI estimation error exists, it is observed that the proposed methods outperform the conventional water filling method which assumes accurate estimation. Also, it is found that both of the proposed methods perform better than the robust method against estimation error, since the channel correlation coefficient is 0.90 in this situation, which has significant influence on the accuracy of CSI. Furthermore, because GA is global searching, the proposed algorithm using GA yields smaller BER than that of the proposed water filling method.

In Fig. 3, we assume a better communication environment. The SN*R*
_{1} is fixed at 20 dB, and the value of *ρ* is assumed as 0.95. Figure 3 shows the BER versus the value of SN*R*
_{2}. It can be seen that the performance of all the methods becomes better as the value of SN*R*
_{2} increases. Compared with Fig. 2, all the methods give smaller BER when SN*R*
_{1}=20 dB and SN*R*
_{2}=10 dB, since the value of *ρ* increases. As SN*R*
_{2} increases beyond 10 dB, the performance of the conventional water filling technique and the robust method against CSI estimation error do not have significant improvement, while BER of the proposed methods continue to decrease.

Figure 4 plots BER versus the value of *L* with SN*R*
_{1} and SN*R*
_{2} equal to 10 dB. It is natural to see that the BER of all the methods decrease as the number of antennas equipped on relay increases. Furthermore, it is observed that the increase of *L* favors the improvement of performance of the conventional water filling technique and the robust algorithm against CSI estimation error. The gap between these two methods and the proposed methods becomes smaller.

In Fig. 5, plots of BER versus the value of *ρ* are given. BERs of all the methods increase as the value of *ρ* decreases. It is verified that channel correlation has significant compact on the performance of all the methods. The proposed methods outperform the conventional MMSE method for all values of *ρ*. When *ρ*=1, performance of the proposed methods and the robust algorithm against CSI estimation error become same.

Figure 6 shows the cumulative distribution function (CDF) of BER with SN*R*
_{1}=SN*R*
_{2}=10 dB, *ρ*=0.90, and *L*=4. It is observed that the proposed methods yield smaller BER at a high probability than that of the other methods.

To compare computational time of the GA and the relaxed solution, 100 independent trials are conducted. The number of antennas at all the nodes are assumed to be the same. From Table 1, it is observed that the relaxed solution cost far less computational time than that of the GA.

## Conclusions

In this paper, we consider the robust design of an AF relay system against CSI mismatch. A relay system equipped with multiple antennas is considered. The conditional expectation of the mean square error is minimized with respect to precoding matrix and equalization coefficient. Computer simulations show that when the estimated CSI is different from the real CSI, the proposed strategy outperforms the MMSE strategies with perfect CSI assumption. The advantages include lower BER with different values of *ρ*, SN*R*
_{1}, SN*R*
_{2}, and *L*. Also, the proposed strategy is more likely to give small BER compared to the conventional MMSE strategy. Therefore, the proposed strategy is more reliable to be used in real applications.

## Appendix A

###
*Proof of the Theorem*

Assume \(\mathbf {Q}=\mathbf {V}_{2}\boldsymbol {\Phi }\mathbf {U}_{1}^{H}\), where \(\boldsymbol {\Phi }=\left (\begin {array}{cc} \boldsymbol {\Phi }_{1} & \boldsymbol {\Phi }_{2} \\ \boldsymbol {\Phi }_{3} & \boldsymbol {\Phi }_{4} \end {array} \right)\) and \(\boldsymbol {\Phi }_{1}\in \mathcal {C}^{N_{s}\times N_{s}}\).

From (4), we have

Substituting ** Φ** and (24) into \(E\left [\mathbf {HH}^{H}|\mathbf {\widehat {H}}_{1},\mathbf {\widehat {H}}_{2}\right ]\) yields (25) shown as follows,

where \(\boldsymbol {\Phi }_{i}^{T}\) denotes the *i*
^{th} row of ** Φ**.

Similarly, we have

Define

where \(\mathbf {z}\in \mathcal {C}^{N_{s}\times 1}\).

From (26) and (27), it is straightforward to show that

Suppose **A**, **B**, and ** △** are non-negative matrices, then

Therefore, (28) and (29) gives the result that

reaches its maximum when ** Φ**=

*Φ*_{0}.

Also, it is easy to show that ** Φ**=

*Φ*_{0}helps to save relay transmission power. Therefore, (11) has a solution as \(\mathbf {Q}=\mathbf {V}_{2}\boldsymbol {\Phi }_{0}\mathbf {U}_{1}^{H}\). □

## Appendix B

###
*Proof of Lemma 1*

From the theorem, we have

where \(b=\sigma _{h_{2}|\widehat {h}_{2}}^{2}\sum _{i=1}^{N_{s}}\left (\lambda _{1,i}^{2}+\sigma _{h_{1}|\widehat {h}_{1}}^{2}\right)\mid \phi _{i}\mid ^{2}\), and **u**
_{2,i
} denotes the *i*
^{th} column of **U**
_{2}.

Similarly, it can be derived that

where \(c=\sigma _{h_{2}|\widehat {h}_{2}}^{2}\sum _{i=1}^{N_{s}}\mid \phi _{i}\mid ^{2}\).

Furthermore, we have

## References

- 1
JN Laneman, DNC Tse, GW Wornell, Cooperative diversity in wireless networks: efficient protocols and outage behavior. IEEE Trans. Inf. Theory.

**50**(12), 3062–3080 (2004). - 2
CS Patel, GL Stuber, TG Pratt, Statistical properties of amplify and forward relay fading channels. IEEE Trans. Veh. Technol.

**55**(1), 1–9 (2006). - 3
G Zheng, Collaborative-relay beamforming with perfect CSI: optimum and distributed implementation. IEEE Signal Proc. Lett.

**16**(4), 257–260 (2009). - 4
B Wang, J Zhang, A Host-Madsen, On the capacity of MIMO relay channels. IEEE Trans. Inf. Theory.

**51**(1), 29–43 (2005). - 5
C-C Chien, H-J Su, H-J Li, Joint beamforming and power allocation for MIMO relay broadcast channel with individual SINR constraints. IEEE Trans. Veh. Technol.

**63**(4), 1660–1677 (2014). - 6
V Havary-Nassab, S Shahbazpanahi, A Grami, Z-Q Luo, Distributed beamforming for relay networks based on second-order statistics of the channel information. IEEE Trans. Signal Process.

**56**(9), 4306–4316 (2008). - 7
O Munoz, J Vidal, A Augustin, Linear transceiver design in non-regenerative relays with channel state information. IEEE Trans. Signal Process.

**55**(6), 2593–2604 (2007). - 8
X Tang, Y Hua, Optimal design of non-regenerative MIMO wireless relays. IEEE Trans. Wirel. Commun.

**6**(4), 1398–1407 (2007). - 9
Y Rong, X Tang, Y Hua, A unified framework for optimizing linear non-regenerative multicarrier MIMO relay communication systems. IEEE Trans. Signal Process.

**57**(12), 4837–4851 (2009). - 10
Y Rong, in

*2010 Australian Communications Theory Workshop*. Robust design for linear non-regenerative MIMO relays (Canberra, 2010), pp. 87–92. - 11
H Shen, J Wang, Robust optimization for amplify-and-forward MIMO relaying from a worst-case perspective. IEEE Trans. Signal Process.

**61**(21), 5458–5471 (2013). - 12
Z Zhang, J Xu, L Qiu, Robust ARQ precoder optimization for AF MIMO relay systems with channel estimation errors. IEEE Trans. Wirel. Commun.

**12**(10), 5236–5247 (2013). - 13
C Xing, S Ma, YC Wu, Robust joint design of linear relay precoder and destination equalizer for dual-hop amplify-and-forward MIMO relay systems. IEEE Trans. Signal Process.

**58**(4), 2273–2283 (2010). - 14
Z Wang, W Chen, J Li, Efficient beamforming for MIMO relaying broadcast channel with imperfect channel estimation. IEEE Trans. Signal Process.

**61**(1), 419–426 (2012). - 15
B Zhang, Z He, K Niu, L Zhang, Robust linear beamforming for MIMO relay broadcast channel with limited feedback. IEEE Signal Proc. Lett.

**17**(2), 209–212 (2010). - 16
Y Cai, RC de Lamare, LL Yang, M Zhao, Robust MMSE precoding based on switched relaying and side information for multiuser MIMO relay systems. IEEE Trans. Veh. Technol.

**64**(12), 5677–5687 (2015). - 17
W Xu, X Dong, WS Lu, MIMO relaying broadcast channels with linear precoding and quantized channel state information feedback. IEEE Trans. Signal Process.

**58**(10), 5233–5245 (2010). - 18
HA Suraweera, TA Tsiftsis, Effect of feedback delay on amplify-and-forward relay networks with beamforming. IEEE Trans. Veh. Technol.

**60**(3), 1265–1271 (2011). - 19
W Guan, H Luo, Joint MMSE transciever design in non-regenerative MIMO relay systems. IEEE Commun. Lett.

**12**(7), 517–519 (2008). - 20
R Mo, YH Chew, C Yuen, Information rate and relay precoder design for amplify-and-forward MIMO relay networks with imperfect channel state information. IEEE Trans. Veh. Technol.

**61**(9), 3958–3968 (2012). - 21
AP Liavas, Tomlinson-Harashima precoding with partial channel knowledge. IEEE Trans. Commun.

**53**(1), 5–9 (2005).

## Acknowledgements

This work was supported by grants from the Central University Research Founding (ZYGX2014J018), National Science Foundation of China (61671121), and UESTC Research Start-up funding (ZYGX2016KYQD106).

### Competing interests

The authors declare that they have no competing interests.

## Author information

### Affiliations

### Corresponding author

## Rights and permissions

**Open Access** This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

## About this article

### Cite this article

Zhang, Y., Zhao, D., Zhao, H. *et al.* Optimization of a MIMO amplify-and- forward relay system with channel state information estimation error and feedback delay.
*EURASIP J. Adv. Signal Process. * **2016, **117 (2016). https://doi.org/10.1186/s13634-016-0400-9

Received:

Accepted:

Published:

### Keywords

- Channel state information mismatch
- Conditional expectation
- Multiple-input multiple-output
- Minimum mean square error
- Relay network