从风速数据中删除异常值
此示例显示如何从公共通道读取数据、修改数据并显示选择元素。在示例中,您修改了 MATLAB Analysis 和 MATLAB Visualization App 提供的代码模板之一。该示例使用来自 ThingSpeak 通道 12397 的数据,该数据从位于马萨诸塞州纳蒂克的 Arduino® 气象站收集天气数据。
从模板代码创建 MATLAB Analysis 脚本
为了检测并消除纳蒂克气象站风速数据中的离群值,您可以使用代码模板编写 MATLAB® 脚本。
转至 ThingSpeak 中的 Apps 选项卡,然后选择 MATLAB Analysis。点击 New,选择 Remove outliers from wind speed data,然后点击 Create。
分析您的数据
MATLAB Code 字段已预先填充代码,用于检测并删除过去六小时内风速数据中的离群值。
1) 设置与 ThingSpeak 通信的变量。readChannelID
是从气象站收集数据的公共通道的通道 ID。windSpeedFieldID
是通道中包含风速值的字段。仅当您从私有通道读取数据时才为 readAPIkey 分配一个值。气象站是公共的,因此对于此示例,不要设置 readAPIkey。
readChannelID = 12397;
windSpeedFieldID = 2;
readAPIKey = '';
2) 使用 thingSpeakRead
函数读取过去六小时的风速值和时间戳。
[windSpeed,timeStamp] = thingSpeakRead(readChannelID,'fields',windSpeedFieldID,'NumMinutes',360,'ReadKey',readAPIKey);
3) 使用 MATLAB isoutlier
函数检查风速数据中的离群值。使用默认设置,此函数计算某个值是否与输入数据集的中位数相差超过三个缩放中位数的绝对偏差。您可以调整输入参量来自定义离群结果。识别离群数据点的索引和干净数据点的索引。使用这些索引来选择与离群数据和干净数据相对应的数据点和时间戳。
outlierDataIndex = isoutlier(windSpeed); cleanDataIndex = ~outlierDataIndex; outlierData = windSpeed(outlierDataIndex); cleanData = windSpeed(cleanDataIndex); outlierTimeStamps = timeStamp(outlierDataIndex); cleanTimeStamps = timeStamp(cleanDataIndex);
4) 使用离群数据和干净数据创建时间表,并显示离群数据点的结果。
outlierDataTable = timetable(outlierTimeStamps,outlierData);
cleanDataTable = timetable(cleanTimeStamps,cleanData);
display(outlierDataTable,'Outlier data');
3×1 timetable outlierTimeStamps outlierData ____________________ ___________ 03-Jun-2019 10:20:54 17 03-Jun-2019 13:26:14 16.6 03-Jun-2019 13:39:33 16.8
点击 Save and Run 来执行您的代码。Output 字段显示您的结果。
将数据写入通道
1) 通过将干净的数据结果写入私有通道来存储它们。要创建 ThingSpeak 通道,转至 Channels 选项卡并选择 My Channels。点击 New Channel。选中相应的复选框,并输入这些通道设置值:
Name -
Cleaned Wind Speed Measurements
Field 1 -
Wind speed (mph)
点击 Save Channel。
2) 在 MATLAB Code 字段中,设置用于写入私有通道的变量。将 writeChannelID
和 writeAPIKey
的给定值替换为您的值。您可以在页面右侧的 Channel Info 面板中找到通道 ID 和 API 密钥。
% Replace with the ID of the channel to write data to. writeChannelID = 17504; % Enter the write API key between the ''. writeAPIKey = '23ZLGOBBU9TWHG2H';
3) 将干净的风速读数及其各自的时间戳写入您的通道。
thingSpeakWrite(writeChannelID,cleanData,'timestamp',cleanTimeStamps,'Writekey',writeAPIKey);
4) 点击 Save and Run 来执行您的代码。ThingSpeak 通道中的图填充了风速时间序列数据,但没有计算出离群值。您可以通过点击页面右侧 Channel Info 面板中的通道链接来访问您的通道。
要以 CSV 格式下载数据,请点击 Data Export 按钮或 Data Import / Export 选项卡。要清除通道中所有已保存的数据,请点击 Channel Settings 选项卡。
计算移动平均值
您可以在模板中添加代码来进一步分析风速数据。除了删除离群值之外,平滑数据集的另一种方法是计算移动平均值。在这种方法中,一组局部数据点的平均值是通过整个数据集的滑动窗口来计算的。使用带有五分钟滑动窗口的 MATLAB movmean
函数来平滑风速数据。此部分未包含在代码模板中。调用 thingSpeakRead.
后,您可以将其包含在代码中
smoothData = movmean(windSpeed,minutes(5),'SamplePoints',timeStamp);
要保存您的新数据,请将其写入您的 ThingSpeak 通道。注释掉现有的 thingSpeakWrite
函数并将新的平滑时间序列数据保存到您的通道。
thingSpeakWrite(writeChannelID,smoothData,'timestamp',timeStamp,'Writekey',writeAPIKey);
要更新结果,请再次点击 Save and Run。
另请参阅
函数
thingSpeakRead
|thingSpeakWrite
|isoutlier
(MATLAB) |movmean
(MATLAB)