[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN111241158A - Anomaly detection method and device for aircraft telemetry data - Google Patents

Anomaly detection method and device for aircraft telemetry data Download PDF

Info

Publication number
CN111241158A
CN111241158A CN202010016293.8A CN202010016293A CN111241158A CN 111241158 A CN111241158 A CN 111241158A CN 202010016293 A CN202010016293 A CN 202010016293A CN 111241158 A CN111241158 A CN 111241158A
Authority
CN
China
Prior art keywords
node
data
detection algorithm
current moment
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
Application number
CN202010016293.8A
Other languages
Chinese (zh)
Other versions
CN111241158B (en
Inventor
詹亚锋
万鹏
曾冠铭
陈曦
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tsinghua University
Original Assignee
Tsinghua University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tsinghua University filed Critical Tsinghua University
Priority to CN202010016293.8A priority Critical patent/CN111241158B/en
Publication of CN111241158A publication Critical patent/CN111241158A/en
Application granted granted Critical
Publication of CN111241158B publication Critical patent/CN111241158B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/24Querying
    • G06F16/245Query processing
    • G06F16/2458Special types of queries, e.g. statistical queries, fuzzy queries or distributed queries
    • G06F16/2474Sequence data queries, e.g. querying versioned data
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/21Design, administration or maintenance of databases
    • G06F16/215Improving data quality; Data cleansing, e.g. de-duplication, removing invalid entries or correcting typographical errors
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/28Databases characterised by their database models, e.g. relational or object models
    • G06F16/284Relational databases
    • G06F16/285Clustering or classification

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Databases & Information Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Quality & Reliability (AREA)
  • Fuzzy Systems (AREA)
  • Mathematical Physics (AREA)
  • Probability & Statistics with Applications (AREA)
  • Software Systems (AREA)
  • Computational Linguistics (AREA)
  • Selective Calling Equipment (AREA)
  • Arrangements For Transmission Of Measured Signals (AREA)

Abstract

The invention provides an anomaly detection method and device for aircraft telemetering data, which relate to the technical field of data processing and comprise the steps of obtaining pilot telemetering data sent by an aircraft, and determining whether the pilot telemetering data is a stable time sequence or not through stationarity detection; after stationarity detection is finished, acquiring telemetering data of a node at the current time, which is sent by an aircraft; if the pilot telemetering data is a stable time sequence, detecting the telemetering data of the node at the current moment by using a boundary detection algorithm, and determining whether the telemetering data of the node at the current moment is abnormal data; if the pilot telemetering data is not the stationary time sequence, detecting the telemetering data of the node at the current moment by using a preset detection algorithm, and determining whether the telemetering data of the node at the current moment is abnormal data, so that the technical problems that the detection step for detecting the abnormal data of the telemetering data is complex and the detection accuracy is low in the prior art are solved.

Description

一种飞行器遥测数据的异常检测方法和装置Anomaly detection method and device for aircraft telemetry data

技术领域technical field

本发明涉及数据处理技术领域,尤其是涉及一种飞行器遥测数据的异常检测方法和装置。The present invention relates to the technical field of data processing, in particular to a method and device for detecting abnormality of telemetry data of an aircraft.

背景技术Background technique

对于卫星遥测而言,空间环境瞬态剧烈变化以及设备线路噪声将会给信号采集带来噪声干扰,从而使得遥测数据中掺杂随机分布的异常数据,异常数据对深空探测任务及空间信息网络的发展带来了巨大的困难。For satellite telemetry, the rapid transient changes in the space environment and equipment line noise will cause noise interference to signal acquisition, so that the telemetry data is mixed with randomly distributed abnormal data. development has brought great difficulties.

现有技术中,一般采用边界检测算法,趋势预测算法,速率约束算法,密度(距离)检测算法等算法对遥测数据进行异常检测,但是上述的方法检测步骤复杂,且检测准确率较低。In the prior art, algorithms such as boundary detection algorithm, trend prediction algorithm, rate constraint algorithm, density (distance) detection algorithm, etc. are generally used to detect abnormality of telemetry data, but the above-mentioned methods have complicated detection steps and low detection accuracy.

针对上述问题,还未提出有效的解决方案。For the above problems, no effective solutions have been proposed yet.

发明内容SUMMARY OF THE INVENTION

有鉴于此,本发明的目的在于提供一种飞行器遥测数据的异常检测方法和装置,以缓解了现有技术中对遥测数据的进行异常数据检测的检测步骤复杂,且检测准确率较低的技术问题。In view of this, the purpose of the present invention is to provide an abnormality detection method and device for aircraft telemetry data, so as to alleviate the complicated detection steps and low detection accuracy in the prior art for abnormal data detection of telemetry data. question.

第一方面,本发明实施例提供了一种飞行器遥测数据的异常检测方法,包括:获取飞行器发送的先导遥测数据,并通过平稳性检测确定所述先导遥测数据是否为平稳时间序列;在完成平稳性检测之后,获取飞行器发送的当前时刻节点的遥测数据;若所述先导遥测数据为平稳时间序列,则利用边界检测算法对当前时刻节点的遥测数据进行检测,确定所述当前时刻节点的遥测数据是否为异常数据;若所述先导遥测数据不是平稳时间序列,则利用预设检测算法对所述当前时刻节点的遥测数据进行检测,确定所述当前时刻节点的遥测数据是否为异常数据,其中,所述预设检测算法包括以下至少之一:三点集同构映射检测算法,左右双陪集赋权映射检测算法。In a first aspect, an embodiment of the present invention provides an abnormality detection method for aircraft telemetry data, including: acquiring pilot telemetry data sent by an aircraft, and determining whether the pilot telemetry data is a stationary time series through stationarity detection; After performance detection, obtain the telemetry data of the node at the current moment sent by the aircraft; if the pilot telemetry data is a stationary time series, use the boundary detection algorithm to detect the telemetry data of the node at the current moment, and determine the telemetry data of the node at the current moment Whether it is abnormal data; if the pilot telemetry data is not a stationary time series, a preset detection algorithm is used to detect the telemetry data of the node at the current moment to determine whether the telemetry data of the node at the current moment is abnormal data, wherein, The preset detection algorithm includes at least one of the following: a three-point set isomorphic map detection algorithm, and a left and right dual coset weighted map detection algorithm.

进一步地,所述方法还包括:对所述异常数据进行清洗。Further, the method further includes: cleaning the abnormal data.

进一步地,对所述异常数据进行清洗,包括:确定出第一目标数据和第二目标数据,其中,所述第一目标数据为所述当前时刻节点之前的第一个正常遥测数据,所述第二目标数据为所述当前时刻节点之后的第一个正常遥测数据;计算所述第一目标数据和所述第二目标数据的均值,并将所述异常数据替换为所述均值。Further, cleaning the abnormal data includes: determining first target data and second target data, wherein the first target data is the first normal telemetry data before the node at the current moment, and the The second target data is the first normal telemetry data after the node at the current moment; the average value of the first target data and the second target data is calculated, and the abnormal data is replaced with the average value.

进一步地,利用预设检测算法对所述当前时刻节点的遥测数据进行检测,确定所述当前时刻节点的遥测数据是否为异常数据,包括:确定所述飞行器的性能类型,其中,所述性能类型包括:第一性能类型和第二性能类型,所述第一性能类型的飞行器的在轨处理能力低于所述第二性能类型的飞行器的在轨处理能力;基于所述性能类型,确定出所述预设检测算法中的目标检测算法;利用所述目标检测算法,对所述当前时刻节点的遥测数据进行检测,确定所述当前时刻节点的遥测数据是否为异常数据。Further, using a preset detection algorithm to detect the telemetry data of the node at the current moment, and determining whether the telemetry data of the node at the current moment is abnormal data, includes: determining the performance type of the aircraft, wherein the performance type Including: a first performance type and a second performance type, the on-orbit processing capability of the aircraft of the first performance type is lower than the on-orbit processing capability of the aircraft of the second performance type; based on the performance type, determine the The target detection algorithm in the preset detection algorithm is used; the target detection algorithm is used to detect the telemetry data of the node at the current moment, and determine whether the telemetry data of the node at the current moment is abnormal data.

进一步地,基于所述性能类型,确定出所述预设检测算法中的目标检测算法,包括:若所述飞行器的性能类型为第一性能类型,则所述目标检测算法为所述三点集同构映射检测算法;若所述飞行器的性能类型为第二性能类型,则所述目标检测算法为所述左右双陪集赋权映射检测算法。Further, determining the target detection algorithm in the preset detection algorithm based on the performance type, including: if the performance type of the aircraft is the first performance type, the target detection algorithm is the three-point set isomorphic mapping detection algorithm; if the performance type of the aircraft is the second performance type, the target detection algorithm is the left and right dual coset weighted mapping detection algorithm.

进一步地,若所述目标检测算法为所述三点集同构映射检测算法;利用所述目标检测算法,对所述当前时刻节点的遥测数据进行检测,确定所述当前时刻节点的遥测数据是否为异常数据,包括:确定出所述当前时刻节点的邻居节点,其中,所述邻居节点包括:第一邻居节点和第二邻居节点,所述第一邻居节点为所述当前时刻节点之前的邻居节点,所述第二邻居节点为所述当前时刻节点之后的邻居节点;利用所述当前时刻节点和所述第一邻居节点,构建第一矢量线段,以及利用所述当前时刻节点所述第二邻居节点构建的第二矢量线段,并确定出所述第一矢量线段与所述第二矢量线段构成的矢量夹角的角度;若所述角度小于预设阈值,则所述当前时刻节点的遥测数据为异常数据。Further, if the target detection algorithm is the three-point set isomorphic mapping detection algorithm; use the target detection algorithm to detect the telemetry data of the node at the current moment to determine whether the telemetry data of the node at the current moment is not. is abnormal data, including: determining the neighbor nodes of the node at the current moment, wherein the neighbor nodes include: a first neighbor node and a second neighbor node, and the first neighbor node is the neighbor before the node at the current moment node, the second neighbor node is the neighbor node after the current moment node; use the current moment node and the first neighbor node to construct a first vector line segment, and use the current moment node to construct the second The second vector line segment constructed by the neighbor node, and determine the angle of the vector angle formed by the first vector line segment and the second vector line segment; if the angle is less than the preset threshold, the telemetry of the node at the current moment The data is abnormal data.

进一步地,若所述目标检测算法为所述左右双陪集赋权映射检测算法;利用所述目标检测算法,对所述遥测数据进行检测,确定出所述遥测数据中的异常数据,包括:构建所述当前时刻节点的k元双陪集;计算所述k元双陪集的距离阈值,以及为所述k元双陪集中的节点赋值,得到节点赋值;确定出所述当前时刻节点的紧密因子比判阈值;结合所述k元双陪集的距离阈值和所述节点赋值,计算出所述当前时刻节点的紧密度;Further, if the target detection algorithm is the left and right bi-coset weighted mapping detection algorithm; using the target detection algorithm to detect the telemetry data to determine abnormal data in the telemetry data, including: Construct the k-ary double coset of the node at the current moment; calculate the distance threshold of the k-ary double coset, and assign values to the nodes in the k-ary double coset, and obtain the node assignment; determine the current moment of the node. Tightness factor comparison threshold; Combining the distance threshold of the k-ary double coset and the node assignment, calculate the tightness of the node at the current moment;

若所述紧密度小于紧密因子比判阈值,则所述当前时刻节点的遥测数据为异常数据。If the closeness is less than the closeness factor comparison threshold, the telemetry data of the node at the current moment is abnormal data.

进一步地,利用边界检测算法对当前时刻节点的遥测数据进行检测,确定所述当前时刻节点的遥测数据是否为异常数据,包括:确定所述当前时刻节点的遥测数据的遥测参数是否处于预设边界范围内;若否,则所述当前时刻节点的遥测数据为异常数据。Further, using a boundary detection algorithm to detect the telemetry data of the node at the current moment, and determining whether the telemetry data of the node at the current moment is abnormal data, includes: determining whether the telemetry parameter of the telemetry data of the node at the current moment is in a preset boundary If not, the telemetry data of the node at the current moment is abnormal data.

第二方面,本发明实施例还提供了一种飞行器遥测数据的异常检测装置,包括:第一获取单元,第二获取单元,第一检测单元和第二检测单元,其中,所述第一获取单元用于获取飞行器发送的先导遥测数据,并通过平稳性检测确定所述先导遥测数据是否为平稳时间序列;所述第二获取单元用于在完成平稳性检测之后,获取飞行器发送的当前时刻节点的遥测数据;所述第一检测单元用于在所述先导遥测数据为平稳时间序列的情况下,利用边界检测算法对当前时刻节点的遥测数据进行检测,确定所述当前时刻节点的遥测数据是否为异常数据;所述第二检测单元用于在所述先导遥测数据不是平稳时间序列的情况下,则利用预设检测算法对所述当前时刻节点的遥测数据进行检测,确定所述当前时刻节点的遥测数据是否为异常数据,其中,所述预设检测算法包括以下至少之一:三点集同构映射检测算法,左右双陪集赋权映射检测算法。In a second aspect, an embodiment of the present invention further provides an abnormality detection device for aircraft telemetry data, including: a first acquisition unit, a second acquisition unit, a first detection unit, and a second detection unit, wherein the first acquisition unit The unit is used to acquire the pilot telemetry data sent by the aircraft, and determine whether the pilot telemetry data is a stationary time series through stationarity detection; the second acquisition unit is used to acquire the current moment node sent by the aircraft after the stationarity detection is completed The first detection unit is used to detect the telemetry data of the node at the current moment by using the boundary detection algorithm when the pilot telemetry data is a stationary time series, and determine whether the telemetry data of the node at the current moment is not. is abnormal data; the second detection unit is configured to use a preset detection algorithm to detect the telemetry data of the node at the current moment when the pilot telemetry data is not a stationary time series, and determine the node at the current moment Whether the telemetry data is abnormal data, the preset detection algorithm includes at least one of the following: a three-point set isomorphic mapping detection algorithm, and a left and right dual coset weighted mapping detection algorithm.

第三方面,本申请实施例还提供了一种具有处理器可执行的非易失的程序代码的计算机可读介质,所述程序代码使所述处理器执行第一方面所述的飞行器遥测数据的异常检测方法。In a third aspect, an embodiment of the present application further provides a computer-readable medium having a non-volatile program code executable by a processor, the program code causing the processor to execute the aircraft telemetry data described in the first aspect anomaly detection method.

在本发明实施例中,首先,获取飞行器发送的先导遥测数据,并通过平稳性检测确定先导遥测数据是否为平稳时间序列;然后,在完成平稳性检测之后,获取飞行器发送的当前时刻节点的遥测数据;如果先导遥测数据为平稳时间序列,则利用边界检测算法对当前时刻节点的遥测数据进行检测,确定当前时刻节点的遥测数据是否为异常数据;如果先导遥测数据不是平稳时间序列,则利用预设检测算法对当前时刻节点的遥测数据进行检测,确定当前时刻节点的遥测数据是否为异常数据,其中,预设检测算法包括以下至少之一:三点集同构映射检测算法,左右双陪集赋权映射检测算法。In the embodiment of the present invention, first, the pilot telemetry data sent by the aircraft is obtained, and the stationarity detection is used to determine whether the pilot telemetry data is a stationary time series; then, after the stationarity detection is completed, the telemetry of the node at the current moment sent by the aircraft is obtained If the pilot telemetry data is a stationary time series, the boundary detection algorithm is used to detect the telemetry data of the node at the current moment to determine whether the telemetry data of the node at the current moment is abnormal data; if the pilot telemetry data is not a stationary time series, the prediction Suppose the detection algorithm detects the telemetry data of the node at the current moment, and determines whether the telemetry data of the node at the current moment is abnormal data, wherein the preset detection algorithm includes at least one of the following: a three-point set isomorphic mapping detection algorithm, left and right double cosets Weighted map detection algorithm.

在本发明实施例中,通过对先导遥测数据进行平稳性检测,确定出先导遥测数据是否为平稳时间序列,并根据检测结果选择不同的检测算法对当前时刻节点的遥测数据进行检测,从而确定出当前时刻节点的遥测数据是否为异常数据,达到了对遥测数据进行异常检测的目的,进而解决了现有技术中对遥测数据的进行异常数据检测的检测步骤复杂,且检测准确率较低的技术问题,从而实现了简化了遥测数据的异常数据检测的检测步骤,提高了检测准确率的技术效果。In the embodiment of the present invention, by performing stationarity detection on the pilot telemetry data, it is determined whether the pilot telemetry data is a stationary time series, and different detection algorithms are selected according to the detection results to detect the telemetry data of the node at the current moment, so as to determine whether the pilot telemetry data is a stationary time series. Whether the telemetry data of the node at the current moment is abnormal data, the purpose of abnormality detection of telemetry data is achieved, and the technology of abnormal data detection of telemetry data in the prior art is solved with complicated detection steps and low detection accuracy. Therefore, the detection steps of abnormal data detection of telemetry data are simplified, and the technical effect of improving detection accuracy is realized.

本发明的其他特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。Other features and advantages of the present invention will be set forth in the description which follows, and in part will be apparent from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the description, claims and drawings.

为使本发明的上述目的、特征和优点能更明显易懂,下文特举较佳实施例,并配合所附附图,作详细说明如下。In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, preferred embodiments are given below, and are described in detail as follows in conjunction with the accompanying drawings.

附图说明Description of drawings

为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。In order to illustrate the specific embodiments of the present invention or the technical solutions in the prior art more clearly, the following briefly introduces the accompanying drawings that need to be used in the description of the specific embodiments or the prior art. Obviously, the accompanying drawings in the following description The drawings are some embodiments of the present invention. For those of ordinary skill in the art, other drawings can also be obtained based on these drawings without creative efforts.

图1为本发明实施例提供的一种飞行器遥测数据的异常检测方法的流程图;1 is a flowchart of a method for detecting anomalies in aircraft telemetry data according to an embodiment of the present invention;

图2为本发明实施例提供的另一种飞行器遥测数据的异常检测方法的流程图;FIG. 2 is a flowchart of another abnormal detection method for aircraft telemetry data provided by an embodiment of the present invention;

图3为本发明实施例提供的一种利用三点集同构映射检测算法检测异常数据的流程图;3 is a flowchart of detecting abnormal data by using a three-point set isomorphic mapping detection algorithm according to an embodiment of the present invention;

图4为本发明实施例提供的一种利用左右双陪集赋权映射检测算法检测异常数据的流程图;4 is a flowchart of detecting abnormal data by utilizing the left and right dual coset weighted mapping detection algorithm provided by an embodiment of the present invention;

图5为本发明实施例提供的一种飞行器遥测数据的异常检测装置的示意图。FIG. 5 is a schematic diagram of an abnormality detection apparatus for aircraft telemetry data according to an embodiment of the present invention.

具体实施方式Detailed ways

为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。In order to make the purposes, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are part of the embodiments of the present invention, but not all of them. example. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.

目前的异常检测方法包括边界检测算法、趋势预测算法、速率约束算法、密度(距离)检测算法等算法,其中:The current anomaly detection methods include boundary detection algorithm, trend prediction algorithm, rate constraint algorithm, density (distance) detection algorithm and other algorithms, among which:

1)边界检测算法:边界检测算法是卫星任务中应用最为广泛的异常检测方法,通常由飞行器控制人员根据地面测试阶段获得的样本以及以往卫星任务的经验,提供各个遥测参数工作的正常边界,当数据超过上下界范围则认为是异常数据;1) Boundary detection algorithm: Boundary detection algorithm is the most widely used anomaly detection method in satellite missions. Usually, the aircraft controller provides the normal boundaries of each telemetry parameter based on the samples obtained during the ground test phase and the experience of previous satellite missions. If the data exceeds the upper and lower bounds, it is considered to be abnormal data;

2)趋势预测算法:由于时间序列通常具有一定趋势,趋势预测算法通过估计时间序列的行为趋势,对后续数据行为进行预测并给出自适应检测边界,当数据超过预测范围则认为是异常数据。2) Trend prediction algorithm: Since the time series usually has a certain trend, the trend prediction algorithm predicts the subsequent data behavior and provides an adaptive detection boundary by estimating the behavior trend of the time series. When the data exceeds the prediction range, it is considered as abnormal data.

3)速率约束算法:由于时间序列具有一维时变特性,速率约束算法基于速率的清洗方法采用了前后点之间速度范围(即斜率)受限的思想,若速率过大就说明当前点具有异变特性,认为是异常数据。3) Rate constraint algorithm: Since the time series has one-dimensional time-varying characteristics, the rate-based cleaning method of the rate constraint algorithm adopts the idea that the speed range (that is, the slope) between the front and rear points is limited. If the rate is too large, it means that the current point has Mutation characteristics, considered abnormal data.

4)密度(距离)检测算法:基于密度(距离)的异常检测算法在多维大数据处理中应用广泛,即给定一个距离阈值,若该范围内某点的邻居数目少于给定检测阈值k,则认为是异常数据;换算为单位距离,该方法同构于在单位空间范围内用户密度检测。4) Density (distance) detection algorithm: Density (distance)-based anomaly detection algorithms are widely used in multi-dimensional big data processing, that is, given a distance threshold, if the number of neighbors of a point in the range is less than the given detection threshold k , it is considered to be abnormal data; converted to unit distance, this method is isomorphic to the detection of user density in the unit space.

但是,上述的算法存在以下缺点:However, the above algorithm has the following disadvantages:

1)边界检测算法:采用常用的边界检测方法,各参数检测边界固定,仅能发现少数极端突变数据,对于短期行为变化并不敏感,需研究适用于局部异常检测的算法,当短期遥测数据存在异常跳变,但是并未超出边界范围,因此,利用边界检测算法将无法检测出该异常跳变,导致出现存在漏检的情况。1) Boundary detection algorithm: The commonly used boundary detection method is adopted. The detection boundary of each parameter is fixed, and only a few extreme mutation data can be found. It is not sensitive to short-term behavior changes. It is necessary to study algorithms suitable for local anomaly detection. When short-term telemetry data exists Abnormal jump, but does not exceed the boundary range, therefore, the boundary detection algorithm will not be able to detect the abnormal jump, resulting in a situation of missed detection.

2)趋势检测算法:由于高实时性处理要求与受限的在轨资源不允许长期观测来获得统计规律,遥测时间序列往往具有非理想趋势,短期统计结果并不稳定,需研究一种基于少量数据的轻量级算法。2) Trend detection algorithm: Due to high real-time processing requirements and limited on-orbit resources, long-term observation is not allowed to obtain statistical laws, telemetry time series often have non-ideal trends, and short-term statistical results are not stable. Lightweight algorithms for data.

3)速率检测范围:当时间序列呈现出随机特性时,前后数据之间的速率变化无规律且速率跳变过大,速率检测无法正常工作,需要根据遥测行为类型进行特殊处理;此外,速率检测边界也是固定的,对于短期速率超出边界的正常数据存在误判的可能性,需研究适用于局部异常检测的算法。3) Rate detection range: When the time series exhibits random characteristics, the rate changes between the data before and after are irregular and the rate jump is too large, the rate detection cannot work normally, and special processing is required according to the type of telemetry behavior; in addition, the rate detection The boundary is also fixed. For normal data whose short-term rate exceeds the boundary, there is a possibility of misjudgment. It is necessary to study an algorithm suitable for local anomaly detection.

4)密度(距离)检测算法:传统的基于密度(距离)的算法并未考虑前后序列权重对于评估当前点紧密程度的影响,对于时变信源而言,虽然前后数据之间具有一定的相关性,但是这种相关性会随着间隔时间的增长而越来越弱。若忽略这种效应,则会出现异常点数值与较晚时刻数据数值相当的情况,存在误判的可能,需研究符合时间序列相关性特征的算法。4) Density (distance) detection algorithm: The traditional density (distance)-based algorithm does not consider the influence of the weights of the sequence before and after on evaluating the tightness of the current point. For time-varying sources, although there is a certain correlation between the data before and after. , but this correlation becomes weaker as the interval grows. If this effect is ignored, the value of the abnormal point will be equivalent to the value of the data at a later time, and there is a possibility of misjudgment. It is necessary to study an algorithm that conforms to the correlation characteristics of time series.

本申请提出以下实施例以解决上述确定,具体实施例说明如下:The application proposes the following embodiments to solve the above determination, and the specific embodiments are described as follows:

实施例一:Example 1:

根据本发明实施例,提供了一种飞行器遥测数据的异常检测方法的实施例,需要说明的是,在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行,并且,虽然在流程图中示出了逻辑顺序,但是在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤。According to an embodiment of the present invention, an embodiment of an abnormality detection method for aircraft telemetry data is provided. It should be noted that the steps shown in the flowchart of the accompanying drawings may be executed in a computer system such as a set of computer-executable instructions. and, although a logical order is shown in the flowcharts, in some cases the steps shown or described may be performed in an order different from that herein.

图1是根据本发明实施例的一种飞行器遥测数据的异常检测方法的流程图,如图1所示,该方法包括如下步骤:FIG. 1 is a flowchart of a method for detecting anomalies in telemetry data of an aircraft according to an embodiment of the present invention. As shown in FIG. 1 , the method includes the following steps:

步骤S102,获取飞行器发送的先导遥测数据,并通过平稳性检测确定所述先导遥测数据是否为平稳时间序列;Step S102, obtaining the pilot telemetry data sent by the aircraft, and determining whether the pilot telemetry data is a stationary time series through stationarity detection;

需要说明的是,上述的先导遥测数据为在当前时刻节点之前获取到的遥测数据,且先导遥测数据的数据量较小,能够便捷的判断出先导遥测数据是否为平稳时间序列。It should be noted that the above-mentioned pilot telemetry data is the telemetry data obtained before the node at the current moment, and the data amount of the pilot telemetry data is small, which can conveniently determine whether the pilot telemetry data is a stationary time series.

步骤S104,在完成平稳性检测之后,获取飞行器发送的当前时刻节点的遥测数据;Step S104, after completing the stationarity detection, obtain the telemetry data of the node at the current moment sent by the aircraft;

步骤S106,若所述先导遥测数据为平稳时间序列,则利用边界检测算法对当前时刻节点的遥测数据进行检测,确定所述当前时刻节点的遥测数据是否为异常数据;Step S106, if the pilot telemetry data is a stationary time series, use a boundary detection algorithm to detect the telemetry data of the node at the current moment to determine whether the telemetry data of the node at the current moment is abnormal data;

具体的,利用边界检测算法对当前时刻节点的遥测数据进行检测的步骤如下:Specifically, the steps of using the boundary detection algorithm to detect the telemetry data of the node at the current moment are as follows:

确定当前时刻节点的遥测数据的遥测参数是否处于预设边界范围内,如果当前时刻节点的遥测数据的遥测参数不处于预设边界范围内,则当前时刻节点的遥测数据为异常数据。It is determined whether the telemetry parameter of the telemetry data of the node at the current moment is within the preset boundary range, and if the telemetry parameter of the telemetry data of the node at the current moment is not within the preset boundary range, the telemetry data of the node at the current moment is abnormal data.

步骤S108,若所述先导遥测数据不是平稳时间序列,则利用预设检测算法对所述当前时刻节点的遥测数据进行检测,确定所述当前时刻节点的遥测数据是否为异常数据,其中,所述预设检测算法包括以下至少之一:三点集同构映射检测算法,左右双陪集赋权映射检测算法。Step S108, if the pilot telemetry data is not a stationary time series, use a preset detection algorithm to detect the telemetry data of the node at the current moment to determine whether the telemetry data of the node at the current moment is abnormal data, wherein the The preset detection algorithm includes at least one of the following: a three-point set isomorphic map detection algorithm, and a left and right dual coset weighted map detection algorithm.

在本发明实施例中,通过对先导遥测数据进行平稳性检测,确定出先导遥测数据是否为平稳时间序列,并根据检测结果选择不同的检测算法对当前时刻节点的遥测数据进行检测,从而确定出当前时刻节点的遥测数据是否为异常数据,达到了对遥测数据进行异常检测的目的,进而解决了现有技术中对遥测数据的进行异常数据检测的检测步骤复杂,且检测准确率较低的技术问题,从而实现了简化了遥测数据的异常数据检测的检测步骤,提高了检测准确率的技术效果。In the embodiment of the present invention, by performing stationarity detection on the pilot telemetry data, it is determined whether the pilot telemetry data is a stationary time series, and different detection algorithms are selected according to the detection results to detect the telemetry data of the node at the current moment, so as to determine whether the pilot telemetry data is a stationary time series. Whether the telemetry data of the node at the current moment is abnormal data, the purpose of abnormality detection of telemetry data is achieved, and the technology of abnormal data detection of telemetry data in the prior art is solved with complicated detection steps and low detection accuracy. Therefore, the detection steps of abnormal data detection of telemetry data are simplified, and the technical effect of improving detection accuracy is realized.

在本发明实施例中,如图2所示,所述方法还包括如下步骤:In this embodiment of the present invention, as shown in FIG. 2 , the method further includes the following steps:

步骤S110,对所述异常数据进行清洗。Step S110, cleaning the abnormal data.

在本发明实施例中,当检测出当前时刻节点的遥测数据为异常数据之后,首先,确定出所述第一目标数据为所述当前时刻节点之前的第一个正常遥测数据(即,第一目标数据),所述第二目标数据为所述当前时刻节点之后的第一个正常遥测数据(即,第二目标数据)。In this embodiment of the present invention, after detecting that the telemetry data of the node at the current moment is abnormal data, first, it is determined that the first target data is the first normal telemetry data before the node at the current moment (that is, the first target data), the second target data is the first normal telemetry data (ie, the second target data) after the node at the current moment.

然后,求解出第一目标数据与第二目标数据的平均值,然后用该平均值替换异常数据,从而达到对异常数据进行清洗的技术效果。Then, the average value of the first target data and the second target data is obtained, and then the abnormal data is replaced by the average value, so as to achieve the technical effect of cleaning the abnormal data.

通过清洗异常数据,提高了遥测数据的准确性,将完成清洗的遥测数据作为后续特征提取与弹性压缩的输入数据,能够为深空探测任务及空间信息网络的发展提供技术支撑。By cleaning abnormal data, the accuracy of telemetry data is improved, and the cleaned telemetry data is used as input data for subsequent feature extraction and elastic compression, which can provide technical support for the development of deep space exploration missions and spatial information networks.

在本发明实施例中,步骤S108还包括如下步骤:In this embodiment of the present invention, step S108 further includes the following steps:

步骤S11,确定所述飞行器的性能类型,其中,所述性能类型包括:第一性能类型和第二性能类型,所述第一性能类型的飞行器的在轨处理能力低于所述第二性能类型的飞行器的在轨处理能力;Step S11, determining the performance type of the aircraft, wherein the performance type includes: a first performance type and a second performance type, and the on-orbit processing capability of the aircraft of the first performance type is lower than that of the second performance type the on-orbit processing capability of the aircraft;

步骤S12,基于所述性能类型,确定出所述预设检测算法中的目标检测算法;Step S12, determining a target detection algorithm in the preset detection algorithm based on the performance type;

步骤S13,利用所述目标检测算法,对所述当前时刻节点的遥测数据进行检测,确定所述当前时刻节点的遥测数据是否为异常数据。In step S13, the target detection algorithm is used to detect the telemetry data of the node at the current moment to determine whether the telemetry data of the node at the current moment is abnormal data.

在本发明实施例中,为了简化异常数据检测的检测步骤,以及提高检测准确性,在对异常数据进行检测时,首先,需要确定出发送当前时刻节点的遥测数据的飞行器的性能类型,并根据性能类型选择对应的异常数据检测算法,对当前时刻节点的遥测数据进行检测。In the embodiment of the present invention, in order to simplify the detection steps of abnormal data detection and improve the detection accuracy, when detecting abnormal data, first, it is necessary to determine the performance type of the aircraft that sends the telemetry data of the node at the current moment, and according to Select the corresponding abnormal data detection algorithm for the performance type to detect the telemetry data of the node at the current moment.

如果飞行器的在轨处理能力较弱,则采用三点集同构映射检测算法对当前时刻节点的遥测数据进行检测。If the on-orbit processing capability of the aircraft is weak, the three-point set isomorphic mapping detection algorithm is used to detect the telemetry data of the node at the current moment.

如果飞行器的在轨处理能力较强,则采用左右双陪集赋权映射检测算法对当前时刻节点的遥测数据进行检测。If the on-orbit processing capability of the aircraft is strong, the left and right double coset weighted mapping detection algorithm is used to detect the telemetry data of the node at the current moment.

在本发明实施例中,如图3所示,利用三点集同构映射检测算法检测异常数据包括如下步骤:In the embodiment of the present invention, as shown in FIG. 3 , detecting abnormal data by using a three-point set isomorphic mapping detection algorithm includes the following steps:

步骤S21,确定出所述当前时刻节点的邻居节点,其中,所述邻居节点包括:第一邻居节点和第二邻居节点,所述第一邻居节点为所述当前时刻节点之前的邻居节点,所述第二邻居节点为所述当前时刻节点之后的邻居节点;Step S21: Determine the neighbor node of the node at the current moment, wherein the neighbor node includes: a first neighbor node and a second neighbor node, the first neighbor node is the neighbor node before the node at the current moment, so The second neighbor node is the neighbor node after the current moment node;

步骤S22,利用所述当前时刻节点和所述第一邻居节点,构建第一矢量线段,以及利用所述当前时刻节点所述第二邻居节点构建的第二矢量线段,并确定出所述第一矢量线段与所述第二矢量线段构成的矢量夹角的角度;Step S22, use the current moment node and the first neighbor node to construct a first vector line segment, and use the current moment node to construct a second vector line segment constructed by the second neighbor node, and determine the first vector line segment. The angle between the vector line segment and the vector angle formed by the second vector line segment;

步骤S23,若所述角度小于预设阈值,则所述当前时刻节点的遥测数据为异常数据。Step S23, if the angle is smaller than a preset threshold, the telemetry data of the node at the current moment is abnormal data.

在本发明实施例中,如图4所示,利用左右双陪集赋权映射检测算法检测异常数据包括如下步骤:In the embodiment of the present invention, as shown in FIG. 4 , detecting abnormal data by using the left and right bi-coset weighted mapping detection algorithm includes the following steps:

步骤S31,构建所述当前时刻节点的k元双陪集;Step S31, constructing the k-ary double coset of the node at the current moment;

步骤S32,计算所述k元双陪集的距离阈值,以及为所述k元双陪集中的节点赋值,得到节点赋值;Step S32, calculating the distance threshold of the k-yuan double coset, and assigning a value to the node in the k-yuan double coset to obtain the node assignment;

步骤S33,确定出所述当前时刻节点的紧密因子比判阈值;Step S33, determining the tightness factor comparison threshold of the node at the current moment;

步骤S34,结合所述k元双陪集的距离阈值和所述节点赋值,计算出所述当前时刻节点的紧密度;Step S34, calculating the tightness of the node at the current moment in combination with the distance threshold of the k-yuan double coset and the node assignment;

步骤S35,若所述紧密度小于所述紧密因子比判阈值,则所述当前时刻节点的遥测数据为异常数据。Step S35, if the closeness is less than the closeness factor comparison threshold, the telemetry data of the node at the current moment is abnormal data.

下面结合图3和图4对上述两种检测算法的检测过程进行详细说明:The detection process of the above two detection algorithms will be described in detail below in conjunction with Figure 3 and Figure 4:

一,利用三点集同构映射检测算法检测异常数据:First, use the three-point set isomorphic mapping detection algorithm to detect abnormal data:

由根据局部三点数据平滑线性假定和3σ原则可知,正常数据不应偏离理论值的3倍子集内部差分均值(即预设阈值为72°)。According to the assumption of smooth linearity of local three-point data and the 3σ principle, the normal data should not deviate from the average value of the internal difference of the subsets that is three times the theoretical value (that is, the preset threshold is 72°).

首先,需要确定出当前时刻节点之前的邻居节点(即第一邻居节点),以及确定出当前时刻节点之后的邻居节点(即第二邻居节点)。First, it is necessary to determine the neighbor node (ie, the first neighbor node) before the node at the current moment, and determine the neighbor node (ie, the second neighbor node) after the node at the current moment.

然后,连接当前时刻节点和所述第一邻居节点,构建第一矢量线段,以及链接当前时刻节点所述第二邻居节点构建的第二矢量线段,从而构成矢量夹角。Then, connect the node at the current moment and the first neighbor node to construct a first vector line segment, and link the node at the current moment and the second vector line segment constructed by the second neighbor node to form a vector angle.

接着,确定出该矢量夹角的角度,如果该角度小于72°,则当前时刻节点的遥测数据为异常数据。Next, the angle of the included angle of the vector is determined. If the angle is less than 72°, the telemetry data of the node at the current moment is abnormal data.

如果该角度大于72°,则当前时刻节点的遥测数据为正常数据。If the angle is greater than 72°, the telemetry data of the node at the current moment is normal data.

二,利用左右双陪集赋权映射检测算法检测异常数据:Second, use the left and right double coset weighted mapping detection algorithm to detect abnormal data:

1,构建当前时刻节点的左/右k元双陪集,其中,k为正整数。1. Construct the left/right k-ary double coset of the node at the current moment, where k is a positive integer.

2,计算出左陪集的距离阈值Dl和右陪集的距离阈值Dr,具体计算算是如下:

Figure BDA0002357943540000111
2. Calculate the distance threshold D l of the left coset and the distance threshold D r of the right coset. The specific calculation is as follows:
Figure BDA0002357943540000111

3,将左/右陪集内各个时刻节点按远近效应以幂指数(αlr)特性赋权映射到当前时刻节点,每个时刻节点赋予权值,pj=(αl)j,pi=(αr)i,即:时间上接近的时刻节点具有较高权值,时间上远离的时刻节点具有低权值(即,远近效应)。3. Map each moment node in the left/right coset to the current moment node according to the far-near effect with power exponent (α lr ) characteristics, each moment node is assigned a weight, p j =(α l ) j , p i =(α r ) i , that is, time nodes that are close in time have higher weights, and time nodes that are far away in time have low weights (ie, far-near effects).

4,确定出当前时刻节点的紧密因子比判阈值T,T的取值取决于前时刻节点的邻域范围内可能存在的异常数据的数目,异常数据的数量为1时,T为1/2,异常数据的数量为2时,T为1/2+1/4,依次类推。一般情况下默认为异常数据的数量为1,即T=1/2。4. Determine the close factor ratio threshold T of the node at the current moment. The value of T depends on the number of abnormal data that may exist in the neighborhood of the node at the previous moment. When the number of abnormal data is 1, T is 1/2 , when the number of abnormal data is 2, T is 1/2+1/4, and so on. In general, the default is that the number of abnormal data is 1, that is, T=1/2.

5,在距离阈值(Dl,Dr)范围内对左右陪集权值进行求和运算得到紧密因子C(k,Dl,Dr),即统计前后2k个点数据集中距离当前点xn的1-范数不大于距离阈值的权值之和作为紧密度衡量指标“紧密因子”,如下式所示:

Figure BDA0002357943540000112
5. In the range of the distance threshold (D l , D r ), the left and right coset weights are summed to obtain the compact factor C (k, D l , D r ), that is, the distance from the current point x in the 2k point data set before and after statistics The sum of the weights whose 1-norm of n is not greater than the distance threshold is used as the closeness measurement index "closeness factor", as shown in the following formula:
Figure BDA0002357943540000112

根据上述公式,计算出当前时刻节点的紧密度。According to the above formula, the tightness of the node at the current moment is calculated.

6,如果上述的紧密度小于紧密因子比判阈值T,则当前时刻节点的遥测数据为异常数据,如果上述的紧密度大于或等于紧密因子比判阈值T,则当前时刻节点的遥测数据为正常数据。6. If the above-mentioned tightness is less than the tightness factor ratio judgment threshold T, then the telemetry data of the node at the current moment is abnormal data; if the above-mentioned tightness is greater than or equal to the tightness factor ratio judgment threshold T, then the telemetry data of the node at the current moment is normal. data.

实施例二:Embodiment 2:

本发明还提供了一种飞行器遥测数据的异常检测装置的实施例,该装置用于执行本发明实施例上述内容所提供的遥测数据的异常检测方法,以下是本发明实施例提供的遥测数据的异常检测装置的具体介绍。The present invention also provides an embodiment of an abnormality detection device for aircraft telemetry data, which is used to execute the method for detecting abnormality in telemetry data provided by the above content of the embodiments of the present invention. The following is an example of the telemetry data provided by the embodiments of the present invention. The specific introduction of abnormal detection device.

如图5所示,上述的遥测数据的异常检测装置包括:第一获取单元10,第二获取单元20,第一检测单元30和第二检测单元40。As shown in FIG. 5 , the above-mentioned abnormality detection apparatus for telemetry data includes: a first acquisition unit 10 , a second acquisition unit 20 , a first detection unit 30 and a second detection unit 40 .

所述第一获取单元10用于获取飞行器发送的先导遥测数据,并通过平稳性检测确定所述先导遥测数据是否为平稳时间序列;The first acquisition unit 10 is configured to acquire pilot telemetry data sent by the aircraft, and determine whether the pilot telemetry data is a stationary time series through stationarity detection;

所述第二获取单元20用于在完成平稳性检测之后,获取飞行器发送的当前时刻节点的遥测数据;The second obtaining unit 20 is configured to obtain the telemetry data of the node at the current moment sent by the aircraft after the stationarity detection is completed;

所述第一检测单元30用于在所述先导遥测数据为平稳时间序列的情况下,利用边界检测算法对当前时刻节点的遥测数据进行检测,确定所述当前时刻节点的遥测数据是否为异常数据;The first detection unit 30 is configured to use a boundary detection algorithm to detect the telemetry data of the node at the current moment when the pilot telemetry data is a stationary time series, and determine whether the telemetry data of the node at the current moment is abnormal data. ;

所述第二检测单元40用于在所述先导遥测数据不是平稳时间序列的情况下,则利用预设检测算法对所述当前时刻节点的遥测数据进行检测,确定所述当前时刻节点的遥测数据是否为异常数据,其中,所述预设检测算法包括以下至少之一:三点集同构映射检测算法,左右双陪集赋权映射检测算法。The second detection unit 40 is configured to use a preset detection algorithm to detect the telemetry data of the node at the current moment when the pilot telemetry data is not a stationary time series, and determine the telemetry data of the node at the current moment. Whether it is abnormal data, wherein the preset detection algorithm includes at least one of the following: a three-point set isomorphic mapping detection algorithm, and a left and right dual coset weighted mapping detection algorithm.

在本发明实施例中,通过对先导遥测数据进行平稳性检测,确定出先导遥测数据是否为平稳时间序列,并根据检测结果选择不同的检测算法对当前时刻节点的遥测数据进行检测,从而确定出当前时刻节点的遥测数据是否为异常数据,达到了对遥测数据进行异常检测的目的,进而解决了现有技术中对遥测数据的进行异常数据检测的检测步骤复杂,且检测准确率较低的技术问题,从而实现了简化了遥测数据的异常数据检测的检测步骤,提高了检测准确率的技术效果。In the embodiment of the present invention, by performing stationarity detection on the pilot telemetry data, it is determined whether the pilot telemetry data is a stationary time series, and different detection algorithms are selected according to the detection results to detect the telemetry data of the node at the current moment, so as to determine whether the pilot telemetry data is a stationary time series. Whether the telemetry data of the node at the current moment is abnormal data, the purpose of abnormality detection of telemetry data is achieved, and the technology of abnormal data detection of telemetry data in the prior art is solved with complicated detection steps and low detection accuracy. Therefore, the detection steps of abnormal data detection of telemetry data are simplified, and the technical effect of improving detection accuracy is realized.

优选地,所述装置还包括:数据清洗单元,用于对所述异常数据进行清洗。Preferably, the device further includes: a data cleaning unit, configured to clean the abnormal data.

优选地,所述数据清洗单元用于:确定出第一目标数据和第二目标数据,其中,所述第一目标数据为所述当前时刻节点之前的第一个正常遥测数据,所述第二目标数据为所述当前时刻节点之后的第一个正常遥测数据;计算所述第一目标数据和所述第二目标数据的均值,并将所述异常数据替换为所述均值。Preferably, the data cleaning unit is configured to: determine first target data and second target data, wherein the first target data is the first normal telemetry data before the node at the current moment, and the second target data is the first normal telemetry data before the node at the current moment. The target data is the first normal telemetry data after the node at the current moment; the average value of the first target data and the second target data is calculated, and the abnormal data is replaced with the average value.

优选地,所述第二检测单元用于:确定所述飞行器的性能类型,其中,所述性能类型包括:第一性能类型和第二性能类型,所述第一性能类型的飞行器的在轨处理能力低于所述第二性能类型的飞行器的在轨处理能力;基于所述性能类型,确定出所述预设检测算法中的目标检测算法;利用所述目标检测算法,对所述当前时刻节点的遥测数据进行检测,确定所述当前时刻节点的遥测数据是否为异常数据。Preferably, the second detection unit is used to: determine a performance type of the aircraft, wherein the performance type includes: a first performance type and a second performance type, and the on-orbit processing of the aircraft of the first performance type The on-orbit processing capability of the aircraft whose capability is lower than the second performance type; the target detection algorithm in the preset detection algorithm is determined based on the performance type; the target detection algorithm is used to detect the node at the current moment The telemetry data of the node is detected to determine whether the telemetry data of the node at the current moment is abnormal data.

优选地,所述第二检测单元用于:若所述飞行器的性能类型为第一性能类型,则所述目标检测算法为所述三点集同构映射检测算法;若所述飞行器的性能类型为第二性能类型,则所述目标检测算法为所述左右双陪集赋权映射检测算法。Preferably, the second detection unit is configured to: if the performance type of the aircraft is the first performance type, the target detection algorithm is the three-point set isomorphic mapping detection algorithm; if the performance type of the aircraft is the second performance type, the target detection algorithm is the left and right dual coset weighted mapping detection algorithm.

优选地,若所述目标检测算法为所述三点集同构映射检测算法,所述第二检测单元用于:确定出所述当前时刻节点的邻居节点,其中,所述邻居节点包括:第一邻居节点和第二邻居节点,所述第一邻居节点为所述当前时刻节点之前的邻居节点,所述第二邻居节点为所述当前时刻节点之后的邻居节点;利用所述当前时刻节点和所述第一邻居节点,构建第一矢量线段,以及利用所述当前时刻节点所述第二邻居节点构建的第二矢量线段,并确定出所述第一矢量线段与所述第二矢量线段构成的矢量夹角的角度;若所述角度小于预设阈值,则所述当前时刻节点的遥测数据为异常数据。Preferably, if the target detection algorithm is the three-point set isomorphic mapping detection algorithm, the second detection unit is configured to: determine the neighbor nodes of the node at the current moment, wherein the neighbor nodes include: a neighbor node and a second neighbor node, the first neighbor node is the neighbor node before the current moment node, and the second neighbor node is the neighbor node after the current moment node; using the current moment node and The first neighbor node constructs a first vector line segment, and uses the second vector line segment constructed by the second neighbor node of the current moment node, and determines that the first vector line segment and the second vector line segment constitute If the angle is smaller than the preset threshold, the telemetry data of the node at the current moment is abnormal data.

优选地,若所述目标检测算法为所述左右双陪集赋权映射检测算法,所述第二检测单元用于:构建所述当前时刻节点的k元双陪集;计算所述k元双陪集的距离阈值,以及为所述k元双陪集中的节点赋值,得到节点赋值;确定出所述当前时刻节点的紧密因子比判阈值;结合所述k元双陪集的距离阈值和所述节点赋值,计算出所述当前时刻节点的紧密度;若所述紧密度小于紧密因子比判阈值,则所述当前时刻节点的遥测数据为异常数据。Preferably, if the target detection algorithm is the left and right dual coset weighted mapping detection algorithm, the second detection unit is used for: constructing a k-ary dual coset of the node at the current moment; calculating the k-ary dual coset The distance threshold of the coset, and assigning values to the nodes in the k-ary double coset to obtain the node assignment; determining the closeness factor ratio judgment threshold of the node at the current moment; combining the distance threshold of the k-ary double coset and all The node assignment is performed to calculate the tightness of the node at the current moment; if the tightness is less than the tightness factor comparison threshold, the telemetry data of the node at the current moment is abnormal data.

优选地,第一检测单元用于:确定所述当前时刻节点的遥测数据的遥测参数是否处于预设边界范围内;若否,则所述当前时刻节点的遥测数据为异常数据。Preferably, the first detection unit is configured to: determine whether the telemetry parameter of the telemetry data of the node at the current moment is within a preset boundary range; if not, the telemetry data of the node at the current moment is abnormal data.

本申请实施例还提供了一种具有处理器可执行的非易失的程序代码的计算机可读介质,所述程序代码使处理器执行上述实施例一所述的飞行器遥测数据的异常检测方法。The embodiment of the present application also provides a computer-readable medium having a non-volatile program code executable by a processor, the program code enables the processor to execute the abnormality detection method for aircraft telemetry data described in the first embodiment.

另外,在本发明实施例的描述中,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。In addition, in the description of the embodiments of the present invention, unless otherwise expressly specified and limited, the terms "installed", "connected" and "connected" should be understood in a broad sense, for example, it may be a fixed connection or a detachable connection , or integrally connected; it can be a mechanical connection or an electrical connection; it can be a direct connection, or an indirect connection through an intermediate medium, or the internal communication between the two components. For those of ordinary skill in the art, the specific meanings of the above terms in the present invention can be understood in specific situations.

在本发明的描述中,需要说明的是,术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一”、“第二”、“第三”仅用于描述目的,而不能理解为指示或暗示相对重要性。In the description of the present invention, it should be noted that the terms "center", "upper", "lower", "left", "right", "vertical", "horizontal", "inner", "outer", etc. The indicated orientation or positional relationship is based on the orientation or positional relationship shown in the accompanying drawings, which is only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying that the indicated device or element must have a specific orientation or a specific orientation. construction and operation, and therefore should not be construed as limiting the invention. Furthermore, the terms "first", "second", and "third" are used for descriptive purposes only and should not be construed to indicate or imply relative importance.

在本申请所提供的几个实施例中,应该理解到,所揭露的系统、装置和方法,可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,又例如,多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些通信接口,装置或单元的间接耦合或通信连接,可以是电性,机械或其它的形式。In the several embodiments provided in this application, it should be understood that the disclosed system, apparatus and method may be implemented in other manners. The apparatus embodiments described above are only illustrative. For example, the division of the units is only a logical function division. In actual implementation, there may be other division methods. For example, multiple units or components may be combined or Can be integrated into another system, or some features can be ignored, or not implemented. On the other hand, the shown or discussed mutual coupling or direct coupling or communication connection may be through some communication interfaces, indirect coupling or communication connection of devices or units, which may be in electrical, mechanical or other forms.

所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。The units described as separate components may or may not be physically separated, and components displayed as units may or may not be physical units, that is, may be located in one place, or may be distributed to multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution in this embodiment.

另外,在本发明各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。In addition, each functional unit in each embodiment of the present invention may be integrated into one processing unit, or each unit may exist physically alone, or two or more units may be integrated into one unit.

最后应说明的是:以上所述实施例,仅为本发明的具体实施方式,用以说明本发明的技术方案,而非对其限制,本发明的保护范围并不局限于此,尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化,或者对其中部分技术特征进行等同替换;而这些修改、变化或者替换,并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应所述以权利要求的保护范围为准。Finally, it should be noted that the above-mentioned embodiments are only specific implementations of the present invention, and are used to illustrate the technical solutions of the present invention, but not to limit them. The protection scope of the present invention is not limited thereto, although referring to the foregoing The embodiment has been described in detail the present invention, those of ordinary skill in the art should understand: any person skilled in the art who is familiar with the technical field within the technical scope disclosed by the present invention can still modify the technical solutions described in the foregoing embodiments. Or can easily think of changes, or equivalently replace some of the technical features; and these modifications, changes or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present invention, and should be covered in the present invention. within the scope of protection. Therefore, the protection scope of the present invention should be based on the protection scope of the claims.

Claims (10)

1. A method of anomaly detection of aircraft telemetry data, comprising:
acquiring pilot telemetering data sent by an aircraft, and determining whether the pilot telemetering data is a stationary time sequence or not through stationary detection;
after stationarity detection is finished, acquiring telemetering data of a node at the current time, which is sent by an aircraft;
if the pilot telemetering data is a stable time sequence, detecting the telemetering data of the node at the current moment by using a boundary detection algorithm, and determining whether the telemetering data of the node at the current moment is abnormal data;
if the pilot telemetering data is not a stationary time sequence, detecting the telemetering data of the current time node by using a preset detection algorithm, and determining whether the telemetering data of the current time node is abnormal data, wherein the preset detection algorithm comprises at least one of the following: a three-point set isomorphic mapping detection algorithm and a left and right double-coset empowerment mapping detection algorithm.
2. The method of claim 1, further comprising:
and cleaning the abnormal data.
3. The method of claim 2, wherein cleansing the anomaly data comprises:
determining first target data and second target data, wherein the first target data is first normal telemetering data before the node at the current moment, and the second target data is first normal telemetering data after the node at the current moment;
calculating a mean value of the first target data and the second target data, and replacing the abnormal data with the mean value.
4. The method of claim 1, wherein detecting the telemetry data of the current time node by using a preset detection algorithm to determine whether the telemetry data of the current time node is abnormal data comprises:
determining a performance type of the aircraft, wherein the performance type comprises: a first performance type and a second performance type, the in-orbit handling capability of the aircraft of the first performance type being lower than the in-orbit handling capability of the aircraft of the second performance type;
determining a target detection algorithm in the preset detection algorithms based on the performance type;
and detecting the telemetering data of the current time node by using the target detection algorithm, and determining whether the telemetering data of the current time node is abnormal data.
5. The method of claim 4, wherein determining a target detection algorithm of the predetermined detection algorithms based on the performance type comprises:
if the performance type of the aircraft is a first performance type, the target detection algorithm is the three-point set isomorphic mapping detection algorithm;
and if the performance type of the aircraft is a second performance type, the target detection algorithm is the left and right double-coset empowerment mapping detection algorithm.
6. The method of claim 5, wherein if the target detection algorithm is the three-point set isomorphic mapping detection algorithm;
detecting the telemetering data of the current time node by using the target detection algorithm, and determining whether the telemetering data of the current time node is abnormal data, wherein the method comprises the following steps:
determining a neighbor node of the current time node, wherein the neighbor node comprises: the first neighbor node is a neighbor node before the node at the current moment, and the second neighbor node is a neighbor node after the node at the current moment;
constructing a first vector line segment by using the current time node and the first neighbor node, constructing a second vector line segment by using the current time node and the second neighbor node, and determining the angle of a vector included angle formed by the first vector line segment and the second vector line segment;
and if the angle is smaller than a preset threshold value, the telemetering data of the node at the current moment is abnormal data.
7. The method according to claim 5, wherein if the target detection algorithm is the left-right double coset empowerment mapping detection algorithm;
detecting the telemetering data by using the target detection algorithm to determine abnormal data in the telemetering data, wherein the method comprises the following steps:
constructing a k-element double coset of the node at the current moment;
calculating a distance threshold of the k-element double coset, and assigning a node in the k-element double coset to obtain a node assignment;
determining a tightness factor comparison threshold value of the node at the current moment;
calculating the compactness of the node at the current moment by combining the distance threshold of the k-element double coset and the node assignment;
and if the compactness is smaller than the compactness factor comparison threshold, the telemetering data of the node at the current moment is abnormal data.
8. The method of claim 1, wherein detecting telemetry data of a node at a current time using a boundary detection algorithm to determine whether the telemetry data of the node at the current time is abnormal data comprises:
determining whether the telemetry parameters of the telemetry data of the node at the current moment are within a preset boundary range;
and if not, the telemetering data of the node at the current moment is abnormal data.
9. An anomaly detection device for aircraft telemetry data, comprising: a first acquisition unit, a second acquisition unit, a first detection unit and a second detection unit, wherein,
the first acquisition unit is used for acquiring pilot telemetering data sent by an aircraft and determining whether the pilot telemetering data is a stationary time sequence or not through stationarity detection;
the second acquisition unit is used for acquiring the telemetering data of the node at the current moment sent by the aircraft after the stationarity detection is finished;
the first detection unit is used for detecting the telemetering data of the node at the current moment by using a boundary detection algorithm under the condition that the pilot telemetering data is a stable time sequence, and determining whether the telemetering data of the node at the current moment is abnormal data;
the second detection unit is configured to, if the pilot telemetry data is not a stationary time series, detect the telemetry data of the current time node by using a preset detection algorithm, and determine whether the telemetry data of the current time node is abnormal data, where the preset detection algorithm includes at least one of: a three-point set isomorphic mapping detection algorithm and a left and right double-coset empowerment mapping detection algorithm.
10. A computer-readable medium having non-volatile program code executable by a processor, the program code causing the processor to perform the method of anomaly detection of aircraft telemetry data as claimed in any one of claims 1 to 8.
CN202010016293.8A 2020-01-07 2020-01-07 Anomaly detection method and device for aircraft telemetry data Active CN111241158B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010016293.8A CN111241158B (en) 2020-01-07 2020-01-07 Anomaly detection method and device for aircraft telemetry data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010016293.8A CN111241158B (en) 2020-01-07 2020-01-07 Anomaly detection method and device for aircraft telemetry data

Publications (2)

Publication Number Publication Date
CN111241158A true CN111241158A (en) 2020-06-05
CN111241158B CN111241158B (en) 2020-09-18

Family

ID=70877672

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010016293.8A Active CN111241158B (en) 2020-01-07 2020-01-07 Anomaly detection method and device for aircraft telemetry data

Country Status (1)

Country Link
CN (1) CN111241158B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112506908A (en) * 2020-12-10 2021-03-16 云南电网有限责任公司玉溪供电局 Electric energy metering data cleaning method and system
CN112526559A (en) * 2020-12-03 2021-03-19 北京航空航天大学 System relevance state monitoring method under multi-working-condition
CN114896237A (en) * 2022-06-02 2022-08-12 西安航天动力研究所 Liquid rocket engine flight telemetry data cleaning method and device

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102928655A (en) * 2012-11-26 2013-02-13 慈溪市供电局 Method for detecting power abnormal data
CN103646167A (en) * 2013-11-22 2014-03-19 北京空间飞行器总体设计部 Satellite abnormal condition detection system based on telemeasuring data
CN103793765A (en) * 2014-02-12 2014-05-14 航天东方红卫星有限公司 Satellite telemetering data predicting method based on Kalman smoothing
CN104899327A (en) * 2015-06-24 2015-09-09 哈尔滨工业大学 Method for detecting abnormal time sequence without class label
CN104915846A (en) * 2015-06-18 2015-09-16 北京京东尚科信息技术有限公司 Electronic commerce time sequence data anomaly detection method and system
CN106055919A (en) * 2016-08-09 2016-10-26 航天东方红卫星有限公司 Satellite abnormity detection method based on telemetry data wavelet transformation
US10073753B2 (en) * 2016-02-14 2018-09-11 Dell Products, Lp System and method to assess information handling system health and resource utilization
CN110460622A (en) * 2019-09-12 2019-11-15 贵州电网有限责任公司 A kind of network anomaly detection method based on Situation Awareness prediction technique

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102928655A (en) * 2012-11-26 2013-02-13 慈溪市供电局 Method for detecting power abnormal data
CN103646167A (en) * 2013-11-22 2014-03-19 北京空间飞行器总体设计部 Satellite abnormal condition detection system based on telemeasuring data
CN103793765A (en) * 2014-02-12 2014-05-14 航天东方红卫星有限公司 Satellite telemetering data predicting method based on Kalman smoothing
CN104915846A (en) * 2015-06-18 2015-09-16 北京京东尚科信息技术有限公司 Electronic commerce time sequence data anomaly detection method and system
CN104899327A (en) * 2015-06-24 2015-09-09 哈尔滨工业大学 Method for detecting abnormal time sequence without class label
US10073753B2 (en) * 2016-02-14 2018-09-11 Dell Products, Lp System and method to assess information handling system health and resource utilization
US20190057006A1 (en) * 2016-02-14 2019-02-21 Dell Products, Lp System and method to assess information handling system health and resource utilization
CN106055919A (en) * 2016-08-09 2016-10-26 航天东方红卫星有限公司 Satellite abnormity detection method based on telemetry data wavelet transformation
CN110460622A (en) * 2019-09-12 2019-11-15 贵州电网有限责任公司 A kind of network anomaly detection method based on Situation Awareness prediction technique

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
DU YING 等: "Anomaly detection of orbit satellite telemetry sequence based on two-window mode", 《THE 30TH CHINESE CONTROL AND DECISION CONFERENCE》 *
YI PENG 等: "Telemetry fault detection for meteorological satellite based on PCA", 《INTERNATIONAL SYMPOSIUM ON COMMUNICATIONS AND INFORMATION TECHNOLOGIES》 *
彭喜元 等: "航天器遥测数据异常检测综述", 《仪器仪表学报》 *
杨甲森: "基于遥测数据相关性的航天器异常检测", 《仪器仪表学报》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112526559A (en) * 2020-12-03 2021-03-19 北京航空航天大学 System relevance state monitoring method under multi-working-condition
CN112506908A (en) * 2020-12-10 2021-03-16 云南电网有限责任公司玉溪供电局 Electric energy metering data cleaning method and system
CN114896237A (en) * 2022-06-02 2022-08-12 西安航天动力研究所 Liquid rocket engine flight telemetry data cleaning method and device

Also Published As

Publication number Publication date
CN111241158B (en) 2020-09-18

Similar Documents

Publication Publication Date Title
CN111241158B (en) Anomaly detection method and device for aircraft telemetry data
JP7010641B2 (en) Abnormality diagnosis method and abnormality diagnosis device
CN102186239B (en) Radio frequency (RF) fingerprint positing method, device and network equipment
JP5437328B2 (en) Observation value reliability evaluation apparatus, observation value reliability evaluation method, and observation value reliability evaluation program
CN110392046B (en) Method and device for detecting abnormity of network access
WO2018142703A1 (en) Anomaly factor estimation device, anomaly factor estimation method, and program
WO2011112562A2 (en) Efficient channel search with sequential probability ratio testing
CN105451260A (en) Network request method, and method and device for measuring network fluctuation
CN109921925A (en) A kind of dial testing method and device
CN109862598A (en) Residence reselecting, device and computer readable storage medium
CN106603324A (en) Training set acquisition device and training set acquisition method
CN110380925B (en) Port selection method and system in network equipment detection
CN111682988A (en) Remote control method, device, storage medium and processor
CN112214677A (en) A point of interest recommendation method, device, electronic device and storage medium
CN110502432B (en) Intelligent test method, device, equipment and readable storage medium
Karimi et al. Recursive dynamic compressive sensing in smart distribution systems
CN109560978B (en) Network flow detection method, device and system and computer readable storage medium
Tsang et al. Nonparametric internet tomography
CN109491702A (en) A kind of prioritization scheme determines method, apparatus, terminal device and storage medium
EP3128728B1 (en) Method and device for detecting a broken binder
US20230106935A1 (en) Network probe placement optimization
US20240338596A1 (en) Signalling model performance feedback information (mpfi)
CN110225540A (en) A kind of fault detection method towards centralization access net
Steinert et al. Long-term adaptation and distributed detection of local network changes
CN116390139A (en) Equipment detection method and device based on neural network, computer equipment and medium

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