CN102981167B - 一种gps/北斗***双模测时完好性监测方法 - Google Patents
一种gps/北斗***双模测时完好性监测方法 Download PDFInfo
- Publication number
- CN102981167B CN102981167B CN201210450203.1A CN201210450203A CN102981167B CN 102981167 B CN102981167 B CN 102981167B CN 201210450203 A CN201210450203 A CN 201210450203A CN 102981167 B CN102981167 B CN 102981167B
- Authority
- CN
- China
- Prior art keywords
- gps
- centerdot
- clock
- delta
- epsiv
- 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.)
- Expired - Fee Related
Links
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明所涉及的是一种测时完好性监测技术,具体涉及一种用户位置未知的情况下应用于测时服务GPS/北斗***双模测时完好性监测方法。本发明包括如下步骤:获取双模卫星导航接收机的GPS/北斗双模观测量;建立基于GPS***误差时和北斗***误差时的线性模型;统一标准差;获取钟差估计值;计算测时完好性风险;检测当前调节因子τ是否能够满足当前完好性要求;完成钟差估计;遍历调节因子。本发明利用全球卫星导航***GPS和北斗卫星导航***的伪距观测量,获得关于用户钟差的冗余估计值,由于用户钟漂与钟差之间的线性关系,使该方法能够直接检测用户端的异常钟漂,也就是间接监测用户端钟差,提高了监测性能。
Description
技术领域
本发明所涉及的是一种测时完好性监测技术,具体涉及一种用户位置未知的情况下应用于测时服务全球定位***GPS/北斗***双模测时完好性监测方法。
背景技术
目前很多工程应用依赖于卫星导航***(如美国的全球定位***GPS、我国的北斗卫星导航***等)作为它们的时间源。一般而言,用户端的时钟源一般由石英晶振体等组成,而温度、湿度、压力和振动等环境状况的变化都可使晶体振荡器的谐振频率随之发生改变,当谐振频率畸变到一定程度时,不仅测时精度变差,甚至会导致接收端所需输出原始观测量(特别是多普勒观测量)的畸变。应对这种信号畸变监测的解决思路之一是在用户端配置高精度的时钟源(如原子钟)作为参考,但是这种方法代价昂贵,并不适用于一般用户。解决思路之二是利用完好性监测方法,即要求用户位置精确已知,构建对应多颗卫星的伪距关于时钟误差的观测量,通过一致性检测技术对定位解在精度、完好性、可用性和连续性四个方面进行实时监测。但是实际应用中要求测时服务的用户其位置可能不能精确已知,因此其应用也具有较大的局限性。
由上述分析可知,在用户位置未知的情况下,将传统的完好性监测算法应用于测时服务,技术难点之一在于如何构建关于时钟的冗余观测量。此外,由于目前尚无针对测时的所需导航性能(Required Navigation Performance,RNP)要求,因而无法建立起衡量测时完好监测性能的指标。为此,本设计提出采用用户自定义指标的GPS/北斗双模测时完好性监测算法,双模星座不仅可以克服单星座无法提供钟差冗余观测量的天然缺陷,而且随着可用卫星数的增加,钟差的估计精度也可随之增高。用户自定义的完好性风险指标评估完好性监测性能,可以满足不同类型用户的测时等级需求,从而扩展本设计应用对象范围。
发明内容
本发明的目的在于提供一种克服钟差冗余观测量、提高精度的GPS/北斗***双模测时完好性监测方法。
本发明的目的是这样实现的:
一种GPS/北斗***双模测时完好性监测方法,包括如下步骤:
(1)获取双模卫星导航接收机的GPS/北斗双模观测量,其中GPS/北斗双模观测量的线性模型为:
该模型的等效表示为Y=GX+ε,下标GPS和BD代表分别对应GPS和北斗的变量,ρ表示伪距观测量,α、β、γ为对应卫星和用户的方向余弦矢量,xu、yu和zu为规定坐标系下的用户坐标,c为光速时间常量,δt代表用户钟差项,ε为观测量噪声,m为北斗卫星的个数和n为GPS卫星的个数;
(2)建立基于GPS***误差时和北斗***误差时的线性模型:
其中, 为***时钟误差观测量,k代表第k个历元,δfu为频漂,Ts为接收机观测量输出周期,σGPS和σBD分别为基于GPS***和北斗***的伪距观测量观测噪声标准差;
(3)统一标准差:完好性监测方法要求各观测量具有统一的统计特性,统一以σBD为标准差表示基于GPS***误差时和北斗***误差时的线性模型:
上述线性模型的等效形式为z=Hξ+ε′;
(4)获取钟差估计值: 其中次优加权矩阵w为 τ为调节因子以控制时钟估计引入的完好性风险,当且仅当τ=1时,w为最优加权矩阵,否则w为次优加权矩阵;
(5)计算测时完好性风险:
其中λ为由预先设置的漏检率所决定的非中心χ2分布中的非中心参数,τ为调节因子;
(6)检测采用当前调节因子τ的完好性风险是否能够满足当前完好性要求Preq,若满足要求,即满足Pt,ir<Preq,执行步骤(7),若不满足要求执行步骤(8);
(7)完成钟差估计,获得用户测时误差的最终解Δtu,final为:
(8)遍历调节因子,若存在次优加权矩阵,重新执行步骤(5);若遍历所有调节因子后仍不能满足完好性风险要求,向用户告警,结束该进程。
所述X的最优加权最小二乘解为:
本发明的有益效果在于:本发明利用全球卫星导航***GPS和北斗卫星导航***的伪距观测量,获得关于用户钟差的冗余估计值,由于用户钟漂与钟差之间的线性关系,使该方法能够直接检测用户端的异常钟漂,也就是间接监测用户端钟差,提高了监测性能。同时,利用次优加权矩阵损失一定的估计精度,提高了可用性水平。此外,针对不同类型用户对测时完好性不同等级需求,利用调节因子调节次优加权矩阵在估计精度和满足用户测时完好性需求方面达到最优平衡。
附图说明
图1是本发明的执行流程图。
具体实施方式
下面结合附图对本发明做进一步描述:
假设有m颗北斗卫星和n颗GPS卫星可用,则基于伪距观测量的线性模型可表示为:
其中下标GPS和BD代表分别对应GPS和北斗的变量,ρ表示伪距观测量,α、β、γ为对应卫星和用户的方向余弦矢量,xu、yu和zu为规定坐标系下的用户坐标,c为光速时间常量,δt代表用户钟差项,ε为观测量噪声。由式(1)可知,完成基于GPS/北斗双模定位模式下的用户位置解算,必须满足(m+n)>=5。为以下叙述方便,式(1)可等效表示为:
Y=GX+ε (2)
则可得关于X的最优加权最小二乘解可表示为:
从式(3)可以获取钟差的最优估计值,其中最优加权矩阵W取为
其中,σGPS和σBD分别为基于GPS和北斗的伪距观测量观测噪声标准差。和通过可以得到接收机关于GPS***时和北斗***时的估计值分别为和(以米为单位)。由于接收机时钟误差本质上是由于接收机本身的时钟源(如石英晶体、原子钟等)振荡以复制“重现”卫星发送的信号而产生,因此不同的***时误差都是基于同一个时钟频率基准,因此对这个公共的时间源基准进行监测将对实时性要求较高或基于授时服务需求来说,利用时钟完好性监测技术保证接收机钟差估计的稳定性成为有效的监测手段之一。
由于用户接收机钟差Δbu和频漂δfu之间可视为积分的关系,表示为
Δbu=Tsδfu (5)
其中Ts为接收机观测量输出周期。由式(5)可知频漂估计的精度和完好性性决定了钟差估计的完好性,而且对钟差估计的监测也可间接通过对频漂的监测得以实现。在多模接收机中,由于接收机时钟误差具有不同类型的***时钟误差观测量(如和),这种观测量的冗余性也就为实现频漂的完好性监测提供了理论基础。基于GPS***误差时和北斗***误差时在第k历元时刻的线性模型可表示为:
其中,
完好性监测算法要求各观测量具有统一的统计特性,因此可将式(6)统一以σBD为标准差可表示为:
为后续叙述方便,记并将式(7)等效描述为:
z=Hξ+ε′ (8)
为在估计精度和完好性监测水平之间取得最佳的平衡,采用以次优加权矩阵提升完好性的方法,加权钟漂估计值为:
次优加权矩阵w取为
其中τ为调节因子以控制时钟估计引入的完好性风险。当且仅当τ=1时,w为最优加权矩阵,否则w为次优加权矩阵。与水平/垂直保护水平的推导相类似,可以推得时钟保护水平TPL为
其中,λ为有预先设置的漏检率所决定的非中心χ2分布中的非中心参数。由于目前各类所需导航性能(Required Navigation Performance,RNP)标准中尚无关于时钟的完好性监测性能水平要求,因此在这里给出关于时钟完好性风险的计算公式,使得不同类型的用户可通过预先定义的相关完好性风险指标Pt,ir确定调节因子τ以满足可用性要求。时钟完好性风险可表示为:
其中ξ代表真实的频漂值。由前述假设可知估计误差值因此完好性风险可表示为:
其中为使所设计的TPL满足完好性风险要求,可调节τ通过损失一定的估计精度满足完好性监测性能的要求。由于难以得到关于τ的解析解,因此可以迭代搜索法获取τ的数值解,用户可以通过控制搜索步长以维持计算量和数值解精度之间的平衡。
当搜索得到满足用户自定义的完好性风险指标之后,用户测时误差的最终解Δtu,final取为:
当经过τ的遍历搜索之后,所得到的完好性风险仍然无法满足用户自定义的完好性监测性能指标,则算法向用户发出告警信号,提示用户此时***的测时精度置信度水平难以满足所需指标要求。
整个设计实现的流程如图1所示,其中GPS/北斗双模观测量由双模卫星导航接收机提供,钟差估计值可由式(9)获取,测时完好性风险可由式(13)计算获得。如果采用当前调节因子τ能够满足当前完好性要求,则可利用由调节因子构成的加权矩阵(如式(10))与式(9)完成钟差估计,否则需要改变调节因子直至满足完好性要求。当遍历所有调节因子后仍不能满足完好性风险要求,则需要向用户告警。
Claims (2)
1.一种GPS/北斗***双模测时完好性监测方法,其特征在于,包括如下步骤:
(1)获取双模卫星导航接收机的GPS/北斗双模观测量,其中GPS/北斗双模观测量的线性模型为:
该模型的等效表示为Y=GX+ε,下标GPS和BD代表分别对应GPS和北斗的变量,ρ表示伪距观测量,α、β、γ为对应卫星和用户的方向余弦矢量,xu、yu和zu为规定坐标系下的用户坐标,c为光速时间常量,δt代表用户钟差项,ε为观测量噪声,m为北斗卫星的个数和n为GPS卫星的个数;
(2)建立基于GPS***误差时和北斗***误差时的线性模型:
其中, 为***时钟误差观测量,k为第k个历元,δfu为频漂,Ts为接收机观测量输出周期,σGPS和σBD分别为基于GPS***和北斗***的伪距观测量观测噪声标准差;
(3)统一标准差:完好性监测方法要求各观测量具有统一的统计特性,统一以σBD为标准差表示基于GPS***误差时和北斗***误差时的线性模型:
上述线性模型的等效形式为z=Hξ+ε′
(4)获取钟差估计值:其中次优加权矩阵w为τ为调节因子以控制时钟估计引入的完好性风险,当且仅当τ=1时,w为最优加权矩阵,否则w为次优加权矩阵;
(5)计算测时完好性风险:
其中dv,λ为由预先设置的漏检率所决定的非中心χ2分布中的非中心参数,τ为调节因子;
(6)检测采用当前调节因子τ的完好性风险是否能够满足当前完好性要求Preq,若满足要求,即满足Pt,ir<Preq,执行步骤(7),若不满足要求执行步骤(8);
(7)完成钟差估计,获得最优加权最小二乘解用户测时误差的最终解Δtu,final为:
(8)遍历调节因子,若存在次优加权矩阵,重新执行步骤(5);若遍历所有调节因子后仍不能满足完好性风险要求,向用户告警,结束该进程。
2.根据权利要求1所述的一种GPS/北斗***双模测时完好性监测方法,其特征在于:所述X的最优加权最小二乘解为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210450203.1A CN102981167B (zh) | 2012-11-12 | 2012-11-12 | 一种gps/北斗***双模测时完好性监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210450203.1A CN102981167B (zh) | 2012-11-12 | 2012-11-12 | 一种gps/北斗***双模测时完好性监测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102981167A CN102981167A (zh) | 2013-03-20 |
CN102981167B true CN102981167B (zh) | 2014-08-20 |
Family
ID=47855371
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210450203.1A Expired - Fee Related CN102981167B (zh) | 2012-11-12 | 2012-11-12 | 一种gps/北斗***双模测时完好性监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102981167B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103454650B (zh) * | 2013-08-20 | 2015-06-24 | 北京航空航天大学 | 一种视觉辅助卫星完好性监测方法 |
CN104280756A (zh) * | 2014-10-30 | 2015-01-14 | 中国科学院国家天文台 | 基于接收机钟差广义延拓逼近法的卫星定位增强方法 |
CN107831648B (zh) * | 2017-09-06 | 2018-09-28 | 和芯星通科技(北京)有限公司 | 一种基于全球卫星导航***的高精度授时方法 |
CN108761498B (zh) * | 2018-03-13 | 2021-08-10 | 南京航空航天大学 | 一种针对高级接收机自主完好性监测的位置估计优化方法 |
CN108594275A (zh) * | 2018-04-26 | 2018-09-28 | 桂林电子科技大学 | 一种北斗+gps双模单点定位方法 |
CN108931790B (zh) * | 2018-05-17 | 2021-06-11 | 西安空间无线电技术研究所 | 一种导航卫星钟慢变故障的卫星自主完好性监测方法 |
CN109581445B (zh) * | 2018-11-01 | 2021-03-19 | 北京航空航天大学 | 一种基于北斗星座的araim子集选择方法及*** |
CN111142124B (zh) * | 2018-11-02 | 2023-04-28 | 千寻位置网络有限公司 | 全球卫星导航***状态空间表达方式完好性监测方法及装置 |
CN111413719B (zh) * | 2020-03-21 | 2022-07-15 | 哈尔滨工程大学 | 一种基于神经网络的北斗实时精密时钟预测方法 |
CN111983656B (zh) * | 2020-08-17 | 2023-10-27 | 电信科学技术第五研究所有限公司 | 一种gnss时统设备中对gnss时间进行最优估计的方法、装置及计算机介质 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7286082B2 (en) * | 2004-02-26 | 2007-10-23 | Raytheon Company | Method and apparatus for wide area augmentation system having GEO uplink subsystem with enhanced clock steering |
CN101799552B (zh) * | 2010-03-11 | 2012-11-21 | 北京航空航天大学 | 双***组合卫星导航接收机定位方法 |
CN102176031B (zh) * | 2011-01-06 | 2013-02-20 | 中国科学院国家授时中心 | 双模导航***中基于***时差接收机完好性故障检测方法 |
-
2012
- 2012-11-12 CN CN201210450203.1A patent/CN102981167B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN102981167A (zh) | 2013-03-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102981167B (zh) | 一种gps/北斗***双模测时完好性监测方法 | |
CN107193029B (zh) | 北斗三频信号的网络rtk基准站间模糊度快速确定方法 | |
Jakowski et al. | Total electron content models and their use in ionosphere monitoring | |
CN109459778B (zh) | 基于抗差方差分量估计的码伪距/多普勒联合测速方法及其应用 | |
CN101971047B (zh) | 卫星导航***的实时完整性监测设备及方法 | |
CN109709591A (zh) | 一种面向智能终端的gnss高精度定位方法 | |
CN103453906B (zh) | 卫星轨道的预测方法 | |
Viandier et al. | GNSS performance enhancement in urban environment based on pseudo-range error model | |
CN101680943A (zh) | 实时动态(rtk)定位中距离相关的误差减轻 | |
CN106324629A (zh) | 一种bds_gps_glonass融合精密单点定位方法 | |
CN102565821B (zh) | 多普勒频偏辅助的卫星导航信号载波周跳实时探测与修复的方法 | |
CN113848577A (zh) | 一种基于动态分区的大规模gnss网并行解算方法及*** | |
CN104597465A (zh) | 一种提高gps与glonass组合精密单点定位收敛速度的方法 | |
CN106842236B (zh) | Gnss接收机周跳探测与修复处理装置 | |
CN113253303A (zh) | 一种用于实时监测单频星基增强***性能的方法 | |
CN105044738A (zh) | 一种接收机自主完好性监视的预测方法及预测*** | |
CN111596315A (zh) | 一种用于实时监测双频多星座星基增强***性能的方法 | |
JP2013534623A (ja) | グローバル・ナビゲーション衛星システム−gnss及び/又は擬似衛星に基づく構造物の地震動又は振動を測定するシステム | |
CN102486540B (zh) | 一种应用于全球卫星定位与导航***中的快速定位方法 | |
CN104050389A (zh) | 一种实时在线评估导航***精确度和完好性的方法 | |
CN108181632B (zh) | 基于模糊度全微分的gnss单频数据周跳探测与修复方法 | |
CN103576167A (zh) | 基于hht和支持向量机的周跳探测与修复方法 | |
CN105509770A (zh) | 一种gnss和mems组合导航***中的气压计在线校正方法 | |
Seepersad et al. | Integrity monitoring in precise point positioning | |
CN106842242B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140820 Termination date: 20191112 |