删除数据中的异常值

此示例使用:

此示例显示如何使用 Hampel 滤波器检测和删除 ThingSpeak™ 通道数据中的离群值。该滤波器与中位数滤波器密切相关，它有助于在不过度平滑数据的情况下消除信号中的离群值。

使用 thingSpeakRead 函数读取存储在通道 22641 中的风速。

data = thingSpeakRead(12397,'NumPoints',80,'Fields',2,'outputFormat','table');

通过应用 hampel 函数删除原始数据中的离群值。将窗口大小指定为 6，或测量窗口两侧约三分钟的数据。此设置允许提供足够的数据来判断每个点是否是离群。由于数据是动态的，您的结果可能会根据当前天气而变化。将标准差数指定为 2。如果数据点与局部中位数的差异超过该数字，则将其替换为局部中位数。

windowSize = 6;
numMedians = 2;
[filteredData,outliers]=hampel(data.WindSpeedmph,windowSize,numMedians);

使用 thingSpeakWrite 函数将检测到的离群值数量发送到 ThingSpeak 通道。更改 channelID 和 writeAPIKey 以将数据发送到您的通道。

channelID=17504;
writeAPIKey='23ZLGOBBU9TWHG2H';
totalOutliers = sum(outliers)
thingSpeakWrite(channelID,totalOutliers,'WriteKey',writeAPIKey);

totalOutliers =

     8

另请参阅