基于SSLLE和自适应阈值的导弹舵面故障检测方法

技术领域

本申请涉及导弹故障检测领域，尤其是一种基于SSLLE和自适应阈值的导弹舵面故障检测方法。

背景技术

随着军事斗争的不断升级，导弹成为了现代化战争中负责远程精确打击和战略威慑任务的主要武器之一，尤其是巡航导弹，由于具有机动性能好、命中精度高等优势，在局部战争中发挥的影响力正在不断提升。

然而在目前现代导弹自身结构日趋复杂、作战任务难度不断提高的主流趋势下，导弹由于外部各种干扰的作用，其在执行任务时发生故障的可能性也在直线上升。一旦导弹出现故障，不仅会直接影响导弹的飞行动力学特性，导致导弹飞行失稳，使得控制性能大幅下降，还会导致任务执行失败。因此及时有效的针对导弹系统进行状态监控并建立完备的故障诊断体系是确保导弹正常执行作战任务的关键，为保障导弹的可靠性与安全性提供可能。

发明内容

本申请人针对上述问题及技术需求，提出了一种基于SSLLE和自适应阈值的导弹舵面故障检测方法，本申请的技术方案如下：

一种基于SSLLE和自适应阈值的导弹舵面故障检测方法，该方法包括在导弹系统飞行过程中的任意第t检测时刻，t为参数且t≥1：

提取导弹系统在第t检测时刻对应的时间窗口内的飞行状态数据的数据波动特征得到高维数据集，覆盖时长为T的时间窗口包括第(t-T)检测时刻至第t检测时刻之间的若干个检测时刻；

利用SSLLE算法通过映射矩阵将高维数据集映射到低维空间得到低维数据集，映射矩阵基于导弹系统的正常稳态数据计算得到，正常稳态数据是正常工作的导弹系统的飞行状态数据；

基于低维数据集提取得到第(t-T)检测时刻对应的故障检测指标值；

基于漂移阈值方法确定第(t-T)检测时刻对应的自适应阈值；

当第(t-T)检测时刻对应的故障检测指标值超过对应的自适应阈值时，确定导弹舵面在第t检测时刻存在故障，否则确定导弹舵面在第t检测时刻正常工作。

其进一步的技术方案为，提取低维数据集的故障检测指标值的方法包括：

计算低维数据集与导弹系统的正常稳态数据的聚类中心的欧式距离，将计算得到的欧式距离作为第(t-T)检测时刻对应的故障检测指标值。

其进一步的技术方案为，确定第(t-T)检测时刻对应的自适应阈值：

确定第(t-T)检测时刻对应的自适应阈值

其进一步的技术方案为，该方法还包括：

确定

其中，

其进一步的技术方案为，该方法还包括：

对导弹系统的正常稳态数据提取数据波动特征得到高维数据集X＝[x

对高维数据集X进行标准化处理得到

计算标准化处理后的高维数据集

基于主成分分析法确定使得

其进一步的技术方案为，提取得到高维数据集的方法包括：

对第t检测时刻对应的时间窗口内的飞行状态数据提取得到表征数据波动特征的若干项统计信息，构建得到高维数据集，统计信息包括标准差、极差、范数和范数变化率。

其进一步的技术方案为，提取的导弹系统的飞行状态数据包括导弹系统的角加速度、加速度、姿态角、迎角和侧滑角。

其进一步的技术方案为，导弹舵面存在的故障的类型包括舵面卡死故障、舵面松浮故障和舵面缺损故障。

本申请的有益技术效果是：

本申请公开了一种基于SSLLE和自适应阈值的导弹舵面故障检测方法，该方法基于故障发生后，导弹系统稳定性被破坏，动力学与运动学发生突变的特点，利用滑动窗口和波动特征对导弹系统飞行状态数据的变化情况进行描述构建高维数据集，然后依据SSLLE算法将高维数据集映射到低维空间，剔除冗余特征，放大数据间的差异性，采用聚类定义故障检测指标值并基于DSPOT计算自适应阈值，以此判断导弹的故障状态，能够更加准确的描述导弹系统尤其是导弹舵面发生故障后的情况，从而可以准确实时的检测到导弹舵面故障，有利于故障的预防与解决，减小故障对导弹系统的危害后果。

附图说明

图1是本申请一个实施例中的导弹舵面故障检测方法的方法流程图。

图2是本申请一个实例中的导弹系统的结构示意图。

图3是一个仿真实例中，在导弹系统正常工作时采集到的迎角的波形图。

图4是图3所示的仿真实例中，在导弹系统正常工作时采集到的导弹系统的飞行弹道图。

图5是一个仿真实例中，在80s处注入舵面松浮故障后时，整个仿真时长内采集到的迎角的波形图。

图6是图5所示的仿真实例中，在80s处注入舵面松浮故障后时，整个仿真时长内采集到的导弹系统的飞行弹道图。

图7是图5和6的仿真实例中，利用本申请的方法计算得到的各个检测时刻的故障检测指标值的曲线图和自适应阈值。

图8是图7的局部放大示意图。

具体实施方式

下面结合附图对本申请的具体实施方式做进一步说明。

本申请公开了一种基于SSLLE和自适应阈值的导弹舵面故障检测方法，该方法包括如下步骤，请参考图1所示的流程图：

在导弹系统飞行过程持续执行如下过程进行导弹舵面故障检测，对于任意第t检测时刻，t为参数且t≥1，执行的方法包括：

步骤S1，提取导弹系统在第t检测时刻对应的时间窗口内的飞行状态数据的数据波动特征得到高维数据集，覆盖时长为T的时间窗口包括第(t-T)检测时刻至第t检测时刻之间的若干个检测时刻。

在一个实施例中，在每个检测时刻提取的导弹系统的飞行状态数据包括导弹系统的角加速度、加速度、姿态角、迎角和侧滑角。

在采集到第t检测时刻对应的时间窗口内的飞行状态数据后，对时间窗口内的飞行状态数据提取得到表征数据波动特征的若干项统计信息从而构建得到高维数据集X＝[x

步骤S2，利用SSLLE算法通过映射矩阵将高维数据集X＝[x

SSLLE算法中本征维数d和映射矩阵的选取决定了降维的效果，本申请中，映射矩阵基于导弹系统的正常稳态数据计算得到，正常稳态数据是正常工作的导弹系统的飞行状态数据。在一个实施例中，确定本征维数d和映射矩阵的方法包括：

获取正常工作的导弹系统在个不同的时刻的飞行状态数据构建得到高维数据集，与步骤S1获取实际在线飞行状态数据的方法类似，因此同样将高维数据集表示为X＝[x

由此可以通过映射矩阵A＝YX

步骤S3，基于低维数据集提取得到第(t-T)检测时刻对应的故障检测指标值。

通过KMEANS可以聚类得到导弹系统的正常稳态数据的聚类中心，然后计算在线获取到的低维数据集与导弹系统的正常稳态数据的聚类中心的欧式距离，将计算得到的欧式距离作为第(t-T)检测时刻对应的故障检测指标值。

步骤S4，基于漂移阈值方法确定第(t-T)检测时刻对应的自适应阈值。

确定第(t-T)检测时刻对应的自适应阈值

和σ

其中，

步骤S5，当第(t-T)检测时刻对应的故障检测指标值超过对应的自适应阈值时，确定导弹舵面在第t检测时刻存在故障，否则确定导弹舵面在第t检测时刻正常工作。也即对于任意第t个检测时刻有

在一个实施例中，利用本申请的方法可以检测的导弹舵面存在的故障的类型包括舵面卡死故障、舵面松浮故障和舵面缺损故障。

在一个实施例中，根据导弹系统的模型可以确定舵面分配矩阵，从而确定在各种类型的故障下的表现形式，舵面分配矩阵表征导弹舵面的舵偏角与舵控信号之间的关系。

比如在一个实例中，导弹系统的模型选取BGM-109战斧巡航导弹作为对象，主要包含弹体动力学&运动学、发动机、执行器(舵面)、大气环境系统、控制系统、导航系统以及制导系统等模块，导弹系统运动方程组由动力学、运动学、质量关系、几何关系与控制关系方程组构成导弹系统模型。导弹系统舵面采用十字型气动外形布局如图2所示，由四个舵机共同完成对于导弹滚转、俯仰和偏航的控制。基于此模型，导弹系统的舵面分配矩阵可以表示为：

其中，δ

(1)舵面卡死故障，也即导弹系统的某一舵面偏转功能失灵，无法继续偏转到目标角度而维持在一固定位置，可以表示为

(2)舵面松浮故障，也即导弹系统的某一舵面无法正常响应舵控信号，无法提供任何力矩而处于松浮状态，可以表示为

(3)舵面缺损故障，也即导弹系统的某一舵面的增益无法达到预期效果，导致执行效能有所下降，可以表示为

在上述三种情况中，t

重复上述过程，在每个检测时刻对导弹舵面进行故障检测，直到因为完成检测时长或者导弹系统的飞行过程结束而到达最后一个检测时刻。

为了说明本申请的方法的有效性，进行如下仿真试验，设置导弹系统的初始飞行速度为240m/s，初始发射位置为[0,0,800]，导航点设置为[5000,0,1000]、[10000,0,1000]、[15000,0,1000]、[20000,0,1000]、[25000,0,1000]，单位均为米，如图4和图6中的圆圈表示。最大仿真时长为120s。

在导弹系统正常工作的状态下，在120s的仿真时长内的各个采样时刻对导弹系统采集到的各类飞行状态数据，其中采集到的迎角的波形图如图3所示，其他各类飞行状态数据不再图示。导弹系统的飞行弹道图如图4所示。

在一个仿真实例中，在t＝80s时，对导弹系统的2号舵面注入舵面松浮故障，持续时间到仿真结束。则在该仿真实例中，在120s的仿真时长内的各个采样时刻对导弹系统采集到的迎角的波形图如图5所示、导弹系统的飞行弹道图如图6所示。

对比图3和图5可以看出，当导弹舵面在80s处出现了舵面松浮故障时，采样到的迎角明显偏离了稳态值，导致导弹系统数据处于异常波动状态。同样的，采样到的导弹系统的加速度、角加速度、姿态角和和侧滑角也都明显偏离了稳态值。对比图4和图6可以看出，当导弹舵面在80s处出现了舵面松浮故障时，导弹系统迅速下坠并可能导致坠毁等严重后果。

基于本申请的方法进行故障检测，使用的时间窗口覆盖的时长T＝3s，各个不同的检测时刻的故障检测指标值的曲线图和自适应阈值的曲线图如图7所示，图7中的局部示意图如图8所示，当导弹系统正常工作、采集到的飞行状态数据在正常范围内波动时，各个检测时刻的自适应阈值会及时更新以保证故障检测准确性。在77.05s之前的仿真时间内，各个检测时刻的故障检测指标值都小于对应的自适应阈值，在77.05s及之后的仿真时间内，各个检测时刻的故障检测指标值都超过对应的自适应阈值，则会在80.05s处检测到dis(t-3)≥th(t-3)，从而检测到导弹舵面存在故障，检测耗时很短，可以及时发现舵面故障。

而且该方法对于各种类型的舵面故障都可以快速准确的检测，在一个仿真实例中，对各种不同类型的舵面故障的故障检测数据如下表所示：

以上所述的仅是本申请的优选实施方式，本申请不限于以上实施例。可以理解，本领域技术人员在不脱离本申请的精神和构思的前提下直接导出或联想到的其他改进和变化，均应认为包含在本申请的保护范围之内。