Research

# Training design for precoded BICM-MIMO systems in block-fading channels

Zohreh Andalibi12*, Ha H Nguyen12 and Joseph E Salt12

Author Affiliations

2 Department of Electrical and Computer Engineering, University of Saskatchewan 57 Campus Drive, Saskatoon, SK S7N5A9, Canada

For all author emails, please log on.

EURASIP Journal on Wireless Communications and Networking 2012, 2012:80  doi:10.1186/1687-1499-2012-80

 Received: 7 August 2011 Accepted: 4 March 2012 Published: 4 March 2012

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

### Abstract

In order to improve bandwidth efficiency and error performance, a new training scheme is proposed for bit-interleaved-coded modulation in multiple-input multiple-output (BICM-MIMO) systems. Typically, in a block-fading channel, the training overhead used for obtaining channel knowledge is proportional to a power of 2 of the number of transmit antennas. However, this overhead can be reduced by embedding pilot symbols within data symbols before precoding. The values, positions, and the number of pilot symbols are found by minimizing the Cramer-Rao bound on the channel estimation error. Computer simulations are presented to demonstrate the advantage of the proposed scheme over other training methods, in terms of both the mean-square-error of the channel estimation and the system's frame-error-rate.

##### Keywords:
BICM-MIMO; block fading; channel estimation; training design; pilot symbols; Cramer-Rao bound; iterative receiver

### 1 Introduction

The pioneering work on multiple-input multiple-output (MIMO) systems [1] shows that a MIMO system can provide a multiplexing gain and accordingly high spectral efficiency over slow fading channels. On the other hand, to achieve a high diversity order, space-time transmission techniques can be implemented at the transmitter [2,3]. To achieve both high diversity order and coding gain in coded modulation systems, the concept of space-time transmission has also been applied [4,5]. In such systems, space-time transmission is typically implemented using a linear space-time matrix, or equivalently a linear precoder, so that a single modulation symbol is efficiently transmitted across multiple transmit antennas. Among many research works on precoder design for coded modulation systems with multiple antennas, the design that considers all the relevant components of the transmitter, namely precoding, modulation, and interleaver, can be found in [5-7]. Specifically, a full-rate precoder with any size and for any number of transmit antennas is designed in [6] to maximize the achievable diversity order and coding gain in MIMO block-fading channels.

It is shown in [6] that the maximum achievable diversity order can be realized by an iterative receiver that employs a soft-input soft-output detector [5] and under the assumption of having the perfect channel state information (CSI) at the receiver. In practice, however, CSI has to be estimated using a channel estimator and it is never perfect. Two types of channel estimators have been used for MIMO block-fading channels in coded modulation systems, i.e., training-based and semi-blind channel estimators [8,9]. In both types of channel estimators, known signals are used to estimate the CSI at the first iteration of the iterative receiver.

Conventionally, for block-fading channels, known signals or the training sequence is included at the beginning of each data block, which is called time-multiplexed training or pilot symbol-assisted modulation (PSAM) scheme [10]. This scheme however reduces bandwidth efficiency of MIMO systems, since the amount of training overhead needed is at least a power of 2 of the number of transmit antennas [11] to ensure the identifiability of the MIMO channel. A straightforward application of the PSAM scheme to a BICM-MIMO system would be time-multiplex data information with the training information after the precoder.

As an alternative to the above conventional PSAM scheme, a potential benefit can be sought by time-multiplexing data information with the training information before the precoder in the transmitter. This new approach shall reduce the required training overhead compared to the conventional PSAM, since the transmitted training symbols are spread over more time periods; thanks to the precoder. This approach shall be referred to as precoded PSAM (PPSAM). Investigating power and time allocations of the training symbols in PPSAM scheme is the main objective of this article.

Moreover, by multiplexing the training sequence before precoder, training symbols can be exploited in both the initialization and iteration phases of the iterative channel estimation process. This is different from a conventional iterative channel estimator using PSAM scheme, in which training sequence is only used at the initialization phase. A natural question is whether the optimal training design for the initialization phase using PPSAM scheme is still optimal for subsequent iterations of an iterative channel estimator. On the one hand, the channel estimation error at the initialization phase translates to an SNR shift in the BER performance [8]. On the other hand, the channel estimation error from the last iteration of the iterative estimator has a strong impact on the error floor of the BER performance [12]. Therefore, optimal training sequence should be designed carefully that considers both initialization and iteration phases.

One of different criteria that have been used to design training sequences is the minimization of the Cramer-Rao bound (CRB) of the channel estimation error [10]. This criterion shall be used in this article due to two main reasons. First, it is directly related to the channel estimation error. Second, since the CRB is a lower bound on the mean-squared-error (MSE) of any unbiased estimator, designing training sequences using this criterion would be applicable to many estimation algorithms. Other design criteria, such as maximizing the channel capacity [8] and minimizing the outage probability [13], are based on some specific channel estimation algorithms.

The article is organized as follows. The system model of BICM-MIMO is presented in Section 2. In Section 3 a lower bound on the MSE of the channel estimator is obtained and the training sequence is designed by minimizing this bound. Section 4 provides numerical results and comparisons. Section 5 concludes the article.

### 2 System model

Figure 1 shows the block diagram of a BICM-MIMO system under consideration. At the transmitter, a channel encoder with a rate-r error-correcting code converts the vector of information bits b into a codeword c. The coded bits are then interleaved by a random interleaver as described in [6] to produce the interleaved codeword c ̃ . The interleaved codewords are segmented into groups of (Nnt - Np) × m bits, where N is the spreading factor of the precoder, nt is the number of transmit antennas, Np is the number of pilot symbols in Nnt precoded symbols and m is the number of bits carried by one symbol of a QAM constellation whose size is |Ω| = 2m. Next, the coded bits are mapped to (Nnt - Np) QAM constellation points. In this step, Np known pilot symbols are inserted in every segmented group of (Nnt - Np) data symbols to produce N super-symbols. Here, each super-symbol refers to a group of nt consecutive symbols. Investigating the positions and the number of pilot symbols (i.e., Np) to be used in each Nnt symbols is the main objective of this article.

Figure 1. Block diagram of a BICM-MIMO system with a linear precoder and proposed training insertion.

Every group of N super-symbols is then spread over N time periods using a linear precoder G. The Nnt × Nnt matrix G multiplies a vector of Nnt QAM symbols at the precoder input, and generates Nnt symbols to be transmitted over nt antennas, over N time periods.

This is illustrated in Figure 2. Let x k = [ x ( k - 1 ) N n t + 1 , x ( k - 1 ) N n t + 2 , , x ( k - 1 ) N n t + N n t ] be the kth vector to be precoded. Then, xkG gives the precoded symbols. Here, xi's are complex data or pilot symbols belonging to the 2m-QAM constellation Ω. It is assumed that the data symbols xi's are i.i.d with variance σ x 2 . After precoding, precoded symbols are transmitted through nt transmit antennas over a block-fading channel.

Figure 2. Spreading a precoded symbol over nt antennas and N time periods-denoted by Z k = [ Z 1 , k , Z 2 , k , , Z n t , k ] T .

With nt transmit antennas and nr receive antennas, the channel is modeled by an nt × nr matrix. For frequency-flat Rayleigh fading, coefficients of the channel matrix are i.i.d. zero-mean circularly symmetric complex Gaussian random variables with variance σ h 2 . The channel is assumed to be block fading with nc different channel realizations during each codeword. For the kth symbol to be precoded, xk, the Nnt × Nnr extended channel matrix, Hk, can be written as

H k = diag H k [ 1 ] , , H k [ 1 ] , N / n s H k [ 2 ] , , H k [ 2 ] , , H k [ n s ] , , H k [ n s ] , (1)

where ns is the number of distinct channel realizations during N time periods of each codeword. To simplify the notation it is also assumed(a) that ns divides N. For example, if the length of a codeword is 64 and nc = 32, then choosing N = 2 would make ns = 1, whereas choosing N = 4 gives ns = 2. Notation H k [ t ] refers to the nt × nr complex matrix k that defines the tth channel realization included in ns channel realizations. The extended channel input/output relationship is expressed by

y k = x k G H k + w k (2)

where y k = [ y ( k - 1 ) N n r + 1 , y ( k - 1 ) N n r + 2 , , y ( k - 1 ) N n r + N n r ] is the received vector at the kth precoding time period and wk is the noise vector with size 1 × Nnr whose components are i.i.d zero-mean circularly symmetric Gaussian random variables with variance N0. It is noted from (2) that although both data and pilot symbols are precoded, the part of the precoder that multiplies the pilot symbols depends on the positions of the pilot symbols in xk. Equivalently, the design of the pilot symbols is governed by the properties of the precoder used. Since this study adopts the transmission framework and precoder design in [6], it is useful to review the properties of the precoder proposed in [6].

In general, the properties of the precoder in [6] are established by the maximum-likelihood decoding analysis and an assumption of ideal channel interleaving. Specifically, this linear precoder which achieves full diversity order and maximum coding gain satisfies the following two conditions:

• A genie condition, which guarantees orthogonal and equal norm sub-rows in the linear precoding matrix. Each sub-row has size nt in a precoding matrix with size Nnt × Nnt.

• Dispersive nucleo algebraic (DNA) condition, which is based on Proposition 2 in [6], forces null and orthogonal nucleotides with size s' = N/ns. Nucleotides refer to subparts of sub-rows with size s'.

A linear precoder that satisfies the above two sets of conditions is called DNA-cyclo precoder and has the best performance in terms of achieving diversity and coding gains with low complexity receiver when N nt. It is suggested in [6] that to generate one class of such a precoder, a Ns' × Ns' cyclotomic rotator, denoted by Φ, that satisfies the genie condition is first selected. Then the orthogonal nucleotides are placed inside an Nnt × Nnt matrix and they are separated with null nucleotides. Therefore, the DNA-cyclo precoder matrix can be expressed by subparts of a cyclotomic rotator as follows:

G = I n t / s Φ [ 1 ] [ 1 ] I n t / s Φ [ N ] [ 1 ] I n t / s Φ [ 1 ] [ 2 ] I n t / s Φ [ N ] [ 2 ] I n t / s Φ [ 1 ] [ N s ] I n t / s Φ [ N ] [ N s ] (3)

where Φ[i ] [j] is the ith sub-row of the jth row of Φ with size 1 × s', In is an identity matrix with size n × n and ⊗ denotes the Kronecker product.

The properties that shall be useful for the problem considered in this article, which are implied directly from the genie and DNA conditions, are ΦΦH = INs ' and Φ [ i ] [ t ] ( Φ [ j ] [ t ] ) H = 1 N δ ( i - j ) . It is also useful to point out that each component of Φ has an exponential form with a scaling factor of 1 N s .

The iterative receiver is also shown in Figure 1. The channel estimator produces an estimate of the channel using the minimum MSE (MMSE) criterion based on the training sequence. Details about channel estimation with the proposed method of inserting training sequence shall be given in Section 3. After channel estimation is performed using the training signal, the soft-input soft-output demodulator uses the MMSE criterion to demodulate the data. The soft-output MMSE demodulator computes the extrinsic information for the interleaved bits, { Λ ext ( c ̃ l ) } , from the received symbols. To obtain Λ-values, the demodulator exploits the a priori information of the coded bits coming from the decoder, { Λ ap ( c ̃ l ) } , and the channel estimate H ^ k . In the first iteration, the demodulator assumes that the a priori Λ-values are zero, except for the pilot symbols. For the corresponding bits of the pilot symbols, the demodulator uses a large number, say ±100 as their a priori Λ-values. The de-interleaved outputs, i.e., { Λ ap ( c l ) } , become the a priori Λ-values used in the channel decoder shown in Figure 1 after removing the information of pilot symbols. The channel decoder uses the maximum a posteriori probability (MAP) algorithm to compute the extrinsic Λ-values { Λ ext ( c l ) } . for all coded bits, which are used again in the next iteration in the demodulator. In subsequent iterations, soft information from the decoder is used to improve the performance of the channel estimator. The detailed operation of the iterative channel estimator is discussed in the following sections.

### 3 Training design and channel estimator

As discussed before, the criterion used for training design in this article is the CRB on the channel estimation error. The bound states that the MSE of any unbiased estimator is lower bounded by the trace of inverse of complex Fisher information matrix (FIM) [14]. To derive FIM, the relation between the channel input and channel output during one block-length, i.e., N/ns time periods, whose corresponding channel matrix is H k [ t ] , is of interest. In the following, index k is omitted, since it suffices to consider the transmission of a single precoded symbol for the purpose of channel estimation. With the previously described structure of the precoder, the channel output during one super-symbol time is given by

y [ i , t ] = I n r τ = 1 N s x [ τ ] Φ [ i , t ] [ τ ] h [ t ] + w [ i , t ] ; t = 1 , , n s , i = 1 , , s (4)

where y[i,t] = y[(t-1)s'+i] represents the ((t - 1)s' + i)th received symbol during N time periods, with size nr × 1. Moreover, h[t] is the column vector formed by vertically stacking the columns of an nt × nr channel realization matrix H[t] and x[τ]'s are constructed by splitting x in Ns' sub-vectors with size 1 × nt/s'. In the following, we call these sub-vectors x[τ]'s nucleo symbols.

It is quite obvious from (4) that, to have all the received super-symbols, y[i,t], contain training information, there should be at least one pilot nucleo (i.e., nt/s' pilot symbols) in each group of Ns' nucleos to be precoded.

With the above structure of the proposed training sequence, the number of pilot symbols in Nnt transmitted symbols would be Np = np × nt/s', where np nucleo symbols in a symbol to be precoded are assigned to training sequence. Therefore, (4) can be rewritten as

y [ i , t ] = I n r τ I d x d [ τ ] Φ d [ i , t ] [ τ ] + τ I p x p [ τ ] Φ p [ i , t ] [ τ ] h [ t ] + w [ i , t ] , (5)

where I d and I p are sets of indexes from {1, . . . , Ns'}, that are assigned to data and pilot nucleos, respectively, and | I d | + | I p | = ( N s - n p ) + n p = N s . Note that the subscripts "d" and "p" are used to differentiate between data and pilot nucleos. For convenience, the notations Φ p [ i , t ] [ τ ] and Φ d [ i , t ] [ τ ] are used to refer to sub-rows of Φ that are multiplied by pilot and data nucloes, i.e., x p [ τ ] and x d [ τ ] , respectively. Furthermore, in the following the notation T[i,t] is used for I n r τ I p x p [ τ ] Φ p [ i , t ] [ τ ] .

The derivation of FIM is given in the next section. Pilot symbols are exploited at the initialization phase and in subsequent iterations considering the special structure of the training sequence. In general, training design can be investigated for these two phases separately. However, for the precoder adopted in this article, the optimal training design obtained for the initialization phase turns out to also be optimal for the iteration phase. Nevertheless, the optimal numbers of pilot nucleos in these two phases of channel estimation are not the same.

#### 3.1 Fisher information matrix

The key steps in deriving the FIM in the initialization phase are now given. Without loss of generality we drop superscript t in (5) and perform all the derivations for the first block period (i.e., t = 1). Collecting all the observations during the first block period of length s' in a vector φ, the FIM for the channel estimation problem at the initialization phase is defined and computed as

FI M init ( n p , x p , I p ) = E φ , h ln p ( φ , h ) h *   ln p ( φ , h ) h * H = E h E φ ln p ( φ | h ) h *   ln p ( φ | h ) h * H h + E h ln p ( h ) h *   ln p ( h ) h * H (6)

where F I M i n i t ( n p , x p , p ) shows the dependency of FIM on those parameters of interest. Using the i.i.d. assumption on noise and data, p(φ|h) can be approximated as a complex normal distribution with mean μ = [ μ 1 T , , μ s T ] T and covariance Rφ = diag[R1, . . . , Rs ']. Moreover, it follows from (5) that μi = Eφ {y[i]|h} = T[i]h and

R i = H σ x 2 I n t / s ( ( Φ d [ i ] ) T ( Φ d [ i ] ) * ) H H + N 0 I n r (7)

where H = ( H [ 1 ] ) T and Φ d [ i ] is the ith sub-matrix of Φ with size (N s' - np) × s' that is assigned to data symbols.

The i.i.d. assumptions on noise and data make the FIM additive. Specifically, FI M init ( n p , x p , I p ) = i = 1 s FIM i init . The quantity FIM i init is obtained as follows:

FIM i init = E h E y ln p ( y | h ) h * ln p ( y | h ) h * H h + σ h - 2 I n t n r .

We know that

ln p ( y | h ) = Constant - ln | R i | - ( y - μ i ) H R i - 1 ( y - μ i ) . (8)

and ln | R i | h l * = trace R i - 1 R i h l * . Therefore,

R i h l * = H σ x 2 I n t / s ( ( Φ d [ i ] ) T ( Φ d [ i ] ) * ) Σ l T (9)

where ∑l is an nr × nt null matrix with only a single element of 1 at position l - 1 n t + 1 , l - 1 mod n t + 1 . The derivative of the third term in (8) is

( y - μ i ) H R i - 1 ( y - μ i ) h l * = - μ i H h l * R i - 1 ( y - μ i ) + ( y - μ i ) H R i - 1 h l * ( y - μ i )

where R i - 1 h l * = - R i - 1 R i h l * R i - 1 and R i h l * is given by (9). In addition,

μ i H h l * = h H h l * ( T [ i ] ) H = e l T ( T [ i ] ) H

where el is an ntnr × 1 null vector with a single element 1 at position l.

Using all the above equations and after some manipulations, one has

( FIM i init ) l , j = E h { e l T ( T [ i ] ) H R i - 1 T [ i ] e j + tr ( R i - 1 H A [ i ] Σ l T R i - 1 Σ j ( A [ i ] ) H H H ) } + σ h - 2 δ ( l - j ) ,

where A [ i ] σ x 2 I n t / s ( ( Φ d [ i ] ) T ( Φ d [ i ] ) * ) .

Using the fact that tr (ABC) = tr (CAB) and summing over s' quantities FIM i init , the total FIM is given by,

FI M init ( n p , x p , I p ) = E h i = 1 s R i - 1 ( ( X p [ i ] ) H X p [ i ] ) + R i - 1 Q i + s σ h - 2 I n t n r (10)

where X p [ i ] = τ I p x p [ τ ] Φ p [ i ] [ τ ] , and Q i = ( A [ i ] ) T H T R i - 1 H * ( A [ i ] ) * .

For designing training sequence, (10) can be simplified further using numerical calculation. Using numerical calculation, it is observed that for a Rayleigh-distributed channel, the matrix E h { R i - 1 } in (10) is approximately a diagonal matrix(b), α I n r . This observation means that Eh {Qi} can be approximated by n r α σ h 2 ( A [ i ] ) T ( A [ i ] ) * . Then, by performing the expectation operation and using the factorization property of the Kronecker product, (10) can be represented as

E h i = 1 s α I n r ( X p [ i ] ) H X p [ i ] + α I n r n r α σ h 2 ( A [ i ] ) τ ( A [ i ] ) * + s σ h - 2 I n r I n t = I n r i = 1 s α ( X p [ i ] ) H X p [ i ] + n r α 2 σ h 2 ( A [ i ] ) τ ( A [ i ] ) * + s σ h - 2 I n t

Moreover, using the property of the Kronecker product (A B)(C D) = (AC) ⊗ (BD), it follows that ( X p [ i ] ) H X p [ i ] = τ I p τ I p ( ( x p [ τ ] ) H x p [ τ ] ) ( ( Φ p [ i ] [ τ ] ) H Φ p [ i ] [ τ ] ) . Therefore (10) can be further simplified to

FI M init ( n p , x p , I p ) = I n r α τ I p τ I p ( x p [ τ ] ) H x p [ τ ] i = 1 s ( Φ p [ i ] [ τ ] ) H Φ p [ i ] [ τ ] + n r α 2 σ h 2 i = 1 s ( A [ i ] ) τ ( A [ i ] ) * + s σ h - 2 I n t (11)

In general, the second term in (11) depends on I p , but not on the training symbols, whereas the first term depends on both xp and I p . Although both terms depend on np, how FIMinit depends on np is determined by I p . Therefore, in the following I p and xp are first optimized. Then np is determined for the optimized I p .

For the iteration phase, specifically the last iteration, estimation and detection are implemented using information about the data symbols as well as the pilot symbols. Thus, the parameter of interest in deriving FIM is θ = [ h T x d ] T . Moreover, μ i = E φ { y [ i ] | θ } = I n r τ x [ τ ] Φ [ i ] [ τ ] h and R i = N 0 I n r . By replacing θ in (6) for h and after some manipulations, the FIM for channel estimation in the iteration phase is given by

FI M iter ( n p , x p , I p ) = N 0 - 1 I n r ( N s - n p N σ x 2 I n t + τ I p τ I p ( x p [ τ ] ) H x p [ τ ] i = 1 s ( Φ p [ i ] [ τ ] ) H Φ p [ i ] [ τ ] ) + s σ h - 2 I n t n r . (12)

#### 3.2 Optimization of training symbols and their positions

This section is first concerned with minimizing the CRB expression for the initialization phase. The minimization is under a constraint on the power budget for the training sequence. Such a constraint is expressed as

τ I p x p [ τ ] Φ p [ i ] [ τ ] ( x p [ τ ] ) H ( Φ p [ i ] [ τ ] ) H P t . (13)

Using the properties of the precoder employed in this study, the above constraint can be simplified to s N τ I p x p [ τ ] ( x p [ τ ] ) H P t . The other obvious constraint is that the training symbols should be selected from QAM constellation Ω. Then, the training symbols, xp's and their positions, specified by I p , are obtained by solving the following constrained optimization problem:

min x p , I p CR B init ( n p , x p , I p ) = min x p , I p tr ( FI M init ( n p , x p , I p ) - 1 ) s . t . s N τ I p x p [ τ ] ( x p [ τ ] ) H P t ( x p [ τ ] ) j Ω , j = 1 , , n t / s , τ I p (14)

where ( x p [ τ ] ) j is the jth pilot symbol in the τth pilot nucleo and the FIM is given in (11).

To proceed, lets consider two separate cases for problem (14): np = 1 and np ≥ 2. Case 1 (np = 1): In this case the FIM is simplified to

I n r α ( x p [ τ ] ) H x p [ τ ] i = 1 s ( Φ p [ i ] [ τ ] ) H Φ p [ i ] [ τ ] + n r α 2 σ h 2 i = 1 s ( A [ i ] ) τ ( A [ i ] ) * + s σ h - 2 I n t , (15)

Because of the shift-invariant property of (15) with respect to τ, τ can be any value in the set {1, 2, . . . , Ns'}. For simplicity, set τ = 1 and the superscript τ is omitted. Using the fact that if X > 0 then tr (X-1) ≥ ∑i 1/(X)i,i, the original optimization problem is simplified by minimizing the lower bound of the objective function.

On the other hand, i = 1 s Φ p [ i ] [ τ ] H Φ p [ i ] [ τ ] = 1 N I s , i = 1 s ( A [ i ] ) T ( A [ i ] ) * = σ x 4 s ( N s - 1 N ) 2 + ( 1 N ) 2 I n t and the constraint is s N x p x p H = s N j = 1 n t / s | ( x p ) j | 2 . Therefore, it is not hard to see that the solution of the simplified optimization problem is | ( x p ) 1 | 2 = | ( x p ) 2 | 2 = = | ( x p ) n t / s | 2 = N P t n t . It means that all pilot symbols should have the same power. For example, one can select corner points of the QAM constellations for the training symbols.

Case 2 (np ≥ 2): In this case there are two options for the placements of pilot nucleos. The first option is to group all pilot nucleos in one single cluster and the second option is to spread pilot nucleos. It can be shown that the CRB is invariant with respect to a shift of the placements of pilot nucleos in both options. Therefore, it suffices to select one cluster or one spread placement. However, the precoder has been designed such that the soft-output demodulator works with uncorrelated inputs and putting pilot nucleos between data nucleos may violate this condition. That condition is satisfied when A[i] has a diagonal form. The implication of this property is to place pilot nucloes equi-spaced in xk and I p = { i 0 + k n ; k = 0 , , n p - 1 } , where n = Ns'/np and i0 ∈ {1, . . . , n}, which leads to A [ i ] = σ x 2 N s - n p N s I s . In this selection it is supposed that np is divisible by Ns'.

Then the FIM in (11) can be represented by

I n r 1 N α τ I p ( x p [ τ ] ) H x p [ τ ] I s + n r α 2 σ x 4 σ h 2 1 s ( N s - n p N ) 2 I n t + s σ h - 2 I n t (16)

To obtain the above expression of the objective function, the following property has been used:

i = 1 s ( Φ p [ i ] [ τ ] ) H Φ p [ i ] [ τ ] n , l = 1 N , τ = τ ; n = l 0 , otherwise (17)

Moreover, the only term that depends on the training symbols is τ I p ( x p [ τ ] ) H x p [ τ ] in (16). Finally, using the constraint on training power, which can be written as

s N τ I p j = 1 n t / s | ( x p [ τ ] ) j | 2 P t , (18)

the solution is given by τ | ( x p [ τ ] ) j | 2 = N P t n t ; j = 1 , , n t / s .

Now consider the training design for the iteration phase. Observe that all the terms in (12) have diagonal forms with equal diagonal elements, except τ I p τ I p ( x p [ τ ] ) H x p [ τ ] i = 1 s ( Φ p [ i ] [ τ ] ) H Φ p [ i ] [ τ ] . This means that the solution of problem (14), but with FI M init ( n p , x p , I p ) replaced by FI M iter ( n p , x p , I p ) , is to choose equal diagonal elements for this term. Therefore, the training sequence designed for the initialization is also optimal for the iteration phase.

In summary, by selecting pilot nucleos such that the sum of the powers of their corresponding pilot symbols with the same indexes are equal, the bound on CRB is minimized. The above condition can give different selections for pilot symbols from a two-dimensional constellation. It should be pointed out, however, that not all selections guarantee that pilot symbols belong to standard QAM constellations.

#### 3.3 Determination of the number of the training symbols

For block-fading channels, the number of pilot nucleos, i.e., np, should be as small as possible that meets the power constraint. Using a larger value for np wastes bandwidth and does not change the system performance.

The optimum numbers of the training symbols in the initialization phase and iteration phase are not the same. This is explained as follows. At the initialization, by looking at (7), it is observed that the first term in (11) is an increasing function of np. However, the second term is a decreasing function of np that is multiplied by nr. Therefore, np that minimizes the CRB are determined by the summation of these two terms, which is also determined by the value of nr. Table 1 gives several examples of optimal np for different sets of nt, nr and N. For the iteration phase, the expression in (12) means that the CRB in the iteration phase always increases by increasing np. Since it is assumed that there is perfect information about the data symbols in the iteration phase, which is not the case in reality, it is most appropriate to select np considering only the initialization phase.

Table 1. Optimum np for several sets of parameters {nt, nr, N}

To demonstrate the optimal training design, Figure 3 shows a graphical structure for a simple example, where P t = 4 σ x 2 , n p = 2 , N = 2, nt = 4 and nr = 2. In this example, ns = 1. Then the size of pilot nucleos should be nt/s' = 2, where s' = N/ns = 2.

Figure 3. Structure of the proposed scheme for the training sequence-when N = 2, nt = 4, nr = 2 and np = 2

#### 3.4 Channel estimation

For the channel estimation task, one can view the received vector during one block length as φ [ t ] = [ ( y [ 1 , t ] ) T , ( y [ 2 , t ] ) T , , ( y [ s , t ] ) T ] T .

At the initialization, the mean and covariance matrix of this vector are given in Section 3.1. By treating the data symbols as nuisance parameters, the MMSE channel estimate can be found as [14]

h ^ [ t ] = σ h 2 T H ( σ h 2 T T H + R φ [ t ] ) φ [ t ] (19)

where T = [(T1)T, . . . , (T[s'])T]T.

In the subsequent iterations, soft information from the decoder is used to improve the performance of the channel estimator. The channel estimator uses such information to compute new estimates of the channel coefficients using expected values of the data symbols. Therefore, the interleaved { Λ ext ( c l ) } from the decoder are fed back to the estimator to calculate the expected values of the data symbols, i.e., E{xd}. The entries of E{xd} are calculated using { Λ ap ( c ̃ l ) } at each iteration by E{(xd)i} = ∑x∈Ω x · p((xd)i = x). The detailed derivations of the probability p((xd)i = x) from Λ-values are given in [15] (note that the calculation depends on the mapping rule in Ω).

To verify the results obtained in this section, Section 4 compares numerically the MSE performance of the above channel estimator obtained with the optimal and suboptimal training sequences.

### 4 Illustrative results

In this section, the frame-error-rate (FER) and MSE performances of BICM-MIMO systems using a MMSE iterative channel estimator are presented. The space-time precoder is the DNA-cyclo precoder that satisfies the properties outlined in Section 2. We consider quadrature phase-shift keying (QPSK) modulation with Gray mapping.

The MSE performance of a BICM-MIMO for a codeword length of 4 × 1024 bits is shown for a 4 × 2 block-fading MIMO channel in Figure 4, when nc = 2. In this figure, Eb is the energy per information bit. The code used is the 16-state convolutional code with generator polynomials (23, 35) in octal form. In Figure 4, the MSE curves are obtained after 1 and 5 iterations of the iterative channel estimation/demodulation/decoding, with the following cyclotomic rotator [16]:

Figure 4. Comparison of MSE performance obtained with the optimal PPSAM and the sub-optimal PPSAM-over a 4 × 2 block-fading channel with nc = 2, when N = 2 and np = 2 after 1 and 5 iterations of iterative channel-estimation/demodulation/decoding.

Φ = 1 2 1 1 e j 6 π / 15 - e j 6 π / 15 e j 2 π / 15 j e j 2 π / 15 - e j 8 π / 15 j e j 8 π / 15 e j 4 π / 15 - e j 4 π / 15 e j 10 π / 15 e j 10 π / 15 e j 6 π / 15 - j e j 6 π / 15 - e j 12 π / 15 - j e j 12 π / 15

and when the setting for N, ns, np and Pt in Figure 3 are used. The channel is generated randomly and is assumed to be Rayleigh distributed. For the purpose of comparison, the results for MSE performances of the optimal PPSAM, denoted by O-PPSAM and the suboptimal PPSAM, denoted by SO-PPSAM as well as the CRB are shown in Figure 4. For SO-PPSAM, two pilot nucleos are inserted as one cluster in front of data nucleos in a symbol to be precoded. In contrast, in the case of O-PPSAM, the optimized training sequence embeds the pilot nucleos at the first and third positions of Ns' = 4 positions for nucleos. The MSE curves show that the performance of the optimal scheme is better than the sub-optimum scheme for the first iteration (i.e., initialization). In fact the MSE performance of the proposed scheme closely approaches the CRB at high Eb/N0 after 5 iterations.

In Figure 5, the FER performance of the system with the PPSAM schemes is compared with the conventional PSAM training scheme for the same system parameters as in Figure 4. The top curve is the FER performance of the system with the conventional PSAM training scheme. Note that for a fair comparison, the training scheme in PSAM also meets the training power constraint as trace ( X p X p H ) = P t , where Xp is the training matrix placed at the beginning of each block of the precoded symbols. The optimal option for PSAM scheme in terms of minimizing the FER as proposed in [11] is to select Xp to have orthogonal columns. The simplest option is 2 × σ x 2 / n t I n t = σ x 2 I n t , which results in the same power budget as that of the proposed scheme.

Figure 5. Comparison of FER performance obtained with the optimal PPSAM, sub-optimal PPSAM and PSAM scheme-over a 4 × 2 block-fading channel with nc = 2, when N = 2 and np = 2 after 1 and 5 iterations of iterative channel-estimation/demodulation/decoding.

As can be seen from Figure 5, the O-PPSAM scheme offers 0.5 dB performance gain as compared to the SO-PPSAM scheme at FER = 10-2. In comparison with PSAM, the performance of the PSAM scheme is about 0.5-1.5 dB worse than the proposed scheme depending on Eb/N0 after 5 iterations. This is expected because the pilot information is embedded in the precoded symbols for the proposed scheme and not for the PSAM scheme. In this way, the demodulator can also make use of this information. Note, however, that for the first iteration, since there is no information about data, PSAM works the best. More importantly, while the proposed scheme uses a little bandwidth for training information (for the system considered in this figure the training overhead is np × nt/s' = 4), the training overhead of PSAM scheme is nt × nt = 16, which is quadruple. To investigate the effect of the number of transmit antennas, two different systems, one with 2 × 2 channel and one with 4 × 2 MIMO channel, are compared in Figures 6 and 7 in terms of MSE and FER, respectively. For both channels, np = 2 and the optimum scheme are used when N = 2, while other system parameters are the same as those used for Figure 4. As can be seen from Figure 6, the MSE of the channel estimation increases when increasing the number of transmit antennas. This is expected because there are more channels to be estimated for the same amount of training information and power as done in the comparison. Nevertheless, the gain in diversity by using more antennas can still improve the overall FER performance as seen in Figure 7.

Figure 6. Comparison of MSE performance obtained with the optimal PPSAM for 2 × 2 and 4 × 2 block-fading channels with nc = 2, when N = 2 and np = 2 after 5 iterations of iterative channel-estimation/demodulation/decoding.

Figure 7. Comparison of FER performance obtained with the optimal PPSAM for 2 × 2 and 4 × 2 block-fading channels with nc = 2, when N = 2 and np = 2 after 5 iterations of iterative channel-estimation/demodulation/decoding.

### 5 Conclusion

In this article, a new training design for a BICM-MIMO system over a block-fading channel has been proposed. The design inserts pilot symbols into the data symbols before precoding. The new training sequence improves bandwidth efficiency as compared to the conventional PSAM scheme and can also be used by the demodulator in the receiver. In order to design the optimal training symbols and their positions, the CRB on the channel estimations at the initialization and at the iteration phases are minimized. Compared to PSAM, performance improvement achieved with the proposed training is about 1.5 dB at a FER level of 10-2.

### Endnotes

aIn practice, since ns is typically an approximated value over some range and since N can be selected, such an assumption can be fulfilled. bUsing the matrix inversion lemma, one has R i - 1 = ( H A [ i ] H H + N 0 I n r ) - 1 = N 0 - 1 I n r + N 0 - 2 H A [ i ] H H ( I n r + N 0 - 1 H A [ i ] H H ) - 1 . Therefore, for high SNR, E { R i - 1 } can be approximated by N 0 - 1 I n r .

### Competing interests

Zohreh Andalibi has received funding from TRLabs of Saskatchewan. This organization partially is financing this manuscript.

### References

1. G Caire, S Shamai, On the achievable throughput of a multiantenna Gaussian broadcast channel. IEEE Trans Inf Theory 49(7), 1691–1706 (2003). Publisher Full Text

2. SM Alamouti, A simple transmit diversity technique for wireless communications. IEEE J Sel Areas Commun 16(8), 1451–1458 (1998). Publisher Full Text

3. V Tarokh, N Seshadri, AR Calderbank, Space-time codes for high data rate wireless communication: performance criterion and code construction. IEEE Trans Inf Theory 44(2), 744–765 (1998). Publisher Full Text

4. J Boutros, E Viterbo, Signal space diversity: a power and bandwidth eficient diversity technique for the Rayleigh fading channel. IEEE Trans Inf Theory 44(4), 1453–1467 (1998). Publisher Full Text

5. J Boutros, N Gresset, L Brunel, Turbo coding and decoding for multiple antenna channels. International Symposium on Turbo Codes and Related Topics (Brest, France, 2003), pp. 1–8

6. N Gresset, L Brunel, J Boutros, Space-time coding techniques with bit-interleaved coded modulations for MIMO block-fading channels. IEEE Trans Inf Theory 54(5), 2156–2178 (2008)

7. N Gresset, JJ Boutros, L Brunel, Optimal linear precoding for BICM over MIMO channels. ISIT, 66 (Chicago, IL, 2004)

8. M Coldrey, P Bohlin, Training-based MIMO systems, Part I: performance comparison. IEEE Trans Signal Process 55(11), 5464–5476 (2007)

9. M Nicoli, S Ferrara, U Spagnolini, Soft-iterative channel estimation: methods and performance analysis. IEEE Trans Signal Process 55(6), 2993–3006 (2007)

10. M Dong, L Tong, BM Sadler, Optimal insertion of pilot symbols for transmissions over time-varying flat fading channels. IEEE Trans Signal Process 52(5), 1403–1418 (2004). Publisher Full Text

11. G Taricco, E Biglieri, Space-time decoding with imperfect channel estimation. IEEE Trans Wirel Commun 4(4), 1874–1888 (2005)

12. Y Huang, JA Ritcey, Joint iterative channel estimation and decoding for bit-interleaved coded modulation over correlated fading channels. IEEE Trans Wirel Commun 4(5), 2549–2558 (2005)

13. P Piantanida, SM Sadough, On the outage capacity of a practical decoder accounting for channel estimation inaccuracies. IEEE Trans Commun 57(5), 1341–1350 (2009)

14. SM Kay, Fundamentals of Statistical Signal Processing: Estimation Theory (Prentice-Hall PTR, New Jersey, 1993)

15. MA Khalighi, JJ Boutros, Semi-blind channel estimation using the EM algorithm in iterative MIMO APP detectors. IEEE Trans Wirel Commun 5(11), 3165–3173 (2006)

16. GM Kraidy, P Rossi, Full-diversity iterative MMSE receivers with space-time precoders over block-fading MIMO channels. Proc IEEE Int Conf Wireless Commun and Signal Processing (Suzhou, 2010), pp. 1–5