[发明专利]一种多星座组合导航系统的故障识别方法有效
申请号: | 201010617250.1 | 申请日: | 2010-12-31 |
公开(公告)号: | CN102135621A | 公开(公告)日: | 2011-07-27 |
发明(设计)人: | 陆伟宁;王千喜;翟羽佳;杨晓昆;胡强;李秋凤;刘岩 | 申请(专利权)人: | 中国航天科工信息技术研究院 |
主分类号: | G01S19/20 | 分类号: | G01S19/20 |
代理公司: | 北京中海智圣知识产权代理有限公司 11282 | 代理人: | 齐晓静 |
地址: | 100048 北京市*** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 本发明公开了一种多星座组合导航系统的故障识别方法,包括进行多星座系统时和空间坐标的统一、故障检测可用性判定、故障判定、识别并排除故障卫星等步骤,能够对多颗卫星同时发生故障进行监测和识别,通过对多星座系统时和空间坐标的统一及多次检测和判断,大大提高了识别和排除故障卫星的成功率,增强了多星座组合导航系统导航定位精度、定位性能和可靠性,因而提高了导航服务性能。 | ||
搜索关键词: | 一种 星座 组合 导航系统 故障 识别 方法 | ||
【主权项】:
1.一种多星座组合导航系统的故障识别方法,其特征在于,包括以下步骤:1)进行多星座系统时和空间坐标的统一1.1)进行多星座系统时的统一:根据每一星座的系统时与UTC之间存在的固定转换关系,将不同星座的系统时进行粗略的解算,统一所有星座系统的系统时间;各星座的系统时与UTC之间的钟差由下述的伪距观测方程得到:ρ i g = [ ( x - X i g ) 2 + ( y - Y i g ) 2 + ( z - Z i g ) 2 ] 1 / 2 + cδ T r g ]]>δ T r g = [ [ ( x - X i g ) 2 + ( y - Y i g ) 2 + ( z - Z i g ) 2 ] 1 / 2 - ρ i g ] / c ]]> 其中,上标g表示系统序号,i为观测值序号(i=1,2,3,4,5);(x,y,z)为接收机在选定坐标系下的坐标,(Xg,Yg,Zg)为利用坐标转换公式将各个星座中的观测卫星转换到选定坐标系下的坐标;为接收机与各星座间的钟差(r=1,2...N),是每颗可视卫星的伪距观测值,c是光速;对于N个星座,接收机有N+3个需要测量的未知量,需要N+3个伪距观测方程来求解,得到各星座的系统时与UTC之间的钟差1.2)进行多星座空间坐标的统一空间坐标的统一采用如下坐标转换公式完成:x y z sys 1 = Δx Δy Δz + ( 1 + m ) 1 - θ z θ y θ z 1 - θ x - θ y θ x 1 x y z sys ]]> 其中,(Δx,Δy,Δz)为地球中心偏移量,θx,θy,θz坐标轴旋转角,m为比例因子,(x,y,z)sys为待转换坐标系中的坐标,(x,y,z)sys1为目标坐标系的坐标;2)故障检测可用性判定2.1)判断可视卫星个数,如果可视卫星个数<N+3,说明故障检测无法进行,系统进行完好性报警,否则继续;2.2)可用性判定:首先求出故障检测门限值σT,公式如下:σ T = σ 0 × T / n - 4 ]]> 其中,σ0为伪距测量误差的方差值;n是可视卫星的个数;门限值T由下式确定:Pr ( SSE / σ 0 2 < T 2 ) = ∫ 0 T 2 f χ ( n - 4 , λ ) 2 ( x ) dx = 1 - P FA ]]> 其中,PFA为可容忍的虚警概率;为自由度为n-4的χ2分布的概率密度函数;为残差统计平方和小于门限值T的统计概率;式中:SSE = v T Wv / σ 0 2 ]]> v=(I-G(GTWG)-1GTW)ε其中,G是由各卫星至接收机的方向余弦向量构成的线性化矩阵,ε是伪距误差向量,W为n×n维观测伪距权矩阵,I为单位阵;假设第i颗卫星存在故障,其偏差为bi,SSE则服从自由度为n-4的非中心化χ2分布,非中心参数λ可以由下式得到:λ = E ( v T Wv ) / σ 0 2 = RPE i 2 / σ 0 2 δHDOP i 2 ]]> 其中,RPE i = ( A 1 i 2 + A 2 i 2 ) W ii 2 b i ]]>δ HDOP i = HDOP i - HDOP = A 1 i 2 + A 2 i 2 Q v ii ]]> A=(GTWG)-1GTQv=W-1-G(GTWG)-1GTHDOP表示所有观测卫星的水平定位精度因子,HDOPi表示去掉第i颗卫星后的水平定位精度因子;计算δHDOPmax;故障检验前,实时就算出各卫星对应的HDOPi,去并取其中的最大值为δHDOPmax;计算水平定位误差保护限值HPL,公式如下;HPL = δHDOP max × σ 0 × λ ]]> 将HPL与水平误差保护限值HAL比较,如果超限,则故障检测不可用,系统进行完好性报警,否则继续;3)故障判定将实际观测伪距的误差的方差与步骤2)中求得的故障检测门限值σT比较,如果则表示检测到故障,工作继续,否则工作结束;4)识别并排除故障卫星4.1)判断可视卫星个数,如果可视卫星个数<N+4,故障识别无法进行,系统进行完好性报警,否则继续;4.2)识别故障卫星:故障卫星的识别采用特征线偏差线法,首先将观测系数矩阵G进行QR分解,得到矩阵QT:Q T = Q X Q P ]]> 其中,QX为QT前4行,QP为QT剩下的n-4行;奇偶残差矢量p为 p=QPy=QP(Gx+ε)=QPε计算每颗卫星的特征偏差线Kcha:K cha = Q p ( 1 , i ) Q p ( 2 , i ) ]]> 计算奇偶空间矢量的特征偏差斜率Kp:Kp=p1/p2其中,p1与p2是奇偶矢量p的元素,如果第i颗卫星的Kcha与Kp非常接近,则第i颗卫星被识别为故障卫星;4.3)排除故障卫星:将分步骤4.2)中识别的故障卫星排除;5)重复步骤1)至步骤4);6)重复步骤2)和步骤3),如果没有故障,则证明已经成功排除故障卫星,工作结束,否则认为此次测量数据无法完成自主完好性测试,系统进行完好性报警。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国航天科工信息技术研究院,未经中国航天科工信息技术研究院许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201010617250.1/,转载请声明来源钻瓜专利网。
- 同类专利
- 专利分类
G01 测量;测试
G01S 无线电定向;无线电导航;采用无线电波测距或测速;采用无线电波的反射或再辐射的定位或存在检测;采用其他波的类似装置
G01S19-00 卫星无线电信标定位系统;利用这种系统传输的信号确定位置、速度或姿态
G01S19-01 .传输时间戳信息的卫星无线电信标定位系统,例如,GPS [全球定位系统]、GLONASS[全球导航卫星系统]或GALILEO
G01S19-38 .利用卫星无线电信标定位系统传输的信号来确定导航方案
G01S19-39 ..传输带有时间戳信息的卫星无线电信标定位系统,例如GPS [全球定位系统], GLONASS [全球导航卫星系统]或GALILEO
G01S19-40 ...校正位置、速度或姿态
G01S19-42 ...确定位置
G01S 无线电定向;无线电导航;采用无线电波测距或测速;采用无线电波的反射或再辐射的定位或存在检测;采用其他波的类似装置
G01S19-00 卫星无线电信标定位系统;利用这种系统传输的信号确定位置、速度或姿态
G01S19-01 .传输时间戳信息的卫星无线电信标定位系统,例如,GPS [全球定位系统]、GLONASS[全球导航卫星系统]或GALILEO
G01S19-38 .利用卫星无线电信标定位系统传输的信号来确定导航方案
G01S19-39 ..传输带有时间戳信息的卫星无线电信标定位系统,例如GPS [全球定位系统], GLONASS [全球导航卫星系统]或GALILEO
G01S19-40 ...校正位置、速度或姿态
G01S19-42 ...确定位置