Abstract
Due to the constraint of cost and size for mobile wireless communication terminals, many orthogonal frequency division multiplexing (OFDM)based systems required the same crystal driving the sampling and the channel frequencies, which leads to the challenge of a more comprehensive sampling clock synchronization scheme needed. In this article, the effect of sampling clock error on the system performance was analyzed by dividing it into sampling clock frequency offset (SFO) and sampling timing error (STE) firstly. After that, we proposed a twostage scheme of sampling clock synchronization based on theoretical derivation: the preliminary SFO was jointly acquired with the carrier frequency offset by using the improved preambleaided algorithm firstly, the timing drift resulted from residual SFO and STE was tracked based on a phase looped lock in the second stage. The deviation properties of the estimation were achieved theoretically, which reveals that both the estimating variances of SFO and timing drift are in inverse proportion to signaltonoise ratio and grow linearly by the number of total subcarrier. The results of the simulation show that the proposed synchronization scheme can introduce preferable tracking and robust synchronizing performance for this kind of OFDMbased system.
Keywords:
System synchronization; Sampling clock error; Joint estimation; PLL; Timing drift1. Introduction
Underwater acoustic channels are characterized by severe bandwidth limitations, long intersymbol interference (ISI) spans, and large Doppler spreads, which lead to significant challengers for reliable communications. Owning to the advantages of high spectral efficiency and robustness against channel fading, orthogonal frequency division multiplexing (OFDM) has been applied in modern communication system extensively such as WLAN, WiMAX, LTE, and DVBT. OFDM is also a good alternative transmission scheme for underwater communication that both remedies the problem of ISI and provides low complexity solutions. However, OFDM is sensitive to intercarrier interference (ICI) caused by channel variations. Underwater channels vary fast due to the large ratio of the platform motion relative to the sound propagation speed. Even with stationary transmitters and receivers, significant ICI could still exist due to wave action and water motion. One of the solutions for eliminating ICI is to estimate and correct the frequency shift adaptively in the synchronization procedure. Therefore, the synchronization scheme for OFDM system used in underwater communications is a critical issue [1].
The sampling clock synchronization in an OFDM system is to mitigate the sampling clock errors due to the mismatch of the crystal oscillators between the transmitter and the receiver. The sampling clock error can essentially be divided into two parts: sampling clock frequency offset (SFO) and sampling timing error (STE). SFO means the offset of sampling frequency between transmitter and receiver. STE implies that the sampling does not align at the central of the samples, which was also named sampling clock phase offset in some literatures [2,3]. The sampling clock error will cause ICI, and a drift in the symbol timing and further worsen ISI [3,4]. The effects of SFO on the system performance are analyzing in terms of BER degradation and ISI in [4,5], respectively.
In the aspect of sampling clock synchronization, there are two different kinds of methods: synchronous sampling and asynchronous sampling [6]. Synchronous sampling methods [7,8] have large timing fluctuation due to highlevel phase noise when compared with asynchronous sampling methods [5,9,10], but the asynchronous sampling methods traditionally needs the interpolation in time domain, which is computationally complex and needs more processing time. Therefore, it cannot be applied to the area where there is a strict requirement on estimation delay time. Moreover, asynchronous sampling method is more sensitive to carrier frequency offset (CFO). A combined SFO and CFO estimation algorithm is presented in [11,12]. However, it has such prerequisite that the timing offset and the initial CFO should been corrected ideally.
In many wireless communication systems, such as WiMAX, there are the constraints of low cost and miniaturization for mobile communication device, and the specification requires that the same crystal must be used to drive the sampling and the channel frequency, which is adopted for many personal handheld terminals (e.g., Smart phone) especially. However, the same reference clock property introduces a new challenge for the design of joint carrier and sampling clock synchronization, where the frequency offsets have relationship with the sampling frequency offset besides the Doppler shift. It is seen that most research on joint sampling and frequency synchronization [1115] have only considered the sampling clock frequency synchronization but ignored the STE. In addition, they view the synchronization of sampling clock frequency all at once whereas neglecting the permanent drift of sampling clock frequency.
In this article, an asynchronous scheme for sampling clock synchronization in an OFDM receiver system with the same driving clock source is proposed and analyzed: First, the preliminary SFO is jointly acquired with carrier frequency offset by taking the benefit of the same crystal for both sampling and channel frequencies. Second, by introducing a phase looped lock (PLL) with dynamic adjustable parameters, the timing drift resulted from the residual SFO and STE is tracked and periodically corrected. The proposed scheme has not used the interpolation and consequently lowered the computation complexity and processing time consumption. The requirements of ideal timing offset and the initial CFO are not needed previously because the estimation and correction of SFO was jointly deal with CFO and timing error simultaneously.
This article is organized as follows. First, we, respectively, analyzed the effect of SFO and STE on the performance of OFDM system in Section 2. After that, we have derived the theoretical solution for the estimation of SFO and timing error in Section 3. Based on the derivation result, a practical estimation and correction of sampling clock error scheme is proposed in Section 4. The simulation results and discussion are given in Section 5. Section 6 concludes the article briefly.
2. Analysis of SFO and STE effects on system performance
This section presents the system model and then emphasizes on the analysis of the effect of SFO and STE on the system performance.
The OFDM system model and signal flow are illustrated in Figure 1. It needs to point out the same crystal drives both the sampling and carrier frequency mixing here. As the output of D/A converting, the transmitted OFDM signal during one OFDM symbol duration can be expressed as
where A_{k} denotes the OFDM symbol on the carrier k, N is the number of total carrier, and T_{s} is the sampling interval for OFDM symbol.
Figure 1. OFDM system model with the unified driving crystal.
In the following analysis, for the simplicity we assume the transfer function of channel satisfies H(f) = 1. The s(t) successively passes through the procedures: upper converting onto carrier frequency, additive white Gaussian noise channel, downconverting, then can be written as
In (2), Δf_{c} is the carrier frequency offset between receiver and transmitter, which is caused by Doppler shift and the mismatch between the clock of receiver and transmitter.
Here, we define the A/D sampling interval offset as ΔT = T_{r} – T_{s} (T_{r} is the sampling interval of A/D converter), and the STE as Δt, Δt ≤ T_{s}/2. The SFO obviously equals to Δf_{s} = 1/T_{r} – 1/T_{s}. During an OFDM symbol, the sampling time of A/D converter at receiver is t = nT_{r} – Δt, n = 0,1,…,N – 1. The sampled signal can be expressed as
and is the sample of channel noise.
Due to the same driving crystal for sampling and channel frequencies, the relationship between SFO and the carrier frequency offset can be denote as
where f_{s} equals to 1/T_{s}, ε_{c} is the relative carrier frequency offset and defined as the ratio of Δf_{s} to subcarrier spacing. ε_{s} is SFO normalized by f_{s}. In (4), the effect of Doppler shift on carrier frequency offset has been neglected.
2.1. STE
We first analyze the effect of STE on the received signal. For the convenience of the following derivation, we assume Δf_{c} = 0 and the symbol timing synchronization is already completed. Formula (3) is changed into
where n = 0,1,…,N – 1. Let . The signal is then transformed into frequency domain
where is the interference portion arisen from channel noise. By using the summation of geometric sequence, (6) is changed to be
The R(k) can be further divided into two parts in terms of m = k and m ≠ k:
where
When the effect of STE is only concerned, we can let sampling interval offset ΔT = 0. According to that sin x/x approach to 0 with x → 0, the first term in (8) comes to A_{k}exp{−j2πkΔt/(NT_{s})}, and the second term in (8) equals to 0. Consequently, R(k) is expressed as
From (9), it can be seen that the STE will not affect the amplitude of the useful signal, neither introduce ICI. However, the received signal has a phase rotation due to STE, and the angle of phase rotation increases linearly with subcarrier index k. When considering Δt ≤ T_{s}/2, the phase rotation effect is slighter than the effect of symbol timing error. Both the phase rotation from STE and symbol timing error can be corrected in synchronization or channel equalization procedure.
2.2. SFO
In the same way, let Δt = 0 when we analyze the effect of SFO on the received signal, formula (8) can be rewritten as
The first term in (10) denotes the useful component, which is attenuated and rotated due to the SFO. The second term represents the interference to the kth subcarrier from other subcarrier, the ICI destroys the orthogonality of OFDM signal and degrade the signaltonoise ratio (SNR) consequently. The third term is the noise.
For the aims of analyzing the effect of sampling clock synchronization detail, we will extend the above results to the case of successive OFDM symbols. The nth sample of the lth OFDM symbol can be deduced from (3)
where , and the sampling time t = T_{r}{(N + N_{cp})l + n}.
Similarly, the received OFDM signal in frequency domain is
where
and
In (12), exp{−2πk(1 + ΔT/T_{s})(N + N_{cp})l} in R_{1}(l,k) denotes the phase rotation of the kth subcarrier, which have relation with both subcarrier index and OFDM symbol index. Hence, if the SFO is not compensated during several accumulated symbol periods, the degree of phase rotation will increase in the symbol sequence. When the sampling interval offset satisfies ΔT < 0, A/D converter will produce an extra sample for every –T_{s}/ΔT samples. When the relative sampling interval offset satisfies ΔT > 0, A/D converter will lost a sample for every T_{s}/ΔT samples, which means the ISI will arise at a later time.
3. Design and theoretical analysis for sampling clock error estimation
Based on the analysis of the effect of sampling clock offset, this section will elaborate on how to achieve the accurate estimation of SFO and STE firstly, then analyze the deviation property of estimation result. Next, the timing drift caused by residual SFO and STE is derived together with performance analysis. On the basis of this, we propose a complete sampling clock synchronization scheme that composed by acquisition stage and tracking stage to realize the estimation and correction of sampling clock error.
3.1. Design of training sequence
The preambleaided estimation is adopted in this article. The preamble is the first symbol of transmission frame, of which the subcarriers are modulated using a boosted BPSK modulation with a specific pseudo noise (PN) code.
Figure 2 gives the structure of preamble sequence in frequency domain. Only each third of subcarriers (M = N/3) are modulated with PN sequence, the remainder subcarriers are filled with zero. It can be proved that this structure of the preamble in frequency domain means three repletion parts in time domain:
where we assume M = N/3 is integer, C(k) is the PN sequence.
Figure 2. The structure of preamble sequence in frequency domain.
The frame period is defined as L_{s} in sample at transmitter, it means the preamble symbol will appear for every L_{s} samples. The measuring of the frame period at receiver is used in the estimation of timing error in the tracking stage.
3.2. Sampling frequency offset estimation
From (4), the estimation of sampling frequency offset can be translated into the estimation of carrier frequency offset, and the result of carrier frequency estimation can be used to compensate both the sampling frequency offset and carrier frequency offset.
In general, the relative carrier frequency offset is divided into fractional part and integer part and estimated separately
where ε_{F} represents the fractional carrier frequency offset (FCFO) and ε_{I} is the integer carrier frequency offset(ICFO).
3.2.1. Fractional carrier frequency offset
Consider the preamble in time domain there are three identical parts, except a phase shift between the adjacent parts caused by the carrier frequency offset. Therefore, estimation of FCFO is based on the idea: if the conjugate of the first part is multiplied by the second part, the effect of channel can be eliminated, and the result is exactly the phase shift.
When neglecting the cyclic prefix, the received OFDM symbol in time domain can be expressed as
and n = 0,1,…,N – 1. C(k) is the preamble sequence in time domain, and the H(k) is the transfer function on the kth subcarrier. Here, we use ε_{F} to represent the overall carrier frequency offset. It will be proved that the estimation by the method actually is the fractional frequency offset later.
Let define the P(M) as the correlation function of r(n)
By substituting
From (17), it can be found frequency offset is equal to
The estimation range of frequency offset by (18) is (−3π/2, 3π/2), which means that the offset beyond the range cannot to be calculated using the aforementioned method. Therefore, the fractional carrier frequency offset is mainly achieved here.
The derivation in (16) has not considered the channel noise. If the channel noise exists, then the received signal in (15) will be rewritten as
w_{n} is the complex additive white Gaussian noise with the zero mean and the variance σ_{n}^{2}. Similarly in [14], it shows the maximumlikelihood estimation (MLE) of the fractional frequency offset is
When using the method in [14], the variance of ε_{F} is
The ratio of signal powertonoise power (SNR) on an OFDM symbol is defined as SNR = R_{h}(0)A/σ_{n}^{2}. Actually, because the fractional frequency offset is the MLE by (20), it states that the CramerRao bounds are almost met by the estimation with high SNR in [14]. Here, the phase shift is calculated over all of the subcarriers in (16), which usually can satisfy the condition of the high average effective SNR.
3.2.2. Integer carrier frequency offset
The integer carrier frequency offset should be estimated in frequency domain for accurate acquisition. We still achieve the estimation by taking the use of the property of PN sequence and preamble structure.
First, it is known that the PN sequence used for the preamble have the autocorrelation property
We rewrite the value of the kth subcarrier value in received signal in frequency domain as
Here W(k) is the FFT conversion of w(n), which can be proved that W(k) is still a Gaussian variable with the zero mean and variance σ_{n}^{2}. We use C(k + ε_{I}) to denote that the signal has shift ε_{I} subcarriers compared to local preamble sequence.
When shifting the local preamble sequence with i subcarriers, then multiplying with R(k):
As to the second term in (24), since the C(k + i) is a BPSK modulated symbol with value +1 or −1, the product will still hold as the same property of Gaussian variable W(k). Especially, when the shift number i equals to the existed frequency offset ε_{I}, the (24) will be
Here, we have used the result C(k)^{2} = 1 and eliminated the C(k) element in the result.
Next, the correlation function of Y(k) is defined as
Based on (26), we define a metric of integer frequency offset estimation as
It is not difficult to deduce that the value of R_{y}(3) on the condition of i = ε_{I} is more large than that while i ≠ ε_{I}. It can be explained as follows. When i ≠ ε_{I} the first term of Y(k) in (24) will be the product of H(k) and random +1/–1, which produces the counteract within the sum of R_{y}(3). For R_{y}(0), it actually represents the total power of signal and noise, which is used in (27) for normalizing the M(i).
Figure 3 shows the metric M(i) at different values of shift i, while the SNR is 0 dB. The received signal is assumed without carrier frequency offset. It is showed that the maximum metric value is achieved at i = 0 as the predicting, which means the maximum metric appears when there are no subcarrier offset between the received preamble sequence and local sequence.
Figure 3. M( i ) versus shift I, no frequency offset.
Therefore, the estimation of integer carrier frequency offset is to find the value of i when satisfying
The following will elaborate on the statistical property of M(i). First, we rewrite the (26) while i = ε_{I}
Each term in (29) can be viewed as the arithmetic mean for discrete sequence, when the accumulate number N is large enough, the arithmetic mean is approximately the statistical expectation. Since the result of linear operation on Gaussian random variable is still Gaussian, the second and third terms in (29) are equal to the mean of Gaussian variable, which both come out zero. The fourth term is the autocorrelation of white Gaussian noise, whose result is an impulse function of τ. Then, we simplify (29) into
where
From (30), we have
When the coherence bandwidth of channel is greater than the range of three subcarriers, the transfer function on the adjacent three subcarriers is viewed correlatively; therefore, there is the approximation
Regarding the transmitted preamble sequence, the power of signal can be calculated
Taking (31)–(34) together, we have
It can be deduced that the metric approaches to 1 as SNR increases. The metric value is greater than 0.9 when SNR is above 5 dB, which could be used for verifying the result of (28).
3.2.3. SFO
Once the estimation of fractional carrier frequency offset and integer carrier frequency offset are achieved, the sampling clock frequency can be calculated using the linear relationship with the carrier frequency offset as (4).
The method in (36) ignores the existence of Doppler shift in the CFO, so the estimation of SFO is not strict. However, we will deal with amend the estimation in tracking stage.
By now we have already achieved the joint estimation of carrier frequency offset and SFO, but the correction for them should be operated independently. The FCFO offset should be corrected in time domain before the estimation of ICFO, and the ICFO is estimated in frequency domain then fed back into time domain to correct the preamble in next frame. With consideration of the effect of SFO given in Section 2.3, the estimation of SFO is translated into timing error and corrected in frequency domain.
3.3. Timing drift estimation
Except residual sampling clock error, the fractional symbol timing error that may be ignored also contributes to the timing drift, which can be considered as a uniform adjustable variable. This adjusting variable is derived in frequency domain and then fed back into time domain to adjust digital oscillator, guaranteeing the synchronization stability. First, we focus on the estimation of timing drift.
We assume that the received signal compensated with the estimation ε_{s} still has the timing drift denoted by Δn sampling interval (note it may be decimals), that is
In order to eliminate C(k) in the first term of (37), we similarly define
When τ = 3 the correlation of Y(k) is
When neglecting the noise, (39) only reserves the first term:
if the coherence bandwidth of channel is greater than the range of three subcarriers, the transfer functions on the adjacent three subcarriers are viewed relatively, the phase of R_{h}(3) is approach to 0. Therefore, the timing drift can be estimated as
The ML estimation Δn of is
In the following, we will investigate the performance of the estimator by through finding the statistical distribution of timing drift Δn and its variance.
In (39), since the conjugate of complex Gaussian variable W(k) is still a complex Gaussian with the same expectation and variance, the second and third terms within the summation are independent Gaussian variable. Based on the central limit theorem, the summation of the two terms in terms of k is complex Gaussian process with distribution N(0, 2σ_{s}^{2}σ_{n}^{2}), which is expressed as N_{e} = N_{c} + jN_{s}. The imaginary N_{c} and real part N_{s} are independent Gaussian process with zero mean and variance σ_{s}^{2}σ_{n}^{2}. The summation of the last term is approximately the autocorrelation of Gaussian noise and turn to be 0 while τ ≠ 0. Therefore, R_{y}(3) in (39) can be rewritten in real and imaginary parts separately as
and define
It is clear that A_{c} and A_{s} are subject to Gaussian distribution according to the property of Gaussian process, and their expectation and variance are
The joint probability density function for A_{c} and A_{s} is
Let define
It can be seen if we find the statistical property of φ, we will get the distribution of estimation of Δn. By using Jacobian determinant
The joint PDF of G and φ can be calculated as
Thereby we can get the conditional PDF for φ, especially in the condition of high SNR and while φ is small value, the conditional PDF is
It shows the φ is a Gaussian variable with zero mean and variance σ_{n}^{2}/σ_{s}^{2}. From (54) we can get the variance of Δn
The result reveals that the variance of timing drift is in inverse proportion to SNR and grows linearly by square of the number of total subcarrier, i.e., the timing drift will be obvious in the condition of either bad channel quality or large subcarrier number.
4. PLLbased twostage sampling clock estimation and correction scheme
The proposed sampling clock synchronization scheme composed by acquisition stage and tracking stage is illustrated in Figure 4. The acquisition procedure focuses on the preliminary SFO (carrier frequency offset jointly), of which the fractional part and integer part are handled sequentially. Once the acquisition is completed, it will transfer to tracking stage and no longer return unless the restart of synchronization. In tracking stage, the timing drift is estimated and compensated for each frame by using a PLL. Once locked, the followed samples after correction can be further demodulated.
Figure 4. Practical sampling clock synchronization scheme.
4.1. Acquisition stage
The processing of acquisition stage includes estimation and correction of FCFO, estimation of ICFO, and correcting SFO and CFO at last.
The estimation and correction of fractional carrier frequency offset is implemented in time domain. Before that, it is assumed that a coarse symbol timing synchronization is achieved, and then the start position of a frame in sample is already determined with the toleration of ±1/2 cyclic prefix length. In our implementation, the symbol timing method in [13] by exploiting the correlation property of repetition and cyclic prefix is conveniently adopted here.
As the first step of overall synchronization process, it need to extract one preamble sequence (N samples) from the samples of A/D output based on the result of symbol timing synchronization. Let use P_{s} to denote the starting position of the current frame in received samples. When considering the multipath delay of wireless channel, the actual position should be P_{s} + ρN_{cp}, and 0 ≤ ρ ≤ 1. The optimum ρ is set with respect to the average multipath delay, when ρ is too greater, the delayed copy of last symbol will result in interference on the current symbol. When ρ is less than the optimum value, it will result in the ICI to next symbol. Certainly, the arbitrary timing error can be compensated in tracking stage.
Second step, by using the method given in Section 3.2.1, we can get the estimation of fractional carrier frequency offset ε_{F}, correcting is operated on the extract N samples of the preamble
Third step, the N corrected samples are transferred into frequency domain through FFT operation, R(k), k = 0,…,N – 1. According to the method mentioned in Section 3.2.2, the shift i can be found from the maximum metric M(i), which is the estimation of integer carrier frequency offset. Another modified method is to find both the maximum M(i_{1}) and the secondary M(i_{2}). Once the ratio M(i_{1})/M(i_{2}) exceeds a specific threshold (e.g., 3 in practice) we will take i_{1} as the estimation of integer carrier frequency offset. The advantage of the latter method is able to avoid the impulsive interference.
Fourth step, as the last step of acquisition stage, is to correct integer carrier frequency offset and the SFO. The integer carrier frequency offset estimated in frequency domain should be fed back to time domain to correct
which is similarly with the way of correction for fractional carrier frequency offset, it should be pointed out r'(n) is the time domain preamble after correction based on (56).
In order to correct the derived SFO in frequency domain, the impact stemmed from ε_{s} should be converted into timing error. While the frame length at transmitter in sample is L_{s}, the estimated frame length at receiver as L_{r}, then we can give the relationship between the frame length and sampling clock frequency
and the timing error caused by SFO can be derived as
In the process of extracting the preamble sequence, another two types of timing error also have been introduced. Actually, the starting position of the current frame is determined by the starting position of the last frame and the estimated frame length
Base on (60), the achieved P_{s} may include the fractional part. The integer part of P_{s} is used to extract the preamble sequence, while the fractional part should be compensated
Another type of timing error is caused by arbitrary delaying the starting position inside cyclic prefix, which is expressed as
The three timing errors are corrected together on frequency domain samples
where l is the OFDM symbol index in one frame and the preamble is the first symbol with l = 1. When the tracking stage is not in the locked state, the preamble is only needed to correct as (63) because the other symbols will not be demodulated currently.
4.2. Tracking stage
Once the processing of acquisition stage finished, the system automatically moves on to tracking stage. In tracking stage, a PLL is used for maintaining the receiver stays time and frequency locked to transmitter. Using the relationship between frequency and timing errors, the PLL initializes the expected frame period. This loop runs from acquisition onwards, adjusting frame timing each frame.
Fifth step, following the acquisition stage, the timing drift is recognized by the method in (42). The estimation of timing drift is first used to update the frame period and the starting position of next frame. The compensation scale in one frame should not be too large in case fluctuation happens, the timing drift is used in an alpha (firstorder) PLL to track the starting cycle timing of the next preamble, and a beta (secondorder) PLL to adjust the local measure of frame period in samples. Therefore, the frame period and starting position of next frame are updated
The alpha and beta PLL filter values are dynamically adjusted according to the phase of tracking stage, which will be demonstrated in Section 5.
Furthermore, since the same crystal clock used for channel frequency and sampling, the timing drift corresponds to the offset of carrier frequency. The carrier frequency offset can be updated
The correction will operate on the samples of the next frame based on the updated parameters. The correction process is same as that on the fourth step, except for Δn_{1} on the fourth step should be replaced by Δn.
Sixth step, after tracking the timing drift for several frames, the PLL will transfer into the steady state to allow the corrected data to be further demodulated. The lock metric is defined to indicate whether the PLL is locked or not
where M_{PLL,last} represents the lock metric calculated in last frame, v is the forgetting factor (0 < v ≤ 1). If M_{PLL} exceeds a predefined upper threshold PLL is regarded as locked, and if M_{PLL} is less than a predefined lower threshold it is regarded as losing lock of PLL, then the synchronization system will restart and return back to acquisition stage again. The function of the lock metric is essentially the timing drift passing through a low pass filter with the transfer function
which aims at getting a steady value of lock metric.
The overall sampling clock synchronization scheme is practically described with the aforementioned six steps processing. It is especially pointed out, in this scheme, that the synchronization of both carrier frequency and sampling clock is recovered thoroughly.
5. Simulation and discussion
The simulation is implemented to evaluate the performance of proposed sampling clock synchronization in AWGN channel and multipath fading channel (ITU VA channel [16]), the multipath channel has been developed with Doppler spread of 38.9 kHz. The parameters of OFDM transmission system are as follows: frame length 5 ms, sampling clock 22.4 MHz (8 times interpolation), subcarrier spacing of 10.94 kHz with 256 subcarriers. The OFDM subcarriers are modulated by QPSK, ReedSolomon, and circular convolution code are used as the outer and inner forwarding error correction scheme with code rate 1/2.
In the simulation process, the forgetting factor is set to 0.125. Four groups of PLL parameters are used for specified number of frames sequentially, and the alpha and beta filter values reduce as the group index increasing: the first group with α = 0.1, β = 0, the second group with α = 0.1, β = 0.0025, the third group with α = 0.05, β = 6.25 × 10^{–4} and the fourth group with α = 0.01, β = 2.5 × 10^{–5}. As the alpha and beta filter values reduce the PLL loses its ability to make large adjustments to the timing, but tracks the training symbol to a great resolution.
Figure 5 depicts the tracking process of SFO in AWGN and multipath fading channel. The simulation is performed at SNR 10 dB and normalized SFO 73.728 ppm. In the first 100 number of frames, the PLL is initialized using the first group parameters, of which the beta filter value is set to zero for considering that the timing error is mainly caused by the inaccurate of symbol timing at the beginning of synchronization, so sampling clock frequency remains major error. The next two group parameters are used for 100 number of frame sequentially, but during the third 100 frames if the lock metric exceeds a upper threshold 1 it will retreat running with the second group parameters. From Figure 5, it can be seen that the SFO is nearly eliminated after 300 number of frames. Since then the fourth group parameters are bring into use all the time unless the lock metric exceed the threshold 0.1 and go back using the third group parameters.
Figure 5. SFO tracking process.
Because the CFO synchronization in the proposed synchronization scheme is the fundamental algorithm, which affects the performance of overall the synchronization scheme, the proposed CFO synchronization algorithm is compared with another two typical methods of frequency offset synchronization in simulation, the algorithm (Schmidl’s algorithm for short) by using training sequence in [17] and the algorithm (Beek’s algorithm for short) based on the cyclic prefix [18]. For the convenience of analysis of results, the mean square error (MSE) between actual offset and the estimation result is adopted. Figures 6 and 7 show the comparison of three CFO synchronization schemes in term of MSE norm. it can be seen that the proposed scheme achieves about 10^{–5} MSE from 0 dB SNR, which is obviously better than the other two methods.
The BER performance of the proposed synchronization scheme in multipath channel is given in Figure 8. In the simulation, we adopt the method of symbol timing synchronization proposed in [19]. The BER is figured up after the PLL is locked. The simulation result shows that, when using proposed synchronization scheme together with robust QPSK modulation and RS + CC with 1/2 code rate, the BER approach 10^{–6} when SNR ≥ 4 dB. It proved that the proposed sampling clock synchronization scheme can achieve preferable performance for the OFDM system with unified driving crystal.
Figure 8. BER performance of proposed synchronization scheme.
6. Conclusion
Regarding the requirement of miniaturization and low cost for wireless personal communication device, on which the same crystal to drive the sampling and channel frequency is equipped generally. The proposed scheme of sampling clock synchronization is perfectly appropriate for this kind of OFDMbased communication system. One of the advantages of the scheme is that it needs not an accurate result of symbol timing synchronization (only with the toleration of ±1/2 cyclic prefix length). This scheme also has the advantage of low complexity that has been verified on a multicore DSP platform for a wireless system in author’s research work.
Competing interests
The authors declare that they have no competing interests.
Acknowledgments
This study was supported by the National Science and Technology Major Projects under grant 2012ZX03003011 and 2012ZX03003007, and the BUPT Research Innovation Project under grant 2011RC0111.
References

SF Mason, CR Berger, SL Zhou, P Willett, Detection, synchronization, and Doppler scale estimation with multicarrier waveforms in underwater acoustic communication. IEEE J. Sel. Areas Commun. 26(9), 1638–1649 (2008)

DK Kim, SH Do, HB Cho, HJ Chol, KB Kim, A new joint algorithm of symbol timing recovery and sampling clock adjustment for OFDM systems. IEEE Trans. Consum. Electron. 44(3), 1142–1149 (1998). Publisher Full Text

BG Yang, KB Letaief, S Roger, Timing recovery for OFDM transmission. IEEE J. Sel. Areas Commun. 18(22), 2278–2291 (2000)

T Pollet, P Spruyt, M Moeneclaey, The BER performance of OFDM systems using nonsynchronized sampling. Proceeding of IEEE Globecom’94, vol. 1, San Francisco, 253–257 (1994)

J Kim, EJ Powers, Y Cho, A nonsynchronized sampling scheme. the 36th Asilomar Conference on Signals, Systems and Computers, vol. 2, Pacific Grove, CA, 1900–1904 (2002)

AI Bo, Y Shen, ZD Zhong, BH Zhang, Enhanced sampling clock offset correction based on time domain estimation scheme. IEEE Trans. Consum. Electron. 57(2), 696–704 (2011)

BG Yang, KB Letaief, RS Cheng, ZG Cao, An improved combined symbol and sampling clock synchronization method for OFDM systems. Proceedings of the 1999 IEEE Wireless Communications and Networking Conference, vol. 3, New Orleans, 1153–1157 (1999)

DW Paranchych, NC Beaulieu, Performance of a digital symbol synchronizer in cochannel interference and noise. IEEE Trans. Commun. 48(11), 1945–1954 (2000). Publisher Full Text

BG Yang, MA ZhengXin, CAO ZhiGang, MLoriented DA sampling clock synchronization for OFDM systems. Proceedings of the International Conference on Communication Technology, vol. 1, Beijing, 781–784 (2000)

FJ Harris, M Ric, Multirate digital filters for symbol timing synchronization in software defined radios. IEEE J. Sel. Areas Commun. 19(12), 2346–2357 (2001). Publisher Full Text

D Wang, AQ Hu, A combined residual frequency and sampling clock offset estimation for OFDM systems. IEEE Asia Pacific Conference on Circuits and Systems, vol.1, Singapore, 1184–1187 (2006)

WD Xiang, T Pratt, A simple cascade carrier frequency and sampling clock offsets estimation method for OFDM systems. First IEEE Consumer Communications and Networking Conference, vol. 1, Las vegas, 718–720 (2004)

H Lee, J Lee, Joint clock and frequency synchronization for OFDMbased cellular systems. IEEE Signal Process. Lett. 18(12), 757–760 (2011)

E Del CastilloSanchez, FJ LopezMartinez, E MartosNaya, JT Entrambasaguas, Joint time, frequency and sampling clock synchronization for OFDMbased systems. Wireless Communications and Networking Conference, vol. 1, Budapest, 1–6 (2009)

M Sliskovic, Carrier and sampling frequency offset estimation and correction in multicarrier systems. IEEE Global Telecommunications Conference, vol. 1, San Antonio, 285–289 (2001)

ITUR M.1225 Guidelines for evaluation of radio transmission technologies for IMT2000[S] (http://www, 1997), . itu.int webcite

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

PH Moose, A technique for orthogonal frequency division multiplexing frequency offset correction. IEEE Trans. Commun. 10(42), 2908–2914 (1994)

JJ Van De Beek, M Sandell, PO Borjesson, ML estimation of time and frequency offset in OFDM systems. IEEE Trans. Signal Process. 45(7), 1800–1805 (1997). Publisher Full Text