一种自适应确定聚类数的电力负荷数据加权增量聚类方法

1.本发明涉及电力负荷加权增量聚类方法，具体地说是将自适应确定聚类数的dcs统计模型与自相关加权增量模糊c均值聚类算法相结合的电力负荷增量聚类方法。


背景技术：

2.用户用电负荷数据聚类是用户用电负荷特性建模的核心内容，将电力负荷聚类用于用户细分能够有效解决用户用电负荷特性随机时变性问题和区域分散性问题，开展用户用电负荷数据聚类研究具有重要的现实意义。zhang等人在论文“electricity consumption pattern analysis beyond traditional clustering methods：a novel self-adapting semi-supervised clustering method and application case study”提出一种新的基于自适应度量学习过程的半监督自动聚类方法，通过对负荷字典和用户日常用电量进行分析，为样本数据集中用户构建了完整的用户用电行为模式。cn112215490a公开了一种基于相关性系数改进的k-means的电力负荷聚类分析方法，其利用小波变换进行数据预处理，然后采用主成分分析对数据进行降维，最后利用皮尔逊相关系数对电力负荷曲线进行聚类。此方法存在对时变数据隐含信息难以挖掘的问题。kaur等人在论文“behavior segmentation of electricity consumption patterns：a cluster analytical approach”中，分别从每个用户用电时间序列中提取用户行为特征，然后采用高斯混合聚类算法进行用户用电聚类研究，从而将用户划分为不同的集群并提取用户用电消费模式。cn108805213b公开了一种涉及小波熵降维的电力负荷曲线双层聚类方法，其利用小波熵值对数据进行降维，结合双层谱聚类获得形态相似和精细化的负荷类簇。该方法存在初始聚类数难以设定的问题。liu等人在论文“a moving shape-based robust fuzzy k-modes clustering algorithm for electricity profiles”提出了一种新的基于运动形状的鲁棒模糊k-modes聚类方法，旨在准确识别电力负荷时间序列中的形状模式，实验表明所提方法能够有效捕捉电力用户负荷使用模式。cn111860634a涉及一种基于ochnn-k-means算法的负荷聚类方法，包括数据的采集与处理，利用elbow method和交叉验证获得最佳的聚类个数，通过ochnn方法获得初始类中心，最后进行k-means聚类方法实现负荷数据划分。该方法存在对噪声点敏感性问题。
3.电力信息采集技术的不断创新发展，使得电力系统数据来源和种类多样化，数据量急剧增加，所获得的海量数据多具备时序性、高维性和高波动性等大数据特征。传统电力负荷聚类算法已难以适用，人们尝试将统计分析模型和聚类分析算法相结合。aslan等人在论文“robust-learning fuzzy c-means clustering algorithm with unknown number of clusters”提出一种基于阈值自回归模型的时间序列聚类方法，其主要通过拟合线性自回归模型和非线性阈值自回归模型，对每个序列的数据生成机制(dgm)进行逼近，并将模型估计输出用于特征提取，根据模型估计得到的特征向量进行聚类，可用于确定和监控一组跨时间同步移动的时间序列变量。cn113780343a公开提出一种基于lttb降维的双边斜率dwt距离负荷谱聚类方法，包括数据的采集与处理，利用lttb进行降维，计算双边斜率距离
并作为dwt的相似性度量，对处理后的数据集进行dwt谱聚类算法。即使lttb降维算法优于平均降采样算法，但由于电力负荷数据具有时变性和异方差性，存在难以从不断增长的电力负荷数据量中有效挖掘细节信息。cn109272058b公开了一种集成电力负荷曲线聚类方法，首先利用som神经网络对数据进行粗聚类，随后利用dbscan进行再聚类。存在对时序数据特征提取处理不足、难以根据有效的时序信息进行聚类分析。otranto在论文“fuzzy clustering with the entropy of attribute weights”提出一种基于模型的异方差时间序列三级聚类分析算法，分别根据时间序列的无条件波动率、时变波动率参数以及garch模型相应参数依次分级完成聚类，第一级根据无条件波动率对时间序列进行分类，第二级在第一级分类基础上对具有相似时变波动率的时间序列进行分类，最后在前两级分类基础上区分具有相同garch模型参数的时间序列，从而获得时间序列数据的准确分类。


技术实现要素：

4.本发明所要解决的技术问题是：提供一种自适应确定聚类数的dcs统计模型和自相关加权增量模糊c均值聚类算法相结合的电力负荷增量聚类方法，电力负荷数据是同时具有时变性、高维性以及异方差性等特点的时间序列数据，在自适应确定聚类数的基础上，建立一种基于高斯分布数据观测驱动的动态条件评分(dcs)统计模型，利用dcs统计模型参数的条件矩估计计算得到自相关值数据集，并通过加权模糊c均值聚类算法(irfcm)获得数据权重信息，将增量学习聚类方法与加权模糊c均值聚类算法相结合，再结合不同用电负荷数据流时序特性，构建自相关加权增量模糊c均值聚类算法(awi-fcm)实现聚类分析，能够在不依赖任何聚类算法下，实现电力负荷数据集上最佳聚类数的选取，并且通过对数据点分配恰当的权重值减少噪声点对聚类中心的影响，从而加强聚类的鲁棒性，提高聚类的准确率，能够有效解决电力负荷数据存在的高维性和异方差特性所导致的聚类鲁棒性低、聚类准确率低的问题，同时重点解决新的隐藏信息随着时间的推移不断出现时的有效聚类问题。
5.本发明解决该技术问题所采用的技术方案是：将自适应确定聚类数dcs统计模型与自相关加权增量模糊c均值聚类算法相结合的电力负荷增量聚类方法，针对用户用电负荷时间序列数据的时变性、高维性以及异方差性和其在聚类过程中存在的初始聚类数难设定，对噪声点敏感等问题，采用i-nice算法自适应确定聚类数，建立一种基于高斯分布数据观测驱动的动态条件评分(dcs)统计模型，利用统计模型参数的条件矩估计计算得到自相关值数据集，并通过加权模糊c均值聚类算法(irfcm)获得权重信息，将增量学习聚类方法与加权模糊c均值聚类算法相结合，再结合不同用电负荷数据流时序特性，构建自相关加权增量模糊c均值聚类算法(awi-fcm)，具体步骤如下：
6.第一步，用电用户日负荷数据预处理与分析：
7.获取k个用户的同年同月最大负荷日的日负荷数据集，对数据集进行数据预处理，包括对缺失值的查找与填补和对异常值的检测与修正，得到的新数据集作为电力负荷时间序列模型建立以及聚类分析所用数据，绘制用户日负荷曲线，通过对日负荷曲线的波动情况进行观察分析，对负荷曲线进行初步分类。
8.第二步，自适应确定聚类数：
9.通过i-nice算法确定聚类数并不依赖任何聚类算法，其通过对原始数据集的数据
分布等特征进行分析得到数据集最佳聚类数：
10.第2.1步，设观测点个数为q，观测点集合为p＝{p1，p2，...，pq}，数据集为n行m列的矩阵数据点个数为n，任意观测点到数据集的距离集合为y
p
＝{y1，y2，...yn}，且满足yk≥0，k＝1，2，...n，假设距离分布为混合伽马模型gamma mixture model(gmm)。
11.距离集合y
p
的gmm模型定义为：
[0012][0013]
其中θ为模型的参数向量，θj包含了形状参数αj和尺度参数βj，πj是混合比例，每个分量的密度函数为单伽马分布的密度函数：
[0014][0015]
其中形状参数αj和尺度参数βj应满足αj＞0，βj＞0，混合比例应满足：伽马函数γ(x)的定义为：
[0016][0017]
第2.2步，对于距离集合y
p
的gmm模型，运用em算法通过最大化对数似然函数求解模型参数，其对数似然函数的表达式为：
[0018][0019]
第2.3步，引入随机变量z＝{zi}，zi表示y
p
中的yi被分配给第j个伽马成分，在e步中根据初始化或上一步参数值计算y
p
中yi属于第j个成分的概率p(zi＝j|yi，θn)为：
[0020][0021]
对数似然函数关于随机变量z＝{zi}的期望值计算公式为式(5)，在m步中通过最大化对数似然函数的期望值求解下一次迭代中参数θ
n+1
的估计值，如式(6)所示，
[0022][0023][0024]
第2.4步，运用拉格朗日乘子法求解期望最大化对数似然函数的参数值，得到混合比例的参数估计值其计算公式为：
[0025][0026]
形状参数和尺度参数的估计值和为：
[0027][0028][0029]
其中为digamma函数，式(9)是非线性方程，没有封闭形式的解，运用牛
顿法求解参数的值。
[0030]
第2.5步，每个观测点建立mmax个gmm模型，运用二阶akaik信息准则(aicc)计算指标值进行模型的选择，aicc值达到最大时的gmm成分数即为最佳的聚类数。
[0031]
aicc的计算公式为：
[0032][0033]
其中l(θ
*
)为对数似然函数的最大值，n为数据点个数，q为参数的数量，q＝3m，通过计算aicc的值确定每个观测点处的最佳gmm模型，比较多个观测点的最佳gmm模型得到最终的聚类数。
[0034]
第三步，确定dcs的高斯分布用户假设：
[0035]
随机选取初步分类的每个类别中的用电用户，绘制直方图以及qq图(quantile-quantile plot)进行可视化，由图可知，qq图实际上是一个散点图，能够显示数据集与所选理论分布的拟合程度，若日负荷时间序列数据的散点位于红线附近，则说明观察到的向量近似服从高斯分布，将高斯分布(12)作为用户数据y
t
的基础假设：
[0036][0037]
其中μ
t
为时变参数均值，为时变参数方差；
[0038]
第四步，建立dcs模型，并利用极大似然估计法估计模型参数：
[0039]
建立dcs(p，q)模型，令p＝q＝1，用电用户数据得到基于高斯分布假设条件下的dcs(1，1)模型，也即f
t
第t个实现的更新方程，如式(13)所示：
[0040]ft
＝ω+as
t-1
+bf
t-1
ꢀꢀꢀꢀꢀꢀ
(13)
[0041]
其中时变参数向量常数向量实矩阵实矩阵其中标量参数包括ω
μ
、a
μ
等是静态参数向量θ的函数，对ω，a和b运用极大似然估计方法进行参数估计；
[0042]
第五步，计算模型的相关参数，得到最终的dcs模型：
[0043]
根据t-1时刻的驱动向量s
t-1
及时变参数向量f
t-1
计算得到dcs模型中f
t
，其中驱动向量
[0044]
第5.1步，计算条件得分向量其计算公式如式(14)(15)所示：
[0045][0046][0047]
第5.2步，计算缩放矩阵s
t
，其计算公式如式(16)(17)所示：
[0048][0049][0050]
第5.3步，得到时变参数更新方程，如式(18)(19)所示：
[0051]
[0052][0053]
第5.4步，定义无条件矩为：
[0054]
κ＝(i
n-b)-1
ω
ꢀꢀꢀꢀꢀꢀꢀ
(20)
[0055]
其中in为单位矩阵，将参数b和ω的估计值带入上述公式即得到时变参数均值和方差的无条件矩的估计值。无条件矩表示长期平均值，随着时间的推移，波动将达到一个平均值，通过对无条件矩进行计算能够对时间序列数据进行长期分布的分析。
[0056]
将dcs的参数估计值和相应的参数带入时变参数的更新方程得到时变均值和时变方差的条件矩估计值，根据式(20)得到无条件矩估计值，将条件矩估计值和无条件矩估计值分别作为聚类算法的输入数据集。
[0057]
至此，完成电力负荷数据的自适应确定聚类数的dcs模型建立。
[0058]
时间序列电力负荷数据基于自适应确定聚类数的dcs模型的awi-fcm聚
[0059]
类分析：
[0060]
第六步，根据dcs模型计算时间序列的自相关值：
[0061]
在动态条件评分模型中，对每个模型的第r(r＝1，2，...，r)个条件矩估计进行r级聚类。dcs(1，1)模型意味着r＝2，因此r＝1和r＝2的时间序列的条件矩估计分别意味着dcs(1，1)模型的参数和的估计。根据第r个条件矩估计来获得第r个估计自相关值将时变均值和时变方差的条件矩估计值作为计算自相关值的数据集分别进行计算，也就是计算时间序列y
t
的第r个条件矩在滞后l时的估计自相关值，计算公式如(21)所示：
[0062][0063]
其中是第k个时间序列从时刻t到时刻t-l过程的r阶条件矩的均值，时间序列k和k
′
基于自相关的r阶条件矩之间的距离为：
[0064][0065]
第七步，根据所得的序列自相关估计值数据集进行irfcm聚类：
[0066]
基于rfcm聚类算法进行改进得到加权模糊c均值聚类算法(irfcm)，其同时考虑数据点和数据特征的不同权重，通过加入可调参数α避免数据集的数据个数对权重信息和聚类中心产生影响，具体算法如下：
[0067]
第7.1步，设数据集为k行t列的矩阵数据点个数为k，数据特征个数为t，聚类算法的目标函数和约束条件为：
[0068][0069][0070][0071]
其中，聚类数为c，vi为第i个聚类中心，u
ik
为第k个数据点xk对第i个聚类中心vi的隶属度，r
ij
为第j个数据点对第i个聚类类别的特征权重，w
ik
为第k个数据点xk对第i个聚类类别的数据权重。
[0072]
第7.2步，利用数据特征权重和数据权重的l2范数正则化协调权重分散。
[0073]
其中，η是可调参数，δ和ε是正则化参数，通过调整δ的值使得更多数据特征信息参与聚类中心的形成，使得聚类算法聚类准确性更高，α参数用以避免大数据集的数据个数对权重信息产生影响，当数据集数据个数k≤40时，数据个数对权重信息影响较小，当数据个数k＞40时，数据个数将影响产生权重矩阵的精确度，权重矩阵的分辨率较低，无法进行有效聚类，通过调整α的值提高权重矩阵的分辨率和精确度。
[0074]
第7.3步，利用拉格朗日乘子法求解目标函数和约束条件，得到关于模糊隶属度u
ik
、数据特征权重r
ij
、数据权重w
ik
、聚类中心v
ij
的迭代公式：
[0075][0076][0077][0078][0079]
第八步，根据irfcm聚类得到的数据权重进行基于dcs模型的awi-fcm：
[0080]
将a-wfcm算法作为基础，根据所得自相关估计值数据集进行irfcm聚类，得到数据权重w
ik
，将权重信息w
ik
和数据集作为增量聚类的算法基础，得到了所提出的自适应加权增量模糊c均值聚类算法，具体算法如下：
[0081]
第8.1步，将k个时间序列数据点随机分为p(p＝1，2，...，p)个数据块，每个数据块有k/p个数据点，对第一个数据块p＝1执行a-wfcm聚类。
[0082]
第8.1.1步，根据i-nice算法得到的聚类数c随机初始化隶属度矩阵，使得其满足以下约束条件：
[0083][0084]
第8.1.2步，根据得到的隶属度矩阵计算聚类中心，其计算公式为：
[0085][0086]
第8.1.3步，根据隶属度矩阵和聚类中心计算目标函数j的值，其计算公式如式(31)所示，若前后两次j值之差小于规定的阈值ε，则结束迭代，输出隶属度矩阵和聚类中心，否则进行下一步；
[0087][0088]
第8.1.4步，根据聚类中心重新计算得到隶属度矩阵并返回步骤8.1.2，隶属度矩阵的计算公式为：
[0089][0090]
第8.2步，重新根据式(33)赋予计算得到的数据块聚类中心权重，将该数据块聚类中心和计算得到的权重加入下一个数据块中进行a-wfcm聚类(步骤8.1.1至8.1.4)；
[0091][0092]
其中，n为当前数据块中数据点的个数，j为上一个数据块中聚类中心的个数。
[0093]
第8.3步，重复步骤8.2直到所有数据块完成聚类，最后根据最终的聚类中心结果根据式(32)重新计算所有数据点的隶属度。
[0094]
至此，完成电力负荷数据的基于自适应确定聚类数的dcs模型的自相关加权增量模糊c均值聚类。
[0095]
上述第二步中通过i-nice算法自适应确定聚类数，该算法并不依赖任何聚类算法，其通过对原始数据集的数据分布等特征进行分析得到数据集最佳聚类数。假设数据点到数据集的距离分布为混合伽马模型，通过em算法计算其中参数，运用二阶akaik信息准则(aicc)计算指标值，当aicc值达到最大时的gmm成分数即为最佳的聚类数。有效解决了传统聚类算法存在的难以设定初始化聚类数目的问题。
[0096]
上述第三步中时间序列用户数据y
t
中t的值取决于测量电力负荷数据的分辨率，时变参数是根据每个时间点t以及原始时间序列y
t
进行计算得到，绘制初步分类中典型用户的直方图和qq图来确定dcs模型的用户假设，将高斯分布作为全部用户的dcs模型的概率分布假设，用来描述多种情况下的随机变量，具有较高的适普性。
[0097]
上述第四步中dcs(p，q)模型令p＝q＝1时，已可以在模型复杂度较低的情况下描述时间序列数据，根据高斯分布用户假设的参数确定时变参数向量f
t
的向量表达式，同时确定后续待估计参数的向量表达式。
[0098]
上述第六步中计算时间序列y
t
的第3r个条件矩在滞后l时的估计自相关值时，滞后的阶数l通过对原时间序列的时序图和自相关图显示的特征进行分析得到。
[0099]
上述第七步中利用加权模糊c均值聚类算法(irfcm)，其考虑数据点的不同权重，并加入可调参数α避免数据集的数据个数对权重信息和聚类中心产生影响，利用拉格朗日乘子法得到权重信息。
[0100]
本发明的有益效果是：与现有技术相比，本发明的突出的实质性特点和显著进步如下：
[0101]
(1)本发明方法与cn112215490a相比，本发明方法所具有的优点是能够有效挖掘时间序列电力负荷数据中的时变性来建立统计模型并进行聚类分析。
[0102]
(2)本发明方法与cn108805213b相比，本发明方法所具有的优点是能够自适应确定最佳聚类数目，提高聚类效率。
[0103]
(3)本发明方法与cn111860634a相比，本发明方法所具有的优点是通过对数据点进行一个权重分配，能够降低噪声点对聚类效果的影响。
[0104]
(4)本发明方法与cn113780343a相比，本发明方法所具有的优点是能够针对时间序列时变性和异方差性，重点捕捉时间序列数据的细节信息，能够准确地进行聚类分析。
[0105]
(5)本发明方法与cn109272058b相比，本发明方法所具有的优点是可以充分提取时间序列数据的时变特征，可以对动态更新的用户数据将新增加的用户加入原数据中进行聚类分析。
[0106]
(6)本发明构建了一种自适应确定聚类数的dcs模型，首先利用i-nice算法自适应确定聚类数，其假设数据集的距离分布为混合伽马模型，通过em算法计算其中参数，运用二阶akaik信息准则(aicc)计算指标值，当aicc值达到最大时的gmm成分数即为最佳的聚类
数。其次建立dcs模型，其利用似然函数的标度分数随时间更新参数，利用时间序列的数据特征为在广泛的非线性模型中引入时变参数提供了一个统一的框架。以自适应确定聚类数的dcs模型对时间序列数据进行有效的数据分析，再进行聚类分析，提高了时间序列数据聚类的效率和准确率。
[0107]
(7)本发明方法针对初始聚类数目难以设定的问题，首先自适应确定聚类数，接着根据用户用电日负荷时间序列数据的时变性、高维性以及异方差等特性，建立基于高斯分布数据观测驱动的动态条件评分(dcs)统计模型，为了描述时间序列电力负荷数据的前后关系，节约计算成本，利用dcs模型参数的条件矩估计计算得到自相关值数据集，根据不同的数据点对聚类中心有着不同的重要性，进行irfcm聚类，最后结合增量聚类，构建自相关加权增量模糊c均值聚类算法(awi-fcm)，因此，本发明方法是一种自适应确定聚类数的dcs统计模型与自相关加权增量模糊c均值聚类算法相结合的电力负荷增量聚类方法，能够在不依赖任何聚类算法的情况下，实现电力负荷数据集上最佳聚类数的选取，并且通过对数据点分配恰当的权重值能够减少噪声点对聚类中心的影响，从而加强聚类的鲁棒性，提高聚类的准确率，有效解决电力负荷数据存在的高维性和异方差特性所导致的聚类效率低、聚类准确率低的问题，重点解决新的隐藏信息随着时间的推移不断出现时的有效聚类问题。
附图说明
[0108]
下面结合附图和实施例对本发明进一步说明。
[0109]
图1是本发明自适应确定聚类数的dcs统计模型与自相关加权增量模糊c均值聚类算法相结合的电力负荷增量聚类方法的流程示意框图。
[0110]
图2是本发明实施方案中最大负荷日负荷曲线图。
[0111]
图3是本发明实施方案中用户用电数据集a-fcm聚类的结果图。
[0112]
图4是本发明实施方案中典型用户的直方图和qq图。
[0113]
图5是本发明实施方案中用户用电负荷数据集在无条件矩估计值情况下进行dcs-fcm聚类的结果图。
[0114]
图6是本发明实施方案中用户用电负荷数据集在相同条件均值下进行dcs-a-fcm聚类的结果图。
[0115]
图7是本发明实施方案中用户用电负荷数据集在相同条件方差下进行dcs-a-fcm聚类的结果图。
[0116]
图8是本发明实施方案中用户用电负荷数据集在相同条件均值下聚类的数据权重。
[0117]
图9是本发明实施方案中用户用电负荷数据集在相同条件方差下聚类的数据权重。
[0118]
图10是本发明实施方案中用户用电负荷数据集在相同条件均值dcs-wi-fcm下聚类的结果图。
[0119]
图11是本发明实施方案中用户用电负荷数据集在相同条件方差dcs-wi-fcm下聚类的结果图。
[0120]
图12是本发明实施方案中用户用电负荷数据集聚类性能评价指标值示意图。
[0121]
图13是本发明实施方案中用户用电负荷数据集在两个条件矩上进行dcs-wi-fcm聚类的结果图。
具体实施方式
[0122]
图1所示实施例表明，本发明dcs统计模型与自相关加权增量模糊c均值聚类算法相结合的电力负荷增量聚类方法的流程如下：
[0123]
1.用户日电力负荷数据预处理与分析
→
2.自适应确定聚类数
→
3.确定dcs模型的高斯分布用户假设
→
4.建立dcs模型，并利用极大似然估计法估计模型参数
→
5.计算dcs模型的相关参数，得到最终的dcs模型
→
6.基于dcs模型计算时间序列的自相关值
→
7.根据序列自相关估计值进行irfcm聚类
→
8.根据irfcm聚类得到的数据权重进行基于dcs模型的awi-fcm聚类。
[0124]
实施例1
[0125]
本实施例的dcs统计模型与自相关增量模糊c均值聚类算法相结合的电力负荷增量聚类方法，具体步骤如下：
[0126]
第一步，用户日负荷数据预处理与分析：
[0127]
获取中国k＝11个用电公司的2020年3月最大负荷日的日负荷数据集，对数据集进行数据预处理，包括对缺失值的查找与填补和对异常值的检测与修正，得到的新数据集作为电力负荷时间序列模型建立以及聚类分析所用数据，绘制用户日负荷曲线，通过对日负荷曲线的波动情况进行观察分析，对负荷曲线进行初步分类。
[0128]
第二步，自适应确定聚类数：
[0129]
通过i-nice算法确定聚类数并不依赖任何聚类算法，其通过对原始数据集的数据分布等特征进行分析得到数据集最佳聚类数：
[0130]
第2.1步，设观测点个数为q，观测点集合为p＝{p1，p2，...，pq}，数据集为n行m列的矩阵数据点个数为n，任意观测点到数据集的距离集合为y
p
＝{y1，y2，...yn}，且满足yk≥0，k＝1，2，...n，假设距离分布为混合伽马模型gamma mixture model(gmm)。
[0131]
距离集合y
p
的gmm模型定义为：
[0132][0133]
其中θ为模型的参数向量，θj包含了形状参数αj和尺度参数βj，πj是混合比例，每个分量的密度函数为单伽马分布的密度函数：
[0134][0135]
其中形状参数αj和尺度参数βj应满足αj＞0，βj＞0，混合比例应满足：伽马函数γ(x)的定义为：
[0136][0137]
第2.2步，对于距离集合y
p
的gmm模型，运用em算法通过最大化对数似然函数求解模型参数，其对数似然函数的表达式为：
[0138][0139]
第2.3步，引入随机变量z＝{zi}，zi表示y
p
中的yi被分配给第j个伽马成分，在e步
中根据初始化或上一步参数值计算y
p
中yi属于第j个成分的概率p(zi＝j|yi，θn)为：
[0140][0141]
对数似然函数关于随机变量z＝{zi}的期望值计算公式为式(5)，在m步中通过最大化对数似然函数的期望值求解下一次迭代中参数θ
n+1
的估计值，如式(6)所示，
[0142][0143][0144]
第2.4步，运用拉格朗日乘子法求解期望最大化对数似然函数的参数值，得到混合比例的参数估计值其计算公式为：
[0145][0146]
形状参数和尺度参数的估计值和为：
[0147][0148][0149]
其中为digamma函数，式(9)是非线性方程，没有封闭形式的解，运用牛顿法求解参数的值。
[0150]
第2.5步，每个观测点建立mmax个gmm模型，运用二阶akaik信息准则(aicc)计算指标值进行模型的选择，aicc值达到最大时的gmm成分数即为最佳的聚类数。
[0151]
aicc的计算公式为：
[0152][0153]
其中l(θ
*
)为对数似然函数的最大值，n为数据点个数，q为参数的数量，q＝3m，通过计算aicc的值确定每个观测点处的最佳gmm模型，比较多个观测点的最佳gmm模型得到最终的聚类数。
[0154]
第三步，确定dcs的高斯分布用户假设：
[0155]
随机选取初步分类的每个类别中的用电用户，绘制直方图以及qq图(quantile-quantile plot)进行可视化，由图可知，qq图实际上是一个散点图，若日负荷时间序列数据的散点位于红线附近，则说明观察到的向量近似服从高斯分布，将高斯分布(12)作为用户数据y
t
的基础假设：
[0156][0157]
其中μ
t
为时变参数均值，为时变参数方差；
[0158]
第四步，建立dcs模型，并利用极大似然估计法估计模型参数：
[0159]
建立dcs(p，q)模型，令p＝q＝1，用电用户数据得到基于高斯分布假
设条件下的dcs(1，1)模型，也即f
t
第t个实现的更新方程，如式(13)所示：
[0160]ft
＝ω+as
t-1
+bf
t-1
ꢀꢀꢀꢀꢀꢀ
(13)
[0161]
其中时变参数向量常数向量实矩阵实矩阵其中标量参数包括ω
μ
、aμ等是静态参数向量θ的函数，对ω，a和b运用极大似然估计方法进行参数估计；
[0162]
第五步，计算模型的相关参数，得到最终的dcs模型：
[0163]
根据t-1时刻的驱动向量s
t-1
及时变参数向量f
t-1
计算得到dcs模型中f
t
，其中驱动向量
[0164]
第5.1步，计算条件得分向量其计算公式如式(14)(15)所示：
[0165][0166][0167]
第5.2步，计算缩放矩阵s
t
，其计算公式如式(16)(17)所示：
[0168][0169][0170]
第5.3步，得到时变参数更新方程，如式(18)(19)所示：
[0171][0172][0173]
第5.4步，定义无条件矩为：
[0174]
κ＝(i
n-b)-1
ω
ꢀꢀꢀꢀꢀꢀꢀ
(20)
[0175]
其中in为单位矩阵，将参数b和ω的估计值带入上述公式即得到时变参数均值和方差的无条件矩的估计值。无条件矩表示长期平均值，随着时间的推移，波动将达到一个平均值，通过对无条件矩进行计算能够对时间序列数据进行长期分布的分析。
[0176]
将dcs的参数估计值和相应的参数带入时变参数的更新方程得到时变均值和时变方差的条件矩估计值，根据式(20)得到无条件矩估计值，将条件矩估计值和无条件矩估计值作为聚类算法的输入数据集。
[0177]
至此，完成电力负荷数据的自适应确定聚类数的dcs模型建立。
[0178]
时间序列电力负荷数据基于自适应确定聚类数的dcs模型的awi-fcm聚
[0179]
类分析：
[0180]
第六步，根据dcs模型计算时间序列的自相关值：
[0181]
在动态条件评分模型中，对每个模型的第r(r＝1，2，...，r)个条件矩估计进行r级聚类。dcs(1，1)模型意味着r＝2，因此r＝1和r＝2的时间序列的条件矩估计分别意味着dcs(1，1)模型的参数和的估计。根据第r个条件矩估计来获得第r个估计自相关值将时变均值和时变方差的条件矩估计值作为计算自相关值的数据集分别进行计算，也就是计算时间序列y
t
的第r个条件矩在滞后l时的估计自相关值，计算公式如(21)所示：
[0182][0183]
其中是第k个时间序列从时刻t到时刻t-l过程的r阶条件矩的均值，时间序列k和k
′
基于自相关的r阶条件矩之间的距离为：
[0184][0185]
第七步，根据所得的序列自相关估计值数据集进行irfcm聚类：
[0186]
基于rfcm聚类算法进行改进得到加权模糊c均值聚类算法(irfcm)，其同时考虑数据点和数据特征的不同权重，通过加入可调参数α避免数据集的数据个数对权重信息和聚类中心产生影响，具体算法如下：
[0187]
第7.1步，设数据集为k行t列的矩阵数据点个数为k，数据特征个数为t，聚类算法的目标函数和约束条件为：
[0188][0189][0190][0191]
其中，聚类数为c，vi为第i个聚类中心，u
ik
为第k个数据点xk对第i个聚类中心vi的隶属度，r
ij
为第j个数据点对第i个聚类类别的特征权重，w
ik
为第k个数据点xk对第i个聚类类别的数据权重。
[0192]
第7.2步，利用数据特征权重和数据权重的l2范数正则化协调权重分散。
[0193]
其中，η是可调参数，δ和ε是正则化参数，通过调整δ的值使得更多数据特征信息参与聚类中心的形成，使得聚类算法聚类准确性更高，α参数用以避免大数据集的数据个数对权重信息产生影响，当数据集数据个数k≤40时，数据个数对权重信息影响较小，当数据个数k＞40时，数据个数将影响产生权重矩阵的精确度，权重矩阵的分辨率较低，无法进行有效聚类，通过调整α的值提高权重矩阵的分辨率和精确度。
[0194]
第7.3步，利用拉格朗日乘子法求解目标函数和约束条件，得到关于模糊隶属度u
ik
、数据特征权重r
ij
、数据权重w
ik
、聚类中心v
ij
的迭代公式：
[0195][0196][0197][0198][0199]
第八步，根据irfcm聚类得到的数据权重进行基于dcs模型的awi-fcm：
[0200]
将a-wfcm算法作为基础，根据所得自相关估计值数据集进行irfcm聚类，得到数据权重w
ik
，将权重信息w
ik
和数据集作为增量聚类的算法基础，具体算法如下：
[0201]
第8.1步，将k个时间序列数据点随机分为p(p＝1，2，...，p)个数据块，每个数据块有k/p个数据点，对第一个数据块p＝1执行a-wfcm聚类。
[0202]
第8.1.1步，根据i-nice算法得到的聚类数c随机初始化隶属度矩阵，使得其满足以下约束条件：
[0203][0204]
第8.1.2步，根据得到的隶属度矩阵计算聚类中心，其计算公式为：
[0205][0206]
第8.1.3步，根据隶属度矩阵和聚类中心计算目标函数j的值，其计算公式如式(31)所示，若前后两次j值之差小于规定的阈值ε，则结束迭代，输出隶属度矩阵和聚类中心，否则进行下一步；
[0207][0208]
第8.1.4步，根据聚类中心重新计算得到隶属度矩阵并返回步骤8.1.2，隶属度矩阵的计算公式为：
[0209][0210]
第8.2步，重新根据式(33)赋予计算得到的数据块聚类中心权重，将该数据块聚类中心和计算得到的权重加入下一个数据块中进行a-wfcm聚类(步骤8.1.1至8.1.4)；
[0211][0212]
其中，n为当前数据块中数据点的个数，j为上一个数据块中聚类中心的个数。
[0213]
第8.3步，重复步骤8.2直到所有数据块完成聚类，最后根据最终的聚类中心结果根据式(32)重新计算所有数据点的隶属度。
[0214]
至此，完成电力负荷数据的基于dcs模型的自相关加权增量模糊c均值聚类。
[0215]
图2显示了11家用电公司最大负荷日的负荷曲线图，根据图中各公司的波动情况进行简单分类，具有明显波峰的公司分为一类，波动较平稳的分为一类。
[0216]
图3显示了a-fcm聚类方法对公司数据集进行聚类的得到的聚类结果，11家公司共分成了4类，第一类包括begj12；第二类包括begj3、hgxx、wlym、hpjwsz399、hpjwsz3、szpls和cgc；第三类包括zxdzkj与hzdmt；第四类包括wn公司。
[0217]
图4显示了随机选取4家公司，对其最大负荷日负荷数据进行直方图和qq图的绘制，图像显示最大日负荷时间序列数据的散点位于红线附近，表明观察到的向量近似服从高斯分布，进而将高斯分布作为用户假设。
[0218]
图5显示了11家公司的最大日负荷时间序列的基于dcs模型无条件矩下的聚类结果图，共有4个类别，第一个类别包括begj12、hzdmt、zxdzkj和wn；第二个类别有begj3、hgxx、wlym；公司hpjwsz399、hpjwsz3、szpls属于第三个类别，属于cgc第四个类别。其中公司hgxx有31.0％的概率属于第一个类别，wlym有9.5％的概率属于第一个类别，隶属度的不确定性较高。
[0219]
图6和图7显示了运用dcs-a-fcm聚类方法对公司数据集进行聚类得到的聚类结果，根据条件均值进行聚类的结果是：第一个类别包括begj12和hzdmt；第二个类别有
begj3、hgxx、wlym、hpjwsz399、hpjwsz3、szpls和cgc；wn公司属于第三个类别，zxdzkj属于第四个类别。在相同的条件方差下进行的聚类结果时：第一个类别有begj12、hzdmt、zxdzkj；第二类包括begj3、hgxx、wlym、hpjwsz3和cgc；第三类为wn公司，第四类为hpjwsz399、szpls公司。总的来说，根据条件方差进行的dcs-a-fcm聚类是不确定性水平较高的聚类，因为所有公司对所属类别的隶属度都较小，但基于条件均值的聚类具有更低的不确定性。
[0220]
图8和图9是用户用电负荷数据集在相同条件均值下和在相同条件方差下聚类的数据权重。根据相同条件均值下数据权重信息图可以看出，第一个类别的形成中数据权重较大的有wlym公司，第二个类别中begj3公司对该类形成的影响最大，第三类只有一家公司，第四类中begj12公司对该类形成的影响最大；根据相同条件方差下数据权重信息图可以看出，第一个类别的形成中数据权重较大的有begj12公司，其对第一类别的形成影响最大，第二个类别中begj3公司对该类形成的影响最大，第三类和第四类只有一家公司。根据数据特征权重信息图可以看出，在相同的条件方差下进行聚类时，在时间段00：00
‑‑
16：35的时变参数条件方差对聚类类别中类别二的形成影响较大，其他类别中特征权重均在0.0035上下波动且波动幅度较小。
[0221]
图10和图11显示了运用dcs-wi-fcm聚类方法对公司数据集进行聚类得到的聚类结果，在条件平均聚类中，第一个聚类包含begj12和hzdmt，zxdzkj在第二个聚类中，wn在第三个聚类中；第四个聚类包括begj3、hgxx、wlym、hpjwsz399、hpjwsz3、szpls、cgc；在相同的条件方差下，第一个聚类包括begj12、hzdmt、zxdzkj；第二集群为begj3、hpjwsz3、cgc，wn公司位于第三集群，第四集群为hgxx、wlym、hpjwsz399、szpls。
[0222]
图12显示了针对11家用电公司的最大负荷日数据进行的a-fcm聚类、paa-fcm聚类、dcs-a-fcm(mean)聚类、dcs-a-fcm(var)聚类、dcs-wi-fcm(mean)聚类和dcs-wi-fcm(var)聚类结果进行对比，分别采用sc、chi、dbi、di内部评价指标并绘制折线图进行有效性评价，四个评价指标值表明，基于条件均值的dcs-a-fcm聚类和dcs-wi-fcm聚类结果是相同的。当根据条件方差进行聚类时，通过比较四个评价指标的值，dcs-wi-fcm聚类比dcs-a-fcm聚类好得多。总之，dcs-wi-fcm算法将自适应确定聚类数的dcs统计模型与自相关加权增量模糊聚类相结合，比a-fcm聚类算法和经典的时间序列聚类方法paa-fcm更好、更有效。
[0223]
图13是用户用电负荷数据集在两个条件矩上进行dcs-wi-fcm聚类的结果图。将11家公司中的cgc公司来作为增量验证dcs-wi-fcm聚类算法的有效性，cgc公司按照条件均值进行聚类时属于类别4的概率为99.7％，基于条件方差进行聚类时，cgc属于第二类中的概率为71.8％，属于第4类的概率为17.8％。根据两个条件矩聚类，begj3、hpjwsz3和cgc公司始终处于同一组中，这意味着cgc公司与begj3和hpjwsz3公司共享相同的条件分布。

技术特征：
1.自适应确定聚类数的dcs统计模型与自相关加权增量模糊c均值聚类算法相结合的电力负荷增量聚类方法，其特征在于：针对当前电力负荷时间序列数据存在的高维性、波动性和相关性等特性，并且传统聚类算法存在初始化聚类数目难以确定，对噪声点敏感等问题，为了描述时间序列数据的前后的依赖关系，降低高维时间序列数据下聚类的计算成本，结合不同用电负荷数据流时序特性，构建基于自适应确定聚类数的dcs统计模型的自相关加权增量模糊c均值聚类算法具体步骤如下：第一步，用电用户日负荷数据预处理与分析：获取k个用户的同年同月最大负荷日的日负荷数据集，对数据集进行数据预处理，包括对缺失值的查找与填补和对异常值的检测与修正，得到的新数据集作为电力负荷时间序列模型建立以及聚类分析所用数据，绘制用户日负荷曲线，通过对日负荷曲线的波动情况进行观察分析，对负荷曲线进行初步分类。第二步，自适应确定聚类数：通过i-nice算法确定聚类数并不依赖任何聚类算法，其通过对原始数据集的数据分布等特征进行分析得到数据集最佳聚类数：第2.1步，设观测点个数为q，观测点集合为p＝{p1，p2，...，p
q
}，数据集为n行m列的矩阵数据点个数为n，任意观测点到数据集的距离集合为y
p
＝{y1，y2，...y
n
}，且满足y
k
≥0，k＝1，2，...n，假设距离分布为混合伽马模型gamma mixture model(gmm)。距离集合y
p
的gmm模型定义为：其中θ为模型的参数向量，θ
j
包含了形状参数α
j
和尺度参数β
j
，π
j
是混合比例，每个分量的密度函数为单伽马分布的密度函数：其中形状参数α
j
和尺度参数β
j
应满足α
j
＞0，β
j
＞0，混合比例应满足：伽马函数γ(x)的定义为：第2.2步，对于距离集合y
p
的gmm模型，运用em算法通过最大化对数似然函数求解模型参数，其对数似然函数的表达式为：第2.3步，引入随机变量z＝{z
i
}，z
i
表示y
p
中的y
i
被分配给第j个伽马成分，在e步中根据初始化或上一步参数值计算y
p
中y
i
属于第j个成分的概率p(z
i
＝j|y
i
，θ
n
)为：对数似然函数关于随机变量z＝{z
i
}的期望值计算公式为式(5)，在m步中通过最大化对数似然函数的期望值求解下一次迭代中参数θ
n+1
的估计值，如式(6)所示，
第2.4步，运用拉格朗日乘子法求解期望最大化对数似然函数的参数值，得到混合比例的参数估计值其计算公式为：形状参数和尺度参数的估计值和为：为：其中为digamma函数，式(9)是非线性方程，没有封闭形式的解，运用牛顿法求解参数的值。第2.5步，每个观测点建立mmax个gmm模型，运用二阶akaik信息准则(aicc)计算指标值进行模型的选择，aicc值达到最大时的gmm成分数即为最佳的聚类数。aicc的计算公式为：其中l(θ
*
)为对数似然函数的最大值，n为数据点个数，q为参数的数量，q＝3m，通过计算aicc的值确定每个观测点处的最佳gmm模型，比较多个观测点的最佳gmm模型得到最终的聚类数。第三步，确定dcs的高斯分布用户假设：随机选取初步分类的每个类别中的用电用户，绘制直方图以及qq图(quantile-quantile plot)进行可视化，由图可知，qq图实际上是一个散点图，能够显示数据集与所选理论分布的拟合程度，若日负荷时间序列数据的散点位于红线附近，则说明观察到的向量近似服从高斯分布，将高斯分布(12)作为用户数据y
t
的基础假设：其中μ
t
为时变参数均值，为时变参数方差；第四步，建立dcs模型，并利用极大似然估计法估计模型参数：建立dcs(p，q)模型，令p＝q＝1，用电用户数据得到基于高斯分布假设条件下的dcs(1，1)模型，也即f
t
第t个实现的更新方程，如式(13)所示：f
t
＝ω+as
t-1
+bf
t-1
ꢀꢀꢀꢀꢀꢀꢀ
(13)其中时变参数向量常数向量实矩阵实矩阵其中标量参数包括ω
μ
、α
μ
等是静态参数向量θ的函数，对ω，a和b运用极大似然估计方
法进行参数估计；第五步，计算模型的相关参数，得到最终的dcs模型：根据t-1时刻的驱动向量s
t-1
及时变参数向量f
t-1
计算得到dcs模型中f
t
，其中驱动向量第5.1步，计算条件得分向量其计算公式如式(14)(15)所示：其计算公式如式(14)(15)所示：第5.2步，计算缩放矩阵s
t
，其计算公式如式(16)(17)所示：，其计算公式如式(16)(17)所示：第5.3步，得到时变参数更新方程，如式(18)(19)所示：第5.3步，得到时变参数更新方程，如式(18)(19)所示：第5.4步，定义无条件矩为：κ＝(i
n-b)-1
ω
ꢀꢀꢀꢀꢀꢀꢀ
(20)其中i
n
为单位矩阵，将参数b和ω的估计值带入上述公式即得到时变参数均值和方差的无条件矩的估计值。无条件矩表示长期平均值，随着时间的推移，波动将达到一个平均值，通过对无条件矩进行计算能够对时间序列数据进行长期分布的分析。将dcs的参数估计值和相应的参数带入时变参数的更新方程得到时变均值和时变方差的条件矩估计值，根据式(20)得到无条件矩估计值，将条件矩估计值和无条件矩估计值作为聚类算法的输入数据集。至此，完成电力负荷数据的自适应确定聚类数的dcs模型建立。时间序列电力负荷数据基于自适应确定聚类数的dcs模型的awi-fcm聚类分析：第六步，根据dcs模型计算时间序列的自相关值：在动态条件评分模型中，对每个模型的第r(r＝1，2，...，r)个条件矩估计进行r级聚类。dcs(1，1)模型意味着r＝2，因此r＝1和r＝2的时间序列的条件矩估计分别意味着dcs(1，1)模型的参数和的估计。根据第r个条件矩估计来获得第r个估计自相关值将时变均值和时变方差的条件矩估计值作为计算自相关值的数据集分别进行计算，也就是计算时间序列y
t
的第r个条件矩在滞后l时的估计自相关值，计算公式如(21)所示：其中是第k个时间序列从时刻t到时刻t-l过程的r阶条件矩的均值，时间序列k和k
′
基于自相关的r阶条件矩之间的距离为：
第七步，根据所得的序列自相关估计值数据集进行irfcm聚类：基于rfcm聚类算法进行改进得到加权模糊c均值聚类算法(irfcm)，其同时考虑数据点和数据特征的不同权重，通过加入可调参数α避免数据集的数据个数对权重信息和聚类中心产生影响，具体算法如下：第7.1步，设数据集为k行t列的矩阵数据点个数为k，数据特征个数为t，聚类算法的目标函数和约束条件为：类算法的目标函数和约束条件为：类算法的目标函数和约束条件为：其中，聚类数为c，v
i
为第i个聚类中心，u
ik
为第k个数据点x
k
对第i个聚类中心v
i
的隶属度，r
ij
为第j个数据点对第i个聚类类别的特征权重，w
ik
为第k个数据点x
k
对第i个聚类类别的数据权重。第7.2步，利用数据特征权重和数据权重的l2范数正则化协调权重分散。其中，η是可调参数，δ和ε是正则化参数，通过调整δ的值使得更多数据特征信息参与聚类中心的形成，使得聚类算法聚类准确性更高，α参数用以避免大数据集的数据个数对权重信息产生影响，当数据集数据个数k≤40时，数据个数对权重信息影响较小，当数据个数k＞40时，数据个数将影响产生权重矩阵的精确度，权重矩阵的分辨率较低，无法进行有效聚类，通过调整α的值提高权重矩阵的分辨率和精确度。第7.3步，利用拉格朗日乘子法求解目标函数和约束条件，得到关于模糊隶属度u
ik
、数据特征权重r
ij
、数据权重w
ik
、聚类中心v
ij
的迭代公式：的迭代公式：的迭代公式：的迭代公式：第八步，根据irfcm聚类得到的数据权重进行基于dcs模型的awi-fcm：将a-wfcm算法作为基础，根据所得自相关估计值数据集进行irfcm聚类，得到数据权重w
ik
，将权重信息w
ik
和数据集作为增量聚类的算法基础，，具体算法如下：第8.1步，将k个时间序列数据点随机分为p(p＝1，2，...，p)个数据块，每个数据块有k/p个数据点，对第一个数据块p＝1执行a-wfcm聚类。第8.1.1步，根据i-nice算法得到的聚类数c随机初始化隶属度矩阵，使得其满足以下约束条件：
第8.1.2步，根据得到的隶属度矩阵计算聚类中心，其计算公式为：第8.1.3步，根据隶属度矩阵和聚类中心计算目标函数j的值，其计算公式如式(31)所示，若前后两次j值之差小于规定的阈值ε，则结束迭代，输出隶属度矩阵和聚类中心，否则进行下一步；第8.1.4步，根据聚类中心重新计算得到隶属度矩阵并返回步骤8.1.2，隶属度矩阵的计算公式为：第8.2步，重新根据式(33)赋予计算得到的数据块聚类中心权重，将该数据块聚类中心和计算得到的权重加入下一个数据块中进行a-wfcm聚类(步骤8.1.1至8.1.4)；其中，n为当前数据块中数据点的个数，j为上一个数据块中聚类中心的个数。第8.3步，重复步骤8.2直到所有数据块完成聚类，最后根据最终的聚类中心结果根据式(32)重新计算所有数据点的隶属度。至此，完成电力负荷数据的基于自适应确定聚类数的dcs模型的自相关加权增量模糊c均值聚类。2.根据权利要求1所述自适应确定聚类数的dcs统计模型与自相关加权增量模糊c均值聚类算法相结合的电力负荷增量聚类方法，其特征在于：所述第二步中通过i-nice算法自适应确定聚类数，该算法并不依赖任何聚类算法，其通过对原始数据集的数据分布等特征进行分析得到数据集最佳聚类数。假设数据点到数据集的距离分布为混合伽马模型，通过em算法计算其中参数，运用二阶akaik信息准则(aicc)计算指标值，当aicc值达到最大时的gmm成分数即为最佳的聚类数。有效解决了传统聚类算法存在的难以设定初始化聚类数目的问题。3.根据权利要求1所述自适应确定聚类数的dcs统计模型与自相关加权增量模糊c均值聚类算法相结合的电力负荷增量聚类方法，其特征在于：所述第三步中时间序列用户数据y
t
中t的值取决于测量电力负荷数据的分辨率，时变参数是根据每个时间点t以及原始时间序列y
t
进行计算得到，绘制初步分类中典型用户的直方图和qq图来确定dcs模型的用户假设，将高斯分布作为全部用户的dcs模型的概率分布假设，用来描述多种情况下的随机变量，具有较高的适普性。4.根据权利要求1所述自适应确定聚类数的dcs统计模型与自相关加权增量模糊c均值聚类算法相结合的电力负荷增量聚类方法，其特征在于：所述第四步中dcs(p，q)模型令p＝q＝1已可以在模型复杂度较低的情况下描述时间序列数据，根据高斯分布用户假设的参数确定时变参数向量f
t
的向量表达式，同时确定后续待估计参数的向量表达式。5.根据权利要求1所述自适应确定聚类数的dcs统计模型与自相关加权增量模糊c均值聚类算法相结合的电力负荷增量聚类方法，其特征在于：所述第七步中利用加权模糊c均值
聚类算法(irfcm)，其考虑数据点的不同权重，并加入可调参数α避免数据集的数据个数对权重信息和聚类中心产生影响，利用拉格朗日乘子法得到权重信息。

技术总结
本发明结合了自适应确定聚类数的DCS统计模型与加权增量模糊C均值聚类算法，是针对电力负荷数据存在的高维性和异方差性、传统聚类算法初始化聚类数目难以确定、对噪声点敏感等问题进行的时间序列聚类。首先通过I-nice算法自适应确定最佳聚类数，构建动态条件评分模型，并计算参数序列的自相关值；然后进行自适应加权模糊C均值聚类分析得到数据权重信息；最后结合DCS模型参数数据集和数据权重信息进行增量聚类。本发明能够自适应确定电力负荷数据集上的最佳聚类数，并通过对数据点分配适当权重来减少噪声点对聚类中心的影响，有效解决数据高维性和异方差性导致的聚类效果差的问题，重点解决随时间推移而出现新的隐藏信息的有效聚类的问题。有效聚类的问题。有效聚类的问题。


技术研发人员：张勇 李欣玥 王莉
受保护的技术使用者：天津商业大学
技术研发日：2023.04.19
技术公布日：2023/7/21