- Research Article
- Open Access
- Published:

# Time-Frequency Based Channel Estimation for High-Mobility OFDM Systems–Part I: MIMO Case

*EURASIP Journal on Advances in Signal Processing*
**volume 2010**, Article number: 549197 (2010)

## Abstract

Multiple-input multiple-output (MIMO) systems hold the potential to drastically improve the spectral efficiency and link reliability in future wireless communications systems. A particularly promising candidate for next-generation fixed and mobile wireless systems is the combination of MIMO technology with Orthogonal Frequency Division Multiplexing (OFDM). OFDM has become the standard method because of its advantages over single carrier modulation schemes on multipath, frequency selective fading channels. Doppler frequency shifts are expected in fast-moving environments, causing the channel to vary in time, that degrades the performance of OFDM systems. In this paper, we present a time-varying channel modeling and estimation method based on the Discrete Evolutionary Transform to obtain a complete characterization of MIMO-OFDM channels. Performance of the proposed method is evaluated and compared on different levels of channel noise and Doppler frequency shifts.

## 1. Introduction

The major challenges in future wireless communications systems are increased spectral efficiency and improved link reliability. The wireless channel constitutes a hostile propagation medium, which suffers from fading (caused by destructive addition of multipath components) and interference from other users. Diversity provides the receiver with several (ideally independent) replicas of the transmitted signal and is therefore a powerful means to combat fading and interference and thereby improve link reliability. Common forms of diversity are space-time diversity [1] and space-frequency diversity [2]. In recent years the use of spatial (or antenna) diversity has become very popular, which is mostly due to the fact that it can be provided without loss in spectral efficiency. Receive diversity, that is, the use of multiple antennas on the receiver side of a wireless link, is a well-studied subject [3]. Driven by mobile wireless applications, where it is difficult to deploy multiple antennas in the handset, the use of multiple antennas on the transmitter side combined with signal processing and coding has become known under the name of space-time coding [4] and is currently an active area of research. The use of multiple antennas at both ends of a wireless link (multiple-input multiple-output (MIMO) technology) has been demonstrated to have the potential of achieving extraordinary data rates [5]. The corresponding technology is known as spatial multiplexing [6] or BLAST [7] and yields an impressive increase in spectral efficiency. Most of the previous work in the area of MIMO wireless has been restricted to narrow-band systems. Besides spatial diversity broadband MIMO channels, however, offer to higher capacity and frequency diversity due to delay spread. Orthogonal frequency division multiplexing (OFDM) significantly reduces receiver complexity in wireless broadband systems. The use of MIMO technology in combination with OFDM, that is, MIMO-OFDM [6], therefore seems to be an attractive solution for future broadband wireless systems [8, 9]. However, intercarrier interference (ICI) due to Doppler shifts, phase offset, local oscillator frequency shifts, and multi-path fading severely degrades the performance of OFDM systems [10]. Most of the channel estimation methods assume a linear time-invariant model for the channel, which is not valid for the next-generation, fast-moving environments [11]. Recently a time-frequency varying MIMO-OFDM channel estimation approach is presented where discrete prolate spheroidal sequences are used to obtain a robust time-varying channel estimator that does not require any channel statistics [12]. A time-varying model of the channel can be obtained by employing time-frequency representation methods. Here we present a time-varying MIMO-OFDM channel estimation based on the discrete evolutionary representation of the channel output. The Discrete Evolutionary Transform (DET) [13] provides a time-frequency representation of the received signal by means of which the spreading function of the multi-path, fading, and frequency selective channel can be modeled and estimated.

The rest of the paper is organized as follows. In Section 2, we give a brief summary of the wireless parametric channel model used in our approach and the MIMO-OFDM communication system. Section 3 presents time-varying modeling and estimation of MIMO-OFDM channels via DET. A time-frequency receiver is also given in Section 3 for the detection of data symbols using estimated channel parameters. In Section 4, we give some simulation results to illustrate the performance of our algorithm for different levels of channel noise and Doppler frequency-shifts and compare with other existing methods. Conclusions are drawn in Section 5.

## 2. MIMO-OFDM System Model

In this section we give a brief introduction to the time-varying, parametric communication channel model used in our work and the MIMO-OFDM signal model.

### 2.1. Parametric Channel Model

In wireless communications, the multi-path, fading channel with Doppler frequency-shifts may be modeled as a linear time-varying system with the following impulse response [14–16]:

where are independent Gaussian processes with zero mean and variance, are delay profiles describing the channel dispersion with as the maximum delay, and is the total number of paths. The variance is a measure of the average signal power received at path , characterized by the relative attenuation of that path, . Assuming that the sampling frequency is high enough, so the delays are integer multiples of the sampling interval, a discrete-time channel model can be obtained as [15, 17]

where is the time index, represents the Doppler frequency-shift, is the relative attenuation, and is the delay in path . The Doppler frequency shift , on the carrier frequency , is caused by an object with radial velocity and can be approximated by [15]

where is the speed of light in the transmission medium. In the new generation wireless mobile communication systems, with fast moving objects and high carrier frequencies, Doppler frequency-shifts become significant and have to be taken into consideration. The channel parameters cannot be easily estimated from the impulse response; however the estimation problem can be solved in the time-frequency domain by means of the so called spreading function.

The spreading function is related to the generalized transfer function and the bifrequency function. The generalized transfer function of the linear, time-varying channel is obtained by taking the DFT of with respect to , that is,

where . Now, the channel bi-frequency function is found by computing the DFT of with respect to time variable, :

where . Calculating the inverse DFT of with respect to (or by taking the DFT of in (2) with respect to ), we have the spreading function of the channel as

provided that the Doppler frequency shifts are integer multiples of the frequency sampling interval . displays peaks located at the time-frequency positions determined by the delays and the corresponding Doppler frequencies, with as their amplitudes [17]. In our approach, we extract this information from the spreading function of the received signal and then detect the transmitted data symbol.

### 2.2. MIMO-OFDM Signal Model

In an OFDM communication system, the available bandwidth is divided into subchannels. The input data is also divided into -bit parallel bit streams and then mapped onto some transmit symbols drawn from an arbitrary constellation points where is the time index, and , denotes the frequency or subcarrier index.

Some pilot symbols are inserted at some preassigned positions , known to the receiver: where is the number of pilots, and the integer is the distance between adjacent pilots in an OFDM symbol [10]. The OFDM symbol is obtained by taking -point inverse DFT and then adding a cyclic prefix (CP) of length where is chosen such that , and is the time-support of the channel impulse response. This is done to mitigate the effects of intersymbol interference (ISI) caused by the channel dispersion in time.

A MIMO-OFDM system with transmit and receive antennas is depicted in Figure 1. The incoming bits are modulated to form , where is the transmit antenna index. After parallel-to-serial (P/S) conversion, the signal transmitted by the antenna becomes

where , , and each OFDM symbol is sample long. The channel output suffers from multi-path propagation, fading and Doppler frequency shifts:

where denotes the receiver. is indexing for the wireless time-varying channel between the transmitter and the receiver antennas. The transmit signal is also corrupted by additive white Gaussian noise over the channel. The received signal for the frame can then be written as . The receiver discards the Cyclic Prefix and demodulates the signal using a -point DFT as

where is the Fourier transform of the channel noise at the receiver. If the Doppler effects in all the channel paths are negligible, , then the channel is almost time-invariant within one OFDM symbol. In that case, above equation becomes

where is the frequency response of the channel between the transmitter and receiver antenna, and if there are large Doppler frequency shifts in the channel, then the time-invariance assumption above is no longer valid. Here we consider modeling and estimation of the channels frequency responses and approach the problem from a time-frequency point of view [17].

## 3. Time-Varying Channel Estimationfor MIMO-OFDM Systems

In this section we present a time-frequency procedure to characterize time-varying MIMO-OFDM channels. We also propose a time-frequency receiver that uses the estimated channel fading, delay, and Doppler parameters to recover the transmitted symbols. In the following, we briefly present the Discrete Evolutionary Transform (DET) as a tool for the time-frequency representation of time-varying MIMO-OFDM channels.

### 3.1. The Discrete Evolutionary Transform

Wold-Cramer representation [18] of a nonstationary random signal can be expressed as an infinite sum of sinusoids with random and time-varying amplitudes and phases, or

where is considered a random process with orthogonal increments. This is a generalization of the spectral representation of stationary processes. Priestley's evolutionary spectrum [18, 19] of is given as the magnitude square of the evolutionary kernel . Analogous to the above Wold-Cramer representation, a discrete, time-frequency representation for a deterministic signal with a time-dependent spectrum is possible [13, 20]:

where , is the number of frequency samples, and is a time-frequency evolutionary kernel. A similar representation can be given in terms of the corresponding bi-frequency kernel :

where and are discrete frequencies. Discrete evolutionary transformation (DET) is obtained by expressing the kernels or by means of the signal [13]. Thus, for the representation in (12) the DET that provides the evolutionary kernel , , is given by

where is, in general, a time- and frequency-dependent window.

The DET can be seen as a generalization of the short-time Fourier transform, where the windows are constant. The windows can be obtained from either the Gabor representation that uses nonorthogonal frames or the Malvar wavelet representation that uses orthogonal systems. Details of how the windows can be obtained for the Gabor and Malvar representations are given in [13]. However, for the representation of multipath wireless channel outputs, we consider windows that are adapted to the Doppler frequencies of the channel.

### 3.2. MIMO-OFDM Channel Estimation by Using DET

We will now consider the computation of the spreading function by means of the evolutionary transformation of the received signal. The output of the channel, after discarding the cyclic prefix, for the OFDM symbol can be written using (8) as

where we ignore the additive channel noise for simplicity. The equation above can be rewritten in matrix form as

where

The input-output relation for the whole system results as

where

If is known and then input symbols can be estimated by the following relation,

Now calculating the discrete evolutionary representation of , we get

Comparing the representations of in (21) and (15), we get the kernel as

above relation is also valid at the preassigned pilot positions

where and is a decimated version of the . Note that is again the number of pilots, and is the distance between adjacent pilots. Finally if the pilot symbols are chosen to be orthogonal to each other, the decimated frequency response of the channel between the transmitter and the receiver antennas may be obtained as

Taking the inverse DFT of with respect to and DFT with respect to , we obtain the downsampled spreading function ,

By comparing and , we observe that the channel parameters , , and calculated from can also be estimated from the downsampled spreading function .

In the following, we present a method to estimate the spreading function of the MIMO channel from the received signal. The time-frequency evolutionary kernel of the channel output in the receiver is obtained as

We consider windows of the form , for presented in [17] that depends on the Doppler frequency . This window will give us the correct representation of only when in fact, using the window , the above representation of becomes,

which is the expected result multiplied by . In our estimation procedure, we use windows where is chosen in a discrete set with certain increments, . When coincides with one of the Doppler frequencies in the channel, the spreading function displays a large peak at the time-frequency position , corresponding to delay and Doppler frequency of that transmission path, with magnitude proportional to attenuation . When does not coincide with any of the Doppler frequencies, the spreading function displays a random sequence of peaks spread over all possible delays. Then it is possible to determine a threshold that permits us to obtain the most significant peaks of the spreading function corresponding to possible delays and Doppler frequencies. In our experiments we observed that peaks having amplitudes larger than 65% of the maximum peak are due to an actual transmission; otherwise they are considered as noise. Thus, by searching in the possible Doppler frequency range, we are able to estimate all the parameters of a multi-path, fading, and time-varying MIMO-OFDM channel via the spreading function of the channel.

According to (27), we need the input pilot symbols to estimate the channel frequency response. Here we consider simple, uniform pilot patterns; however improved patterns may be employed as well [11].

### 3.3. Time-Frequency Receiver

After estimating the spreading function and the corresponding frequency response of the channel, data symbols can be detected using a time-frequency receiver given in (20). On the other hand, the channel output in (9) can be rewritten as

where is the bi-frequency function of the channel during OFDM symbol, and the above equation indicates a circular convolution with the data symbols. Based on aforementioned equality it is possible to write the MIMO-OFDM system consisting of transmitter and receiver antennas in a matrix form as

or

where is a matrix; , and are vectors defined by , , and , respectively. Finally, data symbols can be estimated by using a simple time-frequency receiver: .

The exhaustive search for the channel Doppler frequencies may seem to increase the computational cost of the proposed method. However, considering the carrier frequencies and maximum possible velocities in the environment, Doppler frequencies lie in a certain band which can be easily covered by the algorithm. Furthermore, our channel estimation approach does not require any a priori information on the statistics of the channel as in the case of many other channel estimation methods [11].

In the following, we demonstrate time-varying MIMO-OFDM channel estimation performance of our time-frequency-based approach by means of examples.

## 4. Simulations

In the experiments, a 2-input, 2-output MIMO-OFDM system is considered, and the wireless channels are simulated randomly; that is, the number of paths, , the delays, and the Doppler frequency shift of each path are picked randomly. Input data is QPSK coded and modulated onto subcarriers, 16 of which are assigned to pilot symbols. The Signal-to-Noise Ratio (SNR) of the channel noise is changed between 5 and 35 dB, for two different values of the maximum Doppler frequency. Figure 2 depicts an example of the estimated spreading function for a MIMO system, during one OFDM symbol. In our simulations, we use Doppler frequencies normalized by the subcarrier spacing [21]:

The performance of our channel estimation method is investigated and compared with that of a recently proposed time-frequency varying MIMO-OFDM channel estimation approach [12].This method applies the discrete prolate spheroidal sequences to obtain a robust time-varying channel estimator that does not require any channel statistics, similar to our proposed method.

We show the mean square error (MSE) of our channel estimation approach and the 2D Slepian-based two-dimensional channel estimation method [12] for different channel noise levels and for normalized Doppler frequencies and in Figure 3. As seen from the graphs, our method outperforms the Slepian-based approach in terms of both estimation error and robustness against increased Doppler frequencies. Note that our method is capable of estimating and compensating for large Doppler shifts yielding a similar MSE for both and .

We then investigate the bit error rate (BER) versus SNR of the channel noise performance of both channel estimation approaches for different numbers of the pilot symbols in one OFDM block, and show the results in Figure 4. Notice that for a fixed number of pilots, proposed TF-based approach achieves the same BER with about 15 dB less SNR than the Slepian-based method. Note also that increasing the number of pilots improves the BER performance in both methods.

Finally the effect of the number of channel propagation paths is investigated. BER is calculated for both TF and Slepian-based channel estimation approaches by changing the number of paths from to and for and pilots. Results given in Figure 5 show that for the same number of paths, proposed TF-based approach achieves approximately times less BER than the 2D Slepian-based channel estimation method. Also notice that the BER performance goes down rapidly for the number of paths larger than .

## 5. Conclusions

In this work, we present a time-varying estimation of MIMO-OFDM channels for high-mobility communication systems by means of discrete evolutionary transform. The main advantage of the proposed method is that it does not assume any statistics on the communication channel. The parametric channel model used in this approach allows us to obtain a two-dimensional representation for the channel and estimate its parameters from the spreading function. We observe that the method is robust against large variations on the channel frequency response, that is, fast fading. Simulations show that our time-frequency-based method has considerably better channel estimation and bit error performance compared to a similar time-frequency varying channel estimation approach [12].

## References

- 1.
Tarokh V, Seshadri N, Calderbank AR: Space-time codes for high data rate wireless communication: performance criterion and code construction.

*IEEE Transactions on Information Theory*1998, 44(2):744-765. 10.1109/18.661517 - 2.
Cirpan HA, Panayirci E, Doǧan H: Nondata-aided channel estimation for OFDM systems with space-frequency transmit diversity.

*IEEE Transactions on Vehicular Technology*2006, 55(2):449-457. 10.1109/TVT.2005.863427 - 3.
Shah A, Haimovich AM: Performance analysis of optimum combining in wireless communications with Rayleigh fading and cochannel interference.

*IEEE Transactions on Communications*1998, 46(4):473-479. 10.1109/26.664303 - 4.
Alamouti SM: A simple transmit diversity technique for wireless communications.

*IEEE Journal on Selected Areas in Communications*1998, 16(8):1451-1458. 10.1109/49.730453 - 5.
Paulraj AJ, Kailath T: Increasing capacity in wireless broadcast systems using distributed transmission/directional reception. US patent no. 5,345,599, 1994

- 6.
Bölcskei H, Gesbert D, Paulraj AJ: On the capacity of OFDM-based spatial multiplexing systems.

*IEEE Transactions on Communications*2002, 50(2):225-234. - 7.
Foschini GJ, Gans MJ: On limits of wireless communications in a fading environment when using multiple antennas.

*Wireless Personal Communications*1998, 6(3):311-335. 10.1023/A:1008889222784 - 8.
Goldsmith A:

*Wireless Communications*. Cambridge University Press, New York, NY, USA; 2005. - 9.
Zemen T, Mecklenbräuker CF, Wehinger J, Müller RR: Iterative joint time-variant channel estimation and multi-user detection for MC-CDMA.

*IEEE Transactions on Wireless Communications*2006, 5(6):1469-1478. - 10.
Stüber GL, Barry JR, Mclaughlin SW, Li YE, Ingram MA, Pratt TG: Broadband MIMO-OFDM wireless communications.

*Proceedings of the IEEE*2004, 92(2):271-293. 10.1109/JPROC.2003.821912 - 11.
Kang SG, Ha YM, Joo EK: A comparative investigation on channel estimation algorithms for OFDM in mobile communications.

*IEEE Transactions on Broadcasting*2003, 49(2):142-149. 10.1109/TBC.2003.810263 - 12.
Rossi PS, Müller RR: Slepian-based two-dimensional estimation of time-frequency variant MIMO-OFDM channels.

*IEEE Signal Processing Letters*2008, 15: 21-24. - 13.
Suleesathira R, Chaparro LF, Akan A: Discrete evolutionary transform for time-frequency signal analysis.

*Journal of the Franklin Institute*2000, 337(4):347-364. 10.1016/S0016-0032(00)00041-7 - 14.
Bello PA: Characterization of randomly time-variant linear channels.

*IEEE Transactions on Communications*1963, 11(4):360-393. 10.1109/TCOM.1963.1088793 - 15.
Hahm MD, Mitrovski ZI, Titlebaum EL: Deconvolution in the presence of Doppler with application to specular multipath parameter estimation.

*IEEE Transactions on Signal Processing*1997, 45(9):2203-2219. 10.1109/78.622944 - 16.
Shu F, Lee J, Wu L-N, Zhao G-L: Time-frequency channel estimation for digital amplitude modulation broadcasting systems based on OFDM.

*IEE Proceedings on Communications*2003, 150(4):259-264. 10.1049/ip-com:20030215 - 17.
Akan A, Chaparro LF: Modeling and estimation of wireless OFDM channels by using time-frequency analysis.

*Circuits, Systems, and Signal Processing*2006, 25(3):389-403. 10.1007/s00034-005-0210-z - 18.
Priestley MB:

*Non-Linear and Non-Stationary Time Series Analysis*. Academic Press, London, UK; 1988. - 19.
Melard G, Schutter AH: Contributions to evolutionary spectral theory.

*Journal of Time Series Analysis*1989, 10: 41-63. 10.1111/j.1467-9892.1989.tb00014.x - 20.
Akan A, Chaparro LF: Multi-window Gabor expansion for evolutionary spectral analysis.

*Signal Processing*1997, 63(3):249-262. 10.1016/S0165-1684(97)00161-8 - 21.
Tang Z, Cannizzaro RC, Leus G, Banelli P: Pilot-assisted time-varying channel estimation for OFDM systems.

*IEEE Transactions on Signal Processing*2007, 55(5):2226-2238.

## Acknowledgments

This work was partially supported by The Research Fund of The University of Istanbul, Project nos. 6904, 3898, and 4382.

## Author information

### Affiliations

### Corresponding author

## Rights and permissions

**Open Access**
This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (
https://creativecommons.org/licenses/by/2.0
), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## About this article

### Cite this article

Önen, E., Akan, A. & Chaparro, L. Time-Frequency Based Channel Estimation for High-Mobility OFDM Systems–Part I: MIMO Case.
*EURASIP J. Adv. Signal Process.* **2010, **549197 (2010). https://doi.org/10.1155/2010/549197

Received:

Accepted:

Published:

### Keywords

- Orthogonal Frequency Division Multiplex
- Channel Estimation
- Doppler Frequency
- Orthogonal Frequency Division Multiplex System
- Orthogonal Frequency Division Multiplex Symbol