李 強，朱衍波，方 堃（民航數(shù)據(jù)通信有限責任公司，北京 0009；北京航空航天大學，北京 009）

基于碼減載波和B值的GBAS精度性能分析

李 強1，朱衍波1，方 堃2
（1民航數(shù)據(jù)通信有限責任公司，北京 1000191；2北京航空航天大學，北京 100191）

對GBAS精度性能進行了分析。分別按照基于碼減載波 （code minus carrier，CMC）和B值兩種方法對GBAS基準站實際數(shù)據(jù)進行了分析。CMC分析結果和B值的分析結果一致表明，GBAS基準站的精度性能等級約為GAD-A2。CMC分析法與B值分析方法具有一些不同的特點。CMC需要使用雙頻觀測量；B值只需要單頻觀測量。CMC消除了短周期的多路徑效應，主要反映的是長周期的多路徑效應；B值經(jīng)過GBAS基準站接收機平均，反映了不同接收機組合的多路徑綜合效應，更能反映GBAS差分校正誤差的實際影響。除此之外，CMC分析法與B值分析法的天空圖可以直觀地評估GBAS站址的多路徑環(huán)境。

地基增強系統(tǒng)；地基增強系統(tǒng)精度性能評估；碼減載波；B值；多路徑

衛(wèi)星導航地基增強系統(tǒng)（ground based augmenta tion system，GBAS）是滿足民航業(yè)務需求的重要系統(tǒng)之一。GBAS系統(tǒng)向用戶發(fā)播差分校正參數(shù)、校正標準差參數(shù)等。這些參數(shù)的準確度和可信度是用戶保持定位精度和完好性的重要保證。為了保證GBAS正常工作，并提供高精度的差分數(shù)據(jù)和完好性參數(shù)，國際標準組織規(guī)定了GBAS基準站的運行性能標準。GBAS精度性能是其中的重要方面之一[1-3]。

Collins的Mc Graw博士詳細討論了LAAS精度的理論模型和相關假設[4]，并提出了LAAS地面系統(tǒng)精度和機載性能精度等級的參數(shù)化模型。除了理論分析外，GBAS精度性能也可通過長期監(jiān)測其性能獲得其統(tǒng)計數(shù)據(jù)進行評估。Braasch提出分析碼多路徑的方法分析了單接收機的多路徑環(huán)境對接收機碼精度的影響[5]。美國聯(lián)邦航空局FAA的LAAS項目組建立了監(jiān)測接收機長期分析LAAS基準接收機的運行性能[6]。

北京航空航天大學建立了試驗型的GBAS基準站，自主研發(fā)了GBAS基準站軟件，可實時監(jiān)測GPS衛(wèi)星狀態(tài)和GPS差分校正值等，積累了實際數(shù)據(jù)。本文基于GBAS數(shù)據(jù)評估GBAS精度性能，采用了碼減載波（CMC）和B值兩種方法，分析了GBAS精度性能。并根據(jù)CMC和B值分析了GBAS基準站的多路徑環(huán)境影響。

GBAS基準站多路徑環(huán)境影響提供了GBAS基準站選址評估依據(jù)。GBAS精度性能分析可用于評估GBAS接收機的運行等級，提供確定GBAS差分值標準差的依據(jù)之一。

1 北航試驗型GBAS站精度分析方法

1.1 北航試驗型GBAS站

北航建立了試驗型的GBAS基準站，如圖1所示。如心樓GBAS處理器計算本地差分和完好性信息并實時發(fā)送給新主樓的GBAS性能監(jiān)測處理器。

圖1 北航試驗型GBAS基準站系統(tǒng)框圖（左）和天線（右）Fig.1 GBAS prototype at BUAA（left）and GPS receiver（right）

GBAS基準站包括3臺基準型GPS接收機，連接至3副間距約為10 m的天線，分別為接收機A、B、C。新主樓的性能監(jiān)測接收機距基準站距離約為200 m。

通過一段時間運行數(shù)據(jù)的積累，可通過原始數(shù)據(jù)和B值分析GBAS精度的性能。

國際標準組織RTCA根據(jù)GNSS接收機的技術現(xiàn)狀和未來趨勢，結合GBAS性能需求，定義了GBAS地面站系統(tǒng)的GPS性能標準的形式為[5]

其中：M為基準接收機數(shù)目；θi指第i顆測距源的仰角；a0、a1、a2、θ0是精度等級的參數(shù)，精度等級包括A、B、C，定義如表1所示。

表1 GBAS-GPS精度需求參數(shù)Tab.1 GBAS-GPS accuracy requirement parameters

1.2 GBAS精度分析方法

GBAS精度分析采用兩種方法。第1種方法是基于碼減載波進行計算。

GPS偽距和相位的基本公式如下，

偽距

相位

其中：D為衛(wèi)星到接收機幾何距離；I1為載頻L1上的電離層延遲；URE為廣播星歷誤差；MP1為偽距P1的多路徑誤差；N1為相位L1的載波相位模糊度；υ1為相位L1的多路徑誤差，可忽略。

碼減載波（CMC）為CMC1包括電離層和L1相位模糊度的影響，可結合L2相位消除電離層影響

結合CMC1，可獲得消除了相位模糊度參數(shù)和電離層延遲影響的P1上的多路徑MP1估值。

第2種方法為B值分析方法。

GBAS的B值計算公式為

B值統(tǒng)計特征的理論值為

上述公式中的偽距為采用100 s長度的載波平滑偽距，短時間的B值具有相關性，因此在分析中按照200 s間隔采樣，可基本保證B值樣本是獨立的。

2 GBAS數(shù)據(jù)結果與精度分析

采用了2012年6月1日—6月4日之間的數(shù)據(jù)。GBAS基準站在該時間段，接收機C由于故障沒有接收數(shù)據(jù)，GBAS基準站接收和記錄接收機A與接收機B的數(shù)據(jù)。這4天的結果類似，下文主要以6月1日的結果為例說明分析過程和結果。

接收機A的CMC包括了2012年6月1日1天所有可見衛(wèi)星，隨時間和仰角變化如圖2所示，上圖為CMC隨時間變化，下圖為CMC隨仰角變化。可以看出，由于低仰角衛(wèi)星受到多路徑影響大，CMC在低仰角變化幅度較大，在高仰角變化幅度較小。

圖2 接收機A的CMC隨時間變化（上）和隨仰角變化（下）Fig.2 CMC time series（top）and CMC variation with elevation angle（bottom）

按照仰角5°間隔對接收機A的CMC數(shù)據(jù)進行分區(qū)統(tǒng)計，統(tǒng)計結果如圖3所示。上圖中疊加的深實線為標準規(guī)定的GAD-A2曲線，淺實線為GAD-B2曲線。下圖為每個仰角區(qū)間的樣本數(shù)。CMC均值在各仰角區(qū)間約為0。CMC標準差略小于GAD-A2標準，略大于GAD-B2標準。低仰角樣本數(shù)多于高仰角樣本數(shù)。

圖3 接收機A的CMC標準差隨仰角變化（上）和各仰角區(qū)間的樣本數(shù)（下）Fig.3 CMC standard variation with respect to elevation angle（top）and number of independent sample with respect to elevation angle（bottom）

統(tǒng)計2012年6月1日—6月4日接收機A的CMC標準差隨仰角變化，得到表2。從表2中可以看出，CMC標準差在數(shù)日內的變化比較穩(wěn)定。

接收機A的CMC標準差隨仰角變化的1～4日均值與GAD-A2相應值的對比曲線如圖4所示。從圖4中可以看出，CMC標準差的均值均在GAD-A的相應值以下，并留有比較大的冗余。

將接收機A的CMC隨仰角、方位角進行統(tǒng)計，可得CMC隨仰角、方位角的天空圖，如圖5所示。CMC天空圖直觀地顯示了多路徑在不同視角方向的影響，低仰角衛(wèi)星受到的多路徑影響比高仰角的嚴重。

表2 接收機A的CMC標準差（6月1—4日）隨仰角變化Tab.2 CMC standard variation with respect to elevation angle of receiver A in June 1-4 day，2012

圖4 接收機A的CMC標準差均值（6月1—4日）隨仰角變化Fig.4 CMC standard variation with respect to elevation angle of receiver A from in 1-4 day，2012

接收機B的CMC數(shù)據(jù)結果類似于接收機A，CMC標準差略低于GAD-A2標準，高于GAD-B2標準，不再贅述。

GBAS基準站計算的2012年6月1日的接收機A的所有衛(wèi)星B值如圖6所示。B值的標準差隨仰角增大而減小。

按照5°仰角間隔，對B值統(tǒng)計，并根據(jù)B值理論標準差與GBAS精度的關系換算后的B值標準差如圖7所示。上圖顯示，B值標準差同樣位于GAD-B2 和GAD-A2之間；下圖為B值樣本在各仰角區(qū)間的樣本數(shù)量。

圖5 接收機A的CMC天空圖Fig.5 Skymap of CMC statistics at receirer A

圖6 接收機A的B值隨時間變化（上）和隨仰角變化（下）Fig.6 B values time series（top）and B values variation with elevation angle（bottom）

圖7 接收機A的B值標準差隨仰角變化（上）和各仰角區(qū)間的樣本數(shù)（下）Fig.7 B values standard variation with respect to elevation angle（top）and number of independent sample with respect to elevation angle（bottom）

根據(jù)2012年6月1日—6月4日的B值數(shù)據(jù)統(tǒng)計得到標準差隨仰角的變化如表3所示。

表3 B值統(tǒng)計的精度標準差表Tab.3 Accuracy from B values vs GAD-A2

3 結語

GBAS精度性能監(jiān)測對評估GBAS運行性能具有重要意義，本文對此利用實際數(shù)據(jù)進行了分析。數(shù)據(jù)來源于北航試驗型的GBAS基準站，該站對接收的GPS信號進行GBAS處理，記錄原始數(shù)據(jù)和處理數(shù)據(jù)。本文利用2012年6月1日—6月4日的GBAS原始數(shù)據(jù)和處理數(shù)據(jù)分析了GBAS精度性能。

采取基于CMC和B值分析兩種分析方法。這兩種方法獲得GBAS精度性能結果一致，并均能反映GBAS站址的多路徑環(huán)境影響。兩種方法的不同點：

1）B值只需要單頻觀測量，CMC需要使用雙頻觀測量；

2）CMC消除模糊度后主要反映長周期的多路徑效應；

3）B值經(jīng)過更過的GBAS基準站的接收機平均后，綜合了接收機組合后的多路徑效應，更能反映GBAS差分校正誤差的實際大小。

分析結果表明，北航的GBAS基準站接收機的精度性能介于GAD-B2和GAD-A2之間，按照保守取法，可認為精度性能為GAD-A2。

[1]Minimum Operational Performance Specification for Global Navigation Satellite Ground Based Augmentation System Ground Equipment toSupport CategoryⅠOperations，ED-114[S].2003.

[2] NERI P，MACABIAU C，AZOULAI L，et al.GBAS NSE Model for CATⅡ/ⅢAutoland Simulations[C]//Proceedings of IEEE/ION Plans. Indian Wells，CA，2010：694-707.

[3]SHIVELY C A，HSIAOT T.Availability of GAST D GBAS Considering Continuity of Airborne Monitors[C]//Proceedings of the 2010 international technical meeting of the institute of navigation.San Diego，CA，2010：365-375.

[4]GARY A MCGRAW.Development of the LAAS Accuracy Models[C]// ION GPS，2000.

[5]BRAASCH M.Global Positioning System：TheoryandApplications[M]. Washington D C：American Institute of Aeronautics and Astronautics，1996.

[6] FAA.Local Area Augmentation System Performance Analysis and Activities Report[R].

[7]Minimum Aviation System Performance Standards（MASPS）for the Local Area Augmentation System（LAAS），RTCA/DO-245[S].

（責任編輯：黨亞茹）

Performance analysis of GBAS based on code minus carrier and B values

LI Qiang1，ZHU Yan-bo1，F(xiàn)ANG Kun2
（1.Aviation Data Communication Corporation，Beijing 100191，China；2.Beihang University，Beijing 100191，China）

The paper assessed GBAS accuracy performance based on two methods.A GBAS prototype was developed and used to collect and process GPS data.Based on the raw data，two methods called CMC and B values were used to analysis to GBAS accuracy performance.The result of both methods showed consistently that the GBAS prototype accuracy was on the level of GAD-A2 class.Compared with each other，CMC needs to use observations of dual frequencies，however，B values rely on single frequency.Secondly，CMC suppressed short-term multipath errors，and long-term multipath effects remained.B values，however，mainly consist of composite effects of multiple receivers′errors after averaging across reference receivers，which reflect realistic errors of GBAS corrections as both occupied similar processes.Lastly，sky-maps of root mean squares of CMC and B values reflected surrounding multipath effects of GBAS sites.

GBAS；GBAS accuracy performance assessment；code minus carrier；B values；multipath

TN967.1

A

1674－5590（2013）01－0009－04

2012-05-04；

2012-08-10

國家重點基礎研究發(fā)展計劃（973計劃）（2011CB707004）；國家科技支撐計劃（2011BAH24B02）；國家自然科學基金項目（61179054）

李 強（1981—），男，湖北荊門人，工程師，碩士，研究方向為衛(wèi)星導航及其應用.