CN110873887A - Gps和bds紧结合的定位方法及可读存储介质 - Google Patents
Gps和bds紧结合的定位方法及可读存储介质 Download PDFInfo
- Publication number
- CN110873887A CN110873887A CN201911387743.8A CN201911387743A CN110873887A CN 110873887 A CN110873887 A CN 110873887A CN 201911387743 A CN201911387743 A CN 201911387743A CN 110873887 A CN110873887 A CN 110873887A
- Authority
- CN
- China
- Prior art keywords
- gps
- bds
- satellite
- difference
- double
- 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
- G01S19/45—Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement
- G01S19/46—Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement the supplementary measurement being of a radio-wave signal type
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)
Abstract
本公开涉及卫星定位技术领域,公开了一种GPS和BDS紧结合的定位方法,包括以下步骤:选择参考星,分别在系统内和系统间做双差,组GPS和BDS紧结合双差观测方程模型;根据GPS和BDS紧结合双差观测方程模型,采用卡尔曼滤波法进行状态参数估计。本公开的一些技术效果在于:实现不同系统之间的优化融合处理,可提高定位精度、可用性和可靠性。
Description
技术领域
本公开涉及卫星定位技术领域,特别涉及一种GPS和BDS紧结合的定位方法及可读存储介质。
背景技术
随着目前GPS、BDS等系统建设的稳步推进,多GNSS数据处理方法及应用的研究逐渐成为卫星导航领域的研究热点.已有结果表明,综合多个GNSS系统的信号并在观测值层面统一处理能有效提高GNSS定位的可靠性和稳定性,尤其在山区、城市峡谷等信号遮挡严重的地区,可充分有效利用原始观测值。
但是,当前多GNSS融合相对定位,尤其GPS和BDS融合定位,通常只在各个卫星系统内部组双差观测值,然后再联合不同的系统进行组合定位解算,虽然在一般条件下能够达到精度要求,但是存在以下缺陷:①在车载导航特别是树木遮挡较为严重等应用场景下,定位精度无法达到令人满意的精度;②没有充分利用多平多系统的资源和尽可能多的使用观测到的信息,以达到实现不同系统之间的最优融合处理,提高定位精度、可用性和可靠性。
发明内容
为主要解决GPS和BDS融合定位的技术问题,本公开在一方面提出了GPS和BDS紧结合的定位方法,包括以下步骤:选择参考星,分别在系统内和系统间做双差,组GPS和BDS紧结合双差观测方程模型;根据GPS和BDS紧结合双差观测方程模型,采用卡尔曼滤波法进行状态参数估计。
优选地,所述的选择参考星执行如下选择策略:当前历元的GPS和BDS所有卫星,①在没有发生周跳的卫星当中选取高度角最大的卫星作为参考星;②如果所有卫星均发生周跳,则在所有卫星当中选择高度角最大的卫星作为参考星。
优选地,所述的GPS和BDS紧结合双差观测方程为:
其中,各符号的含义:表示双差算子,Δ表示单差算子,S、K表示GPS或BDS系统,i和j表示S和K的频率,P表示伪距观测值,Φ表示相位观测值,dtr、dtS分别表示接收机钟差和卫星钟差,I、T分别表示电离层延迟和对流层延迟,dr、dS分别表示接收机和卫星的伪距硬件延迟,δr、δS分别表示接收机和卫星的载波相位硬件延迟,分别表示接收机和卫星的载波初始相位偏差,ε、分别表示伪距和载波的多路径、噪声以及未模型化的误差。
优选地,所述的GPS和BDS紧结合双差观测方程:在短基线相对定位中时,
①相同频率,同一卫星导航系统:
②相同频率,不同卫星导航系统:
③不同频率,不同卫星导航系统:
优选地,所述的GPS和BDS紧结合双差观测方程:在不同频率组双差时差分系统间偏差计算执行步骤:
优选地,所述的GPS和BDS紧结合双差观测方程:在不同频率组双差时模糊度计算执行步骤:
即在不同频率时,GPS和BDS紧结合的模糊度方程为:
优选地,所述的状态参数估计:设当前历元观测到mG个GPS卫星和nC个BDS卫星,假设参考星为1G,在GPS系统内和GPS/BDS系统间作双差,则状态参数方程如下:
将观测方程线性化之后,以矩阵的形式列出如下:
优选地,所述的状态参数包括:位置参数、差分系统间偏差参数、站间单差模糊度参数。
在另一些方面,本公开提出了一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现所述的GPS和BDS紧结合的定位方法的步骤。
本公开的一些技术效果在于:提供了GPS和BDS紧结合的定位方法,实现不同系统之间的优化融合处理,可提高定位精度、可用性和可靠性。
具体实施方式
下文将对本公开涉及的技术手段或技术效果作进一步的展开描述,显然,所提供的实施例(或实施方式)仅是本公开意旨涵盖的部分实施方式,而并非全部。基于本公开中的实施例的明示或暗示,本领域技术人员在没有作出创造性劳动的前提下所能获得的所有其他实施例,都将在本公开保护的范围之内。
在总体思路上,本公开提出了一种GPS和BDS紧结合的定位方法,包括以下步骤:选择参考星,分别在系统内和系统间做双差,组GPS和BDS紧结合双差观测方程模型;根据GPS和BDS紧结合双差观测方程模型,采用卡尔曼滤波法进行状态参数估计。
本公开提出的GPS和BDS紧结合的定位方法,特点之一是作双差时共用参考星,在系统内和系统间作双差,可以组建系统内和系统间的双差观测值,减少了参考卫星数,增加了额外观测值,提高了定位模型强度,特别是对于车载导航、自动驾驶等观测环境遮挡严重、可用卫星数较少的条件下,松组合模型的性能受到限制,此时可以显著提升浮点模糊度的精度、相关性,从而提高定位的精度和可靠性。
本公开提出的GPS和BDS紧结合的定位方法,可普遍适用于GPS和BDS的融合定位;提出的模型方法和处理过程可使用单频或双频观测值进行定位,有利于降低设备成本;实时选择参考星,避免固定GEO卫星为参考星产生的定位偏差。
具体的,在一些实施例中,选择参考星,分别在系统内和系统间做双差,组GPS和BDS紧结合双差观测方程模型;根据GPS和BDS紧结合双差观测方程模型,采用卡尔曼滤波法进行状态参数估计。
在一些实施例中,选择参考星执行如下选择策略:当前历元的GPS和BDS所有卫星,①在没有发生周跳的卫星当中选取高度角最大的卫星作为参考星;②如果所有卫星均发生周跳,则在所有卫星当中选择高度角最大的卫星作为参考星。传统方法是分系统来分别选择参考星,或固定GEO卫星为多系统的共同参考星,本公开根据每一历元的搜星情况,对当前历元的所有卫星进行分析,依据周跳、高度角因素实时动态的选择GPS和BDS的共同参考星,从而避免参考星选择不当产生的定位偏差。
在一些实施例中,GPS、BDS紧组合相对定位函数模型GNSS非差伪距和载波相位观测方程如下:
在系统内和系统间做差,组GPS和BDS紧组合双差观测方程数学模型GPS和BDS紧结合双差观测方程为:
其中,各符号的含义:表示双差算子,Δ表示单差算子,S、K表示GPS或BDS系统,i和j表示S和K的频率,P表示伪距观测值,Φ表示相位观测值,dtr、dtS分别表示接收机钟差和卫星钟差,I、T分别表示电离层延迟和对流层延迟,dr、dS分别表示接收机和卫星的伪距硬件延迟,δr、δS分别表示接收机和卫星的载波相位硬件延迟,分别表示接收机和卫星的载波初始相位偏差.ε、分别表示伪距和载波的多路径、噪声以及未模型化的误差。
在一些实施例中,GPS和BDS紧结合双差观测方程:在短基线相对定位中时,GPS与BDS短基线紧组合相对定位模型在短基线中,电离层延迟和对流层延迟可以忽略不计,GPS与BDS短基线紧组合相对定位模型总共分三种组合情况,各组合模型如下,其中符号定义与前述相同。
a)相同频率,同一卫星导航系统
设有卫星s和卫星k为同一卫星导航系统卫星,且信号频率相同,即为常规的系统内双差方程,联合多个系统内的双差方程即为目前广泛采用的松组合双差模型(LCM):
b)相同频率,不同卫星导航系统
对于GPS和BDS(主要是北斗3号系统卫星)两个不同的卫星导航系统卫星,北斗3号卫星与GPS卫星有重叠频率的信号,即为不同系统间的重叠频率的紧组合双差模型(TCM)
c)不同频率,不同卫星导航系统
可将GPS与BDS不同频率的双差观测方程转换成如下形式:
再对其进行变换,可以将GPS与BDS不重叠的频率的方程转换成重叠频率的形式的方程:
d)组方程组
假设观测了mG个GPS卫星和nC个BDS卫星,假设参考星为1G,在GPS L1系统内和GPSL1/BDS B1系统间作双差,则双差观测方程如下式所示。
状态参数为:
将观测方程线性化之后,以矩阵的形式列出如下:
在一些实施例GPS和BDS紧结合双差观测方程:在不同频率组双差时差分系统间偏差计算执行步骤:
当由于当前历元的GPS和BDS的卫星频率不同,两者组双差时接收机载波初始相位和硬件延迟无法分离,需要进行参数重组:
在一些实施例中,GPS和BDS紧结合双差观测方程:在不同频率组双差时模糊度计算执行步骤:
为了保证单差模糊度的整数特性,必须从双差观测值中扣除单差偏差项的影响,同时单差偏差项需要满足一定的精度,不影响双差模糊度的固定,即:
在一些实施例中,当前历元的状态参数估计:设当前历元观测到mG个GPS卫星和nC个BDS卫星,假设参考星为1G,在GPS系统内和GPS/BDS系统间作双差,则状态参数方程如下:
将观测方程线性化之后,以矩阵的形式列出如下:
在一些实施例中,GPS与BDS紧组合参数包括位置参数、差分系统间偏差参数、各系统各卫星的站间单差模糊度参数,通常用序惯最小二乘或卡尔曼滤波算法进行参数估计,本公开采用扩展的卡尔曼滤波算法进行参数估计。状态方程和观测方程表示如下:
Xk=Φk,k-1Xk-1+Γk-1wk
Lk=HkXk+vk
上式中,Xk、Xk-1分别表示第k和k-1个历元时刻的状态向量,Φk,k-1表示从k-1历元到第k历元时刻的状态转移矩阵,Γk-1为系统噪声驱动矩阵;Lk为第k历元时刻的观测向量,Hk为对应的系数矩阵,wk和vk分别表示系统状态噪声和观测噪声向量,二者为互不相关的零均值白噪声,对应的随机模型如下:
E(wk)=0,Cov{wk,wj}=Qkδkj
E(vk)=0,Cov{vk,vj}=Rkδkj
Cov{wk,vj}=0
其中Qk、Rk分别表示系统状态向量的方差阵和观测噪声的方差阵,前者也称为卡尔曼滤波的过程噪声。δkj为kronecker函数,定义如下:
由状态更新的结果作为初值,联立观测值组成观测方程按最小方差估计准则得到参数的修正结果,也即测量更新:
Pk=(I-KkHk)Pk,k-1
上式中,Pk,k-1分别表示状态更新的向量及其方差协方差阵,I为单位矩阵,Kk为增益矩阵,它通过状态方差协方差阵和观测方差协方差阵来动态调节状态信息和观测信息对当前历元状态估值的影响,本质上反映的是二者的权比关系。
通常将卡尔曼滤波分为上述的两个计算过程:时间更新(状态更新、一步预测),即上一个历元的状态向量和方差协方差阵通过状态转移矩阵预测得到当前历元的状态向量和方差协方差阵;测量更新,即利用当前历元的观测信息对预测的状态向量及其方差协方差阵进行修正。对应的递推流程如下:对于模糊度固定的处理策略则选择部分模糊度固定,分别从卫星层面、频率层面和模糊度层面等角度选取模糊度子集,对其进行尝试模糊度固定,由固定的模糊度得到固定解,再对固定解进行验证,通过则输出固定解,否则输出浮点解;
在一些实施例中,对下一个历元判断是否更换参考星,判断是否需要重新估计差分系统间偏差参数。如果参考星不变,那么不重新估计差分系统间偏差参数,将上一个历元的相应数值直接传递到下一个历元,如果参考星发生了变化,则重置差分系统间偏差参数,让其重新滤波收敛。
在一些实施例中,状态参数包括:位置参数、差分系统间偏差参数、站间单差模糊度参数。
在一些方面,本公开提出了一种计算机可读存储介质,其上存储有计算机程序,其特征在于:所述计算机程序被处理器执行时实现前述的动态加权方法的步骤。本领域技术人员可以理解的是,实施例中的全部或部分步骤,可以通过计算机程序来指令相关的硬件实现,该程序可以存储于计算机可读介质中,可读介质可以包括闪存盘、移动硬盘、只读存储器、随机存取器、磁盘或光盘等各种可以存储程序代码的介质。
在符合本领域技术人员的知识和能力水平范围内,本文提及的各种实施例或者技术特征在不冲突的情况下,可以相互组合而作为另外一些可选实施例,这些并未被一一罗列出来的、由有限数量的技术特征组合形成的有限数量的可选实施例,仍属于本公开揭露的技术范围内,亦是本领域技术人员结合本文内容所能理解或推断而得出的。
另外,多数实施例的描述是基于不同的重点而展开的,未详述之处,可参见现有技术的内容或本文的其他相关描述进行理解。
再次强调,上文所列举的实施例,为本公开较为典型的、较佳实施例,仅用于详细说明、解释本公开的技术方案,以便于读者理解,并不用以限制本公开的保护范围或者应用。在本公开的精神和原则之内所作的任何修改、等同替换、改进等而获得的技术方案,都应被涵盖在本公开的保护范围之内。
Claims (9)
1.GPS和BDS紧结合的定位方法,其特征在于,包括以下步骤:
选择参考星,分别在系统内和系统间做双差,组GPS和BDS紧结合双差观测方程模型;
根据GPS和BDS紧结合双差观测方程模型,采用卡尔曼滤波法进行状态参数估计。
2.如权利要求1所述的定位方法,其特征在于:所述的选择参考星执行如下选择策略:当前历元的GPS和BDS所有卫星,
①在没有发生周跳的卫星当中选取高度角最大的卫星作为参考星;
②如果所有卫星均发生周跳,则在所有卫星当中选择高度角最大的卫星作为参考星。
8.如权利要求7所述的定位方法,其特征在于,所述的状态参数包括:位置参数、差分系统间偏差参数、站间单差模糊度参数。
9.计算机可读存储介质,其上存储有计算机程序,其特征在于:所述计算机程序被处理器执行时实现如权利要求1~8中任一所述的GPS和BDS紧结合的定位方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911387743.8A CN110873887A (zh) | 2019-12-26 | 2019-12-26 | Gps和bds紧结合的定位方法及可读存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911387743.8A CN110873887A (zh) | 2019-12-26 | 2019-12-26 | Gps和bds紧结合的定位方法及可读存储介质 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110873887A true CN110873887A (zh) | 2020-03-10 |
Family
ID=69718392
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911387743.8A Pending CN110873887A (zh) | 2019-12-26 | 2019-12-26 | Gps和bds紧结合的定位方法及可读存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110873887A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115480274A (zh) * | 2021-05-31 | 2022-12-16 | 千寻位置网络有限公司 | 一种参考星选择方法、装置及接收机 |
CN116106953A (zh) * | 2023-04-14 | 2023-05-12 | 中南大学 | 一种gnss相对定位半参数估计的时间差分载波相位增强方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5451964A (en) * | 1994-07-29 | 1995-09-19 | Del Norte Technology, Inc. | Method and system for resolving double difference GPS carrier phase integer ambiguity utilizing decentralized Kalman filters |
CN106569242A (zh) * | 2016-11-04 | 2017-04-19 | 中国科学院测量与地球物理研究所 | 固定参考星的gnss单差处理方法 |
CN108519614A (zh) * | 2018-03-16 | 2018-09-11 | 东南大学 | 一种gps/bds紧组合载波差分定位方法 |
CN108802780A (zh) * | 2018-03-09 | 2018-11-13 | 东南大学 | 一种gps/bds差分系统间偏差特性分析方法 |
CN109031359A (zh) * | 2018-06-11 | 2018-12-18 | 武汉大学 | 一种基于微波暗室的gnss接收天线绝对相位中心校正方法 |
CN109581455A (zh) * | 2019-01-24 | 2019-04-05 | 东南大学 | 一种bds和gps融合的三频宽巷紧组合定位方法 |
-
2019
- 2019-12-26 CN CN201911387743.8A patent/CN110873887A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5451964A (en) * | 1994-07-29 | 1995-09-19 | Del Norte Technology, Inc. | Method and system for resolving double difference GPS carrier phase integer ambiguity utilizing decentralized Kalman filters |
EP0694791A1 (en) * | 1994-07-29 | 1996-01-31 | Del Norte Technology, Inc. | Method and system for resolving double difference GPS carrier phase integer ambiguity utilizing decentralized kalman filters |
CN106569242A (zh) * | 2016-11-04 | 2017-04-19 | 中国科学院测量与地球物理研究所 | 固定参考星的gnss单差处理方法 |
CN108802780A (zh) * | 2018-03-09 | 2018-11-13 | 东南大学 | 一种gps/bds差分系统间偏差特性分析方法 |
CN108519614A (zh) * | 2018-03-16 | 2018-09-11 | 东南大学 | 一种gps/bds紧组合载波差分定位方法 |
CN109031359A (zh) * | 2018-06-11 | 2018-12-18 | 武汉大学 | 一种基于微波暗室的gnss接收天线绝对相位中心校正方法 |
CN109581455A (zh) * | 2019-01-24 | 2019-04-05 | 东南大学 | 一种bds和gps融合的三频宽巷紧组合定位方法 |
Non-Patent Citations (2)
Title |
---|
苏春循: "基于估计差分系统间偏差的相对定位算法研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
金俭俭 等: "固定GEO卫星为参考的BDS与GPS紧组合定位方法", 《测绘科学》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115480274A (zh) * | 2021-05-31 | 2022-12-16 | 千寻位置网络有限公司 | 一种参考星选择方法、装置及接收机 |
CN116106953A (zh) * | 2023-04-14 | 2023-05-12 | 中南大学 | 一种gnss相对定位半参数估计的时间差分载波相位增强方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AU2018426707A1 (en) | Precise point position and real-time kinematic (PPP-RTK) positioning method and device | |
CN109196379B (zh) | 具有改进的模糊度解算的卫星导航接收器 | |
CN107728180B (zh) | 一种基于多维粒子滤波偏差估计的gnss精密定位方法 | |
CN111308528B (zh) | 一种北斗/gps紧组合虚拟参考站定位方法 | |
CN111913201B (zh) | 一种gnss差分定位方法、装置及计算机可读存储介质 | |
CN110873887A (zh) | Gps和bds紧结合的定位方法及可读存储介质 | |
WO2017160585A1 (en) | Satellite navigation receiver with improved ambiguity resolution | |
CN109597105B (zh) | 一种顾及载波系统间偏差的gps/glonass紧组合定位方法 | |
CN115201870B (zh) | 附先验约束的多频多模gnss非差非组合时间传递方法 | |
CN110824505B (zh) | Gnss卫星接收机的偏差估计方法及系统、定位方法及终端 | |
CN109683182A (zh) | 一种引入系统间偏差的gnss多模组合差分定位方法 | |
CN116148909B (zh) | 多频多模非组合精密单点定位瞬时窄巷模糊度固定方法 | |
CN113204042A (zh) | 一种基于精密单点定位的多星座联合列车定位方法 | |
CN111123322B (zh) | 星载gnss接收机的观测值实时数据预处理方法、系统、介质及设备 | |
CN111381253A (zh) | 用于估计全球导航卫星系统接收器的位置的系统及其方法 | |
CN111736183B (zh) | 一种联合bds2/bds3的精密单点定位方法和装置 | |
CN113933872A (zh) | 多系统差分定位方法及其系统 | |
CN110208841B (zh) | 一种改进的面向非重叠频率的gnss紧组合方法 | |
CN110954927B (zh) | 动态加权方法、装置及可读存储介质 | |
CN113791432A (zh) | 一种提高ais定位精度的组合导航定位方法 | |
CN117955554B (zh) | 基于预报拼接的低轨卫星实时钟差确定方法及系统 | |
CN108627858A (zh) | 一种基于码相位调整量的多径修复方法及装置 | |
CN114415213A (zh) | 顾及信号间偏差稳定性的多频多模gnss精密定位方法和设备 | |
JP7446729B2 (ja) | 測位方法及び測位システム | |
CN117452462B (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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20200310 |
|
WD01 | Invention patent application deemed withdrawn after publication |