基于MHSS的ARAIM算法的详细分析解释
RAIM和ARAIM的区别
完好性支持信息ISM
- 用户测距误差(User Range Error, URE):用来计算精度和连续性的标称误差高斯分布的标准差;
- 用户测距精度(User Range Accuracy, URA):用来计算完好性的标称误差高斯分布的标准差;
- 卫星故障率Psat,星座故障率Pconst;(这些都是先验概率)
- 连续性偏差(bnom):估计精度和连续性的标准状况下典型偏差大小;
- 完好性偏差(bmax):用于无故障条件下完好性评估的最大偏差;
所有需要监测的故障模式(子集)的确定
- 故障数指某个故障模式下(也就是某个子集下)发生故障的卫星/星座数。
- 对于卫星故障或者星座故障,统称为事件,总的卫星数记为,总的星座数记作,那么总的事件数
- 每一种故障模式可以理解为不同的事件组合而成的
这一段是我自己的理解:
因为我们如果考虑可能出现的最大故障数越大,那么子集的数量就越大。但是故障数越大时,就代表其发生概率很低,几乎可以不用去考虑这种情况的发生。因此我们定义了一个最大故障数.
最大故障数的计算方法:
MHSS多假设解分离算法
什么是子集?
即全可见星集合中排除掉某颗或某几颗卫星的卫星集合。
MHSS的基本原理
全可见星解为使用全可见星集合定位的结果,子集解为子集卫星定位的结果。算法认为,在无故障情况下,全可见星解和所有的子集解应聚集一起。而如果某颗卫星存在故障,则使用该故障卫星测量值的全可见星定位解和子集定位解将产生偏移,不含故障卫星的子集定位解将更接近飞机实际位置。通过对比各个子集解和全可见星解的距离即可判断是否存在故障。
若所有子集解距离全可见解小于预设阈值,则认为无故障存在,若有子集解的距离超出阈值,则认为存在故障。从而采用排除算法,实现故障卫星的移除,最后确保至少有一组子集是由全部无故障卫星组成的。
MHSS算法公式推导
线性化伪距观测方程 的推导
设可见星的个数为n,可得以下n几个方程:
..............
1
其中为用户接收机的位置坐标,属于未知量,也是本方程组想要计算的量,为n颗可见性的位置坐标(可从卫星星历计算得到),为各个卫星钟差导致的伪距偏差(属于误差项,也是未知因素),是通过卫星和接收机时钟乘以光速可计算得出的已知量。
由上面的公式,根据最小二乘法可以计算的解——, 即的估计值。(真实值是计算得不到的,只能计算一个估计值)(本文中上面带有尖号的量都是估计量)
得到后,重新构建以下n个方程:
2
公式2可以这么理解,通过公式1的计算, = 右边的都是已知量了,从而可以计算出一个。是通过估计的用户接收机坐标计算到的每一颗可见星的伪距。
定义
,即真实伪距与伪距估计值之差
= 即接收机真实左边减去估计坐标
上面两个都是误差量,属于未知量,我们假设无故障情况下他们都是高斯分布。
根据上面的定义,公式2可以改写成:
3
将公式3进行傅里叶展开,得到n个方程组:
4
公式4用矩阵形式描述为:
5
公式5即为:
6
其中,为测量误差矢量;G为观测矩阵,由接收机与卫星间方向余弦及接收机钟差相关系数组成;y为观测矩阵与线性化伪距预测值之间的差值;x 表示接收机三维位置及时钟与标称值的偏差。
y、x、 都认为是服从高斯分布。
公式6可以理解为伪距观测误差和接收机位置误差相联系起来的方程。
故障识别算法推导
利用最小二乘法,可获得n颗可视卫星下的接收机坐标的估计值:
= 7
=
- 公式7中,为无故障假设下加权最小二乘投影矩阵,可以理解为把 y 投影成 (这个y的含义还不明确);
- 为n×n的加权对角矩阵(只有对角线有元素),其对角线元素是URA的函数,= ,为完好性标称误差对角协方差矩阵,其计算方法如下所示:
8
- 同样也是只有对角线有元素
- 是用于完好性的星历星钟误差的标准差,这个我觉得可以从完好性支持信息ISM的用户测距精度(User Range Accuracy, URA)计算得来
- 和分别是对流层延迟误差和接收机噪声误差。(这个我觉得也是通过一些给定信息计算而来)
相应地,可以获得排除第k颗卫星条件下接收机坐标的估计:
9
- 为第k个对角元素置零的n×n维单位矩阵;
因此,第 k 颗卫星的检验统计量如下:
10
- 每一颗卫星都有自己的检验统计量,用于检验是否有故障
检测门限如下:
11
- ,这里的Sk和S0在公式7、9中,后面的(3,i)我觉得是第三行第i列的意思
- 为垂直方向第 k 颗卫星的检验统计量 ݀ 的标准差,,,其中是n×n的加权对角阵,其对角线元素是URE的函数,,为用于精度和连续性的标称误差对角协方差矩阵,其计算方法如下所示:,是用于精度和连续性的星历星钟误差的标准差,可以从完好性支持信息ISM的用户测距误差(User Range Error, URE)计算得来
- ,为总误警概率(性能要求时给定的),Q为标准正态分布的尾部概率。
从北航《基于子集包含减少 ARAIM 子集数量的方法》文献中看到,检测门限还分东北天三个不同方位的检测门限。三个方位的检测门限在计算时只需将改成各自方向上被分配的误警概率即可。
- 这里的为标称误差,是ISM的完好性偏差信息
对比检验统计量和检测门限,如果,则该卫星存在故障。
求解保护级PL和有效监视门限EMT
.............
SAFETY CRITICAL BOUNDS FOR PRECISE POSITIONING FOR AVIATION AND AUTONOMY 论文
2.3节
- PHMI:危险误导性信息的概率或完整性风险分配
- x:真实的状态
- x尖号:估计状态
- PL:保护级
- q:检验统计量
- T:一个基于误警要求和子集协方差的阈值
本质上,完整性风险是位置误差超过保护水平且测试统计数据低于阈值的概率。
ARAIM方法使用多层次解决方案分离(MHSS)将完整性风险分解为基于故障假设的组件,因此方程可以重写:
4.5节
这里提到的子集过滤器我的理解是一种新的分组策略。
FDE(故障检测排除)原始方法:
第一种方法涉及在故障检测和排除后对AIV(全视图)和子集过滤器进行简单的重新初始化。计算的唯一子集是在当前时期生成保护级别所需的子集。下图的左半部分显示了视图中所有(顶行)和其他子集(第2行之后)的测量使用示意图。
如上图右侧所示。新AIV实际上是原始集合的一个子集,但由于之前未跟踪任何新子集,因此无法使用已跟踪的AIV,必须重新初始化整个过滤器。无法使用以前跟踪的AIV的原因是,如果在给定已收敛AIV的情况下生成新子集,则这些子集将无法容忍AIV解决方案中已包含的测量上的错误。一旦过滤器重新初始化,保护级别将达到峰值,并且由于缺少PRN而导致几何结构减弱,保护级别需要比平常更长的时间才能再次收敛。
FDE(故障检测排除)新方法——one-deck filters:
为了避免在FDE后重新初始化过滤器时保护级别出现大幅跳变,开发了一种新工艺,仅在FDE情况下使用。本质上,如果我们只关注排除1颗卫星的子集,那么新FDE方法的子集将是排除两个卫星的子集,如上图图左侧所示。这些子集不包括在保护级别计算中。这样的话,检测到故障时,过滤器不需要重新初始化,因为新的AIV和子集一直都被跟踪。
新的过滤器由上图右侧的红色框突出显示。此时需要重新初始化新的one-deck上的过滤器,但只要在这些过滤器收敛(我觉得指的是保护级的收敛之前,可以看论文的4.4.1保护级收敛的图)之前未检测到其他故障,这不会影响保护级别。新的AIV是原始子集过滤器之一,新的主要子集已经在one-deck过滤器中收敛。当然,携带甲板上的过滤器在计算上是昂贵的,但是已经探索了子集分组等技术[12],[16],这可以显著降低计算负载。FDE后,由于几何结构减弱,保护级别略有增加,但最重要的是,它们已经收敛。
使用新的FDE方法,注入故障后保护级和定位误差收敛图
基于MHSS的ARAIM算法的详细分析解释相关推荐
- Yolov1目标检测算法详细分析
Yolov1目标检测算法详细分析 Yolov1介绍 这是继RCNN,fast-RCNN 和 faster-RCNN之后,rbg(Ross Girshick)大神挂名的又一大作,起了一个很娱乐化的名字: ...
- ML之NB、LoR:基于NB和LoR算法对Kaggle IMDB影评数据集(国外类似豆瓣电影)情感分析进行分类
ML之NB.LoR:基于NB和LoR算法对Kaggle IMDB影评数据集(国外类似豆瓣电影)情感分析进行分类 注:LoR code waiting 目录 输出结果 核心代码 输出结果 数据集:Dat ...
- yolo算法的优缺点分析_YOLO算法详细解析(一)
目标检测和目标分类 图像识别算法是计算机视觉的基础算法,例如VGG,GoogLeNet,ResNet等,这类算法主要是判断图片中目标的种类. 目标检测算法和图像识别算法类似,但是目标检测算法不仅要识别 ...
- 基于Qt的A*算法可视化分析
代码地址如下: http://www.demodashi.com/demo/13677.html 需求 之前做过一个无人车需要自主寻找最佳路径,所以研究了相关的寻路算法,最终选择A算法,因为其简单易懂 ...
- PDERL:基于DEM的快速精确通视域分析算法介绍(一)
PDERL: 基于DEM的快速精确通视域分析算法(精度与R3相同,速度与XDraw接近) 系列文章: PDERL:基于DEM的快速精确通视域分析算法介绍(一) PDERL:基于DEM的快速精确通视域分 ...
- NLP之TEA之NB/LoR:基于NB和LoR算法对Kaggle IMDB影评数据集(国外类似豆瓣电影)情感分析进行分类
NLP之TEA之NB/LoR:基于NB和LoR算法对Kaggle IMDB影评数据集(国外类似豆瓣电影)情感分析进行分类 注:LoR code waiting 目录 输出结果 核心代码 输出结果 数据 ...
- ML/DL之预测分析类:利用机器学习算法进行预测分析的简介、分析、代码实现之详细攻略
ML/DL之预测分析类:利用机器学习算法进行预测分析的简介.分析.代码实现之详细攻略 目录 机器学习算法进行预测的简介 机器学习算法进行预测的分析 机器学习算法进行预测的代码实现 机器学习算法进行预测 ...
- java消费者模式_基于Java 生产者消费者模式(详细分析)
生产者消费者模式是多线程中最为常见的模式:生产者线程(一个或多个)生成面包放进篮子里(集合或数组),同时,消费者线程(一个或多个)从篮子里(集合或数组)取出面包消耗.虽然它们任务不同,但处理的资源是相 ...
- AMAP-TECH算法大赛开赛!基于车载视频图像的动态路况分析
简介:车载视频图像包含了更多的信息量,给了我们另外一个解决问题的视角.通过视频或图片,可以观察到路面的真实状态,包括机动车数量.道路宽度和空旷度等等.基于车载视频图像可以获取更准确的路况状态,为用户出 ...
- 算法设计与分析 0-1背包问题 动态规划解法【超详细】
0-1背包问题 问题描述 给定i个物品和一个容量为的背包,物品的重量是Wi,其价值为Vi 物品个数为i,背包容量为C. 如何选择装入背包内的物品,使得装入背包中的物品的总价值最大? 其中,每种物品只有 ...
最新文章
- vue怎么引入外网json文件
- 前期绑定 vs 后期绑定
- 数据中心行业在能源转型中将发挥重要作用
- docker基础用法
- poj 1751 Highways 最小生成树Kruskal(、Prim还没写
- DMA(direct memory access)控制方式
- OpenCV-理想高通低通滤波器(C++)
- mysql的cost
- 百度之星冠军分享:AI图像赛事入门
- 计算机网络-名词解释整理
- 支付宝支付接口的使用
- PDF文件在线转换HTML的方法
- 盘点CSV文件在Excel中打开后乱码问题的两种处理方法
- laravel与TP的区别
- mysql拼图游戏代码_HTML+Javascript制作拼图小游戏详解(二)
- Python—网络编程_Mail
- 基于高德导航的Android大作业
- C# 终本案件、综合执行人、裁判文书爬虫
- 服务器宝塔Error: connect ETIMEDOUT
- Visual Studio Codec#
热门文章
- C语言常用库函数(含详细用法))
- 用matlab如何画六边形,matlab怎样直接画出六边形
- a form 出口享惠情况_关税聚焦 | 报关单“出口享惠情况”填报解读
- 无root权限安卓真机抓包HTTPS
- linux ulipad,Python for Ulipad(总结ulipad使用)
- 【python教程入门学习】利用Python自动生成暴力破解的字典
- Android 系统架构图
- 相亲中的最优停止理论-相亲中的数学
- Linux安装nginx做端口映射
- 易语言代码转php,易语言代码转PHP代码有没大佬