CN103823224A - 一种基于北斗卫星导航***的分轨选星方法 - Google Patents
一种基于北斗卫星导航***的分轨选星方法 Download PDFInfo
- Publication number
- CN103823224A CN103823224A CN201410077781.4A CN201410077781A CN103823224A CN 103823224 A CN103823224 A CN 103823224A CN 201410077781 A CN201410077781 A CN 201410077781A CN 103823224 A CN103823224 A CN 103823224A
- Authority
- CN
- China
- Prior art keywords
- satellite
- delta
- star
- satellites
- footstock
- 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
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/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/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/28—Satellite selection
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
本发明公开了一种基于北斗卫星导航***的分轨选星方法。该方法首先通过选取满足高度角预设阈值的卫星作为待选卫星,其次根据待选卫星的轨道类型来划分高轨卫星组和中轨卫星组,在高轨卫星组中确定待选顶座星,在中轨卫星组中选取待选底座星;最后通过计算一颗待选顶座星与三颗待选底座星所构成的定位星座四面体体积来确定当前定位星。本发明针对北斗导航***的空间星座布局角度来设计,在不影响解算精度的前提下有效降低计算复杂度,提高导航解算的实时性。
Description
技术领域
本发明属于卫星导航技术领域,具体地说是一种基于北斗卫星导航***的分轨选星方法。
背景技术
我国正快速发展建设具有自主知识产权的北斗2代全球卫星导航***,根据北斗空间信号接口控制文件(ICDv1.0)的介绍,北斗***全星座建成之后,其空间段由5颗地球静止轨道(GEO)卫星、27颗中圆地球轨道(MEO)卫星和3颗倾斜地球同步轨道(IGSO)卫星组成。目前在轨运行的北斗卫星已经达到16颗,初步形成覆盖亚太地区的定位能力,随着北斗星座的逐步完善,同时观测到的卫星数目也在不断增多。
将所有可见星用于导航解算虽然可以一定程度的减小几何精度因子(GDOP)值,但是会给接收机硬件跟踪通道以及导航解算模块带来负担。因此对当前可见星进行合理优选,在不影响精度的前提下降低计算复杂度,具有非常现实的意义和价值。
传统的选星算法主要包括:最小几何精度因子(GDOP)选星法、最大矢端四面体体积法、最大行列式法等,这些算法均需要对所有可见星进行遍历运算以寻求最佳的选星方案,但在可见星数量较多时计算量巨大,影响算法的实时性。虽然目前已出现了一些以牺牲很小的GDOP值来换取解算实时性的次优选星算法,但迄今为止,尚未见到针对北斗异质星座特性开展的选星方法的报道。
发明内容
本发明解决的技术问题是提供一种针对北斗异质星座特性且无需遍历所有可见星即可获取当前定位星的分轨选星方法。
为解决上述技术问题,本发明一种基于北斗卫星导航***的分轨选星方法,其特征在于,包括以下步骤:
步骤1、选取当前所有可见卫星中高度角大于等于预设阈值的可见卫星作为待选卫星;
步骤2、根据待选卫星所在轨道的类型将其分为高轨卫星组和中轨卫星组;
步骤3、在高轨卫星组中确定待选顶座星;
步骤4、选择中轨卫星组中任意三颗卫星作为待选底座星,计算由上述三颗待选底座星与待选顶座星所构成的定位星座四面体体积;
步骤5、重复步骤4直至遍历选择完中轨卫星组中的所有三颗卫星的组合,将其中四面体体积最大的定位星座所对应的一个待选顶座星和三个待选底座星作为当前定位星。
进一步地优选方案,本发明基于北斗卫星导航***的分轨选星方法中,所述步骤1中可见卫星的高度角按照以下步骤获取,具体为:
步骤1-1、由卫星星历数据计算当前所有可见卫星在地心地固坐标系下的位置坐标,SatPosCar=[xn,yn,zn],其中n=1,2,3…,n表示卫星编号;
步骤1-2、将卫星位置由地心地固坐标系转换到站心坐标系SatPosGeo;
SatPosGeo=S·SatPosCar, λ,φ分别为接收机对应的经度和纬度;
步骤1-3、根据几何关系计算各个可见卫星的高度角θ;
式中, UsrPosCar=[x,y,z]为接收机在地心地固坐标系中的坐标。
进一步地优选方案,本发明基于北斗卫星导航***的分轨选星方法中,步骤1中高度角阈值为5°。
进一步地优选方案,本发明基于北斗卫星导航***的分轨选星方法中,步骤2中将位于地球静止轨道与倾斜地球同步轨道的待选卫星归为高轨卫星组,将位于中圆地球轨道的待选卫星归为中轨卫星组。
进一步地优选方案,本发明基于北斗卫星导航***的分轨选星方法中,所述步骤3中在高轨卫星组中确定待选顶座星,具体为:
判断当前高轨卫星组中的卫星数量,若只存在1颗卫星则将该卫星作为待选顶座星,若至少存在2颗卫星则分别选取高度角最大和最小的两颗卫星作为待选顶座星。
进一步地优选方案,本发明基于北斗卫星导航***的分轨选星方法中,步骤4中计算定位星座四面体体积时,当待选顶座星为两颗时,分别计算这两颗待选顶座星与三个待选底座星所构成的定位星座四面体的体积;
本发明与现有技术相比具备以下显著的优点:1)本发明的分轨选星方法针对北斗导航***的空间星座布局角度来设计;2)本发明在不影响解算精度的前提下有效降低计算复杂度,提高导航解算的实时性,对后续开展的多星座协同导航具有很强的适应性。
附图说明
图1为本发明分轨选星方法的流程图;
图2为本发明分轨选星方法与最小GDOP选星法的仿真结果对比曲线。
具体实施方式
如图1所示,本发明一种基于北斗卫星导航***的分轨选星方法,包括以下步骤:
步骤1、选取当前所有可见卫星中高度角大于等于预设阈值的可见卫星作为待选卫星,其中高度角预设阈值为5°,可见卫星的高度角按照以下步骤获得,具体为:
步骤1-1、由卫星星历数据计算当前所有可见卫星在地心地固坐标系下的位置坐标,SatPosCar=[xn,yn,zn],其中n=1,2,3…,n表示卫星编号;
步骤1-2、将卫星位置由地心地固坐标系转换到站心坐标系SatPosGeo,由于地心地固坐标系下的位置表达不能与实际物理意义相联系,不便于计算卫星的高度角与方位角,因此需要转换到站心坐标系即东北天坐标系中;
SatPosGeo=S·SatPosCar, λ,φ分别为接收机对应的经度和纬度;
步骤1-3、根据几何关系计算各个可见卫星的高度角θ;
式中, UsrPosCar=[x,y,z]为接收机在地心地固坐标系中的坐标;
对于接收机来说,对接收到的卫星信号进行解码,获得当前卫星位置是进行解算的重要前提,高度角计算过程为本领域的常规方法,具体参见《GPS原理与接收机设计》谢钢3.1.5节的相关内容,此处不再赘述!
本发明中选取高度角大于等于5°的卫星是为了降低多路径和大气折射等信号传播路径上的误差;
步骤2、根据待选卫星所在轨道的类型将其分为高轨卫星组和中轨卫星组:
由于地球静止轨道(GEO)与倾斜地球同步轨道(IGSO)卫星的轨道高度均为35786千米,将其归类为高轨卫星组;
由于中圆地球轨道(MEO)卫星的轨道高度为21528千米,将其归类为中轨卫星组;
步骤3、在高轨卫星组中确定待选顶座星:
判断当前高轨卫星组中的卫星数量,若只存在1颗卫星则将该卫星作为待选顶座星,若至少存在2颗卫星则分别选取高度角最大和最小的两颗卫星作为待选顶座星;
步骤4、选择中轨卫星组中任意三颗卫星作为待选底座星,计算由上述三颗待选底座星与待选顶座星所构成的定位星座四面体体积,计算定位星座四面体体积时,当待选顶座星为两颗时,分别计算这两颗待选顶座星与三个待选底座星所构成的定位星座四面体的体积;
其中,定位星座四面体体积计算公式如下:式中,A为待选顶座星在地心地固坐标系下的坐标、B、C、D分别为三颗底座星在地心地固坐标系下的坐标;
步骤5、重复步骤4直至遍历选择完中轨卫星组中的所有三颗卫星的组合,将其中四面体体积最大的定位星座所对应的一个待选顶座星和三个待选底座星作为当前定位星。
如图1所示,基于目前北斗卫星的分布情况,对于步骤3中高轨卫星组中卫星的数目一般均大于等于1,步骤4中的中轨卫星组中卫星数目一般均大于等于3;若出现高轨卫星组中无卫星或中轨卫星组中卫星数目小于3时,则接收机重新搜索卫星。
如图2所示,采用仿真的方式,对使用本发明方法后的北斗卫星导航***的实时GDOP值与传统最小GDOP选星方法进行对比验证。
首先利用卫星仿真工具包软件建立了北斗全面建成后的全星座模型,获取一段时间的精确卫星位置信息,仿真时间为24小时,采样周期为1分钟,接收机位置设为南京的经度纬度高度坐标[118.8°32.03°10m]。
通过编程实现了本发明中的选星算法以及最小GDOP选星算法,两者的实时GDOP对比曲线如图2所示,表1为对图1数据的统计分析结果,由图2和表1可以看出本发明的选星方法所计算出的GDOP值均在3.9以内,达到优良等级,而且与最小GDOP选星法的结果相比,73.1%的数据差值小于0.2,99%的数据差值小于0.5,说明本方法可以正确地选择定位星。
表1:
GDOP值 | 最小GDOP选星法 | 本文算法 |
最大值 | 3.6080 | 3.8366 |
最小值 | 2.1612 | 2.1612 |
均值 | 2.6387 | 2.7687 |
方差 | 0.1009 | 0.0937 |
Claims (6)
1.一种基于北斗卫星导航***的分轨选星方法,其特征在于,包括以下步骤:步骤1、选取当前所有可见卫星中高度角大于等于预设阈值的可见卫星作为待选卫星;
步骤2、根据待选卫星所在轨道的类型将其分为高轨卫星组和中轨卫星组;
步骤3、在高轨卫星组中确定待选顶座星;
步骤4、选择中轨卫星组中任意三颗卫星作为待选底座星,计算由上述三颗待选底座星与待选顶座星所构成的定位星座四面体体积;
步骤5、重复步骤4直至遍历选择完中轨卫星组中的所有三颗卫星的组合,将其中四面体体积最大的定位星座所对应的一个待选顶座星和三个待选底座星作为当前定位星。
2.根据权利要求1所述的基于北斗卫星导航***的分轨选星方法,其特征在于,
步骤1中可见卫星的高度角按照以下步骤获得,具体为:
步骤1-1、由卫星星历数据计算当前所有可见卫星在地心地固坐标系下的位置坐标,SatPosCar=[xn,yn,zn],其中n=1,2,3…,n表示卫星编号;
步骤1-2、将卫星位置由地心地固坐标系转换到站心坐标系SatPosGeo;
SatPosGeo=S·SatPosCar, λ,φ分别为接收机对应的经度和纬度;
步骤1-3、根据几何关系计算各个可见卫星的高度角θ;
式中, UsrPosCar=[x,y,z]为接收机在地心地固坐标系中的坐标。
3.根据权利要求1所述的基于北斗卫星导航***的分轨选星方法,其特征在于,步骤1中高度角预设阈值为5°。
4.根据权利要求1或2所述的基于北斗卫星导航***的分轨选星方法,其特征在于,步骤2中将位于地球静止轨道与倾斜地球同步轨道的待选卫星归为高轨卫星组,将位于中圆地球轨道的待选卫星归为中轨卫星组。
5.根据权利要求4所述的基于北斗卫星导航***的分轨选星方法,其特征在于,步骤3中在高轨卫星组中确定待选顶座星,具体为:
判断当前高轨卫星组中的卫星数量,若只存在1颗卫星则将该卫星作为待选顶座星,若至少存在2颗卫星则分别选取高度角最大和最小的两颗卫星作为待选顶座星。
6.根据权利要求5所述的基于北斗卫星导航***的分轨选星方法,其特征在于,步骤4中计算定位星座四面体体积时,当待选顶座星为两颗时,分别计算这两颗待选顶座星与三个待选底座星所构成的定位星座四面体的体积;
其中,定位星座四面体体积计算公式如下:式中,A为待选顶座星在地心地固坐标系下的坐标、B、C、D分别为三颗底座星在地心地固坐标系下的坐标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410077781.4A CN103823224B (zh) | 2014-03-04 | 2014-03-04 | 一种基于北斗卫星导航***的分轨选星方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410077781.4A CN103823224B (zh) | 2014-03-04 | 2014-03-04 | 一种基于北斗卫星导航***的分轨选星方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103823224A true CN103823224A (zh) | 2014-05-28 |
CN103823224B CN103823224B (zh) | 2016-04-06 |
Family
ID=50758373
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410077781.4A Active CN103823224B (zh) | 2014-03-04 | 2014-03-04 | 一种基于北斗卫星导航***的分轨选星方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103823224B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105005056A (zh) * | 2015-07-24 | 2015-10-28 | 上海微小卫星工程中心 | 一种基于动力学轨道外推的星载辅助gps方法及*** |
CN107064967A (zh) * | 2017-05-23 | 2017-08-18 | 南京航空航天大学 | 一种多星座接收机冷启动可见星搜索方法 |
CN107167824A (zh) * | 2017-07-26 | 2017-09-15 | 天津博创金成技术开发有限公司 | 一种北斗卫星导航***快速选星方法 |
CN107390238A (zh) * | 2017-07-23 | 2017-11-24 | 天津博创金成技术开发有限公司 | 一种北斗导航星座快速选星方法 |
CN107991696A (zh) * | 2017-11-16 | 2018-05-04 | 北京卫星信息工程研究所 | 针对终端载体为高轨卫星的卫星导航信号仿真方法 |
CN113419262A (zh) * | 2021-05-25 | 2021-09-21 | 武汉导航与位置服务工业技术研究院有限责任公司 | 一种全***rtk快速选星方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101571581A (zh) * | 2009-06-15 | 2009-11-04 | 北京航空航天大学 | 全球导航卫星***快速选星方法 |
CN101866010A (zh) * | 2010-05-28 | 2010-10-20 | 北京航空航天大学 | 一种卫星导航***选星方法 |
CN102520422A (zh) * | 2011-12-22 | 2012-06-27 | 成都金本华科技有限公司 | 北斗一号卫星无源定位选星方法 |
CN103185887A (zh) * | 2011-12-27 | 2013-07-03 | 东莞市泰斗微电子科技有限公司 | 一种卫星导航***的选星方法 |
CN103499822A (zh) * | 2013-10-14 | 2014-01-08 | 北京大学 | 一种基于最优gdop和牛顿恒等式的双星座快速选星方法 |
-
2014
- 2014-03-04 CN CN201410077781.4A patent/CN103823224B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101571581A (zh) * | 2009-06-15 | 2009-11-04 | 北京航空航天大学 | 全球导航卫星***快速选星方法 |
CN101866010A (zh) * | 2010-05-28 | 2010-10-20 | 北京航空航天大学 | 一种卫星导航***选星方法 |
CN102520422A (zh) * | 2011-12-22 | 2012-06-27 | 成都金本华科技有限公司 | 北斗一号卫星无源定位选星方法 |
CN103185887A (zh) * | 2011-12-27 | 2013-07-03 | 东莞市泰斗微电子科技有限公司 | 一种卫星导航***的选星方法 |
CN103499822A (zh) * | 2013-10-14 | 2014-01-08 | 北京大学 | 一种基于最优gdop和牛顿恒等式的双星座快速选星方法 |
Non-Patent Citations (2)
Title |
---|
初海彬 等: "基于中轨与静止轨道卫星定位***", 《中国空间科学技术》, no. 2, 30 April 2005 (2005-04-30), pages 52 - 58 * |
张昀申: "北斗导航***多轨道卫星星座分析与设计", 《舰船电子工程》, vol. 33, no. 4, 31 December 2013 (2013-12-31) * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105005056A (zh) * | 2015-07-24 | 2015-10-28 | 上海微小卫星工程中心 | 一种基于动力学轨道外推的星载辅助gps方法及*** |
CN107064967A (zh) * | 2017-05-23 | 2017-08-18 | 南京航空航天大学 | 一种多星座接收机冷启动可见星搜索方法 |
CN107064967B (zh) * | 2017-05-23 | 2020-03-24 | 南京航空航天大学 | 一种多星座接收机冷启动可见星搜索方法 |
CN107390238A (zh) * | 2017-07-23 | 2017-11-24 | 天津博创金成技术开发有限公司 | 一种北斗导航星座快速选星方法 |
CN107167824A (zh) * | 2017-07-26 | 2017-09-15 | 天津博创金成技术开发有限公司 | 一种北斗卫星导航***快速选星方法 |
CN107991696A (zh) * | 2017-11-16 | 2018-05-04 | 北京卫星信息工程研究所 | 针对终端载体为高轨卫星的卫星导航信号仿真方法 |
CN107991696B (zh) * | 2017-11-16 | 2021-08-13 | 北京卫星信息工程研究所 | 针对终端载体为高轨卫星的卫星导航信号仿真方法 |
CN113419262A (zh) * | 2021-05-25 | 2021-09-21 | 武汉导航与位置服务工业技术研究院有限责任公司 | 一种全***rtk快速选星方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103823224B (zh) | 2016-04-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103823224B (zh) | 一种基于北斗卫星导航***的分轨选星方法 | |
CN103760572B (zh) | 一种基于区域cors的单频ppp电离层加权方法 | |
CN102998690B (zh) | 一种基于gps载波双差方程的姿态角直接求解方法 | |
CN1996040B (zh) | 一种用于双星座卫星定位***的选星方法 | |
CN106338753A (zh) | 一种基于地面站/星间链路/gnss联合测量的地球同步轨道星座定轨方法 | |
CN101038337B (zh) | 储存全球导航卫星***的辅助校正数据的方法、装置及全球导航卫星***接收器 | |
Nadarajah et al. | Instantaneous GPS–Galileo attitude determination: single-frequency performance in satellite-deprived environments | |
CN105629263A (zh) | 一种对流层大气延迟误差估计改正方法和改正*** | |
CN108226978B (zh) | 一种基于wgs-84模型的双星定位方法 | |
CN101833080A (zh) | 一种利用gps***附加约束条件的载体姿态测量方法 | |
CN101887128A (zh) | 确定全球卫星导航***导航卫星频间偏差的方法 | |
CN104101888A (zh) | 一种适合于中国区域的gnss电离层延迟精确建模方法 | |
Ma et al. | Estimation of GPS instrumental biases from small scale network | |
CN104133231B (zh) | 一种基于积分多普勒平滑伪距的导航定位方法 | |
CN101561496A (zh) | 一种伪卫星和惯性组合导航***的非线性补偿方法 | |
CN107607971A (zh) | 基于gnss共视时间比对算法的时间频率传递方法及接收机 | |
Xue et al. | Understanding GDOP minimization in GNSS positioning: Infinite solutions, finite solutions and no solution | |
CN105044738A (zh) | 一种接收机自主完好性监视的预测方法及预测*** | |
Soycan et al. | Precise point positioning versus traditional solution for GNSS networks | |
Du et al. | The realization and evaluation of PPP ambiguity resolution with INS aiding in marine survey | |
CN114384564B (zh) | 一种基于多源数据驱动的电离层层析成像方法 | |
CN107782317A (zh) | 一种倒vlbi深空单向导航定位方法 | |
Zheng et al. | Analysis of a HAPS-aided GNSS in urban areas using a RAIM algorithm | |
CN105510946A (zh) | 一种bds卫星载波相位整周模糊度快速解算方法 | |
CN106814376B (zh) | 一种快速精确厘米级单点定位方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20140528 Assignee: Jiangsu the Big Dipper is open to the traffic Electronics Co., Ltd. Assignor: Nanjing University of Aeronautics and Astronautics Contract record no.: 2018320000067 Denomination of invention: Orbit-separating satellite selection method based on Beidou satellite navigation system Granted publication date: 20160406 License type: Exclusive License Record date: 20180320 |
|
EE01 | Entry into force of recordation of patent licensing contract |