A Fast Acquisition Algorithm Overcoming Fuzz Problems for TDDM Spread Spectrum Signal

Abstract

TDDM (time division data modulation) technique will be used in the next generation GNSS (global navigation satellite system) to improve processing performance and to reduce inter-GNSS interference; however, the emergence of TDDM signal causes the estimation frequency and message reversal fuzz problems in the acquisition process of a GNSS receiver. At present, the traditional acquisition methods have some limitations and shortcomings. Therefore, aiming at the unique characteristics of TDDM signal, a fast acquisition algorithm is proposed to overcome these fuzz problems in this paper. In the proposed algorithm, three stages are obtained by some key technologies, which are the I-Q frequency compensation, superposition processing, subsection processing, and reversion position estimation. Besides, the algorithm is simulated from carrier frequency error, code phase error, message inversion error, and processing speed. Theoretical and simulation results show that the new algorithm can quickly overcome the fuzz problems, and the new algorithm is better than the existing algorithm in the speed and accuracy, which demonstrates that this new algorithm is an effective search scheme for the next generation GNSS signals.

1. Introduction

GNSS (global navigation satellite system) [1] with the advantages of high accuracy, high speed, and low cost has become the most widely used application of spread spectrum communication technology [2, 3]. At the same time, TDDM (time division data modulation) [4] technique will be used in the next generation GNSS (global navigation satellite system) to improve processing performance and to reduce inter-GNSS interference. The mechanisms of TDDM modulation are a new spread spectrum technology using long-period code parity. Therefore, the TDDM signal has not only the characteristics of long codes but also the unique characteristics of this modulation method.

In the spread spectrum communication technology research, the signal synchronization research plays a more and more important role in a receiver. Only by owning the signal synchronization ability, we can achieve the purpose of expectant measure, investigation, communication, navigation, and so on. So some synchronization acquisition algorithms are proposed, such as the long code acquisition algorithms [5–7] and the relevant algorithms based on FFT [8, 9]. At the same time, many fast acquisition algorithms [10–13] are also proposed from the perspective of processing speed, and some improved XFAST algorithms [14, 15] have been proposed. Furthermore, some new studies on analysis of the signal synchronization are proposed, such as activation feedback control synchronization scheme [16] and projective synchronization scheme [17]. These algorithms and schemes can reach effective synchronization purpose.

However, considering some synchronization problems, especially acquisition fuzz problems, are caused just because of the TDDM signal comprehensive properties, the above algorithms and schemes cannot solve the acquisition fuzz problems for TDDM signal. At present, the two-channel algorithm [18] is proposed to solve the fuzz problems of TDDM signal. Although the algorithm can deal with fuzz problems, the algorithm precision is not enough, and the processing speed is slow.

In summary, according to the TDDM signal characteristics, combined with the existing research methods, a fast acquisition algorithm is proposed to overcome the fuzz problems for TDDM signal in this paper.

2. The Fuzz Problems of TDDM Signal Acquisition

2.1. The Fuzz Problem Caused by External Factors

Evident frequency error will appear due to relative motion when the signal arrives at the receiver, so the frequency mixing technology of multiple channels is usually used in traditional acquisition algorithms to overcome this problem. However, the residual frequency exists in the mixing signal, which will affect the correlation result. And if the residual frequency satisfies the following equation, the maximum error [19] will appear in the correlation result. Consider 1𝑡×Δ𝑓=2×𝑛(𝑛=1,2⋯+∞), (1) where 𝑡 is the accumulation time, Δ𝑓 is the residual frequency, and 𝑛 is a positive integer.

Because the residual frequency is periodic function, the maximum error will appear in the correlation result when the single cycle of residual frequency satisfies the following equation: 𝑓_𝑒=1,2𝑡 (2) where 𝑓_𝑒 is a cycle of the residual frequency.

Now, the effect of residual frequency is verified by an example. In the condition of 20 ms accumulation time, the theoretical correlation result is shown in Figure 1(a). When the other conditions are the same, but the residual frequency is 25 Hz, the acquisition result was shown in Figure 1(b). These test results show that the correlation result cannot correctly be gotten when the residual frequency satisfies (2). Therefore, the influence of residual frequency is very great which causes fuzz problems of leakage acquisition.

(a) The theoretical correlation result

(b) The acquisition result

Figure 1 

The correlation results in the condition of 20 ms accumulation time.

2.2. The Fuzz Problem Caused by Internal Factors

In TDDM spread spectrum processing, the odd code of PN code is modulated by message, but the even code of PN code is not modulated by message. If the message sign in the received signal is positive and there is no sign inversion, the correct acquisition result can be gotten in the traditional methods. However, if the message sign in the received signal is negative and there is no sign inversion, the test result of traditional acquisition method is shown in Figure 2. If the message sign in the received signal is not only positive but also negative, the test result of traditional acquisition method is shown in Figure 3. These test results showed that the TDDM mechanism caused the acquisition fuzz problems.

(a) The theoretical correlation result

(b) The acquisition result

Figure 2 

The test results in the condition of negative message sign.

(a) The theoretical correlation result

(b) The acquisition result

Figure 3 

The test results in the condition of positive and negative message sign.

3. The Fast Acquisition Algorithm Overcoming Fuzz Problems

Currently, the two-channel algorithm [19] is an effective algorithm for TDDM signal. Firstly, the local code is divided into two channels in the two-channel algorithm. Then, circular correlation operations are completed between each channel signal and received signal of frequency compensation. Finally, each channel maximum peak is calculated, and the message reversion position is calculated. In the process of acquisition, the two-channel algorithm is suitable for the TDDM signal without message inversion, while it is not suitable for the TDDM signal with message inversion. Therefore, this algorithm cannot solve well the fuzz problems caused by TDDM modulation mechanism, and it cannot solve the fuzz problems caused by frequency compensation because the residual frequency influence has not been considered.

Through the above analysis, we can see that the fuzz problems are caused by residual frequency and TDDM mechanism in the process of TDDM signal acquisition. At the same time, the existing algorithm has the limitations of leakage acquisition and erroneous acquisition. Therefore, a fast acquisition algorithm overcoming fuzz problems is proposed, aiming at the unique characteristics of TDDM signal.

3.1. General Scheme

The new algorithm mainly includes search stage, coarse acquisition stage, and precise acquisition stage, and this algorithm general scheme is shown in Figure 4.

Figure 4 

The general scheme of the fast acquisition algorithm.

In the search stage, the search speed is the key problem because of the bigger error, so the 𝐼-𝑄 frequency compensation processing is employed to overcome the frequency fuzz problem; the superposition processing of local code square wave modulation and original code is employed to overcome message reversal fuzz problem. Further, the search speed is enhanced by subsection processing.

In the coarse acquisition stage, the signal phase and frequency error are quickly determined by removing subsection processing and further 𝐼-𝑄 frequency compensation processing.

In the precise acquisition stage, the received signal and local code are divided into two sections, and the length of each section is equal to the duration of half message code. If there is sign reversion in one section of received signal, there is surely no sign reversion in another section. So the technology is used to find reversion section and combined with four channel correlation results to calculate message reversion position.

3.2. Algorithm Principle

The TDDM spread spectrum process is shown in (3), where 𝑑(𝑛) is the message code, 𝐶(𝑛) is the PN code, 𝑆𝑁(𝜔_𝑆𝑛) is the square wave, and 𝜔_𝑆 is the frequency of square wave. If the message sign is positive in accumulation time, the TDDM spread spectrum result is equivalent to the product of message code and PN code. If the message sign is negative, the TDDM spread spectrum result is multiplied by the square wave. Further, the received signal is shown in (4), where 𝜔_𝑅 is the intermediate frequency and 𝜑_𝑅 is the carrier phase. Consider 𝐶_TDDM(𝑛)=(𝑑(𝑛)𝐶(𝑛))|_𝑇=𝑑(𝑛)𝐶(𝑛)(𝑛)𝑑(𝑛)𝐶(𝑛)𝜔,𝑑=1𝑆𝑁_𝑆𝑛,𝑑(𝑛)=−1𝑆_𝑅(𝑛)=𝐴𝐶_TDDM(𝑛)cos𝜔_𝑅𝑛+𝜑_𝑅=𝐴(𝑑(𝑛)𝐶(𝑛))|_𝑇cos𝜔_𝑅𝑛+𝜑_𝑅. (3)(4)

Firstly, the received signal is transmitted into the search stage. In this stage, the parallel 𝐼-𝑄 frequency compensation processes are employed, and the processes are expressed as 𝑆_𝐼(𝑛)=𝑆_𝑅(𝑛)cos𝜔_𝑏𝑛+𝜙_𝑏=𝐴(𝑑(𝑛)𝐶(𝑛))|_𝑇cosΔ𝜔_𝐼𝑆𝑛+Δ𝜑_𝑄(𝑛)=𝑆_𝑅(𝑛)cos𝜔_𝑏+𝜔_𝑒2𝑛+𝜑_𝑏=𝐴(𝑑(𝑛)𝐶(𝑛))|_𝑇cosΔ𝜔_𝑄𝑛+Δ𝜑, (5) where 𝜔_𝑏 is the basal frequency, 𝜑_𝑏 is the basal carrier phase, 𝜔_𝑒 is the residual frequency, Δ𝜑 is the phase after frequency compensation, and Δ𝜔_𝐼 and Δ𝜔_𝑄 are the frequency of 𝐼 channel and 𝑄 channel, respectively.

At the same time, 𝑚 segments local codes are accumulated; each segment length is equal to the length of the received signal, and then the local original code and its square wave modulation results are added, and the process is expressed as 𝐺_𝑖(𝑛)=𝐶_𝑌((𝑖−1)𝑁+𝑛)+𝐶_𝑌((𝑖−1)𝑁+𝑛)𝑆𝑁𝜔_𝑆𝑛,(𝑖=1,2⋯𝑚), (6) where 𝑁 is the point number of correlation processing.

Each segment superposition result is expressed as 𝐺_𝑋(𝑛)=^𝑚_𝑖=1𝐺_𝑖(𝑛)=^𝑚_𝑖=1𝐶_𝑌((𝑖−1)𝑁+𝑛)+𝐶_𝑌((𝑖−1)𝑁+𝑛)𝑆𝑁𝜔_𝑆𝑛. (7)

Then, the circular correlation operation is completed by using the results of (5) and (7), and the process is expressed as 𝑋(𝑛)=^𝑁_𝑚=0𝑆_𝐼(𝑚)+𝑆_𝑄(𝑚)𝐺^∗_𝑋((𝑚−𝑛))_𝑁𝑅_𝑁(𝑛)=^𝑁_𝑚=0𝑆_𝐼(𝑚)𝐺^∗_𝑋((𝑚−𝑛))_𝑁+𝑆_𝑄(𝑚)𝐺^∗_𝑋((𝑚−𝑛))_𝑁×𝑅_𝑁(𝑛)=^𝑁_𝑚=0𝑆_𝐼(𝑚)×^𝑚_𝑖=1𝐶^∗_𝑌(((𝑖−1)𝑁+𝑚−𝑛))_𝑁+𝐶^∗_𝑌(((𝑖−1)𝑁+𝑚−𝑛))_𝑁×𝑆𝑁𝜔_𝑆(𝑚−𝑛)_𝑁+𝑆_𝑄(𝑚)^𝑚_𝑖=1𝐶^∗_𝑌(((𝑖−1)𝑁+𝑚−𝑛))_𝑁+𝐶^∗_𝑌(((𝑖−1)𝑁+𝑚−𝑛))_𝑁×𝑆𝑁𝜔_𝑆(𝑚−𝑛)_𝑁×𝑅_𝑁(𝑛)≈^𝑁_𝑚=0𝐴𝑑(𝑚)cosΔ𝜔_min+𝑚+Δ𝜑_{4𝑚−1𝑁𝑗=1𝑚=0}𝐴𝑑(𝑚)𝐶(𝑚)𝐷_𝑗(𝑚)×cosΔ𝜔_other≈𝑚+Δ𝜑_𝑁𝑚=0𝐴𝑑(𝑚)cosΔ𝜔_min𝑚+Δ𝜑+^4𝑚−1_𝑗=1𝜆_𝑗, (8) where 𝑅_𝑁(𝑛) is the rectangle sequence, 𝐺^∗_𝑋((𝑚−𝑛))_𝑁 is the result of conjugate and circular shifting for 𝐺_𝑋(𝑛), and 𝐷_𝑗(𝑚) is the other PN code.

Because there is no correlation between 𝐷_𝑗(𝑚) and the received signal, the correlation result of 𝐷_𝑗(𝑚) and the received signal is equivalent to noise, and the noise is expressed as 𝜆_𝑖. The maximum proportion peak is calculated using the correlation results, and the process is expressed as 𝑃_𝑉=𝑁⋅max[𝑋(𝑛)]∑^𝑁_𝑛=1[𝑋(𝑛.)] (9)

The frontal operations are redone using new received signal if 𝑃_𝑉 cannot meet the threshold; otherwise, the received signal is transmitted in the coarse acquisition stage.

In the coarse acquisition stage, the estimated frequency value 𝜔_𝑑 is calculated by using the channel of 𝑃_𝑉, and the 𝐼-𝑄 frequency compensation processes for received signal are done anew; they are expressed as 𝑆^′_𝐼(𝑛)=𝑆_𝑅(𝑛)cos𝜔_𝑑𝑛+𝜙_𝑑=𝐴(𝑑(𝑛)𝐶(𝑛))|_𝑇cosΔ𝜔^′_𝐼𝑆𝑛+Δ𝜑_′𝑄(𝑛)=𝑆_𝑅(𝑛)cos𝜔_𝑑+𝜔_𝑒2𝑛+𝜑_𝑑=𝐴(𝑑(𝑛)𝐶(𝑛))|_𝑇cosΔ𝜔^′_𝑄𝑛+Δ𝜑. (10)

The circular correlation between 𝑆^′_𝐼(𝑛), 𝑆^′_𝑄(𝑛), and 𝐺^∗_𝑖((𝑚−𝑛))_𝑁 is expressed as 𝑋_𝐼𝑖(𝑛)=^𝑁_𝑚=0𝑆^′_𝐼(𝑚)𝐺^∗_𝑖((𝑚−𝑛))_𝑁𝑅_𝑁(𝑛)=^𝑁_𝑚=0𝑆^′_𝐼(𝑚)𝐶_𝑌(((𝑖−1)𝑁+𝑚−𝑛))_𝑁+𝐶_𝑌(((𝑖−1)𝑁+𝑚−𝑛))_𝑁×𝑆𝑁𝜔_𝑆(𝑚−𝑛)_𝑁𝑅_𝑁(𝑛)≈^𝑁_𝑚=0𝐴𝑑(𝑚)cosΔ𝜔^′_𝐼+𝑚+Δ𝜑_{4𝑚−1𝑁𝑗=1𝑚=0}𝐴𝑑(𝑚)𝐶(𝑚)𝐷_𝑗(𝑚)×cosΔ𝜔^′_𝐼≈𝑚+Δ𝜑_𝑁𝑚=0𝐴𝑑(𝑚)cosΔ𝜔^′_𝐼𝑚+Δ𝜑+^4𝑚−1_𝑗=1𝜆_𝑗𝑋_𝑄𝑖(𝑛)=^𝑁_𝑚=0𝑆^′_𝑄(𝑚)𝐺^∗_𝑖((𝑚−𝑛))_𝑁𝑅_𝑁(𝑛)=^𝑁_𝑚=0𝑆^′_𝑄(𝑚)𝐶_𝑌(((𝑖−1)𝑁+𝑚−𝑛))_𝑁+𝐶_𝑌(((𝑖−1)𝑁+𝑚−𝑛))_𝑁×𝑆𝑁𝜔_𝑆(𝑚−𝑛)_𝑁𝑅_𝑁(𝑛)≈^𝑁_𝑚=0𝐴𝑑(𝑚)cosΔ𝜔^′_𝑄+𝑚+Δ𝜑_{4𝑚−1𝑁𝑗=1𝑚=0}𝐴𝑑(𝑚)𝐶(𝑚)𝐷_𝑗(𝑚)×cosΔ𝜔^′_𝑄≈𝑚+Δ𝜑_𝑁𝑚=0𝐴𝑑(𝑚)cosΔ𝜔^′_𝑄𝑚+Δ𝜑+^4𝑚−1_𝑗=1𝜆_𝑗. (11)

Further, the maximum proportion peaks are calculated using the above results, as (12). The final estimated frequency is calculated using the channel number of 𝑃_𝑋𝑖, and the code phase is calculated using the position of 𝑃_𝑋𝑖. Through the above operations, the carrier frequency and code phase can be controlled in a certain range. However, the message inversion fuzz problem still exists, so the received signal is transmitted into the precise acquisition stage to estimate message inversion position. Consider 𝑃_𝐼𝑖=𝑁⋅max𝑋_𝐼𝑖(𝑛)∑^𝑁_𝑛=1𝑋_𝐼𝑖(𝑛𝑃_𝑄𝑖)=𝑁⋅max𝑋_𝑄𝑖(𝑛)∑^𝑁_𝑛=1𝑋_𝑄𝑖(𝑛𝑃_𝑋𝑖)=maxmax𝑃_𝐼𝑖,max𝑃_𝑄𝑖. (12)

In the acquisition stage, the received signal is divided into two sections, in which each section’s length is equal to the duration of half message code. The frequency compensation of each section is processed by using the result of coarse acquisition stage, and the processes are expressed as 𝑆₁(𝑛)=𝐴(𝑑(𝑛)𝐶(𝑛))|_𝑇cos𝜔_𝑅𝑛+𝜑_𝑅cos𝜔_𝑔𝑛+𝜑_𝑔=𝐴(𝑑(𝑛)𝐶(𝑛))|_𝑇cosΔ𝜔_𝑔𝑛+Δ𝜑_𝑔,(),𝑆𝑛=1,2⋯𝐿₂(𝑛)=𝐴(𝑑(𝑛)𝐶(𝑛))|_𝑇cos𝜔_𝑅𝑛+𝜑_𝑅cos𝜔_𝑔𝑛+𝜑_𝑔=𝐴(𝑑(𝑛)𝐶(𝑛))|_𝑇cosΔ𝜔_𝑔𝑛+Δ𝜑_𝑔,(𝑛=𝐿+1,⋯2𝐿), (13) where Δ𝜔_𝑔 is the frequency error and 𝐿 is the sampling point.

Then, the local code is generated anew according to the result of coarse acquisition stage, and it is also divided into two sections. The first half of sequence is expressed as in (14) and the second half of sequence is expressed as in (15). Consider 𝐺₁(𝑛)=𝐶_𝑌(𝑛),()𝐺𝑛=1,2⋯𝐿_′1(𝑛)=𝐶_𝑌(𝑛)𝑆𝑁𝜔_𝑆𝑛,(𝑛=1,2⋯𝐿)𝐺₂(𝑛)=𝐶_𝑌(𝑛),()𝐺𝑛=𝐿+1,𝐿+2⋯2𝐿_′2(𝑛)=𝐶_𝑌(𝑛)𝑆𝑁𝜔_𝑆𝑛,(𝑛=𝐿+1,𝐿+2⋯2𝐿). (14)(15)

Furthermore, the correlation process is expressed as 𝑋₁=^𝐿_𝑛−1𝐺₁(𝑛)𝑆₁(𝑛)𝑋^′₁=^𝐿_𝑛−1𝐺^′₁(𝑛)𝑆₁(𝑛)𝑋₂=^2𝐿_𝑛=𝐿+1𝐺₂(𝑛)𝑆₂(𝑛)𝑋^′₂=^2𝐿_𝑛=𝐿+1𝐺^′₂(𝑛)𝑆₂(𝑛). (16)

The maximum of four correlation results is expressed as in (17), and the minimum is expressed as in (18). Consider 𝑋_max=max𝑋₁,𝑋^′₁,𝑋₂,𝑋^′₂𝑋_min=max𝑋₁,𝑋^′₁,𝑋₂,𝑋^′₂. (17)(18)

Let 𝑋_Δ be an intermediate variable; the process methods in the channel of 𝑋_Δ and 𝑋_max are the same; 𝑋_Δ and 𝑋_max are in different section; let 𝜉 be the estimated inversion position of each section. If the sign of received signal is positive and local signal is PN code itself, then 𝑍 is used to express the correlation coefficient. If the sign of received signal is negative and local signal is the modulation result of PN code and square wave, then 𝐹 is used to express the correlation coefficient. 𝑍^′ is used to express the correlation coefficient in other conditions.

If inversion position is in the first section and the sign is from positive to negative, (19) is established, and then it is derived as (20). Furthermore, (21) is established, due to 𝑋_max=𝑋^′₂, 𝑋_min=𝑋₂, and  𝑋_Δ=𝑋^′₁. Consider 𝑋₁()𝑋=𝜉⋅𝑍+𝐿−𝜉⋅𝐹_′1=𝜉⋅𝐹+(𝐿−𝜉)⋅𝑍^′𝑋₂𝑋=𝐿⋅𝐹_′2=𝐿⋅𝑍^′,𝑋𝜉=^′2−𝑋^′₁⋅𝐿𝑋^′₂−𝑋₂,𝑋𝜉=^max−𝑋_Δ⋅𝐿𝑋_max−𝑋_min. (19)(20)(21)

In the same ways, if inversion position is in the first section and the sign is from negative to positive, (22) and (23) are established. If inversion position is in the second section and the sign is from positive to negative, (24) and (25) are established. If inversion position is in the second section and the sign is from negative to positive, (26) and (27) are established. Consider 𝑋₁()𝑋=𝜉⋅𝐹+𝐿−𝜉⋅𝑍_′1=𝜉⋅𝑍^′+()𝑋𝐿−𝜉⋅𝐹₂𝑋=𝐿⋅𝑍_′2𝑋=𝐿⋅𝐹,𝜉=²−𝑋₁⋅𝐿𝑋₂−𝑋^′₂𝑋𝜉=^max−𝑋_Δ⋅𝐿𝑋_max−𝑋_min𝑋₁𝑋=𝐿⋅𝑍_′1𝑋=𝐿⋅𝐹₂()𝑋=𝜉⋅𝑍+𝐿−𝜉⋅𝐹_′2=𝜉⋅𝐹+(𝐿−𝜉)⋅𝑍^′,𝑋𝜉=²−𝑋^′₁⋅𝐿𝑋₁−𝑋^′₁,𝑋𝜉=^Δ−𝑋_min⋅𝐿𝑋_max−𝑋_min𝑋₁𝑋=𝐿⋅𝐹_′1=𝐿⋅𝑍^′𝑋₂()𝑋=𝜉⋅𝐹+𝐿−𝜉⋅𝑍_′2=𝜉⋅𝑍^′+()𝑋𝐿−𝜉⋅𝐹𝜉=^′2−𝑋₁⋅𝐿𝑋^′₁−𝑋₁,𝑋𝜉=^Δ−𝑋_min⋅𝐿𝑋_max−𝑋_min. (22)(23)(24)(25)(26)(27)

Finally, it can be derived from (21), (23), (25), and (27) that the message inversion position is 𝑋𝜉=^max−𝑋_Δ⋅𝐿𝑋_max−𝑋_min,𝑋_max⊂𝑋₂,𝑋^′₂𝑋𝐿+^Δ−𝑋_min⋅𝐿𝑋_max−𝑋_min,𝑋_max⊂𝑋₁,𝑋^′₁. (28)

4. Simulation and Analysis

The new algorithm is simulated and analyzed based on MATLAB simulation platform. The basic parameters are as follows: message rate of 100 Hz, PN code rate of 1 MHz, carrier frequency of 3 MHz, and sampling frequency of 12 MHz.

4.1. Algorithm Simulation

The simulation parameters in the first example are as follows: carrier frequency offset is 845 Hz, PN code offset is 42 ms, accumulation time is 10 ms, SNR is −12 dB, and the message sign is always positive. The simulation parameters in the second example are as follows: carrier frequency offset is 1190 Hz, PN code offset is 24 ms, accumulation time is 10 ms, signal to noise ratio is −12 dB, message sign is from negative to positive, and reversion position is 26.5 ms.

The simulation results of search stage are shown in Figures 5 and 6, and the proportion peak results are shown in Figures 7 and 8. Figures 5 and 7 show that the proportion peak of the ninth channel is maximum value, thus preliminarily estimating carrier frequency as 800 Hz in the first example. Figures 6 and 8 show that the proportion peak of the 13th channel is maximum value, thus preliminarily estimating carrier frequency as 1200 Hz in the second example.

Figure 5 

The simulation result of search stage in the first example.

Figure 6 

The simulation result of search stage in the second example.

Figure 7 

The proportion peak result of search stage in the first example.

Figure 8 

The proportion peak result of search stage in the second example.

In the coarse acquisition stage, the maximum proportion peak results of 𝐼 channel and 𝑄 channel are shown in Figures 9 and 10. Figure 9 shows that the maximum proportion peak of 𝑄 channel is greater than 𝐼 channel; thus further estimated carrier frequency is 825 Hz in the first example. Figure 10 shows that the maximum proportion peak of 𝐼 channel is greater than 𝑄 channel, thus estimating carrier frequency as 1200 Hz in the second example.

Figure 9 

The maximum proportion peak of coarse acquisition stage in the first example.

Figure 10 

The maximum proportion peak of coarse acquisition stage in the second example.

At the same time, the PN code offset is estimated as 42 ms by the position of maximum proportion peak in the first example, and it is shown in Figure 11. The PN code offset is estimated as 24 ms in the second example, and it is shown in Figure 12. Four correlation peaks in the precise acquisition stage are used in (28); thus the message sign is estimated as positive in the first example. The message sign is estimated from negative to positive, and the estimated reversion position is 26.58 ms in the second example.

Figure 11 

The final estimated PN code offset result in the first example.

Figure 12 

The final estimated PN code offset result in the second example.

All simulation results of the first example show that the frequency error is 20 Hz, the PN code phase error is 0, and the message reversion error is 0. The second example results show that the frequency error is 10 Hz, the PN code phase error is 0, and the message reversion error rate is 0.8%.

4.2. Performance Analysis

This chapter aims to analyze this new algorithm adaptability and validity from diverse angles.

The effects of Doppler and message inversion on the search stage are not obvious because this stage is the initial acquisition stage. Therefore, the effects of the segment number selection and SNR are analyzed in this stage. The relationship between the segment number selection and the proportion peak is shown in Figure 13, which shows that the proportion peak decreases gradually along with the increase in segment number. The relationship between SNR and the proportion peak is shown in Figure 14, which shows that the proportion peak decreases gradually along with the decrease of SNR. Because the proportion peak will not meet the threshold requirement when SNR reaches −25 dB, the adaptation SNR of this algorithm is not less than −25 dB.

Figure 13 

The relationship between the segment number and the proportion peak in search stage.

Figure 14 

The relationship between the SNR and the proportion peak in search stage.

The effect of SNR in the coarse acquisition stage is obvious, and the relationship between SNR and the proportion peak is shown in Figure 15, which shows that the proportion peak decreases gradually along with the decrease of SNR. The proportion peak of 𝑄 channel is always greater than 𝐼 channel, so the compensation frequency of 𝑄 channel is more close to the real value in this condition.

Figure 15 

The relationship between the SNR and the proportion peak in coarse acquisition stage.

The effects of Doppler, message inversion, and SNR on the precise acquisition stage are obvious because this stage is a precise stage. The estimation inversion rate result in precise acquisition stage is shown in Figure 16 under the conditions that Doppler error is zero and SNR is zero, and the estimation inversion rate error is shown in Figure 17. The results show that the estimation error is less than 1.5% and it can meet the acquisition requirement. Finally, the relationship between estimation inversion rate error and SNR is shown in Figure 18, which proves that the estimation error increases gradually along with the decrease of SNR, and the effect is irregular. From the previous results, it can be seen that this algorithm estimation inversion rate error is less than 2.5% in the condition of more than −25 dB.

Figure 16 

The estimation inversion rate result in precise acquisition stage.

Figure 17 

The estimation inversion rate error in precise acquisition stage.

Figure 18 

The relationship between estimation inversion rate error and SNR in precise acquisition stage.

4.3. Algorithm Comparison

In order to verify the superiority of the new algorithm, this algorithm is compared with two-channel algorithm on the premise of successful acquisition. The simulation parameters are as follows: carrier frequency offset is 845 Hz, PN code offset is 42 ms, accumulation time is 10 ms, SNR is −12 dB, the message sign is from negative to positive, and the reversion position is 44.5 ms.

With the change of search range, the acquisition time of the two algorithms is shown in Table 1, which shows that the processing speed of the new algorithm increases 9.5 times more than two-channel algorithm. With the change of Doppler, the estimation frequency error of the two algorithms is shown in Table 2, which shows that the new algorithm can obviously reduce the estimation frequency error 6 times less than the two-channel algorithm. Furthermore, with the change of code phase, the estimation code phase error of the two algorithms is shown in Table 3, which shows that the estimation code phase error of the two algorithms is not obvious. With the change of message inversion rate, the estimation inversion error rate of the two algorithms is shown in Table 4, which shows that the new algorithm can obviously reduce the estimation inversion error rate 38 times less than the two-channel algorithm.

Table 1 

The acquisition time of two algorithms with different search range.

Table 2 

The estimation frequency error of two algorithms with different Doppler.

Table 3 

The estimation code phase error of two algorithms with different code phase.

Table 4 

The estimation inversion error rate of two algorithms with different message inversion rate.

5. Conclusion

In this paper, the principle and characteristics of TDDM spread spectrum signal are studied. At the same time, the fuzz problem caused by external factors and internal factors is analyzed, and the two-channel algorithm for TDDM signal is mainly studied. Aiming at the unique characteristics of TDDM signal, a fast acquisition algorithm is proposed to overcome the fuzz problems. The new algorithm is composed of three parts, which are search stage, coarse acquisition stage, and precise acquisition stage. In the search stage, the 𝐼-𝑄 frequency compensation processing is employed to overcome the frequency fuzz problem; the superposition processing of local code and square wave modulation is employed to overcome message reversal fuzz problem. Further, the search speed is enhanced by subsection processing. In the coarse acquisition stage, the signal phase and frequency are quickly determined by removing subsection processing and 𝐼-𝑄 frequency compensation processing. In the precise acquisition stage, the received signal and local code are divided into two sections, and four channel correlation results are combined to calculate message reversion position. Finally, the new algorithm is simulated. The simulation results show that the new algorithm can quickly estimate carrier frequency, code phase, and message inversion position. It can overcome the above fuzz problems, the adaptation SNR of this algorithm is not less than −25 dB, and this algorithm estimation inversion rate error is less than 2.5%. In the same condition, the processing speed of the new algorithm is better than the traditional algorithm. And the estimation frequency error and inversion error of the new algorithm are obviously less than the traditional algorithm.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was supported by the Program for Liaoning Innovative Research Team in University (no. LT2011005), New Century Program for Excellent Talents of Ministry of Education of China (no. NCET-11-1013), Project of Science and Technology Department of Liaoning Province (no. 20121038), Project of Education Department of Liaoning Province (no. L2013085), and the Open Foundation of Key Laboratory of Shenyang Ligong University.