呂志成，劉增軍，王飛雪

（國防科學技術(shù)大學 電子科學與工程學院衛(wèi)星導航研發(fā)中心，湖南 長沙410073）

0 引 言

GPS由美國國防部于1973年開始研發(fā)，經(jīng)過20多年的研制和試驗，于1995年初具備完全的工作能力［1]。經(jīng)過幾十年的不斷發(fā)展，GPS作為新一代的衛(wèi)星導航定位系統(tǒng)已在軍事、交通運輸、測量測繪、高精度時間比對以及資源勘探等諸多領(lǐng)域得到了廣泛的應用。為了進一步鞏固GPS在衛(wèi)星導航領(lǐng)域的領(lǐng)先優(yōu)勢和提高系統(tǒng)服務質(zhì)量，美國政府在1999年提出了GPS現(xiàn)代化計劃，采取的措施包括在導航信號L2上增加第二個民用碼（L2C）和增設(shè)第三個民用導航信號L5（1 176.45MHz）。通過使用所有三個頻率信號（L1C／A，L2C和L5）的載波相位觀測信息以及差分處理技術(shù)，可以快速獲得高精度的解算結(jié)果。主要針對GPS三頻觀測值的組合形式、特征參數(shù)以及應用前景進行研究。

1 載波相位雙差觀測方程

GPS精密定位一般采用雙差載波相位觀測值作為原始觀測數(shù)據(jù)，形成雙差觀測值的過程中可以極大削弱甚至消除許多空間相關(guān)及測站、衛(wèi)星相關(guān)的系統(tǒng)誤差，減少待估參數(shù)的個數(shù)，從而簡化觀測方程的表達形式。雙差載波相位觀測方程可以表示為

式中：▽Δ為雙差算子；φi為第i個頻率載波相位觀測值，單位：周（i＝1，2，3）；λi為第i個頻率載波波長，單位：m（i＝1，2，3）；R 為接收機與衛(wèi)星間幾何距離，單位：m；dR 為衛(wèi)星軌道誤差，單位：m；I1為L1頻率上的電離層時延，單位：m；T為對流層時延，單位：m；Ni為第i個頻率載波相位整周模糊度，單位：周；εΦi為第i個頻率載波相位觀測噪聲（包括接收機噪聲和多徑等），單位：m.

由式（1）可以看到，雙差觀測值中消除了衛(wèi)星鐘差和接收機鐘差的影響，殘余的誤差主要包括衛(wèi)星軌道誤差殘差、電離層時延殘差、對流層時延殘差以及包括接收機噪聲、多徑效應等在內(nèi)的觀測噪聲。多徑誤差和接收機噪聲在時間和空間上均不具備明顯的相關(guān)性，因此，無法通過差分方法進行有效的削弱和消除，在基線解算中通常在隨機模型中進行考慮［2]。

2 載波相位觀測值的線性組合

GPS現(xiàn)代化實現(xiàn)后，將成為擁有三個民用頻率信號的多頻衛(wèi)星導航系統(tǒng)。多個頻率信號的使用，使得載波相位觀測值的線性組合形式更加多樣，合理的選擇線性組合方式和工作順序，將有助于改善整周模糊度的解算成功率和效率［3]，提高粗差探測和周跳檢測與修復的精度和有效性［4－5]。下面詳細分析載波相位觀測值的線性組合形式及其性能，為基線數(shù)據(jù)處理策略的制定提供參考依據(jù)。

2.1 線性組合的一般形式

三頻載波相位觀測量線性組合的一般表達式可寫為

式中：φi（i＝1，2，3）表示第i個載波信號對應的相位觀測值；ki（i＝1，2，3）為線性組合系數(shù)；φLC為線性組合后的載波相位觀測量。φLC對應的整周模糊度NLC，可表示為

式中，Ni（i＝1，2，3）表示第i個載波信號對應的整周模糊度，為了保證NLC的整數(shù)特征，在實際應用中ki（i＝1，2，3）分別取為整數(shù)。

對于三頻雙差觀測方程（見式（1）），其線性組合觀測量可表示為

2.2 線性組合的誤差分析

1）幾何相關(guān)誤差

根據(jù)式（4），雙差觀測值線性組合后幾何相關(guān)誤差可表示為

2）電離層誤差

由式（4）可以看出，雙差觀測值線性組合后電離層誤差為

3）測量噪聲

為了簡化分析，假設(shè)不同頻率上的載波相位觀測誤差相互獨立，以米為單位的測量噪聲均方差等于相應波長的α倍，即σi＝λi·α，則根據(jù)誤差傳播定律，雙差觀測量的觀測噪聲為

線性組合觀測量的觀測噪聲表示為

2.3 電離層無關(guān)線性組合及性能分析

電離層誤差及其殘差是影響基線解算質(zhì)量的重要因素之一，而電離層無關(guān)線性組合可以消除電離層誤差的影響，為整周模糊度的快速、正確搜索提供了可能。電離層無關(guān)線性組合在實際應用中具有重要的意義。下面詳細討論電離層無關(guān)線性組合的形式及性能。

根據(jù)上述對載波相位觀測值線性組合誤差的分析，電離層無關(guān)線性組合應滿足以下條件

上述方程組屬于含不等式約束條件的非滿秩方程組?？紤]到λ1＜λ2＜λ3，現(xiàn)令k2、k3∈Z，k1∈R且限定取值范圍k2＝ ［－10，10]，k3＝ ［－10，10]，得到k1的一組浮點解序列，將k1取整得到其定點解序列「k1」，符號「·」表示對浮點解進行四舍五入運算。理論上，浮點解與定點解的偏差越小，則組合后觀測量包含的電離層誤差越小，偏差最小的五組線性組合及其特征值列于表1。

表1 滿足電離層無關(guān)條件的部分候選組合及其特征值

從表1中可以看到，上述各線性組合除（1，－6，5）外，組合后的波長均小于任意載波波長，較小的波長將對整周模糊度的解算產(chǎn)生不利影響。組合觀測值的電離層誤差（αI和βI）與載波L1相比降低1～2個數(shù)量級，在基線解算中可以忽略不計。與載波L1相比，組合觀測值的觀測噪聲（ασ和βσ）放大數(shù)倍，同時幾何相關(guān)誤差（αG）也有不同程度的增加。線性組合（1，－6，5）在波長和電離層誤差方面較其他組合具有明顯的優(yōu)勢，但以m為單位的觀測噪聲被放大2個數(shù)量級。在實際應用中，需要綜合考慮原始觀測數(shù)據(jù)的質(zhì)量、解算目標和使用的算法等因素，合理選擇線性組合形式。

2.4 幾種特殊的線性組合及性能分析

根據(jù)前面的介紹，通過將多個頻率觀測量進行合理的線性組合，可以構(gòu)造出許多具有長波長、低噪聲等良好特征的組合觀測量，有助于提高整周模糊度解算的成功率和定位精度。選擇不同的系數(shù)（k1，k2，k3），可以得到一系列具有不同特性的組合觀測量，其中比較重要的有超寬巷組合（EWL）、寬巷組合（WL）、中巷組合（ML）和窄巷組合（NL），詳細的參數(shù)列于表2。從表2中可以看到，超寬巷組合波長最長（586.1cm），具有較小的以周為單位的測量噪聲、幾何相關(guān)誤差以及電離層時延誤差，這些特性對于整周模糊度的解算極為有利，但是超寬巷組合觀測量以m為單位的測量噪聲較大，一般不宜直接應用于精密定位計算；寬巷組合和中巷組合波長分別為86.2cm和75.1cm，具有適中的電離層誤差和測量噪聲，有利于整周模糊度的確定，同時與超寬巷組合相比，以m為單位的測量噪聲顯著降低，可以提供更精確的定位結(jié)果；窄巷組合的波長僅為10.7cm，且?guī)缀蜗嚓P(guān)誤差和電離層誤差較大，不利于整周模糊度的確定，但較小的觀測噪聲使其在精密定位解算方面具有明顯的優(yōu)勢；（1，0，0）組合相當于載波 L1單頻觀測，通常使用載波L1觀測量可以得到較好的解算結(jié)果，在實際應用中載波L1的整周模糊度往往也是最關(guān)心的參數(shù)之一。

表2 幾種特殊線性組合及其特征值

分析線性組合序列EWL／WL／ML的特征值的變化可以看到：序貫組合觀測值的波長由5.86 m依次遞減直0.75m，以m為單位的測量噪聲放大系數(shù)由43.56依次遞減至5.58，而以周為單位的測量噪聲放大系數(shù)則保持不變，電離層時延誤差雖然有所增加但幅度有限，這說明采用上述序貫序列組合，在測量噪聲水平一定的條件下，可以逐步提高模糊度解算的成功率和基線解算的精度。

3 線性組合在解算中的應用

GPS現(xiàn)代化的實施，使民用信號的數(shù)量和質(zhì)量得到了不同程度的提高，將不同頻率觀測量進行合理組合可以得到許多具有低噪聲、長波長、電離層無關(guān)和幾何無關(guān)等良好特性的組合觀測量，為提高整周模糊度搜索效率和基線解算質(zhì)量創(chuàng)造了有利條件。如：在電離層誤差修正方面，以前的單頻和雙頻系統(tǒng)只能通過模型或差分方法消除一階電離層時延的影響，誤差改正精度有限，三頻觀測數(shù)據(jù)的組合應用可以將電離層時延消除至高階項，有效分離電離層誤差與對流層誤差、多徑效應等的影響［6]。借助三頻組合觀測量的優(yōu)良特性，可以將觀測數(shù)據(jù)中的粗差和周跳進行放大，更好地進行檢測和修正，提高觀測數(shù)據(jù)質(zhì)量［4]。在整周模糊度解算過程中，還可以充分利用不同組合觀測量各自的特點，在不同的解算階段采用不同的組合形式，進一步提高模糊度搜索的效率和成功率，使GPS在長基線和高動態(tài)實時等應用領(lǐng)域的應用成為可能，類似的方法有 TCAR、CIR、MCAR［7－10]等。

4 結(jié) 論

隨著GPS現(xiàn)代化的逐步實現(xiàn)，對于三頻觀測量的組合形式及其特性的研究逐漸成為GPS應用領(lǐng)域的研究熱點。將多頻觀測數(shù)據(jù)進行合理的組合可以得到許多具有良好特性的組合觀測量，如低噪聲、長波長、電離層無關(guān)和幾何無關(guān)等，這些特征對于特定的應用環(huán)境和系統(tǒng)具有重要的意義，可以有效提高電離層誤差改正、周跳檢測與修正、整周模糊度解算、精密定位等作業(yè)質(zhì)量，使GPS在長基線解算、高動態(tài)實時定位以及RTK測量領(lǐng)域的應用成為可能。

［1]寇艷宏.GPS原理與應用［M].2版.北京：電子工業(yè)出版社，2008.

［2]SATIRAPOD C，WANG Jin－ling，RIZOS C.Simplified MINQUE procedures for the estimation of variance－covariances components of GPS observables［J].Surv.Review，2002，36（286）：582－590.

［3]范建軍.GNSS三頻精密測量定位系統(tǒng)的數(shù)據(jù)處理理論及方法研究［D].長沙：國防科學技術(shù)大學，2007.

［4]王貴文.基于三差觀測量的實時動態(tài)GPS周跳修復方法研究［J].武漢大學學報，2007，32（8）：711－714.

［5]劉旭春.GPS三頻數(shù)據(jù)在周跳和粗差探測與修復中的應用［J].煤炭學報，2006，31（5）：585－588.

［6]伍 岳.GPS現(xiàn)代化后電離層折射誤差高階項的三頻改正方法［J].武漢大學學報，2005，30（7）：601－603.

［7]FORSSELL B，MARTIN NEIRA M，HARRIS R A，et al.Carrier phase ambiguity resolution in GNSS－2［C]／／Proceeding of ION GPS－97，Kansas ，USA，1997：1727－1736.

［8]JUNG J，ENGE P，PERVAN B.Optimization of cascade integer resolution with three civil GPS frequencies［C]／／Proceeding of the ION GPS－00，Salt Lake City，2000：2191－2201.

［9]JUNG J.High integrity carrier phase navigation for future LAAS using multiple civilian GPS signals［C]／／Proceding of ION GPS－99，Nashville，1999：727－736.

［10]WOLFGANG W，JON W.TCAR and MCAR options with Galileo and GPS［C]／／Proceeding of the ION GPS／GNSS－03，Portland，2003.