衛(wèi)星軌道計算從開普勒六根數(shù)到廣播星歷參數(shù)學(xué)習(xí)RTKLIB源碼時decode_eph函數(shù)中那一長串軌道參數(shù)eph-A、eph-e、eph-i0、eph-OMG0……讓人眼花繚亂。這些參數(shù)到底代表什么它們是如何共同確定一顆衛(wèi)星在太空中的位置的本文系統(tǒng)梳理廣播星歷中各個軌道參數(shù)的物理含義及其在衛(wèi)星位置計算中的作用。一、為什么需要廣播星歷廣播星歷是GNSS系統(tǒng)向用戶播發(fā)的衛(wèi)星軌道信息用戶接收機根據(jù)這些參數(shù)計算出任意時刻衛(wèi)星的精確位置。廣播星歷通常包含1個參考時刻、6個開普勒軌道參數(shù)和9個攝動改正參數(shù)。它的核心設(shè)計思路是以開普勒橢圓軌道為基準(zhǔn)再用調(diào)和項修正攝動影響。理解這些參數(shù)首先要從開普勒軌道六要素說起。二、開普勒軌道六要素描述衛(wèi)星軌道的“身份證”描述衛(wèi)星軌道最基本、最經(jīng)典的工具就是開普勒軌道六要素也叫軌道六根數(shù)?？梢园阉鼈兿胂蟪山o一個太空中的橢圓軌道“拍證件照”從六個維度來定義它。這六個參數(shù)可以分為三組第一組決定軌道的大小和形狀參數(shù)符號含義在RTKLIB中的字段軌道半長軸a橢圓軌道長軸的一半決定軌道大小和運行周期eph-A sqrt(a)存儲的是平方根偏心率e描述軌道偏離圓形的程度0為圓形越接近1越扁eph-e半長軸越大衛(wèi)星離地球越遠運行周期越長。GPS衛(wèi)星所在的MEO軌道半長軸約26560公里周期約12小時。GNSS衛(wèi)星的偏心率都極小0.01以內(nèi)軌道非常接近圓形。第二組決定軌道在空間中的朝向參數(shù)符號含義在RTKLIB中的字段軌道傾角i軌道平面與赤道面之間的夾角-eph-i0升交點赤經(jīng)Ω升交點衛(wèi)星從南半球穿入北半球的點與春分點之間的夾角-eph-OMG0這兩個參數(shù)共同確定了軌道平面在三維空間中的方位。軌道傾角決定了軌道面相對赤道面的傾斜程度升交點赤經(jīng)決定了軌道面在赤道面上的“朝向”。第三組決定衛(wèi)星在軌道上的當(dāng)前位置參數(shù)符號含義在RTKLIB中的字段近地點幅角ω從升交點到近地點的角度eph-omg平近點角M描述衛(wèi)星在參考時刻處于軌道上的哪個位置eph-M0近地點幅角決定了橢圓長軸在軌道平面內(nèi)的指向。平近點角則是一個“虛擬角度”通過開普勒方程迭代求解后才能得到衛(wèi)星的真實位置。三、9個攝動改正參數(shù)對抗“軌道漂移”實際衛(wèi)星并非嚴(yán)格運行在理想橢圓軌道上。地球非球形引力、日月引力、太陽光壓等攝動力會使衛(wèi)星軌道不斷偏離參考軌道。為了把精度從公里級提升到米級廣播星歷提供了9個攝動改正參數(shù)。3.1 長期攝動參數(shù)變化率參數(shù)含義在RTKLIB中的字段Δn平均角速度修正值eph-delnΩ?升交點赤經(jīng)變化率eph-OMGd?軌道傾角變化率eph-idot這三個參數(shù)描述軌道根數(shù)隨時間線性變化的速率用于補償長期攝動效應(yīng)。3.2 周期攝動參數(shù)諧波修正參數(shù)含義在RTKLIB中的字段C?c / C?s緯度幅角的余弦/正弦調(diào)和項振幅eph-cuc/eph-cusC?c / C?s軌道半徑的余弦/正弦調(diào)和項振幅eph-crc/eph-crs-C?c / C?s軌道傾角的余弦/正弦調(diào)和項振幅eph-cic/eph-cis這六個參數(shù)用于修正周期性的軌道攝動主要補償約6小時的短周期攝動。它們分別對緯度幅角u、軌道半徑r和軌道傾角i三個方向進行正弦/余弦諧波修正。修正公式以半徑為例r_corrected r Crc * cos(2u) Crs * sin(2u)四、參考時刻TOE與TOC參數(shù)含義在RTKLIB中的字段TOE星歷參考時刻Time of Ephemeriseph-toeTOC鐘差參考時刻Time of Clockeph-tocTOE是開普勒軌道參數(shù)的參考時刻TOC是衛(wèi)星鐘差多項式的參考時刻。兩者可能接近但來源和用途不同——TOE用于計算衛(wèi)星位置TOC用于計算衛(wèi)星鐘差。五、從參數(shù)到坐標(biāo)衛(wèi)星位置的計算流程理解了各個參數(shù)的含義再來看它們?nèi)绾巫罱K確定衛(wèi)星在ECEF坐標(biāo)系中的位置步驟1計算時間差tk t - TOE其中t是信號發(fā)射時刻TOE是星歷參考時刻。步驟2計算平近點角n n? Δn 平均角速度 理論值 修正值 M M? n * tk 平近點角 參考值 角速度 × 時間差步驟3開普勒方程迭代求解偏近點角 EM E - e * sin(E)這個方程沒有解析解需要通過牛頓迭代法求解。步驟4計算真近點角 vv 2 * arctan( sqrt((1e)/(1-e)) * tan(E/2) )步驟5計算緯度幅角 u?u? ω v步驟6加入攝動修正δu Cuc * cos(2u?) Cus * sin(2u?) δr Crc * cos(2u?) Crs * sin(2u?) δi Cic * cos(2u?) Cis * sin(2u?) u u? δu r a * (1 - e * cos(E)) δr i i? ? * tk δi步驟7計算軌道平面坐標(biāo)x r * cos(u) y r * sin(u)步驟8旋轉(zhuǎn)到ECEF坐標(biāo)系X x * cos(Ω) - y * cos(i) * sin(Ω) Y x * sin(Ω) y * cos(i) * cos(Ω) Z y * sin(i)其中Ω Ω? Ω? * tk考慮升交點赤經(jīng)變化率。六、各GNSS系統(tǒng)的差異系統(tǒng)星歷模型主要差異GPS開普勒攝動修正16個參數(shù)2小時更新一次-Galileo開普勒攝動修正與GPS類似3小時更新一次-BDS開普勒攝動修正與GPS兼容但北斗GEO衛(wèi)星有特殊處理-GLONASS初始狀態(tài)向量數(shù)值積分完全不同播發(fā)位置/速度初值用戶通過數(shù)值積分外推-GLONASS不采用開普勒參數(shù)而是直接播發(fā)衛(wèi)星在參考時刻的位置和速度初始條件用戶通過數(shù)值積分計算任意時刻的衛(wèi)星位置-。這也是為什么RTKLIB中GLONASS的星歷解碼函數(shù)是獨立的decode_geph。七、總結(jié)參數(shù)類別參數(shù)作用軌道大小a半長軸決定軌道尺寸和周期軌道形狀e偏心率決定橢圓扁平程度軌道朝向i傾角、Ω升交點赤經(jīng)決定軌道平面在空間中的方位軌道指向ω近地點幅角決定橢圓長軸指向衛(wèi)星位置M平近點角決定衛(wèi)星在軌道上的位置長期攝動Δn、Ω?、?補償軌道根數(shù)的線性漂移周期攝動Cuc/Cus、Crc/Crs、Cic/Cis補償周期性軌道偏差參考時刻TOE、TOC星歷和鐘差的基準(zhǔn)時間一句話總結(jié)6個開普勒參數(shù)定義了理想橢圓軌道3個變化率補償長期漂移6個諧波參數(shù)修正周期攝動——它們共同協(xié)作將衛(wèi)星位置的計算精度從公里級提升到米級。在RTKLIB的decode_eph和eph2pos函數(shù)中這些參數(shù)按照上述流程一步步計算出衛(wèi)星在ECEF坐標(biāo)系中的精確坐標(biāo)。本文基于RTKLIB源碼學(xué)習(xí)、GPS接口規(guī)范IS-GPS-200及公開技術(shù)資料整理如有不當(dāng)之處歡迎指正討論。