一種GNSS測碼偽距絕對定位求解方法*

常占強 ，劉家喜 ，錢淑君 ，張妍欣 ，鄭浩鑫

（1.首都師范大學(xué)資源環(huán)境與旅游學(xué)院，北京 100048；2.三維信息獲取與應(yīng)用教育部重點實驗室，北京 100048；3.山東省第四地質(zhì)礦產(chǎn)勘查院，山東 濰坊 261021）

0 引 言

全球衛(wèi)星導(dǎo)航系統(tǒng)（global navigation satellite system，GNSS）泛指中國的北斗（Beidu）、美國的GPS（globalpositioning system）、俄 羅 斯 的GLONASS（global navigation satellite system）和歐洲的伽利略（Galileo satellite navigation system）.現(xiàn)今，GNSS定位導(dǎo)航功能不僅廣泛應(yīng)用于軍事、測繪、公安交通、數(shù)字測圖、精細(xì)農(nóng)林業(yè)及土地資源調(diào)查與勘測等領(lǐng)域［1-9］，并已滲透到人類日常生活中，如旅游、野外探險等.其中，GNSS測碼偽距絕對定位僅用一臺接收機不需要進行差分處理就可實現(xiàn)全球快速定位導(dǎo)航，在軍事、公安交通、車輛艦船定位導(dǎo)航、精細(xì)農(nóng)林業(yè)、土地資源調(diào)查與勘測等領(lǐng)域中發(fā)揮著重要作用.尤其是精確單點定位技術(shù)（pre‐cise point positioning，PPP）定位精度可達厘米級，受到了越來越廣泛的關(guān)注［10-21］.當(dāng)前，GNSS 測碼偽距絕對定位坐標(biāo)解算方法有常規(guī)經(jīng)典算法和求差法2種［22-26］：常規(guī)經(jīng)典法同時接收處理4顆GNSS衛(wèi)星信號即可實現(xiàn)定位，但需要GNSS接收機初始概略坐標(biāo)值，再用泰勒級數(shù)展開將偽距觀測方程線性化；求差法不需GNSS接收機初始概略坐標(biāo)值，但需多接收處理1顆GNSS衛(wèi)星信號才能實現(xiàn)定位，在一定程度上限制了GNSS的定位區(qū)域，特別是在地形復(fù)雜的山區(qū)峽谷地帶和森林茂密地區(qū)等GNSS衛(wèi)星信號接收困難區(qū).本文將在對2種方法進行系統(tǒng)分析的基礎(chǔ)上，提出一種改進的測碼偽距絕對定位處理方法，并用實例對其進行驗證.

1 GNSS測碼偽距絕對定位算法分析

1.1 常規(guī)經(jīng)典算法

GNSS測碼偽距定位是依據(jù)GNSS衛(wèi)星坐標(biāo)和測碼偽距為已知量建立觀測方程實現(xiàn)的.其測碼偽距觀測方程可表示為

為實現(xiàn)絕對定位，GNSS接收機需同時接收處理4顆以上GNSS衛(wèi)星信號，建立4個以上觀測方程可解算出4個未知數(shù)：GNSS接收機三維坐標(biāo)（X，Y，Z）及接收機時鐘誤差δtT.但方程式（1）為非線性方程，一般用泰勒級數(shù)展開將其線性化.設(shè)GNSS接收機初始概略坐標(biāo)為（X0，Y0，Z0），將式（1）用泰勒級數(shù)展開，略去2次以上項，則

式中：ρj是經(jīng)過電離層、對流層以及衛(wèi)星時鐘誤差改正后的距離，為已知量；ρj0=kj、mj、nj分別是星站歐氏距離在X、Y、Z向的偏導(dǎo)值.

式（2）為線性化后的測碼偽距觀測方程，可直接求取三維坐標(biāo)增量（坐標(biāo)改正數(shù)）δX、δY、δZ.若GNSS接收機同時觀測處理N顆衛(wèi)星信號（N＞4），用最小二乘法解出GNSS接收機三維坐標(biāo)增量，進而解算出GNSS接收機三維坐標(biāo)

再經(jīng)過若干次遞歸迭代，最后得到GNSS接收機精確坐標(biāo).如上所述，用泰勒級數(shù)展開線性化需預(yù)知GNSS 接收機的初始概略坐標(biāo)值（X0，Y0，Z0）.然而，GNSS接收機初始概略坐標(biāo)很難預(yù)知.鑒于此，通常將其設(shè)為（0，0，0）［25-26］.但若所選取接收機初始概略坐標(biāo)值與精確值偏差較大，則泰勒級數(shù)展開略去的2次以上項將對解算結(jié)果產(chǎn)生顯著影響［25-26］，甚至級數(shù)發(fā)散導(dǎo)致迭代求解失敗.顯然，獲取接收機初始概略坐標(biāo)為此方法的關(guān)鍵所在.

1.2 求差法

為了在GNSS接收機初始概略坐標(biāo)未知情況下獲取接收機精確三維坐標(biāo)值，文獻［22-23］提出了求差線性化觀測方程法.其主要過程為：

先對式（1）移項，再兩邊平方得

若GNSS接收機同時接收5顆GNSS衛(wèi)星信號（j=1，2，3，4，5），即建立 5個測碼偽距觀測方程；再將通式（4）中式j(luò)+1兩邊與式j(luò)（j=1，2，3，4）兩邊相減得到4個線性方程

式中，aj、bj、cj、lj為方程組系數(shù)，aj=2(Xj+1-Xj)；bj=

求解式（5）線性方程組可得到GNSS接收機三維精確坐標(biāo)值（X，Y，Z）.可見，求差法不需預(yù)知GNSS接收機的初始概略坐標(biāo)，通過對偽距觀測方程求差處理實現(xiàn)偽距觀測方程的線性化，以解出接收機三維坐標(biāo).如上所述，此方法需額外增加1顆GNSS衛(wèi)星信號可求解出接收機的精確三維坐標(biāo).這對于在衛(wèi)星信號接收困難區(qū)域，如高樓林立大都市與地形復(fù)雜地區(qū)（山谷、盆地、樹木茂密地帶）的定位與導(dǎo)航增加一定難度.

2 對GNSS測碼偽距絕對定位求解法的改進

GNSS接收機在衛(wèi)星信號接收不理想的區(qū)域只能接收4顆GNSS衛(wèi)星信號，且GNSS接收機初始概略坐標(biāo)未知，無論是常規(guī)經(jīng)典算法和求差法都無法求解GNSS接收機三維精確坐標(biāo).本文將提出一種GNSS測碼偽距絕對定位處理方法，在不增加GNSS衛(wèi)星接收個數(shù)的情況下，即僅接收4顆GNSS衛(wèi)星信號，也可解算出GNSS接收機三維精確坐標(biāo).

首先與求差法中式（5）類似，GNSS接收機同時接收4顆GNSS衛(wèi)星信號則可建立4個測碼偽距觀測方程（j=1，2，3，4），再通過式 j+1兩邊與式 j兩邊相減得到以下線性方程

式中：X、Y、Z為GNSS接收機的三維坐標(biāo)；δtT為接收機時鐘誤差；aj、bj、cj、lj為方程組系數(shù)，將式（6）移項得

由于接收機時鐘穩(wěn)定性通常為10-10～10-11量級，而時鐘誤差δtT一般為±10-7s量級，與自由項系數(shù)lj相比非常微小，現(xiàn)暫將其忽略，則式（7）可簡化為

顯然，由式（8）可直接解出GNSS接收機的三維坐標(biāo).由于此前忽略了微小量δtT，因此如此解出的并非GNSS接收機精確坐標(biāo)值，但可作為其三維初始概略坐標(biāo)（X0，Y0，Z0），現(xiàn)將其值代入式（2），可計算出三維坐標(biāo)增量（坐標(biāo)改正數(shù)）（δX，δY，δZ）；再經(jīng)過遞歸迭代，最后獲取到GNSS接收機精確坐標(biāo)值（X，Y，Z）.

3 實驗分析與驗證

為了驗證本文所提出方法的可靠性與精確性，作者用實測數(shù)據(jù)在GNSS接收機的初始概略坐標(biāo)值未知的情況下，用4顆GPS衛(wèi)星信號數(shù)據(jù)，求取GPS接收機精確三維坐標(biāo).

現(xiàn)以三維信息獲取與應(yīng)用教育部重點實驗室1臺GPS接收機2021年10月20日在GPS歷元TOW（time of week）：268008接收的 GPS導(dǎo)航電文數(shù)據(jù)為例，對本方法進行實驗分析驗證.從衛(wèi)星星歷中獲取的4顆GPS衛(wèi)星坐標(biāo)和各顆衛(wèi)星與接收機間觀測偽距以及電離層和對流層延遲差見表1；用本算法求取的GPS接收機概略坐標(biāo)以及由此迭代計算出的精確坐標(biāo)如表2所示.由此可見，在無GPS接收機初始概略坐標(biāo)情況下，本處理方法僅用4顆GPS衛(wèi)星信號數(shù)據(jù)，計算出了GPS接收機精確三維坐標(biāo).

表1 4顆衛(wèi)星的GPS衛(wèi)星坐標(biāo)和主要參數(shù) 單位：m

表2 用實際數(shù)據(jù)計算出的GPS接收機概略坐標(biāo)與精確坐標(biāo) 單位：m

在計算環(huán)境為Intel（R）Core?i7-2600，CPU為3.4G下，對本方法與常規(guī)經(jīng)典法的工作性能進行了對比，二者都是用4顆衛(wèi)星信號數(shù)據(jù)，在用常規(guī)經(jīng)典算法計算中，按照通常做法將初始概略坐標(biāo)設(shè)為（0，0，0）［24-25］，本文所提出的處理方法計算迭代次數(shù)為4，明顯少于常規(guī)經(jīng)典方法的6次，且消耗時間也較短僅為81 ms，而常規(guī)經(jīng)典方法的消耗時間為124 ms.這說明由于求取了GPS接收機概略坐標(biāo)值而加速了級數(shù)收斂，對于減少迭代次數(shù)發(fā)揮了關(guān)鍵作用，并使得實際計算消耗時間較經(jīng)典方法顯著縮短.此外，需要指出的是，若所選取的初始概略坐標(biāo)值偏差較大，即三維坐標(biāo)改正數(shù) δX，δY，δZ 較 大 .此時常規(guī)經(jīng)典方法用泰勒級數(shù)展開線性化測碼偽距觀測方程時，其2次以上級數(shù)余項不再是“微小項”，很可能因為級數(shù)發(fā)散而導(dǎo)致迭代求解失??；而對于本文所提出的方法，則不存在此問題，因此具有較高的魯棒性.

4 結(jié) 論

現(xiàn)有GNSS測碼偽距絕對定位的2種解算方法各具特色：常規(guī)經(jīng)典算法需預(yù)知GNSS接收機的初始概略坐標(biāo)；求差法不需預(yù)知GNSS接收機初始概略坐標(biāo)但需GNSS接收機多接收處理1顆GNSS衛(wèi)星信號才能實現(xiàn)定位.本文提出了一種改進的GNSS絕對定位處理方法，此方法兼具經(jīng)典常規(guī)法和求差法的優(yōu)勢.即使在未知GNSS接收機初始概略坐標(biāo)的情況下，僅接收4顆GNSS衛(wèi)星信號即可精確求取GNSS接收機的三維坐標(biāo)值.這對于GNSS衛(wèi)星信號接收困難區(qū)域，如在地形復(fù)雜的山區(qū)峽谷地帶、森林茂密地區(qū)以及高樓林立的大都市，實現(xiàn)GNSS絕對定位具有較高的理論意義與應(yīng)用價值.