劉 永，張傳定，李?yuàn)檴櫍?琰

（1．信息工程大學(xué)，河南鄭州450001；2．北京衛(wèi)星導(dǎo)航中心，北京100000）

?

基于IGS精密星歷的BERNESE光壓模型參數(shù)分析

劉 永1，張傳定2，李?yuàn)檴?，王 琰1

（1．信息工程大學(xué)，河南鄭州450001；2．北京衛(wèi)星導(dǎo)航中心，北京100000）

摘 要：分析導(dǎo)航衛(wèi)星在光壓的作用下影響較大的光壓模型參數(shù)并尋求其變化規(guī)律。文中根據(jù)自編程序Orbit，對(duì)國(guó)外成熟的GPS光壓模型進(jìn)行實(shí)驗(yàn)分析，利用IGS精密星歷，通過(guò)動(dòng)力學(xué)擬合方法求解一年的BERN0光壓模型的九參數(shù)（D0，Y0，B0，Dc，Yc，Bc，Ds，Ys，Bs）并分析其變化規(guī)律。實(shí)驗(yàn)表明，D0，Y0，B0，Bc，Bs方向的光壓參數(shù)存在周期性變化，其余方向的參數(shù)變化規(guī)律性不強(qiáng)。運(yùn)用相關(guān)性參數(shù)分析B0和Dc；Yc和Bs；Ys和Bc；Ds和Y0 4組光壓參數(shù)強(qiáng)相關(guān)，每組中可考慮其中一個(gè)參數(shù)，結(jié)合5個(gè)參數(shù)的物理規(guī)律，只考慮D0，Y0，B0，Bc，Bs 5個(gè)參數(shù)是合理的，與武大五參數(shù)模型的理論相吻合。

關(guān)鍵詞：BERN0光壓模型；軌道擬合；相關(guān)性分析

太陽(yáng)光壓是GPS衛(wèi)星中最大的非保守力，成為影響其運(yùn)動(dòng)的最大攝動(dòng)力。在GPS衛(wèi)星的運(yùn)行過(guò)程中，為了獲得持續(xù)的能量供自身運(yùn)動(dòng)，在衛(wèi)星的兩側(cè)安裝了太陽(yáng)帆板，太陽(yáng)給衛(wèi)星提供能量的同時(shí)，還對(duì)衛(wèi)星產(chǎn)生一個(gè)力的作用，這個(gè)力就是太陽(yáng)光壓攝動(dòng)。太陽(yáng)光壓受很多因素的影響：太陽(yáng)帆板的材料選取、材料的老化、太陽(yáng)輻射強(qiáng)度、太陽(yáng)帆板受到照射的面積、太陽(yáng)與太陽(yáng)帆板受照面的夾角等。因此可見(jiàn)太陽(yáng)光壓成為衛(wèi)星在定軌過(guò)程中最難確定的非保守力。GPS衛(wèi)星發(fā)射前建立的分析性光壓模型，例如GPS的制造商Rockwell公司聯(lián)合IBM公司針對(duì)不同型號(hào)的衛(wèi)星建立了不同的ROCK光壓模型，伴隨分析數(shù)據(jù)的增多，為消除太陽(yáng)翼的羽流問(wèn)題建立T10／T20模型，并于1996年發(fā)展為T(mén)30光壓模型，之后，國(guó)外學(xué)者根據(jù)衛(wèi)星的長(zhǎng)期觀測(cè)數(shù)據(jù)反推經(jīng)驗(yàn)型的光壓模型，通過(guò)建立HILL方程分析衛(wèi)星運(yùn)動(dòng)過(guò)程中的軌道切向、徑向、法向的變化，得出Colombo模型；德國(guó)的伯爾尼大學(xué)通過(guò)歐洲定軌中心提供的數(shù)據(jù)建立了3個(gè)BERNESE模型，3種經(jīng)驗(yàn)性模型在GPS衛(wèi)星的定軌中應(yīng)用廣泛，效果較好。為彌補(bǔ)經(jīng)驗(yàn)型模型大數(shù)據(jù)問(wèn)題和分析性模型應(yīng)用精度不高的問(wèn)題，建立了JPL模型和Adjustable box－wing模型。我國(guó)目前在建BDS系統(tǒng)，尚未建立適合BDS的光壓模型，而光壓模型的精度對(duì)最終衛(wèi)星定軌的精度影響較大，因此若要提高定軌精度，開(kāi)展光壓模型研究意義重大。

本文借鑒國(guó)外GPS光壓模型研究過(guò)程，通過(guò)對(duì)模型的實(shí)驗(yàn)分析，確定出GPS衛(wèi)星在光壓的作用下，哪幾個(gè)方向參數(shù)影響較大，有規(guī)律可循，進(jìn)而為以后建立北斗衛(wèi)星的光壓模型奠定基礎(chǔ)。

1　常見(jiàn)軟件中使用的光壓模型公式

目前最適合GPS衛(wèi)星定軌的是經(jīng)驗(yàn)型光壓模型：伯爾尼大學(xué)建立的BERNESE模型。BERNESE模型有BERN0、BERN1、BERN2 3種模型，目前的幾個(gè)定軌軟件系統(tǒng)中都是依此模型為基礎(chǔ)變化出來(lái)的，而且各個(gè)分析中心在具體使用時(shí)形式不一樣，通過(guò)研究各個(gè)分析中心的光壓模型的不同性，分析BERNESE九參數(shù)模型在各個(gè)軸上的長(zhǎng)期項(xiàng)以及周期項(xiàng)的變化規(guī)律，最重要的是這9個(gè)參數(shù)沒(méi)有必要全部解算。本文特別針對(duì)經(jīng)驗(yàn)型模型中BERNESE模型進(jìn)行實(shí)驗(yàn)，通過(guò)對(duì)伯爾尼大學(xué)提供的BERNESE九參數(shù)中的BERN0模型，進(jìn)行軌道擬合，反算出接近一年的光壓九參數(shù)，對(duì)光壓九參數(shù)進(jìn)行分析，探究出不同軸方向的變化規(guī)律。在D0，Y0，B0，Bc，Bs 5個(gè)參數(shù)方向上存在周期性的規(guī)律，即這5個(gè)參數(shù)存在物理性的意義。接下來(lái)從相關(guān)性參數(shù)分析角度驗(yàn)證五參數(shù)模型存在的合理性，即B0和Dc，Yc和Bs，Ys和Bc，Ds和Y0 4組光壓參數(shù)存在強(qiáng)相關(guān)性，從數(shù)學(xué)角度驗(yàn)證減少參數(shù)的合理性，即在提高計(jì)算效率，不影響定軌的精度的前提下，為定軌中解算光壓參數(shù)的選取提供了理論基礎(chǔ)。

首先介紹光壓模型基本概念：3軸方向上存在周期性攝動(dòng)，參數(shù)包括3個(gè)軸方向上的輻射壓參數(shù)（D0，Y0，B0）和3個(gè)軸方向上的周期項(xiàng)攝動(dòng)系數(shù)（Dc，Yc，Bc，Ds，Ys，Bs）。式中，SRP（i），（i＝1，2，…，9）為待估參數(shù)，分別對(duì)應(yīng)BERN0光壓模型的九參數(shù)D0，Y0，B0，Dc，Yc，Bc，Ds，Ys，Bs；定義u為衛(wèi)星在軌道平面上距升交點(diǎn)的角度；λ為陰影因子；au代表天文單位的長(zhǎng)度；為加速度分別為慣性系下太陽(yáng)和衛(wèi)星的位置；D0為與衛(wèi)星型號(hào)有關(guān)的面質(zhì)比參數(shù)，其取值與衛(wèi)星型號(hào)以及質(zhì)量有關(guān)。為星固坐標(biāo)系坐標(biāo)軸的單位向量，其中指向地球中心，為太陽(yáng)至衛(wèi)星的單位向量，并且有＝×，＝×，＝×。

1．1 BERNESE

Bernese軟件中使用的光壓模型公式為

其中：

1．2 GAMIT

GAMIT軟件中光壓模型公式為

其中：

1．3 MADOCA

MADOCA軟件中光壓模型的公式為

式中：f為衛(wèi)星到衛(wèi)星在地球完全遮擋的中間位置的角度。

1．4 PANDA

PANDA軟件中使用的光壓模型公式為

其中：

從式（1）～（7）可以看出：

1）MADOCA軟件考慮了9個(gè)參數(shù)，GAMIT軟件也考慮了9個(gè)參數(shù)，但PANDA軟件只考慮了5個(gè)參數(shù)，因此對(duì)比來(lái)看：9個(gè)參數(shù)也不必全部考慮，從公式看出D0，Y0，B0，Bc，Bs 5個(gè)參數(shù)是必須要考慮的。在保證精度不受影響的前提下，參與計(jì)算的參數(shù)越少計(jì)算效率越高。

2）GAMIT軟件把地影因子放到D0方向上，MADOCA光壓模型和PANDA光壓模型把地影因子放到D，Y，B 3個(gè)方向上都考慮。從地影因子的位置看出：太陽(yáng)光壓已不是單純的光壓攝動(dòng)力，還包括其他的力的影響。因此從實(shí)驗(yàn)得到對(duì)光壓系數(shù)的反演，分析光壓系數(shù)的規(guī)律，探究出哪幾個(gè)方向的光壓系數(shù)是不可或缺的，為以后的北斗衛(wèi)星的光壓模型研究奠定理論基礎(chǔ)。因此通過(guò)實(shí)驗(yàn)來(lái)驗(yàn)證只解D0，Y0，B0，Bc，Bs 5個(gè)參數(shù)從物理角度存在的合理性。

2　參數(shù)相關(guān)性檢核

通過(guò)對(duì)一年的光壓九參數(shù)變化規(guī)律分析，從物理角度分析解哪幾個(gè)光壓參數(shù)的合理性，接下來(lái)要對(duì)光壓參數(shù)之間的相關(guān)性進(jìn)行分析，從數(shù)學(xué)角度驗(yàn)證只解算D0，Y0，B0，Bc，Bs五參數(shù)的合理性。參數(shù)相關(guān)性檢核從法方程出發(fā)，計(jì)算權(quán)逆陣，利用權(quán)逆陣輸出

其中：Qij為權(quán)逆陣第i行，第j列元素；Qi，Qj為權(quán)逆陣對(duì)角線第i，j元素。對(duì)于相關(guān)性檢核，從權(quán)逆陣計(jì)算相關(guān)系數(shù)，相關(guān)系數(shù)為1，兩個(gè)光壓參數(shù)線性相關(guān)；相關(guān)系數(shù)為0，兩個(gè)光壓參數(shù)相互獨(dú)立。

3　實(shí)驗(yàn)結(jié)果與分析

3．1 動(dòng)力學(xué)軌道平滑

動(dòng)力學(xué)軌道平滑的計(jì)算過(guò)程如下：

1）將單日IGS精密星歷文件中GPS的PRN02星的坐標(biāo)從ITRF坐標(biāo)系轉(zhuǎn)換到GCRS坐標(biāo)系，共得到96個(gè)歷元的三維坐標(biāo)x，y，z值。

2）設(shè)置衛(wèi)星的初始坐標(biāo)、初始速度和太陽(yáng)光壓參數(shù)（BERNE0九參數(shù)），共15個(gè)參數(shù)作為軌道參數(shù)，以96點(diǎn)GCRS坐標(biāo)作為虛擬的觀測(cè)值，用最小二乘法求解衛(wèi)星軌道參數(shù)。求解過(guò)程中，衛(wèi)星的初始坐標(biāo)和初始速度的近似值用二體問(wèn)題解算，光壓參數(shù)的初始值設(shè)為0。解算出軌道參數(shù)之后，再將其作為參數(shù)近似值，重新平差計(jì)算，直至參數(shù)的變化小于設(shè)定值為止（初始坐標(biāo)變化量小于0．01m）。

3）將軌道擬合后的數(shù)據(jù)進(jìn)行分析，進(jìn)行動(dòng)力學(xué)軌道平滑效果圖的驗(yàn)證（見(jiàn)圖1），通過(guò)點(diǎn)位誤差判斷出PRN02星的軌道擬合狀況，驗(yàn)證PRN02星在GCRS坐標(biāo)系下三維位置誤差及該星的運(yùn)行狀況。

4）將該星的光壓參數(shù)輸出，分析光壓參數(shù)的變化。

5）求解一年的BERN0模型九參數(shù)，分析存在的變化規(guī)律。

圖1為自編程序驗(yàn)證了PRN01－PRN32共32顆星的動(dòng)力學(xué)軌道平滑狀況。橫坐標(biāo)描述的是一天的精密星歷96個(gè)歷元，縱坐標(biāo)描述的是點(diǎn)位的范圍。

圖1　點(diǎn)位誤差比較

3．2 光壓參數(shù)變化規(guī)律分析

利用IGS精密星歷，選取其中的PRN02星，其一年數(shù)據(jù)的時(shí)間范圍為2008年12月14日至2009 年12月16日。利用3．1節(jié)中介紹的分析方法，可以得到9個(gè)方向的光壓參數(shù)一年的隨時(shí)間變化規(guī)律圖像（見(jiàn)圖2～圖10）。圖像中：橫坐標(biāo)是時(shí)間；縱坐標(biāo)是光壓參數(shù)數(shù)值，不同方向的數(shù)值的量級(jí)不一樣，每個(gè)光壓參數(shù)不再是單純的加速度，為方便計(jì)算，包含面值比參數(shù)。

圖2　D0方向變化趨勢(shì)

從圖2中可以看出，D0方向的參數(shù)存在周期變化規(guī)律，變化周期大約為半年，并且在峰頭處存在一個(gè)突變，光壓參數(shù)的變化范圍是（2 193．67，2 253．08）。

圖3　Dc方向變化趨勢(shì)

從圖3可以看出，Dc方向的參數(shù)存在無(wú)規(guī)律的變化，光壓參數(shù)的變化范圍是（－9．59，10．78）。

圖4　Ds方向變化趨勢(shì)

從圖4可以看出，Ds方向的參數(shù)存在無(wú)規(guī)律的變化，光壓參數(shù)的變化范圍是（－10．55，7．96）。

圖5　Y0方向變化趨勢(shì)

從圖5可以看出，Y0方向的參數(shù)存在規(guī)律性的變化，周期大約是一年，光壓參數(shù)的變化范圍是（－6．70，12．18）。

從圖6可以看出，Yc方向的參數(shù)存在無(wú)規(guī)律的變化，光壓參數(shù)的變化范圍是（－6．70，5．31）。

從圖7可以看出，Ys方向的參數(shù)變化規(guī)律性不明顯，光壓參數(shù)的變化范圍是（－1．65，5．12）。

圖6　Yc方向變化趨勢(shì)

圖7　Ys方向變化趨勢(shì)

圖8　B0方向變化趨勢(shì)

圖9　Bc方向變化趨勢(shì)

圖10　Bs方向變化趨勢(shì)

從圖8可以看出，B0方向的參數(shù)存在規(guī)律的變化，周期大約是半年，光壓參數(shù)的變化范圍是（－6．78，28．67）。

從圖9可以看出，Bc方向的參數(shù)存在規(guī)律性變化，光壓參數(shù)的變化范圍是（－48．14，40．73）。

從圖10可以看出，Bs方向的參數(shù)存在規(guī)律性變化，光壓參數(shù)的變化范圍是（－32．14，34．57）。

3．3 相關(guān)系數(shù)驗(yàn)證參數(shù)的相關(guān)性

利用IGS精密星歷（2009年12月15日），選取PRN02星輸出衛(wèi)星光壓參數(shù)的相關(guān)性矩陣見(jiàn)表1。

表1　相關(guān)系數(shù)權(quán)陣

相關(guān)系數(shù)大于0．8的有B0和Dc，Yc和Bs，Ys 和Bc，Ds和Y0 4組光壓參數(shù)。

從圖2～10可以看出：

1）D0方向（即對(duì)日方向）的光壓參數(shù)存在周期性的變化，通過(guò)分析數(shù)據(jù)得出的結(jié)論是在地球繞太陽(yáng)的運(yùn)動(dòng)的近日點(diǎn)時(shí)參數(shù)變化較快，在遠(yuǎn)日點(diǎn)時(shí)參數(shù)變化較慢，但呈現(xiàn)基本接近的周期性變化趨勢(shì)，周期是半年。與地球的運(yùn)動(dòng)速度有關(guān)，即運(yùn)動(dòng)速度快，參數(shù)變化明顯，反之亦然。

2）從光壓參數(shù)的量級(jí)和變化范圍來(lái)看，D0方向的量級(jí)是其它參數(shù)1 000倍，變化幅度較大，所以必須考慮；同時(shí)Y0，B0，Bc，Bs 4個(gè)光壓參數(shù)量級(jí)與剩余的光壓參數(shù)量級(jí)相同，但是變化范圍較大，所以光壓參數(shù)的選取時(shí)有必要考慮D0，Y0，B0，Bc，Bs 5個(gè)參數(shù)。

3）D0，Y0，B0，Bc，Bs數(shù)值是有規(guī)律性變化的，即這5個(gè)光壓參數(shù)存在一定規(guī)律性變化，而其他方向的參數(shù)存在無(wú)規(guī)律的抖動(dòng)。對(duì)光壓參數(shù)的選取有重要的物理意義，此5個(gè)參數(shù)一年的規(guī)律性變化意味著有物理意義。根據(jù)一年光壓模型的九參數(shù)的變化規(guī)律，把光壓模型的參數(shù)固定為包含D0， Y0，B0，Bc，Bs的五參數(shù)。

從表1可以看出：上面4組光壓參數(shù)相關(guān)性較強(qiáng)，每組中可以考慮其中一個(gè)參數(shù)，可以簡(jiǎn)化計(jì)算，結(jié)合物理的變化規(guī)律，從相關(guān)性角度驗(yàn)證了D0，Y0，B0，Bc，Bs 5個(gè)光壓參數(shù)必須考慮的合理性。

4　結(jié) 論

實(shí)驗(yàn)表明：①對(duì)于BERN0光壓模型，D0，Y0，B0，Bc，Bs 5個(gè)光壓參數(shù)存在周期性變化，D0其變化存在半年周期項(xiàng)；且地球在近日點(diǎn)時(shí)衛(wèi)星光壓參數(shù)變化較快，地球在遠(yuǎn)日點(diǎn)時(shí)衛(wèi)星光壓參數(shù)變化較慢，但是其余參數(shù)存在不規(guī)律的變化；②B0和Dc；Yc和Bs；Ys和Bc；Ds和Y0 4組光壓參數(shù)中每組兩參數(shù)間存在強(qiáng)相關(guān)，每組中只需考慮一個(gè)參數(shù)即可。結(jié)合5個(gè)參數(shù)的物理規(guī)律，由此得出只考慮D0，Y0，B0，Bc，Bs 5個(gè)參數(shù)不考慮其他參數(shù)是合理的，與PANDA軟件只考慮五參數(shù)相吻合。由于本文只選擇了單顆GPS衛(wèi)星一年的數(shù)據(jù)，接下來(lái)的工作可以進(jìn)行同類(lèi)型衛(wèi)星光壓參數(shù)分析、計(jì)算效率、模型精度等問(wèn)題展開(kāi)研究。

參考文獻(xiàn)：

［1］ 陳俊平，王解先．GPS定軌中的太陽(yáng)輻射壓模型［J］．天文學(xué)報(bào)，2006，47（3）：310－319．

［2］ DACH R，HUGENTOBLER U，F(xiàn)RIDEZ P，et al．Bernese GPS Software Version 5．0［M］．Astronomical Institute，University of Bern．January 2007．

［3］ 李濟(jì)生．人造衛(wèi)星精密軌道確定［M］．北京：解放軍出版社，1995．

［4］ 王琰．北斗導(dǎo)航衛(wèi)星精密定軌算法與軟件設(shè)計(jì)研究［D］．鄭州：信息工程大學(xué)，2013．

［5］ 吳顯兵．星載GPS低軌衛(wèi)星幾何法定軌及動(dòng)力學(xué)平滑方法研究［D］．鄭州：信息工程大學(xué)，2004．

［6］ 劉基余．GPS衛(wèi)星導(dǎo)航的定位原理與方法［M］．北京：科學(xué)出版社，2003．

［7］ 黃維彬．近代平差理論及其應(yīng)用［M］．北京：解放軍出版社，1990．

［8］ 王琰，宋力杰，黃令勇．GPS衛(wèi)星精密定軌中的攝動(dòng)力分析［J］．測(cè)繪工程，2013，22（5）：16－20．

［9］ 宋小勇，毛悅，賈小林．BERNESE光壓模型參數(shù)的統(tǒng)計(jì)分析［J］．測(cè)繪科學(xué)，2009，34（3）：25－27．

［10］陳秋麗，陳忠貴，王海紅．GPS衛(wèi)星光壓攝動(dòng)建模發(fā)展與啟示［J］．航天器工程，2013（2）：98－103．

［11］焦寧，田龍華，孫秀寧．廣義延拓插值法在IGS精密鐘差插值中的應(yīng)用［J］．測(cè)繪與空間地理信息，2015，38 （4）：144－146．

［12］潘紹林，張顯云，杜寧，等．基于灰色系統(tǒng)模型的IGS精密鐘差預(yù)報(bào)［J］．測(cè)繪工程，2015，24（4）：20－24．

［13］徐克紅，程鵬飛，文漢江．IGS跟蹤站的大地高時(shí)間序列特征分析［J］．測(cè)繪工程，2015，24（5）：19－23．

［責(zé)任編輯：劉文霞］

Analysis of BERN SRP parameters based on IGS ephemeris

LIU Yong1，ZHANG Chuanding2，LI Shanshan1，WANG Yan1

（1．Information Engineering University，Zhengzhou 450001，China；2．Beijing Satellite Navigation Center，Beijing 100000，China）

Abstract：It analyses the foreign GPS SRP（solar radiation parameter）model to seek the directions that SRP parameters have the greatest impact and the parameter variation．And it conducts the orbital dynamics fitting to get the year’s BERN0SRP nine parameters include D0／Y0／B0／Dc／Yc／Bc／Ds／Ys／Bs and the parameter variation with IGS precise ephemeris by the program Orbit．The analysis shows：the SRP parameters in D0／Y0／B0／Bc／Bs directions have cycle regularity but other directions don’t；it conducts the relative analysis of 4pairs of B0／Dc，Yc／Bs，Ys／Bc，Ds／Y0and finds strong correlation．Wuda Model only considers five parameters，D0／Y0／B0／Bc／Bs，as a result of which is reasonable．

Key words：BERN0model；orbital fitting；correlation analysis

作者簡(jiǎn)介：劉 永（1990－），男，碩士研究生．

基金項(xiàng)目：國(guó)家自然科學(xué)基金資助項(xiàng)目（41174026）

收稿日期：2015－06－23

中圖分類(lèi)號(hào)：P228

文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：1006－7949（2016）01－0029－05