CN112529096A - 基于pca的多维航天器遥测数据的故障诊断方法 - Google Patents
基于pca的多维航天器遥测数据的故障诊断方法 Download PDFInfo
- Publication number
- CN112529096A CN112529096A CN202011534550.3A CN202011534550A CN112529096A CN 112529096 A CN112529096 A CN 112529096A CN 202011534550 A CN202011534550 A CN 202011534550A CN 112529096 A CN112529096 A CN 112529096A
- Authority
- CN
- China
- Prior art keywords
- data
- characteristic
- covariance matrix
- spacecraft
- dimension
- 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.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
- G06F18/2135—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on approximation criteria, e.g. principal component analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Abstract
基于PCA的多维航天器遥测数据的故障诊断方法,属于故障诊断领域。本发明为了解决现有采用人工监测并检测航天器的遥测数据,检测效率低的问题。对航天器的遥测数据进行归一化处理得到特征协方差矩阵,根据特征协方差矩阵得到特征协方差矩阵的特征值和特征向量;采用PCA方法对特征协方差矩阵的特征值和特征向量进行降维处理,得到降维数据集的特征值和特征向量;将测试数据向所述特征向量做投影,得到投影后的分量;将投影后的分量与降维数据集的特征值进行一一对比,若投影后的特征值中有一个值大于降维数据集的特征值,则测试数据异常,若投影后的特征值中每一个值均小于降维数据集的特征值,则测试数据正常。它用于检测数据故障。
Description
技术领域
本发明涉及利用航天器遥测数据进行故障诊断技术领域,特别涉及一种基于主成分分析(PCA)的多维航天器遥测数据的故障诊断方法。即应用PCA算法在航天器故障诊断问题,目标为减小航天器的故障率,属于故障诊断领域。
背景技术
航天器作为典型的安全攸关系统,在航天器设计之初已尽可能考虑各种因素,确保高可靠安全运行。随着航天装备需求越来越多样、设计复杂性越来越高,航天器的智能自主功能越来越强,航天器提早发现异常和故障要求越来越高,准确检测和诊断故障的难度也越来越大。随着航天事业的发展,当今发射的航天器数量越来越多,航天器的结构和功能越来越复杂,各种新技术的应用也使得航天器的故障呈现出增多的趋势。航天器故障诊断技术对于提高航天器产品可靠性,减少故障发生概率具有重要的意义。航天器虽然在设计中已充分考虑了各种因素,但仍有许多航天器会发生在轨故障,导致任务失败。航天器飞行任务的扩展以及系统安全可靠性要求的提高促进了航天器故障诊断技术的发展。
航天器的遥测数据量大,维数多,在对其监测过程中,会对其进行检测耗费大量的人力物力。
发明内容
本发明的目的是为了解决现有采用人工监测并检测航天器的遥测数据,检测效率低的问题。现提供基于PCA的多维航天器遥测数据的故障诊断方法。
基于PCA的多维航天器遥测数据的故障诊断方法,所述方法包括以下步骤:
步骤1、对航天器的遥测数据进行归一化处理,得到特征协方差矩阵,根据特征协方差矩阵得到特征协方差矩阵的特征值和特征向量;
步骤2、采用PCA方法对特征协方差矩阵的特征值和特征向量进行降维处理,得到能够反映特征协方差矩阵预设特征范围的数据集,作为降维数据集的特征值和特征向量;
步骤3、将测试数据向所述特征向量做投影,得到投影后的分量;
步骤4、将投影后的分量与降维数据集的特征值进行一一对比,若投影后的特征值中有一个值大于降维数据集的特征值,则测试数据异常,若投影后的特征值中每一个值均小于降维数据集的特征值,则测试数据正常。
优选地,步骤1中,得到特征协方差矩阵的具体过程为:
步骤11、设定航天器的遥测数据的维度n和航天器遥测数据的样本总量m,获取每个维度数据的平均值;
步骤12、采用遥测数据中每个维度下的数据减去对应维度数据的平均值,得到中心化数据;
步骤13、根据中心化数据,得到特征协方差矩阵。
优选地,步骤11中,每个维度数据的平均值,表示为:
式中,i为航天器遥测数据的第i个样本,xij为航天器遥测数据的第i个样本的第j维数据,μj为第j维数据的平均值。
优选地,步骤12中,中心化数据表示为:
优选地,步骤13中,特征协方差矩阵C表示为:
优选地,步骤1中,根据特征协方差矩阵得到特征协方差矩阵的特征值和特征向量,具体为:
特征协方差矩阵的特征值表示为:
Cv=λv 公式4,
式中,λ是特征向量v对应的特征值;
特征协方差矩阵的特征向量表示为:
C=QΣQ-1 公式5,
式中,Q是矩阵C的特征向量组成的矩阵,Σ则是一个对角阵,对角线上的元素就是特征值。
优选地,预设特征范围为降维数据集的特征值和特征向量占特征协方差矩阵的90%。
本发明的有益效果为:
本申请基于PCA的多维航天器遥测数据的故障诊断方法,能够及时有效的监测航天器的运行情况,并且诊断出故障,对于保证航天器的可靠性和降低航天器的故障有较为明显的作用。该方法能够高效的识别出遥测数据的异常,更便于监控航天器的状态。
本申请的优点:
1、根据实际情况对航天器的遥测数据进行主成分分析,减少数据维数,使问题分析更为简单明了。
2、通过主成分分析的到数据的特征向量,在面对新的数据(测试数据)时可以对其是否发生故障进行合理判断。
附图说明
图1为基于PCA的多维航天器遥测数据的故障诊断方法的流程图;
图2为主成分分析方法(PCA)的原理图;
图3为测试数据在主成分方向(特征协方差矩阵的特征向量投影后)的分量图;
图4为测试数据是否为故障数据的判断图;
具体实施方式
具体实施方式一:结合图1至图3说明本实施方式,本实施方式所述的基于PCA的多维航天器遥测数据的故障诊断方法,所述方法包括以下步骤:
步骤1、对航天器的遥测数据进行归一化处理,得到特征协方差矩阵,根据特征协方差矩阵得到特征协方差矩阵的特征值和特征向量;
步骤2、采用PCA方法对特征协方差矩阵的特征值和特征向量进行降维处理,得到能够反映特征协方差矩阵预设特征范围的数据集,作为降维数据集的特征值和特征向量;
步骤3、将测试数据向所述特征向量做投影,得到投影后的分量;
步骤4、将投影后的分量与降维数据集的特征值进行一一对比,若投影后的特征值中有一个值大于降维数据集的特征值,则测试数据异常,若投影后的特征值中每一个值均小于降维数据集的特征值,则测试数据正常。
本实施方式中,采用主成分分析方法对特征协方差矩阵的特征值和特征向量进行降维处理,具体过程为:
对特征协方差矩阵的特征值和特征向量(主成分)进行主成分分析,PCA实质上是一个基变换,使得变换后的数据有最大的方差,也就是通过对坐标轴的旋转和坐标原点的平移使得其中一个轴(主轴)与特征协方差矩阵之间的方差最小,坐标转换后去掉高方差的正交轴,得到降维数据集。
本申请中比如特征协方差矩阵为20维空间矩阵,要将20维空间矩阵采用主成分分析方法转换成4维空间矩阵(降维数据集),具体转换成几维空间矩阵,是根据预设特征范围决定的,比如预设特征范围为90%。我们得到的4维空间矩阵比上20维空间矩阵,等于90%,那么就转换成4维空间矩阵。
此时4维空间矩阵在很大程度上反映原来n个变量的影响。即4维空间矩阵可以90%反映原有的n个变量,PCA的目标是寻找k(k<n)个新变量(k维空间矩阵),使它们反映事物的主要特征,压缩原有数据矩阵的规模,将特征向量的维数降低,挑选出最少的维数来概括最重要特征。每个新变量是原有变量的线性组合,体现原有变量的综合效果,具有一定的实际含义。这k个新变量称为“主成分”,它们可以在很大程度上反映原来n个变量的影响,并且这些新变量是互不相关的,也是正交的。通过主成分分析,压缩数据空间,将多元数据的特征在低维空间里直观地表示出来。
具体实施方式二:本实施方式是对具体实施方式一所述的基于PCA的多维航天器遥测数据的故障诊断方法,步骤1中,得到特征协方差矩阵的具体过程为:
步骤11、设定航天器的遥测数据的维度n和航天器遥测数据的样本总量m,获取每个维度数据的平均值;
步骤12、采用遥测数据中每个维度下的数据减去对应维度数据的平均值,得到中心化数据;
步骤13、根据中心化数据,得到特征协方差矩阵。
具体实施方式三:本实施方式是对具体实施方式二所述的基于PCA的多维航天器遥测数据的故障诊断方法,步骤11中,每个维度数据的平均值,表示为:
式中,i为航天器遥测数据的第i个样本,xij为航天器遥测数据的第i个样本的第j维数据,μj为第j维数据的平均值。
具体实施方式四:本实施方式是对具体实施方式三所述的基于PCA的多维航天器遥测数据的故障诊断方法,步骤12中,中心化数据表示为:
具体实施方式五:本实施方式是对具体实施方式四所述的基于PCA的多维航天器遥测数据的故障诊断方法,步骤13中,特征协方差矩阵C表示为:
具体实施方式六:本实施方式是对具体实施方式五所述的基于PCA的多维航天器遥测数据的故障诊断方法,步骤1中,根据特征协方差矩阵得到特征协方差矩阵的特征值和特征向量,具体为:
特征协方差矩阵的特征值表示为:
Cv=λv 公式4,
式中,λ是特征向量v对应的特征值;
特征协方差矩阵的特征向量表示为:
C=QΣQ-1 公式5,
式中,Q是矩阵C的特征向量组成的矩阵,Σ则是一个对角阵,对角线上的元素就是特征值。
本实施方式中,对于矩阵C,有一组特征向量v,将这组向量进行正交化单位化,就能得到一组正交单位向量。
具体实施方式七:本实施方式是对具体实施方式五所述的基于PCA的多维航天器遥测数据的故障诊断方法,预设特征范围为降维数据集的特征值和特征向量占特征协方差矩阵的90%。
Claims (7)
1.基于PCA的多维航天器遥测数据的故障诊断方法,其特征在于,所述方法包括以下步骤:
步骤1、对航天器的遥测数据进行归一化处理,得到特征协方差矩阵,根据特征协方差矩阵得到特征协方差矩阵的特征值和特征向量;
步骤2、采用PCA方法对特征协方差矩阵的特征值和特征向量进行降维处理,得到能够反映特征协方差矩阵预设特征范围的数据集,作为降维数据集的特征值和特征向量;
步骤3、将测试数据向所述特征向量做投影,得到投影后的分量;
步骤4、将投影后的分量与降维数据集的特征值进行一一对比,若投影后的特征值中有一个值大于降维数据集的特征值,则测试数据异常,若投影后的特征值中每一个值均小于降维数据集的特征值,则测试数据正常。
2.根据权利要求1所述的基于PCA的多维航天器遥测数据的故障诊断方法,其特征在于,步骤1中,得到特征协方差矩阵的具体过程为:
步骤11、设定航天器的遥测数据的维度n和航天器遥测数据的样本总量m,获取每个维度数据的平均值;
步骤12、采用遥测数据中每个维度下的数据减去对应维度数据的平均值,得到中心化数据;
步骤13、根据中心化数据,得到特征协方差矩阵。
6.根据权利要求5所述的基于PCA的多维航天器遥测数据的故障诊断方法,其特征在于,步骤1中,根据特征协方差矩阵得到特征协方差矩阵的特征值和特征向量,具体为:
特征协方差矩阵的特征值表示为:
Cv=λv 公式4,
式中,λ是特征向量v对应的特征值;
特征协方差矩阵的特征向量表示为:
C=QΣQ-1 公式5,
式中,Q是矩阵C的特征向量组成的矩阵,Σ则是一个对角阵,对角线上的元素就是特征值。
7.根据权利要求5所述的基于PCA的多维航天器遥测数据的故障诊断方法,其特征在于,预设特征范围为降维数据集的特征值和特征向量占特征协方差矩阵的90%。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011534550.3A CN112529096A (zh) | 2020-12-22 | 2020-12-22 | 基于pca的多维航天器遥测数据的故障诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011534550.3A CN112529096A (zh) | 2020-12-22 | 2020-12-22 | 基于pca的多维航天器遥测数据的故障诊断方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112529096A true CN112529096A (zh) | 2021-03-19 |
Family
ID=74975843
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011534550.3A Pending CN112529096A (zh) | 2020-12-22 | 2020-12-22 | 基于pca的多维航天器遥测数据的故障诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112529096A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113496261A (zh) * | 2021-07-06 | 2021-10-12 | 中国科学院微小卫星创新研究院 | 一种卫星遥测数据的在轨异常监测方法 |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008079695A1 (en) * | 2006-12-21 | 2008-07-03 | Honeywell International Inc. | Monitoring and fault detection in dynamic systems |
CN103413050A (zh) * | 2013-08-20 | 2013-11-27 | 北京工业大学 | 基于极速学习机的运动想象脑电信号投票策略分类方法 |
CN103471865A (zh) * | 2013-09-12 | 2013-12-25 | 北京交通大学 | 基于线性判别法的列车悬挂系统故障分离方法 |
CN105590026A (zh) * | 2015-12-16 | 2016-05-18 | 上海卫星工程研究所 | 基于主成分分析的卫星遥测回归方法 |
CN106546882A (zh) * | 2016-09-23 | 2017-03-29 | 国网江苏省电力公司南京供电公司 | 一种检测电力变压器内部放电故障的方法 |
CN108052954A (zh) * | 2017-11-01 | 2018-05-18 | 佛山科学技术学院 | 基于多级高维特征的样本空间的故障诊断方法 |
CN108062565A (zh) * | 2017-12-12 | 2018-05-22 | 重庆科技学院 | 基于化工te过程的双主元-动态核主元分析故障诊断方法 |
CN108829978A (zh) * | 2018-06-20 | 2018-11-16 | 哈尔滨工业大学 | 基于统计特征评估与高斯过程回归的直流电源在线异常检测方法 |
CN109711461A (zh) * | 2018-12-25 | 2019-05-03 | 中国人民解放军战略支援部队航天工程大学 | 基于主成分分析的迁移学习图片分类方法及其装置 |
CN111274543A (zh) * | 2020-01-17 | 2020-06-12 | 北京空间飞行器总体设计部 | 一种基于高维空间映射的航天器系统异常检测方法 |
CN111678679A (zh) * | 2020-05-06 | 2020-09-18 | 内蒙古电力(集团)有限责任公司电力调度控制分公司 | 一种基于pca-bpnn的断路器故障诊断方法 |
CN111738309A (zh) * | 2020-06-03 | 2020-10-02 | 哈尔滨工业大学 | 多尺度分析和集成学习的气敏传感器故障模式识别方法 |
-
2020
- 2020-12-22 CN CN202011534550.3A patent/CN112529096A/zh active Pending
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008079695A1 (en) * | 2006-12-21 | 2008-07-03 | Honeywell International Inc. | Monitoring and fault detection in dynamic systems |
CN103413050A (zh) * | 2013-08-20 | 2013-11-27 | 北京工业大学 | 基于极速学习机的运动想象脑电信号投票策略分类方法 |
CN103471865A (zh) * | 2013-09-12 | 2013-12-25 | 北京交通大学 | 基于线性判别法的列车悬挂系统故障分离方法 |
CN105590026A (zh) * | 2015-12-16 | 2016-05-18 | 上海卫星工程研究所 | 基于主成分分析的卫星遥测回归方法 |
CN106546882A (zh) * | 2016-09-23 | 2017-03-29 | 国网江苏省电力公司南京供电公司 | 一种检测电力变压器内部放电故障的方法 |
CN108052954A (zh) * | 2017-11-01 | 2018-05-18 | 佛山科学技术学院 | 基于多级高维特征的样本空间的故障诊断方法 |
CN108062565A (zh) * | 2017-12-12 | 2018-05-22 | 重庆科技学院 | 基于化工te过程的双主元-动态核主元分析故障诊断方法 |
CN108829978A (zh) * | 2018-06-20 | 2018-11-16 | 哈尔滨工业大学 | 基于统计特征评估与高斯过程回归的直流电源在线异常检测方法 |
CN109711461A (zh) * | 2018-12-25 | 2019-05-03 | 中国人民解放军战略支援部队航天工程大学 | 基于主成分分析的迁移学习图片分类方法及其装置 |
CN111274543A (zh) * | 2020-01-17 | 2020-06-12 | 北京空间飞行器总体设计部 | 一种基于高维空间映射的航天器系统异常检测方法 |
CN111678679A (zh) * | 2020-05-06 | 2020-09-18 | 内蒙古电力(集团)有限责任公司电力调度控制分公司 | 一种基于pca-bpnn的断路器故障诊断方法 |
CN111738309A (zh) * | 2020-06-03 | 2020-10-02 | 哈尔滨工业大学 | 多尺度分析和集成学习的气敏传感器故障模式识别方法 |
Non-Patent Citations (4)
Title |
---|
FENG BINGQING等: "Anomaly Detection of Spacecraft Attitude Control System Based on Principal Component Analysis", 《2017 29TH CHINESE CONTROL AND DECISION CONFERENCE(CCDC)》 * |
刘敏: "数据驱动的卫星姿态控制系统微小故障检测与预测方法研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 * |
柴敏等: "面向故障诊断的航天器遥测数据降维分析技术", 《弹箭与制导学报》 * |
许力等: "基于变尺度PCA的电力设备载流故障早期预警", 《电力自动化设备》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113496261A (zh) * | 2021-07-06 | 2021-10-12 | 中国科学院微小卫星创新研究院 | 一种卫星遥测数据的在轨异常监测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20220146462A1 (en) | Method for predicting remaining life of numerical control machine tool | |
Chen et al. | Probability-relevant incipient fault detection and diagnosis methodology with applications to electric drive systems | |
CN107505133B (zh) | 基于自适应mrvm的滚动轴承故障概率性智能诊断方法 | |
US7151976B2 (en) | Multivariate control of semiconductor processes | |
US7243048B2 (en) | Fault detection system and method using multiway principal component analysis | |
Kwan et al. | A novel approach to fault diagnostics and prognostics | |
US7313454B2 (en) | Method and apparatus for classifying manufacturing outputs | |
US20090132626A1 (en) | Method and system for detecting difference between plural observed results | |
CN111274543A (zh) | 一种基于高维空间映射的航天器系统异常检测方法 | |
Yin et al. | On PCA-based fault diagnosis techniques | |
CN112529096A (zh) | 基于pca的多维航天器遥测数据的故障诊断方法 | |
CN108052954B (zh) | 基于多级高维特征的样本空间的故障诊断方法 | |
CN113496261A (zh) | 一种卫星遥测数据的在轨异常监测方法 | |
CN110221590B (zh) | 一种基于判别分析的工业过程多故障诊断方法 | |
Kallas et al. | Fault detection and estimation using kernel principal component analysis | |
Jaffel et al. | Online process monitoring using a new PCMD index | |
CN109213120A (zh) | 基于高维特征表示下多级主元空间的故障诊断方法 | |
CN114611606A (zh) | 基于核混合空间投影的故障检测方法 | |
CN117031294A (zh) | 一种电池多故障检测方法、装置及存储介质 | |
Rahul et al. | A statistical approach for failure analysis involving uncertainty in determining ply orientation | |
Zhang et al. | Improved incipient fault detection using Jensen-Shannon divergence and KPCA | |
Zhang et al. | A novel monitoring strategy combining the advantages of NPE and GMM | |
Nassar et al. | Supervised learning algorithms for spacecraft attitude determination and control system health monitoring | |
Yan et al. | Condition monitoring of the aircraft structural fatigue test based on PCA algorithm | |
Márquez-Vera et al. | Adaptive threshold PCA for fault detection and isolation |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210319 |
|
RJ01 | Rejection of invention patent application after publication |