一种卫星遥测数据的在轨异常监测方法技术

本发明专利技术公开一种卫星遥测数据的在轨异常监测方法，首先对被监测的遥测数据进行中心化处理，得到中心化的被监测遥测数据；然后，采用映射函数及重构函数，将中心化后的被监测的遥测数据进行投影和重构，得到重构的遥测数据；最后，根据重构的遥测数据及中心化的被监测遥测数据，计算重构误差，并与预设阈值进行比较判断被监测的数据是否异常。判断被监测的数据是否异常。判断被监测的数据是否异常。

【技术实现步骤摘要】
一种卫星遥测数据的在轨异常监测方法


[0001]本专利技术涉及航空航天
，特别涉及一种卫星遥测数据的在轨异常监测方法。

技术介绍

[0002]卫星在轨运行期间，其内部的运行状态监测系统获得的传感器参量信息编码后，通过遥测系统传输至地面，此遥测数据是地面了解卫星在轨运行状态的唯一依据。遥测数据量大、维度高、关系复杂、相关性及专业性强，反映了卫星所处的轨道信息、性能变化、工作模式切换以及是否出现故障等。对遥测数据的有效分析和智能计算，能够协助地面运管人员判断卫星的性能，为开展各种运行和维护管理工作提供有效依据。尤其是遥测数据中的异常数据，其不同于正常遥测数据的变化规律或不符合航天器工作模式设定，从而可以反映采集设备失效、传输链路受损、对应设备的性能退化、质量问题、机械及电子故障，或者设计不足等问题。及时而有效地发现遥测数据中存在的异常模式，并对被测设备进行远程指令修复、传输链路修复、软件故障排除、故障预案制定，或者维护服务优化等操作，对于提升地面服务质量，增强航天器设计、研制、生产以及维护各个环节的成熟度、安全性及可靠性具有显著的现实意义。

技术实现思路

[0003]针对现有技术中的部分或全部问题，本专利技术提供一种卫星遥测数据的在轨异常监测方法，包括：
[0004]采用正常历史遥测数据的均值，对被检测的遥测数据进行中心化处理；
[0005]采用映射函数及重构函数，将中心化后的被监测的遥测数据进行投影和重构；以及
[0006]计算重构误差，判断所述被检测的数据是否异常。/>[0007]进一步地，所述映射函数及重构函数通过正常历史遥测数据得到。
[0008]进一步地，所述映射函数及重构函数的学习包括：
[0009]根据中心化的正常历史遥测数据，计算正常遥测数据的协方差矩阵；
[0010]对所述协方差矩阵进行特征值分解，得到特征值及其对应的特征向量；以及
[0011]确定降维重构的准确度的阈值，计算降维后的维度个数k，将特征值按值的大小逐序排序，选前面的k个特征值所相应的特征向量，并按列组成特征矩阵，从而得到映射函数及重构函数。
[0012]本专利技术提供的一种卫星遥测数据的在轨异常监测方法，通过主成分分析(PCA)降维技术，通过学习卫星正常遥测数据，得到高维映射到低维的映射函数以及低维映射回高维的重构函数。在监测过程中，将学习到的映射函数和重构函数应用到被监测的遥测数据上，对被监测的遥测数据进行投影和重构，得到重构后的遥测数据，然后计算重构误差，进而判断所述被监测的遥测数据是否异常。所述方法算法简单，准确率高，可广泛应用于卫星
的在轨异常监测。
附图说明
[0013]为进一步阐明本专利技术的各实施例的以上和其它优点和特征，将参考附图来呈现本专利技术的各实施例的更具体的描述。可以理解，这些附图只描绘本专利技术的典型实施例，因此将不被认为是对其范围的限制。在附图中，为了清楚明了，相同或相应的部件将用相同或类似的标记表示。
[0014]图1示出本专利技术一个实施例的一种卫星遥测数据的在轨异常监测方法的流程示意图。
具体实施方式
[0015]以下的描述中，参考各实施例对本专利技术进行描述。然而，本领域的技术人员将认识到可在没有一个或多个特定细节的情况下或者与其它替换和/或附加方法、材料或组件一起实施各实施例。在其它情形中，未示出或未详细描述公知的结构、材料或操作以免模糊本专利技术的专利技术点。类似地，为了解释的目的，阐述了特定数量、材料和配置，以便提供对本专利技术的实施例的全面理解。然而，本专利技术并不限于这些特定细节。此外，应理解附图中示出的各实施例是说明性表示且不一定按正确比例绘制。
[0016]在本说明书中，对“一个实施例”或“该实施例”的引用意味着结合该实施例描述的特定特征、结构或特性被包括在本专利技术的至少一个实施例中。在本说明书各处中出现的短语“在一个实施例中”并不一定全部指代同一实施例。
[0017]需要说明的是，本专利技术的实施例以特定顺序对工艺步骤进行描述，然而这只是为了阐述该具体实施例，而不是限定各步骤的先后顺序。相反，在本专利技术的不同实施例中，可根据工艺的调节来调整各步骤的先后顺序。
[0018]针对卫星的在轨异常监测需求，本专利技术提供一种卫星遥测数据的在轨异常监测方法。基于主成分分析(PCA)降维技术的异常监测算法是一类无监督算法，其包括学习和监测两个过程。在学习过程中它通过学习卫星正常遥测数据，得到高维映射到低维的映射函数f以及低维映射回高维的重构函数g；以及在监测过程中，将学习到的映射函数f和重构函数g应用到被监测的遥测数据上，对被监测的遥测数据x进行投影和重构，得到重构后的遥测数据然后计算重构误差。如果被监测的遥测数据是正常的，那么重构数据应当接近原始遥测数据，即重构误差接近0；反之，如果被监测的遥测数据是异常的，那么重构数据与原始遥测数据差别会很大，即重构误差比较大。每个遥测数据的重构误差的计算公式为：
[0019][0020]r
i
的大小反映了被监测的遥测数据与学习的正常遥测数据之间的不同程度，当r
i
超过一定值的时候，则说明被监测的遥测数据是异常的。下面结合实施例附图，对本专利技术的方案作进一步描述。
[0021]图1示出本专利技术一个实施例的一种卫星遥测数据的在轨异常监测方法的流程示意图。如图1所示，一种卫星遥测数据的在轨异常监测方法，包括：
[0022]首先，在步骤101，数据中心化处理。采用正常历史遥测数据的均值，对被监测的遥测数据进行中心化处理，得到中心化的被监测遥测数据Y。在PCA中，数据在多维空间中处理，在多维空间中，每个遥测参数代表一个维度，同一时刻的多个遥测参数构成了多维空间中的一个遥测数据点，多个时刻的多个遥测参数构成了多维空间中一系列的遥测数据点。在对卫星遥测数据做PCA降维重构的过程中，每个遥测数据点的向量结构被看作一个行向量。由于在PCA的原理推导中涉及到协方差的计算，而计算协方差需要对数据样本进行中心化处理，在PCA中，对数据进行中心化处理，也可称为去均值，通过如下公式实现：
[0023][0024]其中，x
i
为多维遥测数据点的第i维度上的遥测参数，m为学习的正常多维遥测数据点的个数，以及x
i0
为多维遥测数据点第i维度上中心化后的遥测参数；
[0025]因此，所述中心化的被监测遥测数据Y＝{y
10
，y
20
，
…
，y
m0
}，其中，y
i
为被监测的遥测数据的第i维度上的遥测参数，m为遥测数据的个数；
[0026]接下来，在步骤102，数据投影及重构。采用映射函数f及重构函数g，将中心化后的被监测的遥测数据进行投影和重构。在本专利技术的一个实施例中，所述映射函数f及重构函数g通过正常历史遥测数据得到，包括：
[0027]首先，根据中心化的正常历史遥测数据X0，计算正常遥测数据的协方差矩阵：...

【技术保护点】

【技术特征摘要】
1.一种卫星遥测数据的在轨异常监测方法，其特征在于，包括步骤：对被监测的遥测数据进行中心化处理，得到中心化的被监测遥测数据Y；采用映射函数及重构函数，将中心化后的被监测的遥测数据进行投影和重构，得到重构的遥测数据以及根据所述重构的遥测数据及中心化的被监测遥测数据，计算重构误差R，并将所述重构误差与预设阈值进行比较：若所述重构误差大于预设阈值，则所述被监测的数据异常。2.如权利要求1所述的方法，其特征在于，所述中心化的被监测遥测数据Y＝{y
10
，y
20
，
…
，y
m0
}，其中，y
i
为被监测的遥测数据的第i维度上的遥测参数,m为遥测数据的个数。3.如权利要求1所述的方法，...

【专利技术属性】
技术研发人员：杨琼，张军，冷佳醒，张鸽，王亚宾，曾繁彬，丁澍恺，李国通，林宝军，沈苑，蒋桂忠，刘迎春，王学良，陈志峰，龚文斌，任前义，李锴，
申请(专利权)人：上海微小卫星工程中心，
类型：发明
国别省市：