CN108020847B - 针对高级接收机自主完好性监测中故障模式的确定方法 - Google Patents
针对高级接收机自主完好性监测中故障模式的确定方法 Download PDFInfo
- Publication number
- CN108020847B CN108020847B CN201711204949.3A CN201711204949A CN108020847B CN 108020847 B CN108020847 B CN 108020847B CN 201711204949 A CN201711204949 A CN 201711204949A CN 108020847 B CN108020847 B CN 108020847B
- Authority
- CN
- China
- Prior art keywords
- fault
- probability
- failure mode
- failure
- event
- 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.)
- Active
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/20—Integrity monitoring, fault detection or fault isolation of space segment
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Security & Cryptography (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公布了一种针对高级接收机自主完好性监测中故障模式的确定方法;对高级接收机自主完好性监测中最大故障数、故障模式、监测子集和未监测故障概率的定义,提出了一种概率累加式的故障模式确定方法,将故障模式的确定与最大故障数、监测子集的概率和未监测故障概率的计算结合在一起进行求解,本方法将不可观测的监测子集的判定和排除嵌入到子集的遍历当中,简化了未监测故障概率的计算方法,排除了不可观测子集移除后剩余已监测子集的概率总和低于完好性风险阈值的隐患;同时本方法得到的检测子集数量更少,已监测故障概率可以更为精确的满足完好性风险的监测阈值,缓解了计算负担;此外,本方法省略了最大故障数的计算过程。
Description
技术领域
本发明涉及卫星导航技术领域,特别涉及一种针对高级接收机自主完好性监测中故障模式的确定方法。
背景技术
接收机自主完好性监测(Receiver Autonomous Integrity Monitoring,RAIM)是卫星导航接收机根据冗余的GNSS信息自主进行故障诊断与排除的方法。RAIM算法是包含在接收机里面的,因此被称为自主监测,RAIM也是最直接、最及时、应用最广泛、研究最深入、计算效率最高的完好性监测方法。RAIM技术仅支持侧向导航,无法满足国际民航组织规定的垂直引导航向性能定义的200英尺以下垂直引导航道性能需求(Localized Performancewith Vertical guidance,LPV-200)。高级接收机自主完好性监测技术(Advanced RAIM,ARAIM),是由GNSS发展构架研究小组(GNSS Evolutionary Architecture Study,GEAS)设计主要用于在2030年前为航空LPV-200操作提供完好性监测的解决方案。
ARAIM由空间段、用户段和地面段组成,ARAIM的概念可以大致概括为用户接收机接受空间段多个星座的多频信号测量值,根据地面监测站提供的完好性支持信息(Integrity Support Message,ISM),确定出需要进行监测的故障模式以及对应的监测子集的概率,用户算法计算各个子集的位置估计和完好性边界,从而识别和排除故障测量值,得到定位解的保护水平。ARAIM可以充当一种协调各个星座完好性增强平台,缓解各个星座的差错并能融合他们的完好性方面的性能,对个别星座的负面变化不敏感。
ARAIM的用户接收机算法是基于多假设解分离(Multiple Hypothesis SolutionSeparation,MHSS)的方法,在此之前第一步是根据完好性支持信息ISM和接收机定位解算得到的卫星数、星座数等参数,确定故障模式的数量和对应监测子集的概率。故障模式选择得当,既可以满足完好性风险的阈值要求,又可以尽可能减少多假设解分离的子集数量,缓解计算压力。
高级接收机自主完好性监测基础算法(baseline algorithm)通过数学估算方法得到最大故障数和初始已监测故障概率,然后判断监测子集的可观测性,排除掉不可观测的子集概率后得到最终已监测故障概率。这样的处理流程存在两个问题:一是数学估算方法得到的初始已监测故障概率不够精确,是充分条件而不是必要条件;二是先计算已监测故障概率,后进行可观测性判断,存在去掉不可观测的子集概率后已观测故障概率不能满足完好性风险阈值要求的隐患。
发明内容
本发明所要解决的技术问题是高级通过概率累加的方法解决背景技术中存在的两个问题,且该发明适用于导航接收机自主完好性监测。
本发明为实现上述目的,采用如下技术方案:
针对高级接收机自主完好性监测中故障模式的确定方法,具体包含如下步骤:
步骤1,导航接收机接收卫星信号、完好性支持信息ISM和完好性风险阈值PTHRES;
步骤2,导航接收机根据卫星信号进行定位解算,获取参与定位的卫星数量Nsat和星座数量Nconst;
步骤3,将Nsat颗卫星和Nconst个星座分别发生故障的事件定义为相互不受干扰的独立事件,进而从完好性支持信息ISM中提取对应的故障事件概率Pevent,k,其中1≤k≤(Nsat+Nconst);
步骤5,选择一个故障模式Fevent_temp,生成和故障模式Fevent_temp对应的监测子集;
判断故障模式Fevent_temp是否可观测:如果是,则进入步骤6,反之,则进入步骤9;
步骤6,计算故障模式对应监测子集的概率,并保存至数组pfault中;
步骤7,将步骤6得到的概率累加到已监测故障事件概率Psum中,即Psum=Psum+pfault,则监测子集的数量Nfault加一,同时保存故障模式至序列Idxfault中;
步骤8,判断已监测故障事件概率Psum是否已经超过阈值(1-PTHRES),如果是,则进入步骤11;如果不是,则进入步骤9;
步骤9,判断故障模式选择序列IDX是不是已经遍历完,如果是,则进入步骤10;如果不是,则进入步骤5;
步骤10,故障阶数Nfault_order加一,进入步骤5;
步骤11,计算未监测的故障模式概率Pnot_monitored,将遍历结束时所对应的故障阶数Nfault_order作为最大故障数,进而输出pfault,Nfault,Idxfault和Pnot_monitored。
作为本发明针对高级接收机自主完好性监测中故障模式的确定方法的进一步优选方案,在步骤4中,故障模式选择序列IDX具体表示如下:
作为本发明针对高级接收机自主完好性监测中故障模式的确定方法的进一步优选方案,在步骤4中,故障模式概率由大到小进行排序是在同一故障阶数内,概率大的故障模式在前,概率小的故障模式在后。
作为本发明针对高级接收机自主完好性监测中故障模式的确定方法的进一步优选方案,在步骤5中,从IDX中按照前后顺序选择一个故障模式Fevent_temp。
作为本发明针对高级接收机自主完好性监测中故障模式的确定方法的进一步优选方案,未监测故障模式的概率具体计算如下:
Pnot_monitored=1-Psum
本发明与现有技术相比有益效果体现在:
1、本发明将故障模式、最大故障数、已监测故障概率、未监测故障概率的计算紧密联系在一起,同时求解完成;
2、本发明通过概率累加的方法获得已监测故障概率,每次累加后都与阈值进行判断,得到的已监测故障概率更加精确,相比数学估算法方法也更具物理意义,得到的监测子集更少,缓解了计算负担;
3、本发明将监测子集的可观测性判断嵌入到故障模式的遍历当中,先判断子集是否可观测,通过可观测性判断后再计算子集概率并累加至已监测故障概率中,排除了最后进行子集可观测性判断导致的已监测故障概率不能满足阈值要求的隐患;
4、本发明省略了最大故障数的计算过程,故障模式遍历结束时的故障阶数即为最大故障数。
附图说明
图1是本发明针对接收机自主完好性监测的故障模式确定方法的具体流程图;
图2是本发明故障模式累加示意图。
具体实施方式
下面结合附图对发明的技术方案进行详细说明:
具体实施例如下:如图1所示,针对高级接收机自主完好性监测中故障模式的确定方法,具体需要完成以下工作:
一,基于概率累加的最大故障数计算法:
本发明采用故障累加的方式,按照“故障数由低到高+故障概率由大到低”的顺序,依次累加故障模式的概率,直到满足完好性风险阈值。“故障数由低到高”的意思是从无故障到一个故障到两个故障到更多故障,如图2所示。
“故障概率由大到低”的意思是在同一故障阶数内,先计算故障概率大的事件,后计算故障小的事件。具体方法是将Nsat颗卫星和Nconst个星座分别发生故障的事件定义为相互不受干扰的独立事件,在一个故障阶数内,共有个故障模式,Nfault_order是故障阶数,Nfault_order初始值为1。对个事件的概率按照由大到小顺序进行排序,生成故障事件概率选择序列IDX。
其中,Fmode,i代表第i个故障模式,代表第i个故障模式的概率,Pevent,j表示故障模式中故障事件的概率,i,j都是正整数。按照选择序列IDX,依次选择故障模式进行可观测性判断。当已监测故障概率满足阈值要求时,遍历结束,当时的故障阶数就是最大故障数。
二,监测子集的不可观测性:
通过遍历故障模式得到的监测子集,存在不可观测性,如果子集无法通过可观测性判断,则需要剔除掉,其概率不能累加至已监测故障概率中,不可观测性发生的情况包括:子集里只有三颗或更少的卫星,且这些卫星属于同一个星座;子集里只有四颗或更少的卫星,且这些卫星属于两个或更多的星座;子集里只有五颗或更少的卫星,且这些卫星属于三个或更多的星座;子集里只有六颗或更少的卫星,且这些卫星属于四个星座。
基础算法是通过估算最大故障数,得到未监测多故障概率Pmultiple,not_monitored,然后通过遍历各个监测子集得到不可观测子集并得到它们的概率总和Punobservable,通过两者求和得到不可监测故障概率:
Pnot_monitored=Pmultiple,not_monitored+Punobservable (2)
在保护水平的具体计算过程中,垂直保护水平需要满足公式(3):
水平保护水平需要满足公式(4):
其中q=1,2,最终得到的水平保护水平为:
由公式(3)和公式(4)可知,未监测概率是作为一个整体参与计算的,所以分别求解未监测多故障概率和不可观测故障概率并没有必要,本方法提出的未监测概率求解方法更简单直接。
三,计算故障模式对应监测子集的概率:
每颗星的故障概率和每个星座的故障概率可以从完好性支持信息ISM里面得到。Nsat颗卫星和Nconst个星座分别发生故障的事件为相互不受干扰的独立事件,定义每个事件发生故障的概率为Pevent,k,其中1≤k≤(Nsat+Nconst)。无故障概率为:
故障阶数为1时,监测子集的概率为:
故障阶数为2时,监测子集的概率为:
故障阶数为3或更多时,对应监测子集的概率以此类推。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下还可以作出若干改进,这些改进也应视为本发明的保护范围。
Claims (5)
1.针对高级接收机自主完好性监测中故障模式的确定方法,其特征在于:具体包含如下步骤:
步骤1,导航接收机接收卫星信号、完好性支持信息ISM和完好性风险阈值PTHRES;
步骤2,导航接收机根据卫星信号进行定位解算,获取参与定位的卫星数量Nsat和星座数量Nconst;
步骤3,将Nsat颗卫星和Nconst个星座分别发生故障的事件定义为相互不受干扰的独立事件,进而从完好性支持信息ISM中提取对应的故障事件概率Pevent,k,其中1≤k≤(Nsat+Nconst);
步骤5,选择一个故障模式Fevent_temp,生成和故障模式Fevent_temp对应的监测子集;
判断故障模式Fevent_temp是否可观测:如果是,则进入步骤6,反之,则进入步骤9;
步骤6,计算故障模式对应监测子集的概率,并保存至数组pfault中;
步骤7,将步骤6得到的概率累加到已监测故障事件概率Psum中,即Psum=Psum+pfault,则监测子集的数量Nfault加一,同时保存故障模式至序列Idxfault中;
步骤8,判断已监测故障事件概率Psum是否已经超过阈值(1-PTHRES),如果是,则进入步骤11;如果不是,则进入步骤9;
步骤9,判断故障模式选择序列IDX是不是已经遍历完,如果是,则进入步骤10;如果不是,则进入步骤5;
步骤10,故障阶数Nfault_order加一,进入步骤5;
步骤11,计算未监测的故障模式概率Pnot_monitored,将遍历结束时所对应的故障阶数Nfault_order作为最大故障数,进而输出pfault,Nfault,Idxfault和Pnot_monitored。
3.根据权利要求1所述的针对高级接收机自主完好性监测中故障模式的确定方法,其特征在于:在步骤4中,故障模式概率由大到小进行排序是在同一故障阶数内,概率大的故障模式在前,概率小的故障模式在后。
4.根据权利要求1所述的针对高级接收机自主完好性监测中故障模式的确定方法,其特征在于:在步骤5中,从IDX中按照前后顺序选择一个故障模式Fevent_temp。
5.根据权利要求1所述的针对高级接收机自主完好性监测中故障模式的确定方法,其特征在于:未监测故障模式的概率具体计算如下:
Pnot_monitored=1-Psum。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711204949.3A CN108020847B (zh) | 2017-11-27 | 2017-11-27 | 针对高级接收机自主完好性监测中故障模式的确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711204949.3A CN108020847B (zh) | 2017-11-27 | 2017-11-27 | 针对高级接收机自主完好性监测中故障模式的确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108020847A CN108020847A (zh) | 2018-05-11 |
CN108020847B true CN108020847B (zh) | 2021-05-28 |
Family
ID=62077353
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711204949.3A Active CN108020847B (zh) | 2017-11-27 | 2017-11-27 | 针对高级接收机自主完好性监测中故障模式的确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108020847B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109581445B (zh) * | 2018-11-01 | 2021-03-19 | 北京航空航天大学 | 一种基于北斗星座的araim子集选择方法及系统 |
CN110196434B (zh) * | 2019-03-29 | 2022-05-20 | 南京航空航天大学 | 一种高级接收机自主完好性监测的星座动态选择方法 |
CN110007317B (zh) * | 2019-04-10 | 2022-06-17 | 南京航空航天大学 | 一种选星优化的高级接收机自主完好性监测方法 |
CN112731460B (zh) * | 2019-10-28 | 2022-10-28 | 千寻位置网络有限公司 | 卫星导航系统校准数据的完好性风险的测试方法及其系统 |
CN111060133B (zh) * | 2019-12-04 | 2020-10-20 | 南京航空航天大学 | 一种用于城市复杂环境的组合导航完好性监测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6798377B1 (en) * | 2003-05-31 | 2004-09-28 | Trimble Navigation, Ltd. | Adaptive threshold logic implementation for RAIM fault detection and exclusion function |
CN101806903A (zh) * | 2010-04-01 | 2010-08-18 | 北京航空航天大学 | 一种用于卫星导航系统的接收机自主完好性监测方法 |
CN104267410A (zh) * | 2014-10-10 | 2015-01-07 | 北京航空航天大学 | 机载完好性监测中多故障的排除方法和装置 |
CN104316934A (zh) * | 2014-10-22 | 2015-01-28 | 北京航空航天大学 | Araim地面监测站中伪距加速度故障检测方法和装置 |
CN106646526A (zh) * | 2017-02-09 | 2017-05-10 | 南京航空航天大学 | 一种可同时检测识别多种故障的接收机自主完好性检测方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1637899A1 (en) * | 2004-09-20 | 2006-03-22 | EADS Astrium GmbH | Method and apparatus for providing integrity information for users of global navigation system |
EP2648018B1 (en) * | 2012-04-02 | 2014-10-15 | Airbus DS GmbH | An improved RAIM algorithm |
US9547086B2 (en) * | 2013-03-26 | 2017-01-17 | Honeywell International Inc. | Selected aspects of advanced receiver autonomous integrity monitoring application to kalman filter based navigation filter |
CN105954768B (zh) * | 2016-04-25 | 2018-03-30 | 北京航空航天大学 | 基于araim风险分配优化的ads‑b可用性预测方法和装置 |
-
2017
- 2017-11-27 CN CN201711204949.3A patent/CN108020847B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6798377B1 (en) * | 2003-05-31 | 2004-09-28 | Trimble Navigation, Ltd. | Adaptive threshold logic implementation for RAIM fault detection and exclusion function |
CN101806903A (zh) * | 2010-04-01 | 2010-08-18 | 北京航空航天大学 | 一种用于卫星导航系统的接收机自主完好性监测方法 |
CN104267410A (zh) * | 2014-10-10 | 2015-01-07 | 北京航空航天大学 | 机载完好性监测中多故障的排除方法和装置 |
CN104316934A (zh) * | 2014-10-22 | 2015-01-28 | 北京航空航天大学 | Araim地面监测站中伪距加速度故障检测方法和装置 |
CN106646526A (zh) * | 2017-02-09 | 2017-05-10 | 南京航空航天大学 | 一种可同时检测识别多种故障的接收机自主完好性检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108020847A (zh) | 2018-05-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108020847B (zh) | 针对高级接收机自主完好性监测中故障模式的确定方法 | |
CN110007317B (zh) | 一种选星优化的高级接收机自主完好性监测方法 | |
CN106646526B (zh) | 一种可同时检测识别多种故障的接收机自主完好性检测方法 | |
KR101574819B1 (ko) | 비정상적 의사-거리 측정에 대해 무선 내비게이션 수신기 사용자를 보호하는 방법 | |
CN110515097B (zh) | 应用于基准站的gnss卫星观测粗差剔除方法和装置 | |
CN108761498B (zh) | 一种针对高级接收机自主完好性监测的位置估计优化方法 | |
CN110196434B (zh) | 一种高级接收机自主完好性监测的星座动态选择方法 | |
US20090273511A1 (en) | Method of operating a satellite navigation receiver | |
CN104075734B (zh) | 水下组合导航故障智能诊断方法 | |
EP3012665B1 (en) | Systems and methods for averaging satellite sigmas and readmitting excluded satellite measurements into differential corrections and integrity monitors | |
CN110646820B (zh) | Rtk定位数据的质检方法、装置、设备和存储介质 | |
Feng et al. | Integrity monitoring for carrier phase ambiguities | |
CN104280741B (zh) | 电离层异常检测方法 | |
CN109061683B (zh) | 一种优化水平保护级的h-araim系统 | |
CN113009520B (zh) | 一种卫星导航矢量跟踪环的完好性检测方法 | |
CN112269193B (zh) | 一种监测站伪距观测序列辅助的导航卫星故障检测方法 | |
CN105158778B (zh) | 多系统联合实施载波相位差分故障卫星剔除方法及其系统 | |
CN104459722B (zh) | 一种基于多余观测分量的整周模糊度可靠性检验方法 | |
CN106842236A (zh) | Gnss接收机周跳探测与修复处理方法及装置 | |
CN112130177A (zh) | 一种基于稳定分布的地基增强系统完好性监测方法 | |
CN105511481B (zh) | 一种星载定轨优化方法 | |
CN105204047B (zh) | 一种卫星导航系统中观测量单个粗差的探测与修复方法 | |
CN105549046A (zh) | Gnss接收机周跳探测与修复处理方法 | |
CN113295140A (zh) | 一种基于虚拟参考站检测铁路沉降的方法及系统 | |
Tmazirte et al. | Fast multi fault detection & exclusion approach for GNSS integrity monitoring |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |