基于证据推理的激光惯组最优测试时机确定方法与流程

本发明属于大型工业结构控制系统的惯性器件健康管理与维护领域，基于证据推理有效融合专家知识与监测数据，确定复杂条件下的激光惯组最优测试时机。
背景技术：
：激光惯组作为大型工业结构控制系统的重要导航设备，在火箭飞行过程中建立空间运动的方位基准，实时测得角速度和运动的线加速度等飞行参数，送入弹载计算机转换为模拟量。这些模拟量经放大综合后，作为伺服机构的输入，转换为相应的指令信号驱动喷管摆动，从而控制和稳定火箭的飞行姿态。作为大型工业结构控制系统的核心部件，激光惯组的性能在很大程度上决定了火箭导航系统的可靠性与精度。激光惯组主要包括激光陀螺和石英加速度计，这两大器件也是惯性导航的主要误差源。陀螺仪和加速度计的故障率与精度很大程度上影响了惯性导航系统的可靠性。以激光陀螺为例，角速度增量信号通过a/d转换成为数字脉冲输出，而由于脉冲计数器只能计数整数个脉冲信号，因此在一个采样周期内，整数部分脉冲增量输出，而小数部分则将累计，由此产生量测误差，影响导航精度。通常而言，激光惯组制造工艺精细、造价昂贵，在实际对其进行测试的持续时间短，获取的故障数据较少，且存在各种不确定性。因此，如何有效处理高价值的小样本，解决不确定性带来的影响，是当前面临的突出问题之一。对于纯惯性制导的火箭而言，其惯性器件的误差随时间不断积累，使得导航误差接近70％。目前，在实际应用中，对激光惯组进行测试所获得的是其原始累计脉冲量，无法直接反映出惯性器件的故障状态，使得操作人员难以实时分析和掌握激光惯组的性能。因此，如何有效结合专家知识建立惯组的指标体系，进行融合，分析脉冲量和静态导航精度之间的关系，建立精确的惯组导航模型，以对惯性器件进行最优测试，是当前面临突出问题之二。技术实现要素：本发明的目的是：针对激光惯组测量得到的原始累计脉冲量无法反映惯性器件的故障状态，利用证据推理(evidentialreasoning，简称er)有效融合激光惯组的特征指标，嵌入专家知识，建立起惯组的原始累计脉冲量与静态导航误差之间的非线性解析模型。在er建模过程中，考虑到初始参数(指标权重和参考值)对输出的影响较大这一实际，采用优化算法对初始参数进行优化，得到更为精确的er模型；在建立er模型的基础上，利用维纳过程对未来时刻的导航误差进行预测，为惯组后期的最优测试与健康维护提供借鉴和依据。本发明的技术方案在于：基于证据推理的激光惯组最优测试时机确定方法，该方法如下：融合激光惯组的特征指标与专家知识，建立起以脉冲增量为输入、静态导航误差为输出的指标体系；再根据er模型建立输入与输出关系的非线性模型；基于非线性模型，利用维纳过程进行预测，从而得到激光惯组的最优测试时机。基于证据推理的激光惯组最优测试时机确定方法，该方法如下：(1)脉冲数据预处理：对原始获取的激光惯组的原始累计脉冲量进行差分得到脉冲增量；(2)蒙特卡罗仿真：利用蒙特卡罗算法对脉冲增量进行仿真，并引入随机噪声，代入导航方程，获得对应的静态导航误差为输出；(3)基于er模型的导航模型构建：以脉冲增量为输入，并利用er模型进行多源信息融合，以静态导航误差为输出，建立输入与输出关系的非线性模型；(4)基于维纳过程的最优测试时机确定：利用维纳过程进行预测，结合专家知识确定导航误差阈值，得到激光惯组的最优测试时机。所述脉冲数据预处理过程还包括剔除差分得到的脉冲增量的异常值；利用蒙特卡罗算法对剔除了异常值的脉冲增量进行仿真。利用蒙特卡罗算法对剔除了异常值的脉冲增量转化为模拟量，引入随机噪声，代入导航方程，获得对应的静态导航误差为输出。所述基于er模型的导航模型构建的具体过程为：以陀螺仪三个轴方向的累计脉冲量分别为xw,yw,zw，加速度计三个轴方向的累计脉冲量xa,ya,za，中间层指标陀螺仪为w，加速度计为a构建指标体系；利用er模型，将w与a进行融合，最终得到输出导航误差y，并进行参数优化。本发明的技术效果在于：针对激光惯组维护存在的主观性与不确定性等问题，本发明将专家的知识与监测数据有效融合，推理过程透明，结果具有说服力，为惯组后期的最优测试与健康维护提供准确依据。附图说明图1为本发明的指标体系。图2为本发明的具体实施步骤。图3为本发明的模型测试结果。图4为本发明选择的最优测试时机。具体实施方式融合专家知识与监测数据的激光惯组最优测试时机确定方法，该过程如下。一、脉冲数据预处理：对原始获取的激光惯组的原始累计脉冲量进行差分得到脉冲增量；对原始获取的激光惯组的测试数据为加速度计和陀螺仪的原始累计脉冲量，而实际导航中所需的脉冲输入形式为脉冲增量，因此需要对原始累计脉冲量进行差分得到脉冲增量；对脉冲增量进行统计后，若某一时刻脉冲增量明显高于某个阈值时，判断其为异常值。则需对该异常值进行剔除；其具体过程为：设异常值个数为n，总脉冲数为m，则异常率可计算为：k＝m/n(1)将剔除异常值后得到的脉冲增量作为步骤(2)中导航方程的输入。二、蒙特卡罗仿真：导航方程用于在理想情况下实现激光惯组脉冲增量与静态导航误差之间的转换，而激光惯组脉冲增量通常带有量化误差、机械安装误差以及外界干扰等。因此，导航方程无法真正反映出激光惯组的脉冲增量与静态导航误差之间的复杂动态关系，需要对脉冲增量进行蒙特卡罗仿真，模拟得到激光惯组在不同状态下的静态导航误差；由于脉冲增量为数字量，对其进行蒙特卡罗仿真不具有实际意义，因此需将剔除异常值后得到的脉冲增量转化为对应的模拟量。并引入随机噪声，代入导航方程，该代入导航方程的过程如下：其中，前60000组模拟量用于瞄准，而后每60000组模拟量用于静态导航，取每次静态导航最后时刻对应的静态导航误差作为该组模拟量的输出，共计获得80组输出。三、基于er模型的导航模型构建：(1)构建指标体系：在er模型中，输入指标的选择不仅要符合通用性和完备性原则，还要考虑到激光惯组实际情况，从实际入手，紧贴客观现实。为了有效建立起激光惯组的指标体系，根据相关专家知识我们知道激光惯组是以脉冲量作为输出，通过脉冲量能够反映出载体的速度与角加速度。因此对于激光惯组来说，脉冲量是其最直接同时也是最明显的指标。因此，在方法的最开始，选择激光惯组的原始累计脉冲量作为基础指标。构建指标体系过程中，以去除异常值后的脉冲增量每间隔60000组累加一次作为底层输入指标；即其中，xi表示激光惯组的脉冲增量，xi为累计脉冲量，作为底层指标，加速度计与陀螺仪构成中间层指标，如附图1所示。根据上述构建原则，选取部分输入xi和输出yi，如下表：表1部分输入与输出(2)er模型：在er模型中，当输入的指标确定以后，首先定义这一系列指标为：e＝(e1,e2,…,el)(3)在融合过程中，每个指标会被赋予不同的权重，假设各个指标的权重为w＝(w1,…,wi,…,wl)，归一化得到：定义每个指标的n个参考值为将各个指标的输入利用基于效用的转化规则转化为信度分布：将所有的指标输入的数值都描述为信度分布形式时，根据yang等提出的解析er算法，融合所有指标，得到最后模型的输出结果：在实际应用中，需要的通常是数值型输出，应此需要把输出的信度分布转化为数值，设输出结果的参考值为d＝{d1,…d2,…dn}，则此时的结果输出数值为若即结果不存在无知性，则式(9)可进一步简化为记陀螺仪三个轴方向的累计脉冲量分别为xw,yw,zw，加速度计三个轴方向的累计脉冲量分别为xa,ya,za，记中间层指标陀螺仪为w，加速度计为a，对w与a进行融合，最终得到输出导航误差y。(3)er模型的参数优化：为了准确建立累计脉冲量与导航误差的关系，需要对参数进行优化，其中待优化参数包括6个底层指标权重及对应的参考值(每个指标设立3个参考值)，构成待优化参数向量：p＝[w1,w2,w3,w4,w5,w6,w7,w8,f11,f12,f13,f21,f22,f23,f31,f32,f33,f41,f42,f43,f51,f52,f53,f61,f62,f63]其中，w1,w2,w3,w4,w5,w6为底层指标的权重，w7,w8为中间层指标的权重，f为输入指标的参考值。在得到6个底层指标的输入之后，需要对输入设立参考值及初始权重，输出的导航误差参考值定为0，10，20，不参与优化。根据之前所述的参考值取值范围及权重约束条件，得到参考值及权重的初始给定值如下：表2初始参数设er模型输出的静态导航误差为由导航方程得到的静态导航误差为y，基于均方差最小原则，构建以下优化模型：式(11)实质上是带约束的目标函数求最小值问题，利用matlab进行优化求解。选择前50组数据进行训练，由于待优化参数向量维数较低，选用matlab中fmincon优化工具箱进行优化，最终得到优化后的参数向量取值如下：表3优化参数选取全部数据进行测试，测试效果见附图4，为进一步测试所建立模型的有效性，选取多个误差指标进行检验，所选误差指标包括均方误差mse、平均绝对百分误差mape、平均绝对误差mad这三个值，下表为计算结果：表4误差分析误差rmse(米)mapemad(米)优化后er模型1.55550.18090.9425优化的er模型9.28901.89492.4280。四、基于维纳过程的最优测试时机确定：基于步骤三中得到的静态导航误差输出，利用维纳过程进行预测，结合专家知识确定导航误差阈值，从而确定激光惯组的最优测试时机。(1)维纳过程：设随机过程{x(t)，t≥0}的状态空间是e＝(-∞，+∞)，若满足：a.x(t)是平稳独立增量过程；b.x(s+t)-x(s)～n(0，c2t)，即任一增量服从n(0，c2t)的正态分布；则称x(t)是维纳过程。(2)参数估计：标准维纳过程的表达形式如下：y＝a+λt+σβ(t)(12)其中，a为性能退化初值；λ为一次项系数，σ为漂移系数，β(t)为标准brown运动。由于存在上述未知参数，需要讲行参数估计，采用极大似然法，估计结果如下：其中，{(yi，ti)，i＝1，2，，3，...n}表示观测值与对应的观测时刻。得到α＝-583.33，λ＝7.3839，σ＝1.1878；即y＝-583.33+7.3839t+1.1878β(t)，t＝1，2，.....80(16)。(3)激光惯组最优测试时机确定：为了确定激光惯组的最优测试时机，通过专家知识可知，静态导航误差可以反映出激光惯组的故障状态。当静态导航误差达到不同的阈值时，可以反映出激光惯组的性能状态。因此利用维纳过程对输出的导航误差进行预测，设置不同的阈值y1,y2，当导航误差达到不同阈值时给出不同的测试频率f1,f2，从而确定惯组的最优测试时机。取t＝81，......90，预测结果如下：表5预测结果根据专家知识与实际情况，设立阈值，此时最优测试时机选择见附图4。当前第1页12