22.08.2013 Views

Untitled

Untitled

Untitled

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

where subscript c denotes continuous time, h c(t;T) is the<br />

convolution of the spectral pulse f c (tr) (t), the TV impulse response<br />

f c (ch) (t;τ ), the receive-filter fc (rec) (t), s(l) is the sequence of input<br />

symbols, and vc(t) is the noise process. h c(t;T) can be truncated to<br />

an order L, which is common practice in communication<br />

applications. On approach in characterizing the variation of the<br />

impulse response h(n;l) is to consider it as a stochastic process in<br />

the time index n. In communication, tracking the variations of the<br />

channel taps is of importance. To this end, linear TV channel is<br />

assumed in [6].<br />

III. FADING CHANNELS: BASIS EXPANSION MODELS<br />

In this paper, we model the doubly selective channel h(n;v)<br />

using a basis expansion model (BEM). In this BEM, the channel is<br />

modeled as a TV FIR filter, where each tap is expressed as a<br />

superposition of complex exponential basis functions with the<br />

frequencies on a DFT grid, as described next.<br />

Let us first make the following assumptions:<br />

A1) The delay spread is bounded by τ max.<br />

A2) The Doppler spread is bounded by fmax.<br />

Under assumption A1) and A2), it is possible to accurately<br />

model the doubly selective channel h(n;v) for n∈{0,1,..., N −1}<br />

as<br />

(see Fig.2)<br />

h(<br />

n;<br />

v)<br />

=<br />

L<br />

∑ v−l<br />

Q / 2<br />

j 2πqn<br />

/ N<br />

∑e<br />

hq,<br />

l<br />

l=<br />

0 q=<br />

−Q<br />

/ 2<br />

δ (4)<br />

where L and Q satisfy the following conditions:<br />

C1) LT≥τ max;<br />

C2) Q/(NT) ≥2f max;<br />

h− Q /2,0 h− Q /2,1 h− Q /2,2<br />

Q/2, L<br />

j2 Q/2 n/ N<br />

h Q /2,0 h Q /2,1 h Q /2,2 hQ/2, L e π<br />

h −<br />

Fig. 2. BEM channel.<br />

e<br />

− j2 πQ/2<br />

n/ N<br />

where T is the symbol period. In this expansion model, L<br />

represents the discrete delay spread (expressed in multiples of T,<br />

the delay resolution of the model), and Q/2 represents the discrete<br />

Doppler spread (expressed in multiples of 1/(NT), the Doppler<br />

resolution of the model). Note that the coefficients h q,l remain<br />

invariant for n∈{0,1,..., N −1}<br />

and hence are the BEM<br />

coefficients of interest[7].<br />

IV. BLOCK DATA MODEL.<br />

The discrete-time baseband equivalent model of the OFDM<br />

system under consideration is shown in Fig. 3. In the transmitter,<br />

the transmitted high-speed data is first converted into parallel data<br />

of N subcarriers by the serial-to-parallel converter (S/P). After<br />

each symbol is modulated by the corresponding subcarrier, it is<br />

sampled and D/A converted. The sampled symbols, implemented<br />

by an inverse fast Fourier transform (IFFT), can be expressed as<br />

follows:<br />

x<br />

∑ −<br />

=<br />

=<br />

1 N<br />

n<br />

m 0<br />

X<br />

e<br />

j 2πnm<br />

/ N<br />

m<br />

, 0 ≤ n ≤ N<br />

(5)<br />

where xn represents the nth sample of the output of the IFFT[6].<br />

Input Bit<br />

Stream<br />

Output Bit<br />

Stream<br />

S/P<br />

P/S<br />

X m<br />

~<br />

IFFT<br />

Equalizer<br />

Y m<br />

Xm<br />

FFT<br />

yn<br />

P/S<br />

xn<br />

D/A<br />

S/P A/D<br />

AWGN<br />

wn<br />

Multipath<br />

Channel<br />

Fig. 3. A baseband block diagram for an OFDM system.<br />

Assuming that the multipath fading channel consists of L<br />

discrete paths, the received signal can be written as<br />

L−1<br />

∑<br />

y = h x + ω = h x + h x + ... + h x 1+ ωn<br />

n n, l n−l n n,0 n n,1 n−1 n, L−1 n− L+<br />

l=<br />

0<br />

where hn,l and ω n represent the complex random variable for the<br />

lth path of the CIR and additive white Gaussian noise (AWGN) at<br />

time n, respectively. For reason of simplicity, a cyclic extension of<br />

length N G, used to avoid ISI and to preserve the orthogonality of<br />

subcarriers, is not shown in Fig.4 and not used in this paper’s<br />

simulation. The demodulated signal in the frequency domain is<br />

obtained by taking the FFT of y n as<br />

m<br />

N−1 L−1<br />

∑∑<br />

k= 0 l=<br />

0<br />

k<br />

( m−k) l<br />

− j2 πlk/<br />

N<br />

Y = X H e + W<br />

L−1<br />

⎡ 0 − j2 πlm/<br />

N⎤<br />

∑ He l Xm<br />

l=<br />

0<br />

= ⎢ ⎥<br />

⎣ ⎦<br />

N−1 L−1<br />

∑∑<br />

+ XH e + W<br />

k≠ m l=<br />

0<br />

( m−k) − j2 πlk/<br />

N<br />

k l m<br />

= α X + β + W , 0≤m≤ N −1<br />

m m m m<br />

( m k)<br />

where Wm denotes the FFT of wn . Also, H represents the<br />

l<br />

−<br />

FFT of a time-variant multipath channel hn,l<br />

, which is modeled<br />

using BEM. It can be represented by<br />

1<br />

H h e<br />

N −1<br />

( m−k) − j2 π n( m−k)/ N<br />

l = ∑ n, l<br />

N n=<br />

0<br />

N −1<br />

Q /2 ⎡<br />

∑ ∑<br />

j2 πqn/ N⎤ − j2 πn(<br />

m−k)/ N<br />

hql , e e<br />

n= 0 q=−Q/2 1<br />

= ⎢ ⎥<br />

N ⎣ ⎦<br />

m<br />

+<br />

(6)<br />

(7)<br />

(8)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!