基于自相關(guān)側(cè)峰消除的北斗B1C快速高精度捕獲算法

許 睿，唐瑞琪，羅 凱，曾慶化，劉建業(yè)

（南京航空航天大學(xué) 自動(dòng)化學(xué)院，南京 211106）

我國北斗衛(wèi)星導(dǎo)航系統(tǒng)（BeiDou Navigation Satellite System，BDS）第三代民用信號(hào)新增B1C 信號(hào)[1]，該信號(hào)采用二進(jìn)制偏移載波（Binary Offset Carrier,BOC）調(diào)制方式[2]，偽隨機(jī)碼周期為10 ms[3]，能夠與現(xiàn)存導(dǎo)航信號(hào)共享頻點(diǎn)，又能提高信號(hào)跟蹤精度與定位精度。

然而與傳統(tǒng)的采用相移鍵控調(diào)制方式的導(dǎo)航信號(hào)（如，GPS L1 C/A 信號(hào)）相比，BOC 調(diào)制信號(hào)自相關(guān)函數(shù)具有多峰特性，會(huì)在接收機(jī)信號(hào)捕獲和跟蹤階段帶來模糊性問題[4]，即信號(hào)自相關(guān)主峰旁存在若干副峰，接收機(jī)有可能捕獲或跟蹤在副峰上，造成誤捕和誤鎖，若不進(jìn)行處理會(huì)進(jìn)一步導(dǎo)致定位結(jié)果不穩(wěn)定或出現(xiàn)較大偏差。自2002年Betz.J.W[4]系統(tǒng)地提出BOC 調(diào)制以來，典型的BOC 調(diào)制信號(hào)多峰消除去模糊算法可分為三類。第一類是邊帶濾波法，也稱為BPSK-Like 算法[5]，通過頻域?yàn)V波將BOC 調(diào)制信號(hào)頻譜上兩個(gè)主波瓣提出，再作為類BPSK 信號(hào)頻譜處理，這種方法會(huì)降低BOC 調(diào)制帶來性能優(yōu)勢，且需要額外的濾波器資源。第二類為峰跳法，即Bump-Jump（BJ）方法[6]，該方法需要額外的相關(guān)器通過對(duì)比第一副峰與主峰的大小來判定捕獲的點(diǎn)是否在主峰上，需要經(jīng)過多次采樣計(jì)算，新增的相關(guān)器涉及大量的加乘運(yùn)算增加硬件開銷，且在信噪比較低的情況下會(huì)出現(xiàn)判決錯(cuò)誤的現(xiàn)象。第三類是自相關(guān)副峰消除技術(shù)[7,8]，通過設(shè)計(jì)不同的相關(guān)函數(shù)并按照一定的組合方式進(jìn)行合成，從而消除主峰之外的所有副峰實(shí)現(xiàn)無模糊捕獲，不需增加大量運(yùn)算和新的硬件，靈活性更大。

除多峰性問題，北斗B1C 信號(hào)采用10 ms 周期的偽隨機(jī)碼，碼周期的增加對(duì)應(yīng)信號(hào)跟蹤階段的積分時(shí)間的增長，這有助于減少噪聲，提高跟蹤精度和偽距估計(jì)精度，但同時(shí)跟蹤環(huán)路的帶寬需要相應(yīng)地降低，這又會(huì)導(dǎo)致跟蹤環(huán)路牽引范圍變窄，因此，跟蹤環(huán)路初始值就需要更準(zhǔn)確，這對(duì)提供跟蹤環(huán)路初值的捕獲階段提出了更高精度的要求。此外，較長的碼周期也增加了捕獲運(yùn)算量與耗時(shí)。對(duì)于北斗三號(hào)B1C 接收機(jī)而言，要實(shí)現(xiàn)信號(hào)的快速準(zhǔn)確跟蹤與定位，就必須先進(jìn)行快速捕獲，并得到精度足夠高的載波頻率與碼相位初值。接收機(jī)通常采用先粗捕獲再精捕獲的方法[9]平衡捕獲的精度與速度，粗捕利用FFT 算法實(shí)現(xiàn)對(duì)頻率和碼相位搜索的并行計(jì)算得到粗略的多普勒頻率與碼相位，精捕則可通過對(duì)FFT 輸出序列進(jìn)行SINC 擬合來獲得更精確的多普勒頻率估值[10]。但對(duì)于采用BOC 調(diào)制的信號(hào)來說，由于副載波的存在會(huì)影響SINC 擬合的效果，導(dǎo)致精捕獲得到的頻率誤差較大，從而會(huì)使得信號(hào)跟蹤失敗。鑒于此，本文針對(duì)北斗B1C信號(hào)設(shè)計(jì)了一種快速高精度捕獲算法，先是基于降采樣和多級(jí)搜索的方式，通過減少數(shù)據(jù)點(diǎn)，快速實(shí)現(xiàn)了對(duì)信號(hào)的粗捕獲并將對(duì)多普勒頻移的估計(jì)誤差縮小到一定范圍，大大減少了信號(hào)捕獲的計(jì)算量，有效提高了B1C 信號(hào)的捕獲速度；再通過多項(xiàng)式曲線擬合來估計(jì)頻率相對(duì)于相關(guān)值的峰值點(diǎn)，實(shí)現(xiàn)了多普勒頻移的高精度捕獲，并有效避免了副載波對(duì)精捕獲的影響。最后基于真實(shí)信號(hào)的實(shí)驗(yàn)結(jié)果驗(yàn)證了該算法對(duì)捕獲精度和速度的提升。

1 B1C 信號(hào)與ASPeCT 算法

1.1 B1C 信號(hào)的數(shù)學(xué)模型

北斗B1C 信號(hào)的載波頻率為1575.42 MHz，與GPS L1 和Galileo E1 共享頻點(diǎn)，帶寬為32.736 MHz，包含兩個(gè)正交分量：數(shù)據(jù)分量與導(dǎo)頻分量。數(shù)據(jù)分量由導(dǎo)航電文和測距碼經(jīng)副載波調(diào)制產(chǎn)生，調(diào)制方式為正弦BOC（1,1）調(diào)制；導(dǎo)頻分量由測距碼經(jīng)副載波調(diào)制產(chǎn)生，采用QMBOC（6,1,4/33）調(diào)制，其功率是數(shù)據(jù)分量功率的3 倍。除使用子載波外，導(dǎo)頻分量還使用周期為18 s、碼長度為1800 碼片的子碼來調(diào)制測距主碼，每個(gè)子碼碼元持續(xù)時(shí)間為10 ms，與測距主碼的周期相等，也與數(shù)據(jù)分量中導(dǎo)航電文的一位（bit）所持續(xù)時(shí)間相同。導(dǎo)頻分量采用的QMBOC 調(diào)制方式將BOC（1,1）分量和BOC（6,1）分量按照29:4 的功率比調(diào)制在相互正交的副載波上加以實(shí)現(xiàn)。整個(gè)B1C信號(hào)實(shí)際包含了三個(gè)實(shí)分量[4]，如式(1)所示。

B1C 信號(hào)的BOC（1,1）為低頻分量，是信號(hào)的窄帶部分；BOC（6,1）為高頻分量，是信號(hào)的寬帶部分。BOC（1,1）分量包含完整的信號(hào)結(jié)構(gòu)，功率也較BOC（6,1）分量更大，因此常規(guī)處理GPS L1 C/A 的窄帶接收設(shè)備也可以接收到北斗B1C 信號(hào)。以一顆星信號(hào)為例，北斗B1C 信號(hào)在窄帶接收機(jī)中的中頻信號(hào)可寫為：

根據(jù)式(2)，窄帶接收機(jī)接收到的北斗B1C 信號(hào)可近似認(rèn)為僅包含BOC（1,1）調(diào)制信號(hào)部分。其自相關(guān)函數(shù)表示為：

式中：τ表示以碼片為單位的碼延時(shí)，代表相關(guān)函數(shù)，代表三角函數(shù)。如圖1所示，其自相關(guān)峰主峰較為尖銳，具有良好的測距性能，但存在副峰，會(huì)導(dǎo)致誤捕和誤鎖，需要對(duì)副峰進(jìn)行處理以消除其不利影響。

圖1 BOC(1,1)調(diào)制方式下的自相關(guān)函數(shù)Fig.1 Autocorrelation function in BOC(1,1)modulation

1.2 自相關(guān)側(cè)峰消除捕獲算法

對(duì)于BOC（n，n）類調(diào)制信號(hào)可以使用自相關(guān)側(cè)峰消除（ Autocorrelation Side-Peak Cancellation Technique,ASPeCT）算法消除副峰，其原理是：首先將本地生成的PRN 碼（測距碼）分成兩組，一組使用副載波進(jìn)行調(diào)制得到BOC 信號(hào)，另一組不作處理為原始PRN 信號(hào)，然后將這兩組本地碼分別與經(jīng)過混頻器處理后的信號(hào)作相關(guān)運(yùn)算，對(duì)于BOC（1,1）調(diào)制信號(hào)而言，得到的信號(hào)自相關(guān)函數(shù)如式(3)所示，其與測距碼的互相關(guān)函數(shù)表示為：

再將兩個(gè)相關(guān)函數(shù)進(jìn)一步組合，使用二者的平方差構(gòu)建出一個(gè)新的相關(guān)函數(shù)，即：

圖2 BOC(1,1)信號(hào)的ASPeCT 算法捕獲原理Fig.2 Acquisition principle of ASPeCT algorithm for BOC(1,1)signal

當(dāng)接收的中頻信號(hào)分別與i支路和q支路上的正弦和余弦載波信號(hào)混頻后，將得到的復(fù)數(shù)形式的混頻結(jié)果進(jìn)行FFT 變換，然后將其分別與BOC(1,1)和PRN本地碼FFT 變換的共軛值相乘，再分別作IFFT 變換后取平方差，得到的值即為ASPeCT 算法所構(gòu)造的相關(guān)值，若峰值超過捕獲門限值，則接收機(jī)捕獲到了信號(hào)，并得到該信號(hào)的頻率和碼相位兩個(gè)參數(shù)值。

2 快速高精度捕獲算法

基于ASPeCT 算法可以解決BOC 信號(hào)自相關(guān)多峰特性帶來的誤捕問題，從而能順利實(shí)現(xiàn)北斗B1C 信號(hào)的粗捕獲過程。但粗捕得到的多普勒頻移估計(jì)誤差較大。若要使得接收機(jī)能夠?qū)崿F(xiàn)信號(hào)的快速穩(wěn)定跟蹤，則需要對(duì)載波頻率進(jìn)行較為精確的估計(jì)，即進(jìn)行精捕獲。由于B1C 信號(hào)中副載波調(diào)制的存在，影響了傳統(tǒng)精捕獲算法中SINC 曲線擬合的效果，導(dǎo)致精捕獲得到的頻率誤差較大，從而可能會(huì)導(dǎo)致跟蹤失敗。此外，北斗B1C 信號(hào)的一個(gè)偽碼周期為10 ms，包含10230個(gè)碼片，與傳統(tǒng)GPS 中C/A 碼每周期1 ms 對(duì)應(yīng)1023個(gè)碼片相比，從周期長度和序列長度上都增加到10倍。這給信號(hào)捕獲帶來了更高的計(jì)算量，相同條件下完成衛(wèi)星信號(hào)捕獲勢必要花費(fèi)更多時(shí)間。因此，需要在捕獲算法的精度與效率上同時(shí)改進(jìn)。

鑒于此，本文設(shè)計(jì)了一種快速高精度捕獲算法，首先基于多級(jí)搜索和降采樣策略快速將信號(hào)的頻率捕獲縮小到一個(gè)較小的范圍，再通過多項(xiàng)式曲線擬合算法，實(shí)現(xiàn)對(duì)B1C 信號(hào)多普勒頻移的高精度捕獲。該算法的整體邏輯思路如圖3所示。

圖3 針對(duì)B1C 信號(hào)的快速高精度捕獲算法框圖Fig.3 Block diagram of fast and high-precision acquisition algorithm for B1C signal

2.1 多級(jí)搜索控制

多級(jí)搜索的目的是將對(duì)信號(hào)頻率的搜索快速縮小到一個(gè)較小的范圍之內(nèi)，以便于使用曲線擬合的方法對(duì)頻率進(jìn)行更加精確地估計(jì)。多級(jí)搜索的思路如圖4所示。

圖4 頻率多級(jí)搜索示意圖Fig.4 Schematic diagram of frequency multilevel search

多級(jí)搜索的級(jí)數(shù)與每級(jí)搜索的步長共同影響捕獲算法的運(yùn)算量，為了達(dá)到更好的曲線擬合效果，最終搜索步長設(shè)為5 Hz，以達(dá)到更為精確的頻率估計(jì)值。而一級(jí)搜索步長則需要在確保頻率估計(jì)值落在正確范圍之內(nèi)的基礎(chǔ)上盡量選擇較大的值，以此減少一級(jí)搜索的次數(shù)。為了降低漏警發(fā)生的概率，接收機(jī)通常應(yīng)當(dāng)對(duì)于頻率誤差絕對(duì)值的設(shè)置不得超過表示相干積分時(shí)間，從而得到頻率搜索步長的取值范圍為：

北斗B1C 信號(hào)的偽碼周期為10 ms，這就意味著該信號(hào)在接收機(jī)捕獲環(huán)節(jié)的相干積分時(shí)間最少應(yīng)為偽碼單周期的10 ms，從而得到應(yīng)當(dāng)不超過88.6 Hz，即當(dāng)頻率誤差大于44.3 Hz 時(shí)，漏警事件的概率將會(huì)逐漸提升。故可設(shè)置為80 Hz，相應(yīng)的級(jí)數(shù)與搜索次數(shù)對(duì)應(yīng)關(guān)系如表1所示，搜索范圍為7.4 kHz。

表1 不同步長下的搜索次數(shù)統(tǒng)計(jì)Tab.1 Statistics of search times under different step size

由表1 可得，三級(jí)搜索且二級(jí)搜索步長在30 Hz和20 Hz 時(shí)的總搜索次數(shù)分別為105 次和104 次，而20 Hz 時(shí)三級(jí)搜索的搜索次數(shù)只有5 次，考慮到精捕階段希望有較多數(shù)據(jù)點(diǎn)以供多項(xiàng)式擬合使用，則30 Hz 的二級(jí)搜索步長更為合適，故三級(jí)搜索的步長應(yīng)該分別設(shè)置為80 Hz、30 Hz 和5 Hz。

基于此多級(jí)搜索策略，可以進(jìn)一步使用時(shí)頻域結(jié)合搜索的方法優(yōu)化算法的結(jié)構(gòu)。考慮到三級(jí)搜索每一次均進(jìn)行了頻率與碼相位在頻域的并行搜索，而三次搜索的碼相位基本一致，即一級(jí)搜索的碼相位已較為精確，之后的碼相位搜索工作是重復(fù)的，所以理論上可以省去后兩次的碼相位搜索來減少計(jì)算量。在具體實(shí)現(xiàn)上，將初次搜索得到的碼相位作為后兩次搜索讀入信號(hào)的起始點(diǎn)，將后兩次的頻域并行搜索改為時(shí)域的頻率單獨(dú)搜索，以便極大地減少了后兩次搜索的計(jì)算量。

2.2 降采樣策略

基于降采樣的信號(hào)快速捕獲算法的思路是使用數(shù)據(jù)預(yù)處理的方法，通過縮短輸入數(shù)據(jù)與本地碼的長度來減少計(jì)算量，從而使得信號(hào)捕獲速度加快。首先在本地偽碼生成環(huán)節(jié)，按照前端設(shè)備的采樣頻率生成偽碼數(shù)據(jù)后以固定倍數(shù)進(jìn)行降采樣，即按照每隔固定長度數(shù)據(jù)量取一個(gè)數(shù)據(jù)點(diǎn)的規(guī)律，利用取出的數(shù)據(jù)點(diǎn)生成新的序列，從而實(shí)現(xiàn)偽碼數(shù)據(jù)的降采樣處理。同樣地，在中頻數(shù)據(jù)讀入環(huán)節(jié)也要按照相同的規(guī)律進(jìn)行數(shù)據(jù)處理，從而保證中頻數(shù)據(jù)和本地碼的對(duì)應(yīng)關(guān)系，再利用降采樣后的數(shù)據(jù)進(jìn)行信號(hào)捕獲時(shí)就可以節(jié)省大量的計(jì)算。

降采樣實(shí)質(zhì)是對(duì)已采集信號(hào)的二次抽樣，抽樣后不僅信號(hào)的采樣頻率降低，實(shí)際中頻頻率也對(duì)應(yīng)改變。研究中所采用的前端設(shè)備，采樣頻率16.368 MHz，中頻頻率4.092 MHz，非正交采樣。當(dāng)二次抽樣倍數(shù)為偶數(shù)倍時(shí)，降采樣后的數(shù)據(jù)中頻降至零中頻，非正交采樣導(dǎo)致頻譜混疊，則捕獲算法無法獲得正確的多普勒頻率。實(shí)驗(yàn)也表明，在對(duì)數(shù)據(jù)進(jìn)行二倍降采樣或者四倍降采樣之后，由于偶數(shù)倍降采樣處理后數(shù)據(jù)的固定特征點(diǎn)損失，捕獲結(jié)果中會(huì)出現(xiàn)兩個(gè)可能的峰值，且二者的載波頻率關(guān)于中心頻率對(duì)稱，即多普勒互為相反數(shù)。這就帶來了很大的誤捕風(fēng)險(xiǎn)。因此需要采用奇數(shù)倍的降采樣倍數(shù)，再考慮到采樣頻率過低會(huì)出現(xiàn)信號(hào)失真的問題，故使用三倍降采樣進(jìn)行數(shù)據(jù)預(yù)處理較為合適。圖5 展示的就是對(duì)本地偽碼進(jìn)行三倍降采樣的思路。相應(yīng)地，接收機(jī)讀入的中頻數(shù)據(jù)也要進(jìn)行相同規(guī)律的處理。

圖5 中，變量N 代表的是在前端采樣頻率下測距碼一個(gè)完整周期的采樣點(diǎn)數(shù)。式(2)中的中頻信號(hào)離散化后得到處理，新的中頻信號(hào)為其中通過三倍降采樣 其中m=1,2,3...N/3。降采樣后的數(shù)據(jù)序列長度縮短到原始長度的三分之一。

圖5 偽碼降采樣示意圖Fig.5 Schematic diagram of PRN code downsampling

結(jié)合2.1 節(jié)的多級(jí)與時(shí)頻域切換搜索策略，接收機(jī)輸入的中頻信號(hào)和本地偽碼經(jīng)過同樣的降采樣處理后，首先是粗捕獲通過FFT 在頻域中對(duì)頻率和碼相位進(jìn)行并行搜索，然后控制器記錄一級(jí)搜索得到的碼相位并將相關(guān)器切換到時(shí)域，同時(shí)控制本地載波的頻率搜索范圍和輸入信號(hào)的起點(diǎn)，只進(jìn)行二級(jí)與三級(jí)的頻率搜索。這樣使得整個(gè)捕獲過程中的計(jì)算量大大減少，可將捕獲速度提高到三倍左右。具體思路如圖3所示。

在整個(gè)捕獲過程中，最耗費(fèi)時(shí)間的是搜索過程中的FFT 運(yùn)算。假設(shè)中頻數(shù)據(jù)一個(gè)完整偽碼周期的采樣點(diǎn)數(shù)為N，一個(gè)數(shù)據(jù)點(diǎn)的FFT 運(yùn)算量為M，精捕獲算法和快速精捕獲算法的FFT 運(yùn)算量如表2所示。

表2 運(yùn)算量對(duì)比Tab.2 Comparison of calculation quantity

由表2 可以得到，與單獨(dú)的精捕獲算法相比，使用降采樣與時(shí)頻域切換操作的快速精捕獲算法將占據(jù)主要計(jì)算量的FFT 運(yùn)算相對(duì)減少了70.16%，考慮到算法中的其他運(yùn)算，整體計(jì)算量減少比例會(huì)略低于該值。需要指出的是，降采樣處理可以有效降低運(yùn)算量，但也使得信號(hào)能量有少量損失，使得捕獲靈敏度略有下降。

2.3 多項(xiàng)式曲線擬合

為增加后續(xù)跟蹤環(huán)路的跟蹤成功率，可以進(jìn)一步提高多普勒頻率估計(jì)精度。2.2 節(jié)中的多級(jí)搜索策略可以將載波頻率估計(jì)縮小到一個(gè)較為精確的范圍，而由于相關(guān)結(jié)果包含SINC 項(xiàng)，在小范圍頻率區(qū)域內(nèi)可簡化近似于拋物線，因此可采用一個(gè)二階多項(xiàng)式曲線擬合，以提高頻率估計(jì)精度，在擬合式選取擬合點(diǎn)數(shù)不易過多，以避免邊界效應(yīng)，本文中選取5 個(gè)點(diǎn)。

其中，

最后根據(jù)獲得的擬合拋物線方程，求最大相關(guān)積分值，即曲線頂點(diǎn)所對(duì)應(yīng)的頻率f，該頻率即為捕獲獲得的載波頻率。以一組PRN40 通道的數(shù)據(jù)為例，繪制擬合前后的相關(guān)曲線如圖6所示。

圖6 二次多項(xiàng)式擬合多普勒Fig.6 Doppler frequency estimated by second-order polynomial fitting

由圖6 可以看到擬合曲線與數(shù)據(jù)點(diǎn)連線較為貼切，很好地反映了相關(guān)值的變化趨勢，圖中標(biāo)注的拋物線頂點(diǎn)相對(duì)應(yīng)的頻點(diǎn)即為多普勒頻移的估計(jì)值。

3 實(shí)驗(yàn)驗(yàn)證與分析

為測試所提方法，使用前端采集器（iP-Solutions E2 GPS/Galileo L1 Front end）采集北斗B1C 真實(shí)信號(hào)進(jìn)行測試，數(shù)據(jù)采集地點(diǎn)在無遮擋的樓頂，前端與采集數(shù)據(jù)主要參數(shù)如表3所示。

表3 采集的B1C 信號(hào)主要參數(shù)Tab.3 Main parameters of collected B1C signal

3.1 副峰消除效果分析

圖7所示為ASPeCT 捕獲算法解算B1C 真實(shí)信號(hào)（ PRN40 ） 獲得的碼延時(shí)-多普勒頻率圖（Delay-Doppler Map,DDM），以及多普勒頻率為-219 Hz 時(shí)的相關(guān)值變化曲線。

圖7 BOC(1,1)信號(hào)ASPeCT 算法捕獲結(jié)果（PRN40）Fig.7 Acquisition result of ASPeCT algorithm for BOC(1,1)signal(PRN40)

圖7(a)顯示該算法可以完成真實(shí)信號(hào)的捕獲，與不使用ASPeCT 算法結(jié)果對(duì)比，如圖7(b)，副峰得到了很大程度上的削減，能夠更有效地避免誤捕。

表4 展示了采用ASPeCT 算法前后，從該數(shù)據(jù)中捕獲到的6 顆衛(wèi)星DDM 中的副峰主峰比。

表4 采用ASPeCT 算法前后峰值對(duì)比Tab.4 Comparison of peak values before and after using ASPeCT algorithm

數(shù)據(jù)表明，北斗B1C 真實(shí)信號(hào)自相關(guān)曲線的副峰與主峰的平均比值約為69.30%，經(jīng)過ASPeCT 算法處理之后該比值約降至9.22%，整體的副峰消除效果可以達(dá)到86.70%，消除了信號(hào)誤捕的風(fēng)險(xiǎn)。

3.2 捕獲速度對(duì)比

表5所示為采用本文所提快速精捕獲算法和普通精捕獲算法耗時(shí)對(duì)比表，捕獲耗時(shí)使用MATLAB 統(tǒng)計(jì)，由于北斗衛(wèi)星數(shù)量較多，僅列出捕獲到衛(wèi)星信號(hào)的通道耗時(shí)。

表5 快速精捕獲算法速度對(duì)比表Tab.5 Speed comparison of fast fine acquisition algorithm

從表5 可以看出快速捕獲算法耗時(shí)顯著降低，普通精捕獲算法單通道平均捕獲用時(shí)8.68 s，而使用快速捕獲算法后單通道平均用時(shí)3.11 s，實(shí)際捕獲時(shí)間相對(duì)縮短了64.17%，較2.2 節(jié)分析的減少70.16%，F(xiàn)FT計(jì)算量略有降低，但考慮捕獲還包含擬合等其他運(yùn)算量，實(shí)際改進(jìn)近2/3，仍然顯著。

3.3 多普勒頻率捕獲精度分析

為分析快速精捕獲算法獲得的多普勒頻率精度，利用鎖相環(huán)輸出的載波多頻率頻率擬合多普勒頻率曲線，根據(jù)多普勒近似線性的變化規(guī)律，可以通過直線擬合來近似得到多普勒頻率實(shí)際曲線，如圖8所示。再根據(jù)擬合多普勒頻率曲線獲得的跟蹤初始時(shí)刻多普勒頻率值，即0 s 時(shí)刻的多普勒頻率，作為捕獲時(shí)刻的精確多普勒頻率估計(jì)，則多頻率頻率的捕獲誤差即為捕獲的多頻率頻率與跟蹤初始時(shí)刻多普勒頻率之差。圖8所示為PRN40通道衛(wèi)星多普勒頻率跟蹤結(jié)果，通過上述方法可得信號(hào)初始時(shí)刻的多普勒頻率估計(jì)誤差為0.63 Hz。

圖8 多普勒直線擬合Fig.8 Doppler line fitting

根據(jù)該方法，獲得精捕獲前后各衛(wèi)星信號(hào)多普勒頻率誤差，結(jié)果如表6所示。

表6 快速精捕獲算法精度對(duì)比表Tab.6 Precision comparison of fast fine acquisition algorithm

從表6 可以看出，當(dāng)僅進(jìn)行粗捕獲時(shí)，頻率搜索步長為80 Hz，以滿足捕獲最大步長要求，可以正確捕獲對(duì)應(yīng)的衛(wèi)星信號(hào)，但由于捕獲得到的多普勒誤差范圍相對(duì)較大，雖然滿足捕獲步長的最低要求，卻無法確保跟蹤環(huán)路的初值精度，發(fā)生跟蹤失敗的可能性較高；當(dāng)使用本文提出的快速精捕獲算法時(shí)，上述測試用例均實(shí)現(xiàn)了成功跟蹤，且多普勒頻率的捕獲誤差約0.46 Hz。

4 結(jié)論

本文針對(duì)北斗B1C 信號(hào)，在基于ASPeCT 算法消除副峰影響的基礎(chǔ)上，設(shè)計(jì)了一種快速精捕獲算法策略，該算法通過多級(jí)搜索和降采樣策略實(shí)現(xiàn)了B1C 信號(hào)的快速捕獲，大大減少了捕獲過程中的FFT 計(jì)算量，并快速將頻率估計(jì)縮小到一個(gè)較小范圍，再結(jié)合多項(xiàng)式曲線擬合算法，實(shí)現(xiàn)了對(duì)多普勒的高精度捕獲。經(jīng)過真實(shí)信號(hào)測試，該快速精捕獲算法將平均捕獲時(shí)間相對(duì)縮短了近2/3，多普勒頻率的平均估計(jì)誤差在1Hz以內(nèi)，保證了后續(xù)信號(hào)跟蹤的成功率。