改進LAMBDA算法的雙頻BDS整周模糊度快速固定

龔宵雪，徐愛功，祝會忠，李 博，馬天明，高 猛

(遼寧工程技術大學 測繪與地理科學學院，遼寧 阜新 123000)

改進LAMBDA算法的雙頻BDS整周模糊度快速固定

龔宵雪，徐愛功，祝會忠，李 博，馬天明，高 猛

(遼寧工程技術大學 測繪與地理科學學院，遼寧 阜新 123000)

針對雙頻BDS短基線解算法方程病態嚴重，模糊度的浮點解精度低，導致LAMBDA算法搜索范圍過大、搜索時間過長等問題，提出一種改進的LAMBDA算法：運用Tikhonov正則化理論構造一種具有明確物理意義的正則化矩陣R，以獲得精確的整周模糊度浮點解；同時用均方誤差矩陣RMSEM代替原算法的協方差陣確定模糊度的搜索范圍，縮短初始化時間，提高搜索效率。實驗結果表明：改進后的LAMBDA算法獲得的模糊度浮點解精度顯著提高，且無需初始化時間，能夠正確固定模糊度的整數解，快速實現cm級定位。

BDS；整周模糊度；LAMBDA算法；正則化；浮點解

0 引言

北斗衛星導航系統(Beidou navigation satellite system，BDS)是我國具有獨立自主產權、自主研發的衛星導航定位系統[1]。2012-12-27，BDS空間信號接口控制文件正式版1.0正式公布，BDS已開始向中國及亞太大部分地區正式提供無源定位、導航、授時服務[2]。截至2016-06-13，BDS在軌衛星共計23顆，其中7顆地球靜止軌道(geostationary Earth orbit，GEO)衛星、8顆傾斜地球同步軌道(inclined geo-synchronous orbits，IGSO)衛星及8顆中圓地球軌道(medium Earth orbit，MEO)衛星的定位精度為10 m、測速精度為0.2 m/s、授時精度為10 ns[3]。到2020年左右將全面建成BDS，形成全球服務能力。

BDS現在也已經加入到高精度定位中。在高精度定位中，通常用載波相位差分技術進行快速定位，而載波相位整周模糊度的準確固定是BDS定位的核心問題[4-5]，整周模糊度一旦確定，相位觀測值即轉變為準確的距離觀測值[6]。關于如何更好地確定模糊度，已有許多學者提出了很多方法。文獻[7]針對病態法矩陣，提出基于Tikhonov正則化原理的一種選擇正則化矩陣R的新方法，減弱了法方程的病態性，但沒有明確R的物理意義；文獻[8]將直接算法(direct calculation，DC)推廣為直接解算短基線全球定位系統(global positioning system，GPS)整周模糊度，但DC算法需要保證基線長度的精度優于0.7 m；文獻[9]利用單歷元測碼偽距觀測值和雙頻載波相位觀測值組成雙差觀測方程，解算L1/L2雙差整周模糊度時具有較高的效率和成功率，利用方差對模糊度進行分組降維處理，縮小了模糊度搜索空間；文獻[10]指出模糊度去相關處理(Z變換)有迭代法和聯合去相關法2種算法，通過比較2種方法的條件數和去相關數，得出聯合去相關法的效果略高于迭代法，在整周模糊度搜索前對原始模糊度進行降相關處理，使搜索橢球更接近球體，提高了搜索效率，但未對模糊度浮點解精度做出說明；文獻[11]提出了一種利用多頻觀測值，基于有幾何模型的逐級模糊度固定方法，先通過偽距固定超寬巷模糊度，然后用無模糊度的超寬巷約束2個寬巷模糊度求解，最后使用無模糊度的寬巷觀測值約束窄巷模糊度求解，其模糊度固定達到100 %的成功率。

對于如何快速又精確地固定模糊度的問題，本文通過運用Tikhonov正則化的理論，構造出一種具有明確物理意義的正則化矩陣R，使LAMBDA算法得到改進，以解決法方程的病態性，獲得精確的整周模糊度浮點解；同時用均方誤差矩陣(mean square error matrix，MSEM)RMSEM代替原算法的協方差陣確定模糊度的搜索范圍，縮短初始化時間，提高搜索效率。

1 LAMBDA算法原理

LAMBDA算法的核心問題就是整周模糊度最小二乘估計，主要分為2個組成部分：整數變換和基于喬里斯基分解(LTDL)的整周搜索。整數最小二乘估計的基本方法為

(1)

1.1 整數變換

(2)

(3)

1.2 整周搜索

(4)

求得最優的整數組合z之后，再進行逆變換，即

(5)

式中Z為整數變換矩陣。變換后的參數N滿足式(1)，也就是說，逆變換后得到的參數N就是要尋求的最佳整數模糊度向量。

在LAMBDA法中，是通過LTDL分解對協因數陣QX進行整數變換來提高搜索效率的，其中L為下三角矩陣，D為對角陣。

2 基線解算的一般LAMBDA算法

LAMBDA算法主要是對載波相位雙差觀測方程中的模糊度進行搜索解算，在短基線數據處理時，雙差可以消除衛星鐘差和接收機鐘差，同時也大幅度削弱了衛星軌道誤差、電離層延遲和對流層延遲等誤差。將BDS的B1、B2雙差載波相位觀測值組成寬巷載波相位觀測值，由于B1、B2載波相位寬巷觀測值的波長達84.7 cm，因而很容易準確確定其整周模糊度，然后解算載波相位觀測值的整周模糊度[13]。

設2臺接收機在某歷元可共同觀測k+1顆衛星，則可組成k個相位雙差觀測方程，其線性化后的雙差觀測方程[14]簡寫為

(6)

若2臺GPS接收機連續觀測n1個歷元，則總的雙差觀測方程為

(7)

可以將式(7)簡寫成

(8)

根據最小二乘(least square method，LS)原則，由式(8)組成法方程，得到LS解為

(9)

式中N0=ATPA， 為LS估計的法矩陣。Y的協因數陣為

(10)

在定位過程中，如果觀測歷元數n1較少時，法矩陣N0的條件數一般在106以上，法方程病態嚴重，模糊度的浮點解精度低，導致LAMBDA算法搜索范圍過大、搜索時間過長，甚至不能正確固定整周模糊度。為了彌補LAMBDA算法的不足，本文對LAMBDA算法進行改進，使整周模糊度可以快速、高效地固定。

3 改進的LAMBDA算法

根據Tikhonov正則化原理，求解式(8)采用的估計準則[7]為

(11)

式中：α是平滑因子或正則化參數；R是正則化矩陣，用于減弱法方程病態；Ω(Y)是穩定泛函；‖?‖表示歐氏2-范數。

(12)

式中：R是奇異矩陣；R1是基線分量部分的先驗協因數陣的逆陣。系數1/100 000是根據待定點近似坐標偏差與相位觀測值精度比選取的。

構造正則化矩陣R后，應用L曲線法進行了大量的實際計算，其計算結果表明選取正則化參數α=1時效果最好。則式(11)變為

(13)

結合式(8)，由估計準則求導組成法方程為

(14)

式中R為正則化矩陣。解算式(14)的法方程，得

(15)

相應的均方誤差矩陣RMSEM近似[14]取為

(16)

式中σ0為單位權中誤差。

4 實驗及結果分析

為了驗證改進后的LAMBDA算法的效果[15]，采用BDS/GPS雙系統雙頻接收機對長度為50 m的基線進行了40 min的數據采集，采樣間隔為1 s，衛星截止高度角為15°。實驗對基線數據按照以下2種方案進行解算：1)傳統的一般算法固定整周模糊度；2)通過正則化改進的LAMBDA算法固定整周模糊度。

本文利用Ratio值作為模糊度得到固定解的置信度指標，Ratio值是由次優的模糊度向量的殘差平方和與最優的模糊度向量的殘差平方和的比值計算得到的。本文取Ratio值為2，當Ratio值大于2時，認為模糊度固定正確。利用改進前和改進后2種LAMBDA算法固定整周模糊度的Ratio值，其結果如圖1所示。

由圖1可知，改進后Ratio值小于2的個數遠遠少于改進前的個數。改進后只有1個歷元Ratio值小于2，改進前有498個歷元Ratio值小于2；所以利用改進后的LAMBDA算法固定的整周模糊度是可靠的。

以1號衛星和2號衛星為例，按照2種算法確定的整周模糊度的浮點解的對比結果如圖2所示。因為改進后的LAMBDA算法提高了模糊度浮點解精度，縮小了LAMBDA算法的搜索范圍和搜索時間，使模糊度的固定成功率得到提高；所以在得到模糊度固定解之前，提高模糊度浮點解非常有必要。

由圖2可知，1號衛星改進前的LAMBDA算法需要大約200個歷元才能確定模糊度浮點解，而改進后的LAMBDA算法第1個歷元就可以把模糊度浮點解確定；2號衛星在改進算法前也是大約需要200個歷元才能確定浮點解，同樣改進后的算法在第1個歷元就能確定浮點解。可以得出，改進后的LAMBDA算法明顯提高了模糊度浮點解的精度，使模糊度的固定成功率得到提高。

仍以1號和2號衛星為例，給出LAMBDA算法改進前后模糊度固定的結果，如圖3所示。

由圖3可知，改進前的LAMBDA算法的模糊度浮點解在開始時偏差大，所以整數解沒有正確固定，經過大約200個歷元以后才能固定整數解，而改進后的LAMBDA算法從第1個歷元開始就能固定到正確的整數解。

下面對改進前和改進后的LAMBDA算法初始化時間和模糊度固定個數進行統計，結果如表1所示。

表1 初始化時間和模糊度固定個數

根據表1的統計結果可以看出：改進前的LAMBDA算法固定了2 208個歷元，而改進后的LAMBDA算法固定了2 399個歷元，100 %固定。

在模糊度固定之后，用改進前和改進后的算法解算高精度定位結果，與測站的精確定位結果進行比較，得到東(E)、北(N)、天頂(U)3個方向上的偏差值，如圖4所示。

根據圖4定位結果，再結合前面的Ratio值、模糊度浮點解和固定解的結果可以看出，改進前的LAMBDA算法需要192個歷元才能達到cm級定位精度，而改進后的LAMBDA算法從第1個歷元開始就能達到cm級定位，且坐標偏差值變化范圍在3 cm之內，與整體解算的結果相比較波動不大，說明改進后LAMBDA算法解的精度更好，定位效率更高。

5 結束語

本文研究了雙頻BDS短基線解算問題，首先利用寬巷確定其整周模糊度，然后解算載波相位觀測值的整周模糊度。針對法方程病態嚴重，模糊度的浮動解精度低，導致LAMBDA算法搜索范圍過大、搜索時間過長問題，運用Tikhonov正則化的理論，構造了一種具有明確物理意義的正則化矩陣R，使LAMBDA算法得到改進，從而獲得精確的整周模糊度浮點解，同時用RMSEM代替原算法的協方差陣確定模糊度的搜索范圍，提高搜索效率。改進后的LAMBDA算法所獲得的模糊度浮點解精度顯著提高，可以無需初始化時間，100 %正確固定模糊度的整數解，快速實現cm級定位。

[1] 楊元喜.北斗衛星導航系統的進展、貢獻與挑戰[J].測繪學報，2010，39(1)：1-6.

[2] 高猛，徐愛功，祝會忠.雙導航定位系統偽距單點定位數據處理方法與精度分析[J].導航定位學報，2014，2(2)：83-88.

[3] 施闖，趙齊樂，李敏，等.北斗衛星導航系統的精密定軌與定位研究[J].中國科學：地球科學，2012，42(6)：854-861.

[4] 盧立果，劉萬科，李江衛.降相關對模糊度解算中搜索效率的影響分析[J].測繪學報，2015，44(5)：481-487.

[5] 李豹，許江寧，曹可勁，等.改進LAMBDA算法實現單頻GPS整周模糊度快速解算[J].中國慣性技術學報，2013，21(3)：365-368.

[6] 魏子卿，葛茂榮.GPS相對定位的數學模型[M].北京：測繪出版社，1988：105.

[7] 王振杰，歐吉坤，柳林濤.單頻GPS快速定位中病態問題的解法研究[J].測繪學報，2005，34(3)：196-201.

[8] 邱蕾，花向紅，蔡華，等.GPS短基線整周模糊度的直接解法[J].武漢大學學報(信息科學版)，2009，34(1)：97-99，104.

[9] 祝會忠，高星偉，秘金鐘，等.一種GPS整周模糊度單歷元解算方法[J].測繪科學，2011，36(4)：9-11.

[10]劉鋒，蘭孝奇，陳飛.基于LAMBDA整周模糊度解算方法中的整周Z變換算法[J].測繪與空間地理信息，2014，37(4)：51-53.

[11]劉炎炎，葉世榕，江鵬，等.基于北斗三頻的短基線單歷元模糊度固定[J].武漢大學學報(信息科學版)，2015，40(2)：209-213.

[12]李征航，黃勁松.GPS測量與數據處理[M].武漢：武漢大學出版社，2012：155-158.

[13]祝會忠，徐愛功，高猛.北斗系統載波相位動態差分定位方法[J].測繪科學，2015，40(9)：9-13.

[14]盧秀山，王振杰，歐吉坤，等.單頻GPS短基線快速定位中的少數歷元算法[J].武漢大學學報(信息科學版)，2007，32(12)：1147-1151.

[15]史仲，吳任洪，王以磊.提高GPS基線解算質量的方法研究[J].測繪通報，2015(4)：69-71.

Rapid resolution of dual frequency BDS integer ambiguity by improved LAMBDA algorithm

GONGXiaoxue，XUAigong，ZHUHuizhong，LIBo，MATianming，GAOMeng

(School of Geomatics，Liaoning Technical University，Fuxin，Liaoning 123000，China)

Aiming at the problem that the solution algorithm of the BDS dual frequency short baseline has ill conditioned equation，and the floating point solution of ambiguity has low precision，leading to the large searching scope and long searching time in using LAMBDA algorithm，the paper proposed an improved LAMBDA algorithm：Tikhonov regularization theory was used to construct a regularized matrixRwith a clear physical meaning for obtaining the accurate floating point solution of integer ambiguity；meanwhile a mean square error matrixRMSEMwas used to replace the covariance matrix of the original algorithm in order to determine the searching range of ambiguity，thus the initialization time was reduced and the searching efficiency was improved.Experimental result showed that by the proposed method，the precision of the floating point solution of ambiguity could be increased significantly without initialization time，the integer ambiguity could be fixed correctly and the positioning with the centimeter level accuracy could be rapidly realized.

BDS；integer ambiguity；LAMBDA algorithm；regularization；floating point solution

2016-10-18

國家重點研發計劃項目(2016YFC0803102)；遼寧省高等學校創新團隊項目(LT2015013)。

龔宵雪(1991—)，女，河南安陽人，碩士研究生，研究方向為GNSS高精度定位。

龔宵雪，徐愛功，祝會忠，等.改進LAMBDA算法的雙頻BDS整周模糊度快速固定[J].導航定位學報，2017，5(3)：72-76,83.(GONG Xiaoxue，XU Aigong，ZHU Huizhong，et al.Rapid resolution of dual frequency BDS integer ambiguity by improved LAMBDA algorithm[J].Journal of Navigation and Positioning，2017，5(3)：72-76,83.)

10.16547/j.cnki.10-1096.20170315.

P228

A

2095-4999(2017)03-0072-06