一種rt-ppp星歷改正數(shù)據(jù)優(yōu)選算法
【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及衛(wèi)星導(dǎo)航技術(shù)領(lǐng)域，具體涉及到一種RT-ppp星歷改正數(shù)據(jù)優(yōu)選算法。
【背景技術(shù)】
[0002] 在RT-ppp定位過(guò)程中，IGU星歷特別是其衛(wèi)星鐘差精度往往成為影響精度和收斂 時(shí)間的關(guān)鍵性因素。鑒于目前尚不存在比較好的星歷質(zhì)量評(píng)估定權(quán)方法，提出此解決方案。 IGU的預(yù)報(bào)星歷和觀測(cè)星歷時(shí)間間隔僅僅為3小時(shí)，但是其軌道精度尤其是衛(wèi)星鐘種差精度 卻相差甚遠(yuǎn)，不同型號(hào)GPS衛(wèi)星的質(zhì)量、形狀及光壓模型都存在差異，導(dǎo)致其預(yù)報(bào)軌道精度 不同，從而造成IGU的預(yù)報(bào)軌道短時(shí)間內(nèi)的質(zhì)量具有非常強(qiáng)的相關(guān)性。

【發(fā)明內(nèi)容】

[0003] 針對(duì)現(xiàn)有技術(shù)中RT-PPP定位存在上述的不足，本發(fā)明提供了一種RT-PPP星歷改正 數(shù)據(jù)優(yōu)選算法，能夠?qū)崟r(shí)為ppp解算提供經(jīng)過(guò)自動(dòng)軌道質(zhì)量控制的星歷，提高了實(shí)時(shí)ppp的 定位精度。
[0004] 本發(fā)明的技術(shù)方案為，一種RT-PPP星歷改正數(shù)據(jù)優(yōu)選算法，包括以下步驟：
[0005] 連續(xù)采集IGU精密星歷的數(shù)據(jù)；
[0006] 根據(jù)采集的數(shù)據(jù)對(duì)IGU觀測(cè)星歷和其對(duì)應(yīng)此時(shí)刻的IGU預(yù)報(bào)星歷進(jìn)行求差；
[0007]根據(jù)IGU觀測(cè)星歷和其對(duì)應(yīng)此時(shí)刻的IGU預(yù)報(bào)星歷的差值進(jìn)行第一次定權(quán)并得到 權(quán)陣；
[0008] 根據(jù)上述權(quán)陣通過(guò)軌道松弛法進(jìn)行參數(shù)約束并進(jìn)行偽距定位解算。
[0009] 上述的優(yōu)選算法，其中，所述根據(jù)采集的數(shù)據(jù)對(duì)IGU觀測(cè)星歷和其對(duì)應(yīng)此時(shí)刻的 IGU預(yù)報(bào)星歷進(jìn)行求差的步驟具體為，將IGU觀測(cè)軌道與對(duì)應(yīng)的時(shí)刻之前儲(chǔ)存下來(lái)的預(yù)報(bào)軌 道進(jìn)行比較，并且求差，得到軌道參數(shù)差值 Vl，其中包括歸算到單位為米的衛(wèi)星鐘種差。
[0010] 上述的優(yōu)選算法，其中，所述根據(jù)ICTJ觀測(cè)星歷和其對(duì)應(yīng)此時(shí)刻的IGU預(yù)報(bào)星歷的 差值進(jìn)行第一次定權(quán)并得到權(quán)陣的步驟具體為，取定標(biāo)準(zhǔn)差〇〇，則初始權(quán)陣P = σο/V;偽距觀 測(cè)方程為:L=A*X+1，Ρ，其中L為偽距觀測(cè)值向量，Α為設(shè)計(jì)矩陣，X是坐標(biāo)位置，1為L(zhǎng)-L〇。
[0011] 上述的優(yōu)選算法，其中，所述根據(jù)上述權(quán)陣通過(guò)軌道松弛法進(jìn)行參數(shù)約束并進(jìn)行 偽距定位解算的步驟包括：
[0012] 將權(quán)陣應(yīng)用于實(shí)時(shí)定位的最小二乘定位方程，并最終求得方差陣；
[0013] 將方差陣通過(guò)軌道松弛法進(jìn)行參數(shù)約束并對(duì)有約束的軌道進(jìn)行偽距定位解算。
[0014] 上述的優(yōu)選算法，其中，所述根據(jù)上述權(quán)陣通過(guò)軌道松弛法進(jìn)行參數(shù)約束并進(jìn)行 偽距定位解算的步驟中，軌道松弛約束條件：
其中〇sxl后為驗(yàn)標(biāo)準(zhǔn)差，通過(guò)解偽距 觀測(cè)方程得到，軌道參數(shù)Sxi和k對(duì)應(yīng)t分布，取定k = 2.5達(dá)到98.8%的置信度。
[0015] 上述的優(yōu)選算法，其中，若所述軌道松弛約束條件滿足：
取定k = 2.5達(dá)到 98.8 %的置信度，則更新標(biāo)準(zhǔn)差為：
[0016]上述的優(yōu)選算法，其中，將求得的新標(biāo)準(zhǔn)差歸算到權(quán)陣中重新進(jìn)行偽距定位最小 二乘計(jì)算直到軌道參數(shù)均滿足軌道松弛約束條件為止，其中所述的軌道松弛約束條件為
，其中σδχ1后為驗(yàn)標(biāo)準(zhǔn)差，通過(guò)解偽距觀測(cè)方程得到，軌道參數(shù)δ χ1和k對(duì)應(yīng)t分布，取 定k = 2.5達(dá)到98.8 %的置信度。
[0017]本發(fā)明提供的一種RT-PPP星歷改正數(shù)據(jù)優(yōu)選算法，包括以下步驟:連續(xù)采集I⑶精 密星歷的數(shù)據(jù);根據(jù)采集的數(shù)據(jù)對(duì)IGU觀測(cè)星歷和其對(duì)應(yīng)此時(shí)刻的IGU預(yù)報(bào)星歷進(jìn)行求差； 根據(jù)IGU觀測(cè)星歷和其對(duì)應(yīng)此時(shí)刻的IGU預(yù)報(bào)星歷的差值進(jìn)行第一次定權(quán)并得到權(quán)陣;本發(fā) 明適用于RT-PPP定位，能夠?qū)崟r(shí)為PPP解算提供經(jīng)過(guò)自動(dòng)軌道質(zhì)量控制的星歷，提高了實(shí)時(shí) PPP的定位精度，以及利用短時(shí)間內(nèi)預(yù)報(bào)星歷質(zhì)量的強(qiáng)相關(guān)性進(jìn)行星歷權(quán)重判定，為軌道松 弛提供權(quán)重初始值，并且無(wú)需為接收機(jī)增加新硬件或者服務(wù)器，僅僅需要在軟件層面上可 以更改現(xiàn)有的RT-PPP接收機(jī)的自動(dòng)化軌道質(zhì)量控制。
【附圖說(shuō)明】
[0018] 通過(guò)閱讀參照以下附圖對(duì)非限制性實(shí)施例所作的詳細(xì)描述，本發(fā)明及其特征、外 形和優(yōu)點(diǎn)將會(huì)變得更明顯。在全部附圖中相同的標(biāo)記指示相同的部分。并未刻意按照比例 繪制附圖，重點(diǎn)在于示出本發(fā)明的主旨。
[0019] 圖1為本發(fā)明提供的一種RT-PPP星歷改正數(shù)據(jù)優(yōu)選算法的示意圖。
【具體實(shí)施方式】
[0020] 在下文的描述中，給出了大量具體的細(xì)節(jié)以便提供對(duì)本發(fā)明更為徹底的理解。然 而，對(duì)于本領(lǐng)域技術(shù)人員而言顯而易見(jiàn)的是，本發(fā)明可以無(wú)需一個(gè)或多個(gè)這些細(xì)節(jié)而得以 實(shí)施。在其他的例子中，為了避免與本發(fā)明發(fā)生混淆，對(duì)于本領(lǐng)域公知的一些技術(shù)特征未進(jìn) 行描述。
[0021] 為了徹底理解本發(fā)明，將在下列的描述中提出詳細(xì)的步驟以及詳細(xì)的結(jié)構(gòu)，以便 闡釋本發(fā)明的技術(shù)方案。本發(fā)明的較佳實(shí)施例詳細(xì)描述如下，然而除了這些詳細(xì)描述外，本 發(fā)明還可以具有其他實(shí)施方式。
[0022] 參照?qǐng)D1所示，本發(fā)明提供了一種RT-PPP星歷改正數(shù)據(jù)優(yōu)選算法，包括以下步驟： [0023]步驟S1:連續(xù)采集IGU精密星歷的數(shù)據(jù)，不斷采集IGU超快精密星歷，并將數(shù)據(jù)用于 后續(xù)的計(jì)算。
[0024]步驟S2:根據(jù)采集的數(shù)據(jù)對(duì)IGU觀測(cè)星歷和其對(duì)應(yīng)此時(shí)刻的IGU預(yù)報(bào)星歷進(jìn)行求 差，將IGU觀測(cè)軌道與對(duì)應(yīng)的時(shí)刻之前儲(chǔ)存下來(lái)的預(yù)報(bào)軌道進(jìn)行比較，并且求差，得到軌道 參數(shù)差值^，其中包括歸算到單位為米的衛(wèi)星鐘種差。
[0025]步驟S3:根據(jù)IGU觀測(cè)星歷和其對(duì)應(yīng)此時(shí)刻的IGU預(yù)報(bào)星歷的差值進(jìn)行第一次定權(quán) 并得到權(quán)陣，具體為，取定標(biāo)準(zhǔn)差〇〇,則初始權(quán)陣Ρ = σ〇/ν;偽距觀測(cè)方程為:L = A*X+1，P，其 中L為偽距觀測(cè)值向量，A為設(shè)計(jì)矩陣，X是坐標(biāo)位置，1為L(zhǎng)-Lo,也就是說(shuō)對(duì)最新的ICTJ觀測(cè)星 歷和其對(duì)應(yīng)此時(shí)刻的IGU預(yù)報(bào)星歷進(jìn)行求差，并根據(jù)此差值進(jìn)行第一次定權(quán)，得到權(quán)陣Po。
[0026] 步驟S4:根據(jù)上述權(quán)陣通過(guò)軌道松弛法進(jìn)行參數(shù)約束并進(jìn)行偽距定位解算