1、第 11 卷 第 2 期 导航定位学报 Vol.11，No.2 2023 年 4 月 Journal of Navigation and Positioning Apr.，2023 引文格式：王森,王腾飞,姚铮,等.GNSS/5G 观测量融合定位算法J.导航定位学报,2023,11(2):71-79.（WANG Sen,WANG Tengfei,YAO Zheng,et al.Integrated positioning algorithm based on GNSS/5G observationsJ.Journal of Navigation and Positioning,2023,11(
2、2):71-79.）DOI:10.16547/ki.10-1096.20230208.GNSS/5G 观测量融合定位算法 王 森1，王腾飞2,3，姚 铮2,3，邢建平1（1.山东大学 微电子学院，济南 250100；2.清华大学 北京信息科学与技术国家研究中心，北京 100084；3.清华大学 电子工程系，北京 100084）摘要：针对随着泛在位置服务的需求扩增，卫星导航系统在城市峡谷等环境下的定位效果不连续、定位精度低，无法满足服务需求的问题，提出一种全球卫星导航系统（GNSS）/第五代移动通信技术（5G）观测量融合定位算法：给出基于最小二乘法的 GNSS/5G 观测量融合定位算法模型；并对
3、比分析 GNSS/5G 融合定位与 GNSS 单系统定位效果，通过增加5G 到达时间差（TDOA）观测量，配合 GNSS 观测量参与定位解算，给出定位性能分析结果。实验结果表明，GNSS 信号观测数据与 5G 信号 TDOA 观测量融合定位相比 GNSS 单系统定位在水平方向、三维（3D）方向定位误差优化精度分别约为 64.6%、58.21%，可有效提升定位性能，为可扩展定位需求的可行性相关研究提供参考。关键词：全球卫星导航系统（GNSS）；第五代移动通信技术（5G）；到达时间差（TDOA）；融合定位；最小二乘法 中图分类号：P228 文献标志码：A 文章编号：2095-4999(2023)0
4、2-0071-09 Integrated positioning algorithm based on GNSS/5 5G observations WANG Sen1,WANG Tengfei2,3,YAO Zheng2,3,XING Jianping1（1.School of Microelectronics,Shandong University,Jinan 250100,China;2.Beijing National Research Center for Information Science and Technology,Tsinghua University,Beijing 1
5、00084,China;3.Department of Electronic Engineering,Tsinghua University,Beijing 100084,China）Abstract：Aiming at the problem that with the increase of the demands for ubiquitous location services,the location performance of satellite navigation systems in urban canyons and other scenes is discontinuou
6、s with low position accuracy,which cannot meet the service requirements,the paper proposed an integrated positioning algorithm based on global navigation satellite system(GNSS)/fifth-generation(5G)observations:an integrated positioning algorithm model for GNSS/5G observations based on the least squa
7、res method was given;and the positioning performances of GNSS/5G integration and a single GNSS were comparatively analyzed;by attaching 5G time difference of arrival(TDOA)observations into the GNSS observations to solve the positioning,the analysis results of positioning performance were given then.
8、Experimental result showed that compared with a single GNSS positioning,the positioning error optimization accuracy of the integrated positioning of GNSS observations and 5 G TDOA observations in the horizontal direction and three dimensional directions would be about 64.6%and 58.21%,respectively,wh
9、ich could effectively improve the positioning performance,providing a reference for the related study on the feasibility of scalable positioning requirements.Keywords:global navigation satellite system(GNSS);fifth-generation(5G)mobile communication technology;time difference of arrival(TDOA);integra
10、ted positioning;least squares 收稿日期：2022-11-18 资助项目：国家重点研发计划项目（2021YFA0716603）；国家自然科学基金面上项目（42274018）；第六届中国科协青年人才托举工程项目（YESS20200226）。作者简介：王森（1994），男，山东淄博人，博士研究生，研究方向为 GNSS 信号基带处理算法。72 导航定位学报 2023 年 4 月 0 引言 卫星导航作为时空信息基础设施的重要组成部分，为用户提供全天时、全天候的时空位置服务，对国防军事、生产生活、行业应用等方面发挥着重要的作用和价值，因此，各国都在加紧建设自身的卫星导航系统，
11、不断革新迭代导航技术，提供精准可靠、连续稳定的定位、导航、授时（positioning,navigation and timing，PNT）服务1。受卫星导航系统自身局限性，全球导航卫星系统（global navigation satellite system，GNSS）信号传输过程中易受干扰、遮蔽、多径等影响，在城市峡谷等环境中信号接收强度降低，从而导致出现定位效果不连续、定位精度低，无法满足服务需求等现象2。在此背景下，研究提高复杂环境下的定位能力，满足不同区域的无缝定位需求，提供泛在位置服务，具有重要的作用及意义3。移动通信系统具备的覆盖范围广、建设成本低、可靠性高等优势，使其可提供多场
12、景下的定位服务，有效弥补卫星导航定位服务的不足4。第五代移动通信技术（fifth-generation mobile communication technology，简称 5G）具备的密集组网、大带宽、多天线等有利的定位条件，使其 5G 信号应用于定位服务具有信源密度高、信号功率强、作用距离短、伪距测量精度高、信号带宽资源丰富、信号多径免疫性强等优势，通常具备更高的定位精度5。5G 毫米波信号提供的基于到达时间（time of arrival,TOA）和基于到达时间差（time difference of arrival,TDOA）的定位测量技术，可提供精准的伪距观测量6。在此背景下，将 5
13、G 系统与卫星导航系统融合应用，可实现相互赋能，拓展精准定位服务范围7-8。本文中以 GNSS 信号观测数据和 5G 信号TDOA 观测量作为研究对象，研究基于最小二乘法的 GNSS/5G 观测量融合定位算法模型，并对比分析 GNSS/5G 融合定位与 GNSS 单系统定位效果，以验证 GNSS 与 5G 系统融合定位的可行性。1 定位原理 1.1 基于最小二乘法的 GNSS 定位原理 GNSS 信号伪距观测方程为 ()()()()()uuusususuu=+ddddnnnnnnnnnnnnnRcttITXXYYZZcttIT-+=-+-+-+-+222（1）式中：un为第n颗卫星的伪距观测量
14、；unR为第n颗卫星与接收位置的几何距离；c为光速；udt为接收机钟差；dnt为第n颗卫星的卫星钟差；nI为第n颗卫星电离层误差；nT为第n颗卫星对流层误差；n为第n颗卫星其他观测噪声；()sss,nnnXYZ为第n颗卫星空间位置坐标；()uuu,XY Z为接收机空间位置坐标。在上述伪距观测方程中，将接收机空间位置坐标及接收机钟差参数设为待估计未知参数，满足4颗及以上观测卫星，构建观测方程组，可求得待估计未知参数。观测方程组为()()()()()()()()usususuuusususuu dd ddNNNNNNNNXXYYZZcttITXXYYZZcttIT=-+-+-+|-+|=-+-+-
15、+|-+22211111111222（2）式中 N 为观测卫星总数。将上述观测方程组用向量形式表述为 ()=+Yh XV（3）式中：Tuu,N=Y1为伪距观测量；h 为向量函数；Tuuuu,dXY Zt=X为待估计未知参数；V 为观测噪声。利用最小二乘法中的高斯-牛顿法9，将上述线性化方程在X0处泰勒级数展开，求得线性化误差方程为 d=+VGXb（4）其中：sususuuuusususuuuuNNNNNNXXYYZZRRRXXYYZZRRR-|=|-|G11111111（5）Tuuuud,dXYZt=X（6）第 2 期 王 森，等.GNSS/5G 观测量融合定位算法 73 uuuuuudddd
16、NNNNNRctctITRctctIT-+-|=|-+-b11111（7）式中：G为函数 h的雅克比矩阵，由观测方程对待估计未知参数求偏导的系数矩阵组成；dX为待估计未知参数的校正量；b为观测量与前一历元观测值的差值，经过多次迭代运算，当差值收敛小于阈值，可求得最优解。解算最优解方程为 ()TTd-=XG GG b1（8）d=+XXX0（9）式中X0为前一历元求得的接收机空间位置坐标。在缺少接收机先验信息时，可设T,xy z=X00000，其中，可预设()T,xy z=0000 0 0或接收机先验近似位置坐标。1.2 GNSS/5G 融合定位原理 5G信号TDOA观测方程为10()()()()()(),GGGu,Gu,Gu,Gu,G,u,Gu,GGuGuGuGuGuGu,u,G ddd d kkkkkkkkkkkRctRctctXXYYZZXXYYZcRZt=-=+-=+=-+-+-+-+-+115551115555115522255522211155515 （10）式中：,Gk15为第 k 个 5G基站与 5G参考基站之间的伪距观测量差值；G15为 5G参考基站与接收位置的伪距观测量
