賀 劉，姚 錚，崔曉偉，陸明泉，郭 婧

（1.清華大學(xué) 電子工程系，北京100084；2.中國民航科學(xué)技術(shù)研究院，北京100084）

隨著GPS與GLONASS的現(xiàn)代化，Galileo與我國北斗系統(tǒng)的建設(shè)，全球GNSS系統(tǒng)進(jìn)入了蓬勃發(fā)展時期，多模、多頻聯(lián)合導(dǎo)航已成為導(dǎo)航系統(tǒng)發(fā)展的主要趨勢［1］.由于測距源數(shù)量的倍增，客觀上增大了測距出現(xiàn)故障偏差的概率，熱噪聲、多徑、干擾等復(fù)雜環(huán)境將削弱測距的可靠性，觀測測距的完好性日趨重要.接收機完好性監(jiān)測（receiver autonomous integrity monitoring，RAIM）技術(shù)作為用戶端唯一的完好性監(jiān)測手段，通過故障檢測與排除（fault detect execution，F(xiàn)DE）算法，利用冗余觀測信息進(jìn)行一致性統(tǒng)計校驗，剔除含粗差的故障測距，以提升測距置信度水平，進(jìn)而保證定位估計的準(zhǔn)確性.

由于受測距噪聲與接收設(shè)備熱噪聲的影響，測距觀測精度不完全相同，方差隨著俯仰角度及C／N0等多種因素相應(yīng)變化.在傳統(tǒng)加權(quán)RAIM 算法與模型 中， NIO-RAIM （novel integrity-optimized RAIM）算法［2］針對衛(wèi)星特征斜率提出最優(yōu)加權(quán)，可以有效地提升系統(tǒng)可用性，但會犧牲一定的定位精度，且復(fù)雜度較高，不易在側(cè)重定位域的實時接收機中使 用；廣 域 增 強 系 統(tǒng)（wide area augmentation system，WAAS）加權(quán)模型［3］的適用范圍較窄，僅限于GPS，不適用于多模多頻場景；雙頻加權(quán)模型［4］基于服務(wù)體積模型（service volume model，SVM），涵蓋了殘余電離層延遲、對流層延遲及多徑等主要誤差源，適用于GPS／Galileo雙模雙頻場景，但頻點間加權(quán)因子等權(quán)，帶來了一定的熱噪聲模型偏差.

傳統(tǒng)加權(quán)模型均未考慮信號體制對接收機熱噪聲的影響.在接收機偽碼跟蹤環(huán)路中，復(fù)合二進(jìn)制偏移載波（composite binary offset carrier，CBOC）、交替二進(jìn)制偏移載波（alternate binary offset carrier，AltBOC）等新型調(diào)制方式與二進(jìn)制相移鍵控（binary phase shift keying，BPSK）類傳統(tǒng)調(diào)制方式的碼跟蹤精度不同.尤其在雙頻組合場景中，采用雙頻定位解算方法雖然能夠有效地抑制電離層延遲，但對各信號熱噪聲及多徑延遲有一定的放大作用，這將加大信號體制間觀測量的精度差距，造成測距觀測精度不等權(quán)現(xiàn)象.

傳統(tǒng)加權(quán)模型忽略了信號體制原因帶來的測距精度不等權(quán)問題，該模型假設(shè)與實際不符的情況將對加權(quán)RAIM 技術(shù)性能的評估帶來偏差.本文在噪聲干擾層面，通過仿真分析了幾種典型信號在偽碼跟蹤延遲鎖相環(huán)路（delay-locked loop，DLL）跟蹤精度的差異性，根據(jù)信號調(diào)制方式、接收機設(shè)計參數(shù)和C／N0，建立新型接收機偽距誤差熱噪聲模型，推導(dǎo)了各頻點測距加權(quán)因子的產(chǎn)生方式，進(jìn)而細(xì)化完善了新型雙頻加權(quán)模型.仿真測試表明，該模型能夠更加合理地反映雙頻測距的實際噪聲，優(yōu)化了雙頻加權(quán)RAIM 算法的FDE 性能，提升了測距置信度水平并有效保證了定位精度，適用于不同信號體制共用的雙頻定位域場景.

1 信號體制對碼偽距測量精度的影響

1.1 多信號體制共用的現(xiàn)實趨勢

隨著信號結(jié)構(gòu)的設(shè)計與發(fā)展，出現(xiàn)了新一代信號體制，如在二進(jìn)制偏移載波（binary offset carrier，BOC）信號基礎(chǔ)上發(fā)展出的CBOC、AltBOC 等新型導(dǎo)航信號調(diào)制方式.新信號易于實現(xiàn)頻段資源復(fù)用、利于信號捕獲，具有降低偽碼跟蹤誤差、抑制信號多徑等顯著優(yōu)勢，已廣泛應(yīng)用于GNSS 各頻點中.傳統(tǒng)BPSK-R 信號將在導(dǎo)航系統(tǒng)中繼續(xù)使用，多種導(dǎo)航信號體制的共用與互操作已成為GNSS系統(tǒng)的現(xiàn)實應(yīng)用趨勢.表1給出目前GNSS系統(tǒng)采用的主要信號調(diào)制方式.表中，BPSK-R（n）指具有（n ×1.023） MHz碼片速率的BPSK-R 信 號；BOC（a，b）指一個（a ×1.023） MHz的方波頻率與一個（b ×1.023） MHz碼片速率產(chǎn)生的BOC 調(diào)制信號；CBOC（6，1，1／11）信號是BOC（6，1）與BOC（1，1）通過時分復(fù)用方式，在導(dǎo)頻通道中按照1／10功率比組合而成；類似地，BOC（6，1）與BOC（1，1）在時間復(fù)用二進(jìn)制偏置載波（time-multiplexed BOC，TMBOC）（6，1，4／33）信號導(dǎo)頻通道中的功率分配比是4／33.

表1 GNSS廣泛采用的信號體制Tab.1 Widely adopted signal structures of GNSS

1.2 信號體制間碼跟蹤精度的差異性

為了分析不同信號體制在接收機碼環(huán)中的跟蹤精度，選取廣泛使用的BOC（14，2）、CBOC（6，1，1／11）、AltBOC（15，10）與BPSK-（1）、BPSK-（10）進(jìn)行分析.分析的跟蹤方法使用非相干超前和滯后（noncoherent early-late processing，NELP）算法.

根據(jù)碼跟蹤誤差理論估計理論［5-9］可知，信號通過DLL環(huán)路鑒別器后，輸出熱噪聲鑒別誤差方差為

式中：c為光速，Bfe為接收機前端帶寬，Ss（f）為信號調(diào)制方式的歸一化功率譜密度，T 為預(yù)檢測積分時間，dt為相關(guān)器超前、滯后間隔.信號通過環(huán)路濾波器后，DLL環(huán)路碼跟蹤誤差σDLL為

式中：BDLL為濾波器帶寬，σDLL為信號在單位測距標(biāo)準(zhǔn)差σ下的偽距測量精度.

σDLL與接收機設(shè)計參數(shù)相關(guān)，假設(shè)濾波器帶寬BDLL為5Hz，預(yù)檢測積分時間T 為20ms，AltBOC（15，10）前端帶寬為51 MHz，其余信號前端帶寬為18MHz.通過仿真得到各信號的碼環(huán)跟蹤誤差理論下限值，如圖1所示.可以看出，不同信號調(diào)制方式下的σDLL差別明顯，尤其在C／N0較低的弱信號條件下，40dB-Hz時BPSK-（1）的碼跟蹤誤差與其他調(diào)制信號相比已超出0.5m.

圖1 各信號的碼跟蹤誤差Fig.1 Code tracking error of signal structures

2 新型雙頻加權(quán)模型

通過處理捕獲信號的直接擴頻序列，實現(xiàn)接收機碼片同步，進(jìn)而可以獲得偽距觀測值.碼偽距作為基本觀測量，置信度是影響PVT 定位精度的主要因素.理想測距為衛(wèi)星與用戶間的真實幾何距離，而實際測距含有噪聲誤差.一般地，測距噪聲服從期望為零、方差為σ2的高斯分布.噪聲方差越小，最終的定位結(jié)果越接近用戶的真實位置.測距噪聲方差是多種誤差方差的期望之和，可以表示為

式中：σURA為星歷誤差，σiono為電離層延遲，σtropo為對流層 延 遲，σmpth為 多 徑 誤 差，σrevr為 接 收 機 熱 噪聲誤差.

2.1 傳統(tǒng)等權(quán)值雙頻模型

為了獲得更精確的定位解，可以通過2種方法抑制或部分消除噪聲的影響：一是通過算法或模型消除部分噪聲誤差；二是通過加權(quán)算法，以提高或降低測距權(quán)值的方式歸一化噪聲方差.在雙頻組合中，依靠頻點測距作差能夠消除大部分的電離層延遲，但熱噪聲屬于系統(tǒng)誤差，很難消除.在GPS／Galileo雙模雙頻場景中，傳統(tǒng)加權(quán)模型和各類誤差方差模型［10-12］可以表示為

式中：測 距 間 用 戶 距 離 精 度（user range error，URA）可取1、0.5、0.2m；σ2L1L2為0.04，表示雙頻誤差參數(shù).對流層延遲模型σi，tropo為

表示各頻點無電離層延遲模型，其中各頻點測距加權(quán)因子相等：

熱噪聲模型RMSpr-air與多徑延遲模型與仰角Ei相關(guān)：

2.2 新型接收機偽距誤差熱噪聲模型

在雙頻組合的實際應(yīng)用中，一般可以假設(shè)URA誤差相等，由于頻點間信號經(jīng)過近似相同的對流層路徑，σtropo相等，σrevr及σmpth相互獨立.

根據(jù)1節(jié)的分析可知，不同信號體制下碼跟蹤誤差存在明顯差異，各頻點測距熱噪聲不等.由式（4.d）、（4.e）、（4.f）可知，傳統(tǒng)模型中熱噪聲方差RMSpr＿air與仰角相關(guān)，且假設(shè)各頻點測距σi，L1與σi，L2相等，這與接收機DLL環(huán)路熱噪聲的實際情況存在一定的偏差.傳統(tǒng)模型的頻率方差σ2L1L2為固定值，不能靈活地兼容各系統(tǒng)各頻點中心頻率不同的情況.

由于傳統(tǒng)模型不能兼顧不同調(diào)制方式組合使用的情況，尤其在如今新型信號和傳統(tǒng)信號共用與互操作廣泛應(yīng)用時期.為了改進(jìn)、修正傳統(tǒng)模型中接收機熱噪聲模型的偏差問題，根據(jù)信號體制的差異性，由式（1）、（2）可得基于信號體制的新型熱噪聲方差模型：

新型熱噪聲模型與信號體制功率譜密度、接收機碼環(huán)設(shè)計參數(shù)及C／N0相關(guān)，將該模型替換傳統(tǒng)模型中的RMS2pr＿air與σ2L1L2項，各頻點測距的加權(quán)因子σ′i為

2.3 新型雙頻加權(quán)模型描述

雙頻組合能夠消除大部分電離層誤差，無電離層偽距ρi 的數(shù)學(xué)模型［13］為

式中：f1、f2為各頻點的中心頻率，ρL1、ρL2為各頻點的偽距.對式（8）求期望后，可得雙頻加權(quán)方差為

由式（9）可見，雙頻組合將增大接收機熱噪聲和多徑誤差.當(dāng)σ′i，L1與σ′i，L2一定時，2個頻率的間隔越大，對誤差的放大作用越?。划?dāng)2個頻點的偽距跟蹤精度不同時，高頻頻點測距誤差所占的比重更大.

綜上所述，基于信號體制優(yōu)化的雙頻加權(quán)模型為

σURA使用傳統(tǒng)估計值.

3 新模型對RAIM 技術(shù)的影響

新型雙頻模型作為一項基礎(chǔ)模型，可以廣泛應(yīng)用于加權(quán)定位估計與加權(quán)RAIM 技術(shù).

根據(jù)加權(quán)最小二乘定位解算模型可知，加權(quán)最小二乘估計解［11］為

式中：R 為歸一化噪聲方差加權(quán)矩陣，

加權(quán)矩陣實現(xiàn)了測距誤差的歸一化處理.由式（11）可得歸一化的偽距殘差矢量w 與映射矩陣U：

w 各分量相互獨立，可得w＝URy＝URn，其中n 為噪聲矢量，服從期望為0、方差為1的高斯分布.由于w 的幅值體現(xiàn)了觀測值間的一致性，可將歸一化的偽距殘差平方和作為檢驗統(tǒng)計量［14］，雙模條件下服從自由度為（n-5）的χ2分布.根據(jù)用戶給定的虛警率Pfa能夠確定歸一化的檢驗門限值T，進(jìn)而判決測距是否存在粗差.

4 仿真驗證與分析

通過仿真得到GPS24 星星座及Galileo27／3／1Walker星座，仰角遮蔽角為5°，時長1d，采樣步長為1min，獲得1 440個采樣點.以5°×5°為采樣區(qū)間，遍歷得到全球范圍可見衛(wèi)星數(shù)的分布情況，如圖2所示.圖中，SatN為可見衛(wèi)星數(shù)，Pr為分布概率.

按照航空無線電委員會（Radio Technical Commission for Aeronautics，RTCA）規(guī)定的“最小操作性能標(biāo)準(zhǔn)”［15］要求，遍歷全球24個地點.測距標(biāo)準(zhǔn)差σ統(tǒng)一取2.5m，URA 取1m，Pfa取0.01.假設(shè)各系統(tǒng)雙頻組合為：GPS L1 C／A 碼BPSK-R（1），中心頻點頻率為1 575.42 MHz，L5BPSK-R（10），中心頻點頻率為1 176.45 MHz；Galileo E1 CBOC（6，1，1／11），中心頻點為1 575.42 MHz，E5a AltBOC（15，10），中心頻點為1 191.795 MHz.

圖2 可見衛(wèi)星數(shù)分布Fig.2 Distribution of number of visible satellites

為了仿真對比各模型對全球RAIM 技術(shù)FDE平均性能的影響.設(shè)置3種處理方法，方法a為單頻處理方法，方法b為雙頻等權(quán)處理方法，方法c為本文提出方法.由于雙模或多模場景的漏檢率極低［16］，考慮故障檢測率和虛警率指標(biāo).為了便于直觀顯示，在虛警量級接近的條件下，對新模型與傳統(tǒng)模型進(jìn)行比較.表2給出各場景3～10σ的全球24地點的平均FDE 性能.單頻方法的故障檢測率較低，虛警率極低.雙頻方法的故障檢測率大幅高于單頻方法，本文方法優(yōu)化了加權(quán)RAIM 技術(shù)的FDE性能，c方法的故障檢測率優(yōu)于b方法，尤其在故障較小時更加顯著.

表2 平均FDE性能Tab.2 Average FDE performance

在上述仿真環(huán)境下，經(jīng)過FDE 步驟，分析各方法的定位精度.選擇北京觀測站，隨機選取一個健康衛(wèi)星，在觀測測距的各個采樣歷元中均加入±5 m至±25m 的隨機突變偏差，得到圖3所示的各場景定位 精 度，三 維 均 方 根 誤 差（root mean square，RMS）列于圖3的右下角.圖中，Lat為緯度，Lng為經(jīng)度.由于故障檢測率較低，含故障偏差測距影響了觀測測距置信度，方法1的RMS較高，為4.41m；方法2的RMS為3.37m；經(jīng)過本文提出模型的進(jìn)一步優(yōu)化，方法3的RMS降至2.96m，與雙頻等權(quán)方法相比，定位精度提升約0.4 m.需要指出，方法2、3是經(jīng)過RAIM 故障檢測與排除，獲得較高測距置信度，并對測距噪聲歸一化加權(quán)定位解算共同優(yōu)化的結(jié)果.

為了不失一般性，表3給出各方法在全球24地點的平均RMS結(jié)果.各方法的定位精度與北京觀測站定位精度相近，可以看出，采用本文模型可以有效地保證定位精度.

圖3 方法1～3的定位精度Fig.3 Positioning accuracy of method 1～3

表3 全球24地點平均RMSTab.3 Average RMS of 24worldwide locations m

5 結(jié) 語

本文在分析幾種典型信號調(diào)制方式在接收機DLL環(huán)路碼跟蹤精度的基礎(chǔ)上，提出基于信號體制的接收機熱噪聲加權(quán)模型.結(jié)合傳統(tǒng)雙頻組合各噪聲方差估計，分析新模型對雙頻加權(quán)RAIM 技術(shù)性能及定位估計精度的影響.仿真測試表明，采用本文模型能夠優(yōu)化提高RAIM 技術(shù)FDE 性能，提升各頻點測距置信度水平，并有效地保證定位精度，適合于信號體制共用條件下的定位域應(yīng)用場景.作為基礎(chǔ)加權(quán)模型，該模型能夠適應(yīng)信號體制多樣化的展趨勢，具有廣泛適用性的特點.

（

）：

［1］郭婧.多模多頻衛(wèi)星導(dǎo)航系統(tǒng)RAIM 技術(shù)研究［D］.北京：清華大學(xué)，2011：1-12.GUO Jing.Research on RAIM for multi-constellation and multi-frequency satellite navigation systems ［D］.Beijing：Tsinghua University，2011：1-12.

［2］MADONNA P，VIOLA S，SFARZO L.NIORAIM algorithm applied to a multiconstellation GNSS：analysis of integrity monitoring performances in various phases of flight［C］∥Position Location and Navigation Symposium（PLANS）.Palm Springs：IEEE，2010：1258-1263.

［3］WALTER T，ENGE P.Weighted RAIM for precision approach［C］∥ION GPS 1995.Palm Springs：［s.n.］，1995：1995-2004.

［4］HARTINGER H，BRUNNER F K.Variances of GPS phase observations：the SIGMA-εmodel［J］.GPS Solutions，1999，2（4）：35-43.

［5］BETZ J W，KOLODZIEJSKI K R.Generalized theory of code tracking with an early-late discriminator Part I：lower bound and coherent processing［J］.IEEE Transactions on Aerospace and Electronic Systems，2009，45（4）：1538-1550.

［6］BETZ J W，KOLODZIEJSKI K R.Generalized theory of code tracking with an early-late discriminator part II：noncoherent processing and numerical results［J］.IEEE Transactions on Aerospace and Electronic Systems，2009，45（4）：1557-1564.

［7］何在民.衛(wèi)星導(dǎo)航信號碼跟蹤精度研究［D］.北京：中國科學(xué)院研究生院，2012：35-43.HE Zai-min.Research on code tracking accuracy for satellite navigation signals［D］.Beijing：Chinese Academy of Sciences，2012：35-43.

［8］LIU A，AN J，WANG A.Effect of transform domain interference suppression on PN code tracking loops［C］∥2010 12th IEEE International Conference on Communication Technology.Nanjing：IEEE，2010：1064-1067.

［9］LIU Y，RAN Y，KE T，et al.Code tracking performance analysis of GNSS signal in the presence of CW interference［J］.Signal Processing，2011，91（4）：970-987.

［10］RTCA DO-229D，Minimum operational performance standards for global positioning system／wide area augmentation system airborne equipment［S］.Washington，DC：RTCA，2006.

［11］KAPLAN E D，HEGARTY C J.GPS 原 理 與 應(yīng) 用［M］.寇艷紅，譯.北京：電子工業(yè)出版社，2012：258-265.

［12］LEE Y，BRAFF R，F(xiàn)ERNOW J，et al.GPS and Galileo with RAIM or WAAS for vertically guided approaches［C］∥Proceedings of the Institute of Navigation GNSS.Long Beach：Institute of Navigation，2005：1801-1825.

［13］劉家興.衛(wèi)星定位精度和自主完好性的改進(jìn)［D］.北京：清華大學(xué)，2008：25-27.LIU Jia-xing.Improvement of satellite positioning sccuracy and sutonomous integrity［D］.Beijing：Tsinghua University，2008：25-27.

［14］代長勇，常青，陳媛.多系統(tǒng)衛(wèi)星導(dǎo)航接收機自主完好性監(jiān)測［J］.遙測遙控，2010，31（6）：11-15.DAI Chang-yong， CHANG Qing， CHEN Yuan.Multi-system satellite navigation receiver autonomous integrity monitoring［J］.Journal of Telemetry，Tracking and Command，2010，31（6）：11-15.

［15］RTCA DO-229D，Minimum operational performance standards for airborne supplemental navigation equipment using global positioning system （GPS）.Washington，DC：RTCA，1991.

［16］許龍霞，張慧君，李孝輝.多模衛(wèi)星導(dǎo)航系統(tǒng)的RAIM算法研究［J］.時間頻率學(xué)報，2011，34（2）：131-138.XU Long-xia，ZHANG Hui-jun，LI Xiao-hui.Study on RAIM algorithm of multi-mode navigation system ［J］.Journal of Time and Frequency，2011，34（2）：131-138.