Android智能手机GNSS数据质量分析

张垠，涂锐,3，洪菊，李芳馨，侯福荣

( 1.中国科学院国家授时中心，西安 710600；2.中国科学院大学，北京 100049；3.中国科学院精密导航定位与定时技术重点实验室，西安 710600 )

0 引言

当今时代，手机等智能设备由于具有成本低、便携性好的特点，该设备已经成为人类生活和工作中不可或缺的工具.目前，智能手机基本都搭载了全球卫星导航系统(GNSS)模块，极大地方便了人们的日常生活.各大卫星导航系统逐步发展，可见卫星数进一步增加，也为基于移动位置服务(LBS)的快速发展提供了保证[1].未来公众对LBS服务的需求将会越来越大，对定位精度、稳定性要求也会越来越高.因此，手机等智能设备的高精度定位受到学者的广泛关注，成为当前导航定位领域的一个研究热点.

在早期的Android 系统版本中，GNSS模块只开放了获取GNSS芯片输出的位置、速度等信息的应用程序编程接口(API)，开发者和用户很难得到原始的GNSS观测数据，因此难以开展进一步的研究.2016年5月谷歌公司宣布，从Android N 版本操作系统开始，开发者和用户可以通过API 来得到GNSS的伪距、载波相位、多普勒等原始观测数据，从而具备进行精确计算的条件，这对于研究智能设备的精确定位有着重要的意义.与测量型等专用GNSS接收机不同，手机、平板电脑等这些大众智能设备，往往采用体积小、成本低的线性极化天线[2]和低成本、低功耗的导航芯片，这必将影响智能设备的GNSS原始观测数据质量，而原始观测数据质量好坏是影响定位精度的主要因素.

国内外许多学者对智能设备的GNSS数据质量和定位精度进行了研究和分析.早在谷歌还未开放获取GNSS原始数据的API 时，就已经有学者对此产生了强烈的兴趣.Kirkko-Jaakkola 等[3]和Humphreys T E等[4]使用定制的智能手机诺基亚Lumia 1520和三星Galaxy S5获得了GNSS原始观测数据，他们发现智能手机的GNSS数据噪声和异常值均明显高于低成本接收机.Håkansson M[5]通过Nexus 9平板电脑评估了多路径效应对精确定位的影响.Paziewski J 等[6]使用华为P20和Javad Alpha+GrAnt_G3T 天线评估了智能手机的GNSS数据质量.Riley 等[7]指出手机由于电池容量限制存在占空比(duty cycle)机制，这种机制降低了功耗但会导致载波相位观测数据不连续.LIU 等[8]、陈波等[9]指出手机GNSS数据存在载波相位与伪距不一致的问题，并且通过平滑伪距的方式进行了单点定位及差分定位的测试.薄亚东等[10]对华为Mate30和华为honor V20的GNSS原始观测数据质量和定位结果进行分析，指出智能手机会发生频繁的失锁现象且多路径误差明显，差分定位存在无法固定模糊度的问题.总体来说，目前针对智能手机GNSS定位的研究相对较少，且大部分并未将单双频手机进行对比分析.

本文同时使用单频、双频智能手机采集GNSS原始观测数据，并将其与北斗星通UR4B0-D高性能GNSS接收机的数据作对比，从卫星可见性、载噪比(CNR)、高度角、多路径效应等方面进行分析，评估智能手机GNSS的原始观测数据质量，为后续手机高精度定位应用提供参考.

1 智能手机GNSS原始观测值获取与分析方法

1.1 智能手机GNSS原始观测值获取

自Android 7.0(即N 版本)开始，Android API等级升级至24，开发者可以通过GNSSClock 和GNSSMearement 两类中的字段直接或者间接的获取伪距观测值、载波相位观测值、多普勒观测值和CNR 等.相关字段如表1所示[11].

表1 Android API 时钟和测量类字段信息

GPS时间、伪距观测值、多普勒观测值并不能直接通过Android API获取到，需要通过计算间接得到.其中GPS时(GPST)可由下式计算得到：

式中，c为光速.

其中TimeNanos和TimeOffsetNanos取当前历元值，而FullBiasNanos和BiasNanos取第一个历元从API 获取到的值.这是因为由于手机钟差不稳定，会导致接收机硬件时钟与GPST的差值会随时间变化，以第一个历元的偏差值作为偏差估计值可以保证接收数据的时间间隔一致，避免给后续处理带来麻烦.

多普勒的计算公式如下：

载波相位观测值及CNR 可以分别通过AccumulatedDelataRange和Cn0DbHz 字段得到.

由上可知，Andriod N 版本开始应用程序可以通过相关API获取分析所需的GNSS原始数据.许多公司利用这些应用程序接口，开发了可以输出GNSS原始观测数据的应用程序(APP)，如谷歌公司的发布的GNSSLogger，Geo++公司发布的Geo++RINEX Logger 等.但由于GNSSLogger 无法直接以RINEX格式输出原始数据，不便于后续分析研究，而Geo++RINEX Logger 可以以RINEX 标准格式输出原始数据，因此选择使用Geo++RINEX Logger 采集原始数据.

1.2 多路径误差分析方法

在卫星信号传输过程中，信号会在接收机附近的地面、墙壁等物体表面发生反射现象，反射后的信号与卫星直接发射的信号同时被接收机接收，从而产生多路径误差.降低多路径误差主要通过使用较好的GNSS天线和选择较为开阔的测量环境两个方面进行.但Android 智能手机受到体积限制，只能使用体积很小的线性极化贴片天线，这使得Android 智能手机受到多路径效应的影响较为严重.同时，Android智能手机的使用环境常常位于高楼林立的城市峡谷中，这使得多路径误差进一步增大，成为了手机定位的主要误差项之一.

利用伪距和载波相位观测值的线性组合，再消去双频信号计算得到的电离层延迟，可以得到L1和L5频点的多路径效应，计算公式如式(6)所示：

2 实验简介

本文使用三星S9+(Exynos)(下文简称三星S9+)、华为Mate30、华为P40 Pro进行静态观测实验，并且与测量型接收机北斗星通UR4B0-D的同步观测数据做对比，3款智能手机的具体参数如表2所示.实验地点为中国科学院国家授时中心临潼园区的综合楼楼顶，视野开阔无遮挡物，观测条件良好，采样率为1 s，采样时长为3 h.实验环境如图1所示.

表2 三款智能手机具体参数

图1 实验环境

部分Android 智能手机为了减少电量消耗以达到更长的续航时间，往往采用占空比技术.在占空比技术下，一个脉冲循环内只有一部分时间可以通电，因此会导致载波相位观测值不连续[7].而部分手机可以在开发者选项中手动选择“强制进行全GNSS测量”以关闭占空比技术，如图2所示.本文选用的三星S9+采用该方法关闭占空比技术，而华为P40 Pro和华为Mate30默认关闭占空比技术.

图2 开发者选项开启“强制进行全GNSS测量”

3 实验结果分析

3.1 可见卫星数分析

对手机原始观测数据进行分析后发现，3台Android 手机均可以接收北斗二号(BDS-2)和北斗三号(BDS-3)卫星数据，但2台双频手机接收到的有双频数据的卫星数量都很有限.华为Mate30可以接收到7颗同时具有L1+L5信号的GPS卫星，3颗同时具有L1+L5的准天顶卫星系统(QZSS)卫星，在不同时段最多可以同时接收到7颗具有E1+E5a 信号的Galileo卫星；华为P40 Pro可以接收到4颗同时具有L1+L5信号的GPS卫星，3颗同时具有L1+L5的QZSS卫星，8颗同时具有E1+E5a 信号的Galileo卫星.但部分卫星的双频信号并不稳定，存在大量的数据缺口，导致实际可用的卫星数目远远小于观测到的可见卫星数.

在3 h 的观测时间内，剔除掉未接收到数据的卫星后，3台智能手机和北斗星通UR4B0-D接收机各历元观测卫星数如图3～6所示.可以非常明显的发现，接收机的可见卫星数非常稳定且观测到的卫星数目是最多的，约为45～50颗.3台智能手机中，三星S9+的可见卫星数是3台手机中最少的，约为25～34颗，但相较于另外2部手机卫星数的波动较小，稳定性更好.华为P40 Pro的可见卫星数约为32～38颗，但卫星数的波动是3台手机中最大的.华为Mate30接收到的卫星数目约为35～44颗，可见卫星数的波动较大，但在第6 224历元附近出现了数据采集中断，导致图形出现不正常的跳变.华为P40 Pro和华为Mate30的可见卫星数都有较大的波动，这种波动主要是由于GLONASS和北斗卫星导航系统(BDS)的可见卫星数波动较大导致的，而GPS和Galileo系统的可见卫星数波动较小.

图3 接收机可见卫星数

图4 三星S9+可见卫星数

图5 华为P40 Pro可见卫星数

图6 华为Mate30 可见卫星数

总体来说，相对于接收机而言手机观测能力和稳定性均表现较差，且不同机型之间由于采用不同的GNSS芯片和天线，观测能力也存在着较大的差异.即使是同品牌的手机，这种差异也比较明显.

3.2 CNR 分析

接收机输出的CNR 反映了所跟踪GNSS卫星的观测值精度和接收机天线的噪声密度，是反映测量噪声水平的重要指标[11].CNR 大小和数据质量成正比.从每个卫星系统中选取3颗卫星，计算接收机和3台手机所选卫星所有历元的平均CNR.不同设备在L1频点的CNR 对比如图7所示.

图7 L1频点不同设备不同卫星CNR 对比

由于三星S9+无法接收双频信号，因此在L5频点只进行华为Mate30和华为P40 Pro与接收机的对比.选取GPS的5颗卫星进行对比，如图8所示.

图8 L5频点不同设备不同卫星CNR 对比

由图7、图8可知，Andriod 智能手机的CNR 明显低于测量型接收机，L1频点3台手机的CNR 比测量型接收机低约10～15 dB·Hz.3台手机中只能接受单频数据的三星S9+CNR 最大，大部分卫星平均CNR 达到了35～40 dB·Hz.而华为P40 Pro的CNR最小，平均CNR 约在30 dB·Hz.在L5频点上，接收机的CNR 均在40 dB·Hz 以上，而两台双频手机的CNR 均较低，华为Mate30 CNR 约为29～34 dB·Hz，而华为P40 Pro CNR 仅有22～27 dB·Hz.较低的CNR 说明手机接收的数据中存在着大量的噪声.

3.3 高度角分析

选取G31、G27两颗卫星分析其L1频点CNR和高度角的关系，图9～11分别为接收机、华为Mate30和华为P40 Pro对G31卫星CNR 和高度角分析的结果，图12～14是三者对G27卫星CNR 和高度角分析的结果.华为P40 Pro仅在6 064历元中接收到G31卫星的数据.红色线为高度角，可见对于G31卫星，在所有历元中高度角由大到小变化，接收机的CNR 也随之减小约15 dB·Hz，但两部手机的CNR 没有表现出逐渐减小的趋势，甚至高度角降低时反而有CNR 增大的现象.对G27卫星，高度角由小到大变化.在高度角较小时，并未出现CNR 随高度角减小而一直减小的现象.可见手机的CNR 与卫星高度角的相关性较小.推测是因为手机接收到的数据中已经包含有大量噪声和多路径误差，高度角减小时噪声和多路径误差存在波动的情况，但并不会增大太多.

图9 接收机G31卫星CNR 与高度角

图10 华为Mate30 G31卫星CNR 与高度角

图12 接收机G27 卫星CNR 与高度角

图11 华为P40 Pro G31卫星CNR 与高度角

图13 华为Mate30 G27 卫星CNR 与高度角

图14 华为P40 Pro G27卫星CNR 与高度角

3.4 多路径误差分析

本文分别使用华为Mate30和华为P40 Pro双频GPS信号计算多路径误差，选取G04卫星数据分别对L1和L5频点进行计算.将两台智能手机多路径误差与接收机UR4B0-D进行比较，计算得到G04卫星L1、L5频点的伪距多路径结果如图15～18所示.

图15 华为Mate30和接收机G04卫星L1频段多路径误差

图16 华为Mate30和接收机G04卫星L5频段多路径误差

图17 华为P40 Pro和接收机G04卫星L1频段多路径误差

图18 华为P40 Pro和接收机G04卫星L5频段多路径误差

计算各终端G04卫星多路径误差均方根(RMS)，结果如表3所示.华为Mate30 从5 000历元开始数据噪声突然增大，导致多路径误差计算值增大，因此只分析前5 000历元的多路径误差.由图15～18以及表3可知，智能手机的多路径误差远远大于接收机的多路径误差，而华为Mate30的多路径误差无论是L1频点还是L5频点均小于华为P40 Pro.对比点L1和L5频点的多路径误差发现，接收机和华为Mate30在L5频点的多路径误差明显小于L1频点的多路径误差，华为P40 Pro在L5频点的多路径误差比L1频点的稍大.可见L5频点一定程度上可以起到抑制多路径误差的作用.总体来说，智能手机由于天线特点导致多路径误差很大，是影响智能手机高精度定位的主要误差.

表3 不同终端G04卫星多路径误差m

4 结束语

本文使用三星S9+、华为Mate30和华为P40 Pro进行GNSS原始数据采集，并与北斗星通UR4B0-D接收机同步采集的数据进行对比，从可见卫星数、CNR、高度角、多路径误差等方面分析了智能手机的GNSS数据质量.结果表明：与接收机相比，Andriod智能手机GNSS原始观测数据质量较差，且不同机型之间存在较大差异.Andriod 智能手机的GNSS数据存在着数据缺失现象，导致可进行数据处理的卫星数目小于可见卫星数目.对CNR 和高度角进行分析后发现，智能手机的CNR 较小，且没有表现出随高度角降低而减小的趋势.对多路径误差进行分析，在L1频点上接收机的多路径误差RMS为0.014 4 m，华为Mate30为2.124 1 m，华为P40 Pro为3.227 9 m；而在L5频点上接收机的多路径误差RMS为0.010 8 m，华为Mate30为1.816 7 m，华为P40 Pro为3.424 2 m，可见多路径误差是影响智能手机高精度定位的主要误差，在将来的定位应用研究中应着重进行处理.

致谢：本文的研究得到了国家自然科学基金项目(41674034,41974032)和中组部、中科院高层次青年人才项目以及王宽诚教育基金会的支持.