Bidirectional algorithms for interference suppression in multiuser systems

Abstract

This paper presents adaptive bidirectional minimum mean-square error parameter estimation algorithms for fast-fading channels. The time correlation between successive channel gains is exploited to improve the estimation and tracking capabilities of adaptive algorithms and provide robustness against time-varying channels. Bidirectional normalized least mean-square and conjugate gradient algorithms are devised along with adaptive mixing parameters that adjust to the time-varying channel correlation properties. An analysis of the proposed algorithms is provided along with a discussion of their performance advantages. Simulations for an application to interference suppression in multiuser DS-CDMA systems show the advantages of the proposed algorithms.

Clarke and de Lamare EURASIP Journal on Wireless Communications and Networking (2015) 2015:228DOI 10.1186/s13638-015-0440-0

http://crossmark.crossref.org/dialog/?doi=10.1186/s13638-015-0440-0-x&domain=pdf

Web End = Patrick Clarke1 and Rodrigo C. de Lamare1,2*

Abstract

This paper presents adaptive bidirectional minimum mean-square error parameter estimation algorithms for fast-fading channels. The time correlation between successive channel gains is exploited to improve the estimation and tracking capabilities of adaptive algorithms and provide robustness against time-varying channels. Bidirectional normalized least mean-square and conjugate gradient algorithms are devised along with adaptive mixing parameters that adjust to the time-varying channel correlation properties. An analysis of the proposed algorithms is provided along with a discussion of their performance advantages. Simulations for an application to interference suppression in multiuser DS-CDMA systems show the advantages of the proposed algorithms.

Keywords: Multiuser systems; Interference suppression; Adaptive algorithms; Mobile channels

1 Introduction

Low-complexity reception and interference suppression are essential in multiuser mobile systems if battery power is to be conserved, data-rates improved and quality of service enhanced. Conventional adaptive schemes fulfill many of these requirements and have been a significant focus of the research literature [18]. However, in time-varying fading channels commonly associated with mobile systems, these adaptive techniques encounter tracking and convergence problems. Optimum closed-form solutions can address these problems but their computational complexity is high and CSI is required. Low-complexity adaptive channel estimation can provide CSI but in highly dynamic channels tracking problems exist due to their finite adaptation rate [9]. An alternative statistical approach is to obtain the correlation structures required for optimal minimum mean-square error (MMSE) or least-squares (LS) filtering [10, 11]. Although this relieves the tracking demands placed on the filtering process, in a Rayleigh fading channel, a zero correlator is the result due to the expectation of a Rayleigh fading coefficient, and therefore the cross-correlation vector, equating to zero, i.e., E [h1[n] ] = 0 and

E

b1[n] r[n]

= 0. In slowly fading channels, this problem

*Correspondence: mailto: [email protected]

Web End [email protected] 1Department of Electronics, University of York, York Y010 5DD, UK 2CETUC/PUCRIO, 22453-900 Rio de Janeiro, Brazil

RESEARCH Open Access

Bidirectional algorithms for interference suppression in multiuser systems

may be overcome by using a time averaged approach where the averaging period is equal to or less than the coherence time of the channel. However, in fast fading channels, an averaging period equal to the coherence time of the channel is insufficient to overcome the effects of additive noise and characterize the multiuser interference (MUI) [1].

Furthermore, the use of optimized convergence parameters such as step sizes and forgetting factors into conventional adaptive algorithms extend their fading range and lead to improved convergence and tracking performance [8, 1218]. However, the stability of adaptive step-sizes and forgetting factors can be a concern unless they are constrained to lie within a predefined region [19]. Other alternative schemes include those based on processing the received data in subblocks [2022] and subspace algorithms [2328]. In addition, the fundamental problem of obtaining the unfaded symbols whilst suppressing MUI remains. Consequently, the application of such algorithms is restricted to low and moderate fading rates. The limitations of conventional estimation approaches led to the development of methods that attempt to track the faded symbol, such as the channel-compensated MMSE solution [29, 30]. This removes the burden of fading coefficient estimation from the receive filter. However, a secondary process is required to perform explicit estimation of the fading coefficients in order to perform symbol estimation[31].

2015 Clarke and de Lamare. 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/

Web End =(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.

Clarke and de Lamare EURASIP Journal on Wireless Communications and Networking (2015) 2015:228 Page 2 of 17

Approaches that avoid tracking and estimation of the fading coefficients were proposed in [3133]. Although a channel might be highly time variant, two adjacent fading coefficient will be similar and have a significant level of correlation as studied in [3133]. These properties can then be exploited to obtain a sequence of faded symbols where the primary purpose of the filter is to suppress multiuser interference and track the ratio between successive fading coefficients; thus, not burdening it with estimation of the fading coefficients themselves. However, this scheme has a number of limitations stemming from the use of only one correlation time instant and a single class of adaptive algorithms.

In this work, a bidirectional MMSE-based interference suppression scheme for highly dynamic fading channels is presented. The non-zero correlation between multiple time instants is exploited to improve the robustness, tracking, and convergence performance of existing MMSE schemes. Unlike existing adaptive solutions [8, 3133], which do not fully exploit the fading correlation between multiple successive time instants, the proposed bidirectional approach exploits the correlation and adaptively weighs the output of the receive filter in order to optimize the estimation performance. Normalized least-mean square (NLMS) and conjugate gradient (CG)-type algorithms are presented that overcome a number of problems associated with applying the recursive least-squares (RLS) algorithm to bidirectional problems. Novel mixing strategies that weigh the contribution of the considered time instants and improve the convergence and steady-state performance, increasing the robustness against the channel discontinuities, are also presented. An analysis of the proposed schemes is developed and establishes the mechanisms and factors behind their behavior and expected performance. The proposed schemes are applied to conventional multiuser DS-CDMA [2] and cooperative DSCDMA systems [3, 4] to assess their MUI suppression and tracking capabilities. The application of the proposed scheme and algorithms to multiple-antenna and multicarrier systems are also possible. Simulations show that the algorithms improve upon existing schemes with minimal increase in complexity.

The main contributions of this work can be summarized as:

Bidirectional MMSE-based interference suppression scheme for highly dynamic fading channels.

Bidirectional adaptive parameter estimation algorithms based on NLMS and CG techniques.

An analysis of the convergence and the computational complexity of the proposed algorithms.

A study of the proposed and existing algorithms in DS-CDMA and cooperative DS-CDMA multiuser systems.

This paper is organized as follows. Section 2 briefly details the signal models of a conventional DS-CDMA system and a cooperative DS-CDMA system. Section 3 presents the proposed scheme and its corresponding optimization problems and the motivation behind their development. Switching and mixing strategies that optimize performance are proposed and assessed in Section 4, followed by the derivation of the proposed algorithms in Section 5. An analysis of the proposed algorithms is given in Section 6, whereas performance evaluation results are presented in Section 7. Conclusions are drawn in Section 8.

2 Signal models

In this section, we describe the signal models of a DSCDMA system operating in the uplink and a cooperative DS-CDMA system in the uplink equipped with relays and the amplify-and-forward (AF) cooperation protocol. These systems are employed for testing the proposed algorithms even though that extensions to multiple-antenna and multi-carrier can also be considered with appropriate modifications of the algorithms.

2.1 DS-CDMA signal model

We consider the uplink of a synchronous DS-CDMA system with K users, N chips per symbol, and Lp (Lp < N)

propagation paths for each link. We assume that the delay is a multiple of the chip rate, the channel is constant during each symbol interval and the spreading codes are repeated from symbol to symbol. The received signal after filtering by a chip-pulse matched filter and sampled at chip rate yields the M-dimensional received vector given by

r[i] = A1b1[i] H1[i] c1[i] +

K

+[i] +n[i] ,

(1)

where M = N + L 1, and ck[i] and Ak are the spreading sequence and signal amplitude of the kth user, respectively. The M N channel matrix with L paths is given by

Hk[i] for the kth user, the M 1 vector [i] corresponds to the intersymbol interference and n[i] is the noise vector. Conventional schemes use BPSK modulation and the differential and bidirectional schemes employ differential BPSK where the sequence of data symbols to be transmitted by the kth user are given by bk[i] = ak[i] bk[ i 1]

where ak[i] is the unmodulated baseband data. Assuming that linear receive processing is adopted, the output of the receive filter is given by

x[i] = wH[i] r[i] , (2)

k=2Akbk[i] Hk[i] ck[i]

MUI

Clarke and de Lamare EURASIP Journal on Wireless Communications and Networking (2015) 2015:228 Page 3 of 17

where w[i] is an M-dimensional vector that corresponds to the receive filter.

2.2 Cooperative DS-CDMA signal model

We also consider the uplink of a cooperative DS-CDMA system with K users, Nr relays, N chips per symbol, and Lp (Lp < N) propagation paths for each link. The system shown in Fig. 1 is equipped with an AF protocol at each relay. The received signals at the nth relay and the destination nodes are filtered by a chip-pulse matched filter, sampled at a chip rate to obtain sufficient statistics and organized into M 1 vectors as described by

rsrn[i] =

K

k=1ask[i] bk[i] hsrn[i] ck[i] +nrn[i] , (3)

rrd[i] =

Nr

n=1arn[i] hrnd[i] rsrn[i] +nd[i] (4)

and

K

rrd[i] =

Nr

n=1

k=1ask[i] arn[i] hsrn[i] hrnd[i] ck[i] bk[i]

+

Nr

n=1arn[i] hrnd[i] nrn[i] +nd[i] .

(5)

where hsrn[i] and hrnd[i] are the channel fading channel coefficients between the source and the nth relay, and the nth relay and the destination, respectively, and nrn[i]

and nd[i] are additive white Gaussian noise vectors at the relays and the destination, respectively.

The received data is processed by a linear receive filter, which produces the output given by

x[i] = wH[i] rrd[i] , (6)

where w[i] is an M-dimensional vector that corresponds to the receive filter for the cooperative system.

3 Proposed bidirectional scheme

Adaptive parameter estimation has two primary objectives: estimation and tracking of the desired parameters. When applied to multiuser wireless systems, these translate into recovery of the desired symbol, tracking of channel variations and suppression of MUI. However, in fast fading channels, these objectives place unrealistic demands on conventional filtering and estimation schemes. Differential techniques reduce these demands by relieving adaptive receivers from the task of tracking fading coefficients [31]. This is achieved by posing an optimization problem where the ratio between two successive received samples is the quantity to be tracked. Such an approach is enabled by the presumption that, although the fading is fast, there is correlation between the adjacent channel samples as described by

f1[i] = E

h1[i] h1[ i + 1] 0, (7)

where h1[i] is the channel coefficient of the desired user. The interference suppression of the resulting receive filter is improved in fast fading environments compared to conventional adaptive receivers but only the ratio of adjacent fading samples is obtained. Consequently, differential MMSE schemes are suited to differential modulation where the ratio between adjacent symbols is the data carrying mechanism.

However, limiting the optimization to two adjacent samples exposes these processes to the negative effects of uncorrelated samples

E

h1[i] h1[i + 1]

0, (8)

Fig. 1 Cooperative DS-CDMA system model

Clarke and de Lamare EURASIP Journal on Wireless Communications and Networking (2015) 2015:228 Page 4 of 17

but also does not exploit the correlation that may be present between two or more adjacent samples, i.e.,

f2[i] = E

h1[i] h1[i 1] > 0

f3[i] = E

2

h1[ i + 1] h1[i 1] > 0. (9)

In order to address these weaknesses, we propose a bidirectional MSE cost function based on multiple adjacent samples so that the number of channel scenarios under which the differential MMSE performs beneficial adaptation is substantially increased. Termed the bidirectional MMSE, due to the use of multiple time instants, the motivation behind this proposition is illustrated by the plots of fading/channel coefficients in Fig. 2, where J1 represents the 2 sample differential MMSE. There is a low level of correlation present between samples i and i 1, thus any adaptation of the receive filter will bring little benefit. However, the proposed scheme for 3 time instants operates over J1, J2, and J3; therefore, it can exploit the correlation between i + 1 and i 1 and past data. Figure 2 gives an example of a channel where there is a significant level of correlation between samples.

The proposed bidirectional scheme can be expressed in the form of an optimization problem as described by

wo[i] = arg min

w[i] E

b[i] wH[i] r[i 1] b[ i 1] wH[i] r[i]

...

+

b[i]

wH[i] r[i (D 1)]

b[ i (D 1)] wH[i] r[i]

2

+

b[

i 1] wH[i] r[i 2]

b[ i 2] wH[i] r[i 1]

2

...

+

b[

i 1] wH[i] r[i (D 1)]

b[ i (D 1)] wH[i] r[i 1]

2

+

b[

i (D 2)] wH[i] r[i (D 1)]

b[ i (D 1)] wH[i] r[i (d 2)]

2

,

(10)

Fig. 2 Fading channels: an uncorrelated fading channel on the top figure and a correlated fading channel on the bottom figure

Clarke and de Lamare EURASIP Journal on Wireless Communications and Networking (2015) 2015:228 Page 5 of 17

where D denotes the number of considered time instants. Introducing summations into (10) yields a more concise form

wo[i] = arg min

w[i] E

D2

d=0

D1

l=d+1

b[i

d] wH[i] r[i l]

b[i l] wH[i] r[i d]

2

,

(11)

where an output power constraint is required to avoid the trivial all-zero receive filter solution

E

wH[i] r[i]

4 Switching and weighting strategies

The advantages of a bidirectional scheme operating over 3 time or more time instants have been verified in our studies. However, the performance of the scheme may be degraded when received vectors based on uncorrelated fading coefficients are employed in the update of the receive filter. This is particularly evident from the example with an uncorrelated channel illustrated in Fig. 2, where the contribution to the cost function represented by J3 is unlikely to aid the accurate adaptation of w[i].

To avoid this, we introduce a set of switching or mixing parameters that determine the weighting of the D constituent elements of the bidirectional cost function. The proposed generalized bidirectional cost function with weighting factors is described by

wo[i] = arg min

w[i] E

D2

d=0

2

= 1. (12)

Although the existing differential scheme operates over 2 correlated samples, the proposed scheme is able to exploit the additional correlation present between multiple adjacent samples. Moreover, it is also possible to obtain further gain by weighting the correlation between multiple adjacent samples. However, the benefit of using multiple time instant is dependent on the fading rate of the channel and the related correlation of the channel coefficients. We have investigated the use of multiple time instants and it turns out that a scheme which exploits 3 adjacent samples captures most of the performance benefits. In particular, we have tested the proposed bidirectional scheme and algorithms with various values of adjacent time instants (between 4 and 8) and verified that exploiting extra time instants above 3 does not yield significant gains. In fact, the number of time instants is a parameter to be chosen by the designer.

The optimization problem of the proposed scheme for 3 time instants is given by

wo[i] = arg min

w[i] E

b[i] wH[i] r[i 1] b[i 1] wH[i] r[i]

D1

l=d+1 n[i]

b[i

d] wH[i] r[i l]

b[ i l] wH[i] r[i d]

2

,

(15)

where n = d(D 3) + l + 1. However, we again focus on the case with D = 3 in the remainder of this work. With these modifications, the proposed bidirectional MSE cost function with 3-time instant and weighting factors is given by

wo[i]=arg min

w[i] E

1[i]

b[i]

wH[i] r[i 1] b[i 1] wH[i] r[i]

2

(J1)

+ 2[i]

b[i]

,

(16)

where 0 n 3 for n = 1, 2, 3 are the weighting factors.

The determination of the receive vector samples that correspond to the scenarios depicted in Fig. 2 is essential if correct optimization of the is to be achieved. The use of CSI to achieve this would be an effective but impractical solution due to the difficulty in obtaining CSI; consequently, other methods must be sought. In this section, we propose the use of two alternative metrics: the signal power differential after interference suppression between the considered time instants, and the error between the considered time instants.

Firstly, we consider a switching scheme where 13 = [0, 1] are determined at each time instant based on the following post-filtering power differential metrics:

P1[i] =

w[i]H

r[i]

wH[i] r[i 2] b[i 2] wH[i] r[i]

2

(J2)

+ 3[i]

b[i1]wH[i]

r[i2]b[i2] wH[i]r[i1]

2

(J3)

2

(J1)

,

(13)

where J1J3 equate to those of Fig. 2, and the time instants of interest have been altered to avoid the use of future samples. In addition to (13), an output power constraint is again required to avoid an all-zero trivial solution as given by

E

wH[i] r[i]

+

b[i]

wH[i] r[i 2] b[i 2] wH[i] r[i]

2

(J2)

+

b[i1]

wH[i] r[i2]b[i 2] wHr[i1]

2

(J3)

= 1. (14)

In what follows, we describe switching and weighting strategies to optimize the proposed scheme and obtain further performance gain.

2

P2[i] =

w[i]H

r[i]

2

2

w[i]H

w[i]H

r[i 1]

r[i 2]

2

2

(17)

P3[i] =

w[i]H r[i 1]

2

w[i]H r[i 2]

2

.

Clarke and de Lamare EURASIP Journal on Wireless Communications and Networking (2015) 2015:228 Page 6 of 17

If the power difference for each of J13 exceeds a predefined threshold, the corresponding is set to zero;

therefore, removing the corresponding element of the cost function from the adaptation process at that time instant. For highly dynamic channels, one requires an adaptive threshold which is able to track the changes in the system and determine appropriate time instants based on successive samples. Consequently, for each n a threshold, Tn[i], related to a time-averaged, windowed, root-mean-square of the relevant differential power is used. The value of n is then determined in the following manner:

n[i] =

0 if Pn[i] Tn[i]

1 otherwise , (18)

where

Tn[i] =

e1[i] = b[i] wH[i 1] r[i 1] b[i 1] wH[i 1] r[i] e2[i] = b[i] wH[i 1] r[ i 2] b[i 2] wH[i 1] r[i]

e3[i] = b[i 1] wH[i 1] r[i 2] b[i 2] wH[i 1] r[i 1] .(23)

The forgetting factor, 0 e 1, is user defined and, along with normalization by the total error, eT[i],

and

3n=1n[0] = 1, ensures

3n=1n[i] = 1 and a convex combination at each time instant.

5 Adaptive algorithms

In order to devise low-complexity adaptive algorithms based on the proposed bidirectional schemes, we consider the minimization of the cost function given by

C(w[i]) = E

PPnRMS[i] +(1 P)PnRMS[i] , (19)

b[i] wH[i 1] r[i 1] b[i 1] wH[i] r[i]

2

b[i]

wH[i 2] r[i 2] b[i 2] wH[i] r[i]

2

+

b[i

+ 1] wH[i 2] r[i 2] b[i 2] wH[i] r[i 1]

2

PnRMS[i] =

1m 1

i

l=imPn[l]2, (20)

and is a positive user defined constant greater than unity that scales the threshold. The threshold is set with the help of computer experiments in a similar way as the step size of the NLMS algorithm is tuned. The aim is to scale the threshold such that it will be used to inform the algorithm about the relevant differential power which should be used.

Although the current sample corresponding to Jn may bring little benefit in terms of adaptation, this does not indicate that all previous cost function elements corresponding to Jn should be discarded. An alternative approach is to use a set of convex mixing parameters that are not restricted to 1 or 0. This allows each element of the cost function to be more precisely weighted based on its previous and current values. However, the setting of these mixing parameters is once again problematic if they are fixed. Accordingly, an adaptive implementation that can take account of the time-varying channels and previous values which continue to have an impact on the adaptation of the filter is sought. The errors extracted from the cost function (16) are chosen as the metric for developing algorithms. These provide an input to the weighting factor calculation process that is directly related to the cost function of (16). The time-varying mixing factors are given by

n[i] = en[i 1] +(1 e)

eT[i] |en[i] |eT[i] (21)

where

eT[i] = |e1[i] | + |e2[i] | + |e3[i] |, (22)

and the individual error terms are calculated as

subject to E

w[i]H r[i]

2

= 1.

(24)

This cost function then forms the basis of the adaptive algorithms derived in this section. However, to reduce the complexity of the derivations, enforcement of the nonzero constraint is not included and instead enforced in a stochastic manner at each time instant after the adaptation step is complete [31].

5.1 Normalized least-mean square algorithm

We begin with the low-complexity NLMS implementation that employs an instantaneous gradient in a steepest descent framework. Firstly, the instantaneous gradient of(24) is taken with respect to w[i], yielding

w[i]C(w[i]) = b[i 1] r[i]

b[i] wH[i 1] r[i 1]

b[i 1] wH[i] r[i]

H

b[ i 2] r[i]

b[i] wH[i 2] r[i 2]

b[i 2] wH[i] r[i]

H

b[i 2] r[i 1]

b[i 1] wH[i 2] r[i 2] b[i 2]

wH[i] r[i 1]

H(25)

At this point, in order to improve the convergence performance of the NLMS algorithm, the bracketed error terms of (25) are modified by replacing the receive filters with the most recently calculated one, w[i 1]. The resulting gradient expression is given by

Clarke and de Lamare EURASIP Journal on Wireless Communications and Networking (2015) 2015:228 Page 7 of 17

w[i]C(w[i]) = b[i 1] r[i]

b[i] wH[i 1] r[i 1] b[ i 1] wH[i 1] r[i]

H

e1[i]

H

b[i 2] r[i]

b[i] wH[i 1] r[i 2] b[i 2] wH[i 1] r[i]

. (26)

e2[i]

b[i 2] r[i 1]

b[i 1] wH[i 1] r[i 2] b[i 2] wH[i 1] r[i 1]

H

e3[i]

Placing the above gradient expression in the steepest descent update recursion, we obtain

w[i] = w[i 1] +

M[i] |wH[i 1] r[i 1] |

b[i 1] r[i] e1[i] +b[i 2] r[i] e2[i] +b[i 2] r[ i 1] e3[i]

,

(27)

where is the step-size and the normalization factor, M[i], is given by

M[i] = MM[ i 1] +(1 M)rH[i] r[i] , (28) where M is an exponential forgetting factor [31]. The enforcement of the constraint is performed by the denominator of (27) which ensures that the receive filter w[i]

does not tend towards a zero correlator as the adaptation progresses.

The incorporation of the variable switching and mixing factors of Section 4 has the potential to improve the performance of the above algorithm by optimizing the weighting of the error terms of (26). Integration of the factors given by (18) and (21) yields

w[i] = w[i 1] +

M[i] |wH[i 1] r[i 1] |[1[i] b[i 1] r[i] e1[i] +2[i] b[i 2] r[i] e2[i]

+3[i] b[i 2] r[i 1] e3[i]] (29) as the receive filter update equation.

5.2 Least squares algorithm

To achieve faster convergence and increased robustness to fading, we now pursue a LS-based solution. Firstly, the bidirectional cost function of (24) is cast as an LS problem by replacing the expected value with a weighted summation, as described by

CLS(w[i])=

i

l=1 il

.

(32)

The overall correlation structure is then formed from the summation of the preceding expressions, yielding

R[i] = R1[i] + R2[i] + R3[i] (33)

and

t[i] = t1[i] +t2[i] +t3[i] . (34)

where

w[i] = R1[i] t[i] . (35)

Similarly to the NLMS algorithm, performance improvements can be expected if the variable switching and mixing factors, (18) and (21), are incorporated into the correlation expressions. The resulting expressions are

R[i] = 1[i] R1[i] +2[i] R2[i] +3[i] R3[i] (36)

t[i] = 1[i] t1[i] +2[i] t2[i] +3[i] t3[i] (37)

Introducing the above expression into the RLS framework would lead to a low-complexity algorithm with improved convergence and robustness compared to the NLMS of Section 5.1. This requires the integration of (31) with the matrix inversion lemma [9, 34]. However, the derivation requires an expression with a rank-1 update of the form

R[i] = R[i 1] +r[i] rH[i] (38)

b[i] wH[i 1] r[i 1] b[i 1] wH[i] r[i]

2

+

b[i]

wH[i 2] r[i 2] b[i 2] wH[i] r[i]

2

+

b[i1]

wH[i 2] r[i2] b[i2] wH[i1] r[i 1] |2

,

(30)

where is an exponential forgetting factor. Proceeding as with the conventional LS derivation, and modifying the equivalent error terms in a similar manner to as in (26), we arrive at the following expressions for the component autocorrelation matrices:

R1[i] = R1[i 1] +b[i 1] r[i] rH[i] b[i 1] R2[i] = R2[i 1] +b[i 2] r[i] rH[i] b[i 2]

R3[i] = R3[i 1] +b[i 2] r[i 1] rH[i 1] b[i 2](31)

and the component cross-correlation vectors:

t1[i] = t3[i 1] +b[i 1] r[i] rH[i 1] w[i 1] b[i]

t2[i] = t2[i 1] +b[i 2] r[i] rH[i 2] w[i 1] b[i]

t3[i] = t3[i 1] +b[i 2] r[i 1] rH[i 2] w[i 1] b[i 1]

Clarke and de Lamare EURASIP Journal on Wireless Communications and Networking (2015) 2015:228 Page 8 of 17

for the autocorrelation matrix; a form which (31) is unable to fit into without assumptions that cause a significant performance degradation. Consequently, an alternative low-complexity algorithm to implement the LS solution given by (31)(35) is required.

5.3 Conjugate gradient algorithm

Due to the particular form of the bidirectional LS formulation and the conventional RLS recursion, an alternative low-complexity method is now derived. The CG technique has been chosen to avoid matrix inversions and due to its excellent convergence properties [3537]. We begin the derivation of the proposed CG-type algorithm with the autocorrelation (31) and cross-correlation (32) structures of Subsection 5.2. Inserting them into the standard CG quadratic form yields

J(w) = wH[i] R[i] w[i] tH[i] w[i] . (39)

From [36], the unique minimiser of (39) is also the minimiser of

R[i] w[i] = t[i] . (40)

This shows the suitability of the CG algorithm to the bidirectional problem. At each time instant, a number of iterations of the following method are required to reach an accurate solution, where the iterations are indexed with the variable j. Other single iteration CG methods are available but these depend upon degeneracya term that describes the situation where the successive CG vectors are not orthogonal [38]. Consequently, we employ the conventional method to ensure satisfactory convergence. At the ith time instant, the gradient and direction vectors are initialized as

g0[i] = w[i]CLS(w[i]) = R[i] w0[i] t[i] (41)

and

d0[i] = g0[i] , (42)

respectively, where the gradient expression is equivalent to those used in the derivation of the previous algorithms. The vectors dj[i] and dj+1[i] are R[i] orthogonal with respect to R[i] such that dj[i] R[i] dl[i] = 0 for j = l. At each iteration, the receive filter is updated as

wj+1[i] = wj[i] +j[i] dj[i] (43)

where j[i] is the minimizer of J(wj+1[i]) such that

j =

dHjgj[i]

dHj[i] R[i] dj[i]

and a new conjugate gradient direction vector is obtained as given by

dj+1[i] = gj+1[i] +j[i] dj[i] (46)

where

j[i] =

gHj+1[i] R[i] dj[i]

dHj[i] R[i] dj[i]

(47)

ensures the R[i] orthogonality between dj[i] and dl[i] where j = l. The iterations (43)(47) are then repeated until j = jmax.

The variable switching and mixing factors can be incorporated into the algorithm to improve performance. This is achieved by operating the CG algorithm over the modified correlation structures given by (36) and (37).

6 Analysis

In this section, we analyze the proposed bidirectional algorithms to gain insight of the expected performance but also to obtain further knowledge into the operation of the proposed and existing algorithms. The unconventional form of the proposed cost functions precludes the application of standard MSE analysis. Consequently, we concentrate on the signal-to-interference-plus-noise ratio (SINR) of the proposed algorithms in order to analyze their interference suppression and tracking performance. We firstly study the NLMS algorithm and the features of its weight error correlation matrix in order to arrive at an analytical SINR expression. Following this, we explore the analogy between the form of the bidirectional expression and convex combinations of adaptive receive filters [39, 40].

6.1 SINR analysis

Let us define the instantaneous SINR expression given by

SINRinst [defines] wH[i] RSw[i]wH[i] RIw[i] , (48)

where RS and RI are the signal and interference and noise correlation matrices, into a form amenable to analysis.

The receive filter error weight vector is described by

[i] = w[i] wo[i] , (49)

where wo is the instantaneous standard MMSE receiver.

Let us now describe the expression of the numerator of(48) with the desired signal component:

Sc = H[i] RS[i] + H[i] RSwo[i] +

PS,opt[i]

wHo[i] RSwo[i]

+ wHo[i] RS[i] ,

. (44)

The gradient vector is then updated according to

gj+1[i] = R[i] wj[i] t[i] (45)

(50)

Clarke and de Lamare EURASIP Journal on Wireless Communications and Networking (2015) 2015:228 Page 9 of 17

and of the interference plus noise component described by

Si+n = H[i] RI[i] + H[i] RIwo[i] + wHo[i] RIwo[i]

PI,opt[i]

where F and R are correlation matrices. Specifically, R13 are autocorrelation matrices given by

R1 = E

r[i] b[i 1] rH[i] b[i 1] R2 = E

r[i] b[i 2] rH[i] b[i 2] R3 = E

r[i 1] b[i 2] rH[i 1] b[i 1]

(55)

and F13 cross-time instant correlation matrices, given by

F1 = E

r[i] b[i 1] rH[i 1] b[i] F2 = E

r[i] b[i 2] rH[i 2] b[i] F3 = E

r[i 1] b[i 2] rH[i 2] b[i 1] .

(56)

Using (54) and the independence assumption of E

eo,13[i] [i]

= 0, E

rH[i] r[ i 1] = 0 and E [bk[i] bk[i 1] ] = 0, we arrive at the expression for K[i]

K[i] = [I + F1 R1 + F2 R2 + F3 R3]K[ i 1] [I + F1 R1 + F2 R2

+F3 R3]

+2R1Jmin,1[i]

+2R2Jmin,2[i]

+2R1Jmin,3[i]

(57)

where Jmin,j[i] = |eo,j|2. Following a similar method, an expression for G[i] can also be reached

G[i] = G[i 1] [F1 R1 + F2 R2 + F3 R3] .(58)

At this point, we study the derived expression to gain an insight into the operation of the bidirectional algorithm and the origins of its advantages over the conventional differential scheme. Equivalent expressions for the conventional stochastic gradient scheme are given by

K[i] = [I + F1 R1] K[i 1] [I + F1 R1]+2R1Jmin,1[i]

G[i] = G[i 1] [F1 R1] .(59)

The bidirectional scheme has a number of additional correlation terms compared to the conventional scheme. Evaluating the cross-time instant matrices yields

F1 = |a1|2c1cH1 E h[i] h[i 1]

+ wHo[i] RI[i] .

(51)

Taking the expectation and the trace of Sc and Si+n, defining K[i] = E[ [i] H[i] ] and G[i] = E[wo[i] H[i]], we can define the following SINR expression:

SINR [defines] tr{E[ Sc] } tr{E[ Si+n] }

=

tr[ K[i] RS + G[i] RS + PS,opt[i] +GH[i] RS] tr[ K[i] RI + G[i] RI + PI,opt[i] +GH[i] RI]

.

(52)

From (52), it is clear that we need to pursue expressions for K[i] and G[i] in order to reach an analytical interpretation of the bidirectional NLMS scheme.

Substituting the filter error weight vector into the filter update expression of (27) yields a recursive expression for the receive filter error weight vector described by

[i] =[i 1]

+ I + r[i] b[i 1] rH[i 1] b[i] r[i] b[i 1] rH[n] b[i 1]

+ r[i] b[i 2] rH[i 2] b[i] r[i] b[i 2] rH[n] b[i 2]

+ r[i 1] b[i 2] rH[i 2] b[i 1]

r[i 1] b[i 2] rH[i 1] b[i 2]

[i 1]

+ r[i] b[i 1] eo,1[i]

+ r[i] b[i 2] eo,2[i]

+ r[i 1] b[i 2] eo,3[i]

,

(53)

where the terms eo,13 denote the error terms of (26) when the optimum filter wo is used. Utilizing the direct averaging approach developed by Kushner [41], often invoked in this type of stochastic analysis, the solution to the stochastic difference equation of (53) can be approximated by the solution to a second equation [9, 42], such that

E

I + r[i] b[i 1] rH[i 1] b[i] r[i] b[i 1] rH[n] b[i 1] +r[i] b[i 2] rH[i 2] b[i] r[i] b[i 2]

rH[n] b[i 2] +r[i 1] b[i 2] rH[i 2] b[i 1]

r[i 1] b[i 2] rH[i 1] b[i 2]

= I + F1 R1 + F2 R2 + F3 R3

,

(54)

h[i] h[i 2]

f1[i]

F2 = |a1|2c1cH1 E

. (60)

h[i 1] h[i 2]

f2[i]

F3 = |a1|2c1cH1 E

f3[i]

Clarke and de Lamare EURASIP Journal on Wireless Communications and Networking (2015) 2015:228 Page 10 of 17

From the expression above, it is clear that the underlying factor that governs the SINR performance of the algorithms is the correlation between the considered time instants, f13, data-ruse and the use of f2. Accordingly, it is the additional correlation factors that the proposed bidirectional algorithms possess that enhances its performance compared to the conventional scheme, confirming the initial motivation behind the proposition of the bidirectional approach. Lastly, the f13 expressions of (60)

are the factors that influence the optimum number of considered time instants.

6.2 Combinations of adaptive receive filters

To further our understanding of the bidirectional algorithms, we follow a heuristic and complementary approach that leads to an analogy with a combination of adaptive filters [39]. The bidirectional LS solution given by (35) is made up of 6 constituent correlation structures that result in a filter output of

y[i] =

(1R1[i] +3R2[i] +3R3[i] )1

(1t1[i] +2t2[i] +3t3[i])

H r[i] .

(61)

Decomposing the expression above leads us to an expression where the signal y[i] is formed from the output of 3 individual adaptive receive filters

y[i] =

R1[i] +21 R2[i] +31 R3[i] 1

t1[i]

H r[i]

+ R1[i] +12 R2[i] +32 R3[i] 1

.

(62)

This is equivalent to a convex combination of adaptive receive filters with varying [39, 40], where each of the 3 filters focuses on the correlation between the 2 of the 3 considered time instants. However, the presence of the autocorrelation matrices in the inverses of the expression also indicates that the remaining time instants also influence the structure of each filter. Although the mixing factors are not separable, we can interpret them as a form of weighting that is present in conventional combinations of adaptive filters. This explains in part the additional control and performance they provide.

7 Simulations

In this section, the proposed bidirectional adaptive algorithms are applied to conventional multiuser and cooperative DS-CDMA systems using the signal models described in Section 2. The application of the proposed algorithms

to multiple-antenna and multi-carrier systems is straightforward and requires a change in the signal models. The individual Rayleigh fading channel coefficients, h[i], are generated using Clarkes model [43] where 20 scatterers are assumed. In all simulations, the number of packets is denoted by Np and the fading rate is given by the dimensionless normalized fading parameter, Tsfd, where Ts is the symbol period and fd is the Doppler frequency shift. The convergence parameters of the algorithms have been optimized resulting in step-sizes forgetting factors of 0.1 and0.99, respectively, e = 0.95, M = 0.99, and the number of CG iterations, jmax = 5.

As detailed in Section 6, the proposed algorithms do not minimize the same MSE as a conventional MMSE receiver; therefore, the MSE is not an adequate performance metric. As a result, BER- and SINR-based metrics are chosen for the purpose of comparison between existing algorithms and the optimum MMSE solution. Due to the rapidly fading channel, the instantaneous SNR, SNRi,

is highly variable and so the SINR alone is also not a satisfactory metric. To overcome this, it is normalized by the instantaneous SNR to give SINRSNRi. This value is negative in all simulations and directly reflects the MUI interference suppression and tracking capabilities of the proposed algorithms [31, 32].

7.1 Conventional DS-CDMA

Here we apply the adaptive algorithms of Section 5 to interference suppression in the uplink of a multiuser DSCDMA system described in Section 2. Each simulation is averaged over Np packets and detailed parameters are specified in each plot.

7.1.1 Analytical results

We first assess the analytical expressions derived in Section 6.1 and their agreement with simulated results. Central to the performance of the differential and bidirectional schemes are the correlation factors f13 and the related assumption of h1[i] h1[i 1]. Examining the effect of the fading rate on the value of f13 shows that f1 f2 f3 at fading rates of up to Tsfd = 0.01. Consequently, after a large number of received symbols with high total receive power

3 [I + F1 R1] [I + F1 R1 + F2 R2

+F3 R3] ,

(63)

due to the decreasing significance of the identity matrix. This indicates that the expected value of the SINR, of the bidirectional scheme, once f1 f2 f3 have stabilized, should be similar to the differential scheme. A second implication is that the bidirectional scheme should converge towards the MMSE level due to the equivalence between the differential scheme and the MMSE solution

t2[i]

H r[i]

+ R1[i] +13 R2[i] +23 R3[i] 1

t3[i]

H r[i]

Clarke and de Lamare EURASIP Journal on Wireless Communications and Networking (2015) 2015:228 Page 11 of 17

[31]. Figure 3 illustrates the analytical performance using the expressions given in Section 6.1.

The correlation matrices are calculated via ensemble averages prior to the start of the algorithm and G[ 0] =

K[0] = I. In Fig. 3, one can see the convergence of the simulated schemes to the analytical and MMSE plots, validating the presented analysis. Due to the highly dynamic nature of the channel, using the expected values of the correlation matrix alone cannot capture the true transient performance of the algorithms. However, the convergence period of the analytical plots within the first 200 iterations can be considered to be within the coherence time and therefore give an indication of the transient performance relative to other analytical plots. Using this justification and the aforementioned analysis, advantages should be present in the transient phase due to the additional correlation information supplied by F2 and F3. This conclusion is supported by Fig. 3 and the similar forms of the analytical and simulated schemes relative to each other and their subsequent convergence.

7.1.2 SINR performance

The SINR/SNR performance of the proposed algorithms is given by Figs. 4 and 5. The performance of the CG implementation of the differential algorithm is marginally below that of the RLS during convergence but the bidirectional scheme provides noticeable improvements. The differential and bidirectional algorithms converge close to the MMSE optimum as expected from the previous analysis. The bidirectional NLMS algorithm provides more significant improvements over the differential scheme,

both in the final stages of convergence and steady-state. These differences can be accounted for by the reduced receive signal power; the matrices equivalent to F2 and F3 improving the consistency of the steady-state performance by reducing the impact of weakly correlated samples; and the NLMSs suitability to data reuse as in the affine projection (AP) algorithm. As expected, the conventional adaptive schemes are unable to converge or track the solution due to the more demanding task of tracking both the fading coefficients and suppressing MUI.

The BER performance of the differential and bidirectional schemes is illustrated in Fig. 6, where the system parameters are equal to those of Figs. 4 and 5. The RLS and CG algorithms converge to near the MMSE level with the bidirectional scheme providing a performance improvement over other considered algorithms. The NLMS schemes exhibit slower BER convergence compared to their SINR performance but reach a level where a decision-directed operation can take place in a severely fading channel. Due to the superior performance of the CG- and RLS-based algorithms, we predominantly focus on their performance for the remainder of this section.

Figure 7 illustrates the performance of the proposed CG and RLS algorithms as the fading rate is increased. The conventional schemes with RLS and CG algorithms are unable to cope with fading rates in excess of Tsfd =

0.005 and begin to diverge at the completion of the training sequence. The proposed bidirectional scheme outperforms the differential schemes but the performance begins to decline once fading rates above fdTs = 0.01 are reached.

Fig. 3 SINR performance comparison of simulated and analytical proposed NLMS algorithms over a single path channel

Clarke and de Lamare EURASIP Journal on Wireless Communications and Networking (2015) 2015:228 Page 12 of 17

Fig. 4 SINR performance comparison of proposed CG algorithms over a single path channel where all schemes have been trained with 150 symbols and then switched to decision-directed mode

Once again, the increase in performance of the bidirectional scheme can be accounted for by the increased correlation information supplied by the matrices F2 and F3 and data reuse. The introduction of the mixing factors into the bidirectional algorithm improves performance further, especially at higher fading rates. A first reason for

this is the improvement in consistency as previously mentioned. However, a second more significant reason can be established by referring back to the observations on the correlation factors f13. Although fading rates of 0.01 may be fast, the assumption h[i 2] h[i 1] h[i] is still valid. Consequently, f1 f2 f3 and equal weighting

Fig. 5 SINR performance comparison of proposed NLMS algorithms over a single path channel where all schemes have been trained with 150 symbols and then switched to decision-directed mode

Clarke and de Lamare EURASIP Journal on Wireless Communications and Networking (2015) 2015:228 Page 13 of 17

Fig. 6 BER performance comparison of proposed schemes during training over a single path channel

is adequate. However, as the fading rate increases beyond Tsfd = 0.01, this assumption breaks down and the correlation information requires unequal weighting for optimum performance, a task fulfilled by the adaptive mixing factor.

A more detailed plot illustrating the performance advantages of the CG switching and mixing parameters presented in is given by Fig. 8. The switching approach

provides little improvement over the standard bidirectional scheme due to its discrete and non-adaptive operation. As previously covered, a low instantaneous value of f13, as indicated by a large power differential, does not indicate that all information gathered on f13 is redundant. The mixing parameter implementations address this shortcoming by adaptively setting the parameters via the

Fig. 7 SINR performance versus fading rate of the proposed CG schemes over a single path channel after 200 training symbols

Clarke and de Lamare EURASIP Journal on Wireless Communications and Networking (2015) 2015:228 Page 14 of 17

Fig. 8 SINR performance over a single path channel of the proposed schemes with switching and mixing factors

error weight expression (21) that accurately reflects the averaged correlation factors. At a fading rate of Tsfd =

0.02 the assumption of f1 f2 f3 begins to diminish in accuracy and therefore unequal weighting is required for performance in excess of the standard bidirectional scheme, as previously mentioned and shown in Fig. 8.

The MUI suppression of the proposed and existing schemes is given by Fig. 9. The bidirectional scheme has significantly improved multiuser performance compared to the differential algorithms at low system loads but diminishes as the number of users increases. This behavior supports the analytical conclusions of Section 6.1

Fig. 9 BER performance against system loading after 500 symbols of the proposed schemes over a single-path channel. Schemes are trained with 150 symbols and then switch to decision directed operation

Clarke and de Lamare EURASIP Journal on Wireless Communications and Networking (2015) 2015:228 Page 15 of 17

Fig. 10 SINR performance of the proposed CG schemes during training in a single-path cooperative DS-CDMA system

by virtue of the convergence of the differential and bidirectional schemes and the increasing accuracy of(63) as system loading, and therefore received power, increases.

7.2 Cooperative DS-CDMA

To further demonstrate the performance of the proposed schemes in cooperative relaying systems [5], we apply

them to an AF cooperative DS-CDMA system detailed in Section 2.

Figure 10 shows that the bidirectional scheme obtains performance benefits over the differential schemes during convergence but, as expected, the performance gap closes as steady-sate is reached. The inclusion of variable mixing parameters improves performance but to a lesser extent than non-cooperative networks due to the more

Fig. 11 BER performance of the proposed CG schemes in a single-path cooperative DS-CDMA system

Clarke and de Lamare EURASIP Journal on Wireless Communications and Networking (2015) 2015:228 Page 16 of 17

challenging scenario of compounding highly time-variant channels.

The improvement BER brought about by the bidirectional schemes is evident from Fig. 11. However, the more challenging environment of a cooperative system with compounded rapid fading has impacted on the BER performance of the schemes, as evidenced by the increased performance gap between the proposed schemes and MMSE reception.

8 Conclusions

In this paper, we have presented a bidirectional MMSE framework that exploits the correlation characteristics of rapidly varying fading channels to overcome the problems associated with conventional adaptive interference suppression techniques in such channels.

An analysis of the proposed schemes has been performed and the reasons behind the performance improvements shown to be the additional correlation information, data reuse, and optimized correlation factor weighting. The conditions under which the differential and bidirectional schemes are equivalent have also been established and the steady-state implications of this detailed. Finally, the proposed algorithms have been assessed in standard and cooperative multiuser DS-CDMA systems and shown to outperform both differential and conventional schemes.

Competing interests

The authors declare that they have no competing interests.

Acknowledgements

Part of this manuscript was presented at the International Symposium on Wireless Communications Systems (ISWCS) in 2013. The work of R. C. de Lamare is partly funded by CNPq and FAPERJ.

Received: 2 March 2015 Accepted: 31 August 2015

References

1. S Verdu, Multiuser Detection. (Cambridge University Press, NY, 1998)2. U Madhow, M Honig, MMSE interference suppression for direct-sequence spread spectrum CDMA. IEEE Trans. Commun. 42(12), 31783188 (1994)

3. W Huang, Y Wong, C Kuo, in IEEE Global Communications Conference. Decode-and-forward cooperative relay with multi-user detection in uplink CDMA networks (Phoenix, USA, 2007)

4. RC de Lamare, Joint iterative power allocation and linear interference suppression algorithms for cooperative ds-cdma networks. IET Commun. 6(13), 19301942 (2012)

5. P Clarke, RC de Lamare, Transmit diversity and relay selection algorithms for multirelay cooperative mimo systems. IEEE Trans. Veh. Technol. 61(3), 10841098 (2012)

6. T Wang, RC de Lamare, PD Mitchell, Low-complexity set-membership channel estimation for cooperative wireless sensor networks. IEEE Trans. Veh. Technol. 60(6), 25942607 (2011)

7. RC de Lamare, R Sampaio-Neto, Minimum mean squared error iterative successive parallel arbitrated decision feedback detectors for DS-CDMA, systems. IEEE Trans. Commun. 5(56), 778789 (2008)

8. R de Lamare, P Diniz, Set-membership adaptive algorithms based on time-varying error bounds for cdma interference suppression. IEEE Trans. Veh. Technol. 58(2), 644654 (2009)

9. S Haykin, Adaptive Filter Theory, 4th edn. (Prentice Hall, NJ, 2002)

10. D Sadler, A Manikas, MMSE multiuser detection for array multicarrier DS-CDMA in fading channels. IEEE Trans. Signal Process. 53(7), 23482358 (2005)

11. R Schoder, W Gerstacker, A Lampe, Noncoherent MMSE interference suppression for DS-CDMA. IEEE Trans. Commun. 50(4), 577587 (2002)

12. S Haykin, A Sayed, J Zeidler, P Yee, P Wei, Adaptive tracking of linear time variant systems by extended RLS algorithms. IEEE Trans. Signal Process. 45(5), 11181128 (1997)

13. J Wang, Fast tracking RLS algorithm using novel variable forgetting factor with unity zone. IEEE Electron. Lett. 27(23), 25502551 (1991)

14. R Kwong, E Johnston, A variable step size LMS algorithm. IEEE Trans. Signal Process. 40(7), 16331642 (1992)

15. B Toplis, S Pasupathy, Tracking improvements in fast RLS algorithm using a variable forgetting factor. IEEE Trans. Acoust., Speech, Signal Process. 36(2), 206227 (1988)

16. S-H Leung, C So, Gradient-based variable forgetting factor RLS algorithm in time-varying environments. IEEE Trans. Signal Process. 53(8), 31413150 (2005)

17. S Hyun-Chool, A Sayed, S Woo-Jin, Variable step-size NLMS and affine projection algorithms. IEEE Signal Process. Lett. 11(2), 132135 (2004)

18. Y Zhang, N Li, J Chambers, Y Hao, New gradient-based variable step size LMS algorithms. EURASIP J. Adv. Signal Process, 19 (2008)

19. S Gelfand, Y Wei, J Krogmeier, The stability of variable step-size LMS algorithms. IEEE Trans. Signal Process. 47(12), 32773288 (1999)

20. P Baracca, S Tomasin, L Vangelista, N Benvenuto, A Morello, Per sub-block equalization of very long ofdm blocks in mobile communications. IEEE Trans. Commun. 59(2), 363368 (2011)

21. S Yerramalli, M Stojanovic, U Mitra, Partial fft demodulation: A detection method for highly doppler distorted ofdm systems. IEEE Trans. Signal Process. 60(11), 59065918 (2012)

22. L Li, A Burr, R de Lamare, in IEEE 77th Vehicular Technology Conference (VTC Spring) 2013. Joint iterative receiver design and multi-segmental channel estimation for ofdm systems over rapidly time-varying channels (Dresden, Germany, 2013), pp. 15

23. ML Honig, JS Goldstein, Adaptive reduced-rank interference suppression based on the multistage wiener filter. IEEE Trans. Commun. 6(50), 986994 (2002)

24. Y Sun, V Tripathi, ML Honig, Adaptive, iterative, reducedrank (turbo) equalization. IEEE Trans. Wirel. Commun. 4(6), 27892800 (2005)

25. RC de Lamare, R Sampaio-Neto, Reducedrank adaptive filtering based on joint iterative optimization of adaptive filters. IEEE Signal Process. Lett. 14(12), 980983 (2007)

26. RC de Lamare, R Sampaio-Neto, Reduced-rank spacetime adaptive interference suppression with joint iterative least squares algorithms for spread-spectrum systems. IEEE Trans. Veh. Technol. 59(3), 12171228 (2010)

27. RC de Lamare, R Sampaio-Neto, Adaptive reduced-rank equalization algorithms based on alternating optimization design techniques for MIMO systems. IEEE Trans. Veh. Technol. 60(6), 24822494 (2011)

28. RC de Lamare, R Sampaio-Neto, Adaptive reduced-rank processing based on joint and iterative interpolation, decimation, and filtering. IEEE Trans. Signal Process. 57(7), 25032514 (2009)

29. M Honig, M Shensa, S Miller, L Milstein, in IEEE Vehicular Technology Conf. Performance of adaptive linear interference suppression for DS-CDMA in the presence of flat rayleigh fading (Phoenix, 1997)

30. HV Poor, X Wang, in Annual Allerton Conf. Commun., Control and Computing. Adaptive multiuser detection in fading channels (Monticello, 1996)

31. U Madhow, K Bruvold, LJ Zhu, Differential MMSE: A framework for robust adaptive interference suppression for DS-CDMA over fading channels. IEEE Trans. Commun. 53(8), 13771390 (2005)

32. LJ Zhu, U Madhow, in IEEE Global Telecoms. Conf. Adaptive interference suppression for direct sequence CDMA over severely time-varying channels (Pheonix, 1997)

33. M Honig, S Miller, M Shensa, L Milstein, Performance of adaptive linear interference suppression in the presence of dynamic fading. IEEE Trans. Commun. 49(4), 635645 (2001)

34. PSR Diniz, Adaptive Filtering: Algorithms and Practical Implementation, 3rd edn. (Springer, 2008)

35. C Meyer, Matrix Analysis and Applied Linear Algebra. (SIAM, 2001)

Clarke and de Lamare EURASIP Journal on Wireless Communications and Networking (2015) 2015:228 Page 17 of 17

36. DG Luenburger, Y Ye, Linear and Nonlinear Programming, 3rd edn. (Springer, 2008)

37. J Lee, H Yu, Y Sung, Beam tracking for interference alignment in time-varying mimo interference channels: A conjugate-gradient-based approach. IEEE Trans. Veh. Technol. 63(2), 958964 (2014)

38. PS Chang, AN Wilson, Analysis of conjugate gradient algorithms for adaptive filtering. IEEE Trans. Signal Process. 48(2), 409418 (2000)

39. NJ Bershad, JCM Bermudez, J-Y Yourneret, An affine combination of two LMS adaptive filters - transient mean square analysis. IEEE Trans. Signal Process. 56(8), 18531864 (2008)

40. J Arenas-Garcia, A Figueiras-Vidal, A Sayed, Mean-square performance of a convex combination of two adaptive filters. IEEE Trans. Signal Process. 45(3), 10781090 (2006)

41. H Kushner, Approximation and Weak Convergence Methods for Random Processes with Applications to Stochastic System Theory. (MIT Press, 1984)

42. A Sayed, Adaptive Filters. (John Wiley & Sons, 2008)43. W Jakes, Microwave Mobile Communications. (Wiley-IEE Press, 1994)

Submit your manuscript to a journal and benet from:

7 Convenient online submission7 Rigorous peer review7 Immediate publication on acceptance7 Open access: articles freely available online 7 High visibility within the eld7 Retaining the copyright to your article

Submit your next manuscript at 7 springeropen.com