Adaptive Channel Modeling For MIMO Wireless Communications: Chengjin Zhang Robert R. Bitmead
Adaptive Channel Modeling For MIMO Wireless Communications: Chengjin Zhang Robert R. Bitmead
Adaptive Channel Modeling For MIMO Wireless Communications: Chengjin Zhang Robert R. Bitmead
Communications
Chengjin Zhang
Robert R. Bitmead
Abstract The application of state-space-based subspace system identification methods to training-based estimation for timevarying MIMO frequency-selective channels is explored with the
motivation of possible parsimonious parametrization and direct
model complexity control. The comparison between the statespace-based channel estimation algorithm and the FIR-based
RLS algorithm shows the former is a more robust modeling
approach than the later.
I. I NTRODUCTION
Digital communication using multiple transmit and receive
antennas has been one of the most important technical developments in modern communications. In a rich scattering
environment, multi-input-multi-output (MIMO) systems offer
significant capacity gain at no cost of extra spectrum [1].
Furthermore, space-time channel codes [2] can be applied to
utilize multiple transmit antennas to build spatial redundancy
in the transmitted signal such that maximum spatial diversity
is achieved. Coherent space-time processing schemes assume
the availability of a channel model at the receiver. Therefore,
this model needs to be estimated at the receiver end.
For slowly-varying channels, training-based channel estimation is very common in practice. Most channel estimation
schemes assume an FIR channel model. FIR models for wireless channels have been widely used for their simplicity and
guaranteed stability. An FIR model represents the sub-channels
of a MIMO channel with separately parametrized finite-length
impulse responses. However, when the spatial subchannels
are correlated with each other due to insufficient separation
between adjacent antennas, an FIR model can be very nonparsimonious and contain excessive redundancy. Therefore it
may be beneficial to adopt a state-space model which treats
the whole channel as a single entity and hence captures the
structure in the channel while allowing a more parsimonious
description of the MIMO channel.
It is the subject of this paper to explore the application
of state-space models to represent MIMO frequency-selective
channels in the hope for possible parsimonious parameterization. Since any channel model is essentially an approximation
of the physical channel, the goal of channel modeling is to
find a model that is as close to the real channel as possible
and maintains manageable complexity as well. It is found that
Normalized singular values of the Hankel matrix, 1000 2x2 channel realizations
Normalized singular values of the Hankel matrix, 1000 10x10 channel realizations
1
correlated channel
uncorrelated channel
correlated channel
uncorrelated channel
0.9
0.9
0.8
0.8
0.7
0.7
0.6
0.6
0.5
0.5
0.4
0.4
0.3
0.3
0.2
0.2
0.1
0.1
10
(a) 2 2 channels.
10
15
20
25
30
35
40
45
50
(b) 10 10 channels.
Fig. 1. Hankel singular values of MIMO channels (normalized with respect to the largest value). = 0.8 for the spatially correlated channels; = 0 for
spatially white channels.
correlation increases the number of negligible small Hankel and possess the structure
2 . . . m1
In the simulation, the MIMO channel is assumed to be
1
. . . m2
...
1
be represented by assembling all the taps of the subchannel
impulse responses with the same delay into a matrix H where represents the correlation coefficient between the
fading of two adjacent transmit or receive antennas. The value
and represent the channel as a series of matrices
of is chosen to be 0.8 in the simulation.
.
[H0 , H1 , . . . , HL1 ], H =
channels
..
are computed and averaged over 1000 channel re
.
alizations. [The H -norm of the difference between a MIMO
hp1 ( ) hp2 ( ) . . . hpm ( )
(1) system and its low-order approximation is related to the sum of
the neglected Hankel singular values, see [6] (Theorem 7.11).]
The simulated channels are taken to be square symmetric,
Figure 1 shows the difference between the distributions of
and the length of the impulse responses is set to 6. The channel
Hankel singular values for spatially white channels and spataps are generated as zero-mean circularly-symmetric comtially correlated channels. It can be seen that spatial correlation
plexity Gaussian (ZMCSCG) random variables. The impulse
increases the number of negligible singular values and hence
response of each subchannel is normalized to have unit energy.
allows a low-order state-space model that is a high-quality
The MIMO channel is assumed to the uncorrelated in delay
approximation of the original full-order channel. Furthermore,
but correlated in spatial dimension, i.e.
the larger the size of the MIMO channel is, the lower the order
of a high-quality channel model can be achieved. Therefore,
0, k 6= l
E(vec(Hk )vec(Hl )H ) =
(2) a state-space model is capable of providing good low-order
R, k = l
approximations of large spatially correlated MIMO channels.
where vec() is an operator that stacks the columns of a
xk+1 = Axk + Buk
matrix on top of each other to form a vector. We adopt a
simple model for the spatial correlation structure of H [5]
yk = Cxk + Duk + nk
(6)
H
R
=
=
(3)
(4)
where
uk m 1 input vector at time k;
yk p 1 output vector at time k;
xk n 1 state vector at time k;
nk p 1 white Gaussian noise vector at time k.
On the other hand, the order of an FIR channel can only be
controlled by the length of the impulse response. This is an
inferior approach compared to the model reduction in statespace models, as will be shown later in Section IV.
III. S UBSPACE S YSTEM I DENTIFICATION AND MIMO
C HANNEL E STIMATION
The MIMO FIR model can be fitted by forming a parameter
vector from the columns of H and a regressor vector from
a suitable rearrangement of the input signals at the transmit
antennas. Then Recursive Least Squares is applied to this
vector problem.
However, for state-space models, a special class of estimation algorithms is needed. Subspace system identification
(SSI) refers to a class of recent algorithms, such as MOESP
and N4SID, which apply input-output system identification
methods to determine directly a state-space realization of
a system. The key idea of SSI methods is to estimate the
extended observability matrix through the projection of future
input-output data onto past input-output data based on the
relationship between Hankel matrices of the input and output
given by
Y1,i,M = i X1,M + Hi U1,i,M + N1,i,M , i > n.
(7)
y1
y2
yM
y2
y3
yM+1
Y1,i,M = .
(8)
.
..
..
..
..
.
.
.
yi
yi+1
yi+M1
D
0
0
C
CB
CA
D
0
0
CAB
..
CB D 0 0
i = . , Hi =
.
..
..
..
..
CAi2
.
.
.
.
CAi1
CAi2 B
D
(9)
Then the system matrices A, B, C and D are computed
based on the estimated observability matrix and singular value
decomposition (SVD) algorithm, see [3], [4], [7] for in-depth
treatments.
Based on the channel model given in (6), SSI methods
require the input to satisfy the following requirements for the
channel to be identifiable.
1) The input uk is uncorrelated with the additive Gaussian
white noise nk .
2) The input uk is persistently exciting of order of at least
2 times the maximum order of the channel.
3) The symbols in the input sequence are contiguous and,
for consistency, the number of inputs goes to infinity.
The first assumption is usually satisfied for wireless communication systems. The second one requires the training
AND I NTERPRETATION
2.5
2.5
H error in dB
2
1.5
2
1.5
H error in dB
0.5
0
0.5
0
50
100
150
200
250
50
100
symbol index
= 0.90, 100 data realizations, 20 channel realizations
200
250
2
Recursive MOESP
RLS
steadystate H error
Recursive MOESP
RLS
1.5
1.5
steadystate H error
150
symbol index
0.5
10
15
20
25
30
35
convergence point
40
45
50
0.5
0
20
55
40
60
100
120
convergence point
140
160
180
Estimation error of MOESP and RLS for 2 2 spatially correlated channels with = 0.8.
H error in dB
2
1
0
50
100
150
200
250
symbol index
= 0.93, 100 data realizations, 10 channel realizations
3
steadystate H error
Fig. 2.
80
Recursive MOESP
RLS
2.5
2
1.5
1
0.5
0
10
20
30
40
50
60
convergence point
70
80
90
100
Fig. 3. Estimation error of MOESP and RLS for 3x3 spatially correlated
channels with = 0.8, = 0.93.
V. C ONCLUSION
In this paper, the state-space model is proposed for modeling
MIMO wireless channels with the motivation for a more
parsimonious parameterization. A recursive subspace system
ACKNOWLEDGMENT
This research was supported by CORE Grant No. 02-10109
sponsored by Ericsson.
A PPENDIX
A. SSI for non-contiguous data
Notice that the third assumption of conventional SSI algorithms
(Section III) places limitations on the application of SSI methods
to channel estimation in wireless communication systems where the
training sequences are usually not contiguous in time. Instead, they
lie in the overhead or mid-amble of a frame and are separated by data
symbols. In this circumstance, the state evolution of the received data
must be restarted at the frame boundaries. This is at variance with
the standard formulation of SSI.
A suitable modification of the SSI algorithms for non-contiguous
training data is proposed in [10]. It is shown that the non-contiguousdata approach is similar in estimation power to the contiguous-data
approach when the length of the training sequence is sufficiently
larger than the dimension of the extended observability matrix i
(9).
U (t) u (t + 1)
U (t + 1)
=
Y (t) y (t + 1)
Y (t + 1)
R11 (t)
0
u (t + 1) Q 0
(11)
=
R21 (t) R22 (t) y (t + 1) 0 1
R11 (t + 1) 0mipi
0mi1
=
Q0
R21 (t + 1) R22 (t) (t + 1)
where the pi 1 vector (t + 1) is obtained through a series of
Givens rotation on the large lower triangular matrix. (0, 1] is an
exponential forgetting factor which could be used to adjust the effect
of the past input-output data on the estimate of the current channel
model. R22 (t + 1) can be found by taking an LQ factorization of
[R22 (t) (t + 1)] as follows.
[R22 (t) (t + 1)] = R0 Q0
Therefore, by substituting (12) into (11),
U (t + 1)
R11 (t + 1)
=
Y (t + 1)
R21 (t + 1)
R11 (t + 1)
=
R21 (t + 1)
R11 (t + 1)
=
R21 (t + 1)
(12)
we obtain
0 0
Q0
R0 Q0
0 0 Q1
R0 Q0 Q2
0
Q00
R0
Q1
00
Q =
Q0 Q2
is a matrix with orthogonal rows, i.e. Q00 (Q00 )T = I. Thus
U (t)
R11 (t)mimi
=
Y (t)
R21 (t)pimi
0
Q
R22 (t)pipi
where U (t) and Y (t) are respectively the Hankel matrices of the
input and output data at time t with the form
2
6
6
Y (t) = 6
4
y0
y1
..
.
yi1
y1
y2
..
.
yi
...
...
...
yti+1
yti
..
.
yt
3
7
7
7
5
U (t + 1)
U (t) u (t + 1)
=
(10)
Y (t + 1)
Y (t) y (t + 1)
R22 (t + 1) = R0
In summary, the update from R22 (t) to R22 (t + 1) is a two-step
procedure: computing (t+1) in (11) using Givens rotation followed
by computing R0 in (12).
R EFERENCES
[1] G. J. Foschini and M. J. Gans, On limits of wireless communications in
a fading environment when using multiple antennas, Wireless Personal
Communications, vol. 6, no. 3, pp. 311335, March 1998.
[2] V. Tarokh, H. Jafarkhani, and A. R. Calderbank, Space-time block
coding for wireless communications: Performance results, IEEE J.
Select. Areas Commun., pp. 451460, March 1999.
[3] M. Verhaegen and P. Dewilde, Subspace model identification Part 1.
the output-error state-space model identification class of algorithms,
International Journal of Control, vol. 56, no. 5, pp. 11871210, 1992.
[4] P. Van Overschee and B. De Moor, Subspace identification for linear
systems: theory, implementation, applications.
Dordrecht: Kluwer
Academic Publishers, 1996.
[5] K. Yu and B. Ottersten, Models for mimo propagation channels: a
review, Wireless Communications and Mobile Computing, vol. 2, pp.
273291, Nov. 2002.
[6] K. Zhou and J. C. Doyle, Essentials of robust control, 1st ed. New
Jersey, USA: Prentice Hall, 1997.
[7] L. Ljung, System identification, theory for the user, 2nd ed. New Jersey,
USA: Prentice Hall, 1999.
[8] M. Verhaegen, Identification of the deterministic part of MIMO state
space models given in innovations form from input-output data, Automatica, vol. 30, no. 1, pp. 6174, Jan. 1994.
[9] M. Lovera, T. Gustafsson, and M. Verhaegen, Recursive subspace
identification of linear and non-linear Wiener state-space models,
Automatica, vol. 36, no. 11, pp. 16391650, Nov. 2000.
[10] C. Zhang and R. R. Bitmead, Multiple antenna system equalization
using semi-blind subspace identification methods, Proceedings of 13th
IFAC Symposium on System Identification, Rotterdam, Aug. 2003.