GPS軟件接收機基帶信號處理算法及驗證*

劉瑞華 孔月明 蔚保國 邢兆棟

1.中國民航大學(xué)電子信息工程學(xué)院，天津 300300 2.中國電子科技集團公司第54研究所，石家莊 050081

?

GPS軟件接收機基帶信號處理算法及驗證*

劉瑞華1孔月明1蔚保國2邢兆棟2

1.中國民航大學(xué)電子信息工程學(xué)院，天津 300300 2.中國電子科技集團公司第54研究所，石家莊 050081

主要研究基于數(shù)據(jù)塊的GPS軟件接收機基帶信號處理流程和算法，包括對載波和偽碼的捕獲和跟蹤過程，并解調(diào)出導(dǎo)航電文數(shù)據(jù)。利用實際接收到的衛(wèi)星信號數(shù)據(jù)，結(jié)合MATLAB軟件平臺，對GPS軟件接收機基帶信號處理過程的各個環(huán)節(jié)進行仿真驗證。結(jié)果表明，所研究的基帶信號處理算法能較好的捕獲衛(wèi)星信號，可準(zhǔn)確跟蹤載波信號和偽碼序列，能正確解調(diào)出導(dǎo)航數(shù)據(jù)。同時驗證了該算法在軟件接收機系統(tǒng)中應(yīng)用的可行性，并為開發(fā)實時GPS接收機奠定了基礎(chǔ)。 關(guān)鍵詞 GPS軟件接收機；捕獲；跟蹤；仿真驗證；可行性

隨著現(xiàn)代全球?qū)Ш叫l(wèi)星定位系統(tǒng)的迅速發(fā)展，設(shè)計并實現(xiàn)能靈活充分利用所有導(dǎo)航衛(wèi)星信號的下一代GNSS接收機成為一項極具挑戰(zhàn)性的工程[1]。傳統(tǒng)的GNSS接收機主要基于專用集成電路結(jié)構(gòu)，這種接收機設(shè)計較為固定，一般不可進行重置，用戶很難改變接收機各類參數(shù)以適應(yīng)不同導(dǎo)航信號處理的需求，同時也給基帶信號處理新算法的實現(xiàn)、測試和分析帶來不便[1-3]。而基于軟件無線電思想的GNSS軟件接收機，將硬件作為無線通信基本平臺，用軟件來實現(xiàn)盡可能多的用戶通信功能，只需進行軟件修改便可對接收機實現(xiàn)升級，且這種更改設(shè)置還可以用于處理不同體制的衛(wèi)星信號。

與傳統(tǒng)的GNSS接收機相比，軟件接收機系統(tǒng)設(shè)計簡潔，開發(fā)速度快，成本更低，且有高度可配置性，更具方便性和靈活性，對新一代導(dǎo)航接收機的開發(fā)有著重要意義[2]。本文以接收處理L1波段C/A碼信號的單頻GPS軟件接收機為例，基于GPS中頻信號采樣實際數(shù)據(jù)，并對數(shù)字中頻GPS信號的捕獲和跟蹤過程及其算法進行了分析研究，同時結(jié)合MATLAB軟件，仿真基帶信號數(shù)據(jù)處理流程，最后解調(diào)出導(dǎo)航數(shù)據(jù)，驗證了捕獲和跟蹤算法的正確性及其在接收機系統(tǒng)中應(yīng)用的可行性。

1 GPS信號捕獲算法的分析與仿真

接收機的捕獲是一個粗同步過程，能粗略估計出接收信號的 C/A 碼相位和載波頻率，然后傳遞給跟蹤過程，作為跟蹤環(huán)路的初始化參數(shù)[4-5]。

1.1 GPS信號捕獲算法的原理分析

本文研究的捕獲方法是基于頻域的并行碼相位搜索捕獲算法，通過傅里葉變換實現(xiàn)對C/A碼相位的并行搜索。只需對多個不同的載波頻率點進行捕獲運算，從而將搜索空間消減，加快信號搜索速度[4,6]。并行碼相位搜索捕獲算法的工作原理如圖1所示。

圖1 并行碼相位搜索捕獲算法框圖

采集到包含多顆衛(wèi)星的信號可表示為：

sr(n)=s1(n)+s2(n)+…+sk(n)+e(n)

(1)

其中,sr(n)是采集到的衛(wèi)星信號，sk(n)代表當(dāng)前可見第k顆GPS衛(wèi)星信號，e(n)是高斯白噪聲。對衛(wèi)星信號的搜索就是從這一混合信號中確認(rèn)出可見衛(wèi)星及其信號的C/A碼初始相位和載波粗頻率。

數(shù)字中頻信號表達式為：

sk(n)=Ck(n)Dk(n)cos[(wIF+wd)n+φ]+e(n)

(2)

其中，Ck(n)是第k顆衛(wèi)星的C/A碼序列，Dk(n) 是導(dǎo)航數(shù)據(jù)序列，wIF是載波頻率經(jīng)下變頻后的中頻，wd是載波多普勒頻移，φ是載波初始相位。

輸入的數(shù)字中頻信號與本地載波相乘，分別得到I，Q兩路信號：

(3)

(4)

其中，we是輸入信號與本地載波之間的頻率差，ei(n)，eq(n)是e(n)與本地載波相乘后進行積分的結(jié)果。

由I，Q兩路信號組合得到復(fù)信號，進行歸一化可得：

x(n)=I(n)-jQ(n)

(5)

對x(n)，C(n)分別進行傅里葉變換，結(jié)果為X(k)和C(k)，則：

(6)

式中,m=0,1,2,…,N-1，C*(k)是C(k)的復(fù)共軛。對z(m)取?？傻茫?/p>

(7)

式(7)可用于得到輸入信號與本地信號之間的相關(guān)關(guān)系，最后檢測出最大相關(guān)值，來判斷衛(wèi)星信號是否成功捕獲。對某一頻點的碼相位進行并行搜索，只需要完成2次FFT，一次逐點相乘和一次IFFT運算[7]。

1.2GPS信號捕獲算法的仿真驗證

本文使用美國科羅拉多大學(xué)利用射頻采集前端獲得的數(shù)據(jù)進行MATLAB仿真，驗證GPS軟件接收機中的并行碼相位搜索捕獲算法。

實際采得的GPS數(shù)據(jù)重要參數(shù)為：GPS信號帶寬為19.096MHz；采樣頻率為38.192MHz；中頻頻率為9.548MHz；數(shù)據(jù)時長為37s；以8位的二進制形式存儲。文中采用接收處理L1波段C/A碼信號的單頻GPS軟件接收機：包含8個通道；可捕獲32顆衛(wèi)星；載波頻率搜索范圍為中頻頻率±7KHz；載波頻率的搜索步長取500Hz；C/A碼相位搜索范圍為0～38192；碼相位的搜索步長取1個采樣點；判決捕獲成功與否的相關(guān)峰比值門限設(shè)為2.5。

MATLAB軟件仿真得到32顆衛(wèi)星的捕獲結(jié)果如圖2所示。

圖2 所有衛(wèi)星信號捕獲結(jié)果

該GPS軟件接收機捕獲到PRN21，PRN22，PRN15，PRN18，PRN26，PRN06，PRN09和PRN03這8顆衛(wèi)星的信號，它們的捕獲峰值都大于預(yù)設(shè)門限2.5，相應(yīng)的相關(guān)峰值小于2.5的衛(wèi)星信號沒有被捕獲到。捕獲結(jié)果包括衛(wèi)星對應(yīng)的PRN號、載波頻率、碼相位和多普勒頻移，并按照衛(wèi)星信號強度列于表1。

表1 8個通道的捕獲結(jié)果

結(jié)合表1和圖2可以看出PRN21的衛(wèi)星信號最強，所以選取衛(wèi)星PRN21為例進行分析研究。圖3和4是使用并行碼相位搜索捕獲算法成功捕獲PRN21衛(wèi)星信號的結(jié)果。

圖3 PRN21衛(wèi)星二維捕獲結(jié)果

圖4 PRN21衛(wèi)星三維捕獲結(jié)果

圖3和圖4分別為PRN21衛(wèi)星信號的二維和三維捕獲結(jié)果。循環(huán)互相關(guān)運算完成后，成功捕獲到PRN21衛(wèi)星，其峰值位于頻率9.547MHz和碼相位13404處。

從MATLAB仿真結(jié)果可以看出，GPS軟件接收機利用并行碼相位搜索捕獲算法能成功捕獲到實測數(shù)據(jù)中所包含的8顆衛(wèi)星信號，并粗略估計出各個衛(wèi)星信號對應(yīng)的載波頻率和碼相位參數(shù)，這充分驗證了此算法的正確性及其在GPS接收機系統(tǒng)中應(yīng)用的可行性。

2 GPS信號跟蹤算法的分析與仿真

GPS信號捕獲僅提供了對載波頻率和碼相位參數(shù)的粗略估計。跟蹤的主要目的是使這些估計值精確化，并保持跟蹤，然后從跟蹤到的衛(wèi)星信號中提取各觀測量，解調(diào)出導(dǎo)航數(shù)據(jù)[6]。跟蹤的實質(zhì)就是為了穩(wěn)定的跟蹤信號而采取的一種對環(huán)路參數(shù)的動態(tài)調(diào)整策略。信號跟蹤有2個方面：1)實現(xiàn)對GPS信號中載波分量的跟蹤；2)實現(xiàn)對偽碼分量的跟蹤。

2.1 GPS信號跟蹤算法的原理分析

在GPS軟件接收機內(nèi)部必須有2個跟蹤環(huán)：碼跟蹤環(huán)(DLL)和載波跟蹤環(huán)(PLL)，它們的連接方式如圖5所示，載波跟蹤環(huán)中用來剝離C/A碼的本地碼來自碼跟蹤環(huán)，碼跟蹤環(huán)中用來剝離載波的兩路本地載波信號來自載波跟蹤環(huán)[6,8]。

圖5 組合的DLL和PLL跟蹤環(huán)框圖

PLL與DLL優(yōu)化組合的跟蹤環(huán)工作原理為：

1)載波跟蹤：采用Costas環(huán)實現(xiàn)載波跟蹤環(huán)路，這種鎖相環(huán)對導(dǎo)航數(shù)據(jù)引起的相位翻轉(zhuǎn)不敏感，其目的是將所有的信號能量保留在I支路[9]。本地碼相位精確對齊時，I，Q支路混頻結(jié)果分別為：

(8)

(9)

其中，φe是輸入信號與本地載波的相位差?；祛l后的信號分別經(jīng)過低通濾波器輸出結(jié)果為：

(10)

(11)

輸入信號與本地載波的相位差可由式(12)求得：

(12)

(13)

相位誤差通過載波環(huán)濾波器濾除噪聲后，輸出給數(shù)控振蕩器(NCO)，調(diào)節(jié)振蕩器頻率，使其與輸入信號的頻率保持一致。利用該方式，實現(xiàn)對輸入信號載波的本地復(fù)現(xiàn)，解調(diào)碼跟蹤環(huán)的輸出信號，得到導(dǎo)航電文數(shù)據(jù)。

2)碼跟蹤：碼相位跟蹤的實現(xiàn)通常采用非相干延遲鎖定環(huán)，其運行和控制過程與載波環(huán)很相似，區(qū)別在于DLL需要三路本地偽碼為輸入信號的跟蹤提供延遲依據(jù)，用來判斷本地偽碼與輸入信號之間相關(guān)結(jié)果的最大值。

當(dāng)碼相位需要調(diào)整時，碼相位鑒別器利用E，L支路上的相關(guān)結(jié)果計算出碼相位誤差，經(jīng)濾波后輸入到PRN碼發(fā)生器，調(diào)節(jié)本地偽碼的頻率和相位，使其與輸入信號的偽碼保持對齊。碼相位鑒別器通常使用超前-滯后能量差相位鑒別器，鑒別因子D表示為：

(14)

其中，IE，QE，IL和QL為圖5所示的6個相關(guān)器中的4個輸出。DLL對載波相位不敏感，當(dāng)本地載波存在相位誤差時，可同時利用I，Q兩路信號判斷正確的相關(guān)輸出，從而保證C/A碼相位的正確跟蹤。

2.2GPS信號跟蹤算法的仿真驗證

采用PLL與DLL優(yōu)化組合的跟蹤環(huán)對實測數(shù)據(jù)中的衛(wèi)星信號進行跟蹤處理，跟蹤環(huán)路參數(shù)選擇為：載波環(huán)中，阻尼系數(shù)ξ=0.7，噪聲帶寬BL=25Hz；碼環(huán)中，阻尼系數(shù)ξ=0.7，噪聲帶寬BL=2Hz。MATLAB軟件仿真結(jié)果如圖6所示(取PRN21衛(wèi)星為例)。

圖6反映了PLL和DLL的跟蹤質(zhì)量。圖中左半部分，PLL鑒別器的輸出經(jīng)過短暫的牽引過程后，載波相位誤差抖動逐漸變小，這表明載波環(huán)路維持在良好的跟蹤狀態(tài)。PLL濾波器的輸出隨時間的推移，載波多普勒頻率變化近似呈線性，最終能得到1個穩(wěn)定的信號用來調(diào)節(jié)NCO載波發(fā)生器，生成與輸入信號載波相一致的本地載波。圖中右半部分，DLL鑒別器輸出的碼相位誤差抖動在零值附近相對穩(wěn)定，這表明碼跟蹤環(huán)路維持在良好的跟蹤狀態(tài)。DLL濾波器的輸出隨時間的推移，碼速率變化近似呈線性，最終得到1個穩(wěn)定的信號用來調(diào)節(jié)PRN碼發(fā)生器，生成與輸入信號偽碼完全對齊的本地偽碼。圖6結(jié)果表明PLL與DLL優(yōu)化組合的跟蹤環(huán)能正確、穩(wěn)定的跟蹤輸入信號的載波頻率和C/A碼相位，完成對輸入信號的解調(diào)和解擴。

圖6 PRN21衛(wèi)星載波環(huán)和碼環(huán)鑒別器及環(huán)路濾波器輸出

圖7為PLL處于鎖定狀態(tài)時碼跟蹤環(huán)中6個相關(guān)器的輸出結(jié)果。當(dāng)信號被正確跟蹤時，可以看出信號能量主要集中在I路，從而驗證了本文所設(shè)計的碼跟蹤環(huán)不受載波相位影響的優(yōu)點。

圖7 碼跟蹤環(huán)中同相支路和正交支路的6個相關(guān)器輸出

圖8所示分別為PRN21衛(wèi)星經(jīng)過1s跟蹤積分累加后輸出的Q路和I路分量。Q路輸出為噪聲分量，I路輸出為原始的導(dǎo)航電文數(shù)據(jù)。在整個GPS信號跟蹤環(huán)節(jié)，它們反映出信號的跟蹤質(zhì)量并確定軟件接收機的整體性能。

接收機完成跟蹤后，需要把跟蹤環(huán)I路輸出的積分累加值轉(zhuǎn)換為導(dǎo)航數(shù)據(jù)，理論上1ms輸出1bit數(shù)據(jù)。但實際處理的是帶有噪聲的弱信號，導(dǎo)航電文每一比特持續(xù)20ms，所以應(yīng)每 20ms計算1個平均值，輸出1bit數(shù)據(jù)，并二值化為1和-1[6,10,11]。

圖8 PRN21衛(wèi)星跟蹤環(huán)中Q支路和I支路輸出

在導(dǎo)航數(shù)據(jù)轉(zhuǎn)換時，對數(shù)據(jù)進行簡單的歸一化處理，歸一化公式為：

(15)

根據(jù)導(dǎo)航數(shù)據(jù)位跳變來標(biāo)記每一導(dǎo)航比特的起始位置，相鄰位跳變時刻之差除以20并四舍五入，得到跳變位之間的比特數(shù)。圖9為跟蹤環(huán)輸出對應(yīng)PRN21衛(wèi)星的以比特為單位的部分導(dǎo)航數(shù)據(jù)轉(zhuǎn)換結(jié)果。

圖9 PRN21衛(wèi)星歸一化導(dǎo)航數(shù)據(jù)

由圖9可以看出，利用PLL與DLL優(yōu)化組合的跟蹤環(huán)形式能準(zhǔn)確且穩(wěn)定的跟蹤GPS衛(wèi)星信號，最終得到精確的碼相位值和載波頻率值，正確解調(diào)出各個衛(wèi)星對應(yīng)的導(dǎo)航數(shù)據(jù)，充分驗證了此跟蹤算法的優(yōu)越性及其在GPS軟件接收機系統(tǒng)中應(yīng)用的可行性。

3 結(jié)論

詳細(xì)闡述了GPS軟件接收機對數(shù)字中頻信號的捕獲、跟蹤及解調(diào)過程。主要分析了基于L1波段C/A碼信號的單頻GPS軟件接收機的捕獲和跟蹤算法方案，即在捕獲階段采用基于頻域的并行碼相位搜索捕獲方法，在跟蹤階段采取PLL與DLL優(yōu)化組合的跟蹤環(huán)形式。最后利用MATLAB軟件實現(xiàn)了GPS信號捕獲和跟蹤算法的仿真，解調(diào)出導(dǎo)航數(shù)據(jù)，驗證了信號處理算法在系統(tǒng)中應(yīng)用的正確性和可行性。

[1]OlivierJulien,BoZheng,LeiDong,GérardLachapelle.ACompleteSoftware-BasedifGNSSSignalGeneratorforSoftwareReceiverDevelopment[C].ProceedingofIONGNSS2004,LongBeach,CA, 2004.

[2]MaC,LachapelleG,CannonME.ImplementationofSoftwareGPSReceiver[C].ProceedingofIONGNSS2004 (SessionA3),LongBeach,CA, 2004.

[3]FalinWu,KefeiZhang,andAkioYasuda.GPSSignalAcquisitionandTrackingUsingSoftwareGPSReceiver

[C]. 8thInternationalSymposiumonSignalProcessingandItsApplications,Sydney, 2005.

[4]AkopianD.FastFFTBasedGPSSatelliteAcquisitionMethods[J].Radar,SonarandNavigation, 2005, 152 (4):277-286.

[5]DineshManandhar,RyosukeShibasaki,HideyukiTorimoto.GPSReflectedSignalAnalysisUsingSoftwareReceiver[J].JournalofGlobalPositioningSystems, 2006, 5 (1-2):29-34.

[6]KaiBorre,DennisMAkos,etal.楊東凱，張飛舟，張波，譯.軟件定義的GPS和伽利略接收機[M].北京：國防工業(yè)出版社，2009：63-102. (KaiBorre,DennisMAkos,etal.YangDongkai,ZhangFeizhou,ZhangBo,etal.Translated.ASoftware-DefinedGPSandGalileoReceiver[M].Beijing:NationalDefenseIndustryPress, 2009:63-102.)

[7] 謝剛.GPS原理與接收機設(shè)計[M].北京：電子工業(yè)出版社，2011：372-374. (XieGang.PrinciplesofGPSandReceiverDesign[M].Beijing:PublishingHouseofElectronicsIndustry, 2011:372-374.)

[8] 易維勇，董緒榮，孟凡玉，郭晶，等.GNSS單頻軟件接收機應(yīng)用與編程[M].北京：國防工業(yè)出版社，2010：190-195.(YiWeiyong,DongXurong,MengFanyu,GuoJing,etal.GNSSSingleFrequenceSoftwareReceiverDevelopmentandApplications[M].Beijing:NationalDefenseIndustryPress, 2010:190-195.)

[9]RonaldE.BestPhase-LockedLoops:Design,Simulation,andApplications[M]. 6ed.NewYork:McGraw-HillProfessional, 2004:364-366.

[10]JamesBao-YenTsui.FundamentalsofGlobalPositioningSystemReceiversASoftwareApproach[M].NewYork:JohnWiley&Sons,Inc, 2000:133-144.

[11] 田炳麗，寧春林，丁風(fēng)雷，李建輝.GPS軟件接收機導(dǎo)航電文提取的研究[J].通信技術(shù)，2009，42(3)：212-217.(TianBingli,NingChunlin,DingFenglei,LiJianhui.NavigationDataExtractionBasedonGPSSoftwareReceiver[J].CommunicationsTechnology, 2009, 42 (3):212-217.)

Algorithm and Verification of Baseband Signal Processing for GPS Software Receiver

Liu Ruihua1, Kong Yueming1, Yu Baoguo2, Xing Zhaodong2

1.School of Electronic and Information Engineering, CAUC, Tianjin 300300, China 2.The 54th Research Institute of CETC, Shijiazhuang 050081, China

Inthispaper,thebasebandsignalprocessingflowandalgorithmforGPSsoftwarereceiverbasedonadatablockarestudied,includingtheacquiringandtrackingofcarrierandpseudo-codeaswellasdemodulatingthenavigationdata.TheMATLABsoftwareplatformisusedtosimulateandvalidateeachlinkofthebasebandsignalprocessingforGPSsoftwarereceiverbasedontheactual

satellitesignal.Duringtheexperiment,thesatellitesignalcanbepreferablyacquiredbyusingthisalgorithmresearched,thecarrierandpseudo-codecanbeaccuratelytrackedandthenavigationdatacanbecorrectlydemodulated.Furthermore,thefeasibilityofitsapplicationinthesoftwarereceiversystemisverified,italsocanserveasthereferenceforthedevelopmentofreal-timeGPSsoftwarereceiver.

GPSsoftwarereceiver;Acquisition;Tracking;Simulationandverification;Feasibility

*民航安全能力建設(shè)項目“北斗機載設(shè)備技術(shù)標(biāo)準(zhǔn)規(guī)定與應(yīng)用研究”(AADSA0007)；中央高校基本科研業(yè)務(wù)費中國民航大學(xué)專項“基于北斗的通用航空指揮監(jiān)控系統(tǒng)”(20001006)；中電科54所新技術(shù)研究高校合作項目“通航機載GNSS終端完好性技術(shù)研究”(KX14260033)

2013-03-12

劉瑞華(1965-)，男，陜西人，博士，教授， 研究生導(dǎo)師，主要研究方向為衛(wèi)星導(dǎo)航、 慣性導(dǎo)航和組合導(dǎo)航；孔月明(1989-)，女，山東人，碩士，主要研究方向為衛(wèi)星導(dǎo)航新技術(shù)及應(yīng)用；蔚保國(1966-)，男，內(nèi)蒙古人，博士，研究員，博士生導(dǎo)師，主要研究方向為衛(wèi)星導(dǎo)航總體技術(shù)、航天測控技術(shù)、陣列信號處理技術(shù)和自動測試系統(tǒng)技術(shù)；邢兆棟(1976-)，男，山東人，博士，高級工程師，主要研究方向為衛(wèi)星導(dǎo)航信號處理及航空應(yīng)用。

TN967.1

A

1006-3242(2016)01-0003-06