崔先強(qiáng)，劉 青，賈小林

（1.中南大學(xué) 地球科學(xué)與信息物理學(xué)院，長(zhǎng)沙 410083；2.地理信息工程國(guó)家重點(diǎn)實(shí)驗(yàn)室，西安 710054；3.西安測(cè)繪研究所，西安 710054）

1 引言

在全球衛(wèi)星導(dǎo)航系統(tǒng)（global navigation satellite system，GNSS）中，怎樣快速地捕獲衛(wèi)星信號(hào)、選擇幾何精度因子最佳的衛(wèi)星組合進(jìn)行導(dǎo)航定位，并縮短啟動(dòng)導(dǎo)航定位的時(shí)間，一直是GNSS用戶所關(guān)注的熱點(diǎn)問(wèn)題。衛(wèi)星導(dǎo)航電文提供的歷書參數(shù)信息在接收機(jī)快速地捕獲衛(wèi)星信號(hào)和最佳選星的過(guò)程中起著十分重要的作用。在實(shí)際的導(dǎo)航定位啟動(dòng)過(guò)程中，如果沒有歷書以及相關(guān)的輔助設(shè)施來(lái)提供衛(wèi)星的概略位置，那么接收機(jī)只能進(jìn)行 “滿天搜索”，這將會(huì)大大增加接收機(jī)搜索導(dǎo)航衛(wèi)星的時(shí)間和增大誤搜率（可能搜索到其它衛(wèi)星導(dǎo)航系統(tǒng)的信號(hào)），這也將會(huì)對(duì)接收機(jī)的信號(hào)通道數(shù)提出更高的要求。GNSS系統(tǒng)歷書的主要作用是為用戶提供概略的衛(wèi)星位置，其主要用途表現(xiàn)在4個(gè)方面［1］：（1）有目的地搜索導(dǎo)航衛(wèi)星信號(hào)；（2）找到任意衛(wèi)星的概略多普勒頻移，輔助頻域搜索；（3）精密測(cè)量用戶對(duì)某地某時(shí)的衛(wèi)星可見性進(jìn)行預(yù)報(bào)；（4）用于計(jì)算衛(wèi)星的近似位置。

導(dǎo)航衛(wèi)星的歷書是精度較低的衛(wèi)星星歷，是導(dǎo)航電文所提供的基本數(shù)據(jù)，其對(duì)用戶的重要性卻不容忽視［2］。在2004年發(fā)布的全球定位系統(tǒng)（global positioning system，GPS）信號(hào)接口控制（interface control document，ICD）文檔中，介紹了兩種歷書參數(shù)，即中等精度的歷書（midi almanac）和簡(jiǎn)約歷書（reduced almanac）［3］。其中前者與2000年發(fā)布GPS ICD文檔［4］中介紹的歷書參數(shù)內(nèi)容相同，只是表示的位數(shù)減少了；而簡(jiǎn)約歷書不僅參數(shù)位數(shù)減少，而且參數(shù)構(gòu)成更加簡(jiǎn)單。在我國(guó)北斗衛(wèi)星導(dǎo)航系統(tǒng)（BeiDou navigation satellite system，BDS）建設(shè)過(guò)程中，國(guó)內(nèi)不少學(xué)者都對(duì)GPS的星歷與中等歷書參數(shù)及其用戶算法進(jìn)行了研究分析，據(jù)此設(shè)計(jì)出BDS的星歷與中等歷書參數(shù)及其用戶算法，并以投入使用。對(duì)星歷的研究主要表現(xiàn)在三個(gè)方面：首先，依據(jù)GPS星歷用戶算法導(dǎo)出了BDS中中圓地球軌道（medium earth orbits；MEO）和傾斜地球同步軌道（inclined geosynchronous orbits，IGSO）星歷參數(shù)的生成算法［5-6］；其次，針對(duì) BDS中地球靜止軌道（geostationary earth orbits，GEO）衛(wèi)星軌道傾角接近于零的情況提出了相應(yīng)的解決措施［7-8］；最后，分析了一些影響星歷參數(shù)擬合精度和效率的因素［9-10］。而歷書的研究成果主要是建立BDS歷書參數(shù)的生成算法［1，11］和針對(duì)參數(shù)擬合過(guò)程中存在的一些問(wèn)題提出了相應(yīng)的解決措施［12-13］。目前，所建立的星歷與歷書參數(shù)及其擬合算法已成功的用于我國(guó)BDS的運(yùn)行。

簡(jiǎn)約歷書是在2004年公布的GPS ICD文檔中給出的，目前除了文獻(xiàn) ［2］及文獻(xiàn) ［14］對(duì)其參數(shù)及其用戶算法進(jìn)行分析之外，還未見相關(guān)的文獻(xiàn)報(bào)導(dǎo)。為此，本文將對(duì)GPS簡(jiǎn)約歷書與中等歷書進(jìn)行簡(jiǎn)要地比較分析，并嘗試將其應(yīng)用于我國(guó)的BDS并設(shè)置相應(yīng)的常量參數(shù)，為我國(guó)BDS的簡(jiǎn)約歷書參數(shù)及其算法設(shè)計(jì)提供參考。

2 GPS簡(jiǎn)約歷書及其與中等歷書的比較

GPS的簡(jiǎn)約歷書參數(shù)包括4個(gè)，即歷書參考時(shí)刻toa、長(zhǎng)半軸與參考值之差δA、升交點(diǎn)赤經(jīng)Ω0和參考時(shí)刻的緯度幅角Φ0（Φ0＝M0＋ω），其中需要估計(jì)參數(shù)為后三個(gè)，相應(yīng)的長(zhǎng)半軸參考值為Aref＝26 559 710。相對(duì)于中等歷書而言，缺省參數(shù)的值分別為：e＝0、i＝55°和＝－2.6×10－9半周/秒，相應(yīng)的用戶算法計(jì)算步驟（與文獻(xiàn) ［2］類似）為：（1）改正長(zhǎng)半軸A＝Aref＋δA；（2）計(jì)算衛(wèi)星運(yùn)行的平均角速度（3）計(jì)算時(shí)刻與參考時(shí)刻之差tk＝t－toa；（4）獲取計(jì)算時(shí)刻的緯度幅角Φk＝Φ0＋ntk；（5）獲取計(jì)算時(shí)刻的升交點(diǎn)經(jīng)度（6）計(jì)算衛(wèi)星在地固系中的坐標(biāo)為

顯然，簡(jiǎn)約歷書的用戶算法中并沒有涉及到常量參數(shù)e，也就是說(shuō)e的取值變化并不會(huì)影響簡(jiǎn)約歷書的精度。相對(duì)于中等歷書用戶算法而言，簡(jiǎn)約歷書的用戶算法明顯更加簡(jiǎn)潔，用戶接收機(jī)在使用簡(jiǎn)約歷書參數(shù)計(jì)算衛(wèi)星位置所占用的內(nèi)存和時(shí)間都將會(huì)減少。

與中等歷書參數(shù)相比較，可以發(fā)現(xiàn)（見表1）：（1）相同的參數(shù)，簡(jiǎn)約歷書所占的位數(shù)要少，且其精度也更低，如Ω0；（2）簡(jiǎn)約歷書在衛(wèi)星軌道長(zhǎng)半軸的表示形式上有所變化，即選用了長(zhǎng)半軸與參考值之差，這與GPS 18參數(shù)［15］的表示形式時(shí)一致的；（3）簡(jiǎn)約歷書將近地點(diǎn)幅角ω和參考時(shí)刻的平近點(diǎn)角M0歸化為了一個(gè)參數(shù)Φ0；（4）將參數(shù)e、i和直接以常量的形式給出，不需要在導(dǎo)航電文中進(jìn)行發(fā)播，可以減小衛(wèi)星的內(nèi)存負(fù)荷。

表1 GPS 簡(jiǎn)約歷書與中等歷書對(duì)照比較表

3 BDS簡(jiǎn)約歷書參數(shù)的初步設(shè)計(jì)與生成算法

相對(duì)于美國(guó)的GPS而言，我國(guó)BDS的顯著特點(diǎn)是星座由GEO、IGSO和MEO三類衛(wèi)星共同組成，由于GEO和IGSO衛(wèi)星軌道特性與MEO衛(wèi)星并不相同，尤其是GEO衛(wèi)星軌道面傾角接近于零，導(dǎo)致了GEO衛(wèi)星的廣播星歷用戶算法與IGSO和MEO不同。對(duì)于BDS歷書參數(shù)的生成算法，在不旋轉(zhuǎn)軌道面的條件下通過(guò)實(shí)際計(jì)算可以發(fā)現(xiàn)，GEO的軌道面傾角接近于零并不會(huì)導(dǎo)致歷書參數(shù)擬合發(fā)散，而且其擬合精度損失能夠滿足用戶的實(shí)際需求，具體的計(jì)算結(jié)果見下一節(jié)算例。

為了使我國(guó)BDS能夠更好地與美國(guó)的GPS兼容，初步設(shè)計(jì)與GPS簡(jiǎn)約歷書相類似的參數(shù)，即包括歷書參考時(shí)刻toa、長(zhǎng)半軸與參考值之差δA、升交點(diǎn)赤經(jīng)Ω0和參考時(shí)刻的緯度幅角Φ0（Φ0＝M0＋ω），只是根據(jù)我國(guó)BDS的特殊性，需同時(shí)設(shè)置GEO、IGSO和MEO衛(wèi)星的軌道長(zhǎng)半軸參考值，且MEO衛(wèi)星軌道長(zhǎng)半軸不一定與GPS相同。與BDS的星歷參數(shù)生成算法相同的是，簡(jiǎn)約歷書也同樣涉及到作為常量的軌道面傾角不能直接賦值為零，否則將會(huì)導(dǎo)致GEO衛(wèi)星簡(jiǎn)約歷書生成算法發(fā)散，從而無(wú)法獲得相應(yīng)的簡(jiǎn)約歷書參數(shù)。因此，為了保證簡(jiǎn)約歷書參數(shù)擬合算法的收斂性，將GEO衛(wèi)星軌道面傾角直接賦值為1°，而MEO和IGSO衛(wèi)星軌道面傾角與 《北斗衛(wèi)星導(dǎo)航系統(tǒng)空間信號(hào)接口控制文件》［16］里的說(shuō)明相同，即取為55°。3類衛(wèi)星的偏心率相同，都設(shè)為0。對(duì)于Aref和，則根據(jù)BDS的精密軌道來(lái)擬合中等歷書參數(shù)，并取平均值來(lái)獲得。以GEO衛(wèi)星的長(zhǎng)半軸參考值為例進(jìn)行說(shuō)明：首先使用一周的精密衛(wèi)星軌道分別擬合出每個(gè)衛(wèi)星的中等歷書參數(shù)，然后取五個(gè)衛(wèi)星的長(zhǎng)半軸平均值作為GEO衛(wèi)星的參考值。IGSO、MEO衛(wèi)星的長(zhǎng)半軸參考值以及三類衛(wèi)星的升交點(diǎn)赤經(jīng)的變化率采用同樣的方法即可獲得，所選用的軌道與算例相同。下面將直接給出相應(yīng)的設(shè)定值，如表2所示。

表2 BDS簡(jiǎn)約歷書常量參數(shù)取值

對(duì)于簡(jiǎn)約歷書參數(shù)的生成算法，與廣播星歷和歷書的生成算法是相同的，不同的則是衛(wèi)星位置對(duì)相應(yīng)參數(shù)的偏導(dǎo)數(shù)，這里將只給出依據(jù)前面簡(jiǎn)約歷書的用戶算法推導(dǎo)出的偏導(dǎo)數(shù)表達(dá)式，即

4 計(jì)算與比較

本算例所用數(shù)據(jù)為BDS衛(wèi)星的精密軌道，采樣間隔為5min，所取時(shí)間段為2013－10－01—07，其中03號(hào)星在2013－10－07—09時(shí)有3d機(jī)動(dòng)，缺少7d的軌道，則選擇了2013－09－30—10－06的軌道。我國(guó)BDS在軌衛(wèi)星為14顆，由于提供的精密軌道中MEO 13－14號(hào)星沒有連續(xù)的7d軌道，因此這里只選用了12顆衛(wèi)星進(jìn)行試算，它們分別為：1～5號(hào)是GEO星，6～10號(hào)是IGSO星，11～12號(hào)是MEO星。由于同類衛(wèi)星的計(jì)算結(jié)果類似，這里將只給出1、6和11號(hào)星的擬合殘差曲線圖。曲線圖的橫坐標(biāo)以歷元來(lái)表示，即第一歷元對(duì)應(yīng)于擬合弧段的初始時(shí)刻，最后一個(gè)歷元對(duì)應(yīng)于擬合弧段的最后時(shí)刻。實(shí)際計(jì)算時(shí)，為了便于分析給定的參考值與常量參數(shù)發(fā)生變化對(duì)簡(jiǎn)約歷書的影響，算例中除了使用表2中提供的數(shù)據(jù)外，將對(duì)其進(jìn)行相應(yīng)調(diào)整，并將調(diào)整前后的簡(jiǎn)約歷書擬合誤差作差進(jìn)行分析，調(diào)整后的參數(shù)取值如表3。

表3 參考值與常量參數(shù)調(diào)整后的值

為了便于分析參考值和參數(shù)常量的變化對(duì)簡(jiǎn)約歷書參數(shù)擬合精度損失的影響，接下來(lái)將分別進(jìn)行討論，具體方案及其計(jì)算如果如下。

方案一，中等歷書的計(jì)算結(jié)果，如圖1～圖3。

方案二，表2設(shè)定值的簡(jiǎn)約歷書計(jì)算結(jié)果，如圖4～圖6。

方案三，Aref采用表3數(shù)據(jù)時(shí)與設(shè)定值的計(jì)算結(jié)果之差，如圖7。

方案四，i采用表3數(shù)據(jù)與設(shè)定值的計(jì)算結(jié)果之差，如圖8～圖10。

圖1 01號(hào)GEO衛(wèi)星中等歷書擬合誤差

圖2 061號(hào)IGSO衛(wèi)星中等歷書擬合誤差

圖3 11號(hào)MEO衛(wèi)星中等歷書擬合誤差

圖4 01號(hào)GEO衛(wèi)星簡(jiǎn)約歷書擬合誤差

圖5 06號(hào)IGSO衛(wèi)星簡(jiǎn)約歷書擬合誤差

圖6 11號(hào)MEO衛(wèi)星簡(jiǎn)約歷書擬合誤差

圖7 Aref變化前后擬合誤差之差（所有衛(wèi)星結(jié)果相同）

圖8 01號(hào)GEO衛(wèi)星i變化前后擬合誤差之差

圖9 06號(hào)IGSO衛(wèi)星i變化前后擬合誤差之差

圖10 11號(hào)MEO衛(wèi)星i變化前后擬合誤差之差

圖11 01號(hào)GEO衛(wèi)星變化前后擬合誤差之差

圖12 06號(hào)IGSO衛(wèi)星變化前后擬合誤差之差

圖13 11號(hào)MEO衛(wèi)星變化前后擬合誤差之差

由計(jì)算結(jié)果可以看出：

（1）將方案一（圖1～圖3）與方案二（圖4～圖6）結(jié)果進(jìn)行對(duì)比可以發(fā)現(xiàn)，中等歷書參數(shù)的精度要明顯地優(yōu)于簡(jiǎn)約歷書參數(shù)，原因在于簡(jiǎn)約歷書參數(shù)直接給出的3個(gè)參數(shù)i、e和的值與真實(shí)值之間存在較大的偏差。

（2）根據(jù)表2所設(shè)定的A參考值和常量參數(shù)所計(jì)算出的簡(jiǎn)約歷書參數(shù)的擬合殘差都在600km以內(nèi)，而且MEO衛(wèi)星的結(jié)果（圖6）要明顯地優(yōu)于文獻(xiàn) ［2］的計(jì)算結(jié)果，原因在于文獻(xiàn) ［2］中的簡(jiǎn)約歷書參數(shù)直接來(lái)自于廣播星歷，而本文則是采用最小二乘算法擬合得到，使得部分誤差得到了較好的控制。

（3）從方案三的誤差曲線圖（圖7）可以看出，選擇不同的長(zhǎng)半軸參考值所獲得的簡(jiǎn)約歷書參數(shù)擬合結(jié)果是相同的，即對(duì)長(zhǎng)半軸進(jìn)行調(diào)整前后的擬合殘差之差為0，這說(shuō)明長(zhǎng)半軸參考值設(shè)置的精確與否并不會(huì)影響整個(gè)簡(jiǎn)約歷書參數(shù)的精度。

（4）方案四的計(jì)算結(jié)果（圖8～圖10）顯示出，較小的傾角變化將會(huì)引起簡(jiǎn)約歷書的較大誤差，尤其是針對(duì)GEO衛(wèi)星傾角本就接近于0的情況，因此在設(shè)置簡(jiǎn)約歷書的常量參數(shù)i時(shí)，必須盡可能的逼近于真實(shí)軌道。

（5）從方案五的計(jì)算結(jié)果（圖11～圖13）可以看出，相對(duì)于簡(jiǎn)約的擬合殘差來(lái)說(shuō)，的較小波動(dòng)引起的簡(jiǎn)約歷書擬合殘差的變化并不是很大，尤其是在z方向的影響幾乎接近于0，而且GEO和IGSO衛(wèi)星水平方向的變化具有明顯的系統(tǒng)性傾向。

5 結(jié)束語(yǔ)

本文對(duì)GPS中等歷書和簡(jiǎn)約歷書進(jìn)行了比較分析，并基于此初步設(shè)計(jì)了我國(guó)BDS的簡(jiǎn)約歷書參數(shù)及其相應(yīng)參考值和常量參數(shù)，導(dǎo)出了衛(wèi)星位置對(duì)簡(jiǎn)約歷書參數(shù)的偏導(dǎo)數(shù)，分析了影響B(tài)DS簡(jiǎn)約歷書擬合精度的一些影響因素，可得出如下結(jié)論：

（1）在簡(jiǎn)約歷書的用戶算法中，并沒有涉及到軌道偏心率e，因此直接賦值為0即可。

（2）長(zhǎng)半軸參考值A(chǔ)ref存在一定的誤差并不會(huì)影響簡(jiǎn)約歷書參數(shù)的精度；

（3）軌道面傾角i的微小波動(dòng)可能會(huì)導(dǎo)致簡(jiǎn)約歷書擬合殘差的較大變化，尤其是實(shí)際傾角接近于零的GEO衛(wèi)星，因此在設(shè)計(jì)常量參數(shù)軌道面傾角i時(shí)必須與軌道真實(shí)值一致；

［1］陳留成，胡小工，韓春好，等.導(dǎo)航衛(wèi)星歷書參數(shù)擬合算法［J］.天文學(xué)報(bào)，2008，49（3）：288－296.

［2］陳南，陳大恒，賈小林.GPS新民用導(dǎo)航電文簡(jiǎn)約歷書的算法和性能分析［J］.測(cè)繪工程，2006，15（6）：63－66.

［3］IS－GPS－200，Revision D.Navstar GPS Space Segment/Navigation User Interface［S］.

［4］ICD－GPS－200，Revision C－004.Navstar GPS Space Segment/Navigation User Interface［S］.

［5］崔先強(qiáng)，焦文海，賈曉林，等.GPS廣播星歷參數(shù)擬合算法［J］.測(cè)繪學(xué)院學(xué)報(bào)，2004，21（4）：244－246.

［6］陳留成，韓春好，陳金平.廣播星歷參數(shù)擬合算法研究［J］.測(cè)繪科學(xué)，2007，32（3）：12－14、192.

［7］黃勇，胡小工，王小亞，等.中高軌衛(wèi)星廣播星歷精度分析［J］.天文學(xué)進(jìn)展，2006，24（1）：81－88.

［8］崔先強(qiáng)，楊元喜，吳顯兵.軌道面旋轉(zhuǎn)角對(duì) GEO衛(wèi)星廣播星歷參數(shù)擬合的影響［J］.宇航學(xué)報(bào)，2012，33（5）：590－596.

［9］崔先強(qiáng)，焦文海，秦顯平.GPS廣播星歷參數(shù)擬合算法的探討［J］.測(cè)繪科學(xué)，2006，31（1）：25－26，48.

［10］陳留成，唐波.參考系選擇對(duì) Kepler廣播星歷參數(shù)擬合精度的影響［J］.飛行器測(cè)控學(xué)報(bào)，2006，25（4）：19－25.

［11］崔先強(qiáng)，焦文海，秦顯平.導(dǎo)航衛(wèi)星的歷書參數(shù)及其擬合算法［J］.測(cè)繪科學(xué)技術(shù)學(xué)報(bào)，2006，23（1）：23－25.

［12］陳留成，唐波，梁紅梅.導(dǎo)航衛(wèi)星歷書參數(shù)擬合改進(jìn)算法［J］.天文學(xué)進(jìn)展，2008，26（1）：80－86.

［13］莫中秋，陳留成，董恩強(qiáng)，等.導(dǎo)航衛(wèi)星歷書擬合初值的選取與誤差分析［J］.測(cè)繪科學(xué)技術(shù)學(xué)報(bào)，2008，25（2）：104－107.

［14］陳南，賈小林，崔先強(qiáng).GPS民用導(dǎo)航電文CNAV的特點(diǎn)［J］.全球定位系統(tǒng)，2006（1）：1－6，34.

［15］崔先強(qiáng)，焦文海，賈曉林，等.兩種 GPS廣播星歷參數(shù)算法的比較［J］.空間科學(xué)學(xué)報(bào)，2006，26（5）：382－387.

［16］中國(guó)衛(wèi)星導(dǎo)航系統(tǒng)管理辦公室.北斗衛(wèi)星導(dǎo)航系統(tǒng)空間信號(hào)接口控制文件［EB/OL］.［2013－12－08］.http：//wenku.baidu.com/link？url＝PzhulclXTbE61jIdz6YkAa1B6－IlQE17LTZbngeJRMGQXiaGXhPeFcWmDNbrWd6GaI2oxmYLQx 6TZCxEWKO07Srj62N2i5EYbA2kVJmAFYm.