鄧志鑫，高 宏，王立兵

(1.衛(wèi)星導航系統(tǒng)與裝備技術國家重點實驗室，河北 石家莊 050081； 2.中國人民解放軍63961部隊，北京 100012)

低軌衛(wèi)星導航系統(tǒng)多場景多普勒定位解算方法

鄧志鑫1，高 宏2，王立兵2

(1.衛(wèi)星導航系統(tǒng)與裝備技術國家重點實驗室，河北 石家莊 050081； 2.中國人民解放軍63961部隊，北京 100012)

低軌衛(wèi)星及快速響應衛(wèi)星非常適合作為導航源為地面用戶提供獨立于GNSS的導航定位服務，為了解決用戶在不同低軌衛(wèi)星可見數(shù)量條件下的連續(xù)定位解算問題，并避免傳統(tǒng)偽距定位體制所需要解決的高精度星間時間同步問題，提出多種不同可見星數(shù)量條件下僅利用衛(wèi)星多普勒觀測信息的定位解算方法。仿真結果表明，該方法在大部分區(qū)域內(nèi)都能夠滿足優(yōu)于20 m的三維定位誤差。

低軌衛(wèi)星；多普勒定位；網(wǎng)格搜索；最小二乘

AbstractFast response and low orbit satellites are very suitable for positioning service as navigation sources independent of GNSS for ground users.In order to solve the problem of continuous positioning solution under the condition of different number of low orbit satellites visible to users,and to avoid the increase of system complexity due to high precision synchronization required by traditional pseudo-range positioning system,a positioning method is proposed only using the Doppler satellite observation information under the condition of different visible satellite number.Simulation analysis results show that this method can meet the requirement for three-dimensional positioning error of less than 20 m in most areas.

KeywordsLEO satellite;Doppler locating;grid search;least square

0 引言

在復雜的電磁干擾環(huán)境下,GNSS導航系統(tǒng)中大量的低成本、未加特殊抗干擾措施的導航終端將無法工作[1]。在南北兩極地區(qū)，由于GNSS系統(tǒng)星座采用低傾角設計方案[2]，導致覆蓋性大幅降低。因此，構建GNSS備份無線電導航系統(tǒng)是近年來的研究熱點[3-4]。

低軌衛(wèi)星具有覆蓋區(qū)域廣、可全天候工作以及可觀測性較好等特點[5]，微納衛(wèi)星、快速響應衛(wèi)星還具備建設周期短、成本低廉等優(yōu)勢[6]，因此，低軌衛(wèi)星及快響衛(wèi)星等微小衛(wèi)星平臺非常適合作為導航源為地面用戶提供獨立于GNSS的導航定位服務。

GNSS導航系統(tǒng)用戶通過觀測至少4顆GNSS衛(wèi)星信號的偽距信息進行定位解算[7]。但是該方法應用于低軌衛(wèi)星導航系統(tǒng)存在2個問題：① 由于受到低軌衛(wèi)星數(shù)量限制，大部分區(qū)域內(nèi)的地面用戶無法同時觀測到4顆及以上的低軌衛(wèi)星[8]；② 采用偽距定位解算的前提條件是多顆衛(wèi)星必須要實現(xiàn)嚴格的時間同步[9]，這將極大地增加系統(tǒng)的建設復雜度。而低軌衛(wèi)星由于運動速度很快，具有較好的多普勒觀測性[10]，因此利用多普勒觀測信息進行定位解算是低軌衛(wèi)星導航系統(tǒng)的首選。

近日，美國銥星通信公司宣布其研制的衛(wèi)星授時與定位(STL)系統(tǒng)已經(jīng)具備作為GPS備份或是GPS補充的能力。該系統(tǒng)通過66顆低軌衛(wèi)星及多普勒定位體制為用戶提供比GPS更不易受攻擊的導航服務，但是并未給出具體的定位解算方法。2011年1月，Mark S.Asher等人在國際會議中[11]提出采用單顆LEO衛(wèi)星為地面靜止用戶提供定位手段的思路，但是未給出具體的計算過程，并且其算法需要用戶機發(fā)射信號給低軌衛(wèi)星，這不利于用戶位置的隱蔽。2013年西南交通大學潘磊[12]分析了多普勒定位的原理，并得出了多普勒定位精度與多普勒測量精度、站址位置精度和速度精度有關的結論，但是并未給出適應多種不同可見星數(shù)量條件下的多普勒定位解算方法。

綜上所述，如何利用少量的多普勒觀測信息實現(xiàn)多種不同可見星數(shù)量條件下的連續(xù)定位解算是低軌衛(wèi)星備份導航系統(tǒng)必須解決的關鍵問題。本文提出了一種僅利用衛(wèi)星多普勒觀測信息的精確定位解算方法，給出了計算流程并分析了該方法的仿真試驗性能。

1 多場景條件下的多普勒定位解算方法

1.1 牛頓最小二乘迭代定位解算

由低軌衛(wèi)星備份導航系統(tǒng)的覆蓋性分析可知，地面用戶的可見星數(shù)量由遮蔽角、軌道傾角、軌道高度和每個軌道的衛(wèi)星數(shù)量等因素決定，一般情況下，特定地點的可見衛(wèi)星數(shù)量為1～4顆不等，因此必須采取能夠適應不同可見星數(shù)量的定位解算方法，以滿足不同場景下的連續(xù)導航定位需求。

本文提出能夠應用于多場景下的多普勒定位解算方法，該方法能夠最大程度地滿足不同可見星數(shù)量條件的連續(xù)定位解算需求。每個歷元時刻的定位解算結果由牛頓最小二乘迭代過程獲得，其迭代計算公式為：

[xi+1,yi+1,zi+1]T=[xi,yi,zi]T+[Δxi+1,Δyi+1,Δzi+1]T，

(1)

式中，[xi,yi,zi]T為第i次(i為自然數(shù))迭代后的定位解算結果；[Δxi+1,Δyi+1,Δzi+1]T為第i+1次迭代的定位結果修正量，并有

[Δxi+1,Δyi+1,Δzi+1]T=(GT×G)-1×GT×Y，

(2)

式中，

(3)

Y=

(4)

[x(j),y(j),z(j)]T(j=1，2，…n,n≥3)為第j顆低軌衛(wèi)星或低軌衛(wèi)星第j個觀測點的三維位置坐標；r(j)為第j顆低軌衛(wèi)星或低軌衛(wèi)星第j個觀測點與用戶第i次位置迭代解算結果的觀測距離，即

[Vxj,Vyj,Vzj]T為第j顆低軌衛(wèi)星或低軌衛(wèi)星第j個觀測點當前時刻的三維速度；f0為低軌衛(wèi)星信號的載波頻率；fdj為第j顆低軌衛(wèi)星或低軌衛(wèi)星第j個觀測點的觀測多普勒頻率。當Y的2-范數(shù)小于某固定常數(shù)Dy時，迭代計算結束，得到定位解算結果，否則繼續(xù)進行迭代計算。

1.2 大區(qū)域網(wǎng)格搜索定位解算初值確定

采用牛頓迭代定位解算方法具有計算量小、定位精度高和衛(wèi)星可見星數(shù)量適應性強的特點，但是如果用戶位置初值與真實值距離較遠，如大于100 km以上時，可能會導致迭代計算結果不收斂的情況。為了解決該問題，可采用大區(qū)域網(wǎng)格搜索方法實現(xiàn)用戶在百km量級的粗定位，再將粗定位結果作為牛頓迭代定位解算方法的初值，從而實現(xiàn)大區(qū)域范圍內(nèi)的高精度穩(wěn)健定位解算。

網(wǎng)格搜索定位解算初值確定方法為：首先將大區(qū)域范圍按照經(jīng)緯度劃分為等間距的網(wǎng)格，將每個網(wǎng)格點經(jīng)緯高坐標轉化為ECEF坐標，然后在相同的衛(wèi)星觀測數(shù)量或衛(wèi)星觀測點條件下計算牛頓迭代定位解算方法中的殘差值：

Y=

(5)

式中，[xk,yk,zk]T為第k個網(wǎng)格點的ECEF坐標，求取Y的2-范數(shù)的最小值，其所對應的網(wǎng)格點坐標即可作為牛頓最小二乘多普勒定位解算的迭代初值。

1.3 定位解算工作模式與精度提升

上述多普勒定位解算方法有單歷元定位模式和連續(xù)觀測定位模式2種應用模式。單歷元定位模式是指利用單歷元多普勒觀測信息即可實現(xiàn)用戶的實時定位，一般而言，當用戶可見星數(shù)量≥3時，則可工作于單歷元定位模式。連續(xù)觀測定位模式是指利用多歷元多普勒觀測信息實現(xiàn)靜態(tài)用戶的聯(lián)合定位，一般當用戶可見星數(shù)量<3時，則工作于連續(xù)觀測定位模式。上述為一般情況下的工作模式劃分，有一個特例是如果用戶可見星數(shù)量為2，又具備用戶高程觀測信息，則同樣可工作于單歷元定位模式，即通過2顆低軌衛(wèi)星的多普勒觀測信息和用戶高程觀測信息實現(xiàn)用戶單歷元三維定位。此時的定位解算方法與上面介紹的方法相類似，只不過其中G和Y的計算公式變?yōu)椋?/p>

(6)

Y=

(7)

式中，ah=a+h；bh=b+h；a和b分別為地球的長、短半軸長度；h為用戶的海拔高度估計值。

如果低軌衛(wèi)星可見數(shù)量為1，則只能采用單星多普勒信息連續(xù)觀測的方法實現(xiàn)定位解算，單星連續(xù)觀測的原理如圖1所示。由于低軌衛(wèi)星運行速度快，用戶可觀測時間約為幾分鐘，在有限的觀測時間內(nèi)選取3個及以上的多普勒信息觀測點即可實現(xiàn)用戶的連續(xù)觀測定位解算，觀測點的選取應盡量滿足DOP值最優(yōu)的原則。在此時間內(nèi)用戶本地時鐘的漂移所引入的誤差較小，可忽略不計。

圖1 單星連續(xù)觀測多普勒定位示意

無論是單歷元定位模式還是連續(xù)觀測定位模式，如果用戶處于靜止狀態(tài)，則可通過多組去除野值后的定位解算結果平滑方法進一步提升定位解算精度，此時的定位解算結果[x,y,z]T計算公式為：

(8)

式中，定位解算結果[x,y,z]T為ECEF三維定位坐標；N為定位結果平滑點數(shù)；[xk,yk,zk]T為第k個定位解算結果。

2 定位解算性能仿真驗證

參考美國銥星系統(tǒng)[13]，軌道高度780 km，衛(wèi)星軌道與赤道面傾角為86.4°，星座由6個軌道面組成，每個軌道面11顆星。其全球衛(wèi)星布局如圖2所示。顯然，高緯度地區(qū)的衛(wèi)星可見性較好，赤道地區(qū)的衛(wèi)星可見性最低，因此下面僅仿真分析赤道地區(qū)的星座覆蓋性能。赤道地區(qū)的典型低軌衛(wèi)星布局單元如圖3所示?；静季謫卧目梢娦苑植既鐖D4所示，其中黑色代表可見1顆星、灰色代表可見2顆星、白色代表可見3顆星。

圖2 銥星系統(tǒng)低軌衛(wèi)星全星座布局

圖3 低軌衛(wèi)星赤道地區(qū)典型布局單元

圖4 遮蔽角5°時的基本布局單元可見性分布

設偽衛(wèi)星定軌誤差各方向1 m，多普勒測量誤差1 Hz，衛(wèi)星各方向速度播發(fā)誤差0.1 m/s，用戶高程方差5 m，則可見3顆低軌衛(wèi)星時的多普勒定位誤差分布如圖5所示。圖5中，黑色表示3星定位誤差小于20 m的區(qū)域，約占可定位區(qū)域的95%；灰色表示大于20 m的區(qū)域；白色表示非3星定位區(qū)域。

圖5 三星多普勒定位誤差分布

雙星多普勒及高程聯(lián)合定位誤差分布如圖6所示，其中灰色表示定位誤差小于20 m，約占總定位區(qū)域的85%；白色為大于20 m；黑色為非雙星定位區(qū)域。

圖6 雙星多普勒高程單歷元聯(lián)合定位誤差分布

單星連續(xù)觀測5 min的定位誤差分布如圖7所示。其中灰色表示定位誤差小于20 m，約占可定位區(qū)域的92%；白色為大于20 m；黑色為非單星定位區(qū)域。

圖7 單星連續(xù)觀測5 min多普勒定位誤差分布

采用多普勒觀測量進行定位解算時，每個瞬間都會有若干位置因多普勒觀測值影響導致定位精度較低，但是隨著衛(wèi)星在軌道中運行，定位精度較差的地點快速發(fā)生變化，因此無論哪種多普勒定位解算方法，只要連續(xù)觀測幾個歷元，并進行去除野值后的平滑處理即可得到更為精確的定位精度。經(jīng)過3個歷元平滑后的低軌星座典型布局單元多場景聯(lián)合定位誤差分布如圖8所示，圖8中白色區(qū)域表示定位誤差小于20 m；黑色區(qū)域為定位誤差大于20 m，可知，三維定位誤差優(yōu)于20 m的區(qū)域提升至98.3%，如果進一步增加平滑歷元個數(shù)可進一步提升定位精度。

圖8 多場景多普勒定位誤差小于20 m分布

3 結束語

低軌衛(wèi)星備份導航系統(tǒng)具有覆蓋區(qū)域廣、可全天候工作的優(yōu)勢，采用多普勒定位體制可規(guī)避高精度星間時間同步帶來的系統(tǒng)建設復雜度問題[14]，對于多星間的頻率同步問題，則有多種方法來實現(xiàn)，例如GNSS衛(wèi)星導航校頻、地面雙向校頻和星間校頻等[15]，載荷易于與通信、遙感衛(wèi)星相集成，具有較好的應用可行性。

本文提出了一套可適用于任意可見星數(shù)量的網(wǎng)格搜索與牛頓最小二乘迭代相結合的多普勒定位解算方法，可解決任意低軌星座場景下的自適應定位解算問題，仿真結果表明在典型場景下能夠滿足大于98%的覆蓋區(qū)域內(nèi)實現(xiàn)優(yōu)于20 m的三維定位精度，能夠滿足大部分靜態(tài)及低動態(tài)用戶的導航定位應用需求。

[1] 劉榮，王立平，陳楊.GPS接收機抗干擾性能仿真研究[J].無線電通信技術，2014，40(1)：58-60.

[2] 鄭煒，李斌，任術波，等.GEO衛(wèi)星斜投影多波束的橢圓參數(shù)分析[J].無線電通信技術，2013，39(1)：5-7.

[3] 王玨，何成龍，張金濤，等.基于模糊區(qū)的衛(wèi)星導航空時測量技術研究[J].無線電工程，2015，45(4)：46-49.

[4] 張云，崔曉偉，笪騰飛，等.基于軟件無線電的GNSS干擾和多徑監(jiān)測系統(tǒng)設計[J].電訊技術，2017,57(3)：288-295.

[5] 馮少棟，揭曉，李仰志，等.低軌星座系統(tǒng)業(yè)務建模與仿真[J].宇航學報，2010，31(1)：179-184.

[6] 曾德林.快速響應小衛(wèi)星星座設計及覆蓋性能仿真分析[J].計算機仿真，2014，31(6)：73-77.

[7] 劉競超.復雜環(huán)境下GNSS接收機快速定位研究與分析[D].北京：北京郵電大學，2013.

[8] 楊斌，何鋒，靳瑾，等.LEO衛(wèi)星通信系統(tǒng)覆蓋時間和切換次數(shù)分析[J].電子與信息學報，2014，36(4)：804-809.

[9] 李東新，楊偉才，常青，等.偽衛(wèi)星共視法時間同步模式及精度分析[J].導航定位與授時，2015，2(4)：40-47.

[10] 商春恒.單星無源定位中多普勒信息參數(shù)測量技術研究[D].北京：北京理工大學，2015.

[11] ASHERM S,STAFFORD S J,BAMBERGER R J,et al.Radionavigation Alternatives for US Army Ground Forces in GPS Denied Environments[C]∥International Technical Meeting of the Institute of Navigation,2011：508-532.

[12] 潘磊.基于多普勒頻率差的機載無源定位技術研究[D].成都：西南交通大學，2013.

[13] 潘科炎.銥星星座動力學——系統(tǒng)工程綜合[J].控制工程，2000(1)：29-38.

[14] 徐義.基于多普勒信息的單星無源定位新技術研究[D].長沙：國防科學技術大學，2009.

[15] 吳苗，朱銀兵，李方能，等.無線電導航原理與信號接收技術[M].北京：國防工業(yè)出版社，2015.

Mutiple-sceneDopplerLocatingMethodforLEOSatelliteNavigationSystem

DENG Zhi-xin1,GAO Hong2,WANG Li-bing2

(1.StateKeyLaboratoryofSatelliteNavigationSystemandEquipmentTechnology,ShijiazhuangHebei050081,China; 2.Unit63961,PLA,Beijing100012,China)

TN965.5

A

1003-3106(2017)11-0049-05

鄧志鑫男，(1982—)，博士，高級工程師。主要研究方向：衛(wèi)星導航增強技術。

10.3969/j.issn.1003-3106.2017.11.11

鄧志鑫,高宏,王立兵.低軌衛(wèi)星導航系統(tǒng)多場景多普勒定位解算方法[J].無線電工程,2017,47(11):49-53.[DENG Zhixin,GAO Hong,WANG Libing.Mutiple-scene Doppler Locating Method for LEO Satellite Navigation System[J].Radio Engineering,2017,47(11)：49-53.]

2017-04-18

“十三五”國家重點研發(fā)計劃基金資助項目(2016YFB0502402)。

高宏男，(1979—)，碩士，工程師。主要研究方向：衛(wèi)星導航、大地測量。