# Multiantenna Analog Network Coding for Multihop Wireless Networks

Authors: Ramesh Annavajjala; Amine Maaref; Jinyun Zhang
Journal: International Journal of Digital Multimedia Broadcasting
 (2010)

---

## Abstract

This paper proposes a two-phase minimum mean-square-error bidirectional
amplify-and forward (MMSE-BAF) relaying protocol to allow two sources
exchange independent messages via a relay node equipped with multiple antennas.
MMSE-BAF performs a joint linear MMSE filtering of the received signal after
the multiple access phase before amplifying and forwarding the filtered signal
using a single transmit antenna, possibly through a specific antenna selection
procedure, during the broadcast phase. The proposed protocol extends upon the
so-called analog network coding schemes in the literature in that it inherently
exploits the multiple antennas at the relay station to reduce the noise enhancement
typical of an AF protocol, and can also compensate for link imbalances between
the relay and the sources and is agnostic to sources' modulation and coding
schemes. We derive the instantaneous signal-to-noise ratio expressions for the
received signal by the sources in the downlink and provide extensive link-level
simulations for the MMSE-BAF protocol subject to both frequency flat and
selective fading. Furthermore, we pinpoint the modifications to be incorporated
into the IEEE 802.16e orthogonal-frequency-division multiple access (OFDMA)
cellular standard (mobile WiMax) to enable support of multiantenna bidirectional
communications and show that MMSE-BAF is a viable solution within that
framework.


---

## Body

## 1. Introduction

Half-duplex bidirectional relay systems in which two nodesS1 and S2 wish to exchange independent messages via a third node R, termed relay, give rise to some interesting challenges from a cooperative communications and information-theoretic points of view.  This is especially true when the relay node R is equipped with multiple antennas. Such two-way relay channels have many applications in ad hoc and cellular networks in which all mobile-to-mobile communications have to pass through a common base station. Since full-duplex operation is of little practical interest given current state-of-the-art technology, our focus is on half-duplex nodes, where each active node can either transmit or receive an information message at a given point in time. In particular, without loss of generality, we are interested in the communications part of the problem in a cellular context where two mobile stations wish to exchange data simultaneously via a common base station.The traditional baseline approach for bidirectional communications in half-duplex mode between two sourcesS1 and S2 via a relay station R consists of a 4-phase protocol with a completion time of 4-time-slots (TSs) whereby S1 and S2 send N-bit packets b1 and b2∈{0,1}N to R during TS 1 and TS 2, respectively; R decodes the received packets and then sends b2 to S1 and b1 to S2 during TS 3 and TS 4, respectively. The gist of the 4-phase protocol is to avoid interference by preventing simultaneous transmissions from the sources to the relay and vice versa. However, it was shown in [1] that a three-phase protocol exploiting the so-called network coding concept by decoding packets b1 and b2 at the relay and then broadcasting a single packet b1⊕b2, where ⊕ denotes the bit-wise exclusive-or (XOR) operation, is actually more attractive in terms of achievable throughput, since the desired packet at S1 can be decoded using another XOR operation (and similarly at S1). Better still, a recent concept introduced in [2] and termed analog network coding (ANC) combines the first two phases of the conventional baseline protocol into a single multiple access (MA) phase with simultaneous transmissions from the sources to the relay; the received multiple access signal at the relay is then amplified and broadcast to S1 and S2, thereby yielding the so-called two-phase bidirectional amplify-and-forward (BAF) protocol. A similar concept to ANC, using estimate-and-forward relaying as opposed to AF relaying, has been proposed in [3] under the terminology of physical-layer network coding (PNC). A schematic diagram illustrating the aforementioned bidirectional protocols is illustrated in Figure 1.Bidirectional communication protocols: (a) 4-phase conventional protocol, (b) 3-phase protocol with decode-and-forward network coding, and (c) 2-phase BAF protocol with amplify-and-forward relaying.
(a)(b)(c)Information-theoretic aspects such as bounds on the achievable throughput and the capacity region of the bidirectional relay channel have been investigated in [4–6]. A Markov-chain-based performance analysis for several variants of the BAF protocol was carried out in [7–9]. Linear beamforming filter designs for bidirectional communications with multiantenna relay stations are proposed in [10–12]. In this paper, we propose a simple two-phase minimum-mean-square-error- (MMSE-) BAF protocol which operates by filtering the received signal at the relay station after the MA phase using a specially designed joint linear MMSE filter before amplifying and forwarding the filtered signal during the broadcast phase. Whereas optimal relay beamforming structure for bidirectional multiantenna relay channels is sought in [12], our MMSE-BAF protocol is a simple low-complexity driven approach for bidirectional multiantenna relay channels which exploits linear signal processing on the uplink (from S1 and S2 to R) and transmit antenna selection (TAS) on the downlink (from R to S1 and S2). Besides, MMSE-BAF is a two-phase bidirectional relaying protocol, whereas the multiantenna relaying protocol put forward in [10] is a three-phase protocol which necessitates decoding and reencoding of the received signals at the RS prior to the broadcast phase. Finally, MMSE-BAF differs from the so-called spatial division duplex (SDD) bidirectional relaying scheme proposed in [11] in that it allows to bias the beamforming weights in favor of one of the two source nodes in as much as required to compensate for potential imbalance of the relay-to-source link channel gains or other parameters such as dissimilar signal constellations employed at S1 and S2.The remainder of this paper is structured as follows. Section2 presents the system model for MMSE-BAF relaying. In Section 3, we describe the proposed MMSE-BAF protocol and derive the signal-to-noise ratio (SNR) expressions upon which our extensive simulation results provided in Section 4 are based. Finally, concluding remarks are drawn in Section 5.


## 2. System Setup and Signal Model

The following set of notations is employed throughout this paper: Boldface upper- and lower-case symbols are used to denote matrices and column vectors, respectively.Im denotes the identity matrix of order m. Moreover, (·)*,(·)T,(·)ℋ, and E[·] stand for conjugate, transpose, transpose-conjugate, and expectation operators, respectively.Without loss of generality, and although our proposal and framework are applicable to any type of two-hop bidirectional relay setting, we focus our attention on cellular systems. For that purpose, we consider an infrastructure-based wireless communications system consisting of two mobile stations (MSs),MS1 and MS2, and one base station (BS). A block diagram of the system under consideration is depicted in Figure 2. Both MSs as well as the BS are equipped with multiple antennas for reception with the aim of canceling out potential other cell/user interference but are restricted to transmit through a single transmit antenna. This assumption is dictated by the need to reduce the transmit-power requirements for user terminals and to lower the complexity and cost of a transmission chain at the base station, generally higher than that of a reception chain, especially when accounting for high-cost radio-frequency amplifiers involved in the transmission chain. This is, for instance, the case in current cellular standards such as IEEE 802.16e [13]. Note that our proposed scheme works equally well with both time-division duplex (TDD) as well as frequency-division duplex (FDD) modes of operation. Without loss of generality and for the sake of notational brevity, we focus on the FDD mode in the following analysis. Performance results for both TDD and FDD will be presented in Section 4.Block diagram of the system setup.
(a)
Uplink multiple access phase(b)
Dowlink broadcast phaseComplex baseband transmission is assumed throughout the paper. LetMbs denote the number of receive antennas at the BS, and let h1[n]=[h11[n],h12[n],…,h1Mbs[n]]T and h2[n]=[h21[n],h22[n],…,h2Mbs[n]]T denote the Mbs×1 uplink channels from MS1 and MS2 to the BS, respectively, where n is the discrete-time index. The corresponding time-varying channel vector elements {hij[n]}i=1,2;j=1,…,Mbs are realizations of a zero-mean unit-variance Gaussian wide-sense stationary process. Assuming Mss antennas at each of MS1 and MS2 for downlink reception, we define g1[n]=[g11[n],g12[n],…,g1Mss[n]]T and g2[n]=[g21[n],g22[n],…,g2Mss[n]]T as the downlink channels from BS to MS1 and MS2, respectively. For the special case of Mbs=Mss=1, TDD assumption allows us to set g1[n]=h1T[n] and g2=h2T[n]. Let x1[n] be the signal transmitted from MS1 and intended for MS2 and let x2[n] be the signal transmitted from MS2 and intended for MS1 at time n. Both x1[n] and x2[n] are drawn from two possibly different complex signal constellations with average energies σ12=E[|x1[n]|2] and σ22=E[|x2[n]|2], respectively. Prior to any signal processing at the BS, the Mbs×1-received signal at the end of the MA phase is given by
(1)y[n]=h1[n]x1[n]+h2[n]x2[n]+n[n],
where n[n] is an Mbs×1 additive white Gaussian noise vector at the BS which is modeled as a zero-mean circularly symmetric Gaussian random vector with covariance matrix E[n[n]n[n]ℋ]=σN2IMbs. For notational simplicity, the discrete-time index n is henceforth omitted.


## 3. MMSE-BAF: Description and Analysis

### 3.1. Protocol Description

Figure3 illustrates the physical-layer frame structure for enabling the MMSE-BAF protocol for the uplink MA and the downlink broadcast phases. A frame consists of packets originating from the link-layer whose size depends on the chosen MCS so that each downlink or uplink frame contains a fixed number of symbols. As can be seen in the left-hand side of Figure 3, the uplink frame structure is composed of two parts, one for pilot symbols which are chosen to be orthogonal for MS1 and MS2. Orthogonality of the pilot symbols can be maintained in the time, frequency, or two-dimensional time-frequency grid. The second part is for data symbols. The orthogonal pilot symbols are used to estimate the channels h1 and h2 corresponding to MS1 and MS2. A beamforming weight vector wopt is then computed at the BS based on a joint MMSE criterion to be specified shortly. The BS then estimates an amplification factor β subject to an average power constraint. Likewise, the downlink frame structure contains pilot and data parts. Additionally, it contains a control part consisting of quantized versions of the amplification factor β (a positive scalar value) and two complex scalars values v1:=woptℋh1 and v2:=woptℋh2. The downlink data symbols consist of the beamformed amplified-and-forwarded symbols received on the uplink frame in FDD mode.Block diagram of the system setup.
(a)(b)According to the MMSE-BAF protocol, the BS performs the following set of operations.(1) Jointly minimize the MSE between the received signal y[n] and the transmitted signals x1[n] and x2[n], thus performing a joint linear-MMSE filtering of the received signal, using the following metric:
(2)wopt=argminw∈ℂMbs{δ1E[|x1-wℋy|2∣h1,h2]+δ2E[|x2-wℋy|2∣h1,h2]},
where ℂ is the field of complex numbers and δ1≥0, δ2≥0, δ1+δ2=1, are the two design constants that control the relative weight assigned to the signals of MS1 and MS2. The minimization problem in (2) is a modified Wiener filtering problem whose solution can be easily found using the orthogonality principal in linear mean square estimation and is given by
(3)wopt=(σ12h1h1ℋ+σ22h2h2ℋ+σN2IMbs)-1(δ1σ12h1+δ2σ22h2).
This minimization requires an estimation of both mobile stations’ vector-valued channels h1 and h2.(2) Amplify the linear MMSE-filter output to maintain a constant average transmit power PT which leads to computing the amplification gain factor:
(4)β=PTE[|wℋy|2∣h1,h2]=PTσ12|woptℋh1|2+σ22|woptℋh2|2+σN2∥wopt∥2.(3) Transmit the amplified signal back to the MSs on one of the antennas using an appropriate downlink transmit antenna-selection (TAS) algorithm, based on the uplink channel. One approach inherent to the MMSE-BAF protocol is to select the antenna that has the largest beamformer weight.  We note that, owing to the channel reciprocity, TAS using the largest beamformer weight is applicable for TDD only.


### 3.2. Performance Analysis

Definez:=woptℋy, where y is the uplink received signal (1), as the output of the MMSE filtering operation at the BS. The AF transmitted signal on the downlink is
(5)xr=βz=βwoptℋy.
The received signal on the downlink for MS1 and MS2 is therefore given by
(6)y1=g1xr+n1,(7)y2=g2xr+n2,
where n1 and n2 are the zero-mean AWGN at MS1 and MS2, respectively, with covariance matrix σN2IMss. Without loss of generality, let us focus on the signal received by MS1 (6). A similar signal processing is required at MS2. Incorporating (1) and (5) into (6) yields
(8)y1=g1βwoptℋ(h1x1+h2x2+n)+n1.
Now, assuming that MS1 (i) is able to perfectly estimate its own downlink channel vector g1 owing to the downlink pilot symbols sent by the BS, (ii) knows its own transmitted signal x1, and (iii) is able to extract the value of the amplification factor β as well as the couple (v1,v2), all of which are sent on the downlink control channel as explained in Section 3.1, then MS1 can subtract the so-called self-interference component of the received signal, that is, g1βv1x1, hence winding up with a processed received signal of the form
(9)r1=y1-g1βv1x1=y1-g1βwoptℋh1x1=g1βwoptℋh2x2+g1βwoptℋn+n1︸:=ñ1.
Note that ñ1, defined in the previous equation, is a zero-mean colored noise vector with a conditional covariance matrix given by
(10)Σ1:=E[ñ1ñ1ℋ∣g1,β,wopt](11)=β2g1woptℋE[nnℋ]woptg1ℋ+E[n1n1ℋ](12)=σN2(IMss+β2g1woptℋwoptg1ℋ).
In the absence of knowledge of the conditional covariance matrix Σ1 at the receiver of MS1, an estimate x̂2 of x2 can be obtained as follows:
(13)x̂2=(g1βv2)ℋ(g1βv2)ℋ(g1βv2)r1,
which after simplification becomes
(14)x̂2=x2+1βv2g1ℋñ1g1ℋg1︸:=ň2.
The conditional variance of the newly defined scalar noise term ň2 can be expressed as
(15)Σ1′=E[|ň2|2]=1β2|v2|2g1ℋE[ñ1ℋñ1]g1(g1ℋg1)2=1β2|v2|2g1ℋΣ1g1(g1ℋg1)2.
Now, we are in a position to determine the SNR γ2 at MS1 as
(16)γ2=E[|x2|2]E[|ň2|2](17)=σ22β2|v2|2(g1ℋg1)2g1ℋΣ1g1.
Note that after incorporating (10) into (17) and further simplification, (17) can be expressed as follows:
(18)γ2=σ22β2|v2|2σN2×∥g1∥21+β2∥wopt∥2∥g1∥2.
Similarly, one can evaluate the SNR γ1 for the signal x1 received at MS2 which is found to be
(19)γ1=σ12β2|v1|2(g2ℋg2)2g2ℋΣ2g2,
where Σ2 is a noise covariance matrix (analogous to Σ1) defined as
(20)Σ2:=σN2(IMss+β2g2woptℋwoptg2ℋ).
Again, upon incorporation of (20) into (19), the latter can be simplified to
(21)γ1=σ12β2|v1|2σN2×∥g2∥21+β2∥wopt∥2∥g2∥2.It is worthwhile to mention that the above SNR expressions forγ2 (18) and γ1 (21) have been obtained without exploiting the colored nature of Σ1 and Σ2. Somewhat surprisingly, and as indicated in the appendix, even through whitening of the colored noise, the SNR expressions for γ1 and γ2 remain the same, which is not so bad in itself, as it means that the signal processing cost associated with the whitening operation can be completed eliminated.


## 3.1. Protocol Description

Figure3 illustrates the physical-layer frame structure for enabling the MMSE-BAF protocol for the uplink MA and the downlink broadcast phases. A frame consists of packets originating from the link-layer whose size depends on the chosen MCS so that each downlink or uplink frame contains a fixed number of symbols. As can be seen in the left-hand side of Figure 3, the uplink frame structure is composed of two parts, one for pilot symbols which are chosen to be orthogonal for MS1 and MS2. Orthogonality of the pilot symbols can be maintained in the time, frequency, or two-dimensional time-frequency grid. The second part is for data symbols. The orthogonal pilot symbols are used to estimate the channels h1 and h2 corresponding to MS1 and MS2. A beamforming weight vector wopt is then computed at the BS based on a joint MMSE criterion to be specified shortly. The BS then estimates an amplification factor β subject to an average power constraint. Likewise, the downlink frame structure contains pilot and data parts. Additionally, it contains a control part consisting of quantized versions of the amplification factor β (a positive scalar value) and two complex scalars values v1:=woptℋh1 and v2:=woptℋh2. The downlink data symbols consist of the beamformed amplified-and-forwarded symbols received on the uplink frame in FDD mode.Block diagram of the system setup.
(a)(b)According to the MMSE-BAF protocol, the BS performs the following set of operations.(1) Jointly minimize the MSE between the received signal y[n] and the transmitted signals x1[n] and x2[n], thus performing a joint linear-MMSE filtering of the received signal, using the following metric:
(2)wopt=argminw∈ℂMbs{δ1E[|x1-wℋy|2∣h1,h2]+δ2E[|x2-wℋy|2∣h1,h2]},
where ℂ is the field of complex numbers and δ1≥0, δ2≥0, δ1+δ2=1, are the two design constants that control the relative weight assigned to the signals of MS1 and MS2. The minimization problem in (2) is a modified Wiener filtering problem whose solution can be easily found using the orthogonality principal in linear mean square estimation and is given by
(3)wopt=(σ12h1h1ℋ+σ22h2h2ℋ+σN2IMbs)-1(δ1σ12h1+δ2σ22h2).
This minimization requires an estimation of both mobile stations’ vector-valued channels h1 and h2.(2) Amplify the linear MMSE-filter output to maintain a constant average transmit power PT which leads to computing the amplification gain factor:
(4)β=PTE[|wℋy|2∣h1,h2]=PTσ12|woptℋh1|2+σ22|woptℋh2|2+σN2∥wopt∥2.(3) Transmit the amplified signal back to the MSs on one of the antennas using an appropriate downlink transmit antenna-selection (TAS) algorithm, based on the uplink channel. One approach inherent to the MMSE-BAF protocol is to select the antenna that has the largest beamformer weight.  We note that, owing to the channel reciprocity, TAS using the largest beamformer weight is applicable for TDD only.


## 3.2. Performance Analysis

Definez:=woptℋy, where y is the uplink received signal (1), as the output of the MMSE filtering operation at the BS. The AF transmitted signal on the downlink is
(5)xr=βz=βwoptℋy.
The received signal on the downlink for MS1 and MS2 is therefore given by
(6)y1=g1xr+n1,(7)y2=g2xr+n2,
where n1 and n2 are the zero-mean AWGN at MS1 and MS2, respectively, with covariance matrix σN2IMss. Without loss of generality, let us focus on the signal received by MS1 (6). A similar signal processing is required at MS2. Incorporating (1) and (5) into (6) yields
(8)y1=g1βwoptℋ(h1x1+h2x2+n)+n1.
Now, assuming that MS1 (i) is able to perfectly estimate its own downlink channel vector g1 owing to the downlink pilot symbols sent by the BS, (ii) knows its own transmitted signal x1, and (iii) is able to extract the value of the amplification factor β as well as the couple (v1,v2), all of which are sent on the downlink control channel as explained in Section 3.1, then MS1 can subtract the so-called self-interference component of the received signal, that is, g1βv1x1, hence winding up with a processed received signal of the form
(9)r1=y1-g1βv1x1=y1-g1βwoptℋh1x1=g1βwoptℋh2x2+g1βwoptℋn+n1︸:=ñ1.
Note that ñ1, defined in the previous equation, is a zero-mean colored noise vector with a conditional covariance matrix given by
(10)Σ1:=E[ñ1ñ1ℋ∣g1,β,wopt](11)=β2g1woptℋE[nnℋ]woptg1ℋ+E[n1n1ℋ](12)=σN2(IMss+β2g1woptℋwoptg1ℋ).
In the absence of knowledge of the conditional covariance matrix Σ1 at the receiver of MS1, an estimate x̂2 of x2 can be obtained as follows:
(13)x̂2=(g1βv2)ℋ(g1βv2)ℋ(g1βv2)r1,
which after simplification becomes
(14)x̂2=x2+1βv2g1ℋñ1g1ℋg1︸:=ň2.
The conditional variance of the newly defined scalar noise term ň2 can be expressed as
(15)Σ1′=E[|ň2|2]=1β2|v2|2g1ℋE[ñ1ℋñ1]g1(g1ℋg1)2=1β2|v2|2g1ℋΣ1g1(g1ℋg1)2.
Now, we are in a position to determine the SNR γ2 at MS1 as
(16)γ2=E[|x2|2]E[|ň2|2](17)=σ22β2|v2|2(g1ℋg1)2g1ℋΣ1g1.
Note that after incorporating (10) into (17) and further simplification, (17) can be expressed as follows:
(18)γ2=σ22β2|v2|2σN2×∥g1∥21+β2∥wopt∥2∥g1∥2.
Similarly, one can evaluate the SNR γ1 for the signal x1 received at MS2 which is found to be
(19)γ1=σ12β2|v1|2(g2ℋg2)2g2ℋΣ2g2,
where Σ2 is a noise covariance matrix (analogous to Σ1) defined as
(20)Σ2:=σN2(IMss+β2g2woptℋwoptg2ℋ).
Again, upon incorporation of (20) into (19), the latter can be simplified to
(21)γ1=σ12β2|v1|2σN2×∥g2∥21+β2∥wopt∥2∥g2∥2.It is worthwhile to mention that the above SNR expressions forγ2 (18) and γ1 (21) have been obtained without exploiting the colored nature of Σ1 and Σ2. Somewhat surprisingly, and as indicated in the appendix, even through whitening of the colored noise, the SNR expressions for γ1 and γ2 remain the same, which is not so bad in itself, as it means that the signal processing cost associated with the whitening operation can be completed eliminated.


## 4. Simulation Results and Discussion

In this section, we present some simulation results on the performance of the proposed MMSE-BAF protocol. First, Figure4 shows the MSE performance of MMSE-BAF as a function of the relative loading of user-1 over user-2, δ1, with Mbs=4 antennas at the base-station. Two scenarios are considered: in Figure 4(a), both users average received SNRs at the base-station are set to 10 dB, whereas in Figure 4(b) the average received SNR is set to 20 dB and 40 dB for user-1 and user-2, respectively. The weighted average MSE immediately after the application of MMSE beamforming and the average MSE of user-1 and user-2 after MMSE beamforming followed by AF are obtained by drawing independent channel realizations over 100000 trials. From Figures 4(a) and 4(b), we observe that increasing δ1 minimizes the MSE of user-1 at the expense of an increase in MSE for user-2, whereas an optimum δ1 exists that jointly minimizes the MSE of both users. Interestingly, from an implementation point of view, the range of δ1 is broad to arrive at this optimum overall MSE. Comparing Figures 4(a) and 4(b), we notice that, due to unequal average received SNRs, the individual MSEs as well as the overall MSE are not symmetric functions of δ1. Thus, one should take into account the knowledge of the average uplink SNRs to arrive at an appropriate δ1 to maintain desired MSE levels for each of the two users.Average MSE of MMSE-BAF system with 4 antennas at the base-station.
(a)
Equal average uplink received SNRs(b)
Unequal average uplink received SNRsIn Figure5, the empirical cumulative distribution function (CDF) of two-time-slots-based MMSE-BAF is compared against the four-time-slot-based baseline system. Here, the base-station and the mobile station receivers are each equipped with four receive antennas. The uplink average received SNRs, per antenna, of users 1 and 2 are set to 5 dB and 10 dB, respectively, and the average received SNR per antenna at each mobile station is set to 5 dB. We also assume δ1=δ2. Figure 5 shows that the two-time-slots-based MMSE-BAF system outperforms the baseline system by an order of magnitude.Figure 5
Comparison of empirical capacity CDFs of the MMSE-BAF system against the baseline system with 4 antennas at the base-station as well as the mobile receivers.As argued earlier, the proposed MMSE-BAF protocol is equally attractive to both TDD and FDD systems. Figures6 and 7 show uncoded symbol error rate (SER) performance of MMSE-BAF on block-fading TDD channels when the two users employ dissimilar modulation formats. In both figures, (i) user-1 employs QPSK modulation whereas user-2 employs 16-QAM modulation, (ii) a data frame contains 100 modulation symbols and 20 pilot symbols for channel estimation, and (iii) the channel remains constant over the duration of at least two frames (TDD assumption). In Figure 6, the base-station and the users all have single transmit/receive antennas. With equal average received SNRs at the base-station, Figure 6(a) shows that, at lower average received SNRs, pilot-based channel estimation matches very closely the performance achieved in case of perfect channel knowledge for both users. Since each user has to subtract its own channel-compensated transmitted symbol to decode the other user's modulation symbol, a user transmitting using a higher-order constellation has the potential to generate higher self-interference in the presence of channel estimation errors. Figure 6(a) shows that with equal average received SNRs, the average SER of BPSK exhibits an error floor for an average SNR higher than 30 dB. When user-2 transmits at an SNR that is 20 dB higher than user-1’s SNR, Figure 6(b) shows that the error floor for the BPSK modulation occurs much earlier.Average uncoded SER of MMSE-BAF system with QPSK modulation for user 1 and 16-QAM modulation for user 2.
(a)
Equal average uplink received SNRs(b)
Unequal average uplink received SNRsFigure 7
Performance of transmit antenna selection for TDD-based MMSE-BAF system with 4 antennas at the base-station. User-1 employs QPSK modulation whereas user-2 employs 16-QAM modulation. An uncoded system is considered with realistic channel estimation over block fading channels with 20 pilot and 100 data symbols per fading block.The advantages of transmit antenna selection over transmitting from an arbitrary antenna is investigated in Figure7 when the base-station has four antennas for reception. Exploiting the channel reciprocity of TDD systems, we first compute the elementwise magnitude of the estimated beamformer and downlink transmission is directed from the antenna that has the highest magnitude. It is important to note that once the base-station computes the MMSE beamformer, no additional computation complexity for TAS is required. From Figure 7, we observe that the pilot-based channel estimation has excellent performance in comparison with the ideal performance and our proposed simple TAS yields an impressive gain of close to 3 dB at an average SER of 10-4.We have also investigated the feasibility of bidirectional relaying for OFDM/OFDMA-based 4G cellular standards such as IEEE 802.16e [13]. The IEEE 802.16e system is based on OFDMA physical layer for both uplink and downlink. Current mobile WiMax standard supports various subchannelization procedures, in both uplink and downlink directions, for data transmission in time (OFDM symbols) and frequency (OFDM subcarriers). One such uplink subchannelization procedure is termed partially utilized subchannelization (PUSC) wherein the modulation symbols of a given user are pseudorandomly spread over the frequency band to extract frequency diversity and to average interference across neighboring cells/sectors. Briefly, one slot in UL-PUSC is defined as 48 modulation symbols spanning over three consecutive OFDM symbols (which is a PUSC slot duration). The modulation symbols together with the pilot symbols needed to estimate the uplink channel are sent over 6 tiles distributed over frequency, where a tile is defined as four consecutive subcarriers over three consecutive OFDM symbols. Each tile contains 4 pilot symbols, placed at the corners of the tile, and 8 data symbols. An FEC block in WiMax comprises of a given number of slots and the maximum FEC block size is a function of the modulation order and channel coding rate. The WiMax standard supports 8 modulation order and coding rate combinations. These are QPSK modulation with code rates 1/2 and 3/4, 16-QAM modulation with code rates 1/2 and 3/4, and 64-QAM modulation with rates 1/2, 2/3, 3/4, and 5/6. Figure 8 shows the modified UL-PUSC structure to support bidirectional communications. Each user employs Hadamard sequences as pilot symbols to enable the base-station to estimate the individual channels without interference. For downlink transmission, the base-station can use any subchannelization procedure. However, to render our proposal valid for FDD as well as TDD systems, the downlink subchannel structure is set identical to the uplink one and the broadcast pilots from each tile are used for channel estimation at the mobile stations. For efficient cancellation of self-interference, each mobile station requires the knowledge of uplink channel-related parameters β (a positive scalar value) and the two complex scalars values v1:=woptℋh1 and v2:=woptℋh2. Additional simulation parameters are listed in Table 1.Table 1
Simulation parameters.
ParameterValueBandwidth10 MHzSampling rate11.2 MspsFFT size1024Subcarrier spacing10.9375 KHzUseful symbol duration (Tu)91.4286μsecCyclic prefix (TG)Tu/8Useful subcarriers840Left guard subcarriers92Right guard subcarriers91Channel codingConvolutional turbo coding(with 8 iterations)Carrier frequency2.0 GHzFigure 8
Subchannelization procedure for bidirectional communication using the IEEE 802.16e protocol. In uplink, each user is allocatedK slots spanning Nf subchannels in frequency over Nt slot-durations. Each slot-duration comprises 3 consecutive OFDM symbols, whereas each subchannel contains 6 tiles distributed throughout of the useful frequency band. A tile contains 4 subcarries over a slot-duration with the 4 pilots at the corners of the tile and the remaining 8 tones used for data. In the above figure, there are N tiles with Nt=1, Nf=N/6, and K=NtNf=N/6. The circles filled with violet and red colors indicate the pilot tones of users 1 and 2, respectively, whereas the circles filled with black and blue colors indicate the data tones of users 1 and 2, respectively.In Figures9 and 10, we show the performance of MMSE-BAF when applied to an IEEE 802.16e system in an FDD mode of operation. These simulations are conducted for a base-station with four receive antennas and for a two-antenna mobile station receiver. In these plots, block error rate (BLER) performances of genie-aided perfect channel knowledge are compared against pilot-based realistic channel estimation schemes. For each tile, one channel estimate is obtained by sample averaging the received pilots over that tile. Knowledge of neither the fading statistics nor the delay/Doppler spread is assumed at the receivers. In Figure 9, both users are assumed to encode their data using QPSK modulation with rate-1/2 convolutional turbo coding (CTC). For this MCS, the FEC block size is set to the maximum allowed, which is equal to 10 slots or 10×48×2×1/2=480 information symbols. We have considered ITU Vehicular-A channel model with both low and high Doppler spread values of 32 Hz and 256 Hz, respectively. Figure 9 shows that, under both high and low Doppler scenarios, pilot channel estimation incurs a loss of about 2 dB, and MMSE-BAF works extremely well in supporting data exchanges in high-mobile environments. Figure 10 shows the performance of MMSE-BAF with 64-QAM modulation and a CTC with code rate of 1/2. With this MCS level, supporting a per-user over-the-air spectral efficiency of 3 bits/sec/Hz, we have employed the maximum possible FEC block size of 2 slots, or 2×48×6×1/2=288 information bits. Figure 10(a) shows the block error performance over ITU-Vehicular-A channel with 32 Hz Doppler, whereas Figure 10(b) shows the performance on a Pedestrian-B channel with 6 Hz Doppler. Due to higher-frequency-selectivity of Pedestrian-B channel, compared with the Vehicular-A channel where the channel estimation-based BLER performance (at 1 percent BLER) is about 2 dB away from the ideal performance, the simple sample-average-based channel estimation has a performance degradation of about 3.5 dB compared to the performance with perfect channel knowledge. It is expected that the BLER performance can be significantly improved by incorporating a more complex two-dimensional channel estimation scheme, such as Wiener filtering, which requires knowledge of fading statistics as well as delay and Doppler spread information.Performance of MMSE-BAF on MIMO-OFDMA-based IEEE 802.16e system using uplink partial utilization of subchannelization (PUSC) permutation. Each user employs rate-1/2 convolutional turbo coding (CTC) with QPSK modulation. The FEC block length is 480 information bits which corresponds to 10 slots, with 48 modulation symbols per slot, as per the terminology in [13].
(a)
With 32 Hz Doppler(b)
With 256 Hz DopplerPerformance of MMSE-BAF on MIMO-OFDMA-based IEEE 802.16e system using uplink partial utilization of subchannelization (PUSC) permutation. Each user employs rate-1/2 convolutional turbo coding (CTC) with 64QAM modulation. The FEC block length is 288 information bits which corresponds to 2 slots, with 48 modulation symbols per slot, as per the terminology in [13].
(a)
ITU-Veh-A channel with 32 Hz Doppler(b)
ITU-Ped-B channel with 6 Hz Doppler


## 5. Conclusion

In this paper, we have introduced a so-called MMSE-BAF protocol for bidirectional communications over two-way relay channels with multiantenna relay nodes. The features of this protocol include the usage of analog network coding at the relay node and the evaluation of a receive weight vector for the relay node using a joint linear MMSE filtering operation on the received uplink multiple access signal. Transmit antenna selection using the largest MMSE weight branch on the downlink is also an inherent feature of the proposed protocol in the TDD mode of operation. Extensive link-level simulations have been proposed for both TDD and FDD modes of operations and required modifications to the existing IEEE 802.16e standard have been proposed to accommodate the MMSE-BAF protocol. It has been shown through simulation results that the MMSE-BAF protocol is a simple yet efficient solution to the problem of bidirectional communications in two-way relay channels with multiantenna relays and half-duplex nodes.



---

*Source: 368562-2010-01-19.xml*