一种神经网络辅助粒子滤波的gps接收机自主完好性监测方法
【专利摘要】一种神经网络辅助粒子滤波的GPS接收机自主完好性监测方法，所述方法包括如下步骤：第一步，建立BP神经网络改进粒子滤波算法模型；第二步，建立似然比检验统计量；第三步，使用NNWA-PF和LLR算法对GPS故障卫星进行检测。本发明的优点在于：1、本方法通过神经网络对粒子进行权值分裂和权值调整，增加了有效粒子数目，改善了粒子多样性，抑制了粒子退化和多样性丧失。2、本方法满足非高斯噪声下的GPS接收机自主完好性监测需求。3、本方法采用似然比方法建立检验统计量，利用实测数据进行验证，更符合实际问题，方案的可行性更强。
【专利说明】-种神经网络辅助粒子滤波的GPS接收机自主完好性监测 方法
【技术领域】 [0001] ；本发明设及卫星导航定位结果的可靠性监测技术，是航空卫星导航接 收机定位信息可靠性监测的关键技术，是一种Gl^s接收机自主完好性监测方法，尤其是一 种神经网络辅助粒子滤波的GI^S接收机自主完好性监测方法。
【背景技术】 [0002] ：目前，卫星导航已被应用于军事、航海、航空、测量、交通、勘测等几乎 一切与位置、速度、时间有关的人类活动中。作为世界上最为成熟的Gl^s系统其应用也越 来越广泛，在各种应用领域中，航空卫星导航是对卫星导航系统可靠性要求极为苛刻的高 端应用用户，因此，在航空卫星导航中，完好性监测是航空导航系统不可或缺的一部分，由 于系统故障而引起的导航失败或者故障应该被探测和隔离W保证导航信息可靠性不受影 响。随着全球导航卫星系统（GNSS)的发展W及用户对于GNSS服务性能要求的不断升级， 完好性监测日益受到重视。由于通过卫星导航系统自身的操作控制部分进行卫星故障监测 时，告警时间比较长，通常在15分钟到几小时，不能满足航空导航需求，因此在用户端对卫 星故障进行快速监测，即接收机自主完好性监测（RAIM，Receiver Autonomous Integrity Monitoring)成为研究的热点。
[000引 目前，RAIM算法包括两大类；一类是利用当前伪距观测量的快照（Snapshot)算 法；另一类是滤波方法，包括基于卡尔曼滤波的RAIM算法W及基于粒子滤波的RAIM算法。 该些方法存在的问题主要表现在；（1)要求测量噪声服从高斯分布，而在实际中测量噪声 很难严格服从高斯分布，此时，该算法的性能就会降级。（2)对渐变微小偏差故障检测不灵 敏。（3)非线性非高斯粒子滤波方法存在粒子退化现象和粒子的多样性丧失等问题。
【发明内容】
[0004] ；针对上述现有技术的不足，本发明提供了一种神经网络辅助粒子滤波 的Gl^s接收机自主完好性监测方法。
[0005] 为实现上述目的，本发明采用的技术方案是：一种神经网络辅助粒子滤波的GPS 接收机自主完好性监测方法，所述方法包括如下步骤：
[0006] 第一步，建立BP神经网络改进粒子滤波算法模型；
[0007] 第二步，建立似然比检验统计量；
[000引第S步，使用NNWA-PF和LLR算法对GPS故障卫星进行检测。
[0009] 本发明的优点在于：
[0010] 1、本方法通过神经网络对粒子进行权值分裂和权值调整，增加了有效粒子数目， 改善了粒子多样性，抑制了粒子退化和多样性丧失。
[0011] 2、本方法满足非高斯噪声下的GI^S接收机自主完好性监测需求。
[0012] 3、本方法采用似然比方法建立检验统计量，利用实测数据进行验证，更符合实际 问题，方案的可行性更强。
【具体实施方式】：
[0013] 下面对本发明实施例中的技术方案进行清晰、完整地描述，显然，所描述的实施例 仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通 技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范 围。
[0014] 首先介绍如下的概念与定义：
[0015] 1、接收机自主完好性监测；接收机自主完好性监测是一种基于一致性检验的自主 完好性监测方法，它通过对冗余观测信息的检测来判断测量结果中是否含有较大卫星距离 误差，当系统性能超出指定的容差级时，它能实时有效地给用户提供告警信息。
[0016] 2、粒子滤波：粒子滤波（Particle Filter)算法又称为序贯蒙特卡罗方法，它依 据于大数定理，采用蒙特卡罗方法来求解贝叶斯估计中的积分运算。粒子滤波算法首先依 据系统状态向量的经验分布条件分布在状态空间产生一组随机样本的集合，然后根据观测 量不断地调整粒子的权重和位置，通过调整后粒子的信息修正最初的经验条件分布。当样 本容量很大时，该种蒙特卡罗描述就近似于状态变量真实的后验概率密度函数
[0017] 3、神经网络；BP神经网络是一种有教师信号的按误差逆传播算法训练的多层前 馈网络，是目前应用最广泛的神经网络模型之一。BP网络能学习和存贬大量的输入-输出 模式映射关系，而无需事前揭示描述该种映射关系的数学方程。它的学习规则是使用最速 下降法，通过反向传播来不断调整网络的权值和阔值，使网络的误差平方和最小。BP神经网 络模型拓扑结构包括输入层（input layer)、隐层化ide layer)和输出层（output layer)。 [001引本发明提供的GI^S接收机自主完好性监测方法，具体步骤如下：
[0019] 第一步，建立BP神经网络改进粒子滤波算法模型；
[0020] 假设描述动态系统的状态方程和观测方程为：
[0021] Xk=fk 狂…Vk_i)
[0022] Zk=hk(Xk, rik)
[002引其中，Xk为状态向量，Z巧测量向量，f巧状态转移函数，h巧状态向量和观测向 量之间的传递函数，Vk_i为系统噪声，n k为观测噪声。
[0024] ①粒子集初始化，k = 0:
[0025] 根据先验概率密度pOg产生随机样本，知"，X。?，…X。? (N为随机样本数）。
[0026] ②当k = 1，2，…时，执行W下步骤：
[0027] 1)状态预测；
[002引根据系统的状态方程产生k时刻的先验粒子：
[0029] 找kik-i(i) ;i = 1，2,…，N}?p〇(k|Xk_i)。
[0030] 2)更新；
[0031] 首先，进行权值更新。在获得测量值之后，根据系统的观测方程并利用下式计算粒 子的权值
[003引 0k")= ? k-i")P狂klXk(")，i = 1，…N
[0033] 归一化权值为：
[0034]
【权利要求】
1. 一种神经网络辅助粒子滤波的GPS接收机自主完好性监测方法，所述方法包括如下 步骤： 第一步，建立BP神经网络改进粒子滤波算法模型； 第二步，建立似然比检验统计量； 第三步，使用NNWA-PF和LLR算法对GPS故障卫星进行检测。
2. 如权利要求1所述的一种神经网络辅助粒子滤波的GPS接收机自主完好性监测方 法，其特征在于，所述的第一步建立BP神经网络改进粒子滤波算法模型具体过程为：假设 描述动态系统的状态方程和观测方程为： Xk=fkdVh) Zk=hk(Xk，nk) 其中，Xk为状态向量，Zk为测量向量，fk为状态转移函数，hk为状态向量和观测向量之 间的传递函数，Vlrf为系统噪声，nk为观测噪声； ① 粒子集初始化，k= 0: 根据先验概率密度P(Xci)产生随机样本，Xci'Xtl'…Xtl^ (N为随机样本数）； ② 当k= 1，2，…时，执行以下步骤： 1) 状态预测： 根据系统的状态方程产生k时刻的先验粒子： {Xk丨H(I);i= 1，2，...，N}?P(X1iIXh); 2) 更新： 首先，进行权值更新，在获得测量值之后，根据系统的观测方程并利用下式计算粒子的 权值《kU): ?k(i)= ?H⑴P(Zk|Xk⑴)，i= 1，…N 归一化权值为：
3) 粒子分裂： 将粒子矩阵按权值大小排序，分为高权值矩阵和低权值矩阵，将高权值矩阵中q个权 值非常大的粒子分裂为两个小的、权值减半的粒子，同时将低权值矩阵中的q个权值最小 的粒子舍弃； 4) 权值调整： 将第3)步的粒子权值矩阵进行降序排列，取其中最小的q个粒子，利用BP神经网络调 整粒子的权值；
式中，q为输入/输出层神经元的数量，Zk为教师信号，即系统在该时刻的量测值，yk为 神经网络的输出，输入数据为具有较小权值的粒子，其状态值作为BP神经网络的输入，粒 子权值作为神经网络的初始权值，样本的学习函数为系统量测方程；
5) 重采样： 然后，计算有效粒子数并与设定的有效粒子数阈值Nth_进行比较；如果 则进行重采样，随机设定一个小于1的阈值，通过对粒子权值累加，当累加 值超过这个阈值时就保留该粒子，重复N次，选取了N个粒子，最后将原来的带权样本 丨f=1映射为等权样本试丨，，1匕，否则，执行下面的步骤； 6) 估计： 计算当前时刻系统的状态估计值：
7) 时刻k=k+1，继续计算。
3. 如权利要求1所述的一种神经网络辅助粒子滤波的GPS接收机自主完好性监测方 法，其特征在于，所述的第二步建立似然比检验统计量具体过程为； 将LLR测试定义为各辅助PF和主PF的概率密度函数之比，其中主PF为利用所用的测 量值进行滤波，各辅助PF为分别去掉一个测量值后由其他测量值进行的滤波，表达式如下 式所示：
由于系统状态估计的似然函数可用PF中粒子的归一化权值近似表示，故上式中的pH)和pA(yiIYiJ可表示为：
因此，计算出主PF和各辅助PF中各时刻的粒子归一化权值，即可得到用于一致性检测 的累加LLR，从而对系统是否出现故障进行检测。
4. 如权利要求1所述的一种神经网络辅助粒子滤波的GPS接收机自主完好性监测方 法，其特征在于，所述的第三步使用NNWA-PF和LLR算法对GPS故障卫星进行检测具体过程 为：利用GPS接收机采集卫星位置信息和伪距数据，建立载体状态方程和GPS接收机测量方
误判决（T为判决阈值），如果0 k>T，错误告警设置为ta=t时刻并跳至下面的步骤（1); 如果ek<T，则表示没有错误，跳至步骤（2); ⑴错误隔离：在k>ta下，取出Q颗卫星中累积LLR最大的那个卫星q即为出错的卫 星，则g=q; (2)状态更新：粒子滤波重采样得到更新粒子。
【文档编号】G01S19/20GK104502922SQ201410749171
【公开日】2015年4月8日 申请日期:2014年12月9日 优先权日:2014年12月9日
【发明者】王尔申, 李兴凯, 庞涛 申请人:沈阳航空航天大学