本發(fā)明屬于衛(wèi)星導(dǎo)航定位技術(shù)領(lǐng)域,尤其涉及一種基于低軌衛(wèi)星的星地差分實時精密定位方法。
背景技術(shù):
衛(wèi)星導(dǎo)航定位系統(tǒng)(gnss)能夠提供全球范圍內(nèi)實時的定位服務(wù),已經(jīng)深刻的改變了我們的生活方式,并且在世界各國的許多行業(yè)中得到了廣泛的應(yīng)用。目前主流的全球定位系統(tǒng)有美國的gps,俄羅斯的glonass,中國的beidou系統(tǒng)以及歐洲的galileo系統(tǒng),印度也在積極發(fā)展區(qū)域?qū)Ш较到y(tǒng)irnss。目前標(biāo)準(zhǔn)的gnss定位技術(shù)精度大約是5-10米。對于精度要求更高的應(yīng)用,通常需要使用精密定位的方法。目前精密定定位的方法有差分定位方法和精密單點定位方法。其中差分定位方法又包括局域差分方法和廣域差分方法。廣域差分方法和精密單點定位都需要通過地面監(jiān)測網(wǎng)的數(shù)據(jù)計算出各種精密的信號偏差產(chǎn)品,精密軌道產(chǎn)品等,然后在用戶端使用這些精密產(chǎn)品進(jìn)行改正來提高定位精度。局域差分方法主要是將參考站的觀測數(shù)據(jù)和坐標(biāo)直接播發(fā)給用戶,在用戶端通過觀測值差分的方式消除各種誤差的影響,實現(xiàn)高精度相對定位。
綜上所述,現(xiàn)有技術(shù)存在的問題是:
廣域差分能夠提供大范圍的精密定位服務(wù),但是由于目前缺乏精密電離層模型支持,通常需要20-30分鐘收斂才能獲得厘米級定位結(jié)果。另外廣域差分使用的參數(shù)化方法各不相同,不同的方法之間互不兼容。目前沒有統(tǒng)一的數(shù)據(jù)播發(fā)格式和計算流程,導(dǎo)致不同廣域差分產(chǎn)品之間不能兼容互操作。另外,由于目前對新的gnss系統(tǒng)信號理解不夠深入,導(dǎo)致大部分廣域差分系統(tǒng)只適用于gps信號。局域差分方法模型簡單,收斂快,實時性好,定位精度高,目前基于載波相位的實時動態(tài)定位(rtk)和基于偽距差分的差分gnss技術(shù)(dgnss)長期以來都是使用最為廣泛的精密定位技術(shù)。但是差分技術(shù)也有其自身的局限,即用戶接收機(也稱為流動站)與參考站之間的距離受限制。用戶接收機與參考站之間的距離變遠(yuǎn)導(dǎo)致兩個接收機之間的電離層和對流層延遲相關(guān)性變?nèi)酰@影響定位精度和收斂時間。另外基線變長也會導(dǎo)致用戶接收機與參考站之間的公共可見衛(wèi)星變少,因而長基線數(shù)據(jù)處理通常用在靜態(tài)數(shù)據(jù)處理中。另外長距離相對定位需要構(gòu)建實時通信鏈路,使用地面實時通信鏈路需要大功率無線電發(fā)射設(shè)備或者復(fù)雜的路由轉(zhuǎn)發(fā)網(wǎng)絡(luò),成本昂貴,通信質(zhì)量難以保證。
技術(shù)實現(xiàn)要素:
針對現(xiàn)有技術(shù)存在的問題,本發(fā)明提供了一種基于低軌衛(wèi)星的星地差分實時精密定位方法。
本發(fā)明是這樣實現(xiàn)的,一種基于低軌衛(wèi)星的星地差分實時精密定位方法,所述基于低軌衛(wèi)星的星地差分實時精密定位方法利用低軌衛(wèi)星作為參考站,與地面或近地空間的接收機進(jìn)行差分定位;所述低軌衛(wèi)星利用地面上注的預(yù)報導(dǎo)航衛(wèi)星星歷和鐘差進(jìn)行地面增強的星上精密自主定軌或者只利用自身接收到的觀測值進(jìn)行自主精密定軌;所述低軌衛(wèi)星以一定的時間間隔向地面播發(fā)星上接收機的觀測值和對應(yīng)時刻的衛(wèi)星軌道信息;地面接收機將接收到低軌衛(wèi)星播發(fā)的差分信息與本機接收的導(dǎo)航衛(wèi)星測距信號組成雙差觀測值,進(jìn)行移動基站的rtk或者dgnss定位。作為參考站的低軌衛(wèi)星明顯高于地面用戶站,保證距離較遠(yuǎn)時有足夠的共視衛(wèi)星。使用星地間通信鏈路解決遠(yuǎn)距離實時通信問題,簡化用戶操作。低軌衛(wèi)星可以全球移動,突破了局域差分對作業(yè)距離的限制。
進(jìn)一步,所述差分定位的模式包括:
(1)通過建立實時數(shù)據(jù)鏈,將leo軌道和觀測數(shù)據(jù)播發(fā)給用戶接收機,進(jìn)行實時差分定位計算;
(2)建立星地實時數(shù)據(jù)連接,分別記錄leo觀測數(shù)據(jù)以及地面接收機觀測數(shù)據(jù),進(jìn)行事后處理精密leo定軌和差分定位計算;
(3)不建立星地實時數(shù)據(jù)連接,通過數(shù)傳或測控鏈路定期地將leo觀測數(shù)據(jù)下傳,在地面進(jìn)行l(wèi)eo精密定軌,再使用leo觀測值,leo定軌結(jié)果和地面接收機觀測數(shù)據(jù)進(jìn)行差分精密定位。
進(jìn)一步,所述觀測值的應(yīng)用模式包括:僅使用偽距觀測值的定位模式、使用相位平滑偽距觀測值的定位模式、同時使用偽距和載波相位觀測值的定位模式。
進(jìn)一步,所述leo自主定軌的方式包括:
(1)不使用地面增強信息,僅使用星上觀測值進(jìn)行完全自主精密定軌;
(2)通過星地間上行通信鏈路,將地面的增強信息上傳至低軌衛(wèi)星,進(jìn)行地基增強的星上精密自主定軌,地面增強信息包括但不限于導(dǎo)航衛(wèi)星預(yù)報軌道、預(yù)報鐘差的信息;
(3)通過數(shù)據(jù)中繼衛(wèi)星或者星間鏈路向低軌衛(wèi)星播發(fā)改正數(shù)信息,進(jìn)行星基增強的實時精密定軌,中繼衛(wèi)星增強信息包括但不限于導(dǎo)航衛(wèi)星預(yù)報軌道、預(yù)報鐘差的信息。
進(jìn)一步,差分精密定位的方法基于單個衛(wèi)星導(dǎo)航系統(tǒng)的觀測數(shù)據(jù),或基于多系統(tǒng)觀測值聯(lián)合定位。
差分精密定位的地面用戶坐標(biāo)參考框架包括全球參考框架或使用自定義坐標(biāo)參考框架或經(jīng)過加密處理的全球參考框架。
進(jìn)一步,所述基于低軌衛(wèi)星的星地差分實時精密定位方法具體包括:
1)接收機讀取當(dāng)前歷元的觀測時間、偽距、載波相位、多普勒、載噪比c/n0、鎖定時間的觀測值信息以及當(dāng)前歷元對應(yīng)的導(dǎo)航電文數(shù)據(jù);所述觀測時間、偽距、載波相位、多普勒、載噪比c/n0、鎖定時間的觀測值信息以及當(dāng)前歷元對應(yīng)的導(dǎo)航電文數(shù)據(jù)均由接收機內(nèi)部pvt算法模塊輸出;其中直接用于定位的觀測量為偽距和載波相位觀測值,觀測模型表示為:
pi=ρ+δorb+c(δts-δtr)+ii+δtrop+εpi
φi=ρ+δorb+c(δts-δtr)-ii+δtrop+λini+εφi
式中:pi和φi表示第i個頻點的偽距和載波相位觀測值(以米為單位);
ρ表示衛(wèi)星和接收機間的幾何距離,表示為
ii和δtrop表示電離層和對流層延遲的影響;
ni表示第i個頻點的整周模糊度;
εpi和εφi表示偽距和載波相位的接收機觀測噪聲;
c表示真空中的光速;
λi表示第i個頻點的載波波長,單位為米;
2)計算信號發(fā)射時刻各個導(dǎo)航衛(wèi)星的軌道和鐘差;信號發(fā)射時刻推算方式包括兩種:(1)通過標(biāo)準(zhǔn)單點定位過程迭代求解各個導(dǎo)航衛(wèi)星的精確信號發(fā)射時刻以及對應(yīng)的衛(wèi)星軌道和鐘差,多次迭代后計算出各衛(wèi)星的軌道鐘差以及地面測站的近似坐標(biāo),(2)利用偽距觀測值近似推算信號發(fā)射時刻,推算原理表達(dá)為:
式中ts,tr分別表示信號發(fā)射和信號接收時刻;
3)根據(jù)單點定位計算得到的近似測站坐標(biāo)和各個導(dǎo)航衛(wèi)星的軌道鐘差的信息,計算測站到各個衛(wèi)星之間的幾何距離ρ,衛(wèi)星高度角等信息;進(jìn)一步根據(jù)經(jīng)驗電離層模型,計算電離層對載波相位和偽距的影響ii,根據(jù)經(jīng)驗對流層模型,計算對流層對gnss信號的影響δtrop,修正之后的觀測值中還剩軌道誤差δorb,接收機鐘差δtr和整周模糊度參數(shù)ni的影響;
4)在非差觀測值層面檢測載波相位的周跳;
5)處理接收到的低軌衛(wèi)星播發(fā)的差分觀測數(shù)據(jù);低軌衛(wèi)星播發(fā)的差分觀測數(shù)據(jù)包含信號接收時間,leo收到的觀測值和接受時刻對應(yīng)的低軌衛(wèi)星軌道,及速度信息;首先需要對差分信息進(jìn)行解碼,獲取到必要的信息;
6)計算leo衛(wèi)星觀測時刻對應(yīng)的導(dǎo)航衛(wèi)星信號發(fā)射時刻軌道和鐘差;
7)根據(jù)leo衛(wèi)星的軌道和導(dǎo)航衛(wèi)星的軌道及鐘差對leo星載接收機觀測值進(jìn)行誤差源改正;
8)非差觀測值層面對leo星載接收機的載波相位觀測值進(jìn)行數(shù)據(jù)探測,周跳檢測;對檢測到發(fā)生周跳的觀測值進(jìn)行標(biāo)記;
9)分別對地面接收機和leo星載接收機的數(shù)據(jù)進(jìn)行預(yù)處理后,開始形成雙差觀測值;首先篩選出地面接收機和leo星載接收機的共視衛(wèi)星,然后每個導(dǎo)航系統(tǒng)選一顆參考星,用于形成星間差分觀測值;
10)形成雙差觀測值;根據(jù)步驟9)確定的共視衛(wèi)星和各系統(tǒng)的參考衛(wèi)星,組成雙差觀測值,星地差分的雙差觀測值表達(dá)為:
式中
11)雙差觀測值層面的周跳檢測;采用三差法周跳檢測;如果對非差觀測值的周跳檢測有足夠的信心,可跳過本步驟,直接進(jìn)行下一步;
12)濾波器時間更新;對于實時坐標(biāo)估計,包括:a)使用序貫最小二乘法和卡爾曼濾波法;卡爾曼濾波模型與觀測值選取有關(guān),對于雙頻觀測數(shù)據(jù),有兩種觀測模型:基于無電離層組合觀測值的模型和基于非組合觀測值的模型;
a)基于無電離層組合觀測值的模型通過雙頻觀測值線性組合消除電離層影響的一階項,但是無電離層組合觀測值的觀測值數(shù)量減半,觀測值噪聲是組合前觀測值的3倍;
b)直接使用非組合觀測值進(jìn)行處理,需要每顆衛(wèi)星估計一個電離層參數(shù);站星差分能削弱電離層的影響,按照非組合模型推導(dǎo);對于非組合模型,濾波器狀態(tài)向量選取為:
x=[δr,δtrop,z,i,n]t;
其中δr為3×1坐標(biāo)增量向量,δtrop,z為天頂方向?qū)α鲗友舆t,每個測站估計一個參數(shù);i為各衛(wèi)星視線方向電離層延遲,每顆可見衛(wèi)星估計一個,n為載波相位雙差模糊度參數(shù),每個信號通道估計一個;對于無電離層組合觀測值的情況,濾波器狀態(tài)向量可忽略電離層延遲參數(shù)i;對于靜態(tài)用戶定位,坐標(biāo)增量δr當(dāng)做隨機常數(shù)估計;對于動態(tài)定位,使用隨機游走過程模擬,或附加外部的運動方程對坐標(biāo)增量參數(shù)進(jìn)行約束;δtrop,z使用隨機游走模型模擬,實現(xiàn)方向電離層延遲i或使用先驗信息初始化,然后使用隨機游走過程模擬;載波相位雙差模糊度在連續(xù)跟蹤時使用隨機常數(shù)模型估計,在發(fā)生周跳時需要參數(shù)重置;在發(fā)生參考星換星時,對所有模糊度參數(shù)做特殊處理;
濾波器時間更新方程表示為:
式中
13)濾波器測量更新;根據(jù)步驟10)的分析,用于雙差定位的偽距和載波相位觀測值線性化表示為:
e(y)=(j,m,θ,λ)(δr,δtrop,z,i,n)t;
式中e(·)為數(shù)學(xué)期望運算符,對f個頻點,s顆可見雙差衛(wèi)星的觀測值,對于使用非組合觀測值的數(shù)學(xué)模型而言,j為(2*s*f)×3的jacobbian矩陣,m為(2*s*f)×1的列向量,每個元素對應(yīng)相應(yīng)衛(wèi)星的對流層投影函數(shù)值;θ為(2*s*f)×s的矩陣;λ為(2*s*f)×(f*s)的矩陣;對于使用無電離層組合的觀測模型,除了忽略電離層參數(shù)對應(yīng)的設(shè)計矩陣θ外,所有觀測值線性組合也造成觀測值個數(shù)比非組合的模型減少一半;線性化后,該非組合觀測方程的設(shè)計矩陣定義為:
a=(j,m,θ,λ);
則卡爾曼濾波的測量更新方程表示為:
kt=pt-a(apt-at+rt)-1,
式中kt為t時刻濾波增益矩陣,rt為幾何觀測值方差協(xié)方差矩陣;
14)驗后殘差檢驗;濾波方程的驗后殘差估值表示為:
根據(jù)驗后殘差估值構(gòu)造χ2檢驗量
15)輸出預(yù)報解;如果驗后殘差檢驗未通過,則表示觀測值可能存在粗差;使用數(shù)據(jù)探測法對粗差進(jìn)行檢驗剔除或放棄使用濾波解,改用預(yù)報解
16)模糊度固定;如果濾波解通過驗后殘差檢驗,則利用載波相位模糊度的整數(shù)特性獲得厘米級精度的模糊度固定解;模糊度固定的方法使用最小二乘降相關(guān)算法lambda或者將無電離層組合分解為寬窄巷模糊度的方式進(jìn)行逐級模糊度固定;
將步驟13)中的狀態(tài)向量按照實數(shù)參數(shù)和整數(shù)參數(shù)兩類合并,則觀測方程簡化表示為e(y)=(b,λ)(b,n)t,相應(yīng)的濾波狀態(tài)量
模糊度參數(shù)的濾波解記做
17)模糊度檢驗;由于錯誤的模糊度固定導(dǎo)致定位結(jié)果產(chǎn)生飛點,定位誤差比浮點解還大,對模糊度固定結(jié)果進(jìn)行檢驗;使用ratio檢驗;
18)輸出浮點解;如果ratio檢驗未通過,則認(rèn)為模糊度固定解不可靠,需要使用濾波解
19)輸出固定解;如果ratio檢驗通過,則使用模糊度固定解作為定位結(jié)果輸出;模糊度固定為整數(shù)后,對應(yīng)的實數(shù)參數(shù)
相應(yīng)的方差協(xié)方差矩陣可表示為:
如果模糊度檢驗通過,則使用
本發(fā)明的另一目的在于提供一種基于低軌衛(wèi)星的星地差分實時精密定位系統(tǒng)。
本發(fā)明的另一目的在于提供一種利用上述基于低軌衛(wèi)星的星地差分實時精密定位系統(tǒng)的遠(yuǎn)洋區(qū)域定位系統(tǒng)。
本發(fā)明的另一目的在于提供一種利用上述基于低軌衛(wèi)星的星地差分實時精密定位系統(tǒng)的偏遠(yuǎn)地區(qū)定位系統(tǒng)
本發(fā)明的另一目的在于提供一種利用上述基于低軌衛(wèi)星的星地差分實時精密定位系統(tǒng)的自然災(zāi)害災(zāi)區(qū)定位系統(tǒng)。
本發(fā)明的優(yōu)點及積極效果為:
本發(fā)明提供的利用低軌衛(wèi)星作為參考站,與地面或近地空間的接收機實現(xiàn)差分定位的方法。低軌衛(wèi)星以一定的時間間隔向地面播發(fā)星上接收機的觀測值和對應(yīng)時刻的衛(wèi)星軌道信息,地面接收機將接收到低軌衛(wèi)星播發(fā)的差分信息與本機接收的導(dǎo)航衛(wèi)星測距信號組成雙差觀測值,實現(xiàn)移動基站的rtk或者dgnss定位。該方法使用低軌衛(wèi)星作為參考站,實現(xiàn)了全球范圍的差分定位服務(wù),特別是對遠(yuǎn)洋區(qū)域,偏遠(yuǎn)地區(qū)以及自然災(zāi)害災(zāi)區(qū)實現(xiàn)快速精密定位有重要意義。
本發(fā)明使用局域差分的方法,但是通過移動的參考站平臺提高了參考站的利用率,使得參考站服務(wù)范圍不局限于周圍幾十公里,而是提供全球范圍的差分服務(wù)。如果使用多顆leo衛(wèi)星組成的星座作為參考基站,則可以提供全球范圍不間斷的站星差分定位服務(wù)
本發(fā)明能夠為受自然條件和經(jīng)濟條件影響導(dǎo)致沒有架設(shè)或者無法架設(shè)地面站的區(qū)域,如遠(yuǎn)洋地區(qū),偏遠(yuǎn)地區(qū),地震災(zāi)區(qū)等區(qū)域提供實時差分定位服務(wù)。
本發(fā)明星地差分定位算法模型比精密單點定位(ppp)簡單,通過雙差的方法可以自行消去衛(wèi)星端和接收機端的鐘差和硬件偏差,無需使用外部信號偏差產(chǎn)品進(jìn)行修正。本發(fā)明利用低軌衛(wèi)星平臺,將地面的固定的參考站用移動的低軌衛(wèi)星平臺取代。星地之間構(gòu)成的雙差觀測值對電離層延遲也具有一定的削弱作用。
本發(fā)明涉及的用戶端僅需要一臺能夠接收leo差分信息的接收機即可實現(xiàn)全球范圍內(nèi)的精密定位,星地差分的方法幫助用戶擺脫了對參考站的依賴。無需自行架設(shè)參考站,無需架設(shè)實時通信鏈路。使用方法與ppp完全一樣,但是卻能提供精度和收斂時間都更好的實時差分精密定位服務(wù)。
附圖說明
圖1是本發(fā)明實施例提供的基于低軌衛(wèi)星的星地差分實時精密定位方法流程圖。
圖2是本發(fā)明實施例提供的低軌衛(wèi)星實現(xiàn)星地差分精密定位的原理圖。
具體實施方式
為了使本發(fā)明的目的、技術(shù)方案及優(yōu)點更加清楚明白,以下結(jié)合實施例,對本發(fā)明進(jìn)行進(jìn)一步詳細(xì)說明。應(yīng)當(dāng)理解,此處所描述的具體實施例僅僅用以解釋本發(fā)明,并不用于限定本發(fā)明。
下面結(jié)合附圖對本發(fā)明的應(yīng)用原理作詳細(xì)描述。
本發(fā)明實施例提供的基于低軌衛(wèi)星的星地差分實時精密定位方法,所述基于低軌衛(wèi)星的星地差分實時精密定位方法利用低軌衛(wèi)星作為參考站,與地面或近地空間的接收機進(jìn)行差分定位;所述低軌衛(wèi)星利用地面上注的預(yù)報導(dǎo)航衛(wèi)星星歷和鐘差進(jìn)行地面增強的星上精密自主定軌或者只利用自身接收到的觀測值進(jìn)行自主精密定軌;所述低軌衛(wèi)星以一定的時間間隔向地面播發(fā)星上接收機的觀測值和對應(yīng)時刻的衛(wèi)星軌道信息;地面接收機將接收到低軌衛(wèi)星播發(fā)的差分信息與本機接收的導(dǎo)航衛(wèi)星測距信號組成雙差觀測值,進(jìn)行移動基站的rtk或者dgnss定位。
所述差分定位包括:
(1)通過建立實時數(shù)據(jù)鏈,將leo軌道和觀測數(shù)據(jù)播發(fā)給用戶接收機,進(jìn)行實時差分定位計算;
(2)建立星地實時數(shù)據(jù)連接,分別記錄leo觀測數(shù)據(jù)以及地面接收機觀測數(shù)據(jù),進(jìn)行事后處理leo精密定軌和差分定位計算;
(3)不建立星地實時數(shù)據(jù)連接,通過數(shù)傳或測控鏈路定期地將leo觀測數(shù)據(jù)下傳,在地面進(jìn)行l(wèi)eo精密定軌,再使用leo觀測值,leo定軌結(jié)果和地面接收機觀測數(shù)據(jù)進(jìn)行差分精密定位。
所述觀測值的應(yīng)用模式包括:僅使用偽距觀測值的定位模式、使用相位平滑偽距觀測值的定位模式、同時使用偽距和載波相位觀測值的定位模式。
所述leo自主定軌的方式包括:
(1)不使用地面增強信息,僅使用星上觀測值進(jìn)行完全自主精密定軌;
(2)通過星地間上行通信鏈路,將地面的增強信息上傳至低軌衛(wèi)星,進(jìn)行地基增強的星上精密自主定軌,地面增強信息包括但不限于導(dǎo)航衛(wèi)星預(yù)報軌道、預(yù)報鐘差的信息;
(3)通過數(shù)據(jù)中繼衛(wèi)星或者星間鏈路向低軌衛(wèi)星播發(fā)改正數(shù)信息,進(jìn)行星基增強的實時精密定軌,中繼衛(wèi)星增強信息包括但不限于導(dǎo)航衛(wèi)星預(yù)報軌道、預(yù)報鐘差的信息。
差分精密定位的方法基于單個衛(wèi)星導(dǎo)航系統(tǒng)的觀測數(shù)據(jù),或基于多系統(tǒng)觀測值聯(lián)合定位。
差分精密定位的地面用戶坐標(biāo)參考框架包括全球參考框架或使用自定義坐標(biāo)參考框架或經(jīng)過加密處理的全球參考框架。如itrf2008,cgcs2000參考框架。
星地差分精密定位方法涉及導(dǎo)航衛(wèi)星星座包括,但不限于現(xiàn)有的全球/區(qū)域衛(wèi)星導(dǎo)航系統(tǒng),如gps,glonass,beidou,galileo,irnss,qzss等。
下面結(jié)合具體實施例對本發(fā)明的應(yīng)用原理作進(jìn)一步描述。
本發(fā)明實施例提供的基于低軌衛(wèi)星的星地差分實時精密定位方法如圖1所示,該算法是一個基于擴展的卡爾曼濾波(ekf)的實時定位算法,可以用于靜態(tài)和動態(tài)的數(shù)據(jù)處理。該算法可以在gnss接收機的內(nèi)部的固件算法中實現(xiàn),也可以實現(xiàn)為后處理算法。為了滿足實時定位的需求,以下描述均以接收機內(nèi)部固件算法為參照。具體的實施步驟如下:
1)接收機讀取當(dāng)前歷元的觀測時間,偽距,載波相位,多普勒,載噪比(c/n0),鎖定時間等觀測值信息,以及當(dāng)前歷元對應(yīng)的導(dǎo)航電文數(shù)據(jù),這些信息均可以由接收機內(nèi)部pvt算法模塊輸出。其中直接用于定位的觀測量為偽距和載波相位觀測值,其觀測模型可以表示為:
pi=ρ+δorb+c(δts-δtr)+ii+δtrop+εpi
φi=ρ+δorb+c(δts-δtr)-ii+δtrop+λini+εφi
式中:pi和φi表示第i個頻點的偽距和載波相位觀測值(以米為單位);
ρ表示衛(wèi)星和接收機間的幾何距離,可表示為
ii和δtrop表示電離層和對流層延遲的影響;
ni表示第i個頻點的整周模糊度;
εpi和εφi表示偽距和載波相位的接收機觀測噪聲;
c表示真空中的光速;
λi表示第i個頻點的載波波長(以米為單位);
2)計算信號發(fā)射時刻各個導(dǎo)航衛(wèi)星的軌道和鐘差??紤]到低軌衛(wèi)星的改正數(shù)播發(fā)時有一定的數(shù)據(jù)處理和信號傳播的延遲,需要對信號發(fā)射時刻進(jìn)行準(zhǔn)確推算??紤]到中高軌衛(wèi)星的衛(wèi)星速度,要求信號發(fā)射時刻推算誤差至少小于1μs,才能保證計算軌道誤差小于1cm。信號發(fā)射時刻推算方式有兩種:(1)通過標(biāo)準(zhǔn)單點定位過程迭代求解各個導(dǎo)航衛(wèi)星的精確信號發(fā)射時刻以及對應(yīng)的衛(wèi)星軌道和鐘差,該方法需要多次迭代計算各衛(wèi)星的軌道鐘差以及地面測站的近似坐標(biāo),但是收斂后信號發(fā)射時刻可推算至100ns以內(nèi);(2)利用偽距觀測值近似推算信號發(fā)射時刻,推算原理可表達(dá)為:
式中ts,tr分別表示信號發(fā)射和信號接收時刻。對比偽距碼p1的觀測模型可知,該式?jīng)]有考慮電離層,對流層和接收機鐘差的影響。其中電離層和對流層對測距的影響在幾十米量級,對信號發(fā)射時刻推算的影響在百納秒量級。接收機鐘差的影響能否忽略取決于接收機鐘差的模型。接收機鐘差的模型有實時修正和毫秒修正兩種模型,實時修正模型接收機鐘差穩(wěn)定在幾納秒至幾百納秒的量級,可以滿足推算需求,毫秒修正模型中,接收機鐘差累計達(dá)到1ms或-1ms就會發(fā)生鐘跳,如果接收機鐘使用毫秒修正模型,則必須使用單點定位方法迭代求解導(dǎo)航衛(wèi)星信號發(fā)射時刻。目前大部分接收機使用或者可配置為鐘差實時修正模型。求解出信號發(fā)射時刻后,即可利用廣播星歷計算各個衛(wèi)星的軌道和鐘差。
3)根據(jù)單點定位計算得到的近似測站坐標(biāo)和各個導(dǎo)航衛(wèi)星的軌道鐘差等信息,計算測站到各個衛(wèi)星之間的幾何距離ρ,衛(wèi)星高度角等信息。進(jìn)一步根據(jù)經(jīng)驗電離層模型,如klobuchar模型,計算電離層對載波相位和偽距的影響ii,根據(jù)經(jīng)驗對流層模型,如hopfield模型計算對流層對gnss信號的影響δtrop,修正之后的觀測值中還剩軌道誤差δorb,接收機鐘差δtr和整周模糊度參數(shù)ni的影響。另外,由于幾何距離計算使用的是近似坐標(biāo),修正后的觀測值中還含有幾何距離殘余量的影響,電離層對流層經(jīng)驗?zāi)P鸵仓荒芨恼徊糠钟绊?,?dǎo)航衛(wèi)星鐘差的殘余影響也可達(dá)納秒量級,這些誤差殘余量也會保留在修正后的觀測值中。
4)在非差觀測值層面檢測載波相位的周跳,常見的方法有根據(jù)載波鎖定環(huán)路(pll)輸出的鎖定時間來判斷環(huán)路是否失鎖。一旦環(huán)路失鎖,載波相位就會被重置,此時輸出的載波相位會產(chǎn)生階躍信號,需要進(jìn)行標(biāo)記。另外一種情況就是信號干擾,信號遮擋,或者多路徑的影響導(dǎo)致載波環(huán)路抖動,這種情況一般不會導(dǎo)致環(huán)路失鎖,但是會導(dǎo)致小的周跳,如半周跳或整周跳。半周跳的情況會影響pvt的比特同步和幀同步,但是做極性判斷通常需要一定的時間,因此發(fā)生半周跳后會有幾秒鐘載波環(huán)路和pvt無法判斷是否發(fā)生了半周跳,這種情況數(shù)據(jù)處理中通常都當(dāng)做周跳來處理。對于小的整周跳,需要使用其他方法判斷,但是一般情況出現(xiàn)的概率比較小。判斷小的整周跳常用的方法有無幾何距離法,如果接收機偽距跟蹤噪聲很小,也可以使用mw組合的方法檢測周跳。
5)接下來處理接收到的低軌衛(wèi)星播發(fā)的差分觀測數(shù)據(jù)。低軌衛(wèi)星播發(fā)的差分觀測數(shù)據(jù)包含信號接收時間,leo收到的觀測值和接受時刻對應(yīng)的低軌衛(wèi)星軌道,及速度信息。首先需要對差分信息進(jìn)行解碼,獲取到必要的信息。
6)計算leo衛(wèi)星觀測時刻對應(yīng)的導(dǎo)航衛(wèi)星信號發(fā)射時刻軌道和鐘差。計算方法和步驟102中論述的相同,不同的是leo衛(wèi)星的坐標(biāo)是已知的,如果使用單點定位計算信號發(fā)射時刻是只需要迭代更新衛(wèi)星坐標(biāo),leo衛(wèi)星的坐標(biāo)可以認(rèn)為是真值。考慮到低軌衛(wèi)星的運動速度約為7.6km/s(以500km軌道高度為例),leo觀測值的鐘面時和實際信號接收時刻之差應(yīng)小于1μs。如果leo星載接收機的鐘差使用實時修正模型,則鐘面時與實際信號接收時刻之差量級小于百納秒量級,對leo軌道的影響小于1cm,可忽略不計。leo向地面播發(fā)的軌道信息應(yīng)經(jīng)過適當(dāng)?shù)膬?nèi)插或外推處理來保證leo觀測值的鐘面時和實際信號接收時刻之差量級小于1μs。
7)根據(jù)leo衛(wèi)星的軌道和導(dǎo)航衛(wèi)星的軌道及鐘差對leo星載接收機觀測值進(jìn)行誤差源改正。該步驟執(zhí)行過程與步驟3)類似,不同的是由于leo屬于空間飛行器,不是地面接收機,其受到的電離層和對流層的影響特性與地面接收機有不同。對流層影響通常指地面至30km高度的近地大氣層,因此一般認(rèn)為對流層對leo星載接收機觀測值的影響可忽略,處理時不對對流層延遲進(jìn)行改正。另外電離層對leo星載接收機的觀測值比地面小,根據(jù)chapman理論,低軌衛(wèi)星星載接收機觀測值受到的電離層影響約為地面的30%-50%。對星載接收機觀測值的電離層改正模型也通常使用地面經(jīng)驗?zāi)P统艘阅骋槐壤蜃幼鳛樾禽d接收機觀測值的電離層改正經(jīng)驗?zāi)P?。由于大氣密度剖面曲線和溫度剖面曲線無法精確測定,chapman模型也只是近似的經(jīng)驗?zāi)P?,不能完全消除電離層對觀測值的影響。
8)非差觀測值層面對leo星載接收機的載波相位觀測值進(jìn)行數(shù)據(jù)探測,周跳檢測。處理方法與步驟104相同,不同的是leo衛(wèi)星的運動速度高于地面接收機,采用一些與幾何距離相關(guān)的方法檢測周跳時,需要特殊考慮閾值的選取。對檢測到發(fā)生周跳的觀測值進(jìn)行標(biāo)記。
9)分別對地面接收機和leo星載接收機的數(shù)據(jù)進(jìn)行預(yù)處理后,開始形成雙差觀測值。雙差觀測值通常是先進(jìn)行測站間差分,再進(jìn)行衛(wèi)星間差分。對于星地基線的情況,首先篩選出地面接收機和leo星載接收機的共視衛(wèi)星,然后每個導(dǎo)航系統(tǒng)選一顆參考星,用于形成星間差分觀測值。參考星的觀測值質(zhì)量直接影響所有雙差觀測值的質(zhì)量,因此需要盡量選取觀測值質(zhì)量較好的衛(wèi)星作為參考星,通常情況使用高度角最高的衛(wèi)星或者信號載噪比最好的衛(wèi)星作為參考衛(wèi)星。另一方面,參考星發(fā)生變化會導(dǎo)致所有雙差模糊度參數(shù)不連續(xù)。在換參考星時模糊度參數(shù)需要特殊處理來避免模糊度參數(shù)頻繁地重置,因此選擇參考星時還會盡量保持參考星不變,來保證模糊度參數(shù)的連續(xù)性。對于多系統(tǒng)聯(lián)合定位的情況,需要每個衛(wèi)星系統(tǒng)分別選擇參考星,由于不同導(dǎo)航系統(tǒng)的信號之間硬件延遲有一定差異,而這個差異使用不同系統(tǒng)的衛(wèi)星間差分無法消除。這個系統(tǒng)間偏差可以通過參數(shù)估計的方法消除,更通用的做法是各個導(dǎo)航系統(tǒng)分別選取參考衛(wèi)星來消除接收機鐘差(含系統(tǒng)間偏差)的影響。
10)形成雙差觀測值;根據(jù)步驟9)確定的共視衛(wèi)星和各系統(tǒng)的參考衛(wèi)星,組成雙差觀測值,星地差分的雙差觀測值可表達(dá)為:
式中
11)雙差觀測值層面的周跳檢測。作為非差層面周跳檢測的補充,形成雙差觀測值后還可以在雙差觀測值層面進(jìn)行周跳檢測,典型的方法有三差法周跳檢測。雙差層面的周跳檢測不是必須的步驟,如果對非差觀測值的周跳檢測有足夠的信心,可跳過本步驟,直接進(jìn)行下一步。
12)濾波器時間更新。完成了雙差觀測值的數(shù)據(jù)處理,即可開始參數(shù)估計的過程。對于實時坐標(biāo)估計,可以使用序貫最小二乘法和卡爾曼濾波法??柭鼮V波法適合靜態(tài)和動態(tài)的數(shù)據(jù)處理,序貫最小二乘法通過引入?yún)?shù)消去過程和先驗信息約束也可以實現(xiàn)與卡爾曼濾波等價的效果。本步驟以卡爾曼濾波求解坐標(biāo)參數(shù)為例,闡述站星間差分定位參數(shù)估計的方法??柭鼮V波模型與觀測值選取有關(guān),對于雙頻觀測數(shù)據(jù),有兩種典型觀測模型:基于無電離層組合觀測值的模型和基于非組合觀測值的模型。基于無電離層組合觀測值的模型可以通過雙頻觀測值線性組合消除電離層影響的一階項,但是無電離層組合觀測值的觀測值數(shù)量減半,觀測值噪聲是組合前觀測值的約3倍。另外一種方法是直接使用非組合觀測值進(jìn)行處理,但是需要每顆衛(wèi)星估計一個電離層參數(shù)??紤]到站星差分可以削弱電離層的影響,以下部分按照非組合模型推導(dǎo)。對于非組合模型,濾波器狀態(tài)向量選取為:
x=[δr,δtrop,z,i,n]t;
其中δr為3×1坐標(biāo)增量向量,δtrop,z為天頂方向?qū)α鲗友舆t,每個測站估計一個參數(shù)。i為各衛(wèi)星視線方向電離層延遲,每顆可見衛(wèi)星估計一個(除了參考衛(wèi)星),n為載波相位雙差模糊度參數(shù),每個信號通道估計一個。對于無電離層組合觀測值的情況,濾波器狀態(tài)向量可忽略電離層延遲參數(shù)i。對于靜態(tài)用戶定位,坐標(biāo)增量δr可以當(dāng)做隨機常數(shù)估計,對于動態(tài)定位,可以使用隨機游走過程模擬,也可以附加外部的運動方程對坐標(biāo)增量參數(shù)進(jìn)行約束。δtrop,z通常使用隨機游走模型模擬,實現(xiàn)方向電離層延遲i可以使用先驗信息初始化,然后使用隨機游走過程模擬。載波相位雙差模糊度在連續(xù)跟蹤時使用隨機常數(shù)模型估計,在發(fā)生周跳時需要參數(shù)重置。在發(fā)生參考星換星時,需要對所有模糊度參數(shù)做特殊處理,避免模糊度參數(shù)重置。濾波器時間更新方程可表示為:
式中
13)濾波器測量更新;根據(jù)步驟10)的分析,用于雙差定位的偽距和載波相位觀測值可以線性化表示為:
e(y)=(j,m,θ,λ)(δr,δtrop,z,i,n)t;
式中e(·)為數(shù)學(xué)期望運算符,對f個頻點,s顆可見雙差衛(wèi)星的觀測值(總可見星個數(shù)減去參考星個數(shù)),對于使用非組合觀測值的數(shù)學(xué)模型而言,j為(2*s*f)×3的jacobbian矩陣,m為(2*s*f)×1的列向量,每個元素對應(yīng)相應(yīng)衛(wèi)星的對流層投影函數(shù)值。θ為(2*s*f)×s的矩陣。λ為(2*s*f)×(f*s)的矩陣。對于使用無電離層組合的觀測模型,除了可忽略電離層參數(shù)對應(yīng)的設(shè)計矩陣θ外,所有觀測值線性組合也造成觀測值個數(shù)比非組合的模型減少一半。線性化后,非組合模型觀測方程的設(shè)計矩陣定義為:
a=(j,m,θ,λ);
則卡爾曼濾波的測量更新方程可表示為:
kt=pt-a(apt-at+rt)-1,
pt+=(i-kta)pt-,
式中kt為t時刻濾波增益矩陣,rt為幾何觀測值方差協(xié)方差矩陣。
14)驗后殘差檢驗;濾波方程的驗后殘差估值可表示為:
根據(jù)驗后殘差估值可構(gòu)造χ2檢驗量
15)輸出預(yù)報解;如果驗后殘差檢驗未通過,則表示觀測值可能存在粗差??墒褂脭?shù)據(jù)探測法對粗差進(jìn)行檢驗剔除,也可以放棄使用濾波解,改用預(yù)報解
16)模糊度固定;如果濾波解通過驗后殘差檢驗,則可嘗試?yán)幂d波相位模糊度的整數(shù)特性獲得厘米級精度的模糊度固定解。由于星地差分基線組成雙差觀測值,其定位過程等價于長基線相對定位。模糊度固定的方法可以使用最小二乘降相關(guān)算法(lambda)或者將無電離層組合分解為寬窄巷模糊度的方式進(jìn)行逐級模糊度固定。將步驟13)中的狀態(tài)向量按照實數(shù)參數(shù)和整數(shù)參數(shù)兩類合并,則觀測方程可以簡化表示為e(y)=(b,λ)(b,n)t,相應(yīng)的濾波狀態(tài)量
模糊度參數(shù)的濾波解記做
17)模糊度檢驗;由于錯誤的模糊度固定可能導(dǎo)致定位結(jié)果產(chǎn)生飛點,定位誤差比浮點解還大,因此要對模糊度固定結(jié)果進(jìn)行檢驗。通常使用ratio檢驗。ratio檢驗定義為:
式中
18)輸出浮點解;如果ratio檢驗未通過,則認(rèn)為模糊度固定解不可靠,需要使用濾波解
19)輸出固定解;如果ratio檢驗通過,則使用模糊度固定解作為定位結(jié)果輸出。模糊度固定為整數(shù)后,對應(yīng)的實數(shù)參數(shù)
相應(yīng)的方差協(xié)方差矩陣可表示為:
如果模糊度檢驗通過,則使用
上述步驟可實現(xiàn)地面接收機與低軌衛(wèi)星之間的差分實時精密定位。
本發(fā)明實施例提供的低軌衛(wèi)星實現(xiàn)星地差分精密定位的原理如圖2所示。圖中g(shù)eo/igso/meo指現(xiàn)有的全球或區(qū)域衛(wèi)星導(dǎo)航系統(tǒng)衛(wèi)星,包括但不限于gps、glonass、beidou、galileo、irnss及qzss導(dǎo)航衛(wèi)星系統(tǒng)。leo是指搭載了雙頻gnss信號接收機,具有高精度自主定軌能力,并具有向地面播發(fā)信息能力的低軌衛(wèi)星。地面參考站是指布設(shè)在地球表面的連續(xù)運行參考站,如國際gnss服務(wù)的監(jiān)測網(wǎng)(igs),多gnss試驗網(wǎng)(mgex)等。地面分析中心是指能夠通過網(wǎng)絡(luò)獲取各個地面參考站的觀測數(shù)據(jù),能夠?qū)崟r計算或者預(yù)報導(dǎo)航衛(wèi)星的軌道和鐘差,并且能夠?qū)㈩A(yù)報的導(dǎo)航衛(wèi)星軌道和鐘差上注到低軌衛(wèi)星的地面中心數(shù)據(jù)處理中心。用戶是指具有g(shù)nss信號接收功能,并且能夠接收,解析低軌衛(wèi)星播發(fā)的差分信息,通過差分定位的方式確定自身精密坐標(biāo)的接收設(shè)備。如圖2所示,圖中g(shù)eo/igso/meo導(dǎo)航衛(wèi)星持續(xù)的向地面和低軌衛(wèi)星播發(fā)導(dǎo)航信號,導(dǎo)航信號被低軌衛(wèi)星,用戶和地面參考站接收。地面參考站將接收到的觀測值發(fā)送至地面分析中心,地面分析中心通過綜合計算各個地面站的觀測數(shù)據(jù),確定實時和預(yù)報的精密geo/igso/meo軌道和鐘差信息。考慮到低軌衛(wèi)星測控窗口長度有限,通常需要預(yù)報數(shù)小時的軌道和鐘差,然后將預(yù)報軌道和鐘差信息進(jìn)行數(shù)據(jù)壓縮,在測控窗口內(nèi),將壓縮后的數(shù)據(jù)上傳至低軌衛(wèi)星。低軌衛(wèi)星使用星載接收機的觀測數(shù)據(jù)和接收到的地面增強信息進(jìn)行實時精密軌道確定。然后將自身的軌道和同一時刻對應(yīng)的觀測數(shù)據(jù)進(jìn)行編碼,播發(fā)到地面。地面用戶使用收到的geo/igso/meo衛(wèi)星的觀測值,和低軌衛(wèi)星播發(fā)的差分信息,形成星地差分基線,完成雙差定位。雖然星地基線距離較長,通常有數(shù)百至上千公里,但是形成的雙差觀測值仍然能夠有效地消除接收機鐘差,衛(wèi)星鐘差,削弱電離層影響,有效地簡化精密定位的復(fù)雜度,提高定位效率和精度。另外,本文論述的星地差分系統(tǒng)也可以在自主模式運行,即沒有地面分析中心支持的模式。沒有地面參考站和和地面分析中心支持,星上自主定軌也可以達(dá)到0.5m左右的結(jié)果。直接使用自主定軌的結(jié)果進(jìn)行星地差分,仍然可獲得優(yōu)于1米的實時定位精度。地面分析中心與低軌衛(wèi)星之間的直接通信鏈路也可以通過數(shù)據(jù)中繼衛(wèi)星或者星間鏈路等方式間接的建立。
本發(fā)明能夠提高地面參考站的利用率,將差分定位服務(wù)范圍擴展至全球。
本發(fā)明能夠提供全球范圍的差分定位服務(wù),特別是在遠(yuǎn)洋地區(qū),偏遠(yuǎn)地區(qū),地震災(zāi)區(qū)等沒有條件架設(shè)地面參考站的區(qū)域。
本發(fā)明星地差分定位算法模型比精密單點定位(ppp)簡單,通過雙差的方法可以自行消去衛(wèi)星端和接收機端的鐘差和硬件偏差,因而無需使用外部信號偏差產(chǎn)品進(jìn)行修正。基于星地差分定位系統(tǒng)的定位收斂時間和定位精度方面優(yōu)于精密單點定位。
本發(fā)明能夠提供面向?qū)崟r用戶的高精度定位服務(wù),精密單點定位技術(shù)要依賴多種外部精密數(shù)據(jù)產(chǎn)品對觀測數(shù)據(jù)進(jìn)行改正,部分產(chǎn)品具有一定的延遲,這導(dǎo)致精密單點定位技術(shù)更多地應(yīng)用在后處理和實時性不高的場景。本發(fā)明的方法適合實時高精度定位應(yīng)用。
以上所述僅為本發(fā)明的較佳實施例而已,并不用以限制本發(fā)明,凡在本發(fā)明的精神和原則之內(nèi)所作的任何修改、等同替換和改進(jìn)等,均應(yīng)包含在本發(fā)明的保護(hù)范圍之內(nèi)。