Time-Varying FIR Decision Feedback Equalization of ... - CiteSeerX

11 downloads 0 Views 139KB Size Report
+ HHR−1. ˜CH). −1 . (25). The last line in (25) is obtained by applying the matrix inversion. Lemma. By defining the selection matrix P∆ as: P∆ = ⎡. ⎣0kxv.
Time-Varying FIR Decision Feedback Equalization of Doubly-Selective Channels ∗ Imad Barhumi† , Geert Leus‡ and Marc Moonen K.U.Leuven-ESAT, Kasteelpark Arenberg 10, 3001 Leuven, Belgium Email: {imad.barhumi,geert.leus,marc.moonen}@esat.kuleuven.ac.be

Abstract— In this paper we propose a minimum mean-square error (MMSE) decision feedback (DF) time-varying (TV) finite impulse response (FIR) equalizer for doubly-selective (time- and frequency-selective) channels. We use the basis expansion model (BEM) to approximate the doubly-selective channel and to design the TV FIR DFE. This allows us to turn a large design problem into an equivalent small design problem, containing only the BEM coefficients of the doubly-selective channel and the BEM coefficients of the TV FIR feedforward and feedback equalization filters. Through computer simulations we show that the performance of the proposed TV FIR DF equalizer approaches the performance of the block DF equalizer, while the equalization as well as the design complexity are much lower.

I. I NTRODUCTION

T

HE wireless communication industry has experienced rapid growth in recent years, and digital cellular systems are currently designed to provide high data rates at high speeds. High data rates give rise to intersymbol interference (ISI) due to multipath fading. Such an ISI channel is called frequencyselective. On the other hand, due to mobility and/or frequency offset the received signal is subject to frequency shifts (Doppler shifts). The Doppler effect in conjunction with ISI gives rise to a so-called doubly-selective channel (frequency- and timeselective). Equalization is required to compensate for these channel effects. Equalization can be of different types. Well-known types are i) linear equalization (LE), ii) decision feedback equalization (DFE) iii) maximum likelihood sequence estimation (MLSE). LE is used because of its simplicity and cost effectiveness, but may also sacrifice performance. MLSE is extremely complex, and gives optimum or near optimum performance. DFE on the other hand, provides a good compromise between complexity and performance. The basic idea behind DFE is that once information has been detected, the ISI that it induces on future symbols can be subtracted out before detection of subsequent symbols. DFE has been previously proposed in literature for the case of frequency-selective channels, some based on block processing ∗ This research work was carried out at the ESAT laboratory of the Katholieke Universiteit Leuven, in the frame of the Belgian State, Prime Minister’s Office Federal Office for Scientific, Technical and Cultural Affairs - Interuniversity Poles of Attraction Programme (2002-2007) - IUAP P5/22 (‘Dynamical Systems and Control: Computation, Identification and Modeling’) and P5/11 (‘Mobile multimedia communication systems and networks’), the Concerted Research Action GOA-MEFISTO-666 (Mathematical Engineering for Information and Communication Systems Technology) of the Flemish Government, Research Project FWO nr.G.0196.02 (‘Design of efficient communication techniques for wireless timedispersive multi-user MIMO systems’) and was partially sponsored by IMEC (Flemish Interuniversity Microelectronics Center). The scientific responsibility is assumed by its authors. † Partly supported by the Palestinian European Academic Cooperation in Education (PEACE). ‡ Postdoctoral Fellow of the F.W.O. Vlaanderen.

[1], [2] and others based on serial processing by means of finite impulse response (FIR) filters [3]. In previous work we proposed time-varying (TV) FIR LEs to compensate for the doubly-selective channel effects [4], [5], [6]. In this paper we extend this idea to design a TV FIR DFE, which shows to give better performance and to be more robust against channel imperfections. The proposed TV FIR DFE consists of a TV FIR feedforward filter, a TV FIR feedback filter, and a decision device. We use the basis expansion model (BEM) to model the TV FIR feedforward and feedback filters. The beauty of using the BEM is that it allows us to turn a large design problem into a small design problem, containing only the BEM coefficients of the channel, and the BEM coefficients of the feedforward and feedback TV FIR filters. We also show that this approach extends the results obtained for the case of purely frequency-selective channels. This paper is organized as follows. The system model is described in section II. The block decision feedback equalizer (BDFE) is introduced in section III. The TV FIR DFE equalizer is introduced in section IV. In section V, we show through computer simulations the performance of the proposed equalizer. Finally, our conclusions are drawn in section VI. Notations: We use upper (lower) bold face letters to denote matrices (column vectors). Superscripts ∗ , T , and H represent conjugate, transpose, and Hermitian, respectively. We denote the 1- and 2-dimensional Kronecker delta as δn and δn,m , respectively. We denote the N × N identity matrix as IN and the M × N all-zero matrix as 0M ×N . Finally, diag{x} denotes the diagonal matrix with x on the diagonal. II. S YSTEM M ODEL The system under consideration is depicted in Figure 1. We assume a single-input multiple-output (SIMO) system, where Nr receive antennas are used. Focusing on a baseband-equivalent description, when transmitting a symbol sequence s(n) at rate 1/T and sampling each receive antenna at the same rate 1/T , the received sample sequence at the rth receive antenna can be written as: ∞ X h(r) (n; ν)s(n − ν) + η (r) (n) (1) y (r) (n) = ν=−∞

where η (r) (n) is the additive noise at the rth receive antenna, and h(r) (n; ν) is the doubly-selective channel (time- and frequencyselective) from the transmitter to the rth receive antenna. Here, we model the doubly-selective channel h(r) (n; ν) using the basis expansion model (BEM) [7], [8], [9], [10], which is shown to accurately approximate the well-known Jakes’ model. In the BEM, the channel is modeled as an FIR filter where the

written as:

η (1)(n) y (1)(n)

(1)

h (n; ν)

y(1)

S/P

y(r) = H(r) Tzp s + η (r) ,

(1)

G

N ×1

+

s(n)

η (Nr )(n)

+

y (Nr )(n)

h(Nr )(n; ν)

y(Nr )

S/P

N ×1

˜s

ˆs



G(Nr )

(r)

B

H(r) = (a) Block DFE

y (1)(n)

g (1)(n; ν) +

s(n)

η (Nr )(n) h(Nr )(n; ν)

+

y (Nr )(n)

Q/2 L X X

(r)

hq,l Dq Zl

(5)

l=0 q=−Q/2

η (1)(n) h(1)(n; ν)

(4)

where η is similarly defined as y , Tzp := [IM , 0M ×Lzp ]T , and H(r) is an N × N lower triangular matrix. Using (2), H(r) can be written as: (r)

s˜(n − d)

sˆ(n − d)

where Dq := diag{[1, · · · , ej2πq(N −1)/N ]T }, and Zl is an N × N lower triangular Toeplitz matrix with first column [01×l , 1, 01×(N −l−1) ]T . Substituting (5) in (4), the N ×1 received sample block at the rth receive antenna can be written as:



y(r) =

g (Nr )(n; ν)

Q/2 L X X

(r)

hq,l Dq Zl Tzp s + η (r) .

(6)

l=0 q=−Q/2

b(n; ν) − δd

Stacking the Nr received [y(1)T , . . . , y(Nr )T ]T , we obtain

(b) TV FIR DFE

sample

blocks:

y

:=

H

:=

y = HTzp s + η, Fig. 1.

System Model

taps are expressed as a superposition of complex exponential basis functions. Before we proceed to define the BEM for doublyselective channels, we make the following assumptions: A1) The delay-spread is bounded by τmax ; A2) The Doppler-spread is bounded by fmax . Under assumptions A1) and A2), it is always possible to model the channel between the transmitter and the rth receive antenna in discrete time as: h(r) (n; ν) =

L X l=0

δν−l

Q/2 X

(r)

hq,l ej2πqn/N ,

where η is similarly [H(1)T , . . . , H(Nr )T ]T .

defined

as

y

and

III. B LOCK D ECISION -F EEDBACK E QUALIZATION In this section we review conventional BDFEs. We focus on the minimum mean-square error (MMSE) BDFE. The BDFE consists of a feedforward equalizer represented by the M ×N matrix G(r) on the rth receive antenna, the decision making device, and the feedback equalizer represented by the M ×M matrix B, as shown in Figure 1(a). Hence, a soft estimate of s is computed as:

(2)

˜s =

q=−Q/2

Nr X

G(r) y(r) − Bˆs,

r=1

where N is the block size over which this channel model holds, and L and Q satisfy the following conditions: C1) LT ≥ τmax ; C2) Q/(N T ) ≥ 2fmax , where T is the symbol period. In this expansion model, L represents the delay-spread (expressed in multiples of T , the delay resolution of the model), and Q represents the Doppler-spread (expressed in multiples of 1/(N T ), the Doppler resolution of the (r) model). Note that the coefficients hq,l remain invariant over a period of length N T , and may change from block to block. Substiute (2) in (1), the received sample at the rth receive antenna at time index n can be written as: y (r) (n) =

Q/2 L X X

(r)

ej2πqn/N hq,l s(n − l) + η (r) (n),

(3)

l=0 q=−Q/2

Suppose we want to transmit a symbol burst s(n), where s(n) 6= 0 for n ∈ / {0, 1, . . . , M − 1}, M = N − Lzp and Lzp ≥ L (zp refers to zero-padding). The input-output relation in (3) can also be written in block form. Defining the M × 1 symbol block as s := [s(0), . . . , s(M − 1)]T , the received sample block at the rth receive antenna y(r) := [y (r) (0), . . . , y (r) (N − 1)]T , can be

ˆs =Q (˜s) ,

(7)

where Q (·) is the quantizer used by the decision device. The MMSE BDFE is designed to minimize the error at the decision device. Define the error vector e = ˜s − s. Assuming past decisions are correct (ˆs = s) we get e=

Nr X r=1

  G(r) H(r) Tzp s + η (r) − (B + IM ) s.

(8)

Defining G := [G(1) , . . . , G(Nr ) ], (8) can be written as: e = GHTzp s + Gη − (B + IM ) s.

(9)

Similar to [1], we obtain the MMSE feedforward and the feedback block equalizers under the constraint that (B + IM ) is an upper triangular matrix as: “ ”−1 −1 G = (B + IM ) (HTzp )H R−1 (HTzp )H R−1 η (HTzp ) + Rs η , B =LH − IM ,

(10)

where L is a lower triangular matrix, obtained  by Cholesky −1 = LDLH . factorization of (HTzp )H R−1 η (HTzp ) + Rs The above approach is rather complex in particular for large N . In the next section, we will therefore develop a DFE equalizer,

with TV FIR feedforward and TV FIR feedback filters. Doing this we lower the design and the equalization complexity at the cost of a slightly lower performance.

Hence we can write ˜ s =Rzp

Nr X r

IV. T IME -VARYING FIR D ECISION -F EEDBACK E QUALIZATION

+ Rzp

In this section, we will apply a TV FIR feedforward filter g (n; ν) on the rth receive antenna and a TV FIR feedback filter b(n; ν), as depicted in Figure 1(b) . Hence, a soft estimate of s(n − d) is computed as Nr X ∞ X

g (r) (n; ν)y (r) (n − ν) −

∞ X

Nr X

(11)

Q′ /2

− Rzp

(r)

hq,l Dq Zl Tzp s

q=−Q/2 l=0

L′ X

X

(r)

gq′ ,l′ Dq′ Zl′ η (r)

q ′ =−Q′ /2 l′ =0

Q′′ /2

X

q ′′ =−Q′′ /2

L′′X +∆−1

bq′′ ,l′′ Dq′′ Zl′′ Tzpˆs.

(15)

l′′ =∆

Defining p = q + q ′ and k = l + l′ , and using the property ′ Zl′ Dq = e−j2πql /N Dq Zl′ , (15) can be written as:

b(n; ν)ˆ s(n − ν),

sˆ(n − d) =Q (˜ s(n − d)) ,

Q/2 L X X

(r) gq′ ,l′ Dq′ Zl′

q ′ =−Q′ /2 l′ =0

(Q+Q′ )/2

ν=−∞

r=1 ν=−∞



L X

X

r

(r)

s˜(n − d) =

Q′ /2

′ L+L X

X

˜ s =Rzp

fp,k Dp Zk Tzp s

p=−(Q+Q′ )/2 k=0

where d is the synchronization delay. Since the doubly-selective channel h(r) (n; ν) was modeled by the BEM, it is also convenient to design the TV FIR feedforward and feedback filters, g (r) (n; ν) and b(n; ν), using the BEM. This will allow us to turn a large design problem into an equivalent small design problem, containing only the BEM coefficients of the doubly-selective channel and the BEM coefficients of the TV FIR feedforward and feedback filters. Using the BEM, we design the TV FIR feedforward filter g (r) (n; ν) to have L′ + 1 taps, where the time-variation of each tap is modeled by Q′ + 1 complex exponential basis functions. Similarly we design the TV FIR feedback filter b(n; ν) to have L′′ taps, where the time-variation of each tap is captured by Q′′ + 1 complex exponential basis functions. Hence, we can write the TV FIR feedforward and feedback filters as: Q′ /2



g (r) (n; ν) =

L X

l′ =0

b(n; ν) =

X

δν−l′

L′′X +∆−1

(r)



gq′ ,l′ ej2πq n/N ,

(12a)

q ′ =−Q′ /2 Q′′ /2

X

δν−l′′

l′′ =∆

bq′′ ,l′′ ej2πq

′′

n/N

,

(12b)

q ′′ =−Q′′ /2

where ∆ is chosen to be larger than the synchronization delay; i.e. ∆ ≥ d + 1. since we can only feedback the past decisions and not the future ones. Instead of continuing to work on the sample level, it is easier to switch to the block level at this point. On the block level, (11) corresponds to estimating s as in (7) but with G(r) and B constrained to ′

G

(r)

=Rzp

L X

Q′ /2

X

(r)

gq′ ,l′ Dq′ Zl′ ,

(13)

l′ =0 q ′ =−Q′ /2

B =Rzp

L′′X +∆−1 l′′ =∆

Q′′ /2

X

bq′′ ,l′′ Dq′′ Zl′′ Tzp ,

(14)

q ′′ =−Q′′ /2

where Rzp := [0M ×d IM 0M ×(Lzp −d) ]. It is clear that this requires 0 < d ≤ Lzp . In practice, we will take Lzp = max{L, d}.

+ Rzp

Nr X

Q′ /2

X



L X

(r)

gq′ ,l′ Dq′ Zl′ η (r)

r=1 q ′ =−Q′ /2 l′ =0 Q′′ /2

− Rzp

X

q ′′ =−Q′′ /2

L′′X +∆−1

bq′′ ,l′′ Dq′′ Zl′′ Tzpˆs,

(16)

l′′ =∆

where fp,k :=

Nr X r=1

Q′ /2

X



L X





(r)

(r)

e−j2π(p−q )l /N gq′ ,l′ hp−q′ ,k−l′ . (17)

q ′ =−Q′ /2 l′ =0

We can further rewrite (16) as ˜ + ˜s = (f T ⊗ IM )As

Nr X

˜ T ⊗ IM )Aˆ ˜ (r) − (b ˜s (g(r)T ⊗ IM )Cη

r=1

˜ T ⊗ IM )Aˆ ˜ − (b ˜ s, ˜ + (gT ⊗ IM )(IN ⊗ C)η = (f T ⊗ IM )As r (18) where we have f := [f−Q/2−Q′ /2,0 , . . . , f−Q/2−Q′ /2,L+L′ , . . . , (r) (r) (r) fQ/2+Q′ /2,L+L′ ]T , g(r) := [g−Q′ /2,0 , . . . , g−Q′ /2,L′ , . . . , gQ′ /2,L′ ]T , g := [g(1)T , . . . , g(Nr )T ]T , and the augmented vector ˜ := [0(Q+Q′ −Q′′ )/2×(L+L′ +1) | 0∆×1 bT ′′ 0u×1 , . . . , b −Q /2 0∆×1 bTQ′′ /2 0u×1 | 0(Q+Q′ −Q′′ )/2×(L+L′ +1) ]T , where u = L + L′ − L′′ − ∆ − 1 and bq′′ := [bq′′ ,∆ , · · · , bq′′ ,L′′ +∆−1 ]T , ˜ in a way to match f in dimension. Defining i.e. we construct b the matrices A and C as     D−Q′ /2 Z0 D−Q/2−Q′ /2 Z0     .. ..     . .        A := D−Q/2−Q′ /2 ZL+L′  , C := D−Q′ /2 ZL′  ,     .. ..     . . DQ′ /2 ZL′ DQ/2+Q′ /2 ZL+L′ ˜ and C ˜ in (18) are A ˜ := (I(Q+Q′ +1)(L+L′ +1) ⊗ the matrices A ˜ Rzp )ATzp and C := (I(Q′ +1)(L′ +1) ⊗Rzp )C, respectively. Note that the term in fp,k corresponding to the rth receive antenna is related to a 2-dimensional convolution of the BEM coefficients of the doubly-selective channel for the rth receive antenna and the BEM coefficients of the TV FIR feedforward filter for the

rth receive antenna. This allows us to derive a linear relationship between f and g. We first define the (L′ + 1) × (L′ + L + 1) Toeplitz matrix   (r) (r) hq,0 . . . hq,L 0   (r) .. .. . Tl,L′ +1 (hq,l ) :=  . .   (r) (r) 0 hq,0 . . . hq,L (r)

We then define Hq(r) := Ωq Tl,L′ +1 (hq,l ), where Ωq := diag{[1, ′ e−j2πq/N , . . . , e−j2πqL /N ]T }, and introduce the (Q′ + 1)(L′ + 1) × (Q + Q′ + 1)(L + L′ + 1) block Toeplitz matrix  (r)  (r) H−Q/2 . . . HQ/2 0   .. .. . Tq,Q′ +1 (Hq(r) ) :=  . .   (r) (r) 0 H−Q/2 . . . HQ/2 Introducing the definitions H(r) := Tq,Q′ +1 (Hq(r) ) and H := [H(1)T , . . . , H(Nr )T ]T , we can derive from (17) that f T = gT H.

(19)

Assuming past decisions are correct, the error at the decision device e = ˜s − s can be written as: ˜ ˜ + (gT ⊗ IM )(IN ⊗ C)η e =(g H ⊗ IM )As r ˜ T ⊗ IM )A ˜ + IM )s. − ((b T

(20)

Now we can write the mean-square error M SE = E{eH e} as: ˜ sA ˜ H (HH g∗ ⊗ IM )} M SE = tr{(gT H ⊗ IM )AR ˜ H )(g∗ ⊗ IM )} ˜ η (IN ⊗ C + tr{(gT ⊗ IM )(IN ⊗ C)R r

r

˜ T ⊗ IM )AR ˜ ∗ ⊗ IM )} ˜ sA ˜ H (b + tr{(b ˜ ∗ ⊗ IM )}} ˜ H (b ˜ sA − 2ℜ{tr{(gT H ⊗ IM )AR ˜ s }} − 2ℜ{tr{(gT H ⊗ IM )AR ˜ T ⊗ IM )AR ˜ s }} + tr{Rs }. + 2ℜ{tr{(b

(21)

Introducing the properties tr{(xT ⊗ IM )X} = xT red(X), tr{(xT ⊗ IM )X(x∗ ⊗ IM )} = xT red(X)x∗ , where red(·) splits the matrix up into M × M submatrices and replaces each submatrix by its trace, then (21) can be rewritten as:   M SE = gT HRA˜ HH + RC˜ g∗ ˜ T + eT )R ˜ (b ˜ ∗ + ed ) + (b d A ˜ ∗ } − 2ℜ{gT Hr ˜ }, − 2ℜ{gT HR ˜ b A

A

(22)

˜ η (IN ⊗ ˜ sA ˜ H ), R ˜ := red((IN ⊗ C)R where RA˜ := red(AR r r C H ˜ )) and r ˜ := red(AR ˜ s ). Note also that we have used the C A ˜ s A)e ˜ d , where ed is a (Q + Q′ + fact that tr{Rs } = eTd red(AR ′ 1)(L + L + 1) × 1 unit vector with the 1 in position d(Q + Q′ + 1) + (Q + Q′ )/2 + 1. By solving ∂M SE/∂g = 0, we then obtain: −1  ˜ T + eT )R ˜ HH HR ˜ HH + R ˜ , (23) gT = (b d A C A

T where we used the fact that rH ˜ , substituting (23) in ˜ = ed RA A (22) we get:

˜ T + eT )R⊥ (b ˜ ∗ + ed ), M SE = (b d

(24)

−1  HRA˜ R⊥ = RA˜ − RA˜ HH HRA˜ HH + RC˜ −1  H −1 −1 = RA . ˜ H ˜ + H RC

(25)

where

The last line in (25) is obtained by applying the matrix inversion Lemma. By defining the selection matrix P∆ as: ˛ ˛   ˛ ˛ T ˛ ˛ ˛

˛ P∆ = 0k×v ˛˛˛ ˛ ˛ ˛

I ′′ ⊗ (e∆ ⊗ I ′′ ) Q /2 L

0

0

0

1

0

0

0

T I ′′ ⊗ (e∆ ⊗ I ′′ ) Q /2+1 L

˛ ˛ ˛ ˛ ˛0 ˛ k×v ˛ ˛ ˛ ˛



(26) where k = Q′′ L′′ + 1, v = (Q + Q′ − Q′′ )(L + L′ + 1)/2 and e∆ is an L + L′ + 1 unit vector with the 1 in position (∆ + 1), we can write (24) as: M SE = [bT−Q′′ /2 , · · · , bT−1 , 1, bT0 , · · · , bTQ′′ /2 ] R∆ b∗ , (27) {z } | bT

where R∆ = P∆ R⊥ PH ∆ . The formula in (27) is a quadratic form that is minimized by choosing the TV FIR feedback filter coefficients as: R−1 ex b = T ∆−1 (28) ex R∆ ex where ex is a Q′′ L′′ + 1 unit vector, with the 1 in the position (Q′′ L′′ /2 + 1). By substituting (28) in (23) we can solve for the TV FIR feedforward equalizer coefficients. In comparing the two approaches, BDFE and TV FIR DFE, we consider two types of complexity: design complexity and implementation complexity. The design complexity is the complexity associated with computing the equalizer (feedforward and feedback) coefficieints. To design the BDFE, we need a matrix inversion of size M × M to compute the feedforward matrix G, and a Cholesky factorization of an M × M matrix to compute the feedback matrix B. Both require about O(M 3 ) flops. On the other hand, to design the TV FIR DFE equalizer, we need to compute the inverse of a K × K matrix to obtain the feedforward filter coefficients, where K = (Q + Q′ + 1)(L + L′ + 1), and we need to compute the inverse of a (Q′′ +1)L′′ +1×(Q′′ +1)L′′ +1 matrix to obtain the feedback filter coefficients. So, provided that K and (Q′′ + 1)L′′ + 1 are less than M , the design complexity of the TV FIR DFE is less than the design complexity of the BDFE. In practice and also in our simulations, the TV FIR DFE parameters Q′ , L′ , Q′′ and L′′ always satisfy K < M and (Q′′ + 1)L′′ + 1 ≪ M . The implementation complexity is the complexity associated with estimating a block of M samples. For the BDFE, to estimate a block of M samples we require M N multiply-add (MA) operations per receive antenna coming from the feedforward part, and (M − 1)M/2 MA operations coming from the feedback part. To estimate the same block using the TV FIR DFE we require M (L′ + 1) MA operations per receive antenna coming from the feedforward part and M L′′ MA operations coming from

0

0

10

10 BDFE TV FIR DFE TV FIR LE

BDFE TV FIR DFE TV FIR LE

−1

10 −1

BER

BER

10

−2

10

−2

10

−3

10

−3

10

−4

0

5

10

15

Fig. 2.

20 SNR (dB)

25

30

35

40

BER vs. SNR for N r = 1

the feedback part. So, provided that L′ is less than the block size N , and L′′ is less than M/2, the implementation complexity of the TV FIR DFE is lower than the implementation complexity of the BDFE. In practice and also in our simulations, the TV FIR DFE parameters Q′ , L′ , Q′′ and L′′ always satisfy L′ ≪ N and L′′ ≪ M/2. V. S IMULATIONS In the simulations, we consider a SISO as well as SIMO (Nr = 2) system. The channel is assumed to be doubly-selective, simulated according to the BEM with the following parameters: • transmitted block size N = 800; • symbol period T = 25µsec; • maximum delay spread τmax = 75µsec • maximum Doppler spread fmax = 100Hz; • number of basis functions Q = 2⌈fmax N T ⌉ = 4; • channel order L = ⌈τmax /T ⌉ = 3. Nr 1 2

Q′ 20 12

L′ 20 12

Q′′ 4 4

L′′ 3 3

TABLE I F ILTER COEFFICIENTS FOR DIFFERENT RECEIVE ANTENNAS

The proposed number of equalizer parameters for different receive antennas is listed in table I. In all simulations, QPSK ′ signaling is used, the delay d = ⌊ L+L 2 ⌋ + 1, and ∆ = d + 1. The performance is measured in terms of the BER vs. SNR. In each plot, we draw comparisons between the proposed TV FIR DFE, the TV FIR LE (no feedback case), and the BDFE. As shown in Figure 2 (SISO Nr = 1), the SNR loss for the TV FIR DFE compared to BDFE is less than 1dB at BER = 10−2 . On the other hand, the proposed TV FIR DFE outperforms the TV FIR LE, and the latter suffers from an early error floor. Similarly, in Figure 3 (SIMO Nr = 2), we see that the performance of the proposed TV FIR DFE almost coincides with that of the BDFE, whilst it significantly outperforms the TV FIR LE equalizer.

10

0

2

4

6

8

Fig. 3.

10 SNR (dB)

12

14

16

18

20

BER vs. SNR for N r = 2

VI. C ONCLUSION In this paper, we have proposed a TV FIR DFE to compensate for doubly selective channel distortion. The proposed equalizer provides a good compromise between complexity and performance compared to BDFE and the TV FIR LE. We show that by properly choosing the feedforward and feedback filter coefficients, we can reach the performance of the BDFE. We also show that the TV FIR DFE outperforms the TV FIR LE. R EFERENCES [1] A. Stamoulis, G. B. Giannakis and A. Scaglione, “Block FIR DecisionFeedbak Equalizers for Filterbank Precoded Transmissions with Blind Channel Estimation Capabilities,” IEEE Trans. Signal Processing, vol. 49, pp. 69–83, January 2001. [2] D. Williamson, R. A. Kennedy and G. W. Pulford, “Block Decision Feedback Equalization,” IEEE Trans. Commun., vol. 40, pp. 255–264, February 1992. [3] N. Al-Dhahir and J. M. Cioffi, “MMSE Decision-Feedback Equalizers: Finite Length Results,” IEEE Trans. Inform. Theory, vol. 41, pp. 961–975, July 1995. [4] I. Barhumi, G. Leus and M. Moonen, “Time Varynig FIR Equalization of Doubly-Selective Channels,” in Int. Conf. on Comm., (Anchorage, Alaska USA), May 11-15 2003. accepted. [5] G. Leus, I. Barhumi and M. Moonen, “MMSE-Time Varying FIR Equalization of Doubly-Selective Channels,” in Int. Conf. on ASSP, (Hong Kong), April 6-10 2003. accepted. [6] I. Barhumi, G. Leus and M. Moonen, “Time Varynig FIR Equalization of Doubly-Selective Channels,” IEEE Trans. Wireless Commun., Dec. 2002. submitted. [7] G. B. Giannakis and C. Tepedelenlio˘glu , “Basis Expansion Models and Diversity Techniques for Blind Identification and Equalization of Time Varying Channels,” Proc. IEEE, vol. 86, no. 10, pp. 1969–1986, Oct., 1998. [8] A. M. Sayeed and B. Aazhang, “Joint Multipath-Doppler Diversity in Mobile Wireless Communications,” IEEE Trans. Commun., vol. 47, pp. 123– 132, January 1999. [9] X. Ma, and G. B. Giannakis, “Maximum-Diversity Transmissions over Doubly-Selective Wireless Channels,” IEEE Trans. Inform. Theory, revised, February 2003. see also Proc. of ”WCNC”., vol. 1, pp. 497-501, Orlando, FL, March 17-21, 2002. [10] G. Leus, S. Zhou and G. B. Giannakis, “Orthogonal Multiple Access over Time- and Frequency-Selective Fading,” IEEE Trans. Inform. Theory, revised, June, 2002. see also Proc. of the ACSSC , Pacific Grove, California, November 4-7, 2001.