楊開偉,李娟娟,杜聰慧

(1.中國電子科技集團(tuán)公司第54研究所,石家莊 050081;2.河北省衛(wèi)星導(dǎo)航技術(shù)與裝備工程技術(shù)研究中心,石家莊 050081;3衛(wèi)星導(dǎo)航系統(tǒng)與裝備技術(shù)國家重點(diǎn)實(shí)驗(yàn)室,石家莊 050081)

?

嵌入式GPS/BDS雙系統(tǒng)RTK算法實(shí)現(xiàn)及性能測試

楊開偉1,2,3,李娟娟1,2,3,杜聰慧1,2,3

(1.中國電子科技集團(tuán)公司第54研究所,石家莊 050081;2.河北省衛(wèi)星導(dǎo)航技術(shù)與裝備工程技術(shù)研究中心,石家莊 050081;3衛(wèi)星導(dǎo)航系統(tǒng)與裝備技術(shù)國家重點(diǎn)實(shí)驗(yàn)室,石家莊 050081)

本文基于支持BDS系統(tǒng)和GPS系統(tǒng)的高精度OEM板卡進(jìn)行嵌入式RTK算法設(shè)計(jì),實(shí)現(xiàn)BDS單系統(tǒng)、GPS單系統(tǒng)RTK定位功能,并進(jìn)行室外實(shí)地測試。測試結(jié)果表明，開闊空間靜態(tài)零基線測量北斗系統(tǒng)水平精度亞毫米,垂直1 mm;GPS系統(tǒng)三維統(tǒng)計(jì)精度均優(yōu)于3 mm．開闊空間靜態(tài)10 km基線測量GPS解算結(jié)果較好,水平誤差小于1 cm,垂直誤差大部分小于3 cm,統(tǒng)計(jì)誤差水平4 mm以內(nèi),垂直2 cm以內(nèi);北斗雙頻解算三維誤差大部分在4 cm以內(nèi),統(tǒng)計(jì)精度水平優(yōu)于1.5 cm,垂直優(yōu)于3.5 cm.BDS系統(tǒng)動態(tài)RTK測量同實(shí)際路線符合良好,單歷元中誤差在6 mm以內(nèi)。

RTK;BDS;GPS;高精度定位

0　引　言

隨著我國北斗衛(wèi)星區(qū)域?qū)Ш较到y(tǒng)的正式運(yùn)行,利用北斗衛(wèi)星系統(tǒng)進(jìn)行高精度測量成為導(dǎo)航領(lǐng)域的研究熱點(diǎn)。利用多模多頻高精度OEM板能夠自主、方便地開發(fā)出各種滿足特定需求的北斗導(dǎo)航應(yīng)用系統(tǒng)。將OEM板與計(jì)算機(jī)、通訊技術(shù)相結(jié)合,開發(fā)能夠獲取厘米級實(shí)時動態(tài)定位精度的嵌入式RTK算法,即可廣泛應(yīng)用于傳統(tǒng)測繪、駕校系統(tǒng)、交通運(yùn)輸、變形監(jiān)測、城市管理等領(lǐng)域。

多個頻率多個衛(wèi)星導(dǎo)航系統(tǒng)組合導(dǎo)航定位,能夠在作業(yè)環(huán)境較差的區(qū)域,提供穩(wěn)定、可靠的定位結(jié)果,擴(kuò)展作業(yè)范圍。同時可以提高RTK的作業(yè)距離,提高系統(tǒng)的穩(wěn)定性。鑒于高精度衛(wèi)星導(dǎo)航定位的應(yīng)用需求,研制北斗和GPS的多模多頻高精度OEM板卡,嵌入RTK高精度定位等關(guān)鍵技術(shù),實(shí)現(xiàn)商品化生產(chǎn)受到越來越多的關(guān)注。

1　嵌入式實(shí)時RTK策略設(shè)計(jì)

RTK是一種利用載波相位觀測值進(jìn)行實(shí)時動

態(tài)相對定位的技術(shù)。本文采用基于FPGA原型的多模多頻高精度OEM板,整個系統(tǒng)核心器件是DSP和FPGA.RTK數(shù)據(jù)處理部分是嵌入到DSP上,實(shí)現(xiàn)實(shí)時高精度定位。

嵌入式實(shí)時RTK策略設(shè)計(jì)主要包含RTK任務(wù)與單機(jī)協(xié)同,實(shí)時周跳探測、增減衛(wèi)星條件下的模糊度搜索、以及參考衛(wèi)星變換時的模糊度固定等內(nèi)容。其目的是在復(fù)雜環(huán)境或長時運(yùn)行條件下,在不增加運(yùn)算復(fù)雜度和空間復(fù)雜度的前提下,能夠保證高精度數(shù)據(jù)結(jié)果的穩(wěn)定輸出。策略的成功與否,直接關(guān)系到OEM板卡的穩(wěn)定性以及數(shù)據(jù)結(jié)果的質(zhì)量[1-2]。

1.1RTK定位模型

載波相位動態(tài)相對定位采用消除了衛(wèi)星鐘差、接收機(jī)鐘差,減弱了衛(wèi)星軌道誤差對基線解算的影響的雙差載波相位觀測值。雙差觀測方程為[3-4]:

Δφ(t)· λ=Δρ-ΔN·λ+Δdorb+

ΔIon+ΔTrop+∑δ,

(1)

式中:Δφ(t)為載波相位觀測值; λ為載波波長;Δdorb為軌道誤差;ΔIon為電離層延遲;ΔTrop對流層延遲;Δρ為幾何距離;ΔN為整周模糊度。

設(shè)兩測站同步觀測衛(wèi)星為Sj和Sk,設(shè)Sj為參考衛(wèi)星,則可得雙差觀測方程的線性化形式,即:

Δφk·λ=-[

ΔNk·λ+Δρ,

(2)

式中,下標(biāo)2代表流動站,

Δφk=φk-,

式(2)可改寫為

ΔNk+Δuk(t),

(3)

式中,Δuk(t)=Δρ.

由于流動站的位置是不斷變化的,所以RTK只能采用一個歷元的觀測值,必須預(yù)先解算出整周模糊度,這就是RTK初始化過程。只要同步觀測的衛(wèi)星數(shù)大于4顆,就可以進(jìn)行RTK定位解算。RTK初始化以后,采用最小二乘方法求解。

1.2流程圖

RTK解算主要流程如圖1所示。

圖1　RTK解算流程圖

2　模糊度解算

RTK的關(guān)鍵技術(shù)是在動態(tài)觀測過程中,要求保持對所測衛(wèi)星的連續(xù)跟蹤。一旦發(fā)生失鎖,則要重新進(jìn)行初始化工作。近年來的許多學(xué)者都致力于這一方面的研究和開發(fā)工作。根據(jù)本文RTK數(shù)學(xué)模型,并且采用雙頻觀測值進(jìn)行RTK解算,故需要固定兩個頻點(diǎn)的雙差整周模糊度。本文采用LAMBDA算法進(jìn)行模糊度搜索,統(tǒng)計(jì)檢驗(yàn)采用ratio值檢驗(yàn)法,令ratio值大于3.5為模糊度固定正確,給出固定解[2]。

3　周跳探測

準(zhǔn)確探測與修復(fù)周跳是高精度定位的關(guān)鍵環(huán)節(jié),可以提高定位的連續(xù)性及可靠性。本文使用的DSP,RTK運(yùn)算過程中最多分配12個數(shù)組空間,故需要對PC環(huán)境下編寫的RTK程序進(jìn)行數(shù)組壓縮?？紤]到工程應(yīng)用于實(shí)時性的要求,基于DSP的周跳探測采用電離層殘差法。此方法是利用電離層殘差檢驗(yàn)量歷元間的變化來檢測是否出現(xiàn)周跳。

所用數(shù)據(jù)模型為[5]

L4=λ1φ1-λ2φ2,

(4)

式中,L4組合觀測值消除了站星幾何距離、接收機(jī)鐘差和衛(wèi)星鐘差的影響,消弱了電離層延遲、多路徑效應(yīng)的影響。

第i與第i+1歷元間的L4變化量:

ΔL4(ti,ti+1)=(λ1·φ1(ti+1)-λ2·φ2(ti+1))-(λ1·φ1(ti)-λ2·φ2(ti)).

(5)

當(dāng)電離層比較穩(wěn)定、采樣間隔較短時,ΔL4(ti,ti+1)幾乎為零。當(dāng)載波相位發(fā)生周跳時,電離層殘差不再連續(xù),ΔL4(ti,ti+1)會發(fā)生跳變,很容易探測周跳。

本文采用探測不修復(fù)周跳的方法,如果非參考星周跳,當(dāng)前歷元對周跳的衛(wèi)星重新進(jìn)行模糊度固定;如果參考星發(fā)生周跳,因?yàn)槟：仁撬行l(wèi)星相對于參考星的雙差模糊度,故當(dāng)前歷元需要對所有衛(wèi)星進(jìn)行模糊度搜索。

4　性能測試

為了評定OEM板卡RTK定位精度,本文設(shè)計(jì)了三項(xiàng)測試項(xiàng)目,分別為零基線測試、10 km靜態(tài)基線測試和動態(tài)RTK測量。

4.1靜態(tài)零基線測試

本次試驗(yàn)接收機(jī)天線架在某樓頂基準(zhǔn)點(diǎn),觀測環(huán)境相對比較開闊,天線經(jīng)功分器形成兩路信號,一路送給基準(zhǔn)站接收機(jī),一路送給流動站接收機(jī),形成RTK工作模式。測試采樣率為1 Hz,數(shù)據(jù)采集1 h,設(shè)置衛(wèi)星截止高度角為10°,接收機(jī)實(shí)時處理,流動站按1 Hz實(shí)時上報GGA語句,通過串口連接到PC機(jī)上存儲。零基線靜態(tài)數(shù)據(jù)分別采用單GPS、單北斗進(jìn)行RTK解算。測量結(jié)果與真值“0”之差如圖2～圖3所示,差值統(tǒng)計(jì)(1σ)如表1所示。由圖和表可以看出，北斗雙頻解算結(jié)果N方向、E方向誤差均小于4 mm，U方向誤差大部分小于5 mm，統(tǒng)計(jì)誤差(1 σ )N方向、E方向?yàn)閬喓撩?，U方向?yàn)? mm；GPS解算三維誤差大部分在8 mm以內(nèi)，統(tǒng)計(jì)精度(1 σ )優(yōu)于3 mm.即開闊空間靜態(tài)零基線RTK測量精度可達(dá)到毫米級。

圖2　靜態(tài)零基線GPS解算結(jié)果

圖3　靜態(tài)零基線BDS解算結(jié)果

4.2靜態(tài)10 km基線測試

本次試驗(yàn)基準(zhǔn)站架在某樓頂基準(zhǔn)點(diǎn),流動站約10 km遠(yuǎn)的另一基準(zhǔn)點(diǎn),通過3G模塊進(jìn)行數(shù)據(jù)傳輸,形成RTK工作模式。流動站基準(zhǔn)點(diǎn)相對于基準(zhǔn)站基準(zhǔn)點(diǎn)N/E/U方向真值為8 029.7 816/5 773.905 9/28.400 4.測試采樣率為0.1 Hz,數(shù)據(jù)采集2 h,設(shè)置衛(wèi)星截止高度角為10°.計(jì)算結(jié)果如圖4～圖5所示,與真值差值統(tǒng)計(jì)(1σ)如表2所示。

表1　靜態(tài)零基線解算精度統(tǒng)計(jì) 　(單位：mm)

圖4　靜態(tài)10 km基線GPS解算結(jié)果

圖5　靜態(tài)10 km基線BDS解算結(jié)果

由圖和表可以看出，GPS解算結(jié)果N方向、E方向誤差均小于1 cm，U方向誤差大部分小于3 cm，統(tǒng)計(jì)誤差(1σ)N方向、E方向?yàn)? mm以內(nèi)，U方向2 cm以內(nèi)；北斗解算三維誤差大部分在4 cm以內(nèi)，統(tǒng)計(jì)精度(1σ)N方向、E方向優(yōu)于1.5 cm，U方向優(yōu)于3.5 cm.

表2靜態(tài)10 km基線解算精度統(tǒng)計(jì) (單位：mm)

處理模式固定解比例N方向E方向U方向GPS95.8%3.53.319.6BDS97.1%11.512.731.6

4.3動態(tài)RTK測試

RTK動態(tài)測試試驗(yàn)在市區(qū)公路上進(jìn)行,基準(zhǔn)站架在樓頂基準(zhǔn)點(diǎn)上,流動站在開闊路段上運(yùn)行,通過電臺進(jìn)行數(shù)據(jù)傳輸,形成RTK工作模式,測試采樣率為1 Hz,數(shù)據(jù)采集0.5 h,設(shè)置衛(wèi)星截止高度角為10°,接收機(jī)采用BDS系統(tǒng)雙頻數(shù)據(jù)實(shí)時處理,流動站按1 Hz實(shí)時上報GGA語句,通過串口連接到PC機(jī)上存儲。流動站與基準(zhǔn)站之間的距離為1～2 km,試驗(yàn)持續(xù)時間為1 h.試驗(yàn)數(shù)據(jù)的處理結(jié)果與實(shí)際地圖匹配后如圖6所示。推車先圍繞轉(zhuǎn)盤轉(zhuǎn)行兩圈,然后在東西路上來回,回程過程中沒有數(shù)據(jù)點(diǎn)的兩處是重新開機(jī)引起的,再圍繞轉(zhuǎn)盤轉(zhuǎn)行一圈,然后回到工作樓。單歷元解算中誤差指標(biāo)如圖7所示。由圖可以看出,解算的中誤差基本在6 mm以內(nèi)。

圖6　BDS動態(tài)RTK測試結(jié)果

圖7　BDS動態(tài)RTK測試中誤差

5　結(jié)束語

本文基于多模多頻高精度OEM板卡,討論了嵌入式RTK算法策略,進(jìn)行了BDS單系統(tǒng)、GPS單系統(tǒng)RTK定位功能實(shí)現(xiàn)和室外實(shí)地測試。測試結(jié)果表明:

1) 開闊空間靜態(tài)零基線RTK測量精度可達(dá)到毫米級;

2) 開闊空間靜態(tài)10 km基線測量GPS統(tǒng)計(jì)誤差水平4 mm以內(nèi),垂直2 cm以內(nèi);北斗統(tǒng)計(jì)精度水平優(yōu)于1.5 cm,垂直優(yōu)于3.5 cm;

3) BDS動態(tài)RTK測量同實(shí)際路線符合良好,單歷元平差中誤差在6 mm以內(nèi)。

本文基本實(shí)現(xiàn)了嵌入式RTK基線解算,初步測試性能較好。下一步的工作是:靜態(tài)基線測量精度雖然較好,但需要進(jìn)一步測試穩(wěn)定性及長時間運(yùn)行可行性;動態(tài)測量垂直方向精度不夠高,需要對周跳探測和模糊度解算等方面進(jìn)一步分析研究與完善,以提高解算精度。

[1]KAPLAN E D, HEGARTY C J.GPS原理與應(yīng)用[M].寇艷紅譯.2版.北京:電子工業(yè)出版社,2010.

[2]謝鋼.GPS原理與接收機(jī)設(shè)計(jì)[M].北京:電子工業(yè)出版社,2009.

[3]李征航,黃勁松.GPS測量與數(shù)據(jù)處理[M].武漢:武漢大學(xué)出版社,2011:162-168.

[4]任哲.嵌入式操作系統(tǒng)基礎(chǔ)[M].北京:北京航空航天大學(xué)出版社,2011.

[5]邵儀,基于FPGA的矩陣運(yùn)算固化實(shí)現(xiàn)技術(shù)研究 [D].鄭州:解放軍信息工程大學(xué),2012.

The Realization of the Embedded GPS/BDS RTK Algorithm and the Test of its Performance

YANG Kaiwei1,2,3,LI Juanjuan1,2,3,DU Conghui1,2,3

(1.The54thResearchInstituteofCETC,Shijiazhuang050081,China；2.SatelliteNavigationTechnologyandEquipmentEngineeringTechnologyResearchCenterofHebeiprovince,Shijiazhuang050081,China；3.StateKeyLaboratoryofSatelliteNavigationSystemandEquipmentTechnology,Shijiazhuang050081,China)

Using the high-accuracy OEM board, which contains of BDS and GPS system, this paper realized the embedded algorithm of RTK of BDS、GPS and BDS+GPS. The test of outdoor proved that the accuracy of zero static baseline is sub-mm for the horizontal component and 1 mm for vertical component using BDS. The accuracy of GPS for 10 km static baseline, whose difference of horizontal component is smaller than 1 cm and the vertical component is smaller than 3 cm, is better than BDS, whose difference for three component is better than 4 cm. The statistics accuracy of GPS for 10 km baseline is better than 4 mm for horizontal component and better than 2 cm for vertical component, and the BDS is better than 1.5 cm for horizontal component and better than 3.5 cm for vertical component. For kinematic baseline, the result of integration positioning is accordant to the real topography, and the standard deviation is better than 6 mm.

RTK; BDS; GPS; integration positioning; high-accuracy board

2015-10-09

P228.4

A

1008-9268(2016)03-0035-05

楊開偉(1984-),男,工程師,主要從事衛(wèi)星導(dǎo)航數(shù)據(jù)處理方面的工作。

李娟娟(1984-),女,工程師,主要從事衛(wèi)星導(dǎo)航數(shù)據(jù)處理方面的工作。

杜聰慧(1988-),女,工程師,主要從事衛(wèi)星導(dǎo)航數(shù)據(jù)處理方面的工作。

doi:10.13442/j.gnss.1008-9268.2016.03.008

聯(lián)系人： 楊開偉E-mail:yangkaiwei0603@163.com