#### Abstract

This paper introduces an effective joint scheme of channel estimation and tracking for downlink partial usage of subchannel (DL-PUSC) mode of mobile WiMAX system. Based on the pilot pattern of this particular system, some channel estimation methods including conventional interpolations and a more favorable least-squares line fitting (LSLF) technique are comparatively studied. Besides, channel estimation performance can be remarkably improved by taking advantage of channel tracking derived from the preamble symbol. System performances in terms of packet error rate (PER) and user link throughput are investigated in various channels adopted from the well-known ITU models for mobile environments. Simulation results show a significant performance enhancement when the proposed joint scheme is utilized, at least 5 dB, compared to only commonly used channel estimation approaches.

#### 1. Introduction

Wireless metropolitan area network (Wireless MAN) or worldwide interoperability for microwave access (WiMAX), which is defined in IEEE Std. 802.16d/e [1, 2], is a new technology that provides wireless access in fixed and mobile environments. Some modes in this system utilize orthogonal frequency division multiple access (OFDMA) technique as a modulation method. This technique is adopted from the powerful orthogonal frequency division multiplexing (OFDM) which effectively mitigates the impairment of the time-variant frequency selective fading channel [3, 4]. A typical OFDM system is shown in Figure 1.

DL-PUSC, as specified in [1], is one of the multiple access modes for downlink direction which is popularly used for performance analysis. This scheme divides OFDM symbol into subchannels and assigns them to users/subscribers. Each subchannel is further partitioned into groups of 14 consecutive subcarriers called clusters. Clusters of one user are not continuously connected but are pseudorandomly permuted over OFDM symbol among users so that data of different users are treated equally over the effect of fading channels. Hence, channel estimation and equalization for recovering the original signal of each user must be performed from cluster to cluster.

Pilot-based channel estimation is widely used in OFDM transmission system. By scattering known data called pilots into OFDM symbol at the transmitter, calculating channel values at pilot positions and then interpolating the whole channel values for data subcarriers at the receiver, transmitted information can be recovered. There are many techniques reported for channel estimation; some conventional methods like linear and cubic spline [5–8] interpolations are commonly used due to their low complexity for practical implementation, yet low efficiency. Other methods like transform-domain processing [9, 10] perform better but require higher computation for executing DFT/IDFT. Theoretical optimum method like minimum mean square error estimator (MMSEE) [11–15] gives best performance but is too complicate for practical realization.

Preamble is a special OFDM symbol transmitted at the beginning of transmission frame. Since it contains lots of pilots usually evenly distributed, channel estimation task for preamble symbol is quite easy and accurate. Moreover, as the radio channel is often slowly faded, by some tracking algorithm, the estimated channel can be exploited to enhance the performance of estimation for the subsequent data symbols.

In this paper, a new scheme of channel estimation for DL-PUSC system is proposed. This approach uses LSLF technique combined with channel tracking to form a joint channel estimation scheme. Comparisons between this new method and other conventional approaches such as linear and cubic spline interpolations by simulating a typical 1024-point FFT system profile in different ITU mobile channel models are given to show that a significant improvement in system performance can be achieved.

The following parts of this paper are organized as follows. Section 2 addresses the transmission structure and the signal model. Channel estimation methods such as the commonly used linear and cubic spline interpolations and the new method using LSLF technique are highlighted in Section 3. Joint scheme with tracking algorithm is introduced in Section 4. Simulation setup, results, and discussion are given in Section 5. Finally, Section 6 summarizes and concludes the paper.

#### 2. System Description

##### 2.1. Transmission Structure

A basic transmission structure is shown in Figure 2 in which a preamble symbol starts the frame and data symbols are transmitted right after. Preamble symbol is designed according to different profiles of transmission, that is, in 1024-point FFT mode, there are 284 boosted BPSK pilots, each every three subcarriers, starting from subcarrier index 86 (indexing starts from 0). Subcarriers at other positions are set to 0. Pilot values are generated by a particular Pseudo-Noise code related to IDCell and segment parameters [1].

In DL-PUSC mode, an OFDM symbol is divided into subchannels; each of those is associated to a specific user. Subchannel is further partitioned into clusters; each of which contains a group of 14 consecutive subcarriers. When transmitted, clusters of different users will be permuted among themselves; therefore, they are scattered over the OFDM symbol. The data symbol structure is shown in Figure 3.

Pilot pattern of a cluster in a DL-PUSC frame of data symbols is shown in Figure 4. Pilots are allocated at subcarrier for even symbol and at subcarrier for odd symbol.

##### 2.2. Signal Model

Assume that transmitted frame has *M* OFDM data symbols in which is the symbol at time *m* and *N* is the number of subcarriers in OFDM symbol.

At the receiver, if intersymbol interference is negligible, received signal could be derived as
where **A** is the diagonal matrix whose values are . The channel frequency response is the DFT of the time-varying multipath fading channel impulse response of which a discrete-time version can be obtained as in [9, 10]. **F** is the matrix whose entries are where *L* is the number of channel impulse response taps, and is the additive white Gaussian noise.

In order to recover from , the channel has to be estimated by exploiting the pilots which are located at predefined positions in OFDM symbols. The least-square values of channel frequency response at for pilots are obtained by where denotes the pilot position (in this particular case, ).

The goal is to estimate all the channel values at all data subcarriers from the values of so that should be as much similar as possible to .

Here, we have two kinds of pilot pattern depending on whether they belong to preamble symbol or data symbol. It is obvious that the density of pilot subcarriers in preamble symbol is higher than that in data symbol. Thus, in order to estimate the whole channel to have a reference for channel tracking, conventional method like linear interpolation can be utilized to get a good tradeoff between complexity and performance. On the other hand, for data symbol, in this particular case of DL-PUSC, the interpolation task must be performed from cluster to cluster, and because each cluster contains only 14 consecutive subcarriers, the channel on each cluster can be approximated as a “line”; this fact inspires the idea of using LSLF technique to estimate the partial channel. Therefore, a comparative study is carried out to demonstrate the superiority of this approach to other commonly used methods such as linear [5, 6] and cubic spline [7, 8] interpolations.

Due to the pilot arrangement in data symbols, it is necessary to perform a two-dimension (2D) interpolation scheme or two 1D estimations in cascade. As the number of pilots in time axis is more than that in frequency axis, it is more convenient to estimate first in time and then in frequency.

#### 3. Channel Estimation

##### 3.1. Conventional Methods

###### 3.1.1. Linear Interpolation

Figure 5 shows an example of the linear interpolation technique in time direction for an even number of OFDM symbols in which pilots of even symbols are located at the 5th and 9th locations while those of odd symbols are resided in the 1st and 13th places. Consider **
Then, interpolation in frequency direction can be evaluated as
where *k* denotes the position of channel value inside the interval of two adjacent pilots (,).

###### 3.1.2. Cubic Spline Interpolation

With this method, because we do not have enough pilots in frequency direction, interpolation in time has to be performed first. Moreover, it also requires having at least 8 OFDM symbols to have enough pilot points. Parameters of a third-ordered equation are calculated as indicated in [8] and [16]. Then, the channel values at data positions within the appropriate interval are calculated. After interpolation in time, there are enough points to carry out this task again in frequency direction.

##### 3.2. Channel Estimation Using LSLF

Since a cluster size is small, channel transfer function of that cluster can be considered a “line”. Therefore, the least-squares line that fits all the pilots is thought to be a better approximation of the ideal channel.

###### 3.2.1. Interpolation in Time

Suppose that an even number of OFDM symbols appears in frame and is indexed from 0 to . In time direction, for a cluster, there are four sub-streams according to subcarrier indexes containing channel values at pilots. Stream 1 and stream 13 have pilots at odd locations while stream 5 and stream 9 have pilots at even locations . Define the channel values at pilot positions vector and pilot position vector , in which for stream and for stream ; LSLF technique will find the pair of coefficients to form the line containing the set of points: ; so that the least-squares error is minimized. That means to find a pair of coefficients so that they minimize and so vanish the partial derivatives and . Therefore, from [17], problem turns into solving this system of equations

Channel values at all locations including data and pilots in stream 1, 5, 9, and 13 will be calculated by applying

It is important to derive the maximum number of OFDM symbols *M* so that the fitting by using a “line” is reasonable. By the fact that the fading channel will change in time with a coherent time *T _{c}* as mentioned deeper in next section, it is clear to see that the limit should be . So a rough limit range for

*M*can be and for convenient

*M*should be even number.

###### 3.2.2. Interpolation in Frequency

The same routine as in time axis can be used to interpolate in frequency axis. Now, there is a block of *M* clusters; each cluster contains 14 subcarriers in which 4 locations were estimated values from the previous task. One note is that all clusters now have “pilots” at the same indexes; hence, the complexity is less.

Again, it is crucial to judge the assumption that channel transfer function in a cluster can be viewed as line. The fact that the frequency range of 14 subcarriers or one cluster should be less than the coherent bandwidth of the fading channel [18] gives in which denotes the root mean squared delay spread of the multipath fading channel. Another factor is the bandwidth of the designed system; the performance of this method would degrade when system bandwidth is significantly broader than *B _{c}* so that a small portion as cluster could also be frequency selective.

#### 4. Joint Scheme with Channel Tracking

For preamble symbol, the least-square channel values at pilot positions are

In the case of DL-PUSC mode, pilots in preamble symbol are evenly spaced scattered, one every three subcarriers. The whole channel values can be linearly interpolated (except the two zero-guard interval regions) as where ; ;;

In fact, the channel does not stay the same over time but slowly changes; this is due to the relative movement of all the components influencing the transmission. The most impact factor is the relative speed between mobile station and base station that causes a Doppler frequency shift . The coherent time over which the channel can be viewed as unchanged is in the order of several milliseconds to hundreds of milliseconds. Hence, it is considered “slow” when comparing to an OFDM symbol time slot. As a result, after estimating the whole channel from the preamble symbol, performance of channel estimation for successive data symbols can be enhanced by a tracking algorithm [9].

Suppose that is the estimated channel value at subcarrier *n* of data symbol *m*; it can be recalculated so that some useful information from the preamble symbol can be involved to reduce the distance between it and the real channel, and thus enhance the estimation performance:
where ; .

This tracking algorithm means that the nearer the data symbol is located to the preamble symbol; the more influent it gets from the estimated channel given by the preamble and vice versa.

#### 5. Simulation Results

##### 5.1. Simulation Setup

The typical 1024-point FFT profile, whose parameters are given in Table 1, is chosen for simulation. The number of users requesting service is assumed to be 3.

Channel models are taken from ITU models for mobile environments [19], and carrier frequency is set to be 2.3 GHz: (i)Model 1: ITU_Pedestrian B (Ped.B), speed 6 km/h, and fading frequency Hz, (ii)Model 2: ITU_Vehicular A (Veh.A), speed 30 km/h with Hz, and speed 120 km/h with Hz.

These channel models are time-variant frequency selective channels in Non Line-Of-Sight mobile conditions. Their specific parameters are given in Table 2.

System performance is demonstrated as packet error rate (PER) versus signal-to-noise ratio (SNR) and user link throughput defined as
where *D* is the peak data rate given by
in which *N _{s}, N_{b}, R_{c}, *and

*T*denote the number of subcarriers assigned to a user, number of data bits in a subcarrier, channel coding rate, and the OFDM symbol time, respectively.

_{s}##### 5.2. Simulation Results and Discussion

Simulation results are shown in Figure 6 to Figure 14. The very first notice is that in all channel conditions the LSLF approach always outperforms the other two conventional methods. The improvement varies depending on which modulation mode is used. It is also very clear to see that when channel tracking comes into play, regardless of modulation schemes and channel conditions, the performance is remarkably boosted. The joint scheme of LSLF channel estimation and channel tracking appears to be the best, very robust, and not only highly surpassing other schemes but also able to reach very near to the perfect channel knowledge case.

Figures 6, 7, and 8 show PERs in different channel models. Ped.B channel has quite long delay spread, causing severely frequency-selective faded channel and limiting the performance of channel estimation, particularly in frequency axis. However, due to the slow moving speed, the channel does not change rapidly, giving some favor to estimation in time. One can notice that the effect given by channel tracking in this channel is not as strong as that in Veh.A channels. On the other hand, Veh.A channels have smaller delay spread but higher moving speed, meaning that the channel within a cluster is flat but it changes faster. The coherent time of this channel in case of speed 120 km/h is in order of several milliseconds which can degrade the system performance since it might go below the frame time. However, the LSLF still works properly and with channel tracking; at least 5-dB enhancement can be achieved.

Figures 9, 10, 11, 12, 13, and 14 show the user link throughputs in various channel conditions without and with channel tracking. Obviously, channel tracking give a noticeable improvement and the joint scheme of LSLF channel estimation with tracking significantly increases the link performance.

Another notice is that the performance improvement also depends on modulation modes. The higher modulation modes always suffer higher error in channel estimation, leading to more degradation compared to the perfect channel knowledge case whereas in lower modulation mode, for example, QPSK, the joint scheme is able to reach the ideal case.

Last but not least, it is worth to examine roughly the complexity of LSLF method for practical implementation. From equations (7), it is obvious that the LSLF method needs more computation than linear and cubic spline interpolation but it does not require any complicated process or special design structure. There is no complex operation since the in-phase and quadrature components can be treated separately whereas there are also some terms in (7) that can be reused. Therefore, the superior performance gain obtained by the joint scheme with channel tracking makes this method very promising for realization.

#### 6. Conclusions

This paper has studied a joint channel tracking and estimating scheme which is highly suitable for OFDMA DL-PUSC mode of mobile WiMAX system. System simulation with various standardized channel models for mobile environments showed impressive improvements in both PER and user link throughput. Low complexity and high performance give this joint scheme a high potential for practical implementation.