FIELD OF INVENTION
The present invention is related to communication systems, and more particularly, to the digital modem system using GFSK transmission scheme.
BACKGROUND OF INVENTION
Please refer to FIG. 1, which shows a typical GFSK receiver structure. The antenna 100 receives a GFSK signal 110 propagating through channel, then the RF circuit 120 removes the carrier from the received GFSK signal 110 and obtains a baseband analog complex signal z(t) 130. The purpose of ADC (analog-to-digital converter) 140 is to sample the analog baseband complex signal 130 to a digital baseband complex signal zk,j 150, and the baseband circuit 160 demodulates and processes the digital baseband complex signal 150 and obtains the original binary sequence b(k) 170.
The invention is related to the digital modem system using GFSK scheme to transmit the signal in the baseband. The digital baseband complex signals are demodulated using differential detection. In other word, the decision of a bit is based on the phase difference between the current and its previous received signals. In addition, to demodulate a signal correctly also requires certain mechanism for synchronization. The synchronization tasks include carrier frequency, phase and symbol timing. Generally, a communication system will provide extra information (so-called preamble or training sequence) to aid the receiver to accomplish these tasks. This present invention develops algorithms to perform differential detection of GFSK signal using decision feedback and preamble (training sequence) for frequency, clock and frame synchronization.
GFSK, which employs Gaussian filter for pulse shaping, is an attractive modulation scheme due to its compact spectrum. However, the Inter-Symbol-Interference (ISI) introduced by the Gaussian filter also degrades the bit error rate (BER) performance. Various receiver structures were proposed to improve the BER performance of GFSK owing to the ISI resulted from Gaussian filter. ABBAS et al (reference b.1) proposed a method using differential detection with decision feedback to overcome the ISI issue. In their original paper, they only dealt with GMSK modulation and assumed that clock and frequency have been synchronized perfectly.
As to the synchronization issue, Mehian et al (reference b.2) proposed a method to estimate the symbol timing and frequency offset without training sequence. In their original paper, they only dealt with GMSK modulation and used conventional differential detection.
This present invention modifies, combines, and extends their works from GMSK to GFSK. For a given pre-known preamble, this invention estimates the frequency offset using the preamble as prior information and estimate the symbol timing using the estimated frequency offset. Once the estimated frequency offset and symbol timing are obtained, this information is used to do differential detection incorporating with decision feedback.
SUMMARY OF INVENTION
A method and a circuit of estimating a binary sequence in a GFSK communication system are disclosed in the present invention. First in response to a complex baseband signal and a preamble data t(k), obtain a frequency offset estimation. Then obtain a complex digital decimated signal by estimating a sampling point. Based on the complex digital decimated signal, the frequency offset estimation, using one-bit and two-bit differential detection technique to demodulate the complex digital decimated signal to generate the binary sequence. Finally, in responsive to the binary sequence and the preamble data t(k), obtain a starting bit of the binary sequence.
BRIEF DESCRIPTION OF DRAWINGS
The present invention will become fully understood from the detailed description given herein below with the accompanying drawings, given by way of illustration only and thus not intend to limit the present invention.
FIG. 1 illustrates a typical GFSK receiver structure.
FIG. 2 illustrates the block diagram of present invention.
FIG. 3 illustrates the flowchart of the operation of the carrier synchronization circuit.
FIG. 4 illustrates the mechanism of finding N1, N2, N3, N4.
FIG. 5 illustrates the flowchart of the operation of the clock synchronization circuit.
FIG. 6 illustrates the flowchart of the operation of the demodulation circuit.
FIG. 7 illustrates the flowchart of the operation of the frame synchronization circuit.
DETAILED DESCRIPTION OF PRESENT INVENTION
To describe the invention clearly, a number of definition of terms used herein are given as follows:
The term “symbol” refers to data represented by more than one bit.
The term “preamble” used herein refers to a data string both sender and receiver agree to use as an information header embedded in the transmission signal.
FIG. 2 illustrates the block diagram of the present invention. One aspect of the present invention is to use the baseband circuit for estimating a binary sequence b(k)={{circumflex over (b)}0, {circumflex over (b)}1, {circumflex over (b)}2, . . . , {circumflex over (b)}k−2, {circumflex over (b)}b−1, {circumflex over (b)}bk, {circumflex over (b)}k+1, . . . } in GFSK communication system. In response to a complex baseband signal zk,j 200 and a preamble data t(k) defined in the GFSK communication system, a carrier synchronization circuit 210 generates a frequency offset estimation {circumflex over (Ω)}Δ 220, wherein zk,j 200 represents a j-th sample of a k-th symbol of a received data, 0≦j≦N−1, N represents samples per symbol. A clock synchronization circuit 230 receives the complex digital signal zk,j 200 and the frequency offset estimation {circumflex over (Ω)}Δ 220, estimates a sampling point ĵ and generates a complex digital decimated signal zk,ĵ 240 by using the sampling point ĵ. A demodulation circuit 250 demodulates the complex digital decimated signal zk,ĵ 240 to generate the binary sequence b(k)={{circumflex over (b)}0, {circumflex over (b)}1, {circumflex over (b)}2, . . . , {circumflex over (b)}k−b, {circumflex over (b)}k−1, {circumflex over (b)}k, {circumflex over (b)}k+1, . . . } 260 in response to the complex digital decimated signal zk,ĵ 240 by using the frequency offset estimation {circumflex over (Ω)}Δ 220, a first bit {circumflex over (b)}k−1 270 and a second bit {circumflex over (b)}k−2 275, wherein the first bit {circumflex over (b)}k−1 270 and the second bit {circumflex over (b)}k−2 275 are generated by passing the binary sequence b(k)={{circumflex over (b)}0, {circumflex over (b)}1, {circumflex over (b)}2, . . . , {circumflex over (b)}k−2, {circumflex over (b)}k−1, {circumflex over (b)}k, {circumflex over (b)}k+1, . . . } through a T delay circuit 261 and two T delay circuit (261+262) respectively.
Finally, in FIG. 2, a frame synchronization circuit 280 receives the binary sequence 260 b(k)={{circumflex over (b)}o, {circumflex over (b)}1, {circumflex over (b)}2, . . . , {circumflex over (b)}k−2, {circumflex over (b)}k−1, {circumflex over (b)}k, {circumflex over (b)}k+1, . . . } and uses the preamble data t(k) defined in the GFSK system, to obtain the starting bit 290 of the binary sequence 260 b(k)={{circumflex over (b)}0, {circumflex over (b)}1, {circumflex over (b)}2, . . . , {circumflex over (b)}k−2, {circumflex over (b)}k−1, {circumflex over (b)}k, {circumflex over (b)}k+1, . . . }.
In one exemplar aspect of the present invention, the carrier synchronization circuit 210 generates the frequency offset estimation {circumflex over (Ω)}Δ 220 by performing the following steps shown in FIG. 3. Fist in step 300, the carrier synchronization circuit 210 stores the first L complex baseband signal zk,j 200, where L is a predetermined number, then step 310 the carrier synchronization circuit 210 performs one-bit differential detection over the L symbols in the complex baseband signal zk,j to obtain a parameter ck,j, step 320 performs summation operation of the parameter ck,j, with respect to the L symbols, and performs summation operation with respect to N sampling points to obtain a complex value V. Step 330 estimates a bias Ωd based on the preamble data t(k) defined in the GFSK system, it is noted that step 330 can be executed before step 310 to get the bias Ωd in advance. Step 340 calculates the frequency offset estimation {circumflex over (Ω)}Δ 220 based on the bias Ωd and the angle of the complex value V.
In one preferred embodiment of the present invention, the carrier synchronization circuit 210 stores the complex baseband signal zk,j 200 in step 300. In step 310, calculates ck,j=zk+1,jzk,j*, wherein the zk,j* is the conjugate complex number of zk,j. In step 320, obtain V by summation of ck,j over L and N, wherein L is a predetermined number, step 330 estimate the bias from the preamble data t(k). {circumflex over (Ω)}d is determined by N1, N2, N3, N4 and h, N1 representing number of occurrences of bk=bk+1=1, N2 representing number of occurrences of bk=bk+1=−1, N3 representing number of occurrences of bk=1, bk+1=−1, N4 representing number of occurrences of bk=−1, bk+1=1, bk is a value of a kth symbol of the preamble data t(k), h is the modulation index of the GFSK communication system.
N1, N2, N3, N4 can be obtained as followed, please refer to FIG. 4, store a first L symbols of the preamble data t(k) in a first memory X of L words 400 (L is a predetermined number), provide a second memory N 410 having a first word, a second word, a third word and a fourth word, each word with an initial value of 0 and having an address comprise of base address and an offset address, then repeating the following steps for L-1 times: (1) obtaining s by retrieving last 2 bits of a data in memory X and do “and” operation with “11” (shown in 420); (2) adding 1 to the content of one of the first word, the second word, the third word and the fourth word which its offset address equaling to s (shown in 420); (3) shifting the data in memory X right by one bit (as in 440). After the repeating steps are done, obtain N1, N2, N3, N4 by retrieving contents of the first word, the second word, the third word and the fourth word.
Finally, step 340 obtain the frequency offset estimation {circumflex over (Ω)}Δ 220 by removing {circumflex over (Ω)}d from arg(V) and normalized by T, wherein arg(V) is an angle of V, T is a symbol time.
Therefore, this invention estimates and removes the bias embedded in the preamble data t(k). And this invention can be applied to all of communication systems as long as the preamble introduces a bias.
In one aspect of the present invention, wherein the clock synchronization circuit 230 generates a complex digital decimated signal zk,ĵ 240 by estimating a sampling point ĵ and performing the following steps. Please refer to FIG. 5. Step 510 corrects ck,j based on {circumflex over (Ω)}ΔT to obtain cfk,j, then in Step 520 obtains a first value uj by summation of |Im(cfk,j)| over L, wherein L is a predetermined number. In step 530, chooses the symbol timing point ĵ with the largest uj by
step 540 outputs the complex digital decimated signal zk,ĵ 240.
In the past, the symbol timing is obtained by calculating cfi,j=ck,j·exp(−j·{circumflex over (Ω)}ΔT) at first, wherein ck,j=zk+1,jzk,j*, zk,j* is the conjugate complex number of zk,j. Then obtains uj by summation of |Im(cfk,j)| over L, L is a predetermined number. This present invention simplifies the operation by obtaining a first value A and a second value B by A=C cos({circumflex over (Ω)}ΔT), B=C sin({circumflex over (Ω)}ΔT), wherein C is a constant, the obtain the value uj by summation of |Im(ck,j)*A+Re(ck,j)*B| over L.
Please refer to FIG. 6, in one aspect of the present invention, the demodulation circuit 250 demodulates the complex digital decimated signal zk,ĵ to generate the binary sequence b(k)={{circumflex over (b)}0, {circumflex over (b)}1, {circumflex over (b)}2, . . . , {circumflex over (b)}k−2, {circumflex over (b)}k−1, {circumflex over (b)}k, {circumflex over (b)}k+1, . . . }, and performs the following steps: Step 600 obtains a first angle θ using a first bit {circumflex over (b)}k−1 and the frequency offset estimation {circumflex over (Ω)}Δ 220. Step 610 obtains a second angle Λ using a second bit {circumflex over (b)}k−2 and the frequency offset estimation {circumflex over (Ω)}Δ 220; Perform one-bit differential detection with respect to the first bit {circumflex over (b)}k−1 to obtain a first value, and perform two-bit differential detection with respect to the second bit {circumflex over (b)}k−2 to obtain a second value in step 620, and obtain a complex sequence using the first value, the first angle θ, the second value, and the second angle Λ. Finally, step 630 obtains a binary sequence b(k)={{circumflex over (b)}0, {circumflex over (b)}1, {circumflex over (b)}2, . . . , {circumflex over (b)}k−2, {circumflex over (b)}b−1, {circumflex over (b)}k, {circumflex over (b)}k+1, . . . } using an imaginary part of the complex sequence.
In one preferred embodiment of the invention. Obtaining an estimated binary sequence b(k)={{circumflex over (b)}0, {circumflex over (b)}1, {circumflex over (b)}2, . . . , {circumflex over (b)}k−2, {circumflex over (b)}k−1, {circumflex over (b)}k, {circumflex over (b)}k+1, . . . } can be implemented by the following equations: obtain a first angle θ by Θ=C2{circumflex over (Ω)}ΔT−hk−1·πhδ; obtain a second angle Λ by
wherein the parameter h is the modulation index defined in GFSK system,
which is derived from the a Gaussian function p(t) defined in the GFSK communication system; then obtain the estimated binary sequence b(k)={{circumflex over (b)}0, {circumflex over (b)}1, {circumflex over (b)}2, . . . , {circumflex over (b)}k−2, {circumflex over (b)}k−1, {circumflex over (b)}k, {circumflex over (b)}k+1, . . . } by wherein C1, C2, C3 are constants.
In another aspect of the present invention, the frame synchronization circuit 280 obtains a starting bit 290 of the binary sequence b(k)={{circumflex over (b)}0, {circumflex over (b)}1, {circumflex over (b)}2, . . . , {circumflex over (b)}k−2, {circumflex over (b)}k−1, {circumflex over (b)}k, {circumflex over (b)}k+1, . . . } 260 by performing the following steps shown in FIG. 7. It should be noted that the binary sequence b(k)={{circumflex over (b)}0, {circumflex over (b)}1, {circumflex over (b)}2, . . . , {circumflex over (b)}k−2, {circumflex over (b)}k−1, {circumflex over (b)}k, {circumflex over (b)}k+1, . . . } 260 used herein is limited to unipolar binary sequence having values of {1, 0} in order to perform the XOR operation afterwards. In response to the binary sequence b(k)={{circumflex over (b)}0, {circumflex over (b)}1, {circumflex over (b)}2, . . . , {circumflex over (b)}k−2, {circumflex over (b)}k−1, {circumflex over (b)}k, {circumflex over (b)}k+1, . . . } 260 and the preamble data t(k) in the complex digital signal 200, step 700 obtains a value by performing XOR-operation to the binary sequence b(k)={{circumflex over (b)}0, {circumflex over (b)}1, {circumflex over (b)}2, . . . , {circumflex over (b)}k−2, {circumflex over (b)}k−1, {circumflex over (b)}k, {circumflex over (b)}k+1, . . . } 260 and the preamble data t(k). Step 710 obtains a series of coefficients by generating a sum of the value, then choose a minimum value of the series of coefficients to find out the starting bit 290 of the binary sequence b(k)={{circumflex over (b)}0, {circumflex over (b)}1, {circumflex over (b)}2, . . . , {circumflex over (b)}k−2, {circumflex over (b)}k−1, {circumflex over (b)}k, {circumflex over (b)}k+1, . . . } 260. FIG. 7 can be further implemented by obtaining a series of coefficients coff(k) by summation of b(k+n)(xor)t(k) over M, wherein M is number of bits, wherein (XOR) is an exclusive-or operation; and choose the starting bit 290 of the binary sequence b(k)={{circumflex over (b)}0, {circumflex over (b)}1, {circumflex over (b)}2, . . . , {circumflex over (b)}k−2, {circumflex over (b)}k−1, {circumflex over (b)}k,{circumflex over (b)}k+1, . . . } 260 by choosing the minimum of the series of coefficient by min(coff(k)).
In this example, because the present invention uses (XOR) operation instead of multiplication, thus the invention uses minimum instead of maximum operation in finding the starting bit 290 of the binary sequence b(k)={{circumflex over (b)}0, {circumflex over (b)}1, {circumflex over (b)}2, . . . , {circumflex over (b)}k−2, {circumflex over (b)}k−1, {circumflex over (b)}k, {circumflex over (b)}k+1, . . . } 260. This method can be applied to all the systems as long as the bit stream information is in the form of unipolar binary (0, 1) sequence.
In the foregoing specification the invention has been described with reference to specific exemplar aspects thereof. It will, however, be evident that various modification and changes may be made to thereto without departing from the broader spirit and scope of the invention. The specification and drawings are, accordingly, to be regarded in an illustrative rather than restrictive sense.