文章快速检索  
  高级检索
基于深度森林的卫星ACS执行机构与传感器故障识别
程月华1, 江文建1, 杨浩1, 薛琪1, 廖鹤2     
1. 南京航空航天大学 自动化学院, 南京 211100;
2. 南京航空航天大学 航天学院, 南京 210016
摘要: 针对卫星姿态控制系统(ACS)闭环回路的故障难以辨识的问题,引入深度森林算法,实现执行机构与传感器故障识别。首先针对可获取的少量卫星姿态控制系统遥测数据,结合系统动力学特性,研究合适的特征选择和特征提取方法,再结合深度森林算法进行故障信息学习与辨识,建立故障预测模型,实现执行机构故障与传感器故障的识别。半物理仿真结果表明:在存在气浮台干扰力矩、卫星转动惯量未知、飞轮非线性特性、闭环故障传播等多种不利因素情况下,深度森林算法对于执行机构和传感器故障具有高效的识别能力。
关键词: 深度森林算法    卫星姿态控制系统    执行机构    传感器    故障识别    
Fault identification of actuators and sensors of satellite attitude control systems based on deep forest algorithm
CHENG Yuehua1, JIANG Wenjian1, YANG Hao1, XUE Qi1, LIAO He2     
1. School of Automation Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 211100, China;
2. School of Astronautics, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
Abstract: The difficulty in identifying the faulted sensors and actuators of satellite Attitude Control Systems (ACS) resides in the spreading of faults in the closed loop. The deep forest algorithm is introduced in this study to build a fault prediction model to achieve the isolation of the sensor faults and actuator faults. After collecting healthy ACS telemetry data to group a training set according to the dynamic characteristics of ACS, we apply appropriate feature selection and extraction methods to the training set, obtaining the features of both the sensor faults and actuator faults. The deep forest algorithm, with its strong generalization ability, is then used to learn and identify fault information, thereby establishing a fault prediction model to realize the recognition of actuator and sensor faults. The results of the semi-physical simulation indicate that the proposed method can identify the faults of sensors and actuators effectively in the presence of many uncertain factors, such as the interference moments of air bearing testbeds, unknown moments of the inertia of satellites, nonlinear characteristics of flywheels and fault propagation in the closed loop of ACS.
Keywords: deep forest algorithm    satellite attitude control systems    actuators    sensors    fault identification    

卫星姿态控制系统(Attitude Control System, ACS)执行机构或传感器发生故障,会引发卫星ACS姿态失稳、丢失姿态,从而影响卫星平台的安全稳定运行[1-2]。开展卫星ACS执行机构与传感器故障识别研究,可有效地监测卫星ACS的运行状态,在故障发生时及时识别故障类型,采取措施隔离故障,对提高卫星在轨可靠性和安全性具有重要意义[3]

近年来,国内外学者针对航天器执行机构和传感器的故障诊断问题,开展了诸多的基于模型和基于信号处理的研究工作,并取得了丰硕的研究成果[4-5]。由于卫星ACS组成部件多、空间环境复杂,难以建立精确的数学模型,故基于模型的故障诊断方法应用于卫星ACS具有一定的局限性。近年来,随着计算机与人工智能技术的不断发展以及大量机器学习算法的相继问世,基于数据驱动的智能故障诊断方法被广泛地应用于卫星ACS的故障诊断[6-8]。基于数据驱动的智能故障诊断方法主要通过运用典型的机器学习算法,如人工神经网络(Artificial Neural Network, ANN)[9-10]、支持向量机(Support Vector Machine, SVM)[11-12]和随机森林(Random Forest, RF)[13]等,对卫星执行机构系统的高维采样数据进行表征学习,发现故障数据特征与故障模式之间的联系,从而实现卫星执行机构系统的故障诊断。其中ANN由于其高度的并行处理、联想记忆、自学习以及极强的非线性映射能力而被广泛应用。如文献[9]提出一种基于ANN观测器的陀螺故障诊断方法,利用ACS内的冗余关系设计故障诊断逻辑,同时利用先验模型知识和ANN的非线性建模特性对陀螺故障进行估计,实现了对卫星ACS陀螺的故障诊断。文献[10]利用动态ANN模拟反作用飞轮的运行特性,实现了对卫星ACS的故障诊断。但是,ANN模型的准确性,通常依赖大量的数据作为训练样本。而卫星发生故障的概率低,故障样本少,利用ANN对卫星ACS建立故障模型,其模型准确性势必会受到一定的影响。同时ANN也存在训练时间长、参数调整复杂、易陷入局部最优等问题。针对卫星故障样本少的问题,文献[11]利用SVM对提取到的少量故障特征样本进行分类,实现了卫星ACS执行机构的故障检测与识别。但是,SVM是一种典型的二分类算法,在解决卫星ACS多种故障模式分类问题时存在一定的困难。同时SVM也存在对缺失数据敏感,对参数和核函数的选择敏感等问题。

近年来,随着深度学习(Deep Learning)的概念被提出,深度神经网络(Deep Neural Network, DNN)由于其强大的表征学习能力和拟合能力,成为基于数据驱动的智能故障诊断领域的研究热点之一[14]。但DNN同样存在对训练样本要求高、对超参数敏感和训练时间长等问题。而卫星的故障样本较少,星载计算机性能不高且实时性要求高,因此DNN尚未应用于卫星故障诊断。文献[15-16]在DNN基础上提出了一种基于决策树的集成分类算法,称作深度森林(Deep Forest)算法或多粒度级联森林算法(multi-grained cascade forest, gcForest)。该算法不仅具有强大的表征学习能力,同时也对训练样本量要求不高,训练速度快,超参数少且对超参数的选择不敏感等优势。文献[17]将深度波尔兹曼机(Deep Boltzmann Machine, DBM)与深度森林相结合应用于工业故障诊断。本文针对存在环境干扰力矩、卫星转动惯量未知、飞轮非线性特性、闭环故障传播等诸多不利因素的卫星姿态控制系统发生传感器和执行机构故障情形,利用可获取的少量遥测数据信息,研究合适的特征选择和特征提取方法;利用深度森林强大的泛化能力对故障信息进行学习与辨识,建立故障预测模型,实现传感器与执行机构的故障辨识。

1 问题的提出

卫星ACS结构如图 1所示,由控制器、执行机构、卫星本体以及传感器组成。图中:qd为卫星期望姿态;ωiqi分别为卫星当前角速度和姿态角;ωo为陀螺仪实际输出;qo为姿态敏感器实际输出值;Ms为执行机构实际输出力矩;Md为卫星所受空间环境干扰; fa为执行机构故障模型; fs为传感器故障模型。

图 1 卫星ACS结构框图 Fig. 1 Structure diagram of satellite attitude control systems

定义卫星ACS姿态运动学描述为

$ \mathit{\boldsymbol{\dot q = }}\frac{1}{2}\mathit{\boldsymbol{q}} \odot \mathit{\boldsymbol{\omega = }}\frac{1}{2}\mathit{\boldsymbol{ \boldsymbol{\varXi} }}(\mathit{\boldsymbol{q}})\mathit{\boldsymbol{\omega }} $ (1)
 

式中:ω为角速度;q=[q1:3Tq4]T=[q1q2q3q4]T为航天器姿态四元数,

$ [\mathit{\boldsymbol{q}} \odot ]\mathit{\boldsymbol{ = }}\left[ {\begin{array}{*{20}{c}} {{\mathit{\boldsymbol{q}}_4}{\mathit{\boldsymbol{I}}_3} + [{\mathit{\boldsymbol{q}}_{1:3}} \times ]}&{{\mathit{\boldsymbol{q}}_{1:3}}}\\ { - \mathit{\boldsymbol{q}}_{1:3}^{\rm{T}}}&{{\mathit{\boldsymbol{q}}_4}} \end{array}} \right] = [\mathit{\boldsymbol{ \boldsymbol{\varXi} }}(\mathit{\boldsymbol{q}})\;\;\mathit{\boldsymbol{q}}] $ (2)
 
$ \mathit{\boldsymbol{ \boldsymbol{\varXi} }}(\mathit{\boldsymbol{q}}) = \left[ {\begin{array}{*{20}{c}} {{\mathit{\boldsymbol{q}}_4}{\mathit{\boldsymbol{I}}_3} + [{\mathit{\boldsymbol{q}}_{1:3}} \times ]}\\ { - \mathit{\boldsymbol{q}}_{1:3}^{\rm{T}}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {{\mathit{\boldsymbol{q}}_4}}&{ - {\mathit{\boldsymbol{q}}_3}}&{{\mathit{\boldsymbol{q}}_2}}\\ {{\mathit{\boldsymbol{q}}_3}}&{{\mathit{\boldsymbol{q}}_4}}&{ - {\mathit{\boldsymbol{q}}_1}}\\ { - {\mathit{\boldsymbol{q}}_2}}&{{\mathit{\boldsymbol{q}}_1}}&{{\mathit{\boldsymbol{q}}_4}}\\ { - {\mathit{\boldsymbol{q}}_1}}&{ - {\mathit{\boldsymbol{q}}_2}}&{ - {\mathit{\boldsymbol{q}}_3}} \end{array}} \right] $ (3)
 

本文ACS的执行机构为动量轮,得到的卫星ACS的姿态动力学方程为

$ {\mathit{\boldsymbol{J}}_{\rm{s}}}\mathit{\boldsymbol{\dot \omega + \omega }} \times ({\mathit{\boldsymbol{J}}_{\rm{s}}}\mathit{\boldsymbol{\omega + }}{\mathit{\boldsymbol{H}}_{\rm{w}}}) = {\mathit{\boldsymbol{M}}_{\rm{c}}} + {\mathit{\boldsymbol{M}}_{\rm{d}}} $ (4)
 

式中:Js为卫星转动惯量矩阵;Hw为动量轮角动量;卫星ACS所受总力矩Mq=Mc+Md

1) 卫星运行空间环境复杂,难以建立精确的干扰力矩Md模型。干扰动量轮对卫星ACS的控制力矩,从而影响系统输出。

2) 反作用飞轮作为典型的卫星姿态控制部件,已经成为控制系统不可缺少的配置。但由于飞轮在进行姿态调整过程中,非线性摩擦、电机噪声以及轴承噪声等多种未知因素会带来控制中的不确定干扰力矩,难以精确Mδ的模型。

3) 执行机构发生故障时,其表现主要为执行机构的实际输出力矩会产生偏差Mf,通过动力学与运动学方程,影响传感器测量输出ωoqo。而传感器发生故障时,会通过叠加故障值fs,通过ACS闭环影响控制器的指令输出力矩Mc。ACS中执行机构与传感器故障在系统中传播和相互影响,给故障诊断和辨识带来了困难和挑战。

综上所述,由于卫星空间环境干扰Md和飞轮的非线性以及噪声Mδ的存在,卫星ACS难以用准确的模型描述,其传感器故障模型fs和执行机构故障模型fa识别困难。基于数据驱动的智能模式识别方法的发展,为此提供了一条有效途径。通过采集卫星ACS闭环系统的执行机构与传感器运行数据、以及控制系统中的指令信息,利用智能方法的泛化和学习能力,可实现ACS闭环系统的故障识别问题。

2 卫星ACS执行机构和传感器故障识别策略

深度森林算法是在深度学习理论以及深度神经网络的启发下,以随机森林(Random Forest,RF)算法为基础的一种有监督机器集成学习算法[18-20]。基于深度森林算法的卫星ACS执行机构和传感器故障识别策略如图 2所示。首先,采集卫星ACS测试数据构建数据特征集;然后结合卫星ACS姿态动力学和运动学的分析,对数据特征集进行特征选择以及特征提取;最后通过深度森林算法建立的卫星ACS正常和故障模型对输入特征预测分类结果,实现卫星ACS执行机构与传感器的故障识别。

图 2 卫星ACS故障识别策略 Fig. 2 Fault identification strategy of satellite attitude control systems
3 特征提取与识别方法 3.1 构建原始特征集

记控制指令力矩为Mc、动量轮指令转速增量为Δnc、动量轮实际转速为ns,当前控制周期下动量轮转速的期望增量为Δnc

ACS飞轮采用无刷直流电机控制,其电机的力矩平衡方程为

$ {\mathit{\boldsymbol{M}}_{\rm{s}}} = {\mathit{\boldsymbol{J}}_{\rm{w}}}\frac{{{\rm{d}}{\mathit{\boldsymbol{\omega }}_{\rm{w}}}}}{{{\rm{d}}t}} $ (5)
 

式中:Ms为动量轮实际输出力矩;ωw为动量轮角速度(rad/s);Jw为动量轮的转动惯量。记

$ {\mathit{\boldsymbol{M}}_{\rm{s}}} = \mathit{\boldsymbol{K}}\Delta {\mathit{\boldsymbol{n}}_{\rm{s}}} $ (6)
 

其中:Δns为动量轮实际转速增量;$ \mathit{\boldsymbol{K = }}\frac{{\mathit{\boldsymbol{J}}\pi }}{{30\Delta T}}$,为常值; ΔT为ACS的控制周期。由式(6)可得Ms与Δns成线性关系,因此动量轮指令控制力矩Mc与动量轮指令转速增量Δnc也成线性关系,即

$ {\mathit{\boldsymbol{M}}_{\rm{c}}} = \mathit{\boldsymbol{K}}\Delta {\mathit{\boldsymbol{n}}_{\rm{c}}} $ (7)
 

由于动量轮指令转速增量Δnc反映了飞轮的力矩输出情况,所以将其选择为系统状态空间的特征。

考虑到动量轮存在非线性摩擦力矩以及电机噪声力矩、轴承噪声力矩等不确定干扰力矩,令这些干扰力矩之和为Mδ,则动量轮指令控制力矩Mc与实际控制力矩Ms之间的关系可以表示为

$ {\mathit{\boldsymbol{M}}_{\rm{s}}} = {\mathit{\boldsymbol{M}}_{\rm{c}}} + {\mathit{\boldsymbol{M}}_{\rm{ \mathsf{ δ} }}} + {\mathit{\boldsymbol{f}}_{\rm{a}}} $ (8)
 

式中:fa为执行机构故障产生的偏差力矩,当执行机构无故障时,其值为0。

令当前控制周期动量轮实际转速增量Δns与指令转速增量Δnc的差值为Δrn,其表达式为Δrnns-Δnc,结合式(6)~式(8)可得

$ \Delta {\mathit{\boldsymbol{r}}_n} = \frac{1}{\mathit{\boldsymbol{K}}}({\mathit{\boldsymbol{M}}_{\rm{ \mathsf{ δ} }}} + {\mathit{\boldsymbol{f}}_{\rm{a}}}) $ (9)
 

当执行机构无故障时,$ \Delta {\mathit{\boldsymbol{r}}_n} = \frac{1}{\mathit{\boldsymbol{K}}}{\mathit{\boldsymbol{M}}_\delta }$。而当执行机构发生故障时,fa会在不同故障模式下呈现出不同的变化特点,故Δrn可作为反映状态模式空间的原始特征之一。

令Δθ为当前控制周期气浮台姿态角θ的增量,Δωz为当前控制周期气浮台沿z轴旋转的角速度ωz的增量。由上述分析可知,执行机构故障和传感器故障会使Δrn、Δnc、Δθ和Δωz发生变化,且这些变量在不同的故障模式下会有不同的变化特点。故本文选取Φ=[Δrn, Δnc, Δθ, Δωz]作为反映状态模式空间的特征集。

以上构建的ACS无故障时的数据原始特征集记为Φh,在故障模式f1f2下的数据特征集记为Φf1Φf2,其特征点的空间分布如图 3所示。

图 3 Φf1Φf2Φh的原始特征点空间分布 Fig. 3 Spatial distribution of original feature points of Φf1, Φf2 and Φh
3.2 故障特征提取

特征提取是指提取数据的主要特征,除去数据中包含的冗余信息和错误信息,使得数据的特征更加明显。因为卫星遥测数据的高维性,以及传输过程的损耗,为了进行故障识别有必要进行特征提取。

图 3可以看出,ΦhΦf2的分布较为集中,Φf1分布较为分散,但三者相互之间都有较大的交集,如果直接将原始特征数据集作为输入,训练深度森林以建立故障模式预测分类模型,将影响其分类精度。因此数据特征集Φ进行故障特征提取。

以正常和故障时的Δrn为例来分析故障对特征集Φ产生的影响,如图 4所示, t为时间。对比图 4(a)图 4(b)可以看出,ACS故障情况下Δrn的值整体上升,奇异性成分增多且波动程度增大。图 5给出了不同故障类型下Δrn曲线图。

图 4 Δrn的曲线图 Fig. 4 Graph of Δrn
图 5 不同故障情况下的Δrn曲线图 Fig. 5 Graphs of Δrn under different fault conditions

特征选择方法,不仅要反映出信号的短期的突发特性与变化幅值,还要能反映出信号长期的趋势和统计特性,从而准确区分正常情况与故障情况、以及不同故障模型。

通过分析ACS原始特征集的数据特性,本文经过特征选择后建立的特征集为Φ。其中,均值μ、方差σ2能很好地体现信号的统计学特性,而信息熵H和信号积分I能分别体现信号的变化剧烈程度和变化趋势,计算公式为

$ \left\{ \begin{array}{l} \mu = \frac{1}{N}\sum\limits_{i = 1}^N {{v_i}} \;\;i = 1,2, \cdots ,N\\ {\sigma ^2} = \frac{1}{{N - 1}}\sum\limits_{i = 1}^N {{{({v_i} - \mu )}^2}} \;\;i = 1,2, \cdots ,N\\ H = - \sum\limits_{k = 1}^c {{p_k}{{\log }_2}{p_k}\;\;k = 1,2, \cdots ,c} \\ I = \int {v(t){\rm{d}}t} \end{array} \right. $ (10)
 

式中:v为信号幅值;N为样本个数;c为样本类别数;pk为第k类样本在总样本中所占的比例。

综合以上4个特征指标,对每个子窗口内的所有数据点,计算特征集Φ中各个特征分量的均值Φ1,其中包含μΔrnμΔncμΔθμΔωz;方差Φ2包含σΔrn2σΔnc2σΔθ2σΔωz2;信息熵Φ3包含HΔrnHΔncHΔθHΔωz;积分Φ4包含IΔrnIΔncIΔθIΔωz;由Φ1Φ2Φ3Φ4组成新的特征集Φ′。

在故障模式f1f2情形下,新构建的特征集Φ′f1Φ′f2图 6所示。

图 6 Φ′f1Φ′f2Φ′h的特征点空间分布 Fig. 6 Spatial distribution of feature points of Φ′f1, Φ′f2 and Φ′h

图 6可以看出,Φ′f1Φ′f2Φ′h在三维空间中相互之间的交集明显减小,经过特征选择处理后,特征更加明显,有利于深度森林建立准确的预测分类模型。

构造的基于深度森林的卫星ACS执行机构与传感器故障识别预测分类模型如图 7所示。

图 7 基于深度森林的卫星ACS执行机构与传感器故障识别预测分类模型 Fig. 7 Fault recognition prediction and classification model of actuators and sensors in satellite ACS based on deep forest algorithm
3.3 深度森林故障识别方法

深度森林算法是在深度学习理论以及深度神经网络的启发下,以随机森林算法为基础的一种有监督机器集成学习算法。作为一种具有一定深度的基于决策树的预测算法,深度森林算法将预测过程分为2个阶段:多粒度扫描阶段(multi-grained scanning)和级联森林阶段(cascade forest),如图 8所示。以下将依次对多粒度扫描阶段和级联森林阶段进行分析。

图 8 深度森林算法整体预测过程 Fig. 8 Whole prediction process of deep forest algorithm

为提高深度森林算法中级联森林阶段的预测效果,深度森林算法设置了多粒度扫描阶段来对样本特征进行提取,尽可能地挖掘序列数据特征的顺序关系;通过设置级联森林阶段,以体现其深度学习的过程。级联森林阶段的每一级都由多个不同类型的森林模型组成。深度森林算法利用级联森林阶段对数据特性逐层进行处理,加强了算法的表征学习能力,有利于提高预测精准度。

基于深度森林的卫星ACS执行机构与传感器故障识别方法的流程如图 9所示,整体流程分为训练和测试2个阶段,具体步骤如下:

图 9 基于深度森林的卫星ACS执行机构与传感器故障识别流程 Fig. 9 Flow chart of fault identification of sensor and actuator satellite ACS based on deep forest algorithm

1) 训练阶段

步骤1  建立ACS训练样本集,并对样本数据进行预处理。

步骤2  根据ACS正常和故障样本数据的特点,对训练样本进行特征选择和特征提取,得到训练样本集的特征数据。

步骤3  利用训练样本集的特征数据对深度森林进行训练。按照深度森林算法的超参数设置,对多粒度扫描阶段和级联森林阶段中的森林模型进行构造,确定级联森林的级数,并生成深度森林故障识别模型。

2) 测试阶段

步骤1  建立ACS测试样本集,并对样本数据进行预处理。

步骤2  采用训练样本的特征选择与特征提取方法处理测试样本,得到测试样本集的特征数据。

步骤3  利用深度森林故障识别模型对测试样本集的特征数据进行预测,得到预测分类结果。

4 仿真试验及结果分析 4.1 试验平台介绍

本文开展了半物理仿真试验以验证所提方法的有效性。基于大理石平台的半物理仿真试验如图 10所示,卫星模拟器的系统组成示意图如图 11所示[21]

图 10 某卫星半物理仿真平台示意图 Fig. 10 Schematic diagram of satellite semi-physical simulation platform
图 11 某卫星半物理仿真平台系统组成示意图[21] Fig. 11 Composition of satellite semi-physical simulation platform system[21]

气浮台卫星ACS模拟系统主要包括控制器(星载计算机)、执行机构和传感器3个部分。星载计算机,通过姿态敏感器反馈的当前气浮台姿态信息,结合地面计算机发来的指令以及星载计算机内的控制算法对执行机构施加控制指令;执行机构为反作用飞轮,通过驱动电路输出的电压对气浮台产生相应的控制力矩。传感器包括姿态敏感器和陀螺仪,姿态敏感器用来实时获取当前气浮台的姿态信息,陀螺仪用来实时获取当前气浮台角速度。

试验所用飞轮的输入输出特性如图 12所示。由图 12可看出,飞轮存在物理受限特性及间隙特性,在试验中,将影响指令力矩与实际输出力矩的一致性,从而影响到闭环的控制状态,对执行机构和传感器故障的辨识提出了挑战。

图 12 不同指令力矩下的飞轮输出特性曲线 Fig. 12 Flywheel output characteristic curves under different command torques
4.2 故障模式设置

由于动量轮中的电机转子需要长时间连续不断地做机械运动,因此比较容易发生故障。动量轮的主要故障有动量轮卡死故障、空转故障、驱动电路故障、增益下降故障、摩擦力矩增大故障。其中动量轮卡死故障和空转故障,会对气浮台姿态机动及稳定造成比较大的影响,故障特征明显;摩擦力矩增大故障发生的主要原因包括轴承温度增高、轴承润滑不好、轴承保持架不稳定、壳体泄露导致内部真空环境密封失效、飞轮压力降低,而这些事件发生的概率很小,不考虑其故障情形。

传感器故障主要有卡死故障、偏差故障和噪声故障3种故障模式。其中传感器卡死故障的故障特征明显,易于故障识别,且传感器噪声故障发生的概率相对较大。故障模式如表 1所示。表中:f1为动量轮驱动电路故障; f2为动量轮增益下降故障; f3为动量轮间隙故障; f4为陀螺仪噪声故障; f5为姿态敏感器噪声故障。no为动量轮实际输出转速; ni为动量轮理论输出转速; K1为动量轮驱动电路故障增益系数,0 < K1 < 1。Mo为动量轮实际输出力矩; Mi为动量轮理论输出力矩; K2为动量轮增益下降故障增益系数,0 < K2 < 1。δ(t)为动量轮故障信号; Δt为故障持续时间;ωo为陀螺仪实际输出; ωi为陀螺仪理论输出; ωf为陀螺仪故障信号。θo为姿态敏感器实际输出;θi为姿态敏感器理论输出; θf为姿态敏感器故障信号; tf为故障发生时刻。

表 1 执行机构与传感器故障模式 Table 1 Fault mode of actuators and sensors

故障
类别
故障数学模型 故障原因
f1 $ \begin{array}{l} {n_{\rm{o}}}(t) = \\ \left\{ {\begin{array}{*{20}{c}} {{n_{\rm{i}}}(t)}&{t < {t_{\rm{f}}}}\\ {{K_1}{n_{\rm{i}}}(t)}&{t \ge {t_{\rm{f}}}} \end{array}, {K_1} < 1} \right. \end{array}$ 驱动电路放大器故障
f2 $ \begin{array}{l} {M_{\rm{o}}}(t) = \\ \left\{ {\begin{array}{*{20}{c}} {{M_{\rm{i}}}(t)}&{t < {t_{\rm{f}}}}\\ {{K_2}{M_{\rm{i}}}(t)}&{t \ge {t_{\rm{f}}}} \end{array}, {K_2} < 1} \right. \end{array}$ 驱动电机故障、由于温度的变化, 元器件老化等因素
f3 $ \begin{array}{l} {M_{\rm{o}}}(t) = \\ \left\{ {\begin{array}{*{20}{c}} {{M_{\rm{i}}}(t)}&{t < {t_{\rm{f}}}}\\ {{M_{\rm{i}}}(t) + \delta (t)}&{{t_{\rm{f}}} \le t \le {t_{\rm{f}}} + \Delta t}\\ {{M_{\rm{i}}}(t)}&{t > {t_{\rm{f}}}} \end{array}} \right. \end{array}$ 总线电压故障或是电机存在间歇的时变故障
f4 $ \begin{array}{l} {\omega _{\rm{o}}}(t) = \\ \left\{ {\begin{array}{*{20}{c}} {{\omega _{\rm{i}}}(t)}&{t < {t_{\rm{f}}}}\\ {{\omega _{\rm{i}}}(t) + {\omega _{\rm{f}}}(t)}&{t \ge {t_{\rm{f}}}} \end{array}} \right. \end{array}$ 陀螺线路盒电压输出不稳, 测量噪声增大
f5 $ \begin{array}{l} {\theta _{\rm{o}}}(t) = \\ \left\{ {\begin{array}{*{20}{c}} {{\theta _{\rm{i}}}(t)}&{t < {t_{\rm{f}}}}\\ {{\theta _{\rm{i}}}(t) + {\theta _{\rm{f}}}(t)}&{t \ge {t_{\rm{f}}}} \end{array}} \right. \end{array}$ 存在干扰目标, CCD出现瑕疵点
4.3 故障定位结果及分析

本文采用保持法对深度森林故障识别模型进行准确性评估,即随机选择70%的原始数据集作为训练集,30%的原始数据集作为测试集。

先将训练集进行数据预处理、特征选择、特征提取,并贴上相应类别标签后输入到深度森林中进行模型训练,随后利用深度森林故障识别模型对测试集的特征点进行分类,判别每个特征点所属状态模式。图 13给出所建立的深度森林故障识别模型对某一组测试集的故障识别结果,其中f0为卫星ACS无故障。

图 13 深度森林故障识别模型对某一组测试集的故障识别结果 Fig. 13 Fault identification results of one test set of deep forest algorithm

图 13可以看出,深度森林故障识别模型可实现绝大部分特征点的类别预测,但仍有少部分特征点出现了误判。针对正常和故障采样数据重复多次随机划分训练集和测试集,并对深度森林故障识别模型的识别结果进行统计,得到深度森林故障识别模型预测类别对真实类别的识别准确率如表 2所示,此处的识别准确率表示当样本特征点的真实类别为fi(i=0, 1, …, 5)时,深度森林故障识别模型的预测类别为fj(j=0, 1, …5)的概率。

表 2 深度森林故障识别模型预测类别对真实类别的识别准确率 Table 2 Accuracy rates of predicted category by deep forest fault recognition model to real category
预测类别 识别准确率/%
f0 f1 f2 f3 f4 f5
f0 98.15 0 0 1.23 0 0
f1 0 100 0 0 0 0
f2 0 0 96.55 0 0 0.42
f3 1.85 0 0 96.33 2.25 0.95
f4 0 0 1.65 1.09 97.75 0
f5 0 0 1.8 1.35 0 98.63

对每种状态模式的识别准确率均达到96.33%,如表 2所示。

表 2中可以看出: ①故障f3与故障f4有少部分特征点出现了相互误判。这是因为f3属于动量轮间隙故障,主要通过在一段时间内产生不定大小的动量轮故障偏差力矩δ(t),从而影响姿态角θ和角速度ωz的输出;② f4属于陀螺仪噪声故障,主要通过产生不确定的角速度偏差ωf(t),从而通过闭环回路影响下一周期的系统输出。两者的故障信号都具有一定的随机性,使得两者在故障特征空间上存在一定的交集。当两者在故障特征空间上相交时,造成误判的发生。

从仿真结果看,本文提出的深度森林故障识别模型具有较高的泛化能力,同时也验证了深度森林算法应用于卫星ACS执行机构与传感器故障识别的有效性。

5 结论

本文将深度森林算法应用于卫星ACS执行机构与传感器的故障识别。通过从少量卫星信号数据中提取能反映故障模式特征的信号特征,借助深度森林算法实现由故障特征参量空间向故障模式空间的映射,从而识别故障模式。通过某卫星半物理仿真平台试验,在气浮台干扰力矩、卫星转动惯量未知、飞轮非线性特性、闭环故障传播等多种不定因素影响下,深度森林算法的预测分类准确性达到了96.33%,可有效地实现卫星ACS执行机构与传感器的故障识别。

参考文献
[1] 王巍, 宝音贺喜格, 姜兴渭, 等. 基于状态观测器的卫星执行机构故障诊断方法研究[J]. 黑龙江自动化技术与应用, 1999(4): 9-12.
WANG W, BAOYIN H X G, JIANG X W, et al. Fault diagnosis method of actuators in satellite based on state observer[J]. Techniques of Automation and Applications, 1999(4): 9-12. (in Chinese)
Cited By in Cnki (4) | Click to display the text
[2] LI C L, CHEN X Q, JIA Q X, et al. FDD of the attitude control system of satellite with actuator fault and time delay via two-stage Kalman filter[C]//3rd International Conference on Machinery Electronics and Control Engineering, 2014, 441: 859-863.
[3] 江文建, 姜斌, 廖鹤, 等. 基于ILLE和SVM的卫星执行机构系统故障检测与定位[J]. 航天控制, 2019(3): 4.
JIANG W J, JIANG B, LIAO H, et al. Fault detection and location of actuator system in satellite based on ILLE and SVM[J]. Aerospace Control, 2019(3): 4. (in Chinese)
Cited By in Cnki (2) | Click to display the text
[4] 胡志坤, 孙岩, 姜斌, 等. 一种基于最优未知输入观测器的故障诊断方法[J]. 自动化学报, 2013, 39(8): 1225-1230.
HU Z K, SUN Y, JIANG B, et al. An optimal unknown input observer based fault diagnosis method[J]. Acta Automatica Sinica, 2013, 39(8): 1225-1230. (in Chinese)
Cited By in Cnki (60) | Click to display the text
[5] 樊雯, 程月华, 姜斌. 基于模糊滑模的卫星姿态控制系统故障诊断[J]. 东南大学学报(自然科学版), 2010, 40(S1): 238-243.
FAN W, CHENG Y H, JIANG B. Fault diagnosis in attitude control system of satellite based on fuzzy sliding mode[J]. Journal of Southeast University (Natural Science Edition), 2010, 40(S1): 238-243. (in Chinese)
Cited By in Cnki | Click to display the text
[6] NADERI E, KHORASANI K. A data-driven approach to actuator and sensor fault detection, isolation and estimation in discrete-time linear systems[J]. Automatica, 2017, 85: 165-178.
Click to display the text
[7] WANG K, CHEN J, SONG Z. Concurrent fault detection and anomaly location in closed-loop dynamic systems with measured disturbances[J]. IEEE Transactions on Automation Science and Engineering, 2018.
Click to display the text
[8] CHENG Y, JIANG B, LU N, et al. Incremental locally linear embedding-based fault detection for satellite attitude control systems[J]. Journal of the Franklin Institute, 2016, 353(1): 17-36.
Click to display the text
[9] 李磊, 高永明, 吴止锾, 等. 基于神经网络观测器的卫星姿态控制系统陀螺故障诊断[J]. 计算机测量与控制, 2018, 26(11): 7-11.
LI L, GAO Y M, WU Z H, et al. Gyroscopic fault diagnosis method for satellite attitude control system based on neural network observer[J]. Computer Measurement & Control, 2018, 26(11): 7-11. (in Chinese)
Cited By in Cnki | Click to display the text
[10] MOUSAVI S, KHORASANI K. Fault detection of reaction wheels in attitude control subsystem of formation flying satellites:A dynamic neural network-based approach[J]. International Journal of Intelligent Unmanned Systems, 2014, 2(1): 2-26.
Click to display the text
[11] 陈辛, 魏炳翌, 闻新. 基于支持向量机的卫星执行机构故障诊断研究[J]. 中国空间科学技术, 2018, 38(2): 47-55.
CHEN X, WEI B Y, WEN X. Study of support vector machine based faults diagnosis for satellite's actuators[J]. Chinese Space Science and Technology, 2018, 38(2): 47-55. (in Chinese)
Cited By in Cnki | Click to display the text
[12] JEGADEESHWARAN R, SUGUMARAN V. Fault diagnosis of automobile hydraulic brake system using statistical features and support vector machines[J]. Mechanical Systems and Signal Processing, 2015, 52: 436-446.
Click to display the text
[13] CERRADA M, ZURITA G, CABRERA D, et al. Fault diagnosis in spur gears based on genetic algorithm and random forest[J]. Mechanical Systems and Signal Processing, 2016, 70: 87-103.
Click to display the text
[14] GRAVES A, MOHAMED A, HINTON G. Speech recognition with deep recurrent neural networks[C]//2013 IEEE International Conference on Acoustics, Speech and Signal Processing. Piscataway: IEEE Press, 2013: 6645-6649.
[15] ZHOU Z H, FENG J. Deep forest:Towards an alternative to deep neural networks. arXiv 2017[J]. arXiv preprint arXiv:1702.08835.
Click to display the text
[16] 陈吕鹏, 殷林飞, 余涛, 等. 基于深度森林算法的电力系统短期负荷预测[J]. 电力建设, 2018, 39(11): 42-50.
CHEN L P, YIN L F, YU T, et al. Short-term power load forecasting based on deep forest algorithm[J]. Electric Power Construction, 2018, 39(11): 42-50. (in Chinese)
Cited By in Cnki | Click to display the text
[17] HU G, LI H, XIA Y, et al. A deep Boltzmann machine and multi-grained scanning forest ensemble collaborative method and its application to industrial fault diagnosis[J]. Computers in Industry, 2018, 100: 287-296.
Click to display the text
[18] 沈文博, 孙荣霞, 马少卿, 等. 深度森林与人工神经网络在光伏出力预测的比较[J]. 信息技术与网络安全, 2018, 37(4): 49-51.
SHEN W B, SUN R X, MA S Q, et al. Comparison of deep forest and artificial neural network in prediction of PV output[J]. Information Technology and Network Security, 2018, 37(4): 49-51. (in Chinese)
Cited By in Cnki | Click to display the text
[19] LI M, ZHANG N, PAN B, et al. Hyperspectral image classification based on deep forest and spectral-spatial cooperative feature[C]//International Conference on Image and Graphics. Berlin: Springer, 2017: 325-336.
[20] UTKIN L V, KOVALEV M S, MELDO A A. A deep forest classifier with weights of class probability distribution subsets[J]. Knowledge-Based Systems, 2019, 173: 15-27.
Click to display the text
[21] 康国华, 杨炳辉, 刘瑶, 等. 基于微小卫星编队的聚光操控技术[J]. 航天控制, 2018, 36(2): 42-47.
KANG G H, YANG B H, LIU Y, et al. The technology of light control based on small satellite formation[J]. Aerospace Control, 2018, 36(2): 42-47. (in Chinese)
Cited By in Cnki | Click to display the text
http://dx.doi.org/10.7527/S1000-6893.2019.23778
中国航空学会和北京航空航天大学主办。
0

文章信息

程月华, 江文建, 杨浩, 薛琪, 廖鹤
CHENG Yuehua, JIANG Wenjian, YANG Hao, XUE Qi, LIAO He
基于深度森林的卫星ACS执行机构与传感器故障识别
Fault identification of actuators and sensors of satellite attitude control systems based on deep forest algorithm
航空学报, 2020, 41(S1): 723778.
Acta Aeronautica et Astronautica Sinica, 2020, 41(S1): 723778.
http://dx.doi.org/10.7527/S1000-6893.2019.23778

文章历史

收稿日期: 2019-12-13
退修日期: 2019-12-26
录用日期: 2019-12-31
网络出版时间: 2020-01-03 13:53

相关文章

工作空间