Abstract
In this article, we propose unit vectors in the high dimensional Cartesian coordinate system as the preamble, and then propose a semiblind timing synchronization and channel estimation scheme for orthogonal frequency division multiplexing (OFDM) systems. Due to the lack of useful information in the timedomain, a frequencydomain timing synchronization algorithm is proposed. The proposed semiblind approach consists of three stages. In the first stage, a coarse timing offset related to the delayed timing of the path with the maximum gain in multipath fading channels is obtained. Then, a fine time adjustment algorithm is performed to find the actual delayed timing in channels. Finally, the channel response in the frequencydomain is obtained based on the final timing estimate. Although the required number of additions in the proposed algorithm is higher than those in conventional methods, the simulation results show that the proposed approach has excellent performance of timing synchronization in several channel models at signaltonoise ratio (SNR) smaller than 6 dB. In addition, for a lowdensity paritycheck coded singleinput singleoutput OFDM system, our proposed approach has better biterrorrate performance than conventional approaches for SNR varying from 3 to 8 dB.
Keywords:
Timing synchronization; Fine time adjustment; Channel estimation; Orthogonal frequency division multiplexing (OFDM); Unit vectors; Frequencydomain processingIntroduction
Orthogonal frequency division multiplexing (OFDM) is a promising technology to support highrate wired and wireless applications due to its robustness to multipath delay spread [13]. However, in OFDM systems, synchronization errors can destroy the orthogonality among the subcarriers and result in performance degradation. Thus, timing synchronization in OFDM systems becomes much more challenging due to the increase in the amount of intercarrier interference (ICI) and intersymbol interference (ISI) [1]. Although the soft decoders employing error correction code can improve the system performance at low signaltonoise ratio (SNR), perfect timing synchronization is necessary for the decoder to operate correctly. Therefore, in order to improve the system performance, it is important to find the actual delayed timing in multipath fading channels at the receiver. In addition, channel estimation also plays a crucial role in providing the channel information to the soft decoder and compensating the signal during the demodulation process [4]. Without the knowledge of timing offset and channel information at the receiver, the system will have a poor performance during the entire data transmission.
Various synchronization techniques for orthogonal frequency division multiplexing (OFDM) systems have been developed using welldesigned preambles [514]. Although accurate timing estimation can be achieved, the bandwidth efficiency is also inevitably reduced. In order to reduce the waste of bandwidth, nondata aided synchronization algorithms based on the cyclic prefix (CP) have been proposed [15,16]. However, in some multipath fading channels with nonlineofsight (NLOS) propagation, both dataaided and nondataaided synchronization methods frequently lead to the delayed timing in channels where the delayed path has larger gain than the first path. In this case, the resulting ICI and ISI would degrade the system performance. Also, the channel coding would not perform well because of the synchronization errors. Therefore, in order to solve this problem, a fine time adjustment is needed to modify the frequently delayed timing to the actual delayed timing in channels. In [15], the proposed timing estimator performs well only for the additive white Gaussian noise (AWGN) channels. While the system operates in the multipath fading channels, the proposed algorithm exhibits significantly large fluctuation in the estimated timing offset. In [16], the modified blind timing synchronization method has a good performance in the multipath fading channels with lineofsight (LOS) propagation only when the SNR is greater than 20 dB. In [14], a welldesigned timedomain training sequence is utilized to perform joint timing synchronization and channel estimation. Although the proposed timing estimator has excellent performance at low SNR [14], the power consumption of the proposed preamble is still too large to be adopted in some lowpower wireless applications.
For wireless implantable medical devices, lowpower consumption is necessary in order to prolong the battery operating time. This article develops a semiblind timing synchronization and channel estimation algorithm based on unit vectors, and demonstrates that this algorithm is suitable for multipath fading channels with both LOS and NLOS propagation. Due to the use of unit vectors as the preamble, the power consumption of this preamble at the transmitter is approximately equal to zero. Therefore, the proposed preamble is suitable for any lowpower wireless implantable medical device. In addition, we utilize only one nonzero sample in the training sequence to perform the timing synchronization, and this training sequence definitely lacks useful information at the receiver. Compared with the existing preamblebased methods [59,11,12,14], the number of nonzero elements in the proposed training sequence is the lowest. Thus, the proposed joint approach is called a semiblind method. In this article, we first obtain a coarse timing offset using the crosscorrelation function outputs in the frequencydomain. Then, a fine time adjustment algorithm based on these outputs is applied. Finally, the channel response in the frequencydomain is obtained. Simulation results are represented to verify the effectiveness of our proposed algorithm. This article is an extended version of [13].
This article is organized as follows. Section 2. describes the system and the problem. In Section 3., the proposed semiblind timing synchronization and channel estimation algorithm is presented. Simulation results are provided in Section 4., respectively. Finally, Section 5. concludes this article.
Problem statement
System description
In this article, we consider a trainingsequencebased singleinput singleoutput OFDM system as shown in Figure 1. The training sequence is an unit vector in an Ndimensional Cartesian coordinate system, where N represents the number of subcarriers in the OFDM system. Let ^{pT}=[0…010 …0]={p(n), ∀n∈_{Ω1}} denotes the proposed training sequence, where _{Ω1}={0,1,…,N−1}, p(n)=δ(n−c), c∈{0,1,…,N−1}, the length of ^{pT} is N, and the power of ^{pT} is equal to 1/N. Consider the transmitted packet ^{sT}=[^{pT}^{xT}]={s(n),∀n∈_{Ω2}}, where ^{xT} consists of ℓ OFDM symbols, the length of ^{xT}is ℓ·(N + _{NCP}), _{NCP} denotes the length of cyclic prefix, ℓ is a positive integer, and _{Ω2}={0,1,…,ℓ·(N + _{NCP}) + N−1}. Assume that the cyclic prefix in each OFDM symbol is longer than the maximum delay spread of the channel, and the path delays in the channels are samplespaced. Therefore, the received signal at the receiver can be expressed as
where ∈is the carrier frequency offset (CFO) normalized to the OFDM subcarrier spacing, τ is the timing offset, h(k) represents the kth tap channel impulse response, K is the number of taps in the channel, and w(n) is a complex AWGN sample. After coarse frequency synchronization, the CFOcompensated received signal at the receiver is
where Δ∈ denotes the residual CFO and .
Timing synchronization in the timedomain
For any trainingsequencebased communications system, timing synchronization can be easily achieved based on a welldesigned timing metric in the timedomain. However, in this article, the proposed training sequence is a delta function with unit amplitude. Thus, if we perform the timing synchronization in the timedomain, the correlation function outputs M(d) can be expressed as follows:
where is the estimated timing offset, represents the absolute value of , _{Ω3} is the observation interval, _{Ω3}={0,1,…,D−1}, and D is the length of observation interval. If there is no residual CFO in Equation (2), we rewrite M(d) in Equation (3) as follows:
From Equation (4), it is possible that all correlation function outputs related to the channel, , are smaller than the other elements in N·M(d) at low SNR, where k∈{0,1,…,K−1}. Thus, we will have wrong timing estimates at low SNR as shown in Figure 2. In Figure 2, the delayed timing offset is 65, and c is 31. Then, a maximum correlation function output near the 97th sample is expected. However, in Figure 2, a wrong timing estimate is obtained when SNR is −5 dB.
The proposed approach
Coarse timing synchronization
In order to achieve better timing synchronization performance at low SNR, from Section 2.2, a synchronization method in the timedomain is not suitable for the proposed preamble. However, much more information in the frequencydomain can be utilized to achieve better performance in timing synchronization. Consider there are two unit vectors, _{p1}(n)=δ(n−_{c1}) and _{p2}(n)=δ(n−_{c2}), the crosscorrelation function outputs between these two unit vectors in the frequencydomain are
where ∀_{c1},_{c2}∈{0,1,…,N−1} and m represents the subcarrier index. Therefore, based on Equation (5), a frequencydomain timing synchronization scheme based on the crosscorrelation function outputs is proposed. By employing the crosscorrelation function in the frequencydomain, a timing metric for coarse timing estimation (_{τc}) is given by
where ℜ{u} and ℑ{u} represent the real part and the imaginary part of u, respectively, R(·,·) denotes the Fourier transform of the received signal , ^{b∗}(m) is the complex conjugate of b(m), b(m) denotes the absolute value of b(m), , _{Ωc} is the observation interval, _{Ωc}={0,1,…,L−1}, _{d1} is the time index, _{d1}∈{0,1,…,L−1}, R(_{d1},m) represents the value of the mth subcarrier with respect to _{d1}, U(_{d1}) is the crosscorrelation function output in the frequencydomain, and L is the length of observation interval. In addition, if there is no CFO in Equation (1), _{M1}(_{d1}) in Equation (6) can be further modified to
However, by using both real part and imaginary part of the crosscorrelation function output, more information can be utilized to obtain a better coarse timing estimate.
Assume an unit vector _{pi}(n)=δ(n−_{ci}) is transmitted over a tworay multipath fading channel (_{hi}) without AWGN, a time delay is given by τ, and the power profile of the channel is equal to {0.3, 0.7}, where _{ci}∈{0,1,…,N−1}. Therefore, the received signal is
Consider is equal to [0.3873 + 0.3873j0.5916 + 0.5916j]. Then, the received signal _{ri}(n) is
Based on Equation (6), the crosscorrelation function output (_{M1}(_{d1})) is
Thus, a coarse timing estimate (_{τc}) is
From Equation (10), although _{M1}(_{d1}) gives a maximum value when _{d1}is at the delayed timing of the path with the largest gain in multipath fading channels, the actual delayed timing cannot be obtained.
In addition, for the general channel impulse response h in Equation (1), the received training sequence is
where and E[·] is the expectation operation. Then, the corresponding timing metric _{M1}(_{d1}) is
where k∈{0,…,K−1}, , , and . From Equation (13), we can easily obtain
and
Fine time adjustment
Let us pay attention to Equations (10) and (13). In Equations (10) and (13), two crosscorrelation function outputs related to the multipath fading channel have a strong connection. The correct timing offset can be found using a simple threshold on crosscorrelation outputs. Then, by utilizing the crosscorrelation outputs at two adjacent timing indices, we can obtain the actual delayed timing in the channels. First, a sliding observation vector (SOV) v based on the coarse timing estimate is utilized to perform the fine time adjustment, where
the length of the SOV is V , and _{Ωv}={0,1,…,V−1}. If _{M1}(v(i + 1))>β·_{M1}(v(i)) and _{M1}(v(i + 2))<β·_{M1}(v(i + 1)), the final timing estimate ( ) is v(i + 1), where β is a threshold and i∈_{Ωv}. The detailed procedure of fine time adjustment is described in Algorithm 1
Algorithm 1.
Fine time adjustment.Initial Inputs: _{M1}(v(i)), v
1: fori = 0 toV−1do
2: if_{M1}(v(i + 1))>β·_{M1}(v(i))then
3: u=i + 1
4: else
5: break
6: end if
7: end for
In Algorithm 1, βis utilized to perform the fine time adjustment. Based on Equation (13), _{τc}is approximately equal to the timing index of the path with the largest gain in multipath fading channels. Therefore, the time difference between the timing index of the path with the largest gain and the timing index of the first delayed path in the channels is approximately equal to _{NI}−1, where the actual number of iterations executed in Algorithm 1 is _{NI}and _{NI}<V. Assume the second path has the largest power in the channel. If the correct delayed timing is obtained, and _{M1}(τ−1)<β·_{M1}(τ) must be satisfied. Based on these two conditions, we have
we obtain the bound of the threshold βgiven by
Table 1. The power profiles of different channel models
In general, consider the kth tap has the largest power in the channel. Then, the threshold β in the fine time adjustment can be chosen by satisfying:
where k>0. Moreover, if the first path is the path with the largest gain in the channel, the threshold can be easily set to
where the ^{k″}th tap has the secondlargest power in the channel.
Channel estimation
After the final timing estimate ( ) is found, the channel response in the frequencydomain could be obtained in a simple way. Therefore, the estimated channel response in the frequencydomain is
where is the estimated channel response on the mth subcarrier.
Simulation results
A packetbased lowdensity paritycheck (LDPC) coded singleinput singleoutput (SISO) OFDM system was used for simulations, where each codeword is encoded with code (1600,800) [17] and each packet consists of a training sequence followed by 17 random OFDM data symbols. The structure of OFDM data symbols follows the IEEE 802.11a standard defined in [9], where N=64 and _{NCP}=16. The training sequence of each packet is an unit vector with unit amplitude in the timedomain, where c=31 and the power of the training sequence is 1/64. Quaternary phaseshift keying modulation was adopted in simulations. For each packet transmission, the residual CFO was modeled as a random variable that is uniformly distributed within ±0.1 OFDM subcarrier spacing. In addition, the phase tracker based on the pilots in the frequencydomain is utilized to compensate the phase error [9].
We evaluate the proposed approach and other related schemes [7,9,15] under 6path Rayleigh channels, where the power profiles of their first four taps are described in Table 1 and represents the i + 1th tap power in the channel. A delayed timing offset (τ) is given by 65 samples. Channel Models I and II (CH I and CH II) represent multipath fading channels with NLOS propagation, and Channel Model III (CH III) is a typical multipath fading channel with LOS propagation. For CH I, the power of second tap dominates all channel taps. As for CH II, the third tap has the strongest power in the channel, and it is the worst channel model to evaluate the performance of timing synchronization in this article. Moreover, assume all channels are quasistationary during each packet transmission.
The main motivation of this article is to achieve perfect timing synchronization in very low SNR environments by using unit vectors in the N dimensional Cartesian coordinate system. In Algorithm 1, although the number of iterations is defined by V , the number of iterations actually depends on the comparison between crosscorrelation function outputs. In this article, the actual number of iterations executed in Algorithm 1 (_{NI}) is less than 3. In addition, in CH I and CH II, the thresholds β in the fine time adjustment should be less than 0.5256 and 0.5794, respectively. Therefore, based on the SNR, we employ different thresholds βdefined in fine time adjustment to achieve better performance. For SNR ≤ 0 dB, β is 0.5. As for SNR > 0 dB, β is 0.3. In [7], the timedomain training sequence is generated by a Golay complementary sequence, i.e., ±1, and the length of the timedomain training sequence is N. In addition, the actual threshold in [7] is , where η is a threshold factor and is the strongest channel tap gain estimate. The same criterion for βis applied to η. Moreover, in this article, presimulations and mathematical derivations are not required to choose the threshold for fine time adjustment [7,8,11]. Let L=200 and V=50. Therefore, the length of the interval for finetiming estimation in [7] is also set to 50. For [15], we use four concatenated cyclic prefixes to perform the timing synchronization. The corresponding results are reported in Figure 3. The perfect timing synchronization is defined as the successful acquisition of the position of the first tap in channel models.
In Figure 3, simulation results show that our proposed approach has better performance in timing synchronization at very low SNR. In CH I, the proposed approach achieves perfect timing synchronization when SNR exceeds 1 dB. As for CH II, perfect timing synchronization is achievable using the proposed algorithm when SNR = 6 dB. In CH III, the proposed approach achieves perfect timing synchronization when SNR exceeds −5 dB. Moreover, for c=63, the perfect timing synchronization is achievable at low SNR by only appending one sample to the front of the transmitted packet. The synchronization methods used in IEEE 802.11 standards lead to the delayed path with the maximum gain in channels. Also, the standardized methods are only suitable for the channels with LOS propagation. In CH III, the first tap power is approximately equal to the second tap power. Thus, higher SNR is needed to achieve perfect timing synchronization for [9]. As for [7], the reason why the scheme has poor performance of perfect timing synchronization is that AWGN affects the entire fine time adjustment process at low SNR, especially in CH I and CH II. Therefore, low SNR and wide interval for fine time adjustment significantly degrade the performance in [7,8,11]. For [15], the proposed timing estimator is only suitable in AWGN channel at high SNR. Thus, the time estimator of [15] has poor performance in all channel models. Besides the probability of perfect timing synchronization, we also evaluate the bias and root mean squared error (RMSE) of each approach. The simulation results are listed in Figures 4 and 5. In Figures 4 and 5, our proposed approach has better performance than other methods in [7,9,15]. In Figure 4, our proposed approach performs almost unbiased at any low SNR, and wide interval for fine time adjustment in [7] leads the time estimator to have negative biases. In Figure 5, the proposed approach achieves zero RMSE due to the ability to identify the first arrival path in channels.
In Figure 6, we compare the proposed approach with [7] in terms of biterrorrate (BER). For CH III, zero BER is achievable for the system using our proposed approach when SNR exceeds 6 dB, because there is no timing error to process the received signals. As for CH I and CH II, zero BER is still achievable when SNR exceeds 8 dB. In addition, the BER performance of [7] decreases slowly and still does not reach 1^{0−2} when SNR = 8 dB in CH I and CH II. As for Figure 6, we apply the synchronization methods in [9] to find the timing offset, and then utilize the perfect channel information to compensate the received signal. In Figure 7, zero BER is achievable when SNR exceeds 4 dB in CH III, because perfect timing synchronization is achieved. For CH I and CH II, there are two error floors shown in Figure 7, because timing errors degrade the system performance by introducing ISI and ICI. In addition, from Figure 7, the channel estimation in our proposed algorithm has a 2 dB loss compared to the CramérRao lower bound.
Conclusion
In this article, we have developed a semiblind timing synchronization and channel estimation scheme for OFDM systems based on unit vectors. We also developed a fine time adjustment algorithm to find the actual position of the first arrival path in LOS and NLOS channel models. Based on a simple threshold without any presimulation and theoretical derivation, the proposed fine time adjustment algorithm outperforms conventional schemes even at very low SNR. Simulation results show that there are no timing errors in our proposed time estimator when SNR exceeds 6 dB. In addition, zero BER is also achievable for a LDPC coded SISO OFDM system when SNR exceeds 8 dB.
Competing interests
The authors declare that they have no competing interests.
References

H Sari, G Karam, I Jeanclaude, Transmission techniques for digital terrestrial TV broadcasting. IEEE Commun. Mag 33, 100–109 (1995)

U Reimers, Digitalvideo broadcasting, IEEE Commun. Mag 36, 104–110 (1998)

L Cimini, J Chuang, N Sollenberger, Advanced internet services. IEEE Commun. Mag 36, 150–159 (1998). Publisher Full Text

J Montojo, L Milstein, Channel estimation for nonideal OFDM systems. IEEE Trans. Commun 58(1), 146–156 (2010)

T Schmidl, D Cox, Robust frequency and timing synchronization for OFDM. IEEE Trans. Commun 45(12), 1613–1621 (1997). Publisher Full Text

H Minn, M Zeng, V Bhargava, On timing offset estimation for OFDM systems. IEEE Commun. Lett 4(7), 242–244 (2000)

H Minn, V Bhargava, K Letaief, A robust timing and frequency synchronization for OFDM systems. IEEE Trans. Wirel. Commun 2(4), 822–839 (2003)

H Minn, V Bhargava, K Letaief, A combined timing and frequency synchronization and channel estimation for OFDM. IEEE Trans. Commun 54(3), 416–422 (2006)

IEEE std 802.11, Wireless LAN Medium Access Control (MAC) and Physical Layer (PHY) specifications

Y Mostofi, D Cox, A robust timing synchronization design in OFDM systemspart II: highmobility cases. IEEE Trans. Wirel. Commun 6(12), 4340–4348 (2007)

CL Wang, HC Wang, On joint fine time adjustment and channel estimation for OFDM systems. IEEE Trans. Wirel. Commun 8(10), 4940–4944 (2009)

H AbdzadehZiabari, M Shayesteh, Robust timing and frequency synchronization for OFDM systems. IEEE Trans. Veh. Technol 60(8), 3646–3656 (2011)

TL Kung, K Parhi, Frequency domain symbol synchronization for OFDM systems. IEEE Conference on Electro/Information Technology, 1–5 (2011))

TL Kung, K Parhi, Optimized joint timing synchronization and channel estimation for OFDM systems. IEEE Wirel. Commun. Lett 1(3), 149–152 (2012)

J Beek, M Sandell, P Börjesson, ML estimation of time and frequency offset in OFDM systems. IEEE Trans. Signal Process 45(7), 1800–1805 (1997). Publisher Full Text

S Ma, X Pan, GH Yang, TS Ng, Blind symbol synchronization based on cyclic prefix for OFDM systems. IEEE Trans. Veh. Technol 58(4), 1746–1751 (2009)

R Neal, (http://www), . cs.utoronto.ca/{∼}radford/ldpc.software.html webcite