陳秉柱，谷守周，魏盛桃，秘金鐘

(1.山東科技大學(xué),青島 266590；2.中國(guó)測(cè)繪科學(xué)研究院,北京 100830)

0 引言

智能化的快速發(fā)展，使得大眾對(duì)位置服務(wù)的精度要求越來(lái)越高，而智能手機(jī)的發(fā)展，為用戶提供了更便捷的位置服務(wù)。由于手機(jī)受到其內(nèi)部天線構(gòu)造、全球?qū)Ш叫l(wèi)星系統(tǒng)(Global Navigation Satellite System，GNSS)接收芯片等因素限制，使得手機(jī)原始觀測(cè)值易受外部環(huán)境影響，手機(jī)定位精度較差。自谷歌為移動(dòng)設(shè)備提供了訪問(wèn)GNSS原始觀測(cè)數(shù)據(jù)的接口，研究人員開(kāi)始評(píng)估移動(dòng)終端采集的GNSS原始觀測(cè)數(shù)據(jù)的質(zhì)量，并分析其定位性能[1-2]，實(shí)現(xiàn)手機(jī)高精度定位成為了研究熱點(diǎn)。

目前，對(duì)于手機(jī)高精度定位的研究主要是實(shí)時(shí)動(dòng)態(tài)定位(Real-Time Kinematic, RTK)技術(shù)和精密單點(diǎn)定位(Precise Point Positioning, PPP)技術(shù)等[3-4]。文獻(xiàn)[5]測(cè)試了單頻RTK和雙頻RTK在小米Mi 8上的性能，并證明了在靜態(tài)模式下，單頻RTK和雙頻RTK都可以提供固定的解決方案。文獻(xiàn)[6]利用小米Mi 8獲取到的全球定位系統(tǒng)(Global Positioning System，GPS)/Galileo雙頻觀測(cè)值進(jìn)行了無(wú)電離層組合解算，靜態(tài)模式下雙頻PPP解算能達(dá)到dm級(jí)精度。文獻(xiàn)[7]對(duì)Nexus 9平板電腦的原始觀測(cè)數(shù)據(jù)的特性進(jìn)行了分析，利用PPP模式能大幅提高手機(jī)的定位精度。目前，手機(jī)GNSS接收芯片大多為單頻芯片，且其載波易發(fā)生周跳，因此需要對(duì)單頻觀測(cè)數(shù)據(jù)的周跳情況進(jìn)行探測(cè)。針對(duì)單頻載波觀測(cè)值的周跳探測(cè)方法主要包括：歷元間高次差法、小波分析法、多普勒積分法等，其中，歷元間高次差法算法簡(jiǎn)單，且隨著采樣間隔的減小，探測(cè)精確度越來(lái)越高[8]，常被用于單頻數(shù)據(jù)的周跳探測(cè)。

本文通過(guò)對(duì)手機(jī)觀測(cè)值進(jìn)行數(shù)據(jù)解析，并使用歷元間高次差法對(duì)手機(jī)載波相位觀測(cè)值進(jìn)行了周跳探測(cè)。在此基礎(chǔ)上，對(duì)偽距進(jìn)行載波相位平滑，同時(shí)獲取國(guó)際GNSS服務(wù)(International GNSS Service，IGS)站通過(guò)Ntrip協(xié)議播發(fā)的RTCM3.0格式的狀態(tài)空間表示(State Space Representation，SSR)方法改正信息，并利用開(kāi)發(fā)的手機(jī)實(shí)時(shí)解算軟件進(jìn)行精密單點(diǎn)定位，以達(dá)到實(shí)時(shí)獲取用戶實(shí)時(shí)精密位置的目的。算法流程圖如圖1所示。

圖1 算法流程圖Fig.1 Algorithm flow chart

1 手機(jī)數(shù)據(jù)解析

1.1 手機(jī)數(shù)據(jù)獲取

2016年5月，谷歌為移動(dòng)設(shè)備提供了訪問(wèn)GNSS原始觀測(cè)數(shù)據(jù)的接口后，歐洲GNSS管理局(European GNSS Agency，GSA)發(fā)布了利用安卓設(shè)備原始觀測(cè)值白皮書，提供了手機(jī)原始觀測(cè)數(shù)據(jù)的觀測(cè)字段及其相關(guān)觀測(cè)值的求解方法，文中利用到的部分字段如表1所示[9]。

表1 Android API字段

Android API并不直接提供偽距觀測(cè)值，偽距ρ可由式(1)計(jì)算得到

(1)

式中，tRx表示接收到信號(hào)的時(shí)刻，tTx表示衛(wèi)星信號(hào)發(fā)出的時(shí)刻，c表示光束。其中

tTx=ReceivedSvTimeNanos

(2)

tRx=TimeNanos+TimeOffsetNanos-

(FullBiasNanos+BiasNanos)

(3)

多普勒和載波相位可分別由式(4)和式(5)計(jì)算得到

(4)

(5)

式中，Doppler為多普勒值，Carrier為載波相位，f1為信號(hào)的頻率。

1.2 占空比

由于連續(xù)運(yùn)行的GNSS芯片組會(huì)增加手機(jī)耗電，因此智能手機(jī)廠商采用占空比(duty cycle)技術(shù)來(lái)保持低功耗(圖2)。占空比技術(shù)導(dǎo)致智能手機(jī)的GNSS載波相位觀測(cè)值跟蹤不連續(xù)[10]，在觀測(cè)和定位過(guò)程中，占空比處于關(guān)閉狀態(tài)，可有效減小觀測(cè)過(guò)程中周跳的產(chǎn)生，提高數(shù)據(jù)的可用性[6]。手機(jī)提供了用戶可以自行開(kāi)啟或關(guān)閉占空比的功能模塊，觀測(cè)過(guò)程中關(guān)閉占空比可減小周跳的發(fā)生概率。

圖2 占空比機(jī)制Fig.2 Duty cycle mechanism

2 原理

2.1 高次差法探測(cè)周跳

由于載波相位平滑偽距易受周跳的影響，首先運(yùn)用高次差法對(duì)單頻相位觀測(cè)值的周跳進(jìn)行了探測(cè)，高次差法運(yùn)用的是歷元間差分組合，對(duì)相鄰2個(gè)歷元間的觀測(cè)值依次作差。假設(shè)某顆衛(wèi)星i的載波相位觀測(cè)值為φi，則以求三次差分為例：

歷元間一次差分為

(6)

歷元間二次差分為

(7)

歷元間三次差分為

Δφ″1=Δφ′2-Δφ′1

=φi(k+2)-3φi(k+1)+3φi(k)-φi(k-1)

(8)

在不發(fā)生周跳的情況下，載波相位的變化隨衛(wèi)星與接收機(jī)的距離變化而變化，且是平緩而有規(guī)律的。周跳將破壞這一規(guī)律的特點(diǎn)，利用高次差法對(duì)相位多次差分后周跳會(huì)被放大，從而探測(cè)觀測(cè)值中是否含有周跳。

假設(shè)鐘的短期穩(wěn)定度為5×10-10，考慮到L1載波的頻率f1=1.57542×109Hz，則時(shí)鐘的隨機(jī)誤差給相鄰載波相位造成的周數(shù)影響如下

cycle=s×5×10-10×1.57542×109

(9)

式中，cycle表示周跳數(shù)，s表示采樣間隔。

2.2 廣域精密定位算法原理

廣域增強(qiáng)技術(shù)通過(guò)對(duì)GNSS觀測(cè)量的誤差源加以區(qū)分，然后計(jì)算出每一個(gè)誤差源的誤差修正值(差分改正值)，利用數(shù)據(jù)通信鏈傳輸給用戶，并對(duì)用戶GNSS接收設(shè)備的觀測(cè)值誤差加以改正的方法，以達(dá)到削弱這些誤差源的影響，改善用戶GNSS定位精度的目的。

2.2.1 軌道改正數(shù)

當(dāng)解析出SSR信息中的相應(yīng)數(shù)據(jù)后，可得到k時(shí)刻衛(wèi)星位置改正量為

(10)

(11)

則根據(jù)廣播星歷得到的衛(wèi)星坐標(biāo)和改正量獲得觀測(cè)時(shí)刻k的衛(wèi)星位置為

隨著互聯(lián)網(wǎng)技術(shù)的不斷發(fā)展，中國(guó)各方面水平得到顯著提升，但是在教育行業(yè)有許多班主任的觀念卻沒(méi)有與時(shí)俱進(jìn)，很多班主任受傳統(tǒng)教學(xué)觀念的影響，總是將傳統(tǒng)應(yīng)試教育的教育方式放到新時(shí)代班級(jí)管理中。對(duì)班級(jí)管理的重要性沒(méi)有全面認(rèn)識(shí)，將班主任與學(xué)生的地位理解得不夠現(xiàn)代化，一味地認(rèn)為在校學(xué)生的唯一任務(wù)就是學(xué)習(xí)，還有許多班主任只是將班級(jí)作為教學(xué)組織——教師講課，學(xué)生聽(tīng)，并沒(méi)有按照學(xué)生自己的喜愛(ài)和特長(zhǎng)對(duì)學(xué)生進(jìn)行針對(duì)性教育，更多關(guān)注的是學(xué)生的成績(jī)，不重視培養(yǎng)學(xué)生的綜合能力。這樣的管理現(xiàn)狀對(duì)學(xué)生而言不利于其健康成長(zhǎng)。

(12)

2.2.2 鐘差改正數(shù)

SSR信息中提供的衛(wèi)星鐘差改正數(shù)為：參考時(shí)間t0及計(jì)算衛(wèi)星鐘差多項(xiàng)式的3個(gè)參數(shù)[13-14]，即C0、C1、C2。則IGS實(shí)時(shí)數(shù)據(jù)給出的衛(wèi)星鐘差修正值ΔC為

ΔC=C0+C1(tk-t0)+C2(tk-t0)2

(13)

利用廣播星歷計(jì)算出的衛(wèi)星鐘差和實(shí)時(shí)數(shù)據(jù)的修正值進(jìn)行修正得到的鐘差ΔT為

(14)

式中，T為廣播星歷計(jì)算出的衛(wèi)星鐘差，Vc為光速。

3 實(shí)驗(yàn)分析

實(shí)驗(yàn)以華為P30智能手機(jī)作為智能終端，定位使用GPS單系統(tǒng)，在北京房山人衛(wèi)激光國(guó)家野外科學(xué)觀測(cè)研究站的已知點(diǎn)進(jìn)行數(shù)據(jù)采集和靜態(tài)測(cè)試。在觀測(cè)環(huán)境良好，占空比處于關(guān)閉狀態(tài)的條件下，利用Geo++ RINEX Logger程序采集采樣間隔為1s的P30手機(jī)原始觀測(cè)數(shù)據(jù)文件，以及利用R9接收機(jī)進(jìn)行同步觀測(cè)，觀測(cè)時(shí)段為14:03:54～15:14:50,獲取其觀測(cè)文件，并用于手機(jī)載波相位的周跳分析；將P30手機(jī)安裝自己開(kāi)發(fā)的實(shí)時(shí)偽距單點(diǎn)定位和載波平滑偽距單點(diǎn)定位程序，放置在已知點(diǎn)上進(jìn)行測(cè)量解算，并對(duì)手機(jī)載波相位平滑偽距的可用性進(jìn)行分析；然后，利用開(kāi)發(fā)的廣域精密單點(diǎn)定位程序，采樣間隔設(shè)置為1 s，進(jìn)行實(shí)時(shí)單點(diǎn)解算并記錄生成結(jié)果文件。

3.1 周跳分析

實(shí)驗(yàn)基于P30手機(jī)和R9接收機(jī)采集的原始觀測(cè)數(shù)據(jù)，運(yùn)用高次差法分別對(duì)其載波相位進(jìn)行歷元間四次差處理，利用歷元間差分?jǐn)?shù)據(jù)生成各衛(wèi)星差分值隨歷元的變化圖，并統(tǒng)計(jì)了手機(jī)和接收機(jī)觀測(cè)時(shí)段衛(wèi)星發(fā)生周跳的歷元數(shù)，結(jié)果見(jiàn)表2和表3。

表2 P30手機(jī)周跳統(tǒng)計(jì)

表3 R9接收機(jī)周跳統(tǒng)計(jì)

圖3和圖4所示分別為華為P30手機(jī)和Trim-ble NetR9 GNSS參考站接收機(jī)G3和G25號(hào)衛(wèi)星相位歷元間四次差的周跳探測(cè)圖。由圖3和圖4可以看出：手機(jī)在觀測(cè)條件良好，關(guān)閉占空比的情況下，周跳多是范圍在1～4周的小周跳。表2和表3對(duì)手機(jī)和接收機(jī)觀測(cè)衛(wèi)星的周跳統(tǒng)計(jì)結(jié)果表明：手機(jī)載波相位值在觀測(cè)過(guò)程中，周跳現(xiàn)象比大地型接收機(jī)更為嚴(yán)重，表2中部分衛(wèi)星周跳統(tǒng)計(jì)歷元數(shù)較多，分析原因?yàn)橛^測(cè)過(guò)程中未獲取到手機(jī)的載波相位觀測(cè)值。因?yàn)楦叽尾罘o(wú)法探測(cè)連續(xù)歷元間的周跳情況，所以實(shí)際觀測(cè)中周跳出現(xiàn)的概率要高于統(tǒng)計(jì)值。

圖3 P30手機(jī)G3、G25衛(wèi)星周跳序列Fig.3 G3,G25 satellite cycle slip sequence of P30 phone

圖4 R9接收機(jī)G3、G25衛(wèi)星周跳序列Fig.4 G3,G25 satellite cycle slip sequence of R9 receiver

3.2 載波相位平滑偽距可用性分析

在周跳分析的基礎(chǔ)上，將載波相位平滑窗口設(shè)置為100，進(jìn)行載波相位平滑偽距實(shí)時(shí)單點(diǎn)定位，分析小周跳對(duì)定位精度的影響，驗(yàn)證載波相位平滑偽距的可用性(圖5、圖6)。實(shí)驗(yàn)表明，在手機(jī)關(guān)閉占空比的條件下，雖然周跳相對(duì)于大地型接收機(jī)較多，但采用載波平滑偽距進(jìn)行單點(diǎn)定位依然可以提高定位精度，水平方向精度提高了36.2%，垂直方向精度提高了26.3%，平滑后偽距具有較好的可用性，統(tǒng)計(jì)結(jié)果如表4所示。

圖5 偽距單點(diǎn)定位Fig.5 Pseudorange single point positioning

圖6 載波平滑偽距單點(diǎn)定位Fig.6 Carrier smoothed pseudorange single point positioning

表4 偽距平滑前和平滑后定位精度統(tǒng)計(jì)

3.3 廣域精密定位精度分析

廣域精密單點(diǎn)定位利用載波平滑后偽距和載波相位觀測(cè)值，接入IGS站播發(fā)的廣播星歷和掛載點(diǎn)為IGS03中GPS的實(shí)時(shí)SSR改正信息進(jìn)行解算，測(cè)量時(shí)長(zhǎng)為4.5h，對(duì)實(shí)時(shí)單點(diǎn)定位存儲(chǔ)的結(jié)果文件進(jìn)行處理，并分析其定位精度。圖7所示為定位結(jié)果序列，實(shí)時(shí)定位結(jié)果在E、N、U方向的殘差統(tǒng)計(jì)結(jié)果如表5所示。結(jié)果顯示，精度收斂至1m用時(shí)約20min，收斂后精度可達(dá)到dm級(jí)，相較于原始偽距單點(diǎn)定位精度，在水平方向上提高了66.7%，垂直方向提高了28.9%，水平方向上精度提高更加明顯。

圖7 廣域精密單點(diǎn)定位Fig.7 Wide-area precision single point positioning

表5 偽距單點(diǎn)和廣域精密單點(diǎn)定位精度統(tǒng)計(jì)

4 結(jié)論

本文針對(duì)智能手機(jī)原始觀測(cè)值不能滿足實(shí)時(shí)高精度定位需求的問(wèn)題，采用廣域增強(qiáng)技術(shù)的方法實(shí)現(xiàn)了實(shí)時(shí)高精度定位，分析和實(shí)驗(yàn)結(jié)果表明：

1)華為P30手機(jī)獲取到的單頻相位觀測(cè)值相比大地型接收機(jī)周跳更為嚴(yán)重；

2)在手機(jī)相位觀測(cè)值發(fā)生小周跳的情況下，載波平滑偽距依然可以提高偽距的精確度，進(jìn)而提高定位精度；

3)在靜態(tài)測(cè)試條件下，廣域精密定位可以有效提高手機(jī)的定位精度，水平方向精度提高較為明顯，平面精度可達(dá)1.5m，但高程方向和收斂速度結(jié)果有待提高，需要進(jìn)一步研究；

4)在廣域精密定位可有效提高手機(jī)觀測(cè)值靜態(tài)定位精度的基礎(chǔ)上，進(jìn)行動(dòng)態(tài)測(cè)試將是下一步的工作重點(diǎn)。