对用电负荷时间序列数据进行K-medoids聚类建模和GAM回归|附数据代码

简介: 对用电负荷时间序列数据进行K-medoids聚类建模和GAM回归|附数据代码


通过对用电负荷的消费者进行聚类,我们可以帮助客户提取典型的负荷曲线,提高后续用电量预测的准确性,检测异常或监控整个智能电网(Laurinec等人(2016),Laurinec和Lucká( 2016))点击文末“阅读原文”获取完整代码数据


第一个用例通过K-medoids聚类方法提取典型的电力负荷曲线。

相关视频

image.png

image.png


有50个长度为672的时间序列(消费者),长度为2周的耗电量的时间序列。这些测量数据来自智能电表。

维数太高,会发生维数诅咒。因此,我们必须以某种方式降低维度。最好的方法之一是使用时间序列表示,以减少维数,减少噪声并提取时间序列的主要特征。

对于用电的两个季节性时间序列(每日和每周季节性),基于模型的表示方法是提取典型用电量的最佳方法。

让我们使用一种基于模型的基本表示方法- 平均季节性。在此还有一个非常重要的注意事项,对时间序列进行归一化是对时间序列进行每次聚类或分类之前的必要步骤。我们想要提取典型的消耗曲线,而不是根据消耗量进行聚类。



维数上已大大降低。现在,让我们使用K-medoids聚类方法来提取典型的消耗量。由于我们不知道要选择合适的簇数,即先验信息,因此必须使用验证指数来确定最佳簇数。我将使用Davies-Bouldin指数进行评估,通过Davies-Bouldin指数计算,我们希望找到其最小值。




Kaizong Ye

拓端分析师



我将聚类数的范围设置为2-7。

让我们绘制评估的结果。

74e31087a421a74c417d3cc2e43d186a.png

聚类的“最佳”数目是7。

我们绘制有7个聚类的聚类结果。

4cbc2765ef2a00170420930dbd35db8f.png

我们可以看到5个典型的提取轮廓 (簇的中心)。接下来的两个簇可以称为离群值

现在,让我们尝试一些更复杂的方法来提取季节 GAM回归系数。我们可以提取每日和每周的季节性回归系数 。

## \[1\] 50 53

由于GAM方法中使用样条曲线 。让我们对数据进行聚类并可视化其结果。

让我们绘制 评估的结果。

dd7b8b959c13e72d6d4f1e61832dc265.png

聚类的最佳数目为7。让我们绘制结果。

9a5f15cfe2639a828bdc9eb0131c7e74.png

提取的消费数据比平均季节性数据更平滑。现在,K 中心提取了4个典型的轮廓,并确定了3个簇。

我展示了一些自适应表示的聚类结果,让我们以DFT(离散傅立叶变换)方法为例,并提取前48个DFT系数。

dim(data_dft)
## \[1\] 50 48

让我们绘制评估的结果。


e8e7797feed675d752c2d419f21e8129.png

我们可以在4个簇中看到“肘部”。

94ef6f50ddf006a9cf53b5564b4e074f.png 这些结果可以较好解释。因此,基于模型的时间序列表示在此用例中非常有效 。

建议在每天的时间序列中使用与FeaClip一起的窗口方法。最大的优点是不需要与FeaClip方法一起进行标准化。

 dim(data_feaclip)
## \[1\]  50 112

让我们绘制评估的结果。

0a5534a2b2cc4d3b118ad3d063111fe2.png 我们可以看到现在出现了2个“肘部”。最大的变化是在2到3之间,因此我选择3。

fc5027a757d7f0f5f3cee0c1d73bdb61.png 可分离性好于DFT。但是也可以检查具有不同数量聚类的其他结果。

结论

在本教程中,我展示了如何使用时间序列表示方法来创建用电量的更多特征。然后,用时间序列进行K-medoids聚类,并从创建的聚类中提取典型的负荷曲线。




相关文章
|
6月前
|
机器学习/深度学习 数据可视化 vr&ar
数据分享|SAS与eviews用ARIMA模型对我国大豆产量时间序列预测、稳定性、白噪声检验可视化
数据分享|SAS与eviews用ARIMA模型对我国大豆产量时间序列预测、稳定性、白噪声检验可视化
|
6月前
R语言分布滞后线性和非线性模型DLM和DLNM建模应用| 系列文章
R语言分布滞后线性和非线性模型DLM和DLNM建模应用| 系列文章
|
6月前
|
前端开发 数据建模 计算机视觉
R语言主成分回归(PCR)、 多元线性回归特征降维分析光谱数据和汽车油耗、性能数据
R语言主成分回归(PCR)、 多元线性回归特征降维分析光谱数据和汽车油耗、性能数据
|
6月前
|
计算机视觉
偏最小二乘回归(PLSR)和主成分回归(PCR)分析光谱数据
偏最小二乘回归(PLSR)和主成分回归(PCR)分析光谱数据
|
6月前
|
存储 对象存储
R语言分布滞后线性和非线性模型(DLNM)分析空气污染(臭氧)、温度对死亡率时间序列数据的影响
R语言分布滞后线性和非线性模型(DLNM)分析空气污染(臭氧)、温度对死亡率时间序列数据的影响
|
6月前
|
机器学习/深度学习 算法 数据挖掘
R语言气象模型集成预报:神经网络、回归、svm、决策树用环流因子预测降雨降水数据
R语言气象模型集成预报:神经网络、回归、svm、决策树用环流因子预测降雨降水数据
|
6月前
|
监控 数据可视化 数据挖掘
对用电负荷时间序列数据进行K-medoids聚类建模和GAM回归
对用电负荷时间序列数据进行K-medoids聚类建模和GAM回归
|
6月前
|
监控 数据可视化 数据挖掘
R语言对用电负荷时间序列数据进行K-medoids聚类建模和GAM回归
R语言对用电负荷时间序列数据进行K-medoids聚类建模和GAM回归
|
6月前
|
存储 数据库
R语言用多项式回归和ARIMA模型预测电力负荷时间序列数据
R语言用多项式回归和ARIMA模型预测电力负荷时间序列数据
|
6月前
|
存储
R语言分布滞后线性和非线性模型(DLMs和DLNMs)分析时间序列数据
R语言分布滞后线性和非线性模型(DLMs和DLNMs)分析时间序列数据