劉曉健，李名祺

(上海無線電設(shè)備研究所，上海200090)

1 引 言

基于機(jī)動平臺的時差定位系統(tǒng)是近年來在宇航測控領(lǐng)域中一個受到重視的研究課題［1－2］，它可以飛機(jī)或衛(wèi)星編隊(duì)實(shí)現(xiàn)對目標(biāo)的精確定位和跟蹤，定位的精度很大程度上取決于站間時間同步系統(tǒng)的同步精度，例如，10 ns的時間誤差就會產(chǎn)生約3 m的定位誤差。

通過外部頻率源進(jìn)行時間同步是當(dāng)今應(yīng)用最廣泛的技術(shù)手段。通過衛(wèi)星導(dǎo)航系統(tǒng)，如GPS 和“北斗”進(jìn)行時間同步，因?yàn)槭艿酱髿怆婋x層散射等因素的干擾，其高精度只是統(tǒng)計(jì)意義上的，實(shí)際很難超過20 ns，且野值較多［3－5］。轉(zhuǎn)發(fā)式衛(wèi)星授時方式的精度稍高，然而也只能達(dá)到十幾納秒的水平［6］。其他外部授時方式，如長波、短波、低頻時碼，授時精度約0.1 μs～0.1 s［7］，遠(yuǎn)低于定位精度的要求。

通過自組織的授時網(wǎng)絡(luò)進(jìn)行時間同步是獲得更高精度的有效途徑。由于測量站點(diǎn)間的距離通常較近，空間傳輸造成的誤差遠(yuǎn)小于衛(wèi)星系統(tǒng)，可以達(dá)到優(yōu)于1 ns的時間同步精度。為了實(shí)現(xiàn)這個目標(biāo)，在設(shè)計(jì)和實(shí)現(xiàn)過程中以下幾項(xiàng)技術(shù)得到了綜合運(yùn)用:一是雙向時間同步技術(shù)，它在雙向鐘差測量過程中消除了信道中的公共誤差［8－9］;二是非等量采樣技術(shù)，偽碼跟蹤環(huán)路采用該技術(shù)后，可以通過較低的采樣率得到優(yōu)于碼片寬度百分之一的時間分辨率［10－13］;三是主站持續(xù)廣播偽碼序列，使得接收站能夠始終跟蹤主站的時鐘頻率，提高計(jì)時精度并降低環(huán)路的震蕩;四是優(yōu)化電路設(shè)計(jì)，降低板噪。

理想條件下，以上幾項(xiàng)措施可以在系統(tǒng)層面將時間同步誤差降低到1 ns左右。但實(shí)際情況中，站點(diǎn)之間存在諸如相對運(yùn)動、頻率源差異、環(huán)路抖動等不利影響，如果對這些綜合作用的因素不加應(yīng)對，時間同步的誤差很容易退化到10 ns量級。此外，用戶出于使用需求，常常要求時間同步模塊以遠(yuǎn)高于1 Hz的頻率輸出校時脈沖，這樣會大大壓縮每次獲取站間時差的處理時間，如果授時機(jī)制設(shè)計(jì)不當(dāng)，無法在短時間完成，或者時統(tǒng)精度難以滿足要求。針對上述問題，本文首先對誤差的來源和影響進(jìn)行分析，接著對于其中引起誤差的主要因素給出了具體的控制措施，并對這些措施的效果進(jìn)行了估計(jì)。

2 雙向時間同步原理

雙向時間同步是采用往返計(jì)時消息同步方式來實(shí)現(xiàn)，其基本原理如圖1所示。假設(shè)從站相對主站的鐘差為ε，且上下行傳輸時延相同。從站在1 時刻向主站發(fā)送詢問消息，經(jīng)過傳輸時延tp后在時刻2 到達(dá)主站;在時刻3(與0 時刻的時間差為固定值td)，主站將測得的詢問包的本地到達(dá)時刻TI通過應(yīng)答包向從站發(fā)送，經(jīng)過傳輸時延tp后，在時刻4到達(dá)從站;從站測量應(yīng)答包在本地的到達(dá)時刻TR。由于tp=TI－ε=TR+ε－td，可以得出

從站根據(jù)上式可計(jì)算出相對主站的鐘差ε，并調(diào)整本地的計(jì)時時鐘，以實(shí)現(xiàn)精同步。由式(1)可見，ε 的精度取決于對TI和TR的準(zhǔn)確測量。本系統(tǒng)通過測量偽碼相位的方法實(shí)現(xiàn)該目的。

圖1 往返計(jì)時消息同步方式原理示意圖Fig.1 The schematic diagram of two－way timing scheme

3 誤差來源

鐘差ε 的誤差主要來自TI和TR的測量誤差以及上下行傳輸時延的不對稱兩方面。

3.1 TI 和TR 的測量誤差

3.1.1 主從站的時鐘偏差

偽碼相位的精確測量是高精度授時的基礎(chǔ)，在從站清零時刻、詢問包到達(dá)時刻和應(yīng)答包到達(dá)時刻處，都會通過偽碼相位測出精確的到達(dá)時刻，主站向從站一直發(fā)送偽碼序列，由于主從站的時間偏差會導(dǎo)致主從站產(chǎn)生的偽碼周期發(fā)生變化，從而導(dǎo)致計(jì)時誤差。

主站從站參考頻率都為100 MHz，現(xiàn)假設(shè)主站真實(shí)頻率為f0，從站真實(shí)頻率為f0+Δf，偽碼速率均為10.23 MHz，偽碼跟蹤環(huán)頻率控制字更新時間為一個偽碼周期。偽碼數(shù)控振蕩器(NCO)的頻率控制字為46 b，從站以恢復(fù)出來的主站偽碼為計(jì)時基準(zhǔn)。

主從站初始偽碼頻率控制字

主站產(chǎn)生的偽碼頻率為

從站的初始偽碼頻率為

由于主從站間的頻差導(dǎo)致的在一個偽碼周期內(nèi)的誤差為

假設(shè)主站頻率為100 MHz，主從站間的頻差為10 Hz，由于時鐘不一致在一個偽碼周期內(nèi)引入的測量誤差為0.01 ns。在從站清零時刻，TI和TR處測量偽碼相位時都會引入此誤差，則引入的最大誤差為0.03 ns。

3.1.2 多普勒效應(yīng)

由于多普勒效應(yīng)會導(dǎo)致偽碼周期的變化，它們同樣導(dǎo)致在1、2、3 點(diǎn)處的測量誤差，效果等同于主從站間的頻率偏差導(dǎo)致的誤差。本系統(tǒng)要求的最大相對速度為30 m/s，由此偽碼產(chǎn)生的多普勒頻移為

由于相對運(yùn)動導(dǎo)致的偽碼多普勒頻移很小，一個偽碼周期時間范圍內(nèi)的時間測量誤差為

3.1.3 偽碼環(huán)路跟蹤誤差

(1)采樣率對精度的影響

偽碼跟蹤環(huán)是實(shí)現(xiàn)偽碼相位精確測量的主要工具，而影響其相位分辨率的關(guān)鍵因素之一是采樣速率。Quirk 等［10］人于2002年指出，當(dāng)采樣速率與碼片速率接近互質(zhì)時，可以通過很低的采樣率實(shí)現(xiàn)碼片寬度1%的相位分辨率，這種技術(shù)被稱為非等量采樣技術(shù)(None Commensurate Sampling，NCS)。仿真和試驗(yàn)證明，對于10.23 MHz 的碼片速率，100 MHz的采樣速率可以讓偽碼跟蹤環(huán)路獲得0.07 ns以上的時間測量精度。

(2)碼跟蹤環(huán)熱噪聲誤差

在工程實(shí)際中，延遲鎖定環(huán)不可避免地受到噪聲的干擾。噪聲的來源主要有兩類，一類是隨輸入信號一同進(jìn)入環(huán)路的輸入噪聲和諧波干擾，另一類為跟蹤環(huán)內(nèi)部器件產(chǎn)生的白噪聲。噪聲和干擾的作用勢必會降低跟蹤環(huán)路的性能，使環(huán)路的輸出相位產(chǎn)生抖動。下面在忽略碼環(huán)白噪聲的情況下，僅考慮輸入高斯白噪聲的影響。噪聲是環(huán)路在跟蹤過程達(dá)到穩(wěn)態(tài)后，本地偽碼對接收偽碼存在相位誤差。

設(shè)輸入噪聲是帶限零均值高斯白噪聲，其雙邊功率譜密度是N0/2。當(dāng)無數(shù)據(jù)調(diào)制時，噪聲對回路的影響最嚴(yán)重。此時，噪聲在環(huán)路中引入的均方根跟蹤抖動為［15］

式中，ρL為回路帶寬內(nèi)的信噪比，ρIF為中頻帶寬內(nèi)的信噪比。

其中，Bn為中頻濾波器的單邊等效噪聲帶寬，本地取為10 kHz。當(dāng)環(huán)路帶寬為5 Hz、載噪比C/N0=37 dB－Hz、偽碼速率為10.23 Mchip/s時，熱噪聲誤差為Tdσε=0.49 ns，其中Td為超前偽碼和標(biāo)準(zhǔn)偽碼的間隔，本系統(tǒng)取1/10 碼片。

3.2 上下行傳輸時延的不對稱

在雙向同步的原理介紹中，上下行傳輸試驗(yàn)被假設(shè)為相等。實(shí)際情況中，由于發(fā)射、接收設(shè)備電路的不一致性，和站間相對運(yùn)動的存在，上下行傳輸時延會存在較為明顯的差別，通常范圍在幾十到上百微秒左右。發(fā)射接收設(shè)備的不一致造成的時延誤差近似為恒定值，可以通過多次測量求得后予以消除。而站間相對運(yùn)動造成的傳輸時延不對稱則是一個動態(tài)的量值，假設(shè)上下行鏈路的時延差為tε，則ε =的大小決定于相對運(yùn)動的速度以及詢問和應(yīng)答消息的發(fā)送時間間隔。

根據(jù)裝載平臺的不同，相對運(yùn)動速度為3～30 m/s;詢問和應(yīng)答消息的發(fā)送時間間隔的下限為詢問消息的持續(xù)時間(本系統(tǒng)為16 ms左右)，實(shí)際情況中還要考慮在兩種消息之間加入保護(hù)間隔和應(yīng)答消息報(bào)頭的持續(xù)時間。保護(hù)間隔的取值要考慮站間距離，在站間相隔100 km的情況下，2 ms的保護(hù)間隔通常就已足夠。報(bào)頭的長度與波形設(shè)計(jì)有很大關(guān)系，如果發(fā)射波形能夠讓接收機(jī)的偽碼跟蹤環(huán)持續(xù)并保持穩(wěn)定，則報(bào)頭的長度只需要幾十比特，否則其持續(xù)時間需要數(shù)百毫秒鎖定以保證偽碼跟蹤環(huán)的鎖定。

在本系統(tǒng)中，詢問和應(yīng)答消息的發(fā)送時間間隔定為18 ms，則相對運(yùn)動在上下鏈路中引入的不對稱誤差為

綜上所述，各種因素造成誤差總量約為1.4 ns。

4 誤差控制方法

根據(jù)前文的分析，各種因素對精度影響的程度從大到小依次為上下行傳輸時延的不對稱、偽碼環(huán)路跟蹤誤差和主從站的時鐘偏差，其中，前兩者是最主要的誤差來源。下面重點(diǎn)討論如何對這兩個因素導(dǎo)致的誤差進(jìn)行控制，對其余的誤差控制也將簡要說明。

4.1 補(bǔ)償上下行鏈路的不對稱

補(bǔ)償上下行鏈路不對稱的一種方法是通過測量信號的多普勒頻偏，從而估計(jì)出相對運(yùn)動速度，進(jìn)而計(jì)算出上下行鏈路傳輸時延的補(bǔ)償量。具體做法是通過頻率跟蹤環(huán)路獲取多普勒頻偏，估計(jì)出當(dāng)前的相對運(yùn)動速度，進(jìn)而估算出詢問和應(yīng)答消息的發(fā)送時間間隔內(nèi)主從站之間的距離變化量，從而對上下行時延的不對稱進(jìn)行補(bǔ)償。由于詢問和應(yīng)答消息的發(fā)送時間間隔僅為18 ms，可以近似認(rèn)為在這段時間內(nèi)主從站之間是相對勻速直線運(yùn)動。采用這一方法的一個重要前提是，主從站的頻率源必須具有很高的頻率準(zhǔn)確度與頻率穩(wěn)定度。在目前國內(nèi)的技術(shù)條件下，可以采用經(jīng)過減振設(shè)計(jì)［14］的恒溫晶振。

頻率估計(jì)是通過載波頻率捕獲和跟蹤完成的，其過程分為載波頻率粗捕獲和精跟蹤兩部分。載波頻率的粗捕獲可以由并行頻率捕獲完成或鎖頻環(huán)完成［16］，有效頻率捕獲范圍在Ku 頻段可達(dá)到±100 kHz左右，但是這種方法精度有限，在±1 kHz左右。更高的測量精度需要通過二階鎖相環(huán)路獲得，根據(jù)實(shí)測，其跟蹤精度可以達(dá)到±100 Hz左右，在13 GHz的載頻上相當(dāng)于2.3 m/s的相對運(yùn)動速度。因此，可以將上下行鏈路不對稱的誤差減少到

另一種解決方法是雙程雙向測量體制，即按照圖1進(jìn)行完一次雙向時間同步以后，將主站和從站的角色互換，重新執(zhí)行一次雙向時間同步。如果以上兩次時間同步的間隔很短，可以近似認(rèn)為主從站之間為勻速直線運(yùn)動，這樣兩次求得的鐘差分別為將兩者相加除以2 即可消除上下行鏈路的時延差tε。但是這種方法會增加時統(tǒng)網(wǎng)絡(luò)時序設(shè)計(jì)的復(fù)雜度，在站點(diǎn)數(shù)較多的情況下，效率較低，因此本文采用時頻偏估計(jì)的方法。

4.2 改善偽碼環(huán)路的跟蹤誤差

由式(1)～(2)可見，對偽碼環(huán)路跟蹤誤差影響最大的兩個因素分別是載噪比和環(huán)路帶寬。但是減少環(huán)路帶寬會影響環(huán)路的相應(yīng)速度，嚴(yán)重時會導(dǎo)致發(fā)散，因而減少跟蹤誤差主要考慮如何增加載噪比。其手段主要是提高鏈路增益，在載波頻率確定的前提下，可以通過提高發(fā)射功率、增加天線增益和定時同步算法來實(shí)現(xiàn)。實(shí)際工作中，如能保證載噪比大于45 dB－Hz，就可以使得該項(xiàng)誤差保持在0.1 ns左右。

此外，偽碼通道如果加載數(shù)據(jù)，相鄰比特的極性翻轉(zhuǎn)會導(dǎo)致環(huán)路的振蕩加劇，使得誤差明顯超出理論的估計(jì)值。解決的方法是在將偽隨機(jī)碼和詢問/應(yīng)答消息調(diào)制到同一個頻點(diǎn)上正交的兩路載波(及相差90°的兩路正弦波)并行發(fā)送。這樣接收端的偽碼跟蹤環(huán)路一旦鎖定以后，不會受到消息通道的干擾而出現(xiàn)抖動，同時消息通道的載波因?yàn)榕c偽碼通道的載波具備嚴(yán)格的相位差，所以消息的到達(dá)時刻可以由偽碼跟蹤環(huán)路正確測得。

4.3 頻率源校正

本系統(tǒng)選擇的頻率源為輸出100 MHz的恒溫晶振，其頻率準(zhǔn)確度和頻率穩(wěn)定度分別為10－8和10－10/s 量級。理論上兩站的頻率源誤差在幾赫之間，實(shí)際情況中交付的晶體常常由于電老化不充分，在使用過程中會出現(xiàn)頻率漂移，其范圍可以達(dá)到10 Hz量級。工作溫度的不同也會導(dǎo)致晶體是諧振頻率出現(xiàn)偏差，其值通常在幾百皮秒以內(nèi)。

解決方法之一是為頻率源配置壓控電位計(jì)，當(dāng)發(fā)現(xiàn)主從站的時鐘出現(xiàn)頻率偏差時，可以在線通過程序調(diào)整頻率源的頻率，精度可以達(dá)到1 Hz以內(nèi)。根據(jù)式(2)計(jì)算，可以將由頻率源導(dǎo)致的誤差降低一個量級。

綜上所述，通過綜合應(yīng)用以上幾項(xiàng)措施，可以將期望誤差降低到0.2 ns左右的水平。

5 仿真與實(shí)驗(yàn)

本系統(tǒng)的碼片速率為10.23 MHz，采用QPSK 調(diào)制，滾降系數(shù)設(shè)為0.5，信號3 dB帶寬約為15 MHz，接收機(jī)采樣率為100 MHz。參考文獻(xiàn)［8］和［12］的設(shè)置，信道條件選擇AWGN 信道，載噪比設(shè)為55 dB－Hz，這種情況比較接近編隊(duì)飛行的站點(diǎn)間的信道環(huán)境。采用超前滯后結(jié)構(gòu)的2 階偽碼跟蹤環(huán)，碼環(huán)帶寬選擇5 Hz，載波環(huán)帶寬選擇18 Hz，得到仿真結(jié)果如圖2所示。

圖2 偽碼相位跟蹤精度仿真結(jié)果Fig.2 Simulation result of pseudo code phase－tracking accuracy

從仿真結(jié)果可以看出，當(dāng)碼環(huán)鎖定的時候，因?yàn)樵肼暤拇嬖冢镜貍未a相位調(diào)整量在一個固定值附近波動，在5 s的仿真時間內(nèi)，從1 000點(diǎn)后開始統(tǒng)計(jì)(去掉了前面捕獲到鎖定的過渡階段調(diào)整量)，偽碼相位跟蹤精度可達(dá)0.2 ns左右。

內(nèi)場測試中采用信道模擬器來模擬節(jié)點(diǎn)之間的相對運(yùn)動產(chǎn)生的載波多普勒頻偏，載波頻率為2.5 GHz，I 路為偽碼，Q 路為消息，發(fā)射功率為2 W，主站和各個從站的信號通過射頻線纜接入信道模擬器，而不采用天線輻射的方法，信道的等效衰減在標(biāo)定線纜衰減的基礎(chǔ)上，通過可調(diào)衰減器加以模擬。

時間同步精度的測量通過高速示波器觀察主站和兩個從站輸出的秒脈沖之間的時間差來測量，高速示波器的通道誤差約為0.1 ns。通過對測量結(jié)果進(jìn)行統(tǒng)計(jì)處理，可以得到誤差的統(tǒng)計(jì)結(jié)果。

圖3顯示了內(nèi)場的測試系統(tǒng)組成。

圖3 測試系統(tǒng)組成框圖Fig.3 The block diagram of testing system

測試時分別模擬了3 km和6 km距離的情況，相對運(yùn)動速度為30 m/s。對于每個從站點(diǎn)，讀取100組其與主站點(diǎn)的秒脈沖誤差，然后利用公式計(jì)算其時間同步誤差，其中Δti表示從站與主站的秒脈沖誤差，N 表示測試的次數(shù)(100 次以上)。最終的測試結(jié)果如表1所示。

表1 測試結(jié)果Table 1 Test result

由表1可見，本系統(tǒng)的時間同步精度達(dá)到了0.2 ns的時間同步精度，可以滿足編隊(duì)移動分布式定位系統(tǒng)對于時間同步的技術(shù)要求。

6 結(jié) 論

在機(jī)動平臺之間通過偽碼相位測量完成高精度時間同步，是一項(xiàng)應(yīng)用廣泛的技術(shù)。本文結(jié)合工程研制經(jīng)驗(yàn)，分別從頻率源偏差、多普勒效應(yīng)、環(huán)路抖動和上下行鏈路不對稱方面分析了誤差的來源，并量化了誤差的估計(jì)。在此基礎(chǔ)上，通過綜合利用多種方法，可以有效地在站間機(jī)動條件下取得較為理想的時間同步精度。相比經(jīng)典的雙程雙向測量體制，本方案在多個站點(diǎn)組網(wǎng)的條件下具備反應(yīng)速度快、實(shí)現(xiàn)復(fù)雜度低的優(yōu)點(diǎn)。在今后的研究中，如果考慮結(jié)合載波相位測量技術(shù)，預(yù)計(jì)可以取得更高量級的時間同步精度。

［1］ 胡來招.無源定位［M］.北京:國防工業(yè)出版社，2004.HU Lai－ zhao. Passive Locating［M］. Beijing:National Defense Industry Press，2007.(in Chinese)

［2］ 王旭.TDMA 系統(tǒng)目標(biāo)多站定位理論與算法研究［D］.成都:電子科技大學(xué)，2011.WANG Xu. Research on Muti－ target Locating Theories and Algorithms of TDMA System［D］.Chengdu:University of Electronic Science and Technology of China，2011.(in Chinese)

［3］ 王惠南.GPS 導(dǎo)航原理與應(yīng)用［M］.北京:科學(xué)出版社，2003.WANG Hui－nan.Principles and Applications of GPS Navigation［M］.Beijing:Science Press，2003.(in Chinese)

［4］ 陳仕進(jìn).時間同步方法研究［J］. 測控技術(shù)，2005(4):51－53.CHEN Shi－jin.Research on Time synchronization methods［J］.Measurement and Control Technology，2005(4):51－53.(in Chinese)

［5］ 張國琴，吳玉. 基于GPS 校準(zhǔn)晶振的高精度時鐘的設(shè)計(jì)［J］.儀表技術(shù)，2010(4):23－27.ZHANG Guo－ qin，WU Yu. Design of High Precision Clock Based on the GPS Disciplined Oscillator［J］. Instrumental Technique，2010(4):23－27.(in Chinese)

［6］ 杜曉輝，施滸立，張麗榮，等. 一種轉(zhuǎn)發(fā)式衛(wèi)星授時新方法［J］.天文研究與技術(shù)，2012，1(12):38－42.DU Xiao－h(huán)ui，SHI Hu－li，ZHANG Li－rong，et al.A New Timing Method for Satellite Relay System［J］.Astronomical Research and Technology，2012，1 (12):38－ 42.(in Chinese)

［7］ 王平利，程華軍，林思佳，等.基于GNSS CV 的精密時間服務(wù)系統(tǒng)的設(shè)計(jì)與實(shí)現(xiàn)［J］. 時間頻率學(xué)報(bào)，2012(3):11－16.WANG Ping－li，CHENG Hua－jun，LIN Si－jia，et al.Design and Implementation of Precise Time service system based on GNSS CV［J］.Journal of Time and Frequency，2012(3):11－16.(in Chinese)

［8］ 黃波，胡修林.“北斗”2 導(dǎo)航衛(wèi)星星間測距與時間同步技術(shù)［J］.宇航學(xué)報(bào)，2011，6(7):67－71.HUANG Bo，HU Xiu－ lin. Inter Satellite Ranging and Time Synchronization for BD2［J］. Journal of Astronautics，2011，6(7):67－71.(in Chinese)

［9］ Koppang P，Wheeler P.Working application of TWSTT for high precision remote synchronization［C］ //Proceedings of 1998 the Annual IEEE International Frequency Control Symposium.CA，USA:IEEE，1998:27－29.

［10］ Quirk K J，Srinivasan M. Analysis of sampling and quantization effects on the performance of PN code tracking loops［C］//Proceedings of 2002 IEEE International Conference on Communications. Washington DC:IEEE，2002:1480－1484.

［11］ Quirk K J，Srinivasan M. PN code tracking using noncommensurate sampling［J］.IEEE Transactions on Communications，2006，54(10):1845－1856.

［12］ Ke T，Hu X，Liu Y，et a. l The evaluation methods of sampling rate performance in GNSS receiver［C］//Proceedings of the 3rd International Conference on Space Information Technology.Washington DC:IEEE，2009.

［13］ 柯颋，胡修林，劉禹圻，等. GNSS 接收機(jī)中采樣率對時間鑒別力的影響［J］. 宇航學(xué)報(bào)，2010，31(12):2722－2729.KE Ting，HU Xiu－lin，LIU Yu－qi，et al. Influence of Sampling Rate on Time Discrimination in GNSS Receiver［J］. Journal of Astronautics，2010，31(12):2722－2729.(in Chinese)

［14］ 陳中平，徐淑壹，駕英茜，等.小型抗振晶體振蕩器設(shè)計(jì)［J］.電子科技，2014，24(4):76－78.CHEN Zhong－ ping，XU Shu－ yi，JIA Ying－ qian，et al.Design of A Miniature Anti－Vibration Cristal Oscillator［J］. Electronic Science ＆ Technology，2014，24(4):76－78.(in Chinese)

［15］ 王琦，吳斌.航天測控系統(tǒng)中偽碼測距精度分析［J］.測控遙感與導(dǎo)航定位，2009，39(1):39－42.WANG Qi，WU Bin. Analysis on Precision of Pseudo Noise Code Ranging in Space TT＆C Systems［J］.Measurement and Control of Remote Sensing and Navigation，2009，39(1):39－42.(in Chinese)

［16］ 劉寧，劉雪瑞.GPS 軟件接收機(jī)精確載波頻率捕獲算法分析與實(shí)現(xiàn)［J］.全球定位系統(tǒng)，2011，11(2):36－39.LIU Ning，LIU Xue－rui.Analysis and Implementation of Precise Carrier Frequency Capture of GPS Reciever［J］.Global Positioning System，2011，11 (2):36－ 39.(in Chinese)