北斗衛(wèi)星偽距多路徑系統(tǒng)偏差的修正方法及其對(duì)單頻PPP的影響分析

汪 捷，何錫揚(yáng)

1. 海軍工程大學(xué)，湖北 武漢 430033； 2. 千尋位置網(wǎng)絡(luò)有限公司，上海 200438

?

北斗衛(wèi)星偽距多路徑系統(tǒng)偏差的修正方法及其對(duì)單頻PPP的影響分析

汪 捷1，何錫揚(yáng)2

1. 海軍工程大學(xué)，湖北 武漢 430033； 2. 千尋位置網(wǎng)絡(luò)有限公司，上海 200438

北斗偽距觀測(cè)值存在特有的多路徑系統(tǒng)性偏差，偏差的數(shù)量級(jí)達(dá)到幾個(gè)分米到米。該系統(tǒng)偏差可分為兩類：一類是IGSO/MEO衛(wèi)星隨高度角變化的偽距系統(tǒng)性偏差；另一類是GEO衛(wèi)星(高度角僅微小變化)明顯的偽距系統(tǒng)性偏差。系統(tǒng)性的偽距偏差導(dǎo)致GEO衛(wèi)星MP序列的標(biāo)準(zhǔn)偏差較大，本文針對(duì)GEO衛(wèi)星偽距偏差問題提出了一種基于卡爾曼濾波的修正方法，修正后的GEO衛(wèi)星MP序列的標(biāo)準(zhǔn)偏差下降了10%～16%。基于偽距相位組合的單頻PPP技術(shù)的偽距權(quán)重較大，會(huì)受到北斗偽距偏差的影響，分析表明該系統(tǒng)性偏差將導(dǎo)致單頻PPP定位結(jié)果高程方向產(chǎn)生約1 m的偏差。對(duì)GEO偽距偏差采用提出的卡爾曼濾波修正方法進(jìn)行修正，并應(yīng)用Wanninger和Beer的高度角模型消除IGSO/MEO觀測(cè)值偽距偏差，本文對(duì)修正后的單頻精密單點(diǎn)定位精度進(jìn)行了分析。4個(gè)multi-GNSS experiment(MGEX)站10 d觀測(cè)數(shù)據(jù)的分析結(jié)果表明：僅改正和衛(wèi)星多路徑誤差，高程方向定位結(jié)果精度可改善65%左右；采用本文方法對(duì)GEO衛(wèi)星的多路徑修正后，該方向定位結(jié)果精度改善比例將進(jìn)一步提高至75%左右。

北斗；單頻精密單點(diǎn)定位；濾波；多路徑；偽距；精度分析

我國(guó)于2012年底建成了具備導(dǎo)航定位授時(shí)能力的北斗衛(wèi)星導(dǎo)航定位系統(tǒng)(BeiDou navigation satellite system，BDS)。截至2015年12月，北斗導(dǎo)航衛(wèi)星系統(tǒng)包括5顆Geostationary Orbit(GEO)衛(wèi)星，5顆Inclined Geosynchronous Orbit(IGSO)衛(wèi)星，4顆Medium Earth Orbit(MEO)衛(wèi)星[1]。許多學(xué)者對(duì)北斗系統(tǒng)在標(biāo)準(zhǔn)單點(diǎn)定位[2]、精密單點(diǎn)定位[3-4]、雙頻([5-6])以及多頻相對(duì)定位[7]等多種模式下的定位性能進(jìn)行了分析評(píng)價(jià)，驗(yàn)證了北斗系統(tǒng)的可用性。

近年來(lái)一些研究表明[8-12]，BDS偽距觀測(cè)值和相位觀測(cè)值之間存在系統(tǒng)性差異，該差異在高高度角情況下達(dá)到0.4 m，在低高度角情況下甚至超過1 m。該系統(tǒng)性偏差可分為兩類：一類是IGSO/MEO衛(wèi)星隨高度角變化的偽距系統(tǒng)性偏差；另一類是GEO衛(wèi)星(高度角僅微小變化)明顯的偽距系統(tǒng)性偏差。該系統(tǒng)性偏差對(duì)單點(diǎn)定位、寬巷模糊度固定、單頻PPP定位等都有顯著的影響。針對(duì)此問題，文獻(xiàn)[13]提出了一種基于高度角的偽距改正模型，該模型可應(yīng)用于北斗IGSO衛(wèi)星和MEO衛(wèi)星的偽距觀測(cè)值改正，其測(cè)試結(jié)果表明，該改正模型可顯著改善單頻PPP的定位精度。

對(duì)于GEO衛(wèi)星的系統(tǒng)性偏差[14-15]，因單個(gè)測(cè)站的GEO衛(wèi)星高度角幾乎不變，文獻(xiàn)[13]未對(duì)此建立模型。文獻(xiàn)[16]基于多個(gè)測(cè)站的觀測(cè)值建立了GEO衛(wèi)星的高度角改正模型。對(duì)于單個(gè)靜態(tài)站或低動(dòng)態(tài)觀測(cè)，由于GEO衛(wèi)星高度角幾乎不變，其修正值是恒定的。

文獻(xiàn)[13]及文獻(xiàn)[16]所給出的改正模型，更多地是對(duì)北斗衛(wèi)星端的多路徑系統(tǒng)誤差(星內(nèi)多徑)進(jìn)行建模。本文的分析表明，雖然GEO衛(wèi)星高度角幾乎不變，但是由于受到測(cè)站端多路徑影響較為嚴(yán)重，其MP序列存在若干分米量級(jí)的抖動(dòng)，已有改正模型難以準(zhǔn)確描述衛(wèi)星偽距偏差變化。因此，本文針對(duì)GEO衛(wèi)星偽距偏差問題，提出了一種基于卡爾曼濾波的修正方法。濾波模型兼顧GEO衛(wèi)星觀測(cè)值衛(wèi)星端和接收機(jī)端多路徑的影響，將兩部分誤差統(tǒng)一進(jìn)行處理。

單頻接收機(jī)具有低成本優(yōu)勢(shì)，以單頻接收機(jī)實(shí)現(xiàn)用戶高精度單機(jī)定位已成為近年來(lái)的研究熱點(diǎn)[17]。單頻PPP是其中一種重要的技術(shù)方法，目前已有不少學(xué)者開展了GPS單頻PPP方面的研究工作，取得了豐富的成果[18-20]。基于偽距相位組合的單頻PPP(precise point positioning，PPP)技術(shù)的偽距權(quán)重較大，會(huì)受到上文提到的北斗偽距偏差的影響。本文分析了該偏差對(duì)單頻PPP水平和高程方向精度的影響，并對(duì)比了改正偽距偏差后的效果。

1 BDS偽距多路徑系統(tǒng)誤差

Multipath(MP)組合觀測(cè)值常用于評(píng)估偽距觀測(cè)值的多路徑特征，其計(jì)算公式如下[21-22]

(1)

(2)

MP組合消除了觀測(cè)值的電離層延遲和幾何相關(guān)項(xiàng)(如對(duì)流層延遲、幾何距離、鐘差等)，剩余項(xiàng)包括相位觀測(cè)值模糊度、信號(hào)硬件延遲差、多路徑以及觀測(cè)噪聲。通常情況下，信號(hào)硬件延遲在一定時(shí)間內(nèi)保持穩(wěn)定，而相位觀測(cè)值模糊度在未發(fā)生周跳和失鎖的情況下同樣保持不變，因此B值中包含硬件延遲和模糊度。對(duì)MP組合進(jìn)行長(zhǎng)時(shí)間平均，可確定B值的大小。MP原始組合減去B值后僅剩下觀測(cè)噪聲和多路徑誤差變化部分。由于觀測(cè)噪聲通常表現(xiàn)為高斯白噪聲特性，因此通過MP組合可以分析偽距觀測(cè)值的多路徑系統(tǒng)誤差特性。

圖1給出了基于GMSD測(cè)站數(shù)據(jù)計(jì)算的北斗3種不同類型衛(wèi)星B1頻率偽距觀測(cè)值多路徑系統(tǒng)誤差的典型示例，數(shù)據(jù)觀測(cè)時(shí)間為2016年8月1日—2016年8月2日(DOY213—DOY214)。GEO衛(wèi)星給出了連續(xù)2d的MP序列，IGSO衛(wèi)星和MEO衛(wèi)星則給出了一個(gè)觀測(cè)弧段的MP序列。由圖1可知，3種類型衛(wèi)星觀測(cè)值的多路徑誤差均存在較為明顯的系統(tǒng)性偏差，波動(dòng)量級(jí)從幾個(gè)分米到1m左右。觀測(cè)值的多路徑與信號(hào)高度角存在較強(qiáng)的相關(guān)性，GEO衛(wèi)星的多路徑同時(shí)還呈現(xiàn)出約1d為周期的變化。

IGSO衛(wèi)星和MEO衛(wèi)星偽距觀測(cè)值的系統(tǒng)多路徑誤差可采用Wanninger和Beer的模型進(jìn)行改正。圖2給出了IGSO衛(wèi)星和MEO衛(wèi)星修正多路徑誤差后的MP序列。從圖中可以看到，經(jīng)過模型修正后，IGSO衛(wèi)星和MEO衛(wèi)星MP序列的系統(tǒng)性誤差被明顯削弱或消除。

圖2 修正多路徑誤差后的MP序列Fig.2 MP series after multipath correction

2 GEO衛(wèi)星觀測(cè)值多路徑濾波方法

本節(jié)主要研究如何對(duì)GEO偽距觀測(cè)值多路徑進(jìn)行修正。有關(guān)多路徑的研究方法有多種，如小波變換法[23-24]、恒星日濾波法[25-26]，本文擬采用卡爾曼濾波方法對(duì)多路徑偏差進(jìn)行分析和修正。

2.1 特征分析

MP序列反映了多路徑的量級(jí)，但同時(shí)受到噪聲的影響。本文首先采用卡爾曼濾波對(duì)MP序列進(jìn)行擬合，以去除噪聲的影響，進(jìn)而得到多路徑修正量，對(duì)偽距觀測(cè)值進(jìn)行修正。圖3給出了基于GMSD站2015年8月1日(年積日213)數(shù)據(jù)獲得的GEO C01衛(wèi)星MP組合的原始序列、歷元一次差序列和歷元二次差序列。從圖中可以看出，經(jīng)過一次差或者二次差后的MP序列，趨近于白噪聲。因此，采用常加速模型可表達(dá)出GEO多路徑系統(tǒng)偏差。

圖3 C01觀測(cè)值MP原始序列、一次差序列、二次差序列Fig.3 Original, first-order differenced and second-order differenced MP series for C01 satellite mesurements

2.2 濾波方法

常加速模型的狀態(tài)方程為

(3)

式中，Δt為時(shí)間差。假定加速度是一個(gè)具有常量譜密度的白噪聲過程，系統(tǒng)噪聲(過程噪聲)的協(xié)方差矩陣為

(4)

對(duì)于MP序列，其濾波狀態(tài)為MP系統(tǒng)偏差量、及其隨時(shí)間變化的速度和加速度。采用該濾波模型，可提取出MP序列的多路徑系統(tǒng)偏差，進(jìn)而對(duì)偽距觀測(cè)值進(jìn)行修正。由于該系統(tǒng)偏差存在約一天周期的變化，因此可基于前一天數(shù)據(jù)獲取的修正對(duì)當(dāng)天的偽距觀測(cè)值進(jìn)行修正。

2.3 修正效果

圖4(a)給出了GMSD站C01衛(wèi)星的MP序列濾波示例，數(shù)據(jù)觀測(cè)時(shí)間為2015年8月1日(年積日213)；圖4(b)給出了經(jīng)過修正后的MP序列?？梢钥闯?，本文的濾波方法可以較好地濾波出多路徑系統(tǒng)偏差項(xiàng)，經(jīng)過修正后的MP序列非常接近于白噪聲特性。圖5給出了DOY 213和DOY 214相鄰兩天多路徑偏差項(xiàng)對(duì)比，相鄰兩天的修正量非常類似，進(jìn)一步證實(shí)可基于前一天數(shù)據(jù)獲取的修正量對(duì)當(dāng)天的偽距觀測(cè)值進(jìn)行修正。根據(jù)文獻(xiàn)[27]的研究，GEO衛(wèi)星多路徑周期約為86 160 s，修正時(shí)依據(jù)該周期進(jìn)行。圖6給出了GMSD測(cè)站2015年8月2日4顆GEO衛(wèi)星觀測(cè)值MP序列修正前后的標(biāo)準(zhǔn)偏差統(tǒng)計(jì)值，修正值基于2015年8月1日數(shù)據(jù)計(jì)算得到。經(jīng)過修正后，MP序列的標(biāo)準(zhǔn)偏差下降了10%～16%。

3 修正偽距多路徑系統(tǒng)誤差改善單頻PPP定位精度

本節(jié)分析偽距多路徑偏差對(duì)單頻PPP的影響及修正偽距多路徑偏差后的改善效果。

3.1 單頻PPP模型

單頻PPP的一種形式為偽距相位組合模型[18]，該模型利用偽距和相位一階電離層延遲符號(hào)相反、大小相同的特性，消除了一階電離層延遲的影響；其缺點(diǎn)在于偽距在觀測(cè)方程中占了較大的權(quán)重。當(dāng)偽距噪聲或多路徑較大時(shí)，如上文的BDS偽距多路徑系統(tǒng)偏差，該模型精度將受到影響。

3.2 處理策略

采用GFZ提供的精密星歷和精密衛(wèi)星鐘差產(chǎn)品、以及MGEX的北斗DCB產(chǎn)品進(jìn)行單頻PPP解算。估計(jì)參數(shù)包括測(cè)站坐標(biāo)、接收機(jī)鐘差、天頂對(duì)流層延遲參數(shù)、相位觀測(cè)值模糊度。采用卡爾曼濾波方式進(jìn)行數(shù)據(jù)處理。接收機(jī)鐘差參數(shù)采用白噪聲模型，對(duì)流層參數(shù)初值采用Saastamonien模型確定[28]，估計(jì)參數(shù)采用隨機(jī)游走模型，投影函數(shù)為NMF模型[29]，位置參數(shù)當(dāng)作常數(shù)估計(jì)(靜態(tài)模型)，模糊度參數(shù)在未失鎖和未發(fā)生周跳情況下認(rèn)為不變，發(fā)生周跳或失鎖時(shí)間過長(zhǎng)則重置。偽距與相位驗(yàn)前標(biāo)準(zhǔn)差之比為1∶100。截止高度角取為15°。

采用第2節(jié)的濾波方法消除GEO衛(wèi)星偽距觀測(cè)值的多路徑誤差：基于B1頻率偽距觀測(cè)值和B1和B2頻率的相位觀測(cè)值計(jì)算MP序列，由前一周期提取的多路徑誤差模型，對(duì)當(dāng)前觀測(cè)值進(jìn)行多路徑誤差改正。IGSO和MEO衛(wèi)星的系統(tǒng)性偏差采用Wanninger和Beer的高度角模型進(jìn)行修正。

3.3 數(shù)據(jù)處理與分析

為了測(cè)試改正多路徑系統(tǒng)誤差前后北斗單頻PPP的定位精度情況，選擇MGEX網(wǎng)亞太地區(qū)的4個(gè)測(cè)站進(jìn)行分析，所用數(shù)據(jù)的具體信息見表1。由于需要基于前一天數(shù)據(jù)計(jì)算修正值，同時(shí)獲取了所有測(cè)站2015(212)的觀測(cè)數(shù)據(jù)用于計(jì)算2015(213)的修正值。

表1 測(cè)站信息

圖7給出了GMSD站2015(213)的單頻PPP定位結(jié)果在N、E、U 3個(gè)方向的定位誤差。由圖7可以看出，在收斂之后，修正IGSO和MEO多路徑系統(tǒng)偏差的定位結(jié)果在E方向略差于未修正多路徑系統(tǒng)偏差的定位結(jié)果，在N方向和高程方向則優(yōu)于未修正多路徑系統(tǒng)偏差的定位結(jié)果，尤其在高程方向，定位精度得到了顯著改善。修正GEO、IGSO和MEO多路徑系統(tǒng)偏差的定位結(jié)果則在N、E、U 3個(gè)方向均略優(yōu)于修正IGSO和MEO多路徑系統(tǒng)偏差的定位結(jié)果。對(duì)表1 4個(gè)測(cè)站共39組數(shù)據(jù)(GMSD存在1 d無(wú)法獲取觀測(cè)數(shù)據(jù))的定位結(jié)果進(jìn)行統(tǒng)計(jì)分析表明，大部分測(cè)站均呈現(xiàn)出上述規(guī)律。需要指出的是，存在少數(shù)幾組數(shù)據(jù)的計(jì)算結(jié)果在進(jìn)一步修正GEO衛(wèi)星的系統(tǒng)偏差后，結(jié)果甚至差于未對(duì)GEO衛(wèi)星的系統(tǒng)偏差進(jìn)行修正的結(jié)果，該原因可能是部分?jǐn)?shù)據(jù)GEO衛(wèi)星的多路徑誤差經(jīng)過弧段平均后，其仍存在較大系統(tǒng)性偏差。該系統(tǒng)性偏差將被吸收到式(1)中的常數(shù)項(xiàng)部分中，導(dǎo)致計(jì)算的修正值不準(zhǔn)確。

圖8給出了4個(gè)測(cè)站2015(213—222)所有計(jì)算結(jié)果的坐標(biāo)誤差在各個(gè)方向的坐標(biāo)誤差RMS統(tǒng)計(jì)值。由圖8可以看出，未修正多路徑誤差時(shí)，水平方向定位精度在10 cm以內(nèi)，而高程方向則達(dá)到1 m。修正IGSO和MEO衛(wèi)星觀測(cè)值的多路徑誤差后，雖然水平方向的精度有略微下降，但高程方向的精度得到了極大的改善，從1 m下降到3.5 dm左右，改善了65%。進(jìn)一步修正GEO衛(wèi)星的多路徑誤差后，相比僅修正IGSO和MEO衛(wèi)星觀測(cè)值多路徑誤差的定位結(jié)果，其在平面和高程方向都有一定程度的改善，尤其在高程方向，從3.5 dm下降到了3.0 dm左右，改善了約15%。相比不修正多路徑的情形，高程精度總體改善了75%左右。

本文中對(duì)GEO衛(wèi)星的修正采用的是后處理的方式，且需要用到雙頻的載波相位觀測(cè)值。當(dāng)采用單頻接收機(jī)進(jìn)行實(shí)時(shí)單頻PPP定位時(shí)，該方法并不適用。本文主要側(cè)重于研究GEO衛(wèi)星觀測(cè)值系統(tǒng)多路徑偏差對(duì)單頻PPP定位的影響。

圖1 北斗多路徑誤差特征Fig.1 Characteristics of BeiDou multipath

圖4 C01衛(wèi)星觀測(cè)值MP序列濾波值及修正后序列Fig.4 MP series before and after multipath correction for C01 satellite mesurements

圖5 濾波修正值對(duì)比(DOY 213 VS DOY 214)Fig.5 Comparison of filtered values for two consecutive days (DOY 213 vs DOY 214)

注：圖中百分比表示修正多路徑系統(tǒng)誤差后的MP序列相比未改正的MP序列的標(biāo)準(zhǔn)偏差改進(jìn)程度。圖6 MP序列標(biāo)準(zhǔn)偏差對(duì)比Fig.6 Comparison of Standard deviation for MP series before and after correction

圖7 定位誤差對(duì)比Fig.7 Comparison of positioning errors

圖8 定位誤差RMS值對(duì)比Fig.8 Comparison of RMS values of positioning errors

4 結(jié) 論

北斗偽距觀測(cè)值存在特有的系統(tǒng)性偏差，偏差的數(shù)量級(jí)達(dá)到幾分米到米。該系統(tǒng)偏差可分為兩類：一類是IGSO/MEO衛(wèi)星隨高度角變化的偽距系統(tǒng)性偏差；另一類是GEO衛(wèi)星(高度角僅微小變化)明顯的偽距系統(tǒng)性偏差。

本文對(duì)該系統(tǒng)偏差對(duì)單頻PPP的影響進(jìn)行了分析，表明該系統(tǒng)性偏差將導(dǎo)致單頻PPP定位結(jié)果高程方向產(chǎn)生約1m的偏差。針對(duì)GEO衛(wèi)星偽距偏差，本文提出了一種基于卡爾曼濾波的修正方法，對(duì)衛(wèi)星端和接收機(jī)端多路徑誤差進(jìn)行統(tǒng)一處理，修正后的GEO衛(wèi)星MP序列的標(biāo)準(zhǔn)偏差下降了10～16%。

結(jié)合Wanninger和Beer的高度角模型對(duì)IGSO/MEO偽距偏差進(jìn)行修正，本文分析了同時(shí)修正IGSO/MEO及GEO衛(wèi)星偽距偏差后對(duì)單頻PPP的改善情況。采用4個(gè)MGEX站10 d的數(shù)據(jù)進(jìn)行了處理，分析結(jié)果表明：僅改正和衛(wèi)星多路徑誤差，高程方向定位結(jié)果精度可改善65%左右；采用本文方法對(duì)GEO衛(wèi)星的多路徑修正后，該方向定位結(jié)果精度改善比例將進(jìn)一步提高至75%左右。

[1] 中國(guó)衛(wèi)星導(dǎo)航系統(tǒng)管理辦公室. 北斗衛(wèi)星導(dǎo)航系統(tǒng)空間信號(hào)接口控制文件[EB/OL]. 2.0版. 北京: 中國(guó)衛(wèi)星導(dǎo)航系統(tǒng)管理辦公室, [2016-10-23]. http:∥www.beidou.gov.cn/2013112127.html. CSNO. BeiDou Navigation Satellite System Signal in Space Interface Control Document, Open Service Signal[EB/OL]. (Version 2.0). Beijing: China Satellite Navigation Office[2016-10-23]. http:∥www.beidou.gov.cn/2013112127.html.

[2] 肖國(guó)銳, 隋立芬, 劉長(zhǎng)建, 等. 北斗導(dǎo)航定位系統(tǒng)單點(diǎn)定位中的一種定權(quán)方法[J]. 測(cè)繪學(xué)報(bào), 2014, 43(9): 902-907. DOI: 10.13485/j.cnki.11-2089.2014.0165. XIAO Guorui, SUI Lifen, LIU Changjian, et al. A Method of Determining the Weight Matrix for BeiDou Navigation Satellite System Single Point Positioning[J]. Acta Geodaetica et Cartographica Sinica, 2014, 43(9): 902-907. DOI: 10.13485/j.cnki.11-2089.2014.0165.

[3] 張小紅, 左翔, 李盼, 等. BDS/GPS精密單點(diǎn)定位收斂時(shí)間與定位精度的比較[J]. 測(cè)繪學(xué)報(bào), 2015, 44(3): 250-256. DOI: 10.11947/j.AGCS.2015.20130771. ZHANG Xiaohong, ZUO Xiang, LI Pan, et al. Convergence Time and Positioning Accuracy Comparison Between BDS and GPS Precise Point Positioning[J]. Acta Geodaetica et Cartographica Sinica, 2015, 44(3): 250-256. DOI: 10.11947/j.AGCS.2015.20130771.

[4] 朱永興, 馮來(lái)平, 賈小林, 等. 北斗區(qū)域?qū)Ш较到y(tǒng)的PPP精度分析[J]. 測(cè)繪學(xué)報(bào), 2015, 44(4): 377-383. DOI: 10.11947/j.AGCS.2015.20140082. ZHU Yongxing, FENG Laiping, JIA Xiaolin, et al. The PPP Precision Analysis Based on BDS Regional Navigation System[J]. Acta Geodaetica et Cartographica Sinica, 2015, 44(4): 377-383. DOI: 10.11947/j.AGCS.2015.20140082.

[5] HE Haibo,LI Jinlong, YANG Yuanxi, et al. Performance Assessment of Single-and Dual-Frequency BeiDou/GPS Single-epoch Kinematic Positioning[J]. GPS Solutions, 2013, 18(3): 393-403.

[6] TEUNISSEN P J G,ODOLINSKI R,ODIJK D. Instantaneous BeiDou+GPS RTK Positioning with High Cut-off Elevation Angles [J]. Journal of Geodesy, 2014, 88(4): 335-350.

[7] ZHANG Xiaohong, HE Xiyang. Performance Analysis of Triple-Frequency Ambiguity Resolution with BeiDou Observations[J]. GPS Solutions, 2016, 20(2): 1-13.

[8] HAUSCHILD A, MONTENBRUCK O, SLEEWAEGEN J M, et al. Characterization of Compass M-1 Signals[J]. GPS Solutions, 2012, 16(1): 117-126.

[9] MONTENBRUCK O,HAUSCHILD A, STEIGENBERGER P, et al. A COMPASS for Asia: First Experience with the BeiDou-2 Regional Navigation System[C]∥Poster at IGS Workshop. Olsztyn, Poland, 2012.

[10] MONTENBRUCK O, RIZOS C, WEBER R, et al. Getting a Grip on Multi-GNSS—the International GNSS Service MGEX Campaign[J]. GPS World, 2013, 24(7): 44-49.

[11] WANG G, DE JONG K, LI X, et al. Analysis of Characteristics of BDS Observable Combinations for Wide-lane Integer Ambiguity Resolution[C]∥SUN J, JIAO W, WU H,et al. Proceedings of China Satellite Navigation Conference (CSNC). Berlin, Heidelberg: Springer, 2014.

[12] ZHANG Xiaohong, HE Xiyang, LIU Wanke. Characteristics of Systematic Errors in the BDS Hatch-Melbourne-Wübbena Combination and Its Influence on Wide-lane Ambiguity Resolution[J]. GPS Solutions, 2016, 21(1): 265-277.

[13] WANNINGER L, BEER S. BeiDou Satellite-induced Code Pseudorange Variations: Diagnosis and Therapy[J]. GPS Solutions, 2015, 19(4): 639-648.

[14] WANG Guangxing, DE JONG K, ZHAO Qile, et al. Multipath Analysis of Code Measurements for BeiDou Geostationary Satellites[J]. GPS Solutions, 2014, 19(1): 129-139.

[15] HE Xiyang, ZHANG Xiaohong. Characteristics Analysis of BeiDou Melbourne-Wübbena Combination[C]∥SUN Jiadong, LIU Jingnan, FAN Shiwei, et al. Proceedings of the China Satellite Navigation Conference (CSNC) 2015. Berlin, Heidelberg, Germany: Springer, 2015: 587-593.

[16] LOU Yidong, GONG Xiaopeng, GU Shengfeng, et al. Assessment of Code Bias Variations of BDS Triple-frequency Signals and Their Impacts on Ambiguity Resolution for Long Baselines[J]. GPS Solutions, 2017, 21(1): 177-186.

[17] 張小紅, 李星星, 郭斐, 等. GPS單頻精密單點(diǎn)定位軟件實(shí)現(xiàn)與精度分析[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2008, 33(8): 783-787. ZHANG Xiaohong, LI Xingxing, GUO Fei, et al. Realization and Precision Analysis of Single-frequency Precise Point Positioning Software[J]. Geomatics and Information Science of Wuhan University, 2008, 33(8): 783-787.

[18] CHOY S, ZHANG K, SILCOCK D. An Evaluation of Various Ionospheric Error Mitigation Methods Used in Single Frequency PPP[J]. Journal of Global Positioning Systems, 2008, 7(1): 62-71.

[19] 涂銳, 黃觀文, 張勤, 等. 利用單基準(zhǔn)站改正信息和電離層參數(shù)估計(jì)的單頻PPP算法[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2012, 37(2): 170-173, 182. TU Rui, HUANG Guanwen, ZHANG Qin, et al. PPP Algorithm of Single Frequency Based on Corrections of Single-Base Station and Ionospheric Parameters Estimation [J]. Geomatics and Information Science of Wuhan University, 2012, 37(2): 170-173, 182.

[20] 王利, 張勤, 涂銳, 等. 基于原始觀測(cè)值的單頻精密單點(diǎn)定位算法[J]. 測(cè)繪學(xué)報(bào), 2015, 44(1): 19-25. DOI: 10.11947/j.AGCS.2015.20130315. WANG Li, ZHANG Qin, TU Rui, et al. A Kind of Single-frequency Precise Point Positioning Algorithm Based on the Raw Observations[J]. Acta Geodaetica et Cartographica Sinica, 2015, 44(1): 19-25. DOI: 10.11947/j.AGCS.2015.20130315.

[21] KEE C, PARKINSON B W. Calibration of Multipath Errors on GPS Pseudorange Measurements[C]∥Proceedings of the 7th International Technical Meeting of the Satellite Division of the Institute of Navigation. Salt Lake City, UT: [s.n.], 1994: 353-362.

[22] ESTEY L H, MEERTENS C M. TEQC: the Multi-Purpose Toolkit for GPS/GLONASS Data[J]. GPS Solutions, 3(1): 42-49.

[23] SATIRAPOD C, RIZOS C. Multipath Mitigation By Wavelet Analysis for GPS Base Station Applications[J]. Survey Review, 2005, 38(295): 2-10.

[24] SOUZA E M, MONICO J F G, POLEZEL W G C, et al. An Effective Wavelet Method to Detect and Mitigate Low-Frequency Multipath Effects[C]∥XU P, LIU J, DERMANIS A. Proceedings of International Association of Geodesy Symposia. Berlin, Heidelberg, Germany: Springer, 2008: 413-417.

[25] BORRE K, DE JONG K, PICHOT C. Subsidence Monitoring System Using Real-time GPS Sensors[C]∥Proceedings of ION-GNSS-2001, The Institute of Navigation. Salt Lake City, Utah: [s.n.], 2001: 362-367.

[26] CHOI K, BILICH A, LARSON K M, et al. Modified Sidereal Filtering: Implications for High-Rate GPS Positioning[J]. Geophysical Research Letters, 2004, 31(22): L22608.

[27] MA X, SHEN Y. Multipath Analysis of COMPASS Triple Frequency Observations[C]∥The International Symposium on GPS/GNSS. Xi’an: [s.n.], 2012.

[28] SAASTAMOINEN J. Atmospheric Correction for the Troposphere and Stratosphere in Radio Ranging Satellites[M]∥HENRIKSEN S W, MANCINI A, CHOVITZ B H. The Use of Artificial Satellites for Geodesy. Virginia: AGU, 1972: 274-251.

[29] NIELL A E. Global Mapping Functions for the Atmosphere Delay at Radio Wavelengths[J]. Journal of Geophysical Research: Solid Earth, 1996, 101(B2): 3227-3246.

(責(zé)任編輯：陳品馨)

Correction Model of BeiDou Code Systematic Multipath Errors and Its Impacts on Single-frequency PPP

WANG Jie1,HE Xiyang2

1． Naval University of Engineering,Wuhan 430033，China； 2． Qian Xun Wei Zhi Network Co. Ltd.，Shanghai 200438, China

There are systematic multipath errors on BeiDou code measurements, which are range from several decimeters to larger than 1 meter. They can be divided into two categories, which are systematic variances in IGSO/MEO code measurement and in GEO code measurement. In this contribution, a methodology of correcting BeiDou GEO code multipath is proposed base on Kalman filter algorithm. The standard deviation of GEO MP Series decreases about 10%～16% after correction. The weight of code in single-frequency PPP is great, therefore, code systematic multipath errors have impact on single-frequency PPP. Our analysis indicate that about 1 m bias will be caused by these systematic errors. Then, we evaluated the improvement of single-frequency PPP accuracy after code multipath correction. The systematic errors of GEO code measurements are corrected by applying our proposed Kalman filter method. The systematic errors of IGSO and MEO code measurements are corrected by applying elevation-dependent model proposed by Wanninger and Beer. Ten days observations of four MGEX (Multi-GNSS Experiment) stations are processed. The results indicate that the single-frequency PPP accuracy can be improved remarkably by applying code multipath correction. The accuracy in up direction can be improved by 65% after IGSO and MEO code multipath correction. By applying GEO code multipath correction, the accuracy in up direction can be further improved by 15%.

BeiDou; single-frequency precise point positioning; Kalman filter; multipath; pseudo-range; accuracy analysis

The National Natural Science Foudation of China(Nos.41631072；41504029)

WANG Jie(1973—),male,postgraduate,associate professor,majors in radio and satellite navigation.

汪捷，何錫揚(yáng).北斗衛(wèi)星偽距多路徑系統(tǒng)偏差的修正方法及其對(duì)單頻PPP的影響分析[J].測(cè)繪學(xué)報(bào)，2017,46(7)：841-847.

10.11947/j.AGCS.2017.20170020. WANG Jie,HE Xiyang.Correction Model of BeiDou Code Systematic Multipath Errors and Its Impacts on Single-frequency PPP[J]. Acta Geodaetica et Cartographica Sinica,2017,46(7):841-847. DOI:10.11947/j.AGCS.2017.20170020.

P228

A

1001-1595(2017)07-0841-07

國(guó)家自然科學(xué)基金(41631072；41504029)

2017-04-24

汪 捷(1973—)，男，碩士，副教授，研究方向?yàn)樾l(wèi)星及無(wú)線電導(dǎo)航。

E-mail： wangjie7312@163.com

修回日期： 2017-06-05