#### Abstract

A simple low-complexity algorithm is proposed for blind feedforward symbol timing estimation in time-varying flat-fading channels, assuming constant modulus signalling. Its important properties are brought out using asymptotic analysis. It is shown that at high Signal-to-Noise Ratios (SNRs), the algorithm results in a timing estimate with negligible bias. It also has a Mean-Square-Error (MSE) performance that is superior for short data records and small roll-off factors to that of the Square-Law-(SQL-) and the CycloStationarity-(CS-) based algorithms available to date. Computer simulations are presented to demonstrate effectiveness of the algorithm.

#### 1. Introduction

All communication receivers require accurate symbol timing recovery to provide high data rates. In dynamic environments, use of training sequences for timing estimation is impractical, and blind algorithms [1–7] are preferred. Of these, the SQL algorithm of [4] (which performs poorly for small roll-off factors), the logarithmic nonlinearity-based algorithm (LOGN) of [2], and the recently proposed Constant-Modulus (CM) algorithm [1] are blind feedforward algorithms with low complexity. Further, the complexity involved in the search can be avoided by an approximate Fourier Series (FS) method [1]. In many applications like Wireless Sensor Networks (WSNs), algorithms with* low-complexity* are often desired.

Most blind feedforward timing estimation algorithms proposed to date assume a *deterministic channel* and cannot be used for time-varying flat-fading channels assumed in this paper. This problem is of practical importance these days because of increased mobility. The SQL algorithm is an exception and can be applied to such channels, though its properties in such channels have not been investigated so far. The constant modulus algorithm of [1] can indeed be applied in quasistatic fading channels (though it was proposed for deterministic channels). However, it cannot be applied to time-varying fading channels. The Conditional Maximum Likelihood estimator of [6] is applicable to only to static flat fading channels, and not to time-varying fading channels. This is also true of the algorithm developed in [5]. Despite its practical importance, there has been little research effort to develop blind feedforward algorithms for such links. The* only* exception is the cyclostationary-based algorithm [7], that is specifically designed to work with time-varying fading channels. However, it performs quite poorly with low excess bandwidth signals commonly used today, and for small data records. In this paper we propose a novel low-complexity CM-based blind algorithm for estimation of the timing and demonstrate its effectiveness in such links for small excess bandwith (small roll-off factors), and small data records. We study its asymptotic properties.

#### 2. Problem Formulation

Let denote the i.i.d complex CM () transmitted symbols. In a typical communications receiver, the baseband received signal (after the receiver matched filter) is given by where is the gain, is the (time-varying) fading channel, is the convolution of the transmitter's signaling pulse and the impulse response of the receiver filter, is the symbol interval, is the delay that is to be estimated, and the Gaussian noise after the receiver matched filter. We assume the following.

(**AS1**): the basic pulse used is Nyquist so that and for .

The received signal (after receiver match filtering) sampled with a time offset of is given by

where , is white Gaussian because of (**AS1**), and is the weighted sum of adjacent symbols that arises due to sampling at the improper instant ( denotes the error in timing). Clearly, is approximately (circular) Gaussian by virtue of the central limit theorem. Notice that takes the maximum value of when , or (in which case ). Notice that when the excess bandwidth is large, the variance of is small (the impulse response decays faster for large roll-off factors, as is well known). Also, when the roll-off factor is small, increases sharply with for small . We will cosider the following.

**AS2)**: the channel does not change significantly within two symbol durations ().

With slow channel fading, this is a reasonable assumption and is commonly made.

With the popular Jakes model, , where denotes the doppler frequency, and denotes the zeroth-order Bessel function of the first kind. Assumption (**AS2**) holds when , which is the case with slow fading. Using computer simulations, we establish later in this paper that the algorithm performs well for slow to medium fading channels, but not in the case when the channel changes significantly from one symbol interval to the next. This is not a serious limitation however, and the proposed algorithm will find wide applicability in practical links. To the best knowledge of the author, this is the only work other than [7] that specifically addresses the problem of timing estimation in time-varying fading channels. As noted already, all previous works have assumed deterministic or static fading channels.

The SQL estimator uses the cost function
where is the observation interval. The recent CM algorithm [1] minimizes the cost function
where . In [1], only determinsitic channels were considered. The algorithm is not applicable for time-varying fading channels. Since is CS, autocorrelation is periodic, with the FS coefficient at cycle frequencies given by . It has been shown in [7] (CS algorithm) that the frequency offset can be estimated using the phase of , and the timing estimate subsequently as the phase difference between and at various . This requires the cyclic correlations to converge, *and* the roll-off factor to be large so that the nonzero cycle frequency components are significant. In addition, the algorithm first estimates the frequency offset, and this estimate is used to obtain an estimate of the timing in a two-step procedure. The proposed algorithm on the other hand estimates the timing independent of the CFO.

#### 3. Estimation Algorithm

The following cost function is maximized by the timing estimation algorithm proposed in this paper: In what follows, we briefly study the properties of an estimator based on the above cost function. While a finite sample analysis is analytically difficult, significant insight into properties of the estimator can be obtained from an asymptotic analysis. Denote by . It can be seen that the squared magnitude of the received signal is given by Insight into the asymptotic performance of the squared-law algorithm can be obtained from the above expression. For large , in the absence of noise (where since is assumed to be zero mean). This explains the major failing of the square law algorithm. Clearly, when the excess bandwidth is small, the peak of the cost function will occur at the correct point. However, when the excess bandwidth is large, the baud-rate sampled pulse is aliased. This results in bias for large excess bandwidths. However, is small for large excess bandwidth so that the Mean-Squared Error (MSE) in the estimate is small. For the new cost function, it can be seen that where is the variance of the noise, and the higher-order terms involving and the noise are ignored (assuming high SNR). The second term in the first line of the above equation has a factor of , and this ensures good MSE performance of the algorithm.We assume a small offset , exploit the Gaussian nature of , and write the cost function as First consider the case when there is no noise so that (so that only the first two terms in the above expression are relevant). Notice that the cost function is zero when . To establish that the unbiased nature of the estimator, we only need to verify that it is not greater than zero for other (in other words, it is negative for other ). Since is circular Gaussian as noted already, . From the Cauchy-Schwartz inequality, . Using this, we have . Also, in the range , or , so that . This implies that the second term is much larger in magnitude than the first, so that the cost function is maximum only when , and there is no bias. This implies that a search-based algorithm will always estimate the timing without bias. A careful look at the cost function also reveals that the cost always decreases with an increase in for . This is true because of the fact that the first term is negligible as compared to the second term, and the product of and in the second term is zero at and continuously decreases for in the range. This once again is because of the fact that for all in the range. For this reason, the commonly used first-order method [1] can be used with this algorithm, thereby reducing the complexity, though with loss in performance.

Now consider the case when there is additive noise (). Neglecting the first term involving higher-order products of , we get where (the symbols are of unit variance, and the pulse of unit energy). The above cost function is negative of the product of two terms, one that falls with , and another (within the square brackets) that rises quite sharply with (with a minimum of ), especially at high SNRs. It is easy to verify that the minimum is no longer at , and there is always a bias in the estimation. However, this bias is clearly small because of the high slope of the term in the square brackets. Also, as SNR increases, this bias decreases.

#### 4. Performance Analysis

Performance of proposed algorithm is compared to that of the SQL and CS algorithms (using an oversampling rate of ) in this section. The proposed and SQL algorithms are assumed to use the FS approximation to avoid a search [1]. In all the simulations, the true delay is assumed to be a random variable uniformly distributed between and . Observation windows of length and are used. ( or ). A time-varying flat fading channel with is assumed unless stated otherwise.

Figure 1 depicts the MSE in the estimation of of all three algorithms for various SNRs assuming a roll-off factor of . It can be seen that the proposed algorithm yields performance that is far superior to the other two in all cases, and performance of the CS algorithm is unacceptable for these data record sizes and for the small roll-off factor used (it requires at least symbols, and percent excess bandwidth to yield comparable performance). Figure 2 shows that the proposed estimate has a small bias. Further, this bias decreases with an increase in length of the observation window. Figure 3 depicts the MSE versus the roll-off factor for an SNR of for the proposed and the SQL algorithm, assuming a time-varying fading channel with . We note that the proposed algorithm is more insensitive to roll-off factor than the SQL algorithm, which is an important advantage.

#### 5. Conclusions

A simple low complexity blind feedforward timing estimation algorithm is proposed for time-varying fading channels that exploits the constant modulus property of the transmitted symbols and is advantageous as compared to existing algorithms.

#### Acknowledgment

The work was funded by DST Project no. SR/S3/EECE/031/2008.