CN115792979A - 一种基于pdop贡献度的卫星分步选星方法 - Google Patents
一种基于pdop贡献度的卫星分步选星方法 Download PDFInfo
- Publication number
- CN115792979A CN115792979A CN202211378144.1A CN202211378144A CN115792979A CN 115792979 A CN115792979 A CN 115792979A CN 202211378144 A CN202211378144 A CN 202211378144A CN 115792979 A CN115792979 A CN 115792979A
- Authority
- CN
- China
- Prior art keywords
- satellite
- satellites
- pdop
- contribution degree
- contribution
- 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
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种基于PDOP贡献度的卫星分步选星方法,包括以下步骤:获取待测卫星的位置数据,根据所述位置数据计算所述待测卫星的仰角和方位角;基于所述待测卫星的仰角和方位角,获得基础星座和备选卫星集,基于所述备选卫星集计算卫星的PDOP贡献度,选取卫星加入解算卫星集;设置选取卫星要求,对所述解算卫星集中的卫星进行判断,满足选星数量后,对所选卫星进行解算,获得解算后的位置速度信息。本发明对可见卫星进行选择,选取部分可见星参与定位解算,在保障精度的前提下,减小定位解算时间,降低由于解算延时带来的位置误差。
Description
技术领域
本发明属于全球导航卫星系统技术领域,尤其涉及一种基于PDOP贡献度的卫星分步选星方法。
背景技术
全球导航卫星系统(Global Navigation Satellite System,GNSS)是一个世界范围的定位和定时系统,包含一个或多个卫星星座,随着各国卫星导航系统星座的完善,太空中可供定位的卫星数量不断增加。过多的卫星会导致定位解算时间的增加,产生的延迟误差对于高速运动的飞机而言是不能忽视的。
现有的卫星接收机由于通道数的增加,多是使用所有可以观测到的卫星进行定位解算。但由于过多的卫星带来的延时问题,需要对可观测卫星进行选择。选择的要求是要兼顾定位精度和解算卫星数量两方面的内容。在保证定位精度不受过多影响的情况下,使用较少的卫星进行解算。这就意味着要选出对定位精度提升较高的卫星进行解算。
几何精度因子(Geometric Dilution Of Precision,GDOP)是用来描述GNSS定位精度的物理量,它直接影响到GNSS定位精度。其中位置精度因子(Position Dilution OfPrecision,PDOP)仅与参与解算的可见星位置有关,在伪距测量误差和钟差不变的情况下可见星几何构型直接影响到GNSS定位精度。因此对可见星进行选择时,以较好的PDOP作为选星目标。
发明内容
本发明的目的在于提出一种基于PDOP贡献度的卫星分步选星方法,对可见卫星进行选择,选取部分可见星参与定位解算,在保障精度的前提下,减小定位解算时间,降低由于解算延时带来的位置误差。
为实现上述目的,本发明提供了一种基于PDOP贡献度的卫星分步选星方法,包括以下步骤:
获取待测卫星的位置数据,根据所述位置数据计算所述待测卫星的仰角和方位角;
基于所述待测卫星的仰角和方位角,获得基础星座和备选卫星集,基于所述备选卫星集计算卫星的PDOP贡献度,选取卫星加入解算卫星集;
设置选取卫星要求,对所述解算卫星集中的卫星进行判断,满足选星数量后,对所选卫星进行解算,获得解算后的位置速度信息。
可选的,所述待测卫星的位置数据包括:
卫星位置坐标Xi(t),Yi(t),Zi(t),伪距观测值Si(t),接收机先验位置X,Y,Z。
可选的,基于所述待测卫星的仰角和方位角,获得基础星座和备选卫星集,具体包括:
比较所述待测卫星的仰角,获得基础星座;
将其余卫星中仰角大于35°的卫星纳入备选卫星集。
可选的,所述基础星座包括:天顶星和底角星,
其中,选取仰角最大的卫星作为所述基础星座中的天顶星;
在其余卫星中,根据仰角选取三颗分布最均匀的作为所述基础星座中的底角星。
可选的,基于所述备选卫星集计算卫星的PDOP贡献度,选取卫星加入解算卫星集,具体包括:
对所述备选卫星集解算卫星的PDOP贡献度,获取所述PDOP贡献度最大的卫星;
将所述卫星加入所述解算卫星集,获得新组合PDOP值。
可选的,对所述解算卫星集中的卫星进行判断,获得解算后的位置速度信息,判断方法具体包括:
对所述解算卫星集中的卫星进行位置迭代解算,满足所述选取卫星要求则对所述解算卫星集进行解算获得解算后的位置速度信息;反之则重新计算其余卫星的PDOP贡献度直到满足所述选取卫星要求。
本发明技术效果:本发明公开了一种基于PDOP贡献度的卫星分步选星方法,以PDOP为标准,选取对卫星组合PDOP改善更大的卫星参与解算,提高每一颗参与解算的卫星对定位精度的贡献;控制参与解算的卫星数量,在定位精度和解算时间中取得平衡可以降低卫星接收机芯片的算力要求,降低成本与功耗;本发明拥有较好的实时性和选星质量稳定性,能够有效减小由于计算延时带来的误差。
附图说明
构成本申请的一部分的附图用来提供对本申请的进一步理解,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。在附图中:
图1为本发明实施例基于PDOP贡献度的卫星分步选星方法的流程示意图;
图2为本发明实施例分布选星算法示意图。
具体实施方式
需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本申请。
需要说明的是,在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行,并且,虽然在流程图中示出了逻辑顺序,但是在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤。
如图1所示,本实施例中提供一种基于PDOP贡献度的卫星分步选星方法,包括以下步骤:
S1、获取当前接收机与各可见卫星的距离,各可见卫星当前位置等实验数据;
实验数据为参与解算的数据,每次参与解算的数据包括:
(a)卫星位置坐标Xi(t),Yi(t),Zi(t),其含义为t时刻时卫星i(i为卫星索引号)在地心地固系下的位置坐标。
(b)伪距观测值Si(t),为当前时刻t下,卫星i与接收机的观测距离,由接收机输出。
本发明下述实施例步骤中涉及的伪距观测值为经过电离层误差、对流层误差、卫星钟差修正等修正之后用于解算的伪距值。具体修正内容根据所要求的精度级别而不同。伪距修正步骤为卫星导航解算的必要步骤,具体需参照卫星导航定位算法。本发明适宜使用在伪距观测值修正步骤之后与伪距参与位置解算步骤之前。未经修正的伪距观测值也可使用,但可能将降低发明效果。
(c)接收机先验位置X0(t),Y0(t),Z0(t),代表t时刻接收机在地心地固系下的位置坐标。接收机先验位置可以采用当前时刻预测定位结果,包括但不限于GNSS系统当前时刻实测数据定位结果、其他传感器当前时刻定位结果,GNSS系统当前时刻预测结果等。
本实施例采用上一时刻的卫星定位系统解算的接收机位置为先验位置。
本实施例采用GPS31颗卫星与BeiDou44颗卫星的官方公布历书数据,进行卫星星座仿真。为了表现基于PDOP贡献度的卫星分步选星算法效果,选取民用客机在对GNSS精度要求较高的进近阶段航迹进行仿真,从而验证选星算法对定位精度的改善情况,航迹如图2所示。
S2、计算所有可见星的仰角和方位角;
S3、根据所有卫星的仰角与方位角,挑选四颗星作为基础星座卫星;
S31、比较所有卫星仰角,选择仰角最大的卫星为基础星座选择;
S32、在剩余可见星中,将仰角小于35°的卫星纳入备选卫星集中。在备选卫星集中,选取三颗卫星方位角间距最接近120°的组合。对方位角在0°至120°之间的卫星选取与其方位角差最接近±120°的两颗卫星作为底角星选择。计算三颗卫星间的方位角差值与120°的差,作为该底角星组合的δi。比较所有底角星组合的,选取δi最小的底角星组合作为基础星座的底角星组合。
S4、计算其余卫星的PDOP贡献度,选取贡献度高的卫星添加到解算卫星集中;
设接收机位置为X,Y,Z,对于第i颗卫星,卫星位置坐标为Xi(t),Yi(t),Zi(t),伪距观测值为Si(t)。那么观测矩阵H可写为:
令(HTH)-1=HHH,HHH主对角线上的元素分别为h11、h22、h33、h44,那么:
位置误差系数(Position Dilution Of Precision,PDOP)为:
当参与解算的卫星数量为n时,观测矩阵可改写为:
在增加一颗可见星后,观测矩阵H变为:
其中hn+1=[sinαn+1cosθn+1 cosαn+1cosθn+1 sinθn+1 1]。
根据Sherman-Morrison公式可得:
所以在增加可见星后对DOP值的影响可由Δhn+1对角线上元素大小得出。对角线上元素数值越大,主对角线上元素数值越小,DOP值越小。因此在备选卫星集中选择对DOP值改善最明显的卫星时,可以通过计算各卫星的Δhn+1矩阵得出该卫星的DOP贡献度,选择贡献度最大的卫星加入解算卫星集合。
S5、根据需要设置选星上限,解算卫星数量没有达到上限重复S4,达到上限后进入S6;
根据需要设置选星上限N,在当前选星数量不足N的情况下,回到S4继续选星,在数量达到N时,进入S6进行定位解算。
S6、对解算卫星集中的卫星进行定位迭代解算,输出解算后的位置速度信息;
对于第i颗卫星而言,观测方程如下:
Zi=HiX+εi
其中Zi为第i颗卫星的伪距测量值,Hi为第i颗卫星的观测值,εi为第i颗卫星的观测误差。
那么第i颗卫星观测方程为:
由于Si(t)为非线性函数,需将其泰勒级数展开可得:
其中Si0(t)为Si(t)的泰勒展开零阶常数项,X,Y,Z为当前接收机地心地固坐标系下接收机的位置,Xi(t),Yi(t),Zi(t)为卫星位置坐标,伪距观测值为Si(t)。第i颗卫星的观测向量位置估计量那么观测矩阵H可写为:
令p(t)=P(t)-S0(t)作为位置迭代计算的中间变量,则:
那么:
p(t)=[p1(t) p2(t) … pn(t)]T=HX+Δ
根据观测方程式位置估计量:
在位置迭代解算过程中,经过k次迭代计算后观测矩阵H(k)为:
经过k次位置迭代计算后接收机地心地固坐标系下三维位置X(k)应为:
最终的实验结果如下表1所示。
表1
解算时间(s) | 水平位置误差(m) | 下落阶段垂直位置误差(m) | |
本文方法 | 0.0196 | 12.701 | 7.259 |
全星解算 | 0.0416 | 16.945 | 7.752 |
提升度 | 52.88% | 25.02% | 6.41% |
本发明公开了一种基于PDOP贡献度的卫星分步选星方法,以PDOP为标准,选取对卫星组合PDOP改善更大的卫星参与解算,提高每一颗参与解算的卫星对定位精度的贡献。控制参与解算的卫星数量,在定位精度和解算时间中取得平衡可以降低卫星接收机芯片的算力要求,降低成本与功耗。与别的选星算法相比,拥有较好的实时性和选星质量稳定性,可以有效减小由于计算延时带来的误差。
以上所述,仅为本申请较佳的具体实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应该以权利要求的保护范围为准。
Claims (7)
1.一种基于PDOP贡献度的卫星分步选星方法,其特征在于,包括以下步骤:
获取待测卫星的位置数据,根据所述位置数据计算所述待测卫星的仰角和方位角;
基于所述待测卫星的仰角和方位角,获得基础星座和备选卫星集,基于所述备选卫星集计算卫星的PDOP贡献度,选取卫星加入解算卫星集;
设置选取卫星要求,对所述解算卫星集中的卫星进行判断,满足选星数量后,对所选卫星进行解算,获得解算后的位置速度信息。
2.如权利要求1所述的基于PDOP贡献度的卫星分步选星方法,其特征在于,所述待测卫星的位置数据包括:
卫星位置坐标Xi(t),Yi(t),Zi(t),伪距观测值Si(t),接收机先验位置X,Y,Z。
4.如权利要求1所述的基于PDOP贡献度的卫星分步选星方法,其特征在于,基于所述待测卫星的仰角和方位角,获得基础星座和备选卫星集,具体包括:
比较所述待测卫星的仰角,获得基础星座;
将其余卫星中仰角大于35°的卫星纳入备选卫星集。
5.如权利要求4所述的基于PDOP贡献度的卫星分步选星方法,其特征在于,所述基础星座包括:天顶星和底角星,
其中,选取仰角最大的卫星作为所述基础星座中的天顶星;
在其余卫星中,根据仰角选取三颗分布最均匀的作为所述基础星座中的底角星。
6.如权利要求4所述的基于PDOP贡献度的卫星分步选星方法,其特征在于,基于所述备选卫星集计算卫星的PDOP贡献度,选取卫星加入解算卫星集,具体包括:
对所述备选卫星集解算卫星的PDOP贡献度,获取所述PDOP贡献度最大的卫星;
将所述卫星加入所述解算卫星集,获得新组合PDOP值。
7.如权利要求1所述的基于PDOP贡献度的卫星分步选星方法,其特征在于,对所述解算卫星集中的卫星进行判断,获得解算后的位置速度信息,判断方法具体包括:
对所述解算卫星集中的卫星进行位置迭代解算,满足所述选取卫星要求则对所述解算卫星集进行解算获得解算后的位置速度信息;反之则重新计算其余卫星的PDOP贡献度直到满足所述选取卫星要求。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211378144.1A CN115792979A (zh) | 2022-11-04 | 2022-11-04 | 一种基于pdop贡献度的卫星分步选星方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211378144.1A CN115792979A (zh) | 2022-11-04 | 2022-11-04 | 一种基于pdop贡献度的卫星分步选星方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115792979A true CN115792979A (zh) | 2023-03-14 |
Family
ID=85435606
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211378144.1A Pending CN115792979A (zh) | 2022-11-04 | 2022-11-04 | 一种基于pdop贡献度的卫星分步选星方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115792979A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116299588A (zh) * | 2023-03-17 | 2023-06-23 | 哈尔滨工业大学(深圳) | 一种定位卫星选择方法、装置及存储介质 |
CN117741706A (zh) * | 2023-12-25 | 2024-03-22 | 北京航空航天大学 | 一种基于结果驱动机器学习的无人机诱骗选星方法 |
-
2022
- 2022-11-04 CN CN202211378144.1A patent/CN115792979A/zh active Pending
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116299588A (zh) * | 2023-03-17 | 2023-06-23 | 哈尔滨工业大学(深圳) | 一种定位卫星选择方法、装置及存储介质 |
CN116299588B (zh) * | 2023-03-17 | 2023-11-14 | 哈尔滨工业大学(深圳) | 一种定位卫星选择方法、装置及存储介质 |
CN117741706A (zh) * | 2023-12-25 | 2024-03-22 | 北京航空航天大学 | 一种基于结果驱动机器学习的无人机诱骗选星方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110045407B (zh) | 一种分布式伪卫星/gnss优化定位方法 | |
CN110007317B (zh) | 一种选星优化的高级接收机自主完好性监测方法 | |
CN115792979A (zh) | 一种基于pdop贡献度的卫星分步选星方法 | |
CN108196284B (zh) | 一种进行星间单差模糊度固定的gnss网数据处理方法 | |
CN108845340A (zh) | 基于gnss-rtk的定位方法 | |
Elsobeiey | Precise point positioning using triple-frequency GPS measurements | |
CN112394383B (zh) | 一种卫星与5g基站组合定位方法及装置 | |
CN110146908B (zh) | 一种虚拟参考站观测数据的生成方法 | |
CN109633718B (zh) | 一种归一化加权最小二乘导航定位方法 | |
CN108919316A (zh) | 一种基于局部球对称假设的单站多系统硬件延迟估计方法 | |
CN109212562A (zh) | 一种生成载波伪距观测量的方法 | |
CN106547002A (zh) | 一种提高rtk定位精度的选星方法 | |
MX2011005844A (es) | Un sistema y metodo para determinacion de posicion de referencia. | |
CN104502943B (zh) | 一种室内伪卫星差分相对定位方法 | |
CN105044733B (zh) | 一种高精度的导航卫星tgd参数标定方法 | |
CN113176596B (zh) | 气压高程约束定位方法 | |
CN105842710B (zh) | 一种基于vrs差分原理的低轨双星时差频差精密修正方法 | |
CN112198533B (zh) | 一种多假设下的地基增强系统完好性评估系统及方法 | |
CN111965674B (zh) | 一种基于自适应布谷鸟算法的北斗定位解算方法及系统 | |
CN113419262A (zh) | 一种全系统rtk快速选星方法 | |
CN106814376B (zh) | 一种快速精确厘米级单点定位方法 | |
CN115932920A (zh) | 对流层延迟的插值方法 | |
CN114325779A (zh) | 导航增强系统定位粗差检测方法及装置 | |
CN110045634B (zh) | 一种gnss基准站的非差误差建模方法 | |
CN113465575A (zh) | 一种基于对流层先验信息约束的高落差山区gnss高精度快速网解方法 |
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 |