CN112505735A - 对终端进行定位的方法、装置和存储介质 - Google Patents

对终端进行定位的方法、装置和存储介质 Download PDF

Info

Publication number
CN112505735A
CN112505735A CN202011371254.6A CN202011371254A CN112505735A CN 112505735 A CN112505735 A CN 112505735A CN 202011371254 A CN202011371254 A CN 202011371254A CN 112505735 A CN112505735 A CN 112505735A
Authority
CN
China
Prior art keywords
terminal
satellite
satellites
pseudo
observations
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011371254.6A
Other languages
English (en)
Other versions
CN112505735B (zh
Inventor
苏景岚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tencent Technology Shenzhen Co Ltd
Original Assignee
Tencent Technology Shenzhen Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tencent Technology Shenzhen Co Ltd filed Critical Tencent Technology Shenzhen Co Ltd
Priority to CN202011371254.6A priority Critical patent/CN112505735B/zh
Publication of CN112505735A publication Critical patent/CN112505735A/zh
Application granted granted Critical
Publication of CN112505735B publication Critical patent/CN112505735B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42—Determining position
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/396—Determining accuracy or reliability of position or pseudorange measurements
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/40—Correcting position, velocity or attitude

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Navigation (AREA)

Abstract

本发明提供了对终端进行定位的方法、装置和存储介质。具体地，提供了利用终端上的地图应用程序对终端进行定位的方法。对终端进行定位的方法，包括：获取多个卫星在终端处的多普勒观测值；基于多普勒观测值计算终端的位置；基于计算出的终端的位置计算多个卫星分别到终端的伪距观测值；和基于伪距观测值和终端所在位置的气压值对计算出的终端的位置进行修正。该方法通过融合GNSS伪距观测值和终端的气压值，缩短了对终端地图用户首次定位的时间，提高了定位精度。

Description

对终端进行定位的方法、装置和存储介质
技术领域
本发明涉及终端定位的技术领域，具体地涉及融合了GNSS伪距观测值和气压值的对终端定位方法、装置和存储介质。
背景技术
目前终端（手机、车载终端等）上使用的定位软件（例如地图软件）主要是直接利用终端上的全球导航卫星系统GNSS芯片的定位结果。这种方法仅能够在与卫星信号完全同步后，即终端解算出卫星信号的发射时间后，才能对终端进行定位。在终端首次启动定位功能或GNSS信号较弱时，终端内置的GNSS芯片搜星较慢，无法完全同步卫星信号的发射时刻。这会导致终端无法进行定位或定位较慢。
发明内容
本公开提供了一种可以缓解、减轻或甚至消除上述问题中的一个或多个的用于对终端进行定位的方法、装置和存储介质。
根据本发明的一个方面，提供了一种对终端进行定位的方法。该方法包括：获取多个卫星在终端处的多普勒观测值；基于多普勒观测值计算终端的位置；基于计算出的终端的位置计算多个卫星分别到终端的伪距观测值；和基于伪距观测值和终端所在位置的气压值对计算出的终端的位置进行修正。
在一些实施例中，基于多普勒观测值计算终端的位置包括：根据基准站播发的实时导航星历计算得到多个卫星中每个卫星的位置、速度和卫星钟差变率，基于多普勒观测值、每个卫星的位置、速度和卫星钟差变率计算终端的位置。
在一些实施例中，基于多普勒观测值、每个卫星的位置、速度和卫星钟差变率计算终端的位置包括：
基于终端接收到的多个卫星的多普勒观测值
Figure DEST_PATH_IMAGE001
，n表示卫星的数量，和n颗卫星中的每颗卫星播发信号的波长
Figure 34573DEST_PATH_IMAGE002
，建立卫星多普勒观测值矩阵
Figure DEST_PATH_IMAGE003
，y表示当前时刻终端接收到的卫星多普勒观测值矩阵，
Figure 951713DEST_PATH_IMAGE004
表示第i颗卫星与终端之间的伪距变化率为
Figure DEST_PATH_IMAGE005
的观测值；
将伪距变化率
Figure 43035DEST_PATH_IMAGE005
的数学模型带入多普勒卫星观测矩阵中，得到卫星多普勒观测矩阵
Figure 500561DEST_PATH_IMAGE006
其中
Figure DEST_PATH_IMAGE007
，
Figure 314933DEST_PATH_IMAGE008
表示第i颗卫星的速度，
Figure DEST_PATH_IMAGE009
表示终端的位置，
Figure 950445DEST_PATH_IMAGE010
，
Figure 217478DEST_PATH_IMAGE008
表示卫星s的卫星钟差变率，c为真空中的光速值，
Figure DEST_PATH_IMAGE011
表示终端钟差变率；和
利用高斯-牛顿迭代法对卫星多普勒观测矩阵进行解算，输出计算出的终端的位置
Figure 478695DEST_PATH_IMAGE009
。
在一些实施例中，基于计算出的终端的位置计算多个卫星分别到终端的伪距观测值包括：基于计算出的终端的位置计算多个卫星分别到终端的伪距观测值的整数部分；基于接收卫星信号时刻
Figure 147574DEST_PATH_IMAGE012
和测量到的卫星信号发射时刻
Figure DEST_PATH_IMAGE013
计算多个卫星分别到终端的伪距观测值的小数部分；选择n颗卫星中高度角最大的卫星作为参考卫星，基于该高度角最大的卫星到终端的伪距观测值的整数部分和小数部分，对多个卫星分别到终端的伪距观测值的整数部分进行修正，以消除终端钟差的影响；根据经修正的多个卫星分别到终端的伪距观测值的整数部分、和计算出的多个卫星分别到终端的伪距观测值的小数部分计算多个卫星分别到终端的伪距观测值。
在一些实施例中，基于伪距观测值和终端所在位置的气压值对计算出的终端的位置进行修正包括：基于终端所在位置的气压值计算终端的高度；基于伪距观测值和终端的高度对所计算出的终端的位置进行修正。
在一些实施例中，基于终端所在位置的气压值计算终端的高度包括：基于
Figure 452522DEST_PATH_IMAGE014
计算终端的高度，其中终端所在位置的气压值为
Figure DEST_PATH_IMAGE015
，
Figure 206852DEST_PATH_IMAGE016
，
Figure DEST_PATH_IMAGE017
，
Figure 6181DEST_PATH_IMAGE018
，
Figure DEST_PATH_IMAGE019
，
Figure 342615DEST_PATH_IMAGE020
，g=9.8m/s2，
Figure DEST_PATH_IMAGE021
为计算得到的终端的高度。
在一些实施例中，基于伪距观测值和终端的高度对所计算出的终端的位置进行修正包括：基于伪距观测值的数学模型和终端的高度建立方程式：
Figure DEST_PATH_IMAGE023
其中
Figure 834776DEST_PATH_IMAGE024
和
Figure DEST_PATH_IMAGE025
分别为终端的纬度、经度和高度，
Figure 405564DEST_PATH_IMAGE026
为计算得到的终端的高度，利用高斯-牛顿迭代法对方程式进行解算，输出经修正的终端的位置。
在一些实施例中，伪距观测值的数学模型为
Figure DEST_PATH_IMAGE027
。其中
Figure 680688DEST_PATH_IMAGE028
为卫星i的伪距观测值，
Figure DEST_PATH_IMAGE029
为终端用户位置，
Figure 386476DEST_PATH_IMAGE030
为卫星i的位置，
Figure DEST_PATH_IMAGE031
为终端时钟偏差，
Figure 721642DEST_PATH_IMAGE032
为卫星i的钟差，c为真空中的光速值，
Figure DEST_PATH_IMAGE033
为误差改正数（包括电离层、对流层以及地球自转改正，可由经验模型计算得到），
Figure 263613DEST_PATH_IMAGE034
为卫星i的速度，
Figure DEST_PATH_IMAGE035
为卫星钟差变率可由导航星历计算得到，
Figure 404744DEST_PATH_IMAGE036
表示终端至卫星i的单位观测向量，
Figure DEST_PATH_IMAGE037
表示卫星i所在卫星导航系统与GPS的系统时间偏差，
Figure 902722DEST_PATH_IMAGE038
表示终端接收到卫星信号真实时刻与终端系统时间的偏差。
在一些实施例中，参数ISB (i)满足下述条件：响应于卫星i所在的卫星导航系统为全球定位系统GPS，
Figure DEST_PATH_IMAGE039
；响应于卫星i所在的卫星导航系统为北斗卫星导航系统，
Figure 985953DEST_PATH_IMAGE040
；响应于卫星i所在的卫星导航系统为GLONASS系统，ISB(i)=ISB(GLO)；和响应于卫星i所在的卫星导航系统为伽利略GALILEO系统，
Figure DEST_PATH_IMAGE041
。
在一些实施例中，终端可以为手持终端或车载终端。
在一些实施例中，利用高斯-牛顿迭代法对方程式进行解算，输出经修正的终端的位置包括：将待估计参数
Figure 998909DEST_PATH_IMAGE042
设置为
Figure DEST_PATH_IMAGE043
，其中
Figure 694463DEST_PATH_IMAGE009
为终端的位置，c为真空中的光速值，
Figure 46947DEST_PATH_IMAGE031
为终端时钟偏差，
Figure 989495DEST_PATH_IMAGE038
表示终端接收到卫星信号真实时刻与终端系统时间的偏差；计算参数
Figure 489747DEST_PATH_IMAGE042
的更正量；基于更正量对参数
Figure 910364DEST_PATH_IMAGE042
进行更正；和基于经更正的参数
Figure 382934DEST_PATH_IMAGE042
计算得到经修正的终端的位置。
根据本发明的一个方面，提供了一种对终端进行定位的装置。该装置包括：获取模块、计算模块、伪距观测值计算模块和修正模块。获取模块被配置成获取多个卫星在终端处的多普勒观测值。计算模块被配置成基于多普勒观测值计算终端的位置。伪距观测值计算模块被配置成基于计算出的终端的位置计算多个卫星分别到终端的伪距观测值。修正模块被配置成基于伪距观测值和终端所在位置的气压值对计算出的终端的位置进行修正。
在一些实施例中，修正模块被配置成基于伪距观测值和终端所在位置的气压值对计算出的终端的位置进行修正包括：基于终端所在位置的气压值计算终端的高度；基于伪距观测值和终端的高度对所计算出的终端的位置进行修正。
在一些实施例中，基于多普勒观测值计算终端的位置包括：根据基准站播发的实时导航星历计算得到多个卫星中每个卫星的位置、速度和卫星钟差变率，基于多普勒观测值、每个卫星的位置、速度和卫星钟差变率计算终端的位置。
在一些实施例中，基于计算出的终端的位置计算多个卫星分别到终端的伪距观测值包括：基于计算出的终端的位置计算多个卫星分别到终端的伪距观测值的整数部分；基于接收卫星信号时刻
Figure 807968DEST_PATH_IMAGE044
和测量到的卫星信号发射时刻
Figure 733198DEST_PATH_IMAGE013
计算多个卫星分别到终端的伪距观测值的小数部分；选择n颗卫星中高度角最大的卫星作为参考卫星，基于该高度角最大的卫星到终端的伪距观测值的整数部分和小数部分，对多个卫星分别到终端的伪距观测值的整数部分进行修正，以消除终端钟差的影响；根据经修正的多个卫星分别到终端的伪距观测值的整数部分、和计算出的多个卫星分别到终端的伪距观测值的小数部分计算多个卫星分别到终端的伪距观测值。
根据本发明的另一方面，提供了一种计算设备，包括：处理器；以及存储器，其上存储有指令，指令当在处理器上执行时促使处理器执行上面的方法中的任一个。
根据本发明的另一方面，提供了一种计算机可读存储介质，其上存储有计算机可读指令，计算机可读指令在被执行时实现上面的方法中的任一个。
通过本发明的实施例，提供了基于终端的GNSS伪距观测值和终端所在处的气压值实现对终端快速定位的方法。该方法在终端首次启动定位功能或GNSS卫星信号较弱、终端内置GNSS芯片搜星较慢的情况下，也能实现快速对终端定位，加快了终端用户首次定位的时间。与传统的仅仅依赖GNSS芯片定位的方法相比，定位更合理，提高了终端定位精度，优化了用户体验。
附图说明
在下面结合附图对于示例性实施例的描述中，本发明的更多细节、特征和优点被公开。附图仅用于示出实施方式的目的，而并不认为是对本公开的限制。而且在整个附图中，用相同的附图标记表示相同或相似的要素。在附图中：
图1描绘了GNSS导航芯片结构的原理示意图；
图2描绘了根据本发明一个实施例的终端定位示例场景的示意图；
图3描绘了根据本发明一个实施例的终端设备的示意图；
图4示例性描绘了根据本发明一个实施例的对终端进行定位的方法流程图；
图5描绘了根据本发明一个实施例的基站向终端播发星历的示意图；
图6示例性描绘了根据本发明一个实施例的对终端进行定位的装置框图；
图7a示意性示出了根据相关技术的终端定位结果的示意图；
图7b示意性示出了根据本发明一个实施例的终端定位结果的示意图；和
图8示出了能够实现根据本发明的一些实施例的用于对终端进行定位的方法的计算系统的示意框图。
具体实施方式
下面将参照附图更详细地描述本发明的若干个实施例以便使得本领域技术人员能够理解和实现本发明。然而，本发明可以具体实现为许多不同的形式并且不应被解释为局限于本文所阐述的实施例。相反，本文提供这些实施例以使得本发明全面且完整，并将充分地向本领域技术人员传达本发明的范围。这些实施例用于说明而不是限定本发明。
将理解的是，尽管术语第一、第二、第三等在本文中可以用来描述各种元件、步骤和/或部分，但是这些元件、步骤和/或部分不应当由这些术语限制。这些术语仅用来将一个元件、步骤或部分与另一个元件、步骤或部分相区分。因此，下面讨论的第一元件、步骤或部分可以被称为第二元件、步骤或部分而不偏离本发明的教导。
本文中使用的术语仅出于描述特定实施例的目的并且不意图限制本发明。如本文中使用的，单数形式“一个”、“一”和“该”意图也包括复数形式，除非上下文清楚地另有指示。将进一步理解的是，术语“包括”和/或“包含”当在本说明书中使用时指定特征、整体、步骤、操作、元件和/或部件的存在，但不排除一个或多个其他特征、整体、步骤、操作、元件、部件和/或其群组的存在。如本文中使用的，术语“和/或”包括相关联的列出项目中的一个或多个的任意和全部组合。
除非另有定义，本文中使用的所有术语（包括技术术语和科学术语）具有与本发明所属领域的普通技术人员所通常理解的相同含义。将进一步理解的是，诸如那些在通常使用的字典中定义的术语应当被解释为具有与其在相关领域和/或本说明书上下文中的含义相一致的含义，并且将不在理想化或过于正式的意义上进行解释，除非本文中明确地如此定义。
需要注意的是，在不冲突的前提下，实施例中的特征可以任意组合使用。
在详细介绍本发明的实施例之前，首先对一些相关的概念进行解释：
1、位置服务（LBS，Location Based Services）：位置服务是无线运营公司为用户提供的一种与位置有关的服务。位置服务利用各类型的定位技术来获取定位设备当前的所在位置，通过移动互联网向定位设备提供信息资源和基础服务。LBS服务中融合了移动通讯、互联网络、空间定位、位置信息、大数据等多种信息技术，利用移动互联网络服务平台进行数据更新和交互，使用户可以通过空间定位来获取相应的服务。
2、全球卫星导航系统（GNSS，Global Navigation Satellite System）：全球导航卫星系统是能在地球表面或近地空间的任何地点为用户提供全天候的3维坐标和速度以及时间信息的空基无线电导航定位系统。常见系统有GPS、BDS (BeiDou NavigationSatellite System)、GLONASS和GALILEO四大卫星导航系统。最早出现的是美国的GPS（Global Positioning System），现阶段技术最完善的也是GPS系统。随着近年来BDS、GLONASS系统在亚太地区的全面服务开启，尤其是BDS系统在民用领域发展越来越快。卫星导航系统已经在航空、航海、通信、人员跟踪、消费娱乐、测绘、授时、车辆监控管理和汽车导航与信息服务等方面广泛使用，而且总的发展趋势是为实时应用提供高精度服务。
3、手机GNSS导航芯片：手机导航芯片处理卫星信号并提供用户利用PVT算法估计得到的位置。PVT是基于芯片提供的原始观测值、实时导航星历以及其他信息计算得到的。GNSS导航芯片结构如图1所示。GNSS将通过天线101接收到的卫星信号发送给RF前端模块102进行射频处理，将RF前端模块处理后的信号发送给ADC模数转换模块103进行模拟信号和数字信号的变换。将经过模数变换后的信号发送给带通处理模块104进行频带的滤波，输出原始数据以及导航消息。最后，利用PVT计算模块结合外部传感器信息解算得出PVT统计信息。
4、安卓定位SDK：对于安卓操作系统，平台提供LBS的开发SDK，通过SDK提供的API可以获取得到用户当前的位置、定位精度以及其他相关信息，同时SDK还提供利用PVT的卫星原始观测值以及实时导航星历。
5、高斯一牛顿迭代法（Gauss-Newton iteration method）：是非线性回归模型中求回归参数进行最小二乘的一种迭代方法，该法使用泰勒级数展开式去近似地代替非线性回归模型，然后通过多次迭代，多次修正回归系数，使回归系数不断逼近非线性回归模型的最佳回归系数，最后使原模型的残差平方和达到最小。
图2描绘了根据本发明一个实施例的示例场景的示意图200。通常对终端进行定位需要来自多个卫星的卫星信号。图2中示意性示出了4颗卫星201。通常，对终端进行定位至少需要2颗卫星。如本领域技术人员所理解的，可以利用3颗或更多个卫星。卫星的数量越多，定位精度相应越高。终端设备包括智能手机202和智能车载终端203。终端设备可以是任何类型的移动计算设备，包括移动计算机（例如，Microsoft® Surface®设备、个人数字助理（PDA）、膝上型计算机、笔记本计算机、诸如Apple iPadTM的平板计算机、上网本等）、移动电话（例如，蜂窝电话、诸如Microsoft Windows®电话的智能手机、Apple iPhone、实现了Google® AndroidTM操作系统的电话、Palm®设备，Blackberry®设备等）、可穿戴式计算设备（例如智能手表、头戴式设备，包括智能眼镜，如Google® GlassTM，等）或其他类型的移动设备。在一些实施例中，终端设备也可以是固定式计算设备，例如台式计算机、游戏机等。此外，当媒体内容管理系统包括多个终端设备的情况下，该多个终端设备可以是相同或不同类型的计算设备。各个卫星201与终端设备之间的距离为伪距。
图3描绘了根据本发明一个实施例的终端设备300的示意图。参考图3，设备300具有处理器318。存储器320、显示器314、输入316和气压传感器322与处理器318连接。此外，GPS接收机324与处理器318连接。此外，用于接收GPS信号的天线326与GPS接收机324连接。电源（诸如电池或电池组（未示出））被用于向各种电子组件供应电力。此外，导航设备300可以包括端口，诸如串行数据端口，用于将设备300连接到远程处理器或个人计算机以便将信息（诸如地图信息）上载到设备300，或以便向远程处理器或个人计算机下载信息（诸如路线信息）。可替换地，设备300可以包括无线通信能力，使得数据从远处地点被接收。
图4示例性描绘了根据本发明一个实施例的对终端进行定位的方法400的流程图。该方法400通常是由位于终端（包括智能手机、车载终端等）上的处理器上执行的。如本领域技术人员所应理解的，该方法也可以通过将待处理的数据发送给第三方服务器而在服务器处执行。在步骤401中，获取多个卫星在终端处的多普勒观测值。由于当终端无法完全同步北斗/GPS/GLONASS/伽利略等卫星信号的发射时刻时，无法利用伪距进行GNSS定位。采用多普勒观测值进行用户粗略定位，计算用户粗略位置是为下一步骤的计算伪距观测值做准备。这里，将在t时刻卫星s的多普勒观测值记为
Figure DEST_PATH_IMAGE045
。
在步骤402中，基于多普勒观测值计算终端的位置。在一个实施例中，基于多普勒观测值计算终端的位置包括：根据基准站播发的实时导航星历计算得到多个卫星中每个卫星的位置、速度和卫星钟差变率，基于多普勒观测值、每个卫星的位置、速度和卫星钟差变率计算终端的位置。在另一个实施例中，基于多普勒观测值、每个卫星的位置、速度和卫星钟差变率计算终端的位置包括：
基于终端接收到的多个卫星的多普勒观测值
Figure 957506DEST_PATH_IMAGE001
，n表示卫星的数量，和n颗卫星中的每颗卫星播发信号的波长
Figure 81320DEST_PATH_IMAGE002
，建立卫星多普勒观测值矩阵
Figure 365671DEST_PATH_IMAGE003
，y表示当前时刻终端接收到的卫星多普勒观测值矩阵，
Figure 778198DEST_PATH_IMAGE004
表示第i颗卫星与终端之间的伪距变化率为
Figure 619246DEST_PATH_IMAGE005
的观测值；
将伪距变化率
Figure 800828DEST_PATH_IMAGE005
的数学模型带入多普勒卫星观测矩阵中，得到卫星多普勒观测矩阵
Figure DEST_PATH_IMAGE047
其中
Figure 52818DEST_PATH_IMAGE007
，
Figure 952641DEST_PATH_IMAGE008
表示第i颗卫星的速度，
Figure 518752DEST_PATH_IMAGE009
表示终端的位置，
Figure 132004DEST_PATH_IMAGE010
，
Figure 492579DEST_PATH_IMAGE008
表示卫星s的卫星钟差变率，c为真空中的光速值，
Figure 879698DEST_PATH_IMAGE011
表示终端钟差变率；
利用高斯-牛顿迭代法对卫星多普勒观测矩阵进行解算，输出计算出的终端的位置
Figure 311816DEST_PATH_IMAGE048
。
假设在t时刻卫星s的多普勒观测值为
Figure 202412DEST_PATH_IMAGE045
，卫星s的位置和运行速度为
Figure DEST_PATH_IMAGE049
、
Figure 812516DEST_PATH_IMAGE050
，用户的位置和运动速度为
Figure 421351DEST_PATH_IMAGE048
、
Figure DEST_PATH_IMAGE051
，卫星s和终端u之间的伪距变化率为
Figure 657161DEST_PATH_IMAGE052
，则有如下所示的数学模型：
Figure 667842DEST_PATH_IMAGE054
上其中
Figure DEST_PATH_IMAGE055
表示终端钟差变率，
Figure 953242DEST_PATH_IMAGE056
表示卫星s的卫星钟差变率，c为真空中的光速值，该式是用来估计用户的粗略位置的。
图5描绘了根据本发明一个实施例的基站向终端播发星历的示意图500。如图5中所示，卫星s的位置
Figure 49374DEST_PATH_IMAGE049
、速度
Figure DEST_PATH_IMAGE057
以及卫星钟差变率
Figure 823295DEST_PATH_IMAGE058
是根据基准站播发的实时导航星历计算得到的。基准站502会根据用户终端501（这里以智能手机为例）的请求播发实时导航星历，终端501根据接收到的实时导航星历计算出当前时刻（当前时刻可由移动终端的系统时间获取得到）卫星位置和速度。具体流程如下：
假设在t时刻，终端接收到n颗卫星的多普勒观测值
Figure DEST_PATH_IMAGE059
，则有
Figure DEST_PATH_IMAGE061
上其中
Figure 501532DEST_PATH_IMAGE062
表示卫星i播发信号的波长，y表示当前时刻终端接收到的卫星多普勒观测值矩阵。
Figure DEST_PATH_IMAGE063
表示参数
Figure 640389DEST_PATH_IMAGE052
的观测值。
由于终端用户速度相对于卫星s的速度较小，因此做如下假设：
Figure DEST_PATH_IMAGE065
为了尽量满足上述假设，本方案在提取多普勒观测值y矩阵时，剔除北斗卫星导航系统的地球静止轨道卫星相应的多普勒观测值。
将上述假设代入多普勒观测数学模型中，可得如下近似数学模型：
Figure DEST_PATH_IMAGE067
将参数
Figure 286134DEST_PATH_IMAGE052
的观测值代入上式可得
Figure DEST_PATH_IMAGE069
上述方程其中参数
Figure 847435DEST_PATH_IMAGE048
和
Figure 832708DEST_PATH_IMAGE058
为我们需要估计的参数：用户位置和终端接收机钟差变率，这里设为
Figure 204784DEST_PATH_IMAGE070
；
Figure DEST_PATH_IMAGE071
为卫星i的速度，
Figure 275508DEST_PATH_IMAGE072
为卫星i的位置，
Figure DEST_PATH_IMAGE073
为卫星i的卫星钟差变率，
Figure 141964DEST_PATH_IMAGE074
，它们都是由卫星实时导航星历计算得到的。
该方程式是非线性的，用最小二乘解算精度较差且数值不稳定，因此本方案采用高斯—牛顿迭代法解算上述方程式，具体步骤如下：
1）设置终端用户的初始位置，一般可设置为用户上次结束定位的位置；
2）第k迭代的估计参数为
Figure DEST_PATH_IMAGE075
，计算在
Figure 44061DEST_PATH_IMAGE075
上述方程式的偏导数，即
Figure DEST_PATH_IMAGE077
Figure DEST_PATH_IMAGE079
Figure DEST_PATH_IMAGE081
这里，符号
Figure 836305DEST_PATH_IMAGE082
表示叉乘运算。
3）计算权矩阵W，即
Figure DEST_PATH_IMAGE083
，
Figure 128746DEST_PATH_IMAGE084
，
Figure DEST_PATH_IMAGE085
为卫星i的多普勒观测值的信噪比。
4）计算估计参数的改正量
Figure 48161DEST_PATH_IMAGE086
：
Figure DEST_PATH_IMAGE087
。
5）基于计算得到的改正量更新估计参数：
Figure 555496DEST_PATH_IMAGE088
。
6）结束迭代计算直至满足下列条件，否则从步骤2）再次迭代计算；
Figure DEST_PATH_IMAGE089
。
7）输出所估计的终端位置。
在步骤403中，基于计算出的终端的位置计算多个卫星分别到终端的伪距观测值。在一个实施例中，基于计算出的终端的位置计算多个卫星分别到终端的伪距观测值包括：基于计算出的终端的位置计算多个卫星分别到终端的伪距观测值的整数部分；基于接收卫星信号时刻
Figure 207058DEST_PATH_IMAGE012
和测量到的卫星信号发射时刻
Figure 314691DEST_PATH_IMAGE013
计算多个卫星分别到终端的伪距观测值的小数部分；选择n颗卫星中高度角最大的卫星作为参考卫星，基于该高度角最大的卫星到终端的伪距观测值的整数部分和小数部分，对多个卫星分别到终端的伪距观测值的整数部分进行修正，以消除终端钟差的影响；根据经修正的多个卫星分别到终端的伪距观测值的整数部分、和计算出的多个卫星分别到终端的伪距观测值的小数部分计算多个卫星分别到终端的伪距观测值。
通常情况下，伪距观测值采用下式计算得到：
Figure DEST_PATH_IMAGE091
其中
Figure 21485DEST_PATH_IMAGE092
表示终端接收到卫星信号的时刻，
Figure DEST_PATH_IMAGE093
表示卫星信号在卫星端的发射时刻，参数
Figure 570278DEST_PATH_IMAGE092
和
Figure 720636DEST_PATH_IMAGE093
都是终端用户可获取的观测值，c为真空中的光速值。
由于终端无法完全同步北斗/GPS/GLONASS/伽利略等卫星信号的发射时刻，此时终端获取到的卫星信号的发射时刻是不完整的，即此时卫星信号发射时刻的小数部分是正确的，但是整数部分未知的。因此，不能直接利用
Figure 987670DEST_PATH_IMAGE094
的方法计算伪距观测值，本方案结合3.2.1获得的用户粗略位置，提出一种间接获取伪距观测值的有效方法。
在计算伪距观测值之前，先介绍一个概念：伪距观测值可以表示为光毫秒的单位，用公式表示为：
Figure 186570DEST_PATH_IMAGE096
上其中将伪距表示光毫秒为单位的，z的单位为ms，N为整数部分，
Figure DEST_PATH_IMAGE097
为小数部分，c为真空中的光速值。
1）计算伪距观测值整数部分N的估计值：
Figure DEST_PATH_IMAGE099
其中
Figure 402919DEST_PATH_IMAGE100
表示终端用户与卫星s的几何距离，
Figure DEST_PATH_IMAGE101
为卫星位置，
Figure 724178DEST_PATH_IMAGE009
为终端用户位置，是由步骤S402中获取得到的；
2）计算伪距观测值的小数部分
Figure 478508DEST_PATH_IMAGE097
的估计值：
Figure DEST_PATH_IMAGE103
其中
Figure 521245DEST_PATH_IMAGE012
为接收卫星信号时刻，
Figure 44630DEST_PATH_IMAGE013
为测量到的卫星信号发射时刻，
Figure 474474DEST_PATH_IMAGE104
表示
Figure DEST_PATH_IMAGE105
的整数部分，
Figure 512837DEST_PATH_IMAGE012
和
Figure 53540DEST_PATH_IMAGE013
均是可以直接从终端设备获取到的测量值。
3） 由步骤1）和2），得到n颗卫星的伪距观测值整数部分
4）
Figure 510060DEST_PATH_IMAGE106
和小数部分：
Figure DEST_PATH_IMAGE107
。
5）计算n颗卫星的高度角，选择高度角最高的卫星作为参考卫星，这里假设参考卫星为j。计算卫星i的伪距观测值的整数部分
Figure 907544DEST_PATH_IMAGE108
Figure 636465DEST_PATH_IMAGE110
其中
Figure DEST_PATH_IMAGE111
，
Figure 26864DEST_PATH_IMAGE112
，
Figure DEST_PATH_IMAGE113
和
Figure 524842DEST_PATH_IMAGE114
为卫星j和i的卫星钟差，可由实时导航星历计算得到。5）计算完整的伪距观测值：
Figure 30909DEST_PATH_IMAGE116
在步骤S404中，基于伪距观测值和终端所在位置的气压值对计算出的终端的位置进行修正。在一个实施例中，基于伪距观测值和终端所在位置的气压值对计算出的终端的位置进行修正包括：基于终端所在位置的气压值计算终端的高度；基于伪距观测值和终端的高度对所计算出的终端的位置进行修正。
在一个实施例中，基于终端所在位置的气压值计算终端的高度包括：基于
Figure 43865DEST_PATH_IMAGE014
计算终端的高度，其中终端所在位置的气压值为
Figure 926370DEST_PATH_IMAGE015
，
Figure 278854DEST_PATH_IMAGE016
，
Figure DEST_PATH_IMAGE117
，
Figure 300031DEST_PATH_IMAGE018
，
Figure 737965DEST_PATH_IMAGE019
，
Figure 220899DEST_PATH_IMAGE020
，g=9.8m/s2，
Figure 693469DEST_PATH_IMAGE021
为计算得到的终端的高度。
在一个实施例中，基于伪距观测值和终端的高度对所计算出的终端的位置进行修正包括：基于伪距观测值的数学模型和终端的高度建立方程式：
Figure 541339DEST_PATH_IMAGE118
，
其中
Figure 778154DEST_PATH_IMAGE024
和
Figure 2462DEST_PATH_IMAGE025
分别为终端的纬度、经度和高度，
Figure 329539DEST_PATH_IMAGE021
为计算得到的终端的高度，以及利用高斯-牛顿迭代法对方程式进行解算，输出经修正的终端的位置。
在一个实施例中，伪距观测值的数学模型为：
Figure 676206DEST_PATH_IMAGE027
，其中
Figure 88733DEST_PATH_IMAGE028
为卫星i的伪距观测值，
Figure 398623DEST_PATH_IMAGE029
为终端用户位置，
Figure 580205DEST_PATH_IMAGE030
为卫星i的位置，
Figure 35457DEST_PATH_IMAGE031
为终端时钟偏差，
Figure 263177DEST_PATH_IMAGE032
为卫星i的钟差，c为真空中的光速值，
Figure 829287DEST_PATH_IMAGE033
为误差改正数（包括电离层、对流层以及地球自转改正，可由经验模型计算得到），
Figure 865376DEST_PATH_IMAGE034
为卫星i的速度，
Figure 803114DEST_PATH_IMAGE035
为卫星钟差变率可由导航星历计算得到，
Figure 190233DEST_PATH_IMAGE036
表示终端至卫星i的单位观测向量，
Figure 294455DEST_PATH_IMAGE037
表示卫星i所在卫星导航系统与GPS的系统时间偏差，
Figure 247368DEST_PATH_IMAGE038
表示终端接收到卫星信号真实时刻与终端系统时间的偏差。
在一个实施例中，参数ISB (i)满足下述条件：响应于卫星i所在的卫星导航系统为全球定位系统GPS，
Figure 123051DEST_PATH_IMAGE039
；响应于卫星i所在的卫星导航系统为北斗卫星导航系统，
Figure 731887DEST_PATH_IMAGE040
；响应于卫星i所在的卫星导航系统为GLONASS系统，ISB(i)=ISB(GLO)；和响应于卫星i所在的卫星导航系统为伽利略GALILEO系统，
Figure 905379DEST_PATH_IMAGE041
。
在一个实施例中，利用高斯-牛顿迭代法对方程式进行解算，输出经修正的终端的位置包括：将待估计参数
Figure 712798DEST_PATH_IMAGE042
设置为
Figure 680754DEST_PATH_IMAGE043
，其中
Figure 776886DEST_PATH_IMAGE009
为终端的位置，c为真空中的光速值，
Figure 805934DEST_PATH_IMAGE031
为终端时钟偏差，
Figure 936701DEST_PATH_IMAGE038
表示终端接收到卫星信号真实时刻与终端系统时间的偏差；计算参数
Figure 75559DEST_PATH_IMAGE042
的更正量；基于更正量对参数
Figure 721304DEST_PATH_IMAGE042
进行更正；和基于经更正的参数
Figure 971019DEST_PATH_IMAGE042
计算得到经修正的终端的位置。
具体地，伪距观测值的数学模型为：
Figure 956293DEST_PATH_IMAGE120
其中
Figure 813522DEST_PATH_IMAGE028
为卫星i的伪距观测值，
Figure 884246DEST_PATH_IMAGE029
为终端用户位置，
Figure 937652DEST_PATH_IMAGE030
为卫星i的位置，
Figure 105329DEST_PATH_IMAGE031
为终端时钟偏差，
Figure 585988DEST_PATH_IMAGE032
为卫星i的钟差，c为真空中的光速值，
Figure 878429DEST_PATH_IMAGE033
为误差改正数（包括电离层、对流层以及地球自转改正，可由经验模型计算得到），
Figure 781532DEST_PATH_IMAGE034
为卫星i的速度，
Figure 475819DEST_PATH_IMAGE035
为卫星钟差变率可由导航星历计算得到，
Figure 127380DEST_PATH_IMAGE036
表示终端至卫星i的单位观测向量，
Figure DEST_PATH_IMAGE121
表示卫星i所在卫星导航系统与GPS的系统时间偏差：当i为GPS时，
Figure 235013DEST_PATH_IMAGE039
，当
Figure 630223DEST_PATH_IMAGE122
为北斗卫星导航系统时，
Figure 257644DEST_PATH_IMAGE040
，当i为GLONASS时，
Figure DEST_PATH_IMAGE123
，当i为GALILEO时，
Figure 345686DEST_PATH_IMAGE041
，而在本方案中，我们通过事先采集的GNSS伪距观测值标定
Figure 409457DEST_PATH_IMAGE124
和
Figure DEST_PATH_IMAGE125
，即
Figure 919941DEST_PATH_IMAGE126
可认为是已知的。
Figure 588820DEST_PATH_IMAGE038
表示终端接收到卫星信号真实时刻与终端系统时间的偏差。
n颗卫星组成如下方程式：
Figure 847763DEST_PATH_IMAGE128
移动终端的位置坐标是在ECEF坐标系下的（地心地固坐标系），即
Figure 133251DEST_PATH_IMAGE130
假设移动终端的纬度、经度和高度为
Figure 135842DEST_PATH_IMAGE024
和
Figure 659227DEST_PATH_IMAGE025
，则有如下关系式，即
Figure 636542DEST_PATH_IMAGE132
等式两边求微分，可得
Figure 878167DEST_PATH_IMAGE134
其中
Figure DEST_PATH_IMAGE135
为卯酉圈半径，
Figure 746766DEST_PATH_IMAGE136
为子午圈半径，可由下式计算得到：
Figure 124658DEST_PATH_IMAGE138
Figure 725403DEST_PATH_IMAGE140
Figure 234751DEST_PATH_IMAGE142
Figure 579145DEST_PATH_IMAGE144
其中
Figure DEST_PATH_IMAGE145
和b为地球椭球体的长半径和短半径，e为地球椭球偏心率；
由此可得
Figure DEST_PATH_IMAGE147
假设：
Figure DEST_PATH_IMAGE149
Figure DEST_PATH_IMAGE151
Figure DEST_PATH_IMAGE153
Figure DEST_PATH_IMAGE155
其中
Figure 952488DEST_PATH_IMAGE156
表示移动终端的真实位置，
Figure DEST_PATH_IMAGE157
表示移动终端的估计位置；
代入上式，则气压计高度测量值组成如下方程式：
Figure 520873DEST_PATH_IMAGE023
其中
Figure 737091DEST_PATH_IMAGE024
和
Figure 214338DEST_PATH_IMAGE025
为
Figure 566822DEST_PATH_IMAGE029
对应的纬度、经度和高度。
假设估计参数为
Figure 243791DEST_PATH_IMAGE043
，本方案采用高斯-牛顿迭代法解算上述方程式，具体步骤如下：
1）首先设置估计参数
Figure 478463DEST_PATH_IMAGE042
的初始值，
Figure 164659DEST_PATH_IMAGE029
设置为3.2.1中估计到的用户粗略位置，
Figure 371650DEST_PATH_IMAGE158
设置为零，
Figure 298149DEST_PATH_IMAGE038
设置为零；
2）第k迭代的估计参数为
Figure 223379DEST_PATH_IMAGE075
，计算在
Figure 447687DEST_PATH_IMAGE075
上述方程式的偏导数，即
Figure 837080DEST_PATH_IMAGE160
Figure DEST_PATH_IMAGE161
和
Figure 121431DEST_PATH_IMAGE162
为
Figure DEST_PATH_IMAGE163
对应的纬度、经度和高度
3）计算权矩阵
Figure 579963DEST_PATH_IMAGE164
：
Figure 607962DEST_PATH_IMAGE166
Figure 851861DEST_PATH_IMAGE168
Figure 307114DEST_PATH_IMAGE085
为卫星i的多普勒观测值的信噪比；
Figure DEST_PATH_IMAGE169
为气压计测量高度方差值，本方案中将其设置为
Figure 19986DEST_PATH_IMAGE170
m2；
4）计算估计参数的改正量
Figure DEST_PATH_IMAGE171
：
Figure DEST_PATH_IMAGE173
Figure 648413DEST_PATH_IMAGE174
、
Figure DEST_PATH_IMAGE175
和
Figure 684502DEST_PATH_IMAGE176
为第k次估计值；
5）更新估计参数：
Figure 622240DEST_PATH_IMAGE178
6）结束迭代计算直至满足下列条件，否则从步骤2）再次迭代计算；
Figure DEST_PATH_IMAGE179
；
7）通过更新的参数得到终端的位置。
该方法400在终端首次启动定位功能或GNSS卫星信号较弱、终端内置GNSS芯片搜星较慢的情况下，也能实现快速对终端定位，加快了终端用户首次定位的时间。与传统的仅仅依赖GNSS芯片定位的方法相比，定位更合理，提高了终端定位精度，优化了用户体验。
图6示例性描绘了根据本发明一个实施例的对终端进行定位的装置600的框图。装置600包括获取模块601、计算模块602、伪距观测值计算模块603和修正模块604。获取模块601被配置成获取多个卫星在终端处的多普勒观测值。计算模块602被配置成基于多普勒观测值计算终端的位置。伪距观测值计算模块603被配置成基于计算出的终端的位置计算多个卫星分别到终端的伪距观测值。修正模块604被配置成基于伪距观测值和终端所在位置的气压值对计算出的终端的位置进行修正。
该装置600在终端首次启动定位功能或GNSS卫星信号较弱、终端内置GNSS芯片搜星较慢的情况下，也能实现快速对终端定位，加快了终端用户首次定位的时间。与传统的仅仅依赖GNSS芯片定位的方法相比，定位更合理，提高了终端定位精度，优化了用户体验。
图7a示意性示出了根据相关技术的终端定位结果的示意图。图7b示意性示出了根据本发明一个实施例的终端定位结果的示意图。在该示例中，图7a和图7b是针对相同终端相同位置进行的定位。图7b中的定位比图7a中定位时间快了2s。图7b的定位在道路上，比定位在图7a中道路以外的区域更为合理。可见，本发明提出的解决智能终端定位速度较慢的有效方法，融合智能手机的GNSS伪距观测值和气压计实现用户快速定位，能够有效地辅助地图的定位技术，提高智能终端的定位精度，优化用户体验。特别地，提高了Android智能终端的定位精度。
图8示出了能够实现根据本发明的一些实施例的用于对终端进行定位的方法的计算系统800的示意框图。在一些实施例中，该计算系统800代表终端202或203中的计算系统。
计算系统800可以包括各种不同类型的设备，例如计算设备计算机、客户端设备、片上系统和/或任何其它合适的计算系统或计算系统。
计算系统800可以包括能够诸如通过系统总线814或其他适当的方式连接彼此通信的至少一个处理器802、存储器804、（多个）通信接口806、显示设备808、其他输入/输出（I/O）设备810以及一个或多个大容量存储装置812。
处理器802可以是单个处理单元或多个处理单元，所有处理单元可以包括单个或多个计算单元或者多个核心。处理器802可以被实施成一个或多个微处理器、微型计算机、微控制器、数字信号处理器、中央处理单元、状态机、逻辑电路和/或基于操作指令来操纵信号的任何设备。除了其他能力之外，处理器802可以被配置成获取并且执行存储在存储器804、大容量存储装置812或者其他计算机可读介质中的计算机可读指令，诸如操作系统816的程序代码、应用程序818的程序代码、其他程序820的程序代码等，以实现本发明实施例提供的方法。
存储器804和大容量存储设备812是用于存储指令的计算机存储介质的示例，指令由处理器802执行来实施前面所描述的各种功能。举例来说，存储器804一般可以包括易失性存储器和非易失性存储器二者（例如RAM、ROM等等）。此外，大容量存储设备812一般可以包括硬盘驱动器、固态驱动器、可移除介质、包括外部和可移除驱动器、存储器卡、闪存、软盘、光盘（例如CD、DVD）、存储阵列、网络附属存储、存储区域网等等。存储器804和大容量存储设备812在本文中都可以被统称为存储器或计算机存储介质，并且可以是能够把计算机可读、处理器可执行程序指令存储为计算机程序代码的非瞬时性介质，计算机程序代码可以由处理器802作为被配置成实施在本文的示例中所描述的操作和功能的特定机器来执行。
多个程序模块可以存储在大容量存储设备812上。这些程序包括操作系统816、一个或多个应用程序818、其他程序820和程序数据822，并且它们可以被加载到存储器804以供执行。这样的应用程序或程序模块的示例可以包括例如用于实现本文所提供的方法的计算机程序逻辑（例如，计算机程序代码或指令）。而且，这些程序模块可以分布在不同的物理位置，以实现相应的功能。例如，被描述为由终端202或203中的计算系统执行的方法可以分布在多个计算设备上来完成。
本发明还提供了一种计算机可读存储介质，其上存储有计算机可读指令，计算机可读指令在被执行时实现上述的方法。
虽然在图8中被图示成存储在计算系统800的存储器804中，但是模块814、818、820和822或者其部分可以使用可由计算系统800访问的任何形式的计算机可读介质来实施。如本文所使用的，“计算机可读介质”至少包括两种类型的计算机可读介质，也就是计算机存储介质和通信介质。
计算机存储介质包括通过用于存储信息的任何方法或技术实施的易失性和非易失性、可移除和不可移除介质，信息诸如是计算机可读指令、数据结构、程序模块或者其他数据。计算机存储介质包括而不限于RAM、ROM、EEPROM、闪存或其他存储器技术，CD-ROM、数字通用盘（DVD）、或其他光学存储装置，磁盒、磁带、磁盘存储装置或其他磁性存储设备，或者可以被用来存储信息以供计算系统访问的任何其他非传送介质。
与此相对，通信介质可以在诸如载波或其他传送机制之类的已调数据信号中具体实现计算机可读指令、数据结构、程序模块或其他数据。本文所定义的计算机存储介质不包括通信介质。
计算系统800还可以包括一个或更多通信接口806，以用于诸如通过网络、直接连接等等与其他设备交换数据。通信接口806可以促进在多种网络和协议类型内的通信，其中包括有线网络（例如LAN、电缆等等）和无线网络（例如WLAN、蜂窝、卫星等等）、因特网等等。通信接口806还可以提供与诸如存储阵列、网络附属存储、存储区域网等等中的外部存储装置（未示出）的通信。
在一些示例中，可以包括诸如监视器之类的显示设备808，以用于显示信息和图像。其他I/O设备810可以是接收来自用户的各种输入并且向用户提供各种输出的设备，并且可以包括触摸输入设备、手势输入设备、摄影机、键盘、遥控器、鼠标、打印机、音频输入/输出设备等等。
在本说明书的描述中，术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点被包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
流程图中或在此以其他方式描述的任何过程或方法描述可以被理解为，表示包括一个或多个用于实现定制逻辑功能或过程的步骤的可执行指令的代码的模块、片段或部分，并且本发明的优选实施方式的范围包括另外的实现，其中可以不按所示出或讨论的顺序（包括根据所涉及的功能按基本同时的方式或按相反的顺序）来执行功能，这应被本发明的实施例所属技术领域的技术人员所理解。
此外，在本发明各个实施例中的各功能单元可以集成在一个处理模块中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个模块中。上述集成的模块既可以采用硬件的形式实现，也可以采用软件功能模块的形式实现。集成的模块如果以软件功能模块的形式实现并作为独立的产品销售或使用时，也可以存储在一个计算机可读取存储介质中。
通过研究附图、公开内容和所附的权利要求书，本领域技术人员在实践所要求保护的主题时，能够理解和实现对于所公开的实施例的变型。在权利要求书中，词语“包括”不排除其他元件或步骤，并且不定冠词“一”或“一个”不排除多个。在相互不同的从属权利要求中记载了某些措施的仅有事实并不表明这些措施的组合不能用来获利。

Claims (15)

1.一种对终端进行定位的方法，包括：
获取多个卫星在所述终端处的多普勒观测值；
基于所述多普勒观测值计算所述终端的位置；
基于计算出的所述终端的位置计算所述多个卫星分别到所述终端的伪距观测值；和
基于所述伪距观测值和所述终端所在位置的气压值对计算出的所述终端的位置进行修正。
2.根据权利要求1所述的方法，其中所述基于所述多普勒观测值计算所述终端的位置包括：
根据基准站播发的实时导航星历计算得到所述多个卫星中每个卫星的位置、速度和卫星钟差变率，
基于所述多普勒观测值、所述每个卫星的位置、速度和卫星钟差变率计算所述终端的位置。
3.根据权利要求2所述的方法，其中所述基于所述多普勒观测值、所述每个卫星的位置、速度和卫星钟差变率计算所述终端的位置包括：
基于终端接收到的所述多个卫星的多普勒观测值
Figure DEST_PATH_IMAGE002
，n表示卫星的数量，和所述n颗卫星中的每颗卫星播发信号的波长
Figure DEST_PATH_IMAGE004
，建立卫星多普勒观测值矩阵
Figure DEST_PATH_IMAGE006
，y表示当前时刻终端接收到的卫星多普勒观测值矩阵，
Figure DEST_PATH_IMAGE008
表示第i颗卫星与所述终端之间的伪距变化率为
Figure DEST_PATH_IMAGE010
的观测值；
将伪距变化率
Figure 444073DEST_PATH_IMAGE010
的数学模型带入所述多普勒卫星观测矩阵中，得到卫星多普勒观测矩阵
Figure DEST_PATH_IMAGE012
其中
Figure DEST_PATH_IMAGE014
，
Figure DEST_PATH_IMAGE016
表示第i颗卫星的速度，
Figure DEST_PATH_IMAGE018
表示所述终端的位置，
Figure DEST_PATH_IMAGE020
，
Figure 797432DEST_PATH_IMAGE016
表示卫星s的卫星钟差变率，c为真空中的光速值，
Figure DEST_PATH_IMAGE022
表示终端钟差变率；和
利用高斯-牛顿迭代法对所述卫星多普勒观测矩阵进行解算，输出计算出的所述终端的位置
Figure 390218DEST_PATH_IMAGE018
。
4.根据权利要求1-3中任一项所述的方法，其中所述基于计算出的所述终端的位置计算所述多个卫星分别到所述终端的伪距观测值包括：
基于计算出的所述终端的位置计算所述多个卫星分别到所述终端的伪距观测值的整数部分；
基于接收卫星信号时刻
Figure DEST_PATH_IMAGE024
和测量到的卫星信号发射时刻
Figure DEST_PATH_IMAGE026
计算所述多个卫星分别到所述终端的伪距观测值的小数部分；
选择所述n颗卫星中高度角最大的卫星作为参考卫星，基于该高度角最大的卫星到所述终端的伪距观测值的整数部分和小数部分，对所述多个卫星分别到所述终端的伪距观测值的整数部分进行修正，以消除所述终端钟差的影响；
根据经修正的所述多个卫星分别到所述终端的伪距观测值的整数部分、和计算出的所述多个卫星分别到所述终端的伪距观测值的小数部分计算所述多个卫星分别到所述终端的伪距观测值。
5.根据权利要求1所述的方法，其中所述基于所述伪距观测值和所述终端所在位置的气压值对计算出的所述终端的位置进行修正包括：
基于所述终端所在位置的气压值计算所述终端的高度；
基于所述伪距观测值和所述终端的高度对所计算出的所述终端的位置进行修正。
6.根据权利要求5所述的方法，其中所述基于所述终端所在位置的气压值计算所述终端的高度包括：
基于
Figure DEST_PATH_IMAGE028
计算所述终端的高度，其中终端所在位置的气压值为
Figure DEST_PATH_IMAGE030
，
Figure DEST_PATH_IMAGE032
，
Figure DEST_PATH_IMAGE034
，
Figure DEST_PATH_IMAGE036
，
Figure DEST_PATH_IMAGE038
，
Figure DEST_PATH_IMAGE040
，g=9.8m/s2，
Figure DEST_PATH_IMAGE042
为计算得到的终端的高度。
7. 根据权利要求5所述的方法，其中所述基于所述伪距观测值和所述终端的高度对所计算出的所述终端的位置进行修正包括：
基于伪距观测值的数学模型和所述终端的高度建立方程式
Figure DEST_PATH_IMAGE044
其中
Figure DEST_PATH_IMAGE046
和
Figure DEST_PATH_IMAGE048
分别为所述终端的纬度、经度和高度，
Figure 414456DEST_PATH_IMAGE042
为计算得到的终端的高度，
利用高斯-牛顿迭代法对所述方程式进行解算，输出经修正的所述终端的位置。
8.根据权利要求7所述的方法，其中所述伪距观测值的数学模型为
Figure DEST_PATH_IMAGE050
其中
Figure DEST_PATH_IMAGE052
为卫星i的伪距观测值，
Figure DEST_PATH_IMAGE054
为终端用户位置，
Figure DEST_PATH_IMAGE056
为卫星i的位置，
Figure DEST_PATH_IMAGE058
为终端时钟偏差，
Figure DEST_PATH_IMAGE060
为卫星i的钟差，c为真空中的光速值，
Figure DEST_PATH_IMAGE062
为误差改正数，
Figure DEST_PATH_IMAGE064
为卫星i的速度，
Figure DEST_PATH_IMAGE066
为卫星钟差变率可由导航星历计算得到，
Figure DEST_PATH_IMAGE068
表示终端至卫星i的单位观测向量，
Figure DEST_PATH_IMAGE070
表示卫星i所在卫星导航系统与GPS的系统时间偏差，
Figure DEST_PATH_IMAGE072
表示终端接收到卫星信号真实时刻与终端系统时间的偏差。
9. 根据权利要求8所述的方法，其中参数ISB (i)满足下述条件：
响应于卫星i所在的卫星导航系统为全球定位系统GPS，
Figure DEST_PATH_IMAGE074
；
响应于卫星i所在的卫星导航系统为北斗卫星导航系统，
Figure DEST_PATH_IMAGE076
；
响应于卫星i所在的卫星导航系统为GLONASS系统，ISB(i)=ISB(GLO)；和
响应于卫星i所在的卫星导航系统为伽利略GALILEO系统，
Figure DEST_PATH_IMAGE078
。
10.根据权利要求7所述的方法，其中所述利用高斯-牛顿迭代法对所述方程式进行解算，输出经修正的所述终端的位置包括：
将待估计参数
Figure DEST_PATH_IMAGE080
设置为
Figure DEST_PATH_IMAGE082
，其中
Figure 353462DEST_PATH_IMAGE018
为所述终端的位置，c为真空中的光速值，
Figure 175924DEST_PATH_IMAGE058
为终端时钟偏差，
Figure 505274DEST_PATH_IMAGE072
表示终端接收到卫星信号真实时刻与终端系统时间的偏差；
计算参数
Figure 704175DEST_PATH_IMAGE080
的更正量；
基于所述更正量对参数
Figure 684638DEST_PATH_IMAGE080
进行更正；和
基于经更正的参数
Figure 678002DEST_PATH_IMAGE080
计算得到经修正的所述终端的位置。
11.如权利要求1-3中任一项所述的方法，其中所述终端可以为手持终端或车载终端。
12.一种对终端进行定位的装置，包括：
获取模块，被配置成获取多个卫星在所述终端处的多普勒观测值；
计算模块，被配置成基于所述多普勒观测值计算所述终端的位置；
伪距观测值计算模块，被配置成基于计算出的所述终端的位置计算所述多个卫星分别到所述终端的伪距观测值；和
修正模块，被配置成基于所述伪距观测值和所述终端所在位置的气压值对计算出的所述终端的位置进行修正。
13.根据权利要求12所述的装置，其中所述修正模块被配置成基于所述伪距观测值和所述终端所在位置的气压值对计算出的所述终端的位置进行修正包括：
基于所述终端所在位置的气压值计算所述终端的高度；
基于所述伪距观测值和所述终端的高度对所计算出的所述终端的位置进行修正。
14.一种计算设备，包括存储器和处理器，所述存储器被配置成在其上存储计算机可执行指令，所述计算机可执行指令当在所述处理器上执行时执行权利要求1-12中任一项所述的方法。
15.一种计算机可读存储介质，其上存储计算机可执行指令，所述计算机可执行指令当在处理器上执行时执行权利要求1-12中任一项所述的方法。
CN202011371254.6A 2020-11-30 2020-11-30 对终端进行定位的方法、装置和存储介质 Active CN112505735B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011371254.6A CN112505735B (zh) 2020-11-30 2020-11-30 对终端进行定位的方法、装置和存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011371254.6A CN112505735B (zh) 2020-11-30 2020-11-30 对终端进行定位的方法、装置和存储介质

Publications (2)

Publication Number Publication Date
CN112505735A true CN112505735A (zh) 2021-03-16
CN112505735B CN112505735B (zh) 2022-05-03

Family

ID=74967685

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011371254.6A Active CN112505735B (zh) 2020-11-30 2020-11-30 对终端进行定位的方法、装置和存储介质

Country Status (1)

Country Link
CN (1) CN112505735B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113031037A (zh) * 2021-05-25 2021-06-25 腾讯科技（深圳）有限公司 设备定位方法、装置、电子设备及计算机可读介质
CN115616640A (zh) * 2022-10-28 2023-01-17 北斗天地(北京)科技有限公司 一种智能手机gnss气压计组合定位方法
CN116227237A (zh) * 2023-05-08 2023-06-06 中国空气动力研究与发展中心空天技术研究所 一种航天器飞行中实时位置的精确迭代分析方法及系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4310884A (en) * 1979-05-09 1982-01-12 Chevron Research Company Microcomputer-controlled portable satellite signal simulator for exploration systems, and method of use
CN102486540A (zh) * 2010-12-06 2012-06-06 中国科学院微电子研究所 一种应用于全球卫星定位与导航系统中的快速定位方法
CN104865587A (zh) * 2015-06-15 2015-08-26 中国人民解放军国防科学技术大学 基于geo星座粗定位的北斗接收机快速定位方法
CN107167829A (zh) * 2017-04-06 2017-09-15 广州中硕云空间信息技术有限公司 一种基于mems压力传感器协同卫星定位的方法及系统
CN111158034A (zh) * 2020-01-15 2020-05-15 东方红卫星移动通信有限公司 一种基于低轨卫星多重覆盖场景下的快速定位方法
CN111736185A (zh) * 2020-05-28 2020-10-02 腾讯科技（深圳）有限公司 终端定位方法、装置及计算机可读存储介质和终端设备

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4310884A (en) * 1979-05-09 1982-01-12 Chevron Research Company Microcomputer-controlled portable satellite signal simulator for exploration systems, and method of use
CN102486540A (zh) * 2010-12-06 2012-06-06 中国科学院微电子研究所 一种应用于全球卫星定位与导航系统中的快速定位方法
CN104865587A (zh) * 2015-06-15 2015-08-26 中国人民解放军国防科学技术大学 基于geo星座粗定位的北斗接收机快速定位方法
CN107167829A (zh) * 2017-04-06 2017-09-15 广州中硕云空间信息技术有限公司 一种基于mems压力传感器协同卫星定位的方法及系统
CN111158034A (zh) * 2020-01-15 2020-05-15 东方红卫星移动通信有限公司 一种基于低轨卫星多重覆盖场景下的快速定位方法
CN111736185A (zh) * 2020-05-28 2020-10-02 腾讯科技（深圳）有限公司 终端定位方法、装置及计算机可读存储介质和终端设备

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113031037A (zh) * 2021-05-25 2021-06-25 腾讯科技（深圳）有限公司 设备定位方法、装置、电子设备及计算机可读介质
CN113031037B (zh) * 2021-05-25 2021-08-06 腾讯科技（深圳）有限公司 设备定位方法、装置、电子设备及计算机可读介质
CN115616640A (zh) * 2022-10-28 2023-01-17 北斗天地(北京)科技有限公司 一种智能手机gnss气压计组合定位方法
CN115616640B (zh) * 2022-10-28 2023-10-20 北斗天地(北京)科技有限公司 一种智能手机gnss气压计组合定位方法
CN116227237A (zh) * 2023-05-08 2023-06-06 中国空气动力研究与发展中心空天技术研究所 一种航天器飞行中实时位置的精确迭代分析方法及系统
CN116227237B (zh) * 2023-05-08 2023-07-21 中国空气动力研究与发展中心空天技术研究所 一种航天器飞行中实时位置的精确迭代分析方法及系统

Also Published As

Publication number Publication date
CN112505735B (zh) 2022-05-03

Similar Documents

Publication Publication Date Title
JP7267460B2 (ja) 高完全性衛星測位のためのシステムおよび方法
CN112505735B (zh) 对终端进行定位的方法、装置和存储介质
US11624843B2 (en) Systems and methods for reduced-outlier satellite positioning
US9851448B2 (en) Obtaining pseudorange information using a cellular device
CN103823222B (zh) 通过扩展sps轨道信息进行定位的方法和装置
CN111045034B (zh) 基于广播星历的gnss多系统实时精密时间传递方法及系统
US11885890B2 (en) Method and apparatus of single epoch position bound
JP6279078B2 (ja) 変換装置及びプログラム
WO2022156481A1 (zh) 星历预报方法和装置
CN110361692B (zh) 一种融合定位方法及装置
US11650327B2 (en) Antenna phase center compensation for orbital assistance data
WO2015145719A1 (ja) 測位装置
CN112285749B (zh) 全球导航卫星系统原始观测数据处理方法、装置及存储介质
JP2009025233A (ja) 搬送波位相式測位装置
US20110199257A1 (en) Method and system for updating altitude information for a location by using terrain model information to prime altitude sensors
CN114114334B (zh) Glonass频间偏差标定方法及rtk方法
CN113093250B (zh) 伪距观测数据修复方法、定位信息确定方法及装置
CN112505729B (zh) 确定卫星编号的方法、确定终端位置的方法、设备及介质
US9541650B2 (en) Satellite positioning method, satellite pseudorange calculation apparatus and satellite pseudorange calculation method
CN116577816B (zh) 利用双差定位确定用户站位置的方法、装置及存储介质
CN112987048B (zh) 一种适用于Andriod智能终端的高精度定位方法、智能终端以及存储介质
WO2023216766A1 (zh) 一种定位方法以及相关装置
TWI683122B (zh) 精密單點定位方法及其定位裝置與記錄媒體
CN110531395B (zh) 用于无人车定位的方法、装置和设备
CN117008175A (zh) 整周模糊度的确定方法、装置、电子设备及存储介质

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
REG Reference to a national code

Ref country code: HK

Ref legal event code: DE

Ref document number: 40040656

Country of ref document: HK

GR01 Patent grant
GR01 Patent grant