CN105527639A - 一种基于平滑与外推的卫星定位方法 - Google Patents
一种基于平滑与外推的卫星定位方法 Download PDFInfo
- Publication number
- CN105527639A CN105527639A CN201610052312.6A CN201610052312A CN105527639A CN 105527639 A CN105527639 A CN 105527639A CN 201610052312 A CN201610052312 A CN 201610052312A CN 105527639 A CN105527639 A CN 105527639A
- Authority
- CN
- China
- Prior art keywords
- elevation
- epoch
- extrapolation
- observed reading
- user terminal
- 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.)
- Pending
Links
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/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明涉及一种基于平滑与外推的卫星定位方法。所述方法利用当前历元以前解算出的若干个历元的用户终端所处位置的高程观测值进行外推,得到当前历元的高程预测值,同时解算当前历元的高程观测值,用于外推下一历元的高程预测值。将当前历元的高程预测值作为约束条件求解观测方程组,得到当前历元的用户终端所处位置的坐标。由于外推下一个历元得到的高程观测值比原单点定位得到的高程值更准确,提高了正常情况下的定位精度。当部分卫星信号短暂失锁或被遮挡,以最新的若干个历元的高程观测值的均值作为非完备定位条件下的高程观测值,高程观测值不稳定时,求加权均值,并进行上述相同的外推解算过程,可以完成非完备条件下的三维定位。
Description
技术领域
本发明涉及卫星定位与导航领域,具体涉及一种通过对用户终端所处位置的高程值(即大地高,用户终端至地球基准椭球面的高度)进行平滑与外推,实现非完备定位条件下的卫星定位方法。
背景技术
在城市、峡谷、密林等遮挡物较多的环境下,卫星定位***容易出现可视卫星少于4颗的“非完备定位条件”。此时若不依赖其他辅助手段或信息,用户终端将无法完成正常定位。在非完备条件结束后,仍需要一段时间才能重新捕获和锁定卫星信号进行定位。这种情况影响着卫星定位的连续性和可用性,成为GPS等卫星***在遮挡环境下应用的最大技术瓶颈之一。针对这一问题,常使用组合惯导、伪卫星***及电子地图匹配等方法对卫星星座进行补充,以实现三维定位。但这类通过外界手段辅助的方法,会增加投入与配套设施,且作用范围受到工作方式及布站位置的限制。相比之下,对用户终端坐标位置、测量伪距或接收机钟差进行数学建模处理的辅助定位方法无需额外的硬件设备,使用简单、经济方便。但其中利用位置坐标与伪距进行外推的方法属于矢量参数外推,当用户终端运动方向与速度大小发生改变时,容易产生较大误差。而基于接收机钟差外推的解决方法,因对接收机晶振的要求较高,难以广泛适用于市面上的一般民用接收终端。
另一方面,由于卫星定位***空间几何图形布局的原因,用户终端只能接收到来自天顶方向的卫星信号,而在地底方向没有卫星信号,所以位置精度因子PDOP(positiondilutionofprecision)值较大,从而导致高程方向上的定位误差较大。
发明内容
针对现有技术中存在的上述问题,本发明提出一种基于平滑与外推的卫星定位方法,用于正常定位情况下能够提高定位精度,并能完成非完备定位条件下的三维定位。
现有技术进行卫星定位的方法是:根据球面交会法解算用户终端所处位置的高程观测值(大地高),然后通过求解观测方程组得到用户终端的三维坐标。本发明提出的基于平滑与外推的卫星定位方法是对上述定位方法的改进,具体技术方案如下:
一种基于平滑与外推的卫星定位方法,包括以下步骤:
根据球面交会法解算当前历元n的用户终端所处位置的高程观测值hn,用于外推下一历元n+1的高程预测值hn+1预。对已解算出的历元n-1、历元n-2、……、历元n-m的高程观测值hn-1、hn-2、……、hn-m进行外推,得到当前历元n的高程预测值hn预,将当前历元n的高程预测值hn预作为约束条件求解观测方程组,得到当前历元n的用户终端所处位置的坐标。其中,n=1,2,……;m为外推需要的历元个数,其值由外推精度要求确定。
若部分卫星信号短暂失锁或被遮挡,即在非完备定位条件下,以k个根据球面交会法正常解算出的最新的高程观测值的均值(即进行均值平滑)作为非完备定位条件下的高程观测值hn。高程观测值不稳定时,所述均值为加权均值。
进一步地,对已解算出的高程观测值进行外推采用的模型包括:均值平滑模型,广义延拓插值模型,加权最小二乘模型。
进一步地,所述观测方程组为:
其中,x、y、z与Δt分别为待求的用户终端所处位置的三维坐标与接收机钟差,c为光速;xi、yi和zi为第i颗卫星在地固坐标系下的三维坐标,ρi为第i颗卫星至用户终端的伪距测量值,i=1,2,……;h为用户的高程值;a、b分别为地球基准椭球的长、短半轴。
进一步地,求加权均值时,高程观测值越新(历元序号越大),其权重越大。
进一步地,所述k满足:
式中,σh为不进行均值处理的高程坐标方向的均方误差,为均值处理后所期望的高程坐标方向的均方误差。
与现有技术相比,本发明具有以下有益效果:
本发明所述方法利用当前历元以前解算出的若干个历元的用户终端所处位置的高程观测值进行外推,得到当前历元的高程预测值,同时根据球面交会法解算当前历元的高程观测值,用于外推下一历元的高程预测值。将当前历元的高程预测值作为约束条件求解观测方程组,得到当前历元的用户终端所处位置的坐标。由于卫星***在高程方向上的定位误差较大,外推下一个历元得到的高程观测值比原单点定位得到的高程观测值更准确,因此本发明所述方法提高了正常情况下的定位精度。当部分卫星信号短暂失锁或被遮挡,即在非完备条件下,以正常定位解算出的最新的若干个历元的高程观测值的均值作为所述当前历元的高程观测值,高程观测值不稳定时,求加权均值,并进行上述相同的外推解算过程,可以完成非完备条件下的三维定位。
附图说明
图1为正常定位情况下的定位方法流程图。
具体实施方式
下面结合附图和实施例对本发明做进一步说明。
正常定位情况下定位方法流程图如图1所示,包括以下步骤:
根据球面交会法解算当前历元的用户终端所处位置的高程观测值,用于外推下一历元的高程预测值。利用当前历元以前根据球面交会法解算出的若干个历元的高程观测值(大地高)进行外推,得到当前历元的高程预测值。将当前历元的高程预测值作为约束条件求解观测方程组,得到当前历元的用户终端所处位置的坐标。外推需要的历元个数根据外推精度要求确定。
对已解算出的高程观测值序列进行外推采用的模型包括:均值平滑模型,广义延拓插值模型,加权最小二乘模型。
当用户终端所处位置的高程值h已知后,即相当于一颗卫星处于地球中心点附近,测得的高度值相当于增加了一段伪距值,成为以地心距为半径的另一个交会球,所以能够起到一颗导航卫星的作用。此时的观测方程为:
其中,x、y、z与Δt分别为待求的用户终端所处位置的三维坐标与接收机钟差,c为光速;xi、yi和zi为第i颗卫星在地固坐标系下的三维坐标,ρi为第i颗卫星至用户终端的伪距测量值,i=1,2,……;h为用户终端所处位置的高程值;a、b分别为地球基准椭球的长、短半轴。
上式中的第二个方程为地球椭球约束方程。因此,当h已知时,仅需要三个伪距观测方程即可实现三维定位。
求解上式的方法有多种,其中较为有效的一种算法是:地球椭球方程不直接参与方程组求解,而是单独列出,仅作为变量之间的约束条件。另一种方法是,将(1)式在用户近似位置处进行线性化,按泰勒级数展开为:
式中,cosαi、cosβi和cosγi(i=1,2,3)为第i颗卫星至用户终端的方向余弦;Δx、Δy、Δz与Δt分别为用户终端所处位置三维坐标与接收机钟差未知数的迭代修正值;Δρi(i=1,2,3)为第i颗卫星的伪距迭代修正值。
具体求解步骤如下:
(i)定义初始Δz=0,将用户终端估计位置处的和代入地球椭球方程,得
(ii)由(2)式求得Δx,Δy和Δt;
(iii)由迭代格式得到用户终端坐标x,y和接收机钟差Δt;
(iv)将x和y代入地球椭球方程得到z,并由式计算z坐标增量;
(v)重复上述步骤,直至Δx、Δy、Δz和Δt小于阈值时,迭代终止。
至此,即可解算得到用户终端三维位置坐标。
由于外推下一个历元得到的高程观测值比原单点定位得到的高程值更准确,因此采用本发明所述方法可以提高正常情况下的定位精度。
当部分卫星信号短暂失锁或被遮挡,可视卫星数不足4颗时无法实现定位解算。由于无论在步行还是车载定位中,室外环境下用户终端所处位置的高程值变化一般都比较稳定和平滑,特别是在一般城市环境中,很少发生剧烈、突然的变化。即使在车辆上下高架立交桥的过程中,其高程值的变化一般也在数米以内。所以,在上述非完备定位条件下,可以通过求若干个正常解算出的高程观测值的均值作为高程值,并按照前述外推方法进行处理。具体方法如下:
对根据球面交会法正常定位解算时已得到的最新的k个历元的高程观测值取均值,即进行均值平滑,并将所述均值作为非完备定位条件下的当前历元的高程观测值;高程观测值不稳定时,取所述k个历元高程观测值的加权均值作为非完备定位条件下的当前历元的高程观测值。加权均值的权重按照下面的原则选择:高程值越新(历元序号越大),权重越大。
k满足:
式中,σh为不进行均值处理的高程坐标方向的均方误差,为均值处理后所期望的高程坐标方向的均方误差。
例如,σh约为5m,如果希望达到0.5m,k值最少应该为100。
本发明不限于上述实施方式,本领域技术人员所做出的对上述实施方式任何显而易见的改进或变更,都不会超出本发明的构思和所附权利要求的保护范围。
Claims (5)
1.一种基于平滑与外推的卫星定位方法,其特征在于,所述方法包括以下步骤:
根据球面交会法解算当前历元n的用户终端所处位置的高程观测值hn,用于外推下一历元n+1的高程预测值hn+1预;对已解算出的历元n-1、历元n-2、……、历元n-m的高程观测值hn-1、hn-2、……、hn-m进行外推,得到当前历元n的高程预测值hn预,将当前历元n的高程预测值hn预作为约束条件求解观测方程组,得到当前历元n的用户终端所处位置的坐标;其中,n=1,2,……;m为外推需要的历元个数,其值由外推精度要求确定;
在非完备定位条件下,以k个根据球面交会法正常解算出的最新的高程观测值的均值作为非完备定位条件下的高程观测值hn;高程观测值不稳定时,所述均值为加权均值。
2.根据权利要求1所述的基于平滑与外推的卫星定位方法,其特征在于,对已解算出的高程观测值进行外推采用的模型包括:均值平滑模型,广义延拓插值模型,加权最小二乘模型。
3.根据权利要求1所述的基于平滑与外推的卫星定位方法,其特征在于,所述观测方程组为:
其中,x、y、z与Δt分别为待求的用户终端所处位置的三维坐标与接收机钟差,c为光速;xi、yi和zi为第i颗卫星在地固坐标系下的三维坐标,ρi为第i颗卫星至用户终端的伪距测量值,i=1,2,……;h为用户终端所处位置的高程值;a、b分别为地球基准椭球的长、短半轴。
4.根据权利要求1所述的基于平滑与外推的卫星定位方法,其特征在于,求加权均值时,高程观测值越新,即历元序号越大,其权重越大。
5.根据权利要求1所述的基于平滑与外推的卫星定位方法,其特征在于,所述k满足:
式中,σh为不进行均值处理的高程坐标方向的均方误差,为均值处理后所期望的高程坐标方向的均方误差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610052312.6A CN105527639A (zh) | 2016-01-26 | 2016-01-26 | 一种基于平滑与外推的卫星定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610052312.6A CN105527639A (zh) | 2016-01-26 | 2016-01-26 | 一种基于平滑与外推的卫星定位方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105527639A true CN105527639A (zh) | 2016-04-27 |
Family
ID=55769975
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610052312.6A Pending CN105527639A (zh) | 2016-01-26 | 2016-01-26 | 一种基于平滑与外推的卫星定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105527639A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109752743A (zh) * | 2019-01-17 | 2019-05-14 | 上海华测导航技术股份有限公司 | 卫星导航的二维定位方法 |
CN110109163A (zh) * | 2019-04-09 | 2019-08-09 | 天津大学 | 一种附有高程约束的精密单点定位方法 |
CN111679296A (zh) * | 2020-05-28 | 2020-09-18 | 哈尔滨工业大学 | 一种基于gp迭代外推的单星定位方法 |
CN116893433A (zh) * | 2023-09-08 | 2023-10-17 | 真点科技(北京)有限公司 | 一种实现跟踪站观测值预测的方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110298658A1 (en) * | 2010-06-02 | 2011-12-08 | Qualcomm Incorporated | Position Determination Using Measurements From Past And Present Epochs |
CN103968836A (zh) * | 2014-05-16 | 2014-08-06 | 施浒立 | 一种基于时序伪距差分推算移动目标位置的方法及装置 |
CN104035113A (zh) * | 2014-06-26 | 2014-09-10 | 桂林电子科技大学 | 基于伪距的多模gnss接收机的可靠定位方法 |
CN104280756A (zh) * | 2014-10-30 | 2015-01-14 | 中国科学院国家天文台 | 基于接收机钟差广义延拓逼近法的卫星定位增强方法 |
-
2016
- 2016-01-26 CN CN201610052312.6A patent/CN105527639A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110298658A1 (en) * | 2010-06-02 | 2011-12-08 | Qualcomm Incorporated | Position Determination Using Measurements From Past And Present Epochs |
CN103968836A (zh) * | 2014-05-16 | 2014-08-06 | 施浒立 | 一种基于时序伪距差分推算移动目标位置的方法及装置 |
CN104035113A (zh) * | 2014-06-26 | 2014-09-10 | 桂林电子科技大学 | 基于伪距的多模gnss接收机的可靠定位方法 |
CN104280756A (zh) * | 2014-10-30 | 2015-01-14 | 中国科学院国家天文台 | 基于接收机钟差广义延拓逼近法的卫星定位增强方法 |
Non-Patent Citations (2)
Title |
---|
刘成等: "联合气压测高定位的高程转换问题", 《测绘科学》 * |
杜晓辉等: "基于地面移动通信基站的差分气压测高方法", 《北京航空航天大学学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109752743A (zh) * | 2019-01-17 | 2019-05-14 | 上海华测导航技术股份有限公司 | 卫星导航的二维定位方法 |
CN110109163A (zh) * | 2019-04-09 | 2019-08-09 | 天津大学 | 一种附有高程约束的精密单点定位方法 |
CN110109163B (zh) * | 2019-04-09 | 2022-01-07 | 天津大学 | 一种附有高程约束的精密单点定位方法 |
CN111679296A (zh) * | 2020-05-28 | 2020-09-18 | 哈尔滨工业大学 | 一种基于gp迭代外推的单星定位方法 |
CN111679296B (zh) * | 2020-05-28 | 2023-04-25 | 哈尔滨工业大学 | 一种基于gp迭代外推的单星定位方法 |
CN116893433A (zh) * | 2023-09-08 | 2023-10-17 | 真点科技(北京)有限公司 | 一种实现跟踪站观测值预测的方法及装置 |
CN116893433B (zh) * | 2023-09-08 | 2023-12-05 | 真点科技(北京)有限公司 | 一种实现跟踪站观测值预测的方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110487301A (zh) | 一种雷达辅助机载捷联惯性导航***初始对准方法 | |
CN102928858B (zh) | 基于改进扩展卡尔曼滤波的gnss单点动态定位方法 | |
CN107710017A (zh) | 用于在实时运动模式和相对定位模式之间切换的卫星导航接收器及方法 | |
CN102998690B (zh) | 一种基于gps载波双差方程的姿态角直接求解方法 | |
Georgy et al. | Vehicle navigator using a mixture particle filter for inertial sensors/odometer/map data/GPS integration | |
CN105891863B (zh) | 一种基于高度约束的扩展卡尔曼滤波定位方法 | |
CN105849589A (zh) | 全球导航卫星***、定位终端、定位方法以及记录介质 | |
Mahmoud et al. | Integrated positioning for connected vehicles | |
Elsheikh et al. | Integration of GNSS precise point positioning and reduced inertial sensor system for lane-level car navigation | |
CN103033822B (zh) | 移动信息确定装置、方法以及接收机 | |
US20230403675A1 (en) | Positioning Method, Apparatus, Device and System, and Storage Medium | |
CN105527639A (zh) | 一种基于平滑与外推的卫星定位方法 | |
Yi et al. | Tightly-coupled GPS/INS integration using unscented Kalman filter and particle filter | |
CN105044738A (zh) | 一种接收机自主完好性监视的预测方法及预测*** | |
Dawson et al. | Radar-based multisensor fusion for uninterrupted reliable positioning in GNSS-denied environments | |
CN109613582A (zh) | 一种车载实时单频米级伪距定位方法 | |
Lee et al. | Seamless accurate positioning in deep urban area based on mode switching between DGNSS and multipath mitigation positioning | |
Ye et al. | An improved single-epoch GNSS/INS positioning method for urban canyon environment based on real-time DISB estimation | |
CN104280756A (zh) | 基于接收机钟差广义延拓逼近法的卫星定位增强方法 | |
Li et al. | An indoor and outdoor seamless positioning system for low-cost UGV using PPP/INS/UWB tightly coupled integration | |
Amt et al. | Positioning for range-based land navigation systems using surface topography | |
Georgy et al. | Nonlinear filtering for tightly coupled RISS/GPS integration | |
CN111103603B (zh) | 基于cors***的云定位方法及装置、定位***及云端服务器 | |
Dabove et al. | Mass-market L1 GPS receivers for mobile mapping applications: A novel approach | |
Xu et al. | Multi-Sensor and Analytical Constraints Tightly Augmented BDS-3 RTK for Vehicle-Borne Positioning |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160427 |