CN112000519A - 一种采样数据异常值修复方法及装置 - Google Patents

Abstract

本发明公开了一种采样数据异常值修复方法，通过选取原始采样值序列中采样数据并结合修复输出序列中的采样数据构建数据窗，根据函数极值原理对数据窗内采样数据进行异常数据识别，并结合数据窗内采样数据对疑似异常的数据进行修复，实现了对采样异常数据的剔除。本发明具有简单易实现，延迟小，修复精度高，无需计算门槛值和运算量低等优点，具有良好的应用前景。

Description

一种采样数据异常值修复方法及装置

技术领域

本发明涉及一种对单点采样数据异常值的检测和修复方法与装置，属于数据采集和处理技术领域。

背景技术

随着电子技术的进步，在工业过程控制、精密测量、电力自动化等领域，数字采样系统的应用越来越广泛。对模拟采样的系统而言，由于采样元件的偶发故障或者外部的电气骚扰，会使采样数据出错，对基于数字传输的远程采样系统来说，即使前端采样结果正确，也会因为通信系统的偶发性误码未被系统检出导致错误的数据。根据工程实践，对于雷击浪涌等强干扰导致的异常数据以及由通信系统误码引发的异常数据大多时候偏离正常数据集较远，比较容易识别，而对于磁场骚扰等，采样系统多表现为个别采样点轻微异常，往往会导致识别不出来，另外一方面，工业过程量和电力系统的电气量受负载扰动、系统控制或设备故障等原因往往会伴随着暂态过程，引发测量信号的暂态波动进一步增加异常数据识别难度。

对继电保护来说，电气量采样异常已经成为造成保护误动的主要原因，因此必须对采样值进行实时判断，避免异常数据引起保护装置误动、系统误操作等情况发生。而对测量系统来说，即使较小的异常数据也会影响测量精度，准确快速地剔除异常数据还将进一步提高工业控制系统的控制精度和系统响应时间。对于工业过程控制和电力系统继电保护等系统，要求异常点检测延时必须很短，但是要实现在较短的检测数据窗内准确区分飞点采样数据和系统故障时的正常暂态采样值，有很大的难度。

目前，对于采样异常点的影响已有了一些解决方法：

(1)、滤波法

在采样前端增加模拟滤波器，并降低滤波器带宽和提升阻带衰减，或者将采样数据序列经过数字滤波器，可以根据需要增加滤波器的冲激响应长度，利用数字滤波器来减少或消弱异常点的影响。

(2)、冗余采样方案

系统配置两套前端采样模块，只有两套采样模块数据正常后端应用模块才正常工作，否则进入设备故障应对模式，在单套采样模块异常时可以保证应用系统不发送误动作。

(3)、基于门槛的采样点判决

由于电磁骚扰或其偶发性采样数据异常所持续的时间一般很短，通过设置固定的门槛或动态计算的门槛，只有在采样数据窗内达到门槛的采样数据达到一定数量才可以启动应用模块功能。所以此方案适用于对响应时间不敏感的应用场景。

以上几种方法能在一定程度上减少或避免异常采样数据的影响，均存在一些不足，如针对性不强、灵敏度不高、门槛不易整定、数据窗要求较长、实时性差、系统响应时间长，无法改善异常点对精密测量的影响等问题：

(1)、无论模拟滤波器还是数字滤波器，其实质是将单个异常点的影响分摊到后续的多个采样点中，并没有针对性地剔除异常点，增加传统滤波器的滤波效果是以延长系统响应时间为代价的。此外对于一些精密测量系统来说，滤波器的影响并不能改善异常点对测量精度的影响；

(2)、对于外部电磁骚扰引发的采样异常，即使采样前端冗余也无法解决，因为此时冗余系统采样结果是一致的，而且若冗余系统中的一套系统出现偶发性采样异常必将影响系统整体响应时间。

(3)、基于门槛值的多次判别实际上延长了系统响应时间，严重影响对时间敏感的应用实时性能，如基于采样值策略超快速继电保护、基于采样值控制的电力电子系统等。

综上所述，采样数据的快速异常检测和异常值修复是当前首要解决的问题。

发明内容

本发明的目的在于克服现有技术中的不足，提供一种简单易实现，延迟小，恢复精度高，无需计算门槛值和运算量低且适合硬件加速的采样数据异常值修复方法与装置，解决采样数据的快速异常检测和异常值修复技术问题。

为达到上述目的，本发明是采用下述技术方案实现的：

第一方面，本发明提供了一种采样数据异常值修复方法，所述方法包括以下步骤：

选取原始采样值序列中采样数据并结合修复输出序列中的采样数据构建数据窗；根据函数极值原理对数据窗内采样数据进行异常数据识别；结合数据窗内采样数据对疑似异常的数据进行修复。

进一步的，所述方法还包括以下步骤：对信号进行等间隔连续采样，获得原始采样值序列{x(n)}。

进一步的，所述数据窗的构建方法包括以下步骤：

取原始采样值序列{x(n)}的最新三个采样点以及修复输出序列{y(m)}中的最新两个采样点，构成5个采样数据的数据窗，记为{x₁,x₂,x₃,x₄,x₅}，其中，x₅对应当前时刻的最新采样点x(n)，n为采样序列{x(n)}中当前最新采样点的序号，x₄对应原始采样序列中x(n)的上一采样节拍的值x(n-1)，x₃对应原始采样序列中x(n)的上两个采样节拍的值x(n-2)，而x₂对应当前时刻修复输出序列{y(m)}中上一节拍的修复采样数据y(m-1)，x₁对应修复输出序列{y(m)}中上两个采样节拍的采样数据y(m-2)，m为拟写入到修复输出序列{y(m)}的最新修复输出值序号。

进一步的，识别异常数据的方法包括以下步骤：

若(x₂-x₃)×(x₃-x₄)为非负数，则认为该所述数据窗的采样值正常，否则判别该所述数据窗的采样值为疑似异常。

进一步的，疑似异常数据的修复的方法包括如下步骤：取中间变量u的初值为x₃，发现疑似异常数据后，若(x₃-x₄)×(x₄-x₅)为负数，则按公式(1)重新计算中间变量u；

若(x₃-x₄)×(x₄-x₅)不为负数，则计算表达式[(x₁-x₂)×(x₂-x₃)]的符号，若(x₁-x₂)×(x₂-x₃)为负数，则按公式(2)更新中间变量u，否则按公式(3)更新中间变量u；

更新中间变量u到输出变量y(m)，即y(m)＝u，将y(m)作为修复输出值更新至修复输出序列{y(m)}。

进一步的，所述方法还包括以下步骤：当前数据窗数据修复完成后移动原始采样值序列的读取指针和修复输出序列的读写指针，继续对采样数据进行处理直至原始采样数据被处理完。

第二方面，本发明提供了一种采样数据异常值修复装置，所述装置包括：

数据窗模块：用于选取原始采样值序列中采样数据并结合修复输出序列中的采样数据构建数据窗；

识别模块：用于对数据窗内采样数据进行异常数据识别；

修复模块：用于结合数据窗内采样数据对疑似异常的数据进行修复。

进一步的，还包括采样模块：用于对信号进行等间隔采样，生成原始采样值序列。

第三方面，本发明还提供了一种采样数据异常值修复装置，包括处理器及存储介质和I/O接口设备；所述存储介质用于存储指令；所述I/O接口设备用于获取离线采集的或者在线采集的原始采样值序列；所述处理器用于根据所述指令进行操作以执行根据上述采样数据异常值修复方法的步骤。

第四方面，本发明上述还提供了计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时从I/O接口设备提取原始采样值序列，实现上述采样数据异常值修复方法的步骤。

与现有技术相比，本发明所达到的有益效果：

本发明通过选取原始采样值序列中采样数据并结合修复输出序列中的采样数据构建短数据窗，根据极值原理对数据窗内采样数据进行异常数据识别，并结合数据窗内采样数据对疑似异常的数据进行修复，实现了对采样异常数据的剔除，具有简单易实现，数据缓冲少、延迟小，恢复精度高，无需计算门槛值和运算量低、适合硬件加速等优点，具有良好的应用前景。

附图说明

图1是本发明的数据窗构造示意图；

图2是本发明的异常数据检测和修复流程图；

图3是采样异常数据案例1；

图4是采样异常数据案例2；

图5是采样异常数据案例3；

图6是本发明异常数据剔除和修复效果图。

具体实施方式

下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案，而不能以此来限制本发明的保护范围。

本发明的目的在于克服现有技术中的不足，提供一种简单易实现，延迟小，恢复精度高，无需计算门槛值和运算量低适合硬件加速的采样数据异常值修复方法与装置，解决采样数据的快速异常检测和异常值修复技术问题。

实施例一：

本发明的采样数据异常值修复方法，首先利用连续的三点极值检测发现疑似异常采样数据，随后进一步检测邻近三点极值检测确认异常采样点，最多三次进行疑似异常值识别，再选取邻近正确采样数据进行数据修复。

下面将结合说明书附图，从技术原理和工程实施两个方面对本发明作进一步的说明。

1.技术基础

1.1.异常检测

工程上，异常采样数据是指采样值和系统实际值不符或误差超过一定的范围表现，就局部来说，表现为采样序列的局部极大值(或极小值)，如附图3、附图4和附图5所示。而在数学上，函数极值点可能出现在间断点、边界点或导数为0点，常用判别方法如下：

若函数f(x)可导，f'(x₀)＝0，且

对于

有f'(x₀)＞0(或f'(x₀)＜0)同时

有f'(x₀)＜0(或f'(x₀)＞0)，则x₀是f(x)的极大点(或极小点)，对应的f(x₀)为函数f(x)的极大值(或极小值)。

本发明根据上述数学原理，采用差分替代微分，将采样异常点检测转变为在连续信号的采样序列{x(n)}中寻找符合公式(4)采样数据点：

[x(n-1)-x(n)][x(n)-x(n+1)]＜0 (4)

进一步地，认为x(n)为这3点之中的极大值(或极小值)，同时也是疑似异常采样数据。

1.2.数据修复

根据二阶拉格朗日插值原理，已知信号x(t)在t₀,t₁,t₂时刻的值x(t₀),x(t₁),x(t₂)，对应时刻t₀≤t≤t₂，则信号在时刻t的采样值可以表示为：

x(t)＝k_ax(t₀)+k_bx(t₁)+k_cx(t₂) (5)

其中：

进一步地，本发明在检测出异常采样数据后，剔除当前疑似异常采用数据，同时利用邻近的正常数据，通过公式(5)估计采样异常时刻的系统真实值，并以此替代异常点的异常采样数据，实现数据修复。

2.工程实施

2.1.数据窗构造

本发明实施案例基于5点的数据窗{x₁,x₂,x₃,,x₄,x₅}，其构造方法如附图1所示。从时间上来看，x₁最早，x₂,x₃,x₄次之，x₅最新。从数据来源看x₁,x₂来自修复输出序列{y(m)}，x₃,x₄,x₅对应原始采样序列{x(n)}，其中，x₅对应当前时刻的最新采样点x(n)，n为采样序列{x(n)}中当前最新采样点的序号，x₄对应原始采样序列中x(n)的上一采样节拍的值x(n-1)，x₃对应原始采样序列中x(n)的上两个采样节拍的值x(n-2)，而x₂对应当前时刻修复输出序列{y(m)}中上一节拍的修复采样数据y(m-1)，x₁对应修复输出序列{y(m)}中上两个采样节拍的采样数据y(m-2)，m为拟写入到修复输出序列{y(m)}的最新修复输出值序号。

2.2.异常检测和修复

设系统采样间隔为Ts，附图3～4给出具有一个采样异常点x3(位于时刻t)案例，其中图3a和图4a为单调下降的曲线中出现异常点，图3b和图4b为单调上升的曲线中出现异常点，而附图5a附图5b则分别对应图3a和图3b的移位。本发明根据异常点在数据窗的位置采取不同的异常点识别和修复方法。

场景(1)类似于图3a，在时刻t采样数据出现疑似异常，对应值为x₃，此时x₂,x₃,x₄构成的曲线段在t处出现极小值x₃，而以x₃,x₄,x₅构成的曲线段出现极大值为x₄。对于此种情况，本发明选取数据点x₁,x₂,x₅来估计x₃。

进一步地，参照公式5和附图3，记t₀＝0、t₁＝T_s、t₂＝4T_s，取t＝2T_s，则计算k_a、k_b、k_c得

则疑似异常x₃的估算值

可用下式计算：

场景(2)类似于图5a，在时刻t采样数据出现疑似异常，对应采样值x₃,此时由x₂,x₃,x₄构成的曲线段在t处出现极大值x₃，同时由x₁,x₂,x₃构成的曲线段出现极小值x₂。对于此种情况，本发明选取数据点x₁,x₄,x₅来估计x₃。

进一步地参照公式(5)和附图5，记t₀＝0、t₁＝3T_s、t₂＝4T_s，取t＝2T_s，则计算k_a、k_b、k_c的

则疑似异常x₃的估算值

可用下式计算：

场景(3)对于在时刻t采样数据出现疑似异常，此时由x₂,x₃,x₄构成的曲线段在t处出现极值x₃，同时由x₁,x₂,x₃构成的曲线段和x₃,x₄,x₅构成的曲线段均未出现极值，对于此种情况，本发明选取数据点x₁,x₂,x₄来估计x₃。

进一步地参照公式5和附图5，记t₀＝0、t₁＝T_s、t₂＝3T_s其t＝2T_s，则计算k_a、k_b、k_c的

k_b＝1，

则疑似异常x₃的估算值

可用下式计算：

2.3.技术实现

附图1展示了本发明5点的数据窗{x₁,x₂,x₃,,x₄,x₅}的构造方法，附图2则展示了基于单个异常点的检测与修复流程，本发明的技术实施包括以下步骤：

步骤(1)，以间隔T_s对信号进行连续采样，获得原始采样值序列{x(n)}；

步骤(2)，取信号原始采样值序列{x(n)}的最新三个采样点以及修复输出序列{y(m)}中的最新两个采样点，构成5个采样数据的数据窗，记为{x₁,x₂,x₃,x₄,x₅}，其中，x₅对应当前时刻的最新采样点x(n)，n为采样序列{x(n)}中当前最新采样点的序号，x₄对应原始采样序列中x(n)的上一采样节拍的值x(n-1)，x₃对应原始采样序列中x(n)的上两个采样节拍的值x(n-2)，而x₂对应当前时刻修复输出序列{y(m)}中上一节拍的修复采样数据y(m-1)，x₁对应修复输出序列{y(m)}中上两个采样节拍的采样数据y(m-2)，m为拟写入到修复输出序列{y(m)}的最新修复输出值序号；

步骤(3),取u为中间变量并设初值为x₃，计算表达式[(x₂-x₃)×(x₃-x₄)]的符号，若sgn[(x₂-x₃)×(x₃-x₄)]为-1则转步骤(4),否则直接转步骤(6)，其中sgn[x]表示参数x的正负号，若参数为负则输出-1，若参数为正则输出+1，若参数为0则输出0；

步骤(4),计算表达式[(x₃-x₄)×(x₄-x₅)]的符号，若sgn[(x₃-x₄)×(x₄-x₅)]为-1，则先按公式(9)更新中间变量u后再转步骤(6)，否则直接转步骤(5)：

步骤(5),计算表达式[(x₁-x₂)×(x₂-x₃)]的符号，若sgn[(x₁-x₂)×(x₂-x₃)]为-1，则先按公式(10)更新中间变量u后再转步骤(6)，否则按公式(11)更新中间变量后转步骤(6)；

步骤(6)，更新中间变量u到输出变量y(m)，即y(m)＝u，将y(m)作为修复输出值更新至修复输出序列{y(m)}，并转步骤(2)继续检查新的采样数据。

随着采样的进行，算法将获得原始采样序列的修复输出序列{y(m)},且修复后输出序列{y(m)}在时间上相比原始输入序列{x(n)}滞后两个采样间隔。

附图6展示了本发明一个异常检测与修复的效果，图6a为一个包含异常数据的正弦采样序列，图6b为经过本发明的修复处理之后的正弦采样序列，进一步证明了本发明的实施效果。

综上所述，本发明结合原始采样序列{x(n)}和修复输出序列{y(m)}构建了5点数据窗，通过数据窗内连续3点间的差分符号变化来识别疑似异常点，并合理选择邻近的3个点通过插值完成疑似异常数据点修复，最多三次进行疑似异常值识别。本发明具有简单易实现，同时具有数据缓冲少，输出延迟小，异常数据恢复精度高，无需动态计算门槛，运算量低，适合硬件加速等优点，具有良好的应用前景。

实施例二：

本发明提供了一种采样数据异常值修复装置，能够用于实现实施例一所述的方法步骤，所述装置包括：

数据窗模块：用于选取原始采样值序列中采样数据并结合修复输出序列中的采样数据构建数据窗；

识别模块：用于对数据窗内采样数据进行异常数据识别；

修复模块：用于结合数据窗内采样数据对疑似异常的数据进行修复。

还包括采样模块：用于对信号进行等间隔采样，生成原始采样值序列；

实施例三：

本发明还提供了一种采样数据异常值修复装置，包括处理器、存储介质和和I/O接口设备；所述存储介质用于存储指令；所述I/O接口设备用于获取离线采集的或者在线采集的原始采样值序列；所述处理器用于根据所述指令进行操作以执行根据下述方法的步骤：

步骤(1)，通过I/O接口设备获取离线采集的原始采样值序列{x(n)}，或者以间隔T_s在线对信号进行连续采样，获得原始采样值序列{x(n)}；

步骤(2)，取信号原始采样值序列{x(n)}的最新三个采样点以及修复输出序列{y(m)}中的最新两个采样点，构成5个采样数据的数据窗，记为{x₁,x₂,x₃,x₄,x₅}，其中，x₅对应当前时刻的最新采样点x(n)，n为采样序列{x(n)}中当前最新采样点的序号，x₄对应原始采样序列中x(n)的上一采样节拍的值x(n-1)，x₃对应原始采样序列中x(n)的上两个采样节拍的值x(n-2)，而x₂对应当前时刻修复输出序列{y(m)}中上一节拍的修复采样数据y(m-1)，x₁对应修复输出序列{y(m)}中上两个个采样节拍的采样数据y(m-2)，m为拟写入到修复输出序列{y(m)}的最新修复输出值序号；

步骤(3),取u为中间变量并设初值为x₃，计算表达式[(x₂-x₃)×(x₃-x₄)]的符号，若sgn[(x₂-x₃)×(x₃-x₄)]为-1则转步骤(4),否则直接转步骤(6)，其中sgn[x]表示参数x的正负号；

步骤(4),计算表达式[(x₃-x₄)×(x₄-x₅)]的符号，若sgn[(x₃-x₄)×(x₄-x₅)]为-1，则先按公式(12)更新中间变量u后再转步骤(6)，否则直接转步骤(5)：

步骤(5),计算表达式[(x₁-x₂)×(x₂-x₃)]的符号，若sgn[(x₁-x₂)×(x₂-x₃)]为-1，则先按公式(13)更新中间变量u后再转步骤(6)，否则按公式(14)更新中间变量后转步骤(6)；

步骤(6)，更新中间变量u到输出变量y(m)，即y(m)＝u，将y(m)作为修复输出值更新至修复输出序列{y(m)}，并转步骤(2)继续检查新的采样数据。

随着采样的进行，算法将获得原始采样序列的修复输出序列{y(m)},且修复后输出序列{y(m)}在时间上相比原始输入序列{x(n)}滞后两个采样间隔。

实施例四：

本发明实施例还提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现下述方法的步骤：

步骤(1)，从I/O接口设备提取原始采样值序列{x(n)}；

步骤(2)，取信号原始采样值序列{x(n)}的最新三个采样点以及修复输出序列{y(m)}中的最新两个采样点，构成5个采样数据的数据窗，记为{x₁,x₂,x₃,x₄,x₅}，其中，x₅对应当前时刻的最新采样点x(n)，n为采样序列{x(n)}中当前最新采样点的序号，x₄对应原始采样序列中x(n)的上一采样节拍的值x(n-1)，x₃对应原始采样序列中x(n)的上两个采样节拍的值x(n-2)，而x₂对应当前时刻修复输出序列{y(m)}中上一节拍的修复采样数据y(m-1)，x₁对应修复输出序列{y(m)}中上两个采样节拍的采样数据y(m-2)，m为拟写入到修复输出序列{y(m)}的最新修复输出值序号；

步骤(3),取u为中间变量并设初值为x₃，计算表达式[(x₂-x₃)×(x₃-x₄)]的符号，若sgn[(x₂-x₃)×(x₃-x₄)]为-1则转步骤(4),否则直接转步骤(6)，其中sgn[x]表示参数x的正负号；

步骤(4),计算表达式[(x₃-x₄)×(x₄-x₅)]的符号，若sgn[(x₃-x₄)×(x₄-x₅)]为-1，则先按公式(15)更新中间变量u后再转步骤(6)，否则直接转步骤(5)：

步骤(5),计算表达式[(x₁-x₂)×(x₂-x₃)]的符号，若sgn[(x₁-x₂)×(x₂-x₃)]为-1，则先按公式(16)更新中间变量u后再转步骤(6)，否则按公式(17)更新中间变量后转步骤(6)；

步骤(6)，更新中间变量u到输出变量y(m)，即y(m)＝u，将y(m)作为当前修复输出值更新至修复输出序列{y(m)}，并转步骤(2)继续检查新的采样数据。

随着原始采样数据的读入，算法将获得原始采样序列的修复输出序列{y(m)},且修复后输出序列{y(m)}在时间上相比原始输入序列{x(n)}滞后两个采样间隔。

本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。

本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形也应视为本发明的保护范围。

Claims (10)

1.一种采样数据异常值修复方法，其特征在于，所述方法包括以下步骤：

根据原始采样值序列中采样数据并结合修复输出序列中的采样数据构建数据窗；

根据函数极值原理对数据窗内采样数据进行异常数据识别；

通过数据窗内的采样数据对异常数据进行修复。

2.根据权利要求1所述的一种采样数据异常值修复方法，其特征在于，所述方法还包括以下步骤：

对信号进行等间隔连续采样，获得原始采样值序列{x(n)}。

3.根据权利要求1所述的一种采样数据异常值修复方法，其特征在于，所述数据窗的构建方法包括以下步骤：

取原始采样值序列{x(n)}的最新三个采样点以及修复输出序列{y(m)}中的最新两个采样点，构成5个采样数据的数据窗，记为{x₁,x₂,x₃,x₄,x₅}，其中，x₅对应当前时刻的最新采样点x(n)，n为采样序列{x(n)}中当前最新采样点的序号，x₄对应原始采样序列中x(n)的上一采样节拍的值x(n-1)，x₃对应原始采样序列中x(n)的上两个采样节拍的值x(n-2)，而x₂对应当前时刻修复输出序列{y(m)}中上一节拍的修复采样数据y(m-1)，x₁对应修复输出序列{y(m)}中上两个采样节拍的修复采样数据y(m-2)，m为拟写入到修复输出序列{y(m)}的最新修复输出值序号。

4.根据权利要求3所述的一种采样数据异常值修复方法，其特征在于，识别异常数据的方法包括以下步骤：

若(x₂-x₃)×(x₃-x₄)为负数，则认为所述数据窗内采样值异常，否则判别所述数据窗内采样值正常，并把x₃更新到输出变量y(m)，即y(m)＝x₃，同时将y(m)更新至修复输出序列{y(m)}。

5.根据权利要求4所述的一种采样数据异常值修复方法，其特征在于，异常数据的修复的方法包括如下步骤：

取中间变量u并设置初值为x₃，

发现异常数据后，若(x₃-x₄)×(x₄-x₅)为负数，则按公式(1)重新计算中间变量u；

若(x₃-x₄)×(x₄-x₅)不为负数，则计算表达式[(x₁-x₂)×(x₂-x₃)]的符号，若(x₁-x₂)×(x₂-x₃)为负数，则按公式(2)更新中间变量u，否则按公式(3)更新中间变量u；

更新中间变量u到输出变量y(m)，即y(m)＝u，将y(m)作为修复输出值更新至修复输出序列{y(m)}。

6.根据权利要求1所述的一种采样数据异常值修复方法，其特征在于，所述方法还包括以下步骤：当前数据窗数据修复完成后移动原始采样值序列的读取指针和修复输出序列的读写指针，继续对采样数据进行处理直至原始采样数据被处理完。

7.一种采样数据异常值修复装置，其特征在于，所述装置包括：

数据窗模块：用于选取原始采样值序列中采样数据并结合修复输出序列中的采样数据构建数据窗；

识别模块：用于对数据窗内采样数据进行异常数据识别；

修复模块：用于结合数据窗内采样数据对疑似异常的数据进行修复。

8.根据权利要求7所述的一种采样数据异常值修复装置，其特征在于，还包括采样模块：用于对信号进行等间隔采样，生成原始采样值序列。

9.一种采样数据异常值修复装置，其特征在于，包括处理器及存储介质和I/O接口设备；

所述存储介质用于存储指令；

所述I/O接口设备用于获取离线采集的或者在线采集的原始采样值序列；

所述处理器用于根据所述指令进行操作，从I/O接口设备提取原始采样值序列，以执行根据权利要求1～6任一项所述方法的步骤。

10.计算机可读存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现权利要求1～6任一项所述方法的步骤。

Images