CN115824261A - 一种控制力矩陀螺仪故障检测方法、装置及相关存储介质 - Google Patents
一种控制力矩陀螺仪故障检测方法、装置及相关存储介质 Download PDFInfo
- Publication number
- CN115824261A CN115824261A CN202310060890.4A CN202310060890A CN115824261A CN 115824261 A CN115824261 A CN 115824261A CN 202310060890 A CN202310060890 A CN 202310060890A CN 115824261 A CN115824261 A CN 115824261A
- Authority
- CN
- China
- Prior art keywords
- data
- control moment
- moment gyroscope
- sinc
- fault detection
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 65
- 238000012549 training Methods 0.000 claims abstract description 40
- 238000000034 method Methods 0.000 claims abstract description 37
- 238000005259 measurement Methods 0.000 claims abstract description 8
- 230000015654 memory Effects 0.000 claims description 16
- 230000005856 abnormality Effects 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 10
- 238000000605 extraction Methods 0.000 claims description 5
- 230000011218 segmentation Effects 0.000 claims description 5
- 230000002159 abnormal effect Effects 0.000 abstract description 16
- 238000004364 calculation method Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000013526 transfer learning Methods 0.000 description 5
- 239000002689 soil Substances 0.000 description 4
- 101000794020 Homo sapiens Bromodomain-containing protein 8 Proteins 0.000 description 3
- 101001006782 Homo sapiens Kinesin-associated protein 3 Proteins 0.000 description 3
- 101000615355 Homo sapiens Small acidic protein Proteins 0.000 description 3
- WGKGADVPRVLHHZ-ZHRMCQFGSA-N N-[(1R,2R,3S)-2-hydroxy-3-phenoxazin-10-ylcyclohexyl]-4-(trifluoromethoxy)benzenesulfonamide Chemical compound O[C@H]1[C@@H](CCC[C@@H]1N1C2=CC=CC=C2OC2=C1C=CC=C2)NS(=O)(=O)C1=CC=C(OC(F)(F)F)C=C1 WGKGADVPRVLHHZ-ZHRMCQFGSA-N 0.000 description 3
- 102100021255 Small acidic protein Human genes 0.000 description 3
- 238000013528 artificial neural network Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000007774 longterm Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000013500 data storage Methods 0.000 description 2
- 238000013135 deep learning Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000007787 long-term memory Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 230000000306 recurrent effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000013024 troubleshooting Methods 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Abstract
本发明公开了一种控制力矩陀螺仪故障检测方法、装置及相关存储介质，该方法包括：搭建Sinc‑LSTM网络，采用大规模IMS轴承故障检测数据集和火星车数据集进行两阶段的预训练；获取在轨卫星的控制力矩陀螺仪遥测数据，并进行工况分类预处理；将经工况分类预处理后的控制力矩陀螺仪遥测数据，对所述Sinc‑LSTM网络进行训练，到能够实现控制力矩陀螺仪遥测数据拟合重建的模型；基于所述控制力矩陀螺仪遥测数据拟合重建的模型，对需要检测的控制力矩陀螺仪进行检测；将得到的重建数据与其对应的原始遥测数据作差，阈值判断完成故障检测。该方法提高了控制力矩陀螺仪异常检测的精确度。
Description
技术领域
本发明涉及航天器故障检测技术领域，尤其是涉及基于迁移学习的控制力矩陀螺仪故障检测技术领域，特别涉及一种控制力矩陀螺仪故障检测方法、装置及相关存储介质。
背景技术
控制力矩陀螺仪是航天器姿态调整的惯性执行部件，其稳定工作对航天器的在轨运行具有重要意义。控制力矩陀螺仪的异常检测是卫星在轨状态检测与健康管理技术的基础和重要环节，能够保障卫星在轨的安全运行，以及故障的及时排查，因此具有重大的研究价值。
随着航天技术的发展，航天器系统越来越复杂，数据规模也大幅度提升。对控制力矩陀螺仪等航天器部件来说,传统的依靠人工判别和专家系统判别的异常检测方式会消耗越来越多的人力，并且在海量数据规模下难以提升精度。因此，对于控制力矩陀螺仪的异常检测需要采用深度学习的方式来提升检测效率和精度。
发明内容
本发明的目的在于提供一种控制力矩陀螺仪故障检测方法、装置及相关存储介质，解决了现有技术的人力检测精度低、效率低的问题。
为实现上述目的，本发明采取的技术方案为：
第一方面，本发明实施例提供一种控制力矩陀螺仪故障检测方法，包括以下步骤：
S10、搭建Sinc-LSTM网络，采用大规模IMS轴承故障检测数据集和火星车数据集进行两阶段的预训练；
S20、获取在轨卫星的控制力矩陀螺仪遥测数据，并进行工况分类预处理；
S30、将经工况分类预处理后的控制力矩陀螺仪遥测数据，对所述Sinc-LSTM网络进行训练，得到能够实现控制力矩陀螺仪遥测数据拟合重建的模型；
S40、基于所述控制力矩陀螺仪遥测数据拟合重建的模型，对需要检测的控制力矩陀螺仪进行检测；将得到的重建数据与其对应的原始遥测数据作差，进行自适应的阈值判断完成故障检测。
进一步地，所述步骤S10包括：
S101、搭建Sinc-LSTM网络；
S102、采用大规模IMS轴承故障检测数据集，输入所述Sinc-LSTM网络进行第一阶段预训练；
S103、将火星车数据集输入经过第一阶段预训练的Sinc-LSTM网络中，进行第二阶段预训练。
进一步地，所述步骤S20包括：
获取在轨卫星的控制力矩陀螺仪全周期时序数据进行工况分类，得到不同工况下的控制力矩陀螺仪运行数据。
进一步地，工况分类选用基于密度峰值的自动聚类算法。
进一步地，所述步骤S30中，Sinc-LSTM网络首层利用sinc函数滤波器作为一维卷积核进行特征提取。
进一步地，所述Sinc-LSTM网络中，采用sinc函数的滤波器对LTSM网络进行改进，增强LSTM网络对一维时序数据的拟合重建能力。
进一步地，所述Sinc-LSTM网络中隐藏层的第一层由16×6个SincNet滤波器构成，LSTM网络部分由两层96个单元的单向LSTM网络组成，数据输入全连接层后连接输出层完成训练。
进一步地，所述步骤S40中，将得到的重建数据与其对应的原始遥测数据作差，在差值数据计算过程中，采用指数加权平均的方式对差值序列进行平滑处理；
所述步骤S40中，阈值的设置根据工况分段动态变化，各个工况分段的阈值设置由以下公式决定：
其中，
表示第i个工况类别的异常判断阈值，
表示该工况分段拟合重建数据与原始遥测数据差值的均值，
和
分别表示该工况分段差值的标准差和全寿命周期误差的标准差，
和
表示一对提前设置的常量。
第二方面，本发明实施例还提供一种控制力矩陀螺仪异常检测装置，包括：
遥测数据获取设备，用以采集在轨卫星控制力矩陀螺仪运行数据；
存储器，该存储器中存储有计算机指令；
处理器，与所述运行数据获取设备、存储器数据连接，其通过执行所述计算机指令，从而执行如上述实施例任一项所述的控制力矩陀螺仪异常检测方法，对控制力矩陀螺仪数据进行异常检测。
第三方面，本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有计算机指令，所述计算机指令在被处理器执行时实现如上述实施例任一项所述的控制力矩陀螺仪异常检测方法。
与现有技术相比，本发明具有如下有益效果：
本发明实施例提供的一种控制力矩陀螺仪故障检测方法，包括：搭建Sinc-LSTM网络，采用大规模IMS轴承故障检测数据集和火星车数据集进行两阶段的预训练；获取在轨卫星的控制力矩陀螺仪遥测数据，并进行工况分类预处理；将经工况分类预处理后的控制力矩陀螺仪遥测数据，对所述Sinc-LSTM网络进行训练，得到能够实现控制力矩陀螺仪遥测数据拟合重建的模型；基于所述控制力矩陀螺仪遥测数据拟合重建的模型，对需要检测的控制力矩陀螺仪进行检测；将得到的重建数据与其对应的原始遥测数据作差，进行自适应的阈值判断完成故障检测。该方法提高了控制力矩陀螺仪异常检测的精确度。
附图说明
图1为本发明实施例提供的控制力矩陀螺仪故障检测方法的流程图；
图2为本发明实施例提供的基于迁移学习和工况分类的Sinc-LSTM控制力矩陀螺仪异常检测网络模型示意图；
图3为本发明实施例提供的一组Sinc函数滤波器示意图；
图4为本发明实施例提供的LSTM单元模型示意图；
图5a为本发明实施例提供的原始遥测数据y和预测数据y_hat的对比图；
图5b为本发明实施例提供的差值序列的示意图；
图6为本发明实施例提供的控制力矩陀螺仪异常检测装置结构示意图。
具体实施方式
为使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解，下面结合具体实施方式，进一步阐述本发明。
在本发明的描述中，需要说明的是，术语“上”、“下”、“内”、“外”“前端”、“后端”、“两端”、“一端”、“另一端”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。此外，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性。
在本发明的描述中，需要说明的是，除非另有明确的规定和限定，术语“安装”、“设置有”、“连接”等，应做广义理解，例如“连接”，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。
参照图1所示，本发明提供的一种控制力矩陀螺仪故障检测方法，基于迁移学习与工况分类的Sinc-LSTM控制力矩陀螺仪异常检测网络进行的检测，具体包括如下步骤：
S10、搭建Sinc-LSTM网络，采用大规模IMS轴承故障检测数据集和火星车数据集进行两阶段的预训练；
S20、获取在轨卫星的控制力矩陀螺仪遥测数据，并进行工况分类预处理；
S30、将经工况分类预处理后的控制力矩陀螺仪遥测数据，对所述Sinc-LSTM网络进行训练，得到能够实现控制力矩陀螺仪遥测数据拟合重建的模型；
S40、基于所述控制力矩陀螺仪遥测数据拟合重建的模型，对需要检测的控制力矩陀螺仪进行检测；将得到的重建数据与其对应的原始遥测数据作差，进行自适应的阈值判断完成故障检测。
本实施例中，基于搭建的Sinc-LSTM网络进行预训练，在大规模轴承故障检测数据集IMS以及小规模同类型火星土壤水分主被动卫星及火星车故障数据集进行两阶段的预训练，完成对一维时序数据一般化特征和同类型卫星及火星车特征的学习；对获取的控制力矩陀螺仪全周期数据进行工况自适应聚类，完成对不同时期异常检测阈值的自动区分；利用SincNet网络的Sinc函数滤波器提取控制力矩陀螺仪一维时序数据的特征，作为长短期记忆网络（Long-Short-Term-Memory netwrok，LSTM）的第一层，增强LSTM网络的拟合能力；得到能够实现控制力矩陀螺仪遥测数据拟合重建的模型，基于此，对需要检测的控制力矩陀螺仪进行检测；将得到的重建数据与其对应的原始遥测数据作差，进行自适应的阈值判断完成故障检测，可有效提高控制力矩陀螺仪异常检测的精确度。
下面分别对上述各个步骤进行详细的描述：
步骤1、对网络模型进行两阶段预训练；
对控制力矩陀螺仪异常检测模型分别采用大规模IMS轴承故障检测数据集和同类型火星车（Soil Moisture Active Passive & the Mars Science Laboratory rover,Curiosity）SMAP&MSL故障检测数据集进行两阶段的预训练，具体过程如下：
首先，将来自辛辛那提大学的IMS轴承故障检测数据集输入搭建好的Sinc-LSTM网络中。IMS数据集一共包含三部分，分别是四组轴承在2003-2004年的三个时段做的全寿命周期震动实验数据。将三个数据集合并后输入网络，设置合适的超参数后进行训练。
其次，将来自NASA的SMAP&MSL土壤水分主被动卫星及火星科学实验室火星车数据输入经过第一阶段预训练的Sinc-LSTM网络中。经过第一阶段预训练的Sinc函数滤波器和LSTM单元内部参数作为训练的初始权重参数。
该步骤中，通过预训练，可学习一维时序数据一般化故障和同类型卫星及火星车故障的特征，完成迁移学习的特征空间映射。
步骤2、对控制力矩陀螺仪进行自适应工况分类的预处理；
获取在轨卫星的控制力矩陀螺仪全周期时序数据进行聚类分段，得到不同工况下的控制力矩陀螺仪运行数据，以适应不同时期故障检测阈值变化。
本实施例对控制力矩陀螺仪分段数据进行自动聚类选取的方法是基于密度峰值的自动聚类算法。
基于密度峰值聚类算法是通过计算某个点周围点的密度，以及与周围点的距离来判断属于某一类。该算法认为，对于每个类的中心点，其周围点的密度应当比该类其他点的密度大。同时，本类中心点应该距离其他类中心点要足够远。局部密度的具体公式如下：
假设数据点
的局部密度为
，数据点
到局部密度比它大且距离最近的数据点
的距离为
，则有如下定义：
其中，
表示
和
之间的距离；
表示逻辑判断函数；如果
，则
，其他情况时，
。同时，
为截止距离。
相当于与ｉ点的距离小于
的点的数量。
计算距离具体公式如下。
即，当某点的密度不为最大时，该点不是中心点，则将距离设定为它到离它最近点的距离，当某点的密度最大时，说明该点是中心点，则将该点所对应的距离设为它到离它最远的点的距离。
以密度
作为横轴，距离
作为纵轴，绘制一幅点图，越靠近右上角的点则拥有最大的局部密度和最大的距离，可以认定为簇的中心点。该算法不需要迭代计算就可以将数据划分为合适的簇，同时也不需要预设中心点，具有较少的经验成分。
本实施例中，使用基于密度峰值的自动聚类方法，得到不同工况阶段的分类结果。比如控制力矩陀螺仪的数据采集了控制力矩陀螺仪全寿命周期中均匀分布的155天的数据，将所有共16个通道的数据按小时分段分为3720个数据段，然后采用密度峰值聚类方法对各段数据进行工况分类，在训练过程中相同工况阶段的数据加载入同一个batch输入网络训练。在遥测数据拟合重建完成后根据不同的工况阶段分别进行阈值判断异常检测。
步骤3、将步骤2经过工况分类预处理后的遥测数据输入Sinc-LSTM网络进行拟合重建：
Sinc-LSTM网络结构如图2所示。其中，隐藏层的第一层由16×6 = 96个sinc函数滤波器组成，LSTM网络部分由两层96个单元的单向LSTM网络组成，数据输入全连接层后连接输出层完成训练。
Sinc函数滤波器从SincNet网络中获得。Sinc函数滤波器的结构如图3所示。SincNet最初用于处理语音信号，Sinc函数滤波器的形状类似于一维时序数据中的波峰波谷交替的形式，因此对于语音信号的一维时序数据有更好的拟合能力。Sinc函数滤波器的控制参数主要是滤波器的截止频率，对于语音等信号的深度学习具有一定的可解释性，可以解决深度网络普遍存在的“黑盒子”问题。如图2所示，在隐藏层的Sinc函数滤波器层对控制力矩陀螺仪的16个通道数据分别设置六个具有不同初始截止频率的Sinc函数滤波器，以提取各个通道上不同频段的特征信息。即：SincNet网络首层利用sinc函数滤波器作为一维卷积核进行特征提取。
LSTM神经网络单元结构如图4所示。LSTM长短期记忆网络属于RNN循环神经网络的一种，对于顺序输入的时序数据具有天生的优势。并且通过输入门输出门和遗忘门的设置，改进了RNN会遗忘长期信息的缺点，输入信息长期的特征进行保留，有利于一维时序数据的拟合重建。在两个LSTM层设置了96个LSTM单元，分别对来自sinc函数层各个通道和各个频段的参数进行长短期的特征提取。并在最后输入全连接层完成序列的拟合重建。
拟合重建的方式是对控制力矩陀螺仪各个通道的数据单独建模分析。控制力矩陀螺仪遥测数据，被制作成为了16个数据文档，每个文档包含其余15个通道的数据作为训练集和单一通道的数据作为监督。Sinc-LSTM网络的训练过程中以其余15个通道的信息作为输入训练的数据。训练完成后16个通道的模型参数被分开存放在不同的文件中。
本实施例中，采用SincNet滤波器提取一维时序数据的特征，采用LSTM网络保持对时序的长期记忆能力；采用sinc函数的滤波器对LTSM网络进行改进，增强了LSTM网络对一维时序数据的拟合重建能力，并能够在一定程度上提高信噪比。
步骤4、重建数据与原始遥测数据作差，进行自适应的阈值判断完成故障检测；
根据工况分类结果，不同工况阶段下的阈值设置不同，进行异常检测时需要根据本段数据所处的工况阶段来进行阈值判断。
如图5a-5b所示，重建数据和原始数据作差的图像如图，作差结果超过阈值后，判定为出现异常。
差值数据计算过程中，采用指数加权平均的方式对差值序列进行平滑处理，避免差值计算过程中控制力矩陀螺仪遥测数据突变等因素产生的噪声被误认为异常，产生误报导致精度下降。
同时工况分段阈值的设置是动态变化的，各个工况分段的阈值设置由以下公式决定：
其中，
表示第i个工况类别的异常判断阈值，
代表该工况分段拟合重建数据与原始遥测数据差值的均值，
和
分别代表该工况分段差值的标准差和全寿命周期误差的标准差，
和
是一对提前设置的常量，分别取5%和20%，分别代表能够接受最低的异常检测阈值和由当前工况分段数据波动性决定的异常检测阈值提升。
以下再次通过举例方式对控制力矩陀螺仪Sinc-LSTM异常检测网络的训练及异常检测方法进行说明，该方法包括以下步骤：
步骤（1）对Sinc-LSTM异常检测网络进行两阶段预训练。
首先，将来自辛辛那提大学的IMS轴承故障检测数据集输入搭建好的Sinc-LSTM网络中。IMS数据集一共包含三部分，分别是四组轴承在2003-2004年的三个时段做的全寿命周期震动实验数据。将三个数据集合并后输入网络，设置合适的超参数后进行训练。
其次，将来自NASA的SMAP&MSL土壤水分主被动卫星及火星科学实验室火星车数据输入经过第一阶段预训练的Sinc-LSTM网络中，第一阶段预训练完成后的模型参数保持不变作为第二阶段预训练的初始权重参数。经过第二阶段预训练的Sinc函数滤波器和LSTM单元内部参数作为控制力矩陀螺仪数据训练的初始权重参数。
步骤（2）、获取在轨卫星的控制力矩陀螺仪遥测数据。
遥测数据为某型号卫星所采集的控制力矩陀螺仪运行数据，包括框架电流，电机电流，轴温，航天器姿态信息，框架角度，框架转动指令等，故障及异常信息由专家标注得到。
步骤（3）、对获得的控制力矩陀螺仪遥测数据进行工况分类预处理。
该工况分类预处理具体包括：
首先，将获得的原始遥测数据按小时进行分段切片，然后进行自动聚类，区分控制力矩陀螺仪不同工况阶段的数据。并且将相同工况阶段的数据段合为同一个batch输入Sinc-LSTM网络进行训练。
其中，自动聚类方法选取的是基于密度峰值的自动聚类算法。
步骤（4）、使用完成预处理的控制力矩陀螺仪遥测数据对Sinc-LSTM网络模型进行训练，得到能够实现控制力矩陀螺仪遥测数据拟合重建的模型。
将完成预处理的遥测数据输入到图2所示网络模型中。
使用该模型对控制力矩遥测数据的16个通道完成单独的拟合重建模型的建立。
通过对比重建序列结果与原始遥测数据之间的区别，得到重建误差loss值。
最后，将重建误差loss值传回网络中，进行梯度反馈，修正网络中各个模块单元参数。
通过预定次数的循环计算，使得最终重建误差在设定阈值范围内，结果准确率在期望范围内，得到具有较准确拟合重建遥测数据的网络模型。
最终测试重建结果则依据重建拟合序列和原始遥测数据的差值，计算差值的均值和标准差，以二者的综合判断指标作为拟合重建结果的准确率。
步骤（5）、重建数据与原始遥测数据作差，阈值判断完成故障检测。
差值数据计算过程中，采用指数加权平均的方式对差值序列进行平滑处理，避免差值计算过程中控制力矩陀螺仪遥测数据突变等因素产生的噪声被误认为异常，产生误报导致精度下降。
同时工况分段阈值的设置是动态变化的，各个工况分段的阈值设置由以下公式决定：
其中，
表示第i个工况类别的异常判断阈值，
代表该工况分段拟合重建数据与原始遥测数据差值的均值，
和
分别代表该工况分段差值的标准差和全寿命周期误差的标准差，
和
是一对提前设置的常量，分别取5%和20%，分别代表能够接受最低的异常检测阈值和由当前工况分段数据波动性决定的异常检测阈值提升。
尽管上述实施例以步骤先后顺序描述了上述方法，然而本领域技术人员可以理解，上述各方法的步骤并非严格按照上述先后顺序执行。
训练好的拟合重建模型具有最优组合参数，将遥测数据输入网络模型中，利用已经调整好参数的Sinc函数层和LSTM层对遥测数据进行重建拟合，实现特征提取，最后完成与原始遥测数据的作差和自适应阈值判断。
同时，本发明所采用的网络模型中，sinc函数滤波器由于其可调参数为最高和最低截止频率，因此该部分卷积核具有较少的参数，能够实现更加快速的拟合，减少运行时间。并且，该sinc函数滤波器对一维时序数据有较大的敏感性，其明确的物理含义使得神经网络具有更好的可解释性。
此外，本发明所采用的网络模型中，通过加入LSTM层，可以增强网络对时序数据长短期特征的提取能力，能够更好地拟合长周期的在轨卫星控制力矩陀螺仪遥测数据，从而提升对于异常情况的检测能力。
基于同一发明构思，本发明还提供了一种控制力矩陀螺仪异常检测装置，基于迁移学习和工况分类的控制力矩陀螺仪异常检测装置，可实现在轨卫星控制力矩陀螺仪的异常检测。该异常检测装置可以由软件和/或硬件实现。
如图6所示，该装置包括遥测数据设备、存储器以及处理器。其中，遥测数据获取设备、存储器以及处理器可以通过总线或者其他方式连接。
其中，遥测数据获取设备用以采集电流、轴温、框架角度等控制力矩陀螺仪运行数据，并将该部分数据发送至处理器。
处理器可以为中央处理器(Central Processing Unit，CPU)，或者其他通用处理器、数字信号处理器(Digital Signal Processor，DSP)、专用集成电路(ApplicationSpecific Integrated Circuit，ASIC)、现场可编程门阵列(Field-Programmable GateArray，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等芯片，或者上述各类芯片的组合。
存储器作为一种非暂态计算机可读存储介质，可用于存储非暂态软件程序、非暂态计算机可执行程序以及模块，如上述实施例中的控制力矩陀螺仪异常检测方法对应的程序或指令。
处理器通过运行存储在存储器中的非暂态软件程序或指令，从而执行处理器的各种功能应用以及数据处理，即实现上述方法实施例中的控制力矩陀螺仪异常检测方法。
存储器可以包括程序存储区和数据存储区，其中，程序存储区可存储操作系统、至少一个功能所需要的应用程序；数据存储区可存储处理器所创建的数据等。
此外，存储器可以包括高速随机存取存储器，还可以包括非暂态存储器，例如至少一个磁盘存储器件、闪存器件、或其他非暂态固态存储器件。
在一些方面中，存储器可选地包括相对于处理器远程设置的存储器，这些远程存储器可以通过网络连接至处理器。
可选地，上述网络包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。
Claims (10)
1.一种控制力矩陀螺仪故障检测方法，其特征在于，包括以下步骤：
S10、搭建Sinc-LSTM网络，采用大规模IMS轴承故障检测数据集和火星车数据集进行两阶段的预训练；
S20、获取在轨卫星的控制力矩陀螺仪遥测数据，并进行工况分类预处理；
S30、将经工况分类预处理后的控制力矩陀螺仪遥测数据，对所述Sinc-LSTM网络进行训练，得到能够实现控制力矩陀螺仪遥测数据拟合重建的模型；
S40、基于所述控制力矩陀螺仪遥测数据拟合重建的模型，对需要检测的控制力矩陀螺仪进行检测；将得到的重建数据与其对应的原始遥测数据作差，进行自适应的阈值判断完成故障检测。
2.根据权利要求1所述的一种控制力矩陀螺仪故障检测方法，其特征在于，所述步骤S10包括：
S101、搭建Sinc-LSTM网络；
S102、采用大规模IMS轴承故障检测数据集，输入所述Sinc-LSTM网络进行第一阶段预训练；
S103、将火星车数据集输入经过第一阶段预训练的Sinc-LSTM网络中，进行第二阶段预训练。
3.根据权利要求1所述的一种控制力矩陀螺仪故障检测方法，其特征在于，所述步骤S20包括：
获取在轨卫星的控制力矩陀螺仪全周期时序数据进行工况分类，得到不同工况下的控制力矩陀螺仪运行数据。
4.根据权利要求3所述的一种控制力矩陀螺仪故障检测方法，其特征在于，工况分类选用基于密度峰值的自动聚类算法。
5.根据权利要求1所述的一种控制力矩陀螺仪故障检测方法，其特征在于，所述步骤S30中，Sinc-LSTM网络首层利用sinc函数滤波器作为一维卷积核进行特征提取。
6.根据权利要求5所述的一种控制力矩陀螺仪故障检测方法，其特征在于，所述Sinc-LSTM网络中，采用sinc函数的滤波器对LTSM网络进行改进，增强LSTM网络对一维时序数据的拟合重建能力。
7.根据权利要求6所述的一种控制力矩陀螺仪故障检测方法，其特征在于，所述Sinc-LSTM网络中隐藏层的第一层由16×6个SincNet滤波器构成，LSTM网络部分由两层96个单元的单向LSTM网络组成，数据输入全连接层后连接输出层完成训练。
8.根据权利要求3所述的一种控制力矩陀螺仪故障检测方法，其特征在于，所述步骤S40中，将得到的重建数据与其对应的原始遥测数据作差，在差值数据计算过程中，采用指数加权平均的方式对差值序列进行平滑处理；
所述步骤S40中，阈值的设置根据工况分段动态变化，各个工况分段的阈值设置由以下公式决定：
其中，
表示第i个工况类别的异常判断阈值，
表示该工况分段拟合重建数据与原始遥测数据差值的均值，
和
分别表示该工况分段差值的标准差和全寿命周期误差的标准差，
和
表示一对提前设置的常量。
9.一种控制力矩陀螺仪异常检测装置，其特征在于，包括：
遥测数据获取设备，用以采集在轨卫星控制力矩陀螺仪运行数据；
存储器，该存储器中存储有计算机指令；
处理器，与所述运行数据获取设备、存储器数据连接，其通过执行所述计算机指令，从而执行如权利要求1-8任一项所述的控制力矩陀螺仪故障检测方法，对控制力矩陀螺仪数据进行异常检测。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有计算机指令，所述计算机指令在被处理器执行时实现如权利要求1-8任一项所述的控制力矩陀螺仪故障检测方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310060890.4A CN115824261B (zh) | 2023-01-15 | 2023-01-15 | 一种控制力矩陀螺仪故障检测方法、装置及相关存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310060890.4A CN115824261B (zh) | 2023-01-15 | 2023-01-15 | 一种控制力矩陀螺仪故障检测方法、装置及相关存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115824261A true CN115824261A (zh) | 2023-03-21 |
| CN115824261B CN115824261B (zh) | 2023-05-23 |
Family
ID=85520789
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310060890.4A Active CN115824261B (zh) | 2023-01-15 | 2023-01-15 | 一种控制力矩陀螺仪故障检测方法、装置及相关存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115824261B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116202558A (zh) * | 2023-05-04 | 2023-06-02 | 中国西安卫星测控中心 | 基于增量式数据统计的cmg转动部件工况检测方法 |
| CN117520947A (zh) * | 2023-11-21 | 2024-02-06 | 中国科学院半导体研究所 | 异常监测与故障识别的方法及装置 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105205288A (zh) * | 2015-10-28 | 2015-12-30 | 中国人民解放军国防科学技术大学 | 基于模式演化的卫星长期在轨运行状态的预测方法 |
| CN109934337A (zh) * | 2019-03-14 | 2019-06-25 | 哈尔滨工业大学 | 一种基于集成lstm的航天器遥测数据异常的检测方法 |
| CN110378052A (zh) * | 2019-07-25 | 2019-10-25 | 北京航空航天大学 | 基于循环神经网络考虑未来工况的设备剩余寿命预测方法 |
| CN113392881A (zh) * | 2021-05-27 | 2021-09-14 | 重庆大学 | 一种基于迁移学习的旋转机械故障诊断方法 |
| CN113567858A (zh) * | 2021-07-30 | 2021-10-29 | 北京航空航天大学 | 控制力矩陀螺剩余寿命预测系统 |
| WO2022074052A1 (en) * | 2020-10-06 | 2022-04-14 | Airbus Defence And Space Limited | Fault detection and isolation using a neural network |
| CN115205582A (zh) * | 2022-06-09 | 2022-10-18 | 复旦大学 | 一种拟图智能化航空传感器故障检测与分类方法 |
-
2023
- 2023-01-15 CN CN202310060890.4A patent/CN115824261B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105205288A (zh) * | 2015-10-28 | 2015-12-30 | 中国人民解放军国防科学技术大学 | 基于模式演化的卫星长期在轨运行状态的预测方法 |
| CN109934337A (zh) * | 2019-03-14 | 2019-06-25 | 哈尔滨工业大学 | 一种基于集成lstm的航天器遥测数据异常的检测方法 |
| CN110378052A (zh) * | 2019-07-25 | 2019-10-25 | 北京航空航天大学 | 基于循环神经网络考虑未来工况的设备剩余寿命预测方法 |
| WO2022074052A1 (en) * | 2020-10-06 | 2022-04-14 | Airbus Defence And Space Limited | Fault detection and isolation using a neural network |
| CN113392881A (zh) * | 2021-05-27 | 2021-09-14 | 重庆大学 | 一种基于迁移学习的旋转机械故障诊断方法 |
| CN113567858A (zh) * | 2021-07-30 | 2021-10-29 | 北京航空航天大学 | 控制力矩陀螺剩余寿命预测系统 |
| CN115205582A (zh) * | 2022-06-09 | 2022-10-18 | 复旦大学 | 一种拟图智能化航空传感器故障检测与分类方法 |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116202558A (zh) * | 2023-05-04 | 2023-06-02 | 中国西安卫星测控中心 | 基于增量式数据统计的cmg转动部件工况检测方法 |
| CN116202558B (zh) * | 2023-05-04 | 2023-08-01 | 中国西安卫星测控中心 | 基于增量式数据统计的cmg转动部件工况检测方法 |
| CN117520947A (zh) * | 2023-11-21 | 2024-02-06 | 中国科学院半导体研究所 | 异常监测与故障识别的方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115824261B (zh) | 2023-05-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN115824261B (zh) | 一种控制力矩陀螺仪故障检测方法、装置及相关存储介质 | |
| CN110162018B (zh) | 基于知识蒸馏与隐含层共享的增量式设备故障诊断方法 | |
| CN111340238B (zh) | 一种工业系统的故障诊断方法、装置、设备及存储介质 | |
| CN108334936B (zh) | 基于迁移卷积神经网络的故障预测方法 | |
| CN110135637A (zh) | 基于长短期记忆与自适应提升的微电网短期负荷预测方法 | |
| CN111931806A (zh) | 一种多传感器数据融合的设备故障诊断方法及装置 | |
| CN113992718B (zh) | 一种基于动态宽度图神经网络的群体传感器异常数据检测方法和系统 | |
| CN114363195A (zh) | 面向时间和频谱残差卷积网络的网络流量预测预警方法 | |
| CN113657668A (zh) | 基于lstm网络的电力负荷预测方法及系统 | |
| CN115689008A (zh) | 基于集合经验模态分解的CNN-BiLSTM短期光伏功率预测方法及其系统 | |
| CN111008726A (zh) | 一种电力负荷预测中类图片转换方法 | |
| CN116799796A (zh) | 一种光伏发电功率预测方法、装置、设备及介质 | |
| CN115018170A (zh) | 一种短期风电功率预测方法、装置、设备及存储介质 | |
| CN116702992A (zh) | 一种发电功率预测方法、装置、电子设备及存储介质 | |
| CN115290326A (zh) | 一种滚动轴承故障智能诊断方法 | |
| CN114330814A (zh) | 一种基于vmd分解和改进双层bilstm网络的短期负荷预测方法 | |
| CN116706907B (zh) | 基于模糊推理的光伏发电预测方法和相关设备 | |
| CN114139783A (zh) | 基于非线性加权组合的风电短期功率预测方法及装置 | |
| CN113780644A (zh) | 一种基于在线学习的光伏出力预测方法 | |
| CN115392582A (zh) | 基于增量模糊粗糙集属性约简的作物产量预测方法 | |
| CN115346084A (zh) | 样本处理方法、装置、电子设备、存储介质及程序产品 | |
| CN114648060A (zh) | 基于机器学习的故障信号规范化处理及分类方法 | |
| CN117713688B (zh) | 光伏组件在多朝向及倾角下的低效识别与电量提升方法 | |
| CN116502516B (zh) | 一种航天器部件退化阶段的辨识方法和装置 | |
| CN117056847B (zh) | 一种流式数据的异常检测方法、系统、设备及存储介质 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |