惯性卫星紧组合导航系统的多故障检测装置及检测方法

1.惯性卫星紧组合导航系统的多故障检测方法，其特征在于包括如下步骤：
步骤1，确定惯性卫星紧组合导航系统的观测量矩阵；
步骤2，建立故障检测的检验统计量：构造奇偶空间矩阵并确定观测方程表达式，将观测误差矩阵投影到奇偶空间矩阵得到奇偶矢量，再根据显著性校验原理建立检验统计量；
步骤3，检测识别惯性卫星紧组合导航系统的故障，具体包括如下步骤：
5
步骤3-1，对n颗可见卫星进行子集划分，以每5颗可见卫星为一个子集，共有Cn 个子集；
步骤3-2，在每个子集中，利用步骤2所述的方法建立各个可见卫星的检验统计量，计算各个可见卫星在所属子集中的故障概率；
步骤3-3，叠加可见卫星在所属各子集中的故障概率，对故障概率和做归一化处理得到各个可见卫星的质量；
步骤3-4，判定质量小于规定值的可见卫星为故障卫星；
步骤4，剔除故障卫星，重构观测系统，确定观测方程：
步骤4-1，计算可见卫星数与故障卫星数的差得到可用卫星数：
若可用星数大于或者等于4，进入步骤4-2；
若可用星数小于4，进入步骤4-3；
步骤4-2，删去观测量矩阵中关于故障卫星的数据；
步骤4-3，删去观测量矩阵中质量值小于0.4的故障卫星的数据；
步骤4-4，加权处理修正后的观测量矩阵得到无故障卫星信息的观测量矩阵，卡尔曼滤波器根据无故障卫星信息的观测量矩阵确定观测方程，再将卫星状态信息的误差值反馈给惯性导航装置和GPS接收机。
2.根据权利要求1所述的惯性卫星紧组合导航系统的多故障检测方法，其特征在于所述步骤1的具体实施如下：
步骤1-1，惯导等效伪距测量单元根据惯性导航装置的位置速度信息、原始卫星的位置信息计算惯导与可见卫星的伪距；
步骤1-2，加法器用接收机输出的可见卫星伪距减去步骤1-1得到的惯导与伪距得到原始观测量矩阵；
步骤1-3，惯导辅助装置根据GPS输出的星历文件选择可见卫星之外的一颗卫星为辅助卫星；
步骤1-4，根据惯性导航装置测到的惯导位置信息、GPS接收机提供的原始卫星位置信息、GPS天线与惯导装置中心之间的校正距离，计算辅助伪距向量，将辅助伪距向量加入到步骤1-2得到的原始观测量矩阵；
步骤1-5，利用步骤1-4所得的观测量矩阵计算可见卫星、辅助卫星到GPS接收机的方向余弦确定观测量矩阵。

惯性卫星紧组合导航系统的多故障检测装置及检测方法\n技术领域\n[0001] 本发明涉及惯性卫星紧组合导航系统的多故障检测装置及检测方法，属于卫星定位导航的技术领域。\n背景技术\n[0002] 惯性卫星紧组合导航系统是将惯性导航与卫星导航相结合，采用卫星原始观测信息伪距、伪距率作为观测量的导航系统。相比于松组合，紧组合系统具有较好的精度及动态性能，并且在可见卫星少于4颗的情况下，也能在较短的时间内正常工作。\n[0003] 但由于紧组合系统组合程度较深，其精度及可靠性受到卫星原始观测信息伪距、伪距率质量的影响程度较大。伪距、伪距率容易受到如下因素影响而产生故障：\n[0004] （1）太空中摄动因素、卫星轨道误差、电离层、对流层及空间信号传输畸变等因素影响；\n[0005] （2）卫星导航系统的软、硬件随机故障；\n[0006] （3）人为干扰，特别是敌意干扰，会带来较大误差。\n[0007] 同时，随着各国对卫星导航系统的重视程度和投入力度不断加大，在轨导航卫星数目也在不断增加，多通道接收机输出的原始测量信息同时发生故障的概率也随之增加。\n因此，为保证紧组合导航系统的可靠性及精确性，使系统在故障条件下仍能正常运行，需要检测并有效识别多卫星信息发生故障的情况，对故障及时进行排除、处理，并对系统进行重构。\n[0008] 目前应用于组合导航系统中卫星信息故障检测的算法普遍采用接收机自主完好性监测(RAIM)算法，该算法的基本原理是利用观测量冗余信息的一致性检验来测量观测量偏差情况，并找出故障卫星。常用RAIM算法是基于当前历元观测量的“快照”法，包括最小二乘残差法、奇偶空间法以及距离比较法、解的最大间隔法等，其中前三种方法本质上是等价的。以上算法能够实现在一定虚警率、误警率下系统发生故障的检测，但不能进行故障识别或只能识别单星故障。随着在轨导航卫星数的增多，单星故障识别已经无法满足导航需求，尤其在涉及生命安全等关键应用场合。\n[0009] 多故障识别算法包括n-2法、OWAS法、假设检验法、最优奇偶矢量法等进行观测量子集检验得到，但上述算法的实施对可见卫星数要求在6颗以上；并且以上算法都假设所有观测量一直具有相同的噪声统计特性，而实际上在发生故障以后，非故障观测信息也会受到污染而导致噪声水平加大。\n发明内容\n[0010] 本发明所要解决的技术问题是针对上述背景技术的不足，提供了惯性卫星紧组合导航系统的多故障检测装置及检测方法。\n[0011] 本发明为实现上述发明目的采用如下技术方案：\n[0012] 惯性卫星紧组合导航系统的多故障检测方法，包括如下步骤：\n[0013] 步骤1，确定惯性卫星紧组合导航系统的观测量矩阵；\n[0014] 步骤2，建立故障检测的检验统计量：构造奇偶空间矩阵并确定观测方程表达式，将观测误差矩阵投影到奇偶空间矩阵得到奇偶矢量，再根据显著性校验原理建立检验统计量；\n[0015] 步骤3，检测识别惯性卫星紧组合导航系统的故障，具体包括如下步骤：\n[0016] 步骤3-1，对n颗可见卫星进行子集划分，以每5颗可见卫星为一个子集，共有Cn5个子集；\n[0017] 步骤3-2，在每个子集中，利用步骤2所述的方法建立各个可见卫星的检验统计量，计算各个可见卫星在所属子集中的故障概率；\n[0018] 步骤3-3，叠加可见卫星在所属各子集中的故障概率，对故障概率和做归一化处理得到各个可见卫星的质量；\n[0019] 步骤3-4，判定质量小于规定值的可见卫星为故障卫星；\n[0020] 步骤4，剔除故障卫星，重构观测系统，确定观测方程：\n[0021] 步骤4-1，计算可见卫星数与故障卫星数的差得到可用卫星数：\n[0022] 若可用星数大于或者等于4，进入步骤4-2；\n[0023] 若可用星数小于4，进入步骤4-3；\n[0024] 步骤4-2，删去观测量矩阵中关于故障卫星的数据；\n[0025] 步骤4-3，删去观测量矩阵中质量值小于0.4的故障卫星的数据；\n[0026] 步骤4-4，加权处理修正后的观测量矩阵得到无故障卫星信息的观测量矩阵，卡尔曼滤波器根据无故障卫星信息的观测量矩阵确定观测方程，再将卫星状态信息的误差值反馈给惯性导航装置和GPS接收机。\n[0027] 所述惯性卫星紧组合导航系统的多故障检测方法中，步骤1的具体实施如下：\n[0028] 步骤1-1，惯导等效伪距测量单元根据惯性导航装置的位置速度信息、原始卫星的位置信息计算惯导与可见卫星的伪距；\n[0029] 步骤1-2，加法器用接收机输出的可见卫星伪距减去步骤1-1得到的惯导与伪距得到原始观测量矩阵；\n[0030] 步骤1-3，惯导辅助装置根据GPS输出的星历文件选择可见卫星之外的一颗卫星为辅助卫星；\n[0031] 步骤1-4，根据惯性导航装置测到的惯导位置信息、GPS接收机提供的原始卫星位置信息、GPS天线与惯导装置中心之间的校正距离，计算辅助伪距向量，将辅助伪距向量加入到步骤1-2得到的原始观测量矩阵；\n[0032] 步骤1-5，利用步骤1-4所得的观测量矩阵计算可见卫星、辅助卫星到GPS接收机的方向余弦确定观测量矩阵。\n[0033] 本发明采用上述技术方案，具有以下有益效果：能够对多颗卫星同时发生故障进行检测和识别，根据对卫星观测量的质量评价重构观测系统，提高了故障条件下紧组合导航系统的精度。\n附图说明\n[0034] 图1为惯性卫星紧组合导航系统多故障检测装置的示意图。\n[0035] 图2为检测识别惯性卫星紧组合导航系统的故障的示意图。\n[0036] 图3为故障识别方法得到的故障卫星质量值示意图。\n[0037] 图4为紧组合导航系统未进行观测量矩阵重构时的导航结果示意图。\n[0038] 图5为紧组合导航系统进行观测量矩阵重构时的导航结果示意图。\n具体实施方式\n[0039] 下面结合附图对发明的技术方案进行详细说明：\n[0040] 以紧组合导航系统包含双星故障为例阐述本发明的发明宗旨，惯性卫星紧组合导航系统的仿真采用动态航迹，中等精度惯性器件；GPS伪距残差为20m，不考虑卫星几何精度的影响。\n[0041] 如图1所示的惯性卫星紧组合系统多故障检测装置，包括观测量测量单元、惯导辅助单元、卫星故障评价单元、测量信息筛选单元、卡尔曼滤波器。观测量测量单元包括：惯导等效伪距测量单元和加法器。\n[0042] 惯导等效伪距测量单元的第一输入端与惯性导航装置的输出端连接，惯导等效伪距测量单元的第二输入端与GPS接收机的第一输出端连接，惯导等效伪距测量单元的输出端与加法器的第一输入端连接，GPS接收机的第二输出端与加法器的第二输出端连接，加法器的第一输出端与惯导辅助单元的输入端连接，加法器的第二输出端与测量信息筛选单元的第一输入端连接，惯导辅助单元的输出端与卫星故障评价单元的输入端连接，卫星故障评价单元的输出端与测量信息筛选单元的第二输入端连接，卡尔曼滤波器的第一输入端与测量信息筛选单元的输出端连接，卡尔曼滤波器的第二输入端与惯性导航装置的输出端连接，卡尔曼滤波器的输出端与惯性导航装置的反馈输入端、GPS接收机的反馈输入端连接。\n[0043] 观测量测量单元中：惯导等效伪距测量单元根据惯性导航装置测得的卫星位置信息、GPS接收机提供的原始卫星位置信息计算出可见星和辅助卫星的伪距；再用可见星和辅助卫星的伪距减去接收机提供的可见星伪距得到观测量矩阵。\n[0044] 卫星故障评价单元建立检验统计量，用显著性实验的统计学原理识别故障卫星，测量信息筛选单元剔除故障卫星，根据可用卫星个数修正观测量矩阵，卡尔曼滤波器根据修正后的观测量矢量确定观测系统中无故障卫星的位置信息并将位置信息反馈给惯性导航装置和GPS接收机。\n[0045] 惯性卫星紧组合导航系统的多故障检测方法，包括如下步骤：\n[0046] 步骤1，确定惯性卫星紧组合导航系统的观测方程：\n[0047] 步骤1-1，惯导等效伪距测量单元根据惯性导航装置的位置速度信息、原始卫星的位置信息计算惯导与可见卫星的伪距ρI；\n[0048] 步骤1-2，加法器根据如下公式（1）计算原始测量矢量ρ0：\n[0049] ρ0=ρG-ρI （1）\n[0050] 式中：ρG为接收机输出的原始可见卫星伪距，ρI为惯导与可见卫星伪距；\n[0051] 步骤1-3，惯导辅助装置根据GPS输出的星历文件根据公式（2）选择可见[0052] 卫星之外的一颗\n[0053] 仰角E较大的卫星作为辅助卫星：\n[0054] \n[0055] 式中：φ1为惯导位置经度，φ2为卫星的轨位经度，β为惯导位置纬度；\n[0056] 步骤1-4，根据惯性导航装置测到的惯导位置信息、GPS接收机提供的原始卫星位置信息、GPS天线与惯导装置中心的校正距离d，根据最小二乘法计算辅助伪距向量r，将辅助伪距向量r加入到步骤1-2得到观测量矩阵\n[0057] \n[0058] \n[0059] 式（3）中，(xINS,yINS,zINS)是惯导在ECEF坐标系下的位置信息，(xi,yi,zi)是可见卫星在ECEF坐标系下的位置信息，i=1,2…,n,n为可见卫星的数目；\n[0060] 步骤1-5，确定观测量矩阵\n[0061] 其中，\n[0062] ei1=(x-xi)/ρi，ei2=(y-yi)/ρi，ei3=(z-zi)/ρi，i=1,2…,n（4）；\n[0063] 步骤2，建立故障检测的检验统计量：\n[0064] 步骤2-1，构造奇偶空间矩阵，对观测量矩阵G进行QR分解：\n[0065] G=QR（5）\n[0066] 其中，Q为正交矩阵，R为上三角矩阵，带入观测方程有：\n[0067] QTy=Rx+QTε（6）\n[0068] \n[0069] 其中：y为伪距、伪距率观测量，x为系统状态量，ε为观测误差，矩阵P为奇偶空间矩阵，P的行向量为相互正交的单位向量，同时与G的列向量正交；\n[0070] 步骤2-2，将观测误差矩阵投影到奇偶空间矩阵得到奇偶矢量p，令：p=Py=Pε，奇偶矢量p具有与观测信息故障幅值线性相关的特性；\n[0071] 步骤2-3，根据显著性检验原理建立检验统计量公式（8）：所有观测量都无故障H0(零假设)，至少有一个观测量有故障H1(备择假设)：\n[0072] \n[0073] H0、H1的概率密度函数相同，矢量b表示故障幅度，为未知参数。\n[0074] 检验统计量计算如公式（9）所示：\nT T -1\n[0075] TR(p)=p(PCP) p （9）\n[0076] 式中，C为观测量矩阵的协方差阵。\n[0077] 步骤3，检测识别惯性卫星紧组合导航系统的故障，如图2所示具体包括如下步骤：\n[0078] 步骤3-1，对n颗可见卫星进行子集划分，以每5颗可见卫星为一个子集，共有Cn5个子集；\n[0079] 步骤3-2，在每个子集中，利用步骤2所述的方法计算各个可见卫星的检验统计量，计算假设检验p值，p为检验统计量在零假设下等于或超出给定观测值的概率，利用公式（10）计算各个可见卫星在所属子集中的故障概率：\n[0080] m(p)=F0(TR(p)≥TR(p0)) （10）\n[0081] 其中，F0是零假设(即系统无故障假设)下的概率；\n[0082] 步骤3-3，叠加各个可见卫星在所属子集中的故障概率m(p)，对故障概率和做归一化处理得到各个可见卫星的质量qi；\n[0083] 步骤3-4，判定质量小于0.6的可见卫星为故障卫星；\n[0084] 步骤4，剔除故障卫星，重构观测系统，重构观测方程：\n[0085] 步骤4-1，计算可见卫星数与故障卫星数的差得到可用卫星数：\n[0086] 若可用星数大于或者等于4，进入步骤4-2；\n[0087] 若可用星数小于4，进入步骤4-3；\n[0088] 步骤4-2，删去观测量矩阵中关于故障卫星的数据；\n[0089] 步骤4-3，删去观测量矩阵中质量值小于0.4的故障卫星的数据；\n[0090] 步骤4-4，加权处理修正后的观测量矩阵，表达式如公式（11）所示：\n[0091] \n[0092] 卡尔曼滤波器根据加权处理的观测量矩阵确定观测系统中各无故障卫星的位置信息，再将卫星状态信息的误差值反馈给惯性导航装置和GPS接收机。\n[0093] 在0-2000s内设置卫星故障仿真条件考虑如下因素：可见星数目为9颗，卫星1、卫星3、卫星5同时发生故障，故障幅值分别5σ、8σ、10σ，发生故障时间为950s，故障持续时间50s。仿真结果如图3、图4、图5所示。\n[0094] 图3为本发明提供故障识别方法得到的故障卫星质量值，由本图可以看出，在发生故障的时刻，卫星1、3、5质量值有明显下降，其余卫星质量值无明显下降，但由于故障信息的污染，其质量值也有所减少。\n[0095] 图4为上述故障条件下，紧组合导航系统未进行观测重构时的导航结果。由本图可以看出，导航结果受到严重影响，姿态、位置、速度误差都呈现明显增加。\n[0096] 图5为上述故障条件下，紧组合导航系统进行观测重构后的导航结果。由本图可以看出，在对故障信息排除、非故障信息采用质量值加权的观测重构方法后，导航结果明显提高，恢复到无故障条件下的导航精度水平。\n[0097] 综上所述，本发明所涉及的惯性卫星紧组合导航系统的多故障监测装置能够对多颗卫星同时发生故障进行检测和识别，多故障检测方法根据对卫星观测量的质量评价重构观测系统，提高了故障条件下紧组合导航系统的精度。