劉 成,杜曉輝,王 萌

(1.中國(guó)科學(xué)院國(guó)家天文臺(tái),北京100012;2.中國(guó)科學(xué)院研究生院,北京100190)

0 引 言

2002年,中國(guó)科學(xué)院的科研人員開(kāi)始研發(fā)有別于GPS類(lèi)衛(wèi)星定位系統(tǒng)的、創(chuàng)新的中國(guó)區(qū)域衛(wèi)星導(dǎo)航定位系統(tǒng)CAPS[1-3]。CAPS實(shí)驗(yàn)系統(tǒng)使用C波段商用靜止軌道(GEO)通信衛(wèi)星和小傾角傾斜同步軌道(SIGSO)的通信衛(wèi)星組成轉(zhuǎn)發(fā)式衛(wèi)星導(dǎo)航系統(tǒng),利用通信衛(wèi)星的轉(zhuǎn)發(fā)功能,將地面原子鐘生成的時(shí)間信號(hào)和導(dǎo)航電文通過(guò)衛(wèi)星轉(zhuǎn)發(fā)給用戶(hù),實(shí)現(xiàn)導(dǎo)航定位。

衛(wèi)星星座與布局是一個(gè)衛(wèi)星導(dǎo)航系統(tǒng)的重要基礎(chǔ),在很大程度上決定了系統(tǒng)的定位精度[4-5]。為補(bǔ)充和改善衛(wèi)星星座的分布,提高用戶(hù)的定位精度,CAPS系統(tǒng)創(chuàng)造了應(yīng)用氣壓測(cè)高輔助衛(wèi)星定位[6]的方法。在氣壓測(cè)高技術(shù)的輔助下,即使是衛(wèi)星導(dǎo)航星座全部采用赤道同步軌道(GEO)通信衛(wèi)星的CAPS-I期工程也實(shí)現(xiàn)了三維定位[1,9]。這樣,與發(fā)射專(zhuān)用導(dǎo)航衛(wèi)星的系統(tǒng)相比,建成CAPS系統(tǒng)資金投入較少,具有建設(shè)成本低、周期短和性能優(yōu)越等特點(diǎn)。然而,系統(tǒng)采用的氣壓測(cè)高,本身存在著誤差,也會(huì)給定位結(jié)果帶來(lái)誤差。因此,本文從觀測(cè)方程的解算方法出發(fā),分析了氣壓測(cè)高誤差對(duì)CAPS系統(tǒng)定位精度的影響,并通過(guò)仿真計(jì)算給出了兩者之間的量度關(guān)系,從而為氣壓測(cè)高技術(shù)下的實(shí)際定位工作提供了有意義的參考。

1 氣壓測(cè)高原理

氣壓測(cè)高方法,是利用地球大氣中氣壓隨高度變化的關(guān)系,在測(cè)量大氣壓的情況下獲得該處的地平高度。用戶(hù)接收到CAPS導(dǎo)航電文后,結(jié)合實(shí)測(cè)的氣壓和溫度值,即可求得該點(diǎn)高度值,進(jìn)入測(cè)量方程組進(jìn)行聯(lián)合解算,從而實(shí)現(xiàn)三維定位。

重力場(chǎng)內(nèi)的大氣壓力隨高度增加而減小,通過(guò)使用氣壓傳感器來(lái)測(cè)量氣象參數(shù),即可計(jì)算得到精確的相對(duì)高度[6-7]。

當(dāng)大氣處于流體靜力學(xué)平衡時(shí),按流體力學(xué)原理有:

式中:P為用戶(hù)所測(cè)的大氣壓力;P0為參考點(diǎn)大氣壓;h0為參考點(diǎn)高度;h為用戶(hù)設(shè)備高度;Rd為氣體參數(shù),取287.05J/(kg?K);T為大氣溫度(K);g為重力加速度。

重力加速度g隨高度變化非常緩慢,常作常數(shù)處理。所以上式可化為

由于大氣溫度隨高度的分布很復(fù)雜,難以用函數(shù)關(guān)系表示,因此,對(duì)式(2)直接求積分幾乎是不可能的。為此,通常假設(shè)在同樣的氣壓水平上大氣溫度為常數(shù),并取Tm(K)為P0和P之間大氣層中的平均溫度,則有:

換成以10為底的對(duì)數(shù),并用百分度表示溫度后,有:

式(4)稱(chēng)為L(zhǎng)aplace壓高方程,即為利用氣壓值測(cè)定高度的基本公式。其中,Tm=(T0+TR)/2,T0是基準(zhǔn)點(diǎn)溫度,TR是接收機(jī)測(cè)點(diǎn)測(cè)量值[6-8]。

我國(guó)境內(nèi)及周邊,共有1860個(gè)氣象站對(duì)居地的P0和T0值實(shí)行連續(xù)監(jiān)測(cè),CAPS主控站收集和處理氣象站數(shù)據(jù),以導(dǎo)航電文方式向用戶(hù)廣播。如果用戶(hù)位于氣象站附近,則氣象站的h0、P0和T0可直接用于計(jì)算用戶(hù)設(shè)備所在的高度。如果用戶(hù)設(shè)備離任何一個(gè)氣象站均較遠(yuǎn),則可以利用內(nèi)插法求得h0、P0和T0。

實(shí)驗(yàn)表明:氣象站氣壓高度計(jì)工作情況良好時(shí),中國(guó)大部分地區(qū)高度測(cè)量的精度與用戶(hù)所處高度有關(guān),如表1[9]所示:

表1 中國(guó)大部分地區(qū)高度測(cè)量精度

2 高程約束觀測(cè)方程的解算

CAPS導(dǎo)航定位系統(tǒng)中,氣壓測(cè)高技術(shù)的效果相當(dāng)于增設(shè)了一個(gè)虛擬星座,即相當(dāng)于一顆衛(wèi)星處于地球中心點(diǎn)附近,測(cè)得的高度值相當(dāng)于增加了一段偽距值,成為以地心距為半徑的另一個(gè)交會(huì)球,所以能起到一顆導(dǎo)航衛(wèi)星的作用[5]。這種情況下的觀測(cè)方程為

式中:h為用戶(hù)至地球基準(zhǔn)橢球的高度;a、b分別為地球基準(zhǔn)橢球的長(zhǎng)半軸和短半軸。式(5)中的最后一個(gè)方程式即為地球橢球約束方程。

求解式(5)的方法有多種[10-11],其中較為有效的一種是,地球橢球方程不直接參與方程組求解,而是單獨(dú)列出,僅作為變量之間的約束條件。將式(5)在用戶(hù)近似位置(^x,^y,^z)處經(jīng)線性化,按泰勒級(jí)數(shù)展開(kāi)為

具體求解步驟如下

1)首先可定義初始Δ z=0,將用戶(hù)估計(jì)位置處的^x和^y代入地球橢球方程,得^z;

2)由式(6),求得 Δ x,Δ y和Δ t;

3)由迭代格式X=^X+Δ X,得用戶(hù)坐標(biāo)x,y和接收機(jī)鐘差t;

4)將x和y代入地球橢球方程,得z,并由式Δ z=z-^z計(jì)算z坐標(biāo)增量;

5)重復(fù)上述步驟,直至 Δ x、Δ y、Δ z和 Δ t小于某一規(guī)定閾值[10-12]時(shí),迭代終止。

3 仿真與分析

中國(guó)區(qū)域定位系統(tǒng)CAPS-I中,四顆衛(wèi)星的位置坐標(biāo) Si,j(X,Y,Z)用地心地固直角坐標(biāo)系(ECEF)給出(單位:m),衛(wèi)星具體信息如表2所示。

表2 中國(guó)區(qū)域定位系統(tǒng)CAPS-I衛(wèi)星信息

衛(wèi)星C、D在傾斜軌道上的相位差約為120°(即過(guò)赤道的時(shí)間差約為8 h),仿真中取GM=3.986005e14;地球自轉(zhuǎn)角速度ω=7.2921151467e-5 rad/sec;WGS-84橢球長(zhǎng)半軸a=6378137 m;偏心率e=0.081819190842622。

記等效偽距測(cè)量均方根誤差為σρ,氣高測(cè)高均方根誤差為σh,三維坐標(biāo)位置定位均方根誤差為σu。仿真計(jì)算時(shí),針對(duì)10種不同的衛(wèi)星空間分布情況,在每一種空間分布的偽距值上加上均方根為特定值的隨機(jī)誤差(設(shè)定為2 m,即σρ=2 m),生成多組、大樣本的誤差分布數(shù)據(jù)。

3.1 觀測(cè)衛(wèi)星數(shù)目為四顆

3.1.1 高程約束算法定位結(jié)果

首先考察高程約束算法的定位解算誤差。在仿真中,設(shè)定氣壓測(cè)高的均方根誤差為σh=2 m,利用上述高程約束解算方法進(jìn)行100組絕對(duì)單點(diǎn)定位解算,并將定位誤差結(jié)果與經(jīng)典球面交會(huì)解算方法相比較,如圖1所示。

圖1 四顆衛(wèi)星時(shí)高程約束算法定位結(jié)果示意圖

從圖1中可以看出,高程約束算法明顯地改善了原球面交會(huì)解法的定位精度。統(tǒng)計(jì)圖1中的數(shù)據(jù)可知:原球面交會(huì)解算方法三維定位誤差為1σ=65.3 m,高程約束算法三維定位誤差為1σ=33.1 m。高程約束算法把原定位解算誤差減小了約50%。

3.1.2 氣壓測(cè)高誤差對(duì)定位結(jié)果的影響

實(shí)際定位工作中,由于用戶(hù)儀器設(shè)備精度和所處位置海拔高度不同等原因,測(cè)量得到的氣壓測(cè)高誤差σh也會(huì)不盡相同。為此,在仿真計(jì)算中令σh為各種不同的值(0～20 m),依次進(jìn)行高程約束定位解算,并統(tǒng)計(jì)出每個(gè)σh值所對(duì)應(yīng)的定位誤差均方根值σu,以分析氣壓測(cè)高誤差對(duì)最終定位結(jié)果的影響。如圖2所示。

圖2 四顆衛(wèi)星時(shí)高程誤差與定位解算結(jié)果關(guān)系示意圖

從圖2中可以看出,σh與σu呈現(xiàn)出大致的線性關(guān)系。

這個(gè)結(jié)論說(shuō)明,氣壓測(cè)高誤差σh和由它所產(chǎn)生的那部分定位誤差的量值是基本相同的,并不會(huì)放大最終的坐標(biāo)位置誤差。為建立誤差模型,可以利用樣條插值方法,對(duì)圖中的曲線進(jìn)行一次曲線擬合,得到氣壓測(cè)高誤差與定位誤差的關(guān)系式。一次曲線擬合的表達(dá)式可以寫(xiě)為

式中:定義σh0為σh=0 m時(shí)的三維定位誤差;μ為σh與σu之間的量度關(guān)系系數(shù)。在本情況下,σh0=30.70 m,μ=1.005。σh與σu之間的關(guān)系可表達(dá)為

式(8)即為用戶(hù)能夠觀測(cè)到A、B、C、D四顆衛(wèi)星時(shí),氣壓測(cè)高誤差對(duì)CAPS系統(tǒng)定位精度影響的近似計(jì)算表達(dá)式。在與本仿真情況相近的實(shí)際工作中進(jìn)行綜合定位誤差分析時(shí),μ即可以近似地取為1。

分別統(tǒng)計(jì)各坐標(biāo)方向上的誤差分量大小,結(jié)果如圖3所示。

圖3 四顆衛(wèi)星(σρ=2 m)時(shí)各坐標(biāo)方向上的誤差分量

從圖3中可以看出,在σh逐步增大的過(guò)程中,氣壓測(cè)高誤差對(duì)定位精度所帶來(lái)的影響,在Z坐標(biāo)方向上最大、X坐標(biāo)方向上最小。

3.2 觀測(cè)衛(wèi)星數(shù)目為三顆

3.2.1 高程約束算法定位結(jié)果

當(dāng)用戶(hù)只能觀測(cè)到上述四顆CAPS衛(wèi)星中的三顆時(shí),原本無(wú)法進(jìn)行三維定位工作;但在氣壓測(cè)高技術(shù)的輔助下,起到了一顆虛擬衛(wèi)星的作用,從而可以實(shí)現(xiàn)三維定位。只是星座DOP值較大,定位誤差也會(huì)較大。此時(shí),高程約束算法絕對(duì)單點(diǎn)定位結(jié)果如圖4所示。

圖4 三顆衛(wèi)星時(shí)高程約束算法定位結(jié)果示意圖

統(tǒng)計(jì)圖1中的數(shù)據(jù)可知,此時(shí)高程約束算法三維位置定位誤差約為1σ=53.6 m。

3.2.2 氣壓測(cè)高誤差對(duì)定位結(jié)果的影響

三顆觀測(cè)衛(wèi)星的情況下,同樣在仿真計(jì)算中令σh為各種不同的值(0～20 m),依次進(jìn)行高程約束定位解算,并統(tǒng)計(jì)出每個(gè)σh值所對(duì)應(yīng)的定位誤差均方根值σu,如圖5所示。

圖5 三顆衛(wèi)星時(shí)高程誤差與定位解算結(jié)果關(guān)系示意圖

同樣,通過(guò)一次樣條插值可以得到參數(shù) μ=0.953。此時(shí),σu與σh之間滿(mǎn)足下列誤差模型:

式(10)即為用戶(hù)能夠觀測(cè)到A、B、C、D四顆衛(wèi)星中的三顆時(shí),氣壓測(cè)高誤差對(duì)CAPS系統(tǒng)定位精度影響的近似計(jì)算表達(dá)式。

從圖4、圖5可以看到,即使在只有三顆觀測(cè)衛(wèi)星的情況下,CAPS-I也能通過(guò)氣壓測(cè)高的輔助作用順利實(shí)現(xiàn)三維定位。分別統(tǒng)計(jì)各坐標(biāo)方向上的誤差分量大小,結(jié)果如圖6所示。

可以看出,在σh逐步增大的過(guò)程中,氣壓測(cè)高誤差對(duì)定位精度的所帶來(lái)的影響,依然是在Z坐標(biāo)方向上最大、X坐標(biāo)方向上最小。

圖6 三顆衛(wèi)星(σρ=2 m)時(shí)各坐標(biāo)方向上的誤差分量

從以上各種情況下的仿真結(jié)果可知,無(wú)論觀測(cè)衛(wèi)星的數(shù)目是三顆還是四顆,也無(wú)論偽距測(cè)量值的誤差大小是多少,氣壓測(cè)高誤差總與三維位置定位誤差成近似的一次線性關(guān)系,且兩者之間的線性系數(shù)可以近似取1。換言之,氣壓測(cè)高誤差和由它所產(chǎn)生的那部分定位誤差的大小是基本相同的,并不對(duì)定位誤差產(chǎn)生額外的放大作用。仿真計(jì)算中每種情況下的統(tǒng)計(jì)均為大樣本數(shù)據(jù)統(tǒng)計(jì),所以各統(tǒng)計(jì)結(jié)果具有可靠的統(tǒng)計(jì)特性。

4 結(jié) 論

通過(guò)仿真計(jì)算與統(tǒng)計(jì)分析,我們可以得到以下若干結(jié)論:

1)無(wú)論觀測(cè)衛(wèi)星數(shù)目為三顆還是四顆,即使PDOP值較差,在氣壓測(cè)高的約束和輔助作用下,CAPS-I系統(tǒng)也都能順利地獲得測(cè)站三維位置定位結(jié)果,并且能夠?qū)⒄`差控制在一定范圍之內(nèi)。未來(lái),隨著CAPS系統(tǒng)衛(wèi)星星座的逐步完善與發(fā)展,相信會(huì)得到更好的定位結(jié)果。

2)無(wú)論觀測(cè)衛(wèi)星的數(shù)目為幾顆,也無(wú)論偽距測(cè)量值的誤差大小是多少,氣壓測(cè)高自身誤差對(duì)定位精度的影響總是比較獨(dú)立的,與三維位置定位誤差成近似的一次線性關(guān)系,且兩者之間的線性系數(shù)可以近似取為1。也就是說(shuō),氣壓測(cè)高誤差和由它所產(chǎn)生的那部分定位誤差的量值是基本相同的,并不會(huì)對(duì)最終的坐標(biāo)位置誤差起放大作用。

3)氣壓測(cè)高本身誤差對(duì)定位結(jié)果的影響,大部分體現(xiàn)在Z方向上。

4)在實(shí)際工作中,我們可以根據(jù)不同的衛(wèi)星星座及氣壓測(cè)高情況,統(tǒng)計(jì)和建立如上文所述的先驗(yàn)誤差模型及公式,從而可以指導(dǎo)CAPS等衛(wèi)星系統(tǒng)對(duì)氣壓測(cè)高約束下的定位結(jié)果進(jìn)行適當(dāng)?shù)母恼途?/p>

[1] Ai G X,Shi H L,Wu H T,et al.Positioning system based satellite communication and Chinese Area Positioning System(CAPS)[J].Chin J Astron Astrophys,2008,8(6):611-635.

[2] 艾國(guó)祥,施滸立,吳海濤,等.轉(zhuǎn)發(fā)器衛(wèi)星通信導(dǎo)航定位系統(tǒng),中國(guó):200410046064.1[P].2004-06-06.

[3] 施滸立,孫希延,李志剛.轉(zhuǎn)發(fā)式衛(wèi)星導(dǎo)航原理[M].北京:科學(xué)出版社,2009.

[4] Kaplan E D.Understanding GPS principlesand applications[M].Artech House London,UK.,1996.

[5] 周忠謨,易杰軍,周 琪.GPS衛(wèi)星測(cè)量原理與應(yīng)用[M].北京:測(cè)繪出版社,2004.

[6] 艾國(guó)祥,盛裴軒,杜金林,等.應(yīng)用于CAPS的氣壓測(cè)高虛擬衛(wèi)星[J].中國(guó)科學(xué)G輯:物理學(xué)?力學(xué)?天文學(xué),2008,38(12):1702-1710.

[7] 盛裴軒,毛節(jié)奏,李建國(guó),等.大氣物理學(xué)[M].北京:北京大學(xué)出版社,2003:30-35.

[8] United States committee on Extansion to the Standard Atmosphere.U.S.Standard Atmosphere[M].Washington:US Government Printing Office,1976.

[9] 艾國(guó)祥,施滸立,吳海濤,等.基于通信衛(wèi)星的定位系統(tǒng)原理[J].中國(guó)科學(xué)G輯:物理學(xué)?力學(xué)?天文學(xué),2008,38(12):1615-1633.

[10] 施滸立,裴 軍.轉(zhuǎn)發(fā)式衛(wèi)星導(dǎo)航定位系統(tǒng)量測(cè)方程解[J].中國(guó)科學(xué)G輯:物理學(xué)?力學(xué)?天文學(xué),2008,38(12):1687-1701.

[11] 施滸立,劉 成,耿建平,等.地球高程作為約束時(shí)的衛(wèi)星定位觀測(cè)方程求解方法的探討[G].CAPS內(nèi)部資料,2010.

[12] 紀(jì)元法,孫希延.中國(guó)區(qū)域定位系統(tǒng)的定位精度分析[J].中國(guó)科學(xué)G輯:物理學(xué)?力學(xué)?天文學(xué),2008,38(12):1812-1817.

[13] 施滸立,艾國(guó)祥,韓延本,等.退役衛(wèi)星的多生命周期利用[J].中國(guó)科學(xué)G輯:物理學(xué)?力學(xué)?天文學(xué),2008,38(12):1731-1737.