論精密單點定位整周模糊度解算的不同策略

張寶成，歐吉坤

1.中國科學院測量與地球物理研究所 動力大地測量學重點實驗室，湖北 武漢430077；2.中國科學院 研究生院，北京100049

論精密單點定位整周模糊度解算的不同策略

張寶成1，2，歐吉坤1

1.中國科學院測量與地球物理研究所 動力大地測量學重點實驗室，湖北 武漢430077；2.中國科學院 研究生院，北京100049

由于GPS非差相位觀測值的相位偏差（initial phase biases，IPB）與整周模糊度難以分離，精密單點定位（PPP）估值均為模糊度浮點解。首先對GPS原始觀測方程的秩虧問題進行分析，從參數整合的角度，推導衛星IPB估計的滿秩函數模型，形成一種新的PPP－AR算法。以此為基礎，對已有兩種算法的特點進行對比分析。研究表明，分解法是一種觀測信息的最優利用，且與傳統的星鐘估計方法具有一致性，但未發掘衛星IPB較為穩定的有利約束；非整法對所采用的組合觀測值之間的相關性未加考慮，是一種次優估計，實時性較差，且較依賴于高精度的碼觀測。推導的新算法可有效克服已有算法的不足，便于施加部分參數的合理時變性約束，提高衛星IPB估計的可靠性。

精密單點定位；初始相位偏差；PPP模糊度固定；分解衛星鐘差；非整相位偏差

1 引 言

自1997年Zumberge等人提出PPP概念以來［1］，該技術已廣泛地應用于區域或全球坐標框架維持、高精度導航定位、精密授時以及大氣延遲提取等領域［2］。與聯合多測站GPS觀測數據實施網平差相比，PPP在保證解算精度的同時，實施靈活方便且能有效減少計算負擔。此外，PPP基于非差GPS觀測值，與基于雙差觀測值的數據處理策略相比，其數據利用率更高，計算結果更為可靠［3］。

PPP目前的發展主要受到兩個因素的制約：其一，PPP一般基于IGS等機構發布的衛星軌道和鐘差產品，其實時性實施受到一定程度的影響［4］；其二，非差相位觀測值中的IPBs與整周模糊度難以分離，使得PPP的各類估值均為模糊度浮點解。若能有效利用非差模糊度的整數特性，PPP的估值精度和收斂特性將有進一步改善的空間［5－12］。

當前，伴隨通信技術的發展，基于區域參考站網的GPS觀測數據實時估計、播發以及接收衛星軌道和鐘差的研究已取得較大進展，在充分考慮數據傳輸和處理的時延后，實時PPP的定位精度可達到分米甚至厘米級［4］，有效地滿足了大部分實時應用的需求。

與此同時，PPP模糊度固定（PPP ambiguity resolution，PPP－AR）的研究已經起步，且形成了兩種不同的處理策略：估計分解衛星鐘差（分解法）［6，11－14］或非整相位偏差（非整法）［5，7－10］算法。上述兩種算法均采用區域或全球范圍的GPS網，實時或事后分離衛星IPBs，用于改正PPP相位觀測值，還原非差模糊度的整數特性。與標準PPP相比，PPP－AR可有效改善測站坐標解的東分量精度［7，12］。對于遠離地面基準站的海上動態平臺定位或低軌衛星定軌而言，上述PPP－AR算法還可有效代替目前常用的長基線解算策略［5，9］。

由此可知，為實現PPP－AR，關鍵在于如何利用GPS網估計衛星IPBs，不同網處理策略構成了現有兩種PPP－AR方法的主要差異。此外，第一種方法采用消電離層組合的碼、相位觀測值以及MW組合觀測值，第二種方法對之進一步實施星間單差［6－10，13－14］。上述組合和差分過程“簡單”地消除了未知參數個數，但卻可能模糊對各類線性相關參數所引起的秩虧現象的認識。

為給出完整的PPP－AR概念，筆者從GPS原始觀測方程出發，研究它們的秩虧特性，通過選取適當的基準參數，依次消除了各類秩虧，構建了衛星IPBs估計的滿秩函數模型，形成一種PPP－AR新算法。筆者還分別導出了現有兩種方法各自的數學模型和理論假設，并從觀測數據利用率、實時性以及用戶實施等角度進行了對比分析。

2 衛星IPBs估計

本節推導利用參考站網GPS原始觀測值估計衛星IPBs的函數模型，針對模型中的秩虧問題，采用參數整合的消秩虧策略。

2.1 秩虧函數模型

忽略觀測噪聲和多路徑等誤差，且假定測站和衛星位置已知，則GPS原始觀測方程可表示為［15］

式中，λj表示頻率j相位觀測值的波長。

假定歷元i，n個接收機共觀測到m顆共視衛星，以雙頻觀測數據為例，所有線性化方程可聯合表示為

式中

式中，k＝1，…，n、q＝1，…，m和j＝1，2分別代表測站、衛星和頻率數；式（4）中的Pqk，j（i）和Φqk，j（i）為所有線性化的碼和相位觀測值；式（5）中的子矩陣（從左到右）分別對應于式（6）中的各類參數，含義與式（1）相同（省略了歷元符號i，以簡化表達）。另

式中，?、Im和en分別表示克羅內克積［18］、m維單位矩陣以及各元素均為1的n維列向量；Λ為二維對角陣，其對角元為雙頻GPS觀測值的波長

分析A（i）中后四個子矩陣可知，鐘差、CIDs、IPBs、模糊度以及電離層延遲之間存在線性相關性，導致A（i）存在列秩虧。

2.2 消秩虧策略

秩虧導致參數估計結果非唯一。為消除該秩虧對參數估計的不利影響，提出如下逐次消秩虧的策略：采用參數整合的方法，將式（6）中的部分參數選取為基準［16］，且遵循如下準則：① 基準參數的個數須與秩虧數相等；②衛星IPBs參數不宜被選取為基準，以確保其獨立可估性；③ 最終可估的模糊度須保持其整數特性。具體過程如下：

首先，分析式（8）的D矩陣可知，該矩陣的秩虧源于接收機和衛星的鐘差、CIDs以及IPBs之間的線性相關，將某測站（即基準站，此處假定為第一個參考站）的接收機鐘差dt1、雙頻CIDs b1，j和IPBsφ1，j選取為基準。其余測站和衛星的相應參數可表示為

式中，p＝2，…，n。此時，A（i）中對應的第二和第三個子矩陣可表示為

式中

其次，考察式（11）可知，該矩陣的秩虧數為n＋m－1，且源于CIDs與鐘差之間的線性相關性。將第一個頻率的CIDs（式（10）中的和）選取為基準，可有效消除該秩虧，得到滿秩矩陣和待估參數可分別表示為

第三，單獨考察第j個頻率的衛星IPBs和模糊度的設計矩陣

該矩陣秩虧數為m。解決的辦法是將基準站觀測到的m顆衛星的模糊度Mq1，j作為基準，得到的滿秩矩陣和可估參數分別為

由式（17）可知，經由此步消秩虧，非差整周模糊度被參數化為站間單差的形式。隨后，考察第j個頻率的接收機IPBs和該單差模糊度的設計矩陣

該矩陣秩虧數為n－1?？赏ㄟ^選取某衛星（即基準星，此處假定為第一顆衛星）至所有參考站（除基準站）的單差模糊度作為基準，得到滿秩矩陣和可估參數分別為

式中，t＝2，…，s。至此，站間單差的模糊度被參數化為雙差的形式，其整數特性可在參數估計過程中加以合理的考慮。

最后，經過上述消秩虧過程后，鐘差、CIDs、IPBs以及電離層延遲參數的設計矩陣可聯合表示為

式（21）對應矩陣的秩虧數為n＋m－1，可將第二個頻率的CIDs（Δbp，2和Δbq2）選取為基準，以消除該秩虧。

至此，通過上述的逐項分析，選取適當的基準參數，經參數整合，原矩陣A（i）中的列秩虧可被充分消除，得到的滿秩設計矩陣及部分待估參數分別為

在實際數據處理中，基于上述推導的觀測方程，先固定雙差模糊度，隨后回代，即可求解得到衛星的IPBs和鐘差等參數的模糊度固定解，再播發給PPP用戶使用。

3 現有PPP－AR策略分析

由上節的滿秩函數模型出發，本節推導了兩種現有PPP－AR方法的模型，并分別從觀測信息利用率、實時性以及用戶實施等方面進行了對比分析。

3.1 分解法

分解法［6，11－14］采用消電離層組合的碼和相位觀測值，在觀測域事先消除電離層延遲的影響。通過引入MW組合觀測值，將消電離層組合模糊度轉換為寬巷和窄巷模糊度，同時有效保留原始觀測值中的信息。上述三類組合觀測值與原始碼、相位觀測值之間的轉化關系為

事實上，將列滿秩的轉換矩陣Γ?Inm（rk（Γ?Inm）＝rk（Γ）rk（Inm）＝cl（Γ?Inm）＝3 nm，其中rk（）和cl（）分別表示矩陣的秩和列數［18］）乘以式（22）中的Af（i），且合并雙頻的IPBs和模糊度參數，可得到式（24）中組合觀測值所對應的設計矩陣以及可估參數。具體形式分別為

式中，Ad（i）中不再包含電離層延遲的設計矩陣，其余子矩陣與Af（i）中的對應部分含義相同；和分別為中的站星IPBs；和為的站星CIDs；和分別為和的模糊度參數，其波長因子為其余參數的含義同式（23）。站星鐘差和天頂對流層延遲等參數的形式未受轉換矩陣的影響，此處不再列出。

傳統的衛星鐘差估計算法中，碼和相位觀測方程中包含的可估鐘差參數形式一致［4］，該參數化形式可導致鐘差估值易受碼觀測值的影響，進而出現日間不連續等異?，F象［19］。在分解法的參數估計過程中，碼和相位分別對應不同類型的鐘差參數，具體為將式（26）中的和與對應于的站星鐘差合并，且模型化為具有白噪聲特性的時變參數，由此得到對應于的新站星鐘差分別為

至此，分解法函數模型的設計矩陣可表示為

基于上述推導得到的函數模型，分解法的實施特點可概括如下。

該算法基于上述三類組合觀測值（式（24））實施參數估計，在有效消除電離層延遲的同時，充分利用了原始觀測值中的有效信息。在觀測值加權策略中，若合理考慮由于觀測值組合所引起的數學相關性，則各類參數估計結果與基于原始觀測值等價。

基于卡爾曼濾波算法，分解法可實現各類參數的實時遞歸估計，模糊度參數的整數特性可在參數估計過程中被合理利用，以有效提高其余各類參數的精度和可靠性。

對于用戶而言，PPP－AR的實施與標準PPP類似，區別僅在于消電離層組合的碼、相觀測值需采用分解法計算得到的分別改正（標準PPP算法中，碼、相觀測值均采用改正），同時用戶可選擇MW組合觀測值并施加改正以輔助整周模糊度的快速固定。這意味著目前常用的標準PPP軟件，經過簡單的升級即可滿足模糊度固定的要求，故實施過程將比較經濟、有效。

3.2 非整法

非整法［5，7－10］基于星間單差的和三類觀測值。將差分矩陣I3n?［－em－1Im－1］乘以式（29），可得非整法函數模型的設計矩陣

式中

非整法的實施基于如下的分步數據處理策略。

類似于第一步，將連續弧段內的式（34）經由取平均和取整運算，分離得到（窄巷IPBs）的小數部分至此，可將得到的連同播發給PPP用戶以實現模糊度解算。

與分解法相比，非整法的實施特點可簡要概括如下。

非整法對分解法中所采用的三類組合觀測值進一步實施星間差分運算，以消除與測站有關的偏差影響，進而減少了待估參數的個數。然而，由此所引起的各類觀測值之間的數學相關性在非整法的分步實施過程中未合理考慮：如計算或均基于逐衛星求解的方法，各衛星單差觀測值之間的相關性一般被忽略，導致上述兩類估值的精度評定復雜。

此外，非整法取平均過程使得實時參數估計難以實施，不僅如此，簡單的取平均運算并不能有效消除部分非模型化誤差的影響（如碼觀測值的多路徑等），尤其是當連續弧段的時長較短時，對于較為依賴碼觀測值的求解過程而言，估值結果往往并不可靠。

最后，由于非整法提供給PPP用戶的產品均為星間單差的形式，故用戶需要采用相同的星間單差觀測值，這與目前標準PPP算法中基于非差觀測值的數據處理策略差別較大，用戶的軟件升級所需成本等代價較之分解法更高。

4 結論與建議

本文推導了利用參考站網GPS原始觀測值估計衛星IPBs的滿秩觀測方程。通過卡爾曼濾波，可估計得到衛星IPBs，為PPP－AR提供一種新算法。作為對比，分別導出了目前常用的兩種PPP－AR策略：分解法和非整法的數學模型，并對兩種方法的實施特點進行了對比性的論述。分析表明，分解法的參數估計過程充分利用了觀測值中的有效信息，便于用戶實時或事后實施PPP模糊度解算；非整法采用逐衛星取平均的策略估計星間單差的IPBs，其實時性的實施較為困難，且用戶需采用相同的星間單差運算，增加了可能的計算負擔。

上述兩類方法的本質區別還在于站星IPBs的時變模型選擇：分解法認為兩類參數隨時間變化的特性在參數估計過程中難以有效描述，故該方法將上述兩類偏差的消電離層組合項與對應的站星鐘差合并，且模型化為具有白噪聲性質的參數；與之相反，非整法則認為兩類參數均隨時間緩慢變化，且其隨時間變化的部分可通過簡單的取平均加以消除。

與上述兩種基于組合或差分觀測值的算法不同，本文基于原始觀測值估計衛星IPBs的算法具有如下的優勢：觀測值的組合或差分所引起的數學相關性對于觀測值合理定權的不利影響可以較好地避免；可估電離層延遲的平穩時間變化可作為一種有效的約束以增強各類估值的精度與可靠性；站星IPBs的時變性可采用較為合理的模型加以描述（如隨機游走），以避免分解法中將其模型化為白噪聲所引起的待估參數增多，估值精度不高，同時也可避免非整法中將之過強地約束為時不變參數所引起的可靠性降低等問題。

［1］ ZUMBERGE J，HEFLIN M，JEFFERSON D，et al.Precise Point Positioning for the Efficient and Robust Analysis of GPS Data from Large Networks［J］.Journal of Geophysical Research，1997，102（B3）：5005－5017.

［2］ ZHANG Baocheng，OU Jikun，YUAN Yunbin，et al.Precise Point Positioning Algorithm Based on Original Dualfrequency GPS Code and Carrier－phase Observations and Its Application［J］.Acta Geodaetica et Cartographica Sinica，2010，39（5）：478－483.（張寶成，歐吉坤，袁運斌，等.基于GPS雙頻原始觀測值的精密單點定位算法及應用［J］.測繪學報，2010，39（5）：478－483.）

［3］ BLEWITT G.Carrier Phase Ambiguity Resolution for the Global Positioning System Applied to Geodetic Baselines up to 2000km［J］.Journal of Geophysical Research，1989，94（B8）：10187－10203.

［4］ ZHANG Xiaohong，Li Xinxin，Guo Fei，et al.Servicebased Real－time Precise Point Positioning and Its Application［J］.Chinese Journal of Geophysics，2010，53（6）：1308－1314.（張小紅，李星星，郭斐，等.基于服務系統的實時精密單點定位技術及應用研究［J］.地球物理學報，2010，53（6）：1308－1314.）

［5］ BERTIGER W，DESAI S D，HAINES B，et al.Single Receiver Phase Ambiguity Resolution with GPS Data［J］.Journal of Geodesy，2010，84（5）：327－337.

［6］ COLLINS P，BISNATH S，LAHAYE F，et al.Undifferenced GPS Ambiguity Resolution Using the Decoupled Clock Model and Ambiguity Datum Fixing［J］.Navigation，2010，57（2）：123－135.

［7］ GE M，GENDT G，ROTHACHER M，et al.Resolution of GPS Carrier－phase Ambiguities in Precise Point Positioning（PPP）with Daily Observations［J］.Journal of Geodesy， 2008，82（7）：389－399.

［8］ GENG J，TEFERLE F N，SHI C，et al.Ambiguity Resolution in Precise Point Positioning with Hourly Data［J］.GPS Solutions，2009，13（4）：263－270.

［9］ GENG J，TEFERLE F N，MENG X，et al.Kinematic Precise Point Positioning at Remote Marine Platforms［J］.GPS Solutions，2009，14（4）：343－350.

［10］ GENG J，MENG X，DODSON A H，et al.Integer Ambiguity Resolution in Precise Point Positioning：Method Comparison［J］.Journal of Geodesy，2010，84（9）：569－581.

［11］ LAURICHESSE D，MERCIER F，BERTHIAS JP，et al.Integer Ambiguity Resolution on Undifferenced GPS Phase Measurements and Its Application to PPP and Satellite Precise Orbit Determination［J］.Navigation，2009，56（2）：135－149.

［12］ ZHANG Xiaohong，LI Xinxin.A New Method for Zerodifferenced Ambiguity Resolution and Its Application to PPP［J］.Geomatics and Information Science of Wuhan University，2010，35（6）：657－660.（張小紅，李星星.非差模糊度固定解PPP新方法及實驗［J］.武漢大學學報：信息科學版，2010，35（6）：657－660.）

［13］ COLLINS P.Isolating and Estimating Undifferenced GPS Integer Ambiguities［C］∥Proceedings of ION National Technical Meeting.San Diego：［s.n.］，2008：720－732.

［14］ COLLINS J P.US 2010／0188285A1Decoupled Clock Model with Ambiguity Datum Fixing［P］.Ottawa：Minister of Natural Resources，2010.

［15］ TEUNISSEN P J G，KLEUBERG A.GPS for Geodesy［M］.2nd ed.Berlin：Springer，1998.188－194.

［16］ OU Jikun.Uniform Expression of Solutions of Ill－posed Problems in Surveying Adjustment and the Fitting Method by Selection of the Parameter Weights［J］.Acta Geodaetica et Cartographica Sinica，2004，33（4）：283－288.（歐吉坤.測量平差中不適定問題解的統一表達與選權擬合法［J］.測繪學報，2004，33（4）：283－288.）

［17］ SCHAFRIN B，GRAFAREND E.Generating Classes of Equivalent Linear Models by Nuisance Parameter Elimination，Applications to GPS Observations［J］.Manuscripta Geodetica，1986，11（3）：262－271.

［18］ AMIRI－SIMKOOEI A R.Least－square Variance Component Estimation：Theory and GPS Applications［D］.Delft：Delft University of Technology，2007.

［19］ LAHAYE F，COLLINS P，CERRETTO G，et al.Advances in Time and Frequency Transfer from Dual－frequency GPS Pseudorange and Carrier－phase Observations［C］∥Proceedings of 40th Annual Precise Time and Time Interval（PTTI）Meeting.Virginia：［s.n.］，2008：415－432.

On the Different Strategies for Integer Ambiguity Resolution in Precise Point Positioning

ZHANG Baocheng1，2，OU Jikun1
1.Key Laboratory of Dynamic Geodesy，Institute of Geodesy and Geophysics，Chinese Academy of Sciences，Wuhan 430077，China；2.Graduate University of Chinese Academy of Sciences，Beijing100049，China

Due to the inseparability between the initial phase biases（IPB）and the integer ambiguities in undifferenced（UD）GPS carrier phase observable，the solutions from precise point positioning（PPP）are always ambiguity－float.The rank－defect UD GPS observation equation is firstly investigated，and the full－rank mathematical model for satellite IPB estimation is derived by means of re－parameterization，thus forming a new strategy for PPPAR.Comparative analysis of the characteristics of both existing methods is subsequently conducted based on the derived re－parameterized equations.It is concluded from the analysis that full utilization of information from GPS observations is guaranteed in DSC method.The process of generating DSC is consistent with the standard satellite clocks estimation.However，the stable temporal behaviors of satellite IPB are not reasonably exploited.In contrast，the FPB method doesn’t fully account for the mathematical correlations between the adopted linear combinations of GPS observations，which would result in a sub－optimal satellite IPB.Additionally，the real－time implementation of this method is more troublesome than the DSC，and the GPS code observations with high accuracy are always required.The newly presented strategy here can avoid the shortcomings in both current methods，and reasonable constraints that are converted from the temporal stability of several unknowns can be easily imposed upon the parameter estimation to lead to more reliable satellite IPB.

precise point positioning（PPP）；initial phase biases（IPB）；PPP ambiguity resolution（PPP－AR）；decoupled satellite clocks（DSC）；fractional phase biases（FPB）

ZHANG Baocheng（1985—），male，PhD candidate，majors in study of applications of precise point positioning and network－based RTK.

1001－1595（2011）06－0710－07

P228

A

國家自然科學基金（40874009；41074013）；國家杰出青年科學基金（40625013）；中國科學院研究生科技創新與社會實踐資助專項

叢樹平）

2010－10－11

2011－05－16

張寶成（1985—），男，博士生，研究方向為精密單點定位與網絡RTK技術的應用。

E－mail：b.zhang＠whigg.ac.cn