GPS卫星坐标用计算机计算,GPS卫星位置的计算
【摘 要】針對GPS衛(wèi)星導(dǎo)航電文的特點和里邊各種參數(shù)含義及用法,闡述了利用導(dǎo)航電文提供的數(shù)據(jù)計算衛(wèi)星位置的方法,為進行基線解算時提供準(zhǔn)確的衛(wèi)星坐標(biāo)數(shù)據(jù)。
0前言
全球衛(wèi)星定位系統(tǒng)(GPS)是一個實時的,全天候,全方位的定位系統(tǒng)。設(shè)計此系統(tǒng)的最初目的只是為美國軍事服務(wù),如今卻已擴展到諸多民用領(lǐng)域,包括民用導(dǎo)航與監(jiān)控,交通管理,測繪科學(xué)與技術(shù),土木工程,空間科學(xué)等等。隨著全球定位系統(tǒng)的不斷改進,硬、軟件的不斷完善,應(yīng)用領(lǐng)域正在不斷地開拓。
1 GPS衛(wèi)星導(dǎo)航電文
GPS衛(wèi)星的導(dǎo)航電文是用戶用來定位和導(dǎo)航的數(shù)據(jù)基礎(chǔ)。它主要包括:衛(wèi)星星歷、時鐘改正、電離層時延改正、工作狀態(tài)信息及C/A碼轉(zhuǎn)換到捕獲P碼的信息。這些信息是以二進制碼的形式,按規(guī)定格式組成,按幀向外播送,衛(wèi)星電文又叫數(shù)據(jù)嗎。它的基本單位是長1500bit的一個主幀,傳輸速率是50bit/s,30s傳送完畢一個主幀。一個主幀包括5個子幀,第1、2、3子幀各有10個字碼,每個字碼有30bit;第4、5子幀各有25個頁面,共有37500bit。第1、2、3子幀每30s重復(fù)一次,內(nèi)容每小時更新一次。第4、5子幀是12.5min播完一次,然后再重復(fù)之,其內(nèi)容僅在衛(wèi)星注入新的導(dǎo)航數(shù)據(jù)后才得以更新。
導(dǎo)航電文文件數(shù)據(jù)記錄節(jié)的格式說明如表1,表2,表3所示。
2 GPS衛(wèi)星位置的計算方法
第一步:計算衛(wèi)星運動的平均角速度n。
首先根據(jù)導(dǎo)航文件給出的參數(shù)計算出參考時刻TOE的平均角速度
:
(2.1)
GM=3.9860047×。
然后根據(jù)廣播星歷中給定的攝動參數(shù)n計算觀測時刻衛(wèi)星的平均角速度n
n=+n (2.2)
第二步:計算觀測瞬間衛(wèi)星的平近點角M
(2.3)
式中,為參考時刻TOE時的平近點角,由廣播星歷給出。
第三步:計算偏近點角
用弧度表示的開普勒方程為:
(2.4)
第四步:計算真近點角
(2.5)
(2.6)
(2.7)
第五步:計算升交距角
(2.8)
式中為近地點角距,由導(dǎo)航文件給出。
第六步:計算攝動改正項,,
導(dǎo)航文件中給出了下列6個攝動參數(shù), ,,,,,,據(jù)此可求出由于項引起的升交距角的攝動改正項、衛(wèi)星矢徑的攝動改正項和衛(wèi)星軌道傾角的攝動改正項。計算公式如下:
(2.9)
(2.10)
(2.11)
第七步:對、、進行攝動改正
(2.12)
式中:a為衛(wèi)星軌道的長半徑,a=,由導(dǎo)航文件獲得。為TOE時刻軌道傾角,由導(dǎo)航文件中的開普勒六參數(shù)給出。假設(shè)觀測時刻為2005年9月12日12時12分15秒,則從1981年開始算起,一直到2005年9月12日的天數(shù),用天數(shù)除以7,則余數(shù)即為本周開始時刻的天數(shù),再將余數(shù)與時分全部轉(zhuǎn)換為秒,即為觀測瞬間時刻。
第八步:計算衛(wèi)星在軌道平面坐標(biāo)系中的位置
在軌道平面直角坐標(biāo)系中,衛(wèi)星的平面直角坐標(biāo)為:
(2.13)
第九步:計算觀測瞬間升交點的經(jīng)度
若參考時刻TOE時升交點的赤經(jīng)為,升交點對時間的變化率為,那么觀測瞬間的升交點赤經(jīng)應(yīng)為:
(2.14)
可從導(dǎo)航文件攝動參數(shù)獲得。
設(shè)本周開始時刻(即星期日0時)格林尼治恒星時為,則觀測瞬間的格林尼治恒星時為:
(2.15)
式中:為地球自轉(zhuǎn)角速度,其值為=7.292115×rad/s;
這樣就可求得觀測瞬間升交點的經(jīng)度值為:
(2.16)
令
則有:
(2.17)
注意:導(dǎo)航文件里給出的并不是參考時刻TOE的升交點赤經(jīng),而是該值與本周起始時刻的格林尼治恒星時之差。
第十步:計算衛(wèi)星在瞬時地球坐標(biāo)系中的位置
已知升交點的大地經(jīng)度L以及軌道平面的傾角后,就可通過兩次旋轉(zhuǎn)方便地求得衛(wèi)星在地固坐標(biāo)系中的位置:
(2.18)
式中:
(2.19)
第十一步:計算衛(wèi)星在協(xié)議地球坐標(biāo)系中的位置
觀測瞬間衛(wèi)星在協(xié)議地球坐標(biāo)系中的位置
(2.20)
至此衛(wèi)星位置已求完。
3 結(jié)束語
衛(wèi)星的運動速度和加速度也可用類似方法計算。我們還可從精密星歷中獲得更高精度的衛(wèi)星坐標(biāo)。由于精密星歷文件是給出每15分鐘一個衛(wèi)星坐標(biāo),因此我們可以通過切比雪夫多項式擬合或內(nèi)插法求得觀測瞬間衛(wèi)星位置及運動加速度。
注:文章內(nèi)所有公式及圖表請用PDF形式查看。
總結(jié)
以上是生活随笔為你收集整理的GPS卫星坐标用计算机计算,GPS卫星位置的计算的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 【python二级】红楼梦
- 下一篇: [Luogu P4168] [BZOJ