劉偉平，郝金明，于合理，田英國

信息工程大學(xué)導(dǎo)航與空天目標(biāo)工程學(xué)院，河南 鄭州450001

?

北斗衛(wèi)星多天軌道合成方法及其精度分析

劉偉平，郝金明，于合理，田英國

信息工程大學(xué)導(dǎo)航與空天目標(biāo)工程學(xué)院，河南 鄭州450001

目前，北斗衛(wèi)星多天解軌道主要是通過常規(guī)的數(shù)據(jù)累積解算獲得，處理效率較低。此外，由于開通運(yùn)行不久，北斗衛(wèi)星力學(xué)模型尚未完善，限制了多天解軌道精度的進(jìn)一步提高。針對以上問題，本文提出了附加偽隨機(jī)脈沖的北斗衛(wèi)星多天軌道合成方法，利用法方程疊加實(shí)現(xiàn)北斗衛(wèi)星多天軌道合成，可有效提高軌道解算效率。同時(shí)，利用偽隨機(jī)脈沖對力學(xué)模型的補(bǔ)償作用，進(jìn)一步提高了多天軌道合成的精度水平。實(shí)測數(shù)據(jù)分析表明：在多天軌道合成中，單天邊界附加偽隨機(jī)脈沖對北斗衛(wèi)星軌道精度有明顯改進(jìn)作用。其中，3類衛(wèi)星切向、法向軌道精度的改善程度大于徑向，GEO和IGSO衛(wèi)星軌道精度的提高幅度大于MEO衛(wèi)星。一定范圍內(nèi)，隨著定軌弧段的增長，北斗衛(wèi)星定軌精度能夠得到不同程度的改善，且主要表現(xiàn)在軌道切向。

北斗衛(wèi)星導(dǎo)航系統(tǒng)；精密定軌;法方程疊加;軌道合成;偽隨機(jī)脈沖

北斗衛(wèi)星導(dǎo)航系統(tǒng)是我國自主建設(shè)、獨(dú)立運(yùn)行的全球衛(wèi)星導(dǎo)航系統(tǒng)[1-2]，其軌道精度水平直接影響到系統(tǒng)的拓展應(yīng)用[3-4]。在導(dǎo)航衛(wèi)星精密定軌中，通常以一天為一個(gè)時(shí)段進(jìn)行處理，然而由于單天解軌道存在末端效應(yīng)，即時(shí)段兩端的軌道誤差相對較大，因此，需要將單天解軌道合成為多天解軌道[5-6]。

以往的處理是將幾個(gè)單天時(shí)段的數(shù)據(jù)集合到一起，而后統(tǒng)一處理得到多天解軌道，計(jì)算耗時(shí)耗力，并隨著跟蹤站數(shù)目及處理天數(shù)的增多，對計(jì)算機(jī)性能的要求越來越苛刻。而后發(fā)展起來的法方程疊加方法通過并行處理單天觀測數(shù)據(jù)，再在法方程層面進(jìn)行軌道合成，極大地提高了處理效率[7-8]。國外關(guān)于法方程疊加方法的研究起步較早，早在20世紀(jì)90年代中后期，文獻(xiàn)[5]就研究了應(yīng)用法方程疊加合成多天軌道的方法；隨后，文獻(xiàn)[7]系統(tǒng)探討了法方程疊加原理及其在大地測量與地球動力學(xué)中的應(yīng)用；文獻(xiàn)[9]將法方程疊加原理應(yīng)用到了區(qū)域形變估計(jì)中。國內(nèi)，文獻(xiàn)[10]在深入研究了法方程疊加原理的基礎(chǔ)上，探討了基于法方程疊加的多時(shí)段軌道合成方法；文獻(xiàn)[11]利用法方程疊加原理，分析了多個(gè)子網(wǎng)解軌道合成整網(wǎng)解軌道的處理效果。然而，以上研究大多是針對GPS和GLONASS衛(wèi)星，目前北斗衛(wèi)星采用的多天解軌道主要是依靠傳統(tǒng)的多天數(shù)據(jù)累積解算獲得[12-17]。這種處理模式在精度上與基于法方程疊加的結(jié)果是一致的，但是計(jì)算過程復(fù)雜，無法并行處理各天數(shù)據(jù)，計(jì)算效率較低。此外，在軌道合成中，涉及的主要軌道參數(shù)包括密切元素、力學(xué)參數(shù)和偽隨機(jī)脈沖參數(shù)，已有研究主要集中于密切元素(參考時(shí)刻的6個(gè)開普勒根數(shù))和力學(xué)參數(shù)(如光壓模型參數(shù)等)的合成處理，而對另一類重要參數(shù)——偽隨機(jī)脈沖的處理方法雖有所涉及，但相關(guān)研究還不夠深入。偽隨機(jī)脈沖具有彌補(bǔ)力學(xué)模型精度不足的作用[18]，對當(dāng)前數(shù)據(jù)累積有限、部分?jǐn)z動力難以精確建模的北斗衛(wèi)星精密定軌具有重要意義。

針對以上問題，本文首先給出了法方程疊加的基本原理，并重點(diǎn)推導(dǎo)了軌道合成中主要應(yīng)用的參數(shù)轉(zhuǎn)換方法。而后，在給出密切元素和力學(xué)參數(shù)合成方法的基礎(chǔ)上，著重研究了偽隨機(jī)脈沖的處理方法。據(jù)此，提出了附加偽隨機(jī)脈沖的北斗衛(wèi)星多天軌道合成方法。最后，利用北斗衛(wèi)星實(shí)測數(shù)據(jù)，設(shè)計(jì)相應(yīng)算例，驗(yàn)證了本文方法的應(yīng)用效果。

1 法方程疊加基本原理

在基于法方程疊加的軌道合成中，首先處理單天數(shù)據(jù)，形成單天解法方程并保存，再在法方程水平進(jìn)行參數(shù)轉(zhuǎn)換、疊加等處理，得到多天解軌道。參數(shù)轉(zhuǎn)換的本質(zhì)是對法方程進(jìn)行重新參數(shù)化，其基本思想是通過建立新參數(shù)與原參數(shù)的線性關(guān)系，并利用這一關(guān)系對原法方程系統(tǒng)進(jìn)行轉(zhuǎn)換，形成新的法方程系統(tǒng)，在進(jìn)行法方程疊加之前，通常需要進(jìn)行參數(shù)轉(zhuǎn)換。這里首先推導(dǎo)了法方程轉(zhuǎn)換方法，而后給出了法方程疊加的基本原理。

1.1 參數(shù)轉(zhuǎn)換

設(shè)有誤差方程

(1)

其對應(yīng)的法方程為

(2)

式中，N=BTPB；U=BTPL；P為觀測值的權(quán)陣。

(3)

將式(3)代入式(1)，得

(4)

則新的法方程系統(tǒng)可表示為

(5)

式中

N′=CTBTPBC=CTNC

(6)

U′=CTU-CTNd

(7)

式(6)和式(7)即為對應(yīng)的參數(shù)轉(zhuǎn)換公式。實(shí)際應(yīng)用中，首先建立新參數(shù)與原參數(shù)的線性關(guān)系(如式(3))，即給定C和d的具體形式，而后直接由參數(shù)轉(zhuǎn)換公式獲得新法方程的系數(shù)矩陣，即可建立新的法方程系統(tǒng)。

1.2 法方程疊加

假設(shè)有兩組觀測量，相互之間獨(dú)立，且有觀測方程如下

(8)

(9)

按照最小二乘原理，式(8)和式(9)對應(yīng)的法方程可分別表示為

(10)

(11)

兩組法方程疊加得到如下新的法方程系統(tǒng)

(12)

2 附加偽隨機(jī)脈沖的多天軌道合成

通過法方程疊加進(jìn)行多天軌道合成時(shí)，可將對應(yīng)的單天解法方程中包含的參數(shù)分為兩大類：一類為“感興趣”的軌道參數(shù)；另一類為暫時(shí)“不感興趣”的參數(shù)，如對流層延遲參數(shù)、鐘差參數(shù)等。對于后者，可以在疊加之前首先將其預(yù)消除，待到疊加完成并解得軌道參數(shù)后，再通過回代的方法進(jìn)行解算。關(guān)于參數(shù)的預(yù)消除與回代方法，相關(guān)文獻(xiàn)多有論及[7,10-11]，不再贅述。這里主要關(guān)注軌道參數(shù)，通常包括密切元素和力學(xué)參數(shù)。結(jié)合北斗衛(wèi)星目前數(shù)據(jù)累積有限、部分力學(xué)模型有待進(jìn)一步精化的特點(diǎn)，這里還考慮在軌道合成時(shí)于單天邊界處附加偽隨機(jī)脈沖參數(shù)，以吸收力學(xué)模型誤差。下面首先給出密切元素和力學(xué)參數(shù)的合成方法，而后重點(diǎn)推導(dǎo)偽隨機(jī)脈沖的具體處理方法。

2.1 密切元素和力學(xué)參數(shù)的合成

假定有n天軌道需要合成，那么只要能夠建立第1天的密切元素和力學(xué)參數(shù)與其后n-1天相關(guān)參數(shù)的線性關(guān)系，并利用這一關(guān)系對n-1個(gè)單天解法方程進(jìn)行參數(shù)變換，以統(tǒng)一所有單天解法方程的待估參數(shù)集，最后，根據(jù)法方程疊加原理，即可完成相關(guān)參數(shù)的合成。

如果要將n-1天的軌道通過第1天的軌道參數(shù)表示，就要求相鄰兩天邊界處衛(wèi)星位置與速度具有連續(xù)性，同時(shí)可要求各天的力學(xué)參數(shù)相等。有如下關(guān)系

(13)

利用密切元素與衛(wèi)星位置、速度的對應(yīng)關(guān)系以及式(13)，經(jīng)過推導(dǎo)，可以得到如下線性關(guān)系

(14)

式中，ΔEi、Δqi分別表示第i天(i=2,3,…,n)的密切元素和力學(xué)參數(shù)；ΔE1、Δq1分別表示第1天的密切元素和力學(xué)參數(shù)；Ki,1、Li,1、Mi,1、Ni,1分別表示對應(yīng)的關(guān)系矩陣，具體形式可參閱文獻(xiàn)[7]；I表示單位陣。

利用式(14)建立的線性關(guān)系，參照1.1節(jié)給出的參數(shù)轉(zhuǎn)換方法，對n-1天的法方程系統(tǒng)進(jìn)行參數(shù)轉(zhuǎn)換，統(tǒng)一各天的待估參數(shù)集，而后利用1.2節(jié)給出的法方程疊加方法，即可完成相應(yīng)的合成處理。

2.2 偽隨機(jī)脈沖參數(shù)的合成

為了減弱北斗衛(wèi)星力學(xué)模型精度不足的不利影響，考慮在單天邊界處附加偽隨機(jī)脈沖參數(shù)，其可定義為在給定歷元的瞬間速度增量[19-20]，通過下式表示

vnew(τ)=vold(τ)+so

(15)

式中，s表示偽隨機(jī)脈沖的大小，o表示偽隨機(jī)脈沖的方向，兩者共同定義了偽隨機(jī)脈沖；τ表示偽隨機(jī)脈沖發(fā)生的時(shí)刻；vold、vnew分別表示偽隨機(jī)脈沖發(fā)生前后的衛(wèi)星速度。

在進(jìn)行軌道合成時(shí)，需要考慮當(dāng)前時(shí)刻之前所有偽隨機(jī)脈沖對衛(wèi)星軌道的影響。為便于討論，假定τ為附加偽隨機(jī)脈沖對應(yīng)的參考?xì)v元，ti為第i天法方程對應(yīng)的參考?xì)v元，分兩步處理以上問題：

2.2.1 τ歷元偽隨機(jī)脈沖對τ歷元密切元素的影響

將τ歷元偽隨機(jī)脈沖so在徑向、橫向和軌道面法向的分量分別記為sR、sT、sN。根據(jù)牛頓受攝運(yùn)動方程[21]，引起的密切元素變化量可表示為

(16)

將τ歷元偽隨機(jī)脈沖對τ歷元密切元素的影響以矩陣形式表示，有

ΔEs(τ)=sκs(τ)

(17)

式中，ΔEs(τ)表示偽隨機(jī)脈沖引起的密切元素的變化量；κs(τ)是由式(16)給定的系數(shù)矩陣；s含義與以上相同。

2.2.2 τ歷元密切元素對ti歷元密切元素的影響

根據(jù)衛(wèi)星變分方程，可得狀態(tài)轉(zhuǎn)移矩陣Ms(ti,τ)，其元素為

(18)

則τ歷元密切元素對ti歷元密切元素的影響，可表示為

ΔEs(ti)=Ms(ti,τ)ΔEs(τ)

(19)

綜合式(17)和式(19)，可得τ歷元偽隨機(jī)脈沖對ti歷元密切元素的影響為

ΔEs(ti)=Tss

(20)

并有

Ts=Ms(ti,τ)κs(τ)

(21)

綜上，當(dāng)前時(shí)刻之前所有偽隨機(jī)脈沖對ti歷元密切元素的影響可通過式(20)獲得，表示如下

(22)

對第i天的法方程系統(tǒng)進(jìn)行參數(shù)轉(zhuǎn)換，對照式(3)，有

(23)

di=0

(24)

利用參數(shù)轉(zhuǎn)換公式(6)和式(7)，可得

(25)

(26)

利用以上方法，對各單天解法方程進(jìn)行參數(shù)轉(zhuǎn)換，得到新的法方程系統(tǒng)，而后基于法方程疊加原理，即可完成偽隨機(jī)脈沖參數(shù)的合成。

方便推導(dǎo)起見，以上分別討論了密切元素、力學(xué)參數(shù)和偽隨機(jī)脈沖的合成方法。實(shí)際處理中，3類參數(shù)同時(shí)出現(xiàn)在法方程中，需首先按照處理偽隨機(jī)脈沖的方法對各單天解法方程進(jìn)行參數(shù)轉(zhuǎn)換，而后進(jìn)行密切元素和力學(xué)參數(shù)的參數(shù)轉(zhuǎn)換，最后再進(jìn)行法方程疊加，即可完成多天軌道的合成處理。

3 算例分析

為了驗(yàn)證本文方法的正確性及其在北斗衛(wèi)星精密定軌中的實(shí)際應(yīng)用效果，結(jié)合北斗衛(wèi)星實(shí)測數(shù)據(jù)，首先探討單天邊界附加偽隨機(jī)脈沖對北斗衛(wèi)星軌道精度的改善作用，而后，分析北斗衛(wèi)星多天軌道合成的實(shí)際效果。

3.1 偽隨機(jī)脈沖的作用

采用全球24個(gè)測站2013年6月3日到2013年6月24日共22 d的北斗觀測數(shù)據(jù)進(jìn)行試驗(yàn)分析，其中，12個(gè)測站來自IGS組織的多系統(tǒng)實(shí)驗(yàn)網(wǎng)(multi-GNSS experiment，M-GEX)，12個(gè)測站來自武漢大學(xué)建立的北斗衛(wèi)星觀測實(shí)驗(yàn)網(wǎng)(BeiDou experimental tracking stations,BETS)，測站分布見圖1。

圖1 測站分布圖Fig.1 Station layout

處理中，首先利用單天數(shù)據(jù)以非差動力法進(jìn)行北斗衛(wèi)星單系統(tǒng)精密定軌，并保存單天解法方程，而后，利用本文給出的方法進(jìn)行3d軌道合成。為了方便對比，軌道合成時(shí)采用如下兩種方案，兩種處理方案僅在偽隨機(jī)脈沖的處理上有所區(qū)別，其他處理策略(見表1)相同。

方案1(S1)：3d軌道合成，單天邊界處不附加偽隨機(jī)脈沖。

方案2(S2)：3d軌道合成，單天邊界處附加偽隨機(jī)脈沖。

以3d軌道的重疊弧段軌道差值考察兩種方案的定軌精度，方法如圖2。表2和圖3給出了兩種方案下各衛(wèi)星在所有重疊弧段上徑向(R)、切向(T)、法向(N)和三維(P)的平均RMS情況。

表1 數(shù)據(jù)處理策略

圖2 重疊弧段軌道差值比較方法Fig.2 Method of overlapped-arc comparison

圖3 北斗衛(wèi)星精密定軌精度對比圖Fig.3 Comparison of BeiDou satellites orbit precision

Tab.2 Statistics of BeiDou satellites orbit precision m

由圖3和表2的結(jié)果可見，多天軌道合成時(shí)在單天邊界處附加偽隨機(jī)脈沖，所有北斗衛(wèi)星的精密定軌精度都獲得了明顯改善。平均來講，GEO衛(wèi)星軌道精度在R、T、N方向分別提高了0.062、0.191和0.179 m；IGSO衛(wèi)星分別提高了0.024、0.170和0.109 m；MEO衛(wèi)星分別提高了0.015、0.056和0.032 m。相比R方向，偽隨機(jī)脈沖的引入對3類衛(wèi)星T、N方向軌道精度的改善更為明顯。從三維定軌精度來看，GEO、IGSO和MEO衛(wèi)星分別提高了0.256、0.212和0.071 m，GEO和IGSO衛(wèi)星軌道精度的提高幅度明顯大于MEO衛(wèi)星。偽隨機(jī)脈沖的主要作用是彌補(bǔ)力學(xué)模型精度的不足，以上結(jié)果在一定程度上表明，北斗系統(tǒng)3類衛(wèi)星的力學(xué)模型在T、N方向具有較大的殘余誤差影響，相比MEO衛(wèi)星，GEO和IGSO衛(wèi)星的力學(xué)模型精度具有更大的改進(jìn)空間。

3.2 多天軌道合成效果

為了進(jìn)一步驗(yàn)證北斗衛(wèi)星多天軌道合成的效果，采用2013年6月1日到2013年6月26日共26 d的北斗衛(wèi)星實(shí)測數(shù)據(jù)，所用測站及處理策略與上一算例方案2相同。在附加偽隨機(jī)脈沖的條件下，分別計(jì)算了北斗衛(wèi)星3 d(S1)、5 d(S2)和7 d(S3)解軌道。便于對比起見，僅取多天解軌道的中間一天作為當(dāng)天最終軌道，采用文獻(xiàn)[17]給出的方法，通過統(tǒng)計(jì)連續(xù)兩天最終軌道在00:00:00(或24:00:00)的軌道差異值中位數(shù)，來評定定軌精度。圖4以圖示的形式對比了GEO、IGSO和MEO衛(wèi)星3 d 、5 d和7 d解軌道的三維軌道差異值中位數(shù)，表3則對3種解算方案中各衛(wèi)星徑向(R)、切向(T)、法向(N)和三維(P)的軌道差異值中位數(shù)進(jìn)行了統(tǒng)計(jì)。

圖4 多天軌道合成試驗(yàn)各星軌道精度對比圖Fig.4 Comparison of BeiDou satellites orbit precision in orbit combination

Tab.3 Statistics of BeiDou satellites orbit precision in orbit combination m

由圖4和表3的結(jié)果可見：

(1) 隨著定軌弧段的增長，大部分衛(wèi)星的定軌精度都獲得了不同程度的改善。但是，應(yīng)該引起注意的是，C02星軌道精度隨定軌弧段的增長反而有明顯的降低趨勢，具體原因尚待進(jìn)一步分析。

(2) 從三維軌道精度來看，由3 d定軌弧段增長到5 d時(shí)，GEO、IGSO、MEO平均定軌精度分別提高了23.6%、32.0%和48.6%；由5 d增長到7 d時(shí)，3類衛(wèi)星分別提高了13.0%、17.8%和12.1%。由此可見，定軌弧段由3 d增長到5 d對軌道精度的改善作用更為明顯。這一點(diǎn)在圖4中也有直觀反映，這是由于定軌弧段增長到一定限度對觀測構(gòu)形的改善作用就不再顯著，力學(xué)模型的精度也在一定程度上限制了定軌弧段的長度。此外，定軌弧段由3 d增長到5 d時(shí)，GEO、IGSO和MEO的軌道精度提高幅度依次增大，主要是由于多天數(shù)據(jù)的累積對“運(yùn)動”中的IGSO和MEO(特別是MEO)衛(wèi)星觀測結(jié)構(gòu)改善作用更大，而對“靜止”的GEO衛(wèi)星改善作用稍弱。但是，這種改善對觀測幾何結(jié)構(gòu)較差的GEO衛(wèi)星卻具有重要意義，表現(xiàn)在其軌道精度改進(jìn)的絕對量要明顯大于IGSO和MEO衛(wèi)星。當(dāng)定軌弧段由5 d增長到7 d時(shí)，3類衛(wèi)星軌道精度改善幅度的差異性已經(jīng)不甚明顯。

(3) 以3 d定軌弧長與7 d定軌弧長作對比，在R方向上，GEO、IGSO、MEO衛(wèi)星平均定軌精度分別提高了0.045、0.013和0.019 m，在T方向上，3類衛(wèi)星分別提高了0.215、0.119和0.078 m，在N方向上，3類衛(wèi)星分別提高了0.050、0.084和0.044 m。隨著定軌弧段的增長，3類衛(wèi)星T方向軌道精度改善最為明顯。這主要是由于導(dǎo)航衛(wèi)星的距離觀測量對軌道T方向幾何約束不足，通過多天數(shù)據(jù)累積，有助于改善T方向的觀測幾何約束。

在使用相同的軌道精度評定方法并采用相同長度的定軌弧長的條件下，以上結(jié)果與文獻(xiàn)[17]的定軌精度基本相當(dāng)，本文結(jié)果略優(yōu)。除了附加偽隨機(jī)脈沖對力學(xué)模型缺陷的彌補(bǔ)作用之外，還得益于文中使用了更多的測站數(shù)據(jù)。需要說明的是，以上軌道精度評定方法實(shí)質(zhì)上是以單天邊界處軌道的不連續(xù)性來統(tǒng)計(jì)定軌精度，由于動力法定軌的末端效應(yīng)，定軌弧段兩端軌道精度較低，因此以弧段中間一天邊界處軌道的不連續(xù)性統(tǒng)計(jì)獲得的軌道精度要略優(yōu)于重疊弧段的結(jié)果。此外，隨著定軌弧段的增長，相鄰兩個(gè)定軌弧段公用的觀測數(shù)據(jù)會不斷增多，這在一定程度上也會引起軌道內(nèi)符合精度升高。因此，算例中軌道精度隨定軌弧段增長而不斷提高，除得益于觀測數(shù)據(jù)的累積效應(yīng)之外，以上兩個(gè)因素也具有一定的影響作用。

4 結(jié) 論

針對目前北斗衛(wèi)星多天解軌道計(jì)算依靠常規(guī)方法效率低下的問題，本文深入研究了基于法方程疊加的多天軌道合成方法，同時(shí)，結(jié)合北斗衛(wèi)星部分?jǐn)z動力模型精度尚待進(jìn)一步精化的現(xiàn)實(shí)情況，在研究了密切元素和力學(xué)參數(shù)處理方法的基礎(chǔ)上，重點(diǎn)推導(dǎo)了偽隨機(jī)脈沖參數(shù)的合成方法，提出了附加偽隨機(jī)脈沖的北斗衛(wèi)星多天軌道合成方法。結(jié)合北斗實(shí)測數(shù)據(jù)，對相關(guān)方法進(jìn)行了分析驗(yàn)證，結(jié)果表明：在軌道多天合成中，單天邊界附加偽隨機(jī)脈沖對北斗衛(wèi)星軌道精度有明顯改進(jìn)作用，其中，3類衛(wèi)星T、N方向軌道精度的改善程度大于R方向，GEO和IGSO衛(wèi)星軌道精度的提高幅度大于MEO衛(wèi)星。偽隨機(jī)脈沖的主要作用是彌補(bǔ)力學(xué)模型精度的不足，以上試驗(yàn)結(jié)果在一定程度上也說明：北斗系統(tǒng)3類衛(wèi)星的力學(xué)模型在T、N方向具有較大的殘余誤差影響；相比MEO衛(wèi)星，GEO和IGSO衛(wèi)星的力學(xué)模型精度具有更大的改進(jìn)空間。此外，利用文中方法進(jìn)行北斗衛(wèi)星多天軌道合成，一定范圍內(nèi)，隨著定軌弧段的增長，北斗星定軌精度能夠得到不同程度的改善。定軌弧長由3 d增長到5 d對定軌精度的改善作用要明顯大于由5 d增長到7 d的情況，定軌弧長的增長對北斗3類衛(wèi)星T方向的軌道精度改善作用更大，但是C02星軌道精度隨定軌弧段的增長反而有明顯的降低趨勢，具體原因有待進(jìn)一步分析。

需要指出的是，利用本文方法在附加偽隨機(jī)脈沖條件下獲得的高精度軌道反過來可以用于精化北斗衛(wèi)星力學(xué)模型，其具體方法有待進(jìn)一步深入研究。

致謝：感謝全球連續(xù)監(jiān)測評估系統(tǒng)(iGMAS)信息工程大學(xué)分析中心對本文工作的幫助和支持。

[1] 楊元喜.北斗衛(wèi)星導(dǎo)航系統(tǒng)的進(jìn)展、貢獻(xiàn)與挑戰(zhàn)[J].測繪學(xué)報(bào),2010,39(1): 1-6.

YANG Yuanxi.Progress,Contribution and Challenges of Compass/BeiDou Satellite Navigation System[J].Acta Geodaetica et Cartographica Sinica,2010,39(1): 1-6.

[2] 冉承其.北斗衛(wèi)星導(dǎo)航系統(tǒng)的發(fā)展[R].武漢: 第四屆中國衛(wèi)星導(dǎo)航學(xué)術(shù)年會,2013.RAN Chengqi.Development of BeiDou Navigation Satellite System[R].Wuhan: CSNC 2013.

[3] 中國衛(wèi)星導(dǎo)航系統(tǒng)管理辦公室.北斗衛(wèi)星導(dǎo)航系統(tǒng)發(fā)展報(bào)告2.2[S].北京: 中國衛(wèi)星導(dǎo)航系統(tǒng)管理辦公室,2013: 1-10.China Satellite Navigation Office.Report on the Development of BeiDou Navigation Satellite System V2.2[S].Beijing: China Satellite Navigation Office,2013.

[4] 郝金明,劉偉平,楊力,等.北斗衛(wèi)星導(dǎo)航系統(tǒng)精密定軌技術(shù)研究現(xiàn)狀[J].測繪科學(xué)技術(shù)學(xué)報(bào),2015,32(3): 221-225,230.

HAO Jinming,LIU Weiping,YANG Li,et al.Status of Precise Orbit Determination of BeiDou Navigation Satellite System[J].Journal of Geomatics Science and Technology,2015,32(3): 221-225,230

[5] BEUTLER G,BROCKMANN E,HUGENTOBLER U,et al.Combining Consecutive Short Arcs into Long Arcs for Precise and Efficient GPS Orbit Determination Combining n Consecutive One-Day-Arcs into One n-days-arc[J].Journal of Geodesy,1996,70(5): 287-299.

[6] 姚宜斌.GPS軌道合成原理及其實(shí)現(xiàn)[J].武漢大學(xué)學(xué)報(bào)(信息科學(xué)版),2007,32(6): 510-514.YAO Yibin.Theory and Realization of GPS Orbit Integration[J].Geomatics and Information Science of Wuhan University,2007,32(11): 1487-1492.

[7] BROCKMANN E.Combination of Solutions for Geodetic and Geodynamic Applications of the Global Positioning System (GPS)[D].Berne: University of Bern,1996: 1-7.

[8] 劉偉平.北斗衛(wèi)星導(dǎo)航系統(tǒng)精密軌道確定方法研究[D].鄭州: 信息工程大學(xué),2014: 70-90.

LIU Weiping.Research on Precise Orbit Determination of BeiDou Navigation Satellite System[D].Zhengzhou: Information Engineering University,2014: 70-90.

[9] DONG D,HERRING T A,KING R W.Estimating Regional Deformation from a Combination of Space and Terrestrial Geodetic Data[J].Journal of Geodesy,1998,72(4): 200-214.

[10] 姚宜斌.GPS精密定位定軌后處理算法與實(shí)現(xiàn)[D].武漢: 武漢大學(xué),2004: 1-7.

YAO Yibin.Research on the Algorithm and Realization of Post-processing for GPS Precise Positioning and Orbit Determination[D].Wuhan: Wuhan University,2004: 1-7.

[11] 劉偉平,郝金明,錢龍,等.衛(wèi)星軌道的法方程疊加合成方法[J].測繪科學(xué)技術(shù)學(xué)報(bào),2014,31(1): 23-27.

LIU Weiping,HAO Jinming,QIAN Long,et al.Method of Satellite Orbit Combination Based on Normal Equation Stacking[J].Journal of Geomatics Science and Technology,2014,31(1): 23-27.

[12] 施闖,趙齊樂,李敏,等.北斗衛(wèi)星導(dǎo)航系統(tǒng)的精密定軌與定位研究[J].中國科學(xué): 地球科學(xué),2012,42(6): 854-861.SHI Chuang,ZHAO Qile,LI Min,et al.Precise Orbit Determination of BeiDou Satellites with Precise Positioning[J].Science China Earth Sciences,2012,42(6): 854-861.

[13] ZHAO Qile,GUO Jing,LI Min,et al.Initial Results of Precise Orbit and Clock Determination for COMPASS Navigation Satellite System[J].Journal of Geodesy,2013,87(5): 475-486.

[14] 劉偉平,郝金明,李建文,等.多GNSS融合的北斗衛(wèi)星精密定軌[J].測繪學(xué)報(bào),2014,43(11): 1132-1138.DOI: 10.13485/j.cnki.11-2089.2014.0186.

LIU Weiping,HAO Jinming,LI Jianwen,et al.Multi-GNSS Joint Precise Orbit Determination of BeiDou Navigation Satellite System[J].Acta Geodaetica et Cartographica Sinica,2014,43(11): 1132-1138.DOI: 10.13485/j.cnki.11-2089.2014.0186.

[15] HE Lina,GE Maorong,WANG Jiexian,et al.Experimental Study on the Precise Orbit Determination of the BeiDou Navigation Satellite System[J].Sensors,2013,13(3): 2911-2928.

[16] LIU Yang,LOU Yidong,SHI Chuang,et al.BeiDou Regional Navigation System Network Solution and Precision Analysis[C]∥China Satellite Navigation Conference (CSNC) 2013 Proceedings.Guangzhou: CSNC,2013: 173-186.

[17] STEIGENBERGER P,HUGENTOBLER U,MONTENBRUCK O,et al.Precise Orbit Determination of GIOVE-B Based on the CONGO Network[J].Journal of Geodesy,2011,85(6): 357-365.

[18] BEUTLER G,BROCKMANN E,GURTNER W,et al.Extended Orbit Modeling Techniques at the CODE Processing Center of the International GPS Service for Geodynamics (IGS): Theory and Initial Results[J].Manuscripta Geodaetica,1994,19: 367-386.

[19] 劉偉平,郝金明,田英國,等.一種偽隨機(jī)脈沖的快速參數(shù)估計(jì)方法[J].武漢大學(xué)學(xué)報(bào)(信息科學(xué)版),2015,40(11): 1487-1492.

LIU Weiping,HAO Jinming,TIAN Yingguo,et al.Fast Parameter Estimation for Pseudo-stochastic Pulse[J].Geomatics and Information Science of Wuhan University,2015,40(11): 1487-1492.

[21] 劉林.天體力學(xué)方法[M].南京: 南京大學(xué)出版社,1998: 45-47.

LIU Lin.Method of Celestial Mechanics[M].Nanjing: Nanjing University Press,1998: 45-47.

(責(zé)任編輯：叢樹平)

Solution Method and Precision Analysis of Multi-days Orbit Combination of BeiDou Satellites

LIU Weiping，HAO Jinming，YU Heli，TIAN Yingguo

School of Navigation and Aerospace Engineering，Information Engineering University，Zhengzhou 450001，China

Multi-days orbit of BeiDou satellites is gotten by the routine method of stacking observation data at present,which limits calculation efficiency.And moreover,the dynamic model of BeiDou satellites is not perfect because of the short running time,which also limits the precision of multi-days orbit.The method of multi-days BeiDou satellite orbit combination is present.In this method,normal equation stacking is used to combining some consecutive one-day BeiDou satellite orbits to one multi-days orbit,which can improve calculation efficiency.Meanwhile,pseudo-stochastic pulse is used to compensate the deficiency of dynamic model.The analysis shows that it is the pseudo-stochastic pulse at one-day boundary that improves the precision of BeiDou satellite orbits obviously.And the tangential and normal precision is improved more than the radial.The orbit precision of GEO and IGSO is improved more than that of MEO.And the precision of BeiDou satellite orbit can be improved by extending the orbit arc length with certain limits,which is especially in the tangential direction.

BeiDou navigation satellite system; precise orbit determination; normal equation stacking; orbit combination; pseudo-stochastic pulse

Major Projects of the Second Generation Satellite Navigation System in China (No.GFZX0301040308)

LIU Weiping (1986—)，male，PhD，lecturer，majors in precise satellite orbit determination and positioning.

劉偉平，郝金明，于合理，等.北斗衛(wèi)星多天軌道合成方法及其精度分析[J].測繪學(xué)報(bào)，2016,45(10)：1157-1164.

10.11947/j.AGCS.2016.20160055.

LIU Weiping，HAO Jinming，YU Heli，et al.Solution Method and Precision Analysis of Multi-days Orbit Combination of BeiDou Satellites[J].Acta Geodaetica et Cartographica Sinica,2016,45(10):1157-1164.DOI:10.11947/j.AGCS.2016.20160055.

P228

A

1001-1595(2016)10-1157-08

中國第2代衛(wèi)星導(dǎo)航系統(tǒng)重大專項(xiàng)(GFZX0301040308)

2016-01-31

修回日期：2016-06-21

劉偉平(1986—)，男，博士，講師，研究方向?yàn)樾l(wèi)星精密定軌定位。

E-mail：lwpchxy@sina.com