FPGA Implementable Frame Synchronization Algorithm for Burst Mode GMSK

In time division multiple access (TDMA) communication systems, correctly estimating the synchronization parameters is very important for reliable data transfer. The algorithms used for frequency/phase and symbol timing estimates are generally accepted as knowing the start of signal (SoS) parameter. Therefore, within these parameters, the SoS parameter is of particularly great importance. In this study, a reduced version of the SoS estimation algorithm introduced by Hosseini and Perrins is presented to estimate SoS for Gaussian Minimum Shift Keying (GMSK) modulated signals in burst format over additive white Gaussian noise (AWGN) channels. The reduced algorithm can be implemented on FPGA by using half the number of complex multipliers that would be required by the double correlation method and is robust to carrier frequency/phase errors. Simulations performed under 0.1 normalized frequency offset conditions show that the proposed algorithm has a probability of false lock which is less than 2 7 10 × , even at 0 dB SNR level.


Introduction
The GMSK modulation has many charming advantages.It is efficient in terms of power and bandwidth because GMSK is a type of continuous phase modulation (CPM) [1].In addition, GMSK allows the use of non-linear power amplifiers, which makes GMSK suitable for portable transmitters.Since none of the information is carried as amplitude variations, GMSK is more resilient to noise.As a result of these properties, GMSK is used for TDMA communication systems for data or voice transfer, like GSM, deep space communication missions and wireless body area network [2] [3].Despite its attractive properties, using a maximum likelihood (ML) sequence detector is necessary to optimally demodulate GMSK modulation, which is the main disadvantage of a modulation with memory ( [4], p. 246).
The GMSK receiver becomes computationally complex enough with ML sequence detector.To eliminate GMSK receiver complexity, many MSK-type digital linear receivers and reduced complexity Viterbi receivers are proposed in the literature, and very close performances to optimum receiver are attained.Zhao and Giannakis [5] introduce a reduced complexity Viterbi receiver, Kaleh [6] proposes a linear coherent receiver optimum for all 0 b E N values.Wu and Ng   [7] propose coherent and non-coherent linear receiver structures for GMSK.
However, even with ML sequence detectors, GMSK receivers still fall short in terms of performance.ML sequence detector constructed on frequency/phase and symbol synchronization algorithms gains success, i.e. coherent receiver strip.
Symbol synchronization algorithm estimates the symbol offset for finding the correct sampling point to decrease inter symbol interference.Frequency/phase synchronization algorithm is needed to compensate for the effects of the unstable oscillators used in transmitter and receiver.On the other hand, the GMSK receiver becomes more complicated with frequency/phase and symbol timing estimators in order to achieve performance.
To take the burden of synchronization algorithms away from the hardware, symbol level algorithms rather than sample-based ones can be utilized to make use of the resource sharing advantage of FPGA.GMSK synchronization algorithms can roughly be separated into two groups: data aided and non-data aided algorithms [8].Finding closed loop [9] or open loop [10] symbol-based algorithms is also possible in both groups.Aside from their loop structure, the data aided synchronization algorithms use a known bit sequence for performing estimation [11], while the non-data aided algorithms use only the received waveform itself for estimating the synchronization parameters [12].Regardless of the type of synchronization algorithm used, synchronization algorithms generally assume that the SoS is known.Therefore, SoS parameter should first be estimated with enough accuracy to benefit from synchronization algorithms.Moreover, if a data aided algorithm is selected for synchronization, the SoS parameter of the received signal should be exactly estimated in order to yield low variance synchronization parameter estimations.Although digital synchronization algorithms are by now well established in the literature, the frame synchronization algorithms for GMSK modulation stay in the background, behind the frequency/phase and symbol synchronization algorithms.
In general, there are two common strategies for frame synchronization algorithms: energy detection and correlation detection.Energy-detection-based algorithms are much more sensitive to noise level and finding a detection threshold is not a trivial task at low SNR values [13].However, correlation-based algorithms can allow for frame detection at low SNR values by using a moderately short length preamble sequence.In [14], the introduced two-step algorithm, the preamble detection algorithm (step 1) followed by the preamble estimation algorithm (step 2), showed that 32-bit length preambles can be used to obtain an ac-ceptable level of performance.In this paper, we propose a simplified algorithm that calculates the double correlation metric used for frame synchronization with a remarkable reduction in the number of complex multiplication used by employing a very simple quantization on the expected preamble.The proposed algorithm creates the opportunity for the implementation of additional complex synchronization/detection algorithms of the coherent GMSK receiver strip on the FPGA besides saving in number of complex multiplication.In addition to quantization, by using a 32-symbol length Kasami sequence as preamble sequence of the burst mode transmission, the SoS estimation algorithm is performed in a single step by a simple maximum search operation.A question is, if one can to reduce the number of complex multipliers used in FPGA and the complexity of the SoS estimation algorithm to solve the area efficiency problem.
To answer is yes, but proposed frame synchronization approach should be used.
Performed simulations show that the proposed algorithm has a false lock probability of less than nization algorithm given in [14], which is the same as the ad hoc algorithm introduced in ( [15], p. 487), and is referred to as a double correlation in [16].
In the next section, the GMSK signal model is summarized.Section 3 gives the details of the proposed algorithm.Section 4 presents the simulation conditions and the simulation results obtained.Finally, Section 5 concludes the paper.Now, we proceed with the fundamental properties of the GMSK signal model.

GMSK Signal
The complex envelope representation of the GMSK signal, ( ) b s t , can be ob- tained by using the general binary partial response CPM signal with modulation where b E is the energy per bit, T is the symbol duration, and the phase func- tion, ( ) In Equation ( 2), n α are the equal probable and independent information symbols in the binary alphabet { } 1, 1 + − where the phase pulse, ( ) ( ) f t is the frequency pulse supported over the time interval ( ) 0, LT and has the following properties: For the GMSK signal, the frequency pulse, ( ) where ( ) , and B is the −3 dB bandwidth of the Gaussian frequency pulse (i.e. the communication bandwidth).
( ) By substituting Equation (5) into Equation ( 3), the phase pulse, ( ) g t , can be represented for the GMSK as It can be observed from Equation ( 5) that the frequency pulse duration increases as the bandwidth of the pulse decreases.In practical applications, the pulse is truncated to a specified time interval, ( ) f t is sufficiently well approximated.Throughout this paper, this value is taken to be 4 L = , whereas the time bandwidth product is 0.25 BT = .

Frame Synchronization Algorithm
For inference of the reduced algorithm and simulations for the performance measure, the assumed signal transmission scheme is given by Figure 1.It is assumed that the transmitter starts the burst transmission after a silent period of unknown duration.The transmission starts with a known preamble, and continues with random data.The preamble part of the signal is composed of a Kasami sequence of P N symbol length.Since this work scrutinizes the estimation of the SoS parameter, the length or content of the remaining random data is unimportant.
The received complex baseband signal, ( ) where T is the symbol duration.
The estimation of the SoS is finding the boundary of the known preamble from arbitrarily captured P N samples of the received signal given in Equation ( 8).The decision about the inclusion of preamble fully/partially or not by the captured P N samples is made by using a test function.The result of the test function applied to the captured samples is simply compared with a pre-determined threshold value, which can be determined by a selected probability of false alarm.When the threshold value is exceeded, the preamble is assumed to be located within a particular SoS uncertainty window.The true location of the SoS is determined by simply searching the maximum value of the test function within the SoS uncertainty window.
The ad hoc frame synchronization test function introduced in this work is inspired by the preamble (SoS) detection algorithm derived in [14], which sums the double correlation functions introduced in [16] for different correlation lag values.The test function, ( ) where r represents the captured P N samples of the received signal given in Equation ( 8), [ ] s n is the samples of the expected preamble sequence, and is the design parameter for the performance of the test function.By assuming the products of the samples of the known preamble sequence,

*
s n s n d + , are pre-calculated and stored in a read only memory (ROM), the number of complex multipliers required to calculate the test function given in Equation ( 9) is ( ) s n s nT = .
By using the simple quantization given by Equation ( 10 where [ ] ŝ n represents the quantized P N length preamble sequence.
To further reduce the implementation of the SoS estimation algorithm, the Kasami sequence is selected as the preamble because of its correlation properties.
Therefore, estimating the true location of the SoS can be reduced to searching the maximum value of the test function, Equation (11), within the SoS uncertainty window instead of using a computationally complex SoS estimation algorithm.

Simulations and Results
In the simulation, two different preamble structures are used.One of them is the optimum synchronization sequence for M-ary signals [14].The other is a 32-bit Kasami sequence.
The Kasami sequence used as the preamble is selected from the small set of Kasami sequences generated from the 6 1 x x + + primitive polynomial [17], [18].By using the given primitive polynomial, many different 63-bit length Ka-sami sequences can be generated.The Kasami sequence that we used in the synchronization process is selected as the first 32 consecutive bits of one of the possible 63-bit length sequences.The sequences used in the simulations are presented in Table 1.
There are two synchronizers considered in the simulations.Synchronizer #1 is designed to maximize the test function, Equation (11), in an interval detected by the pre-determined threshold value.Synchronizer # 2 is a combined version of SoS detection, using Equation ( 9), and SoS estimation algorithms proposed by Hosseini and Perrins in [14].Basically, simulation of both Synchronizer   The threshold value, λ , is determined using the Neyman-Pearson criterion for a fixed false alarm probability, FA P .The definition of the false alarm proba- bility is given as such that r consists of only noise samples.Throughout the simulation, all of the threshold values are calculated for For ease of implementation, the D′ value in (11) should be selected to be as small as possible in its defined interval.The tradeoff is that system performance decreases with decreasing D′ .The performance of Synchronizer #1 with respect to D′ is presented in Figure 5.The key point is that the relation between D′ and system performance is not linear.Therefore, we can select D′ in such a way that both system performance and complexity are at an acceptable level.
Under these circumstances, 8 D′ = is a good choice.The effect of quantization on system complexity is mentioned in Section 3. At the expense of computation precision loss, a considerable number of complex multipliers are saved.To examine the trade-off between system performance and quantization, the probability of false lock performance of Synchronizer #1 is compared for the two cases, one in which quantization is utilized and one in which it is not.The comparison is presented in Figure 6.According to the results, the effect of quantization on system performance is negligibly small.We would rather to use quantized values instead of using the true values, which are too close to each other for GMSK modulation see

Conclusion
In this paper, to improve the performance of a coherent receiver, which uses SoS estimation, double correlation algorithm was implemented.Especially, we focused on the SoS estimation problem for burst mode GMSK transmission.Firstly, the number of complex multipliers were halved by quantizing the GMSK modulated preamble data, which is used in the double correlation.Secondly, a single step SoS estimation algorithm which simply estimates SoS as the maximum point of the proposed test function in the uncertainty window was developed by using the Kasami sequence as preamble.Complexity of an algorithm has a prior importance in the sense of its efficiency.The combination of reduction in the number of complex multipliers and a single step SoS estimation reduce the complexity of the frame synchronization process drastically.The simulation results showed that the effect of the quantization on system performance is negligibly small.At low SNR values, the proposed Synchronizer #1 and the Synchronizer # 2 given in [14] have almost the same FL P values.For Thus, the power and area efficiencies of the double correlation algorithm increase while the number of complex multiplier decreases with a reasonable performance loss.The proposed approach can be a practical candidate for frame synchronization algorithms of communication systems using GMSK modulation.
SNR level and 0.1 normalized frequency offset condition.The developed algorithm is a reduced version of the frame synchro-

+
), the real and imaginary parts of the products of the samples of the known preamble sequence, , can take integer values over [ ] 2, 2 − interval.Therefore, the complex multiplication operations turn into simple shift and negate operations.Hence, the required number of complex multiplications for calculating Equation (

Figure 2 .Figure 6 .
Figure 6.Effect of quantization of the Kasami sequence by simulating 7 10 frames with 0.1 normalized frequency offset and 8 D′ = .

Figure 7 .
Figure 7.Comparison of synchronizers by simulating 7 10 frames at 0.1 normalized frequency offset condition.
t is the complex baseband AWGN with zero mean and power spectral density 0 N .The proposed algorithm uses the sampled received signal at the symbol rate, i.e. one sample for each symbol.The sampled received signal, [ ] r t , over an AWGN channel can be represented as ( ) ( ) ( ) 2π e , Figure 1.GMSK transmission scheme in burst mode with preamble.w

Table 1 .
Preamble sequences used in the simulations.