[发明专利]基于互近似熵与聚类的多谐波源谐波责任划分方法

权利要求书

1.一种基于互近似熵与聚类的多谐波源谐波责任划分方法，其特征在于，包括步骤：

(1)通过判断互近似熵值的大小，筛选出有效的谐波电压、电流数据，形成新的数据集A；

(2)调用DBSCAN聚类算法对数据集A进行聚类分析，设置参数ε与MinPts的值，针对不同的系统谐波阻抗划分为不同的类；

(3)根据每个类簇的谐波数据，取出相对应时刻的馈线电流数据，利用最小二乘多元线性回归方法对谐波阻抗进行计算，得到谐波责任划分结果。

2.根据权利要求1所述的基于互近似熵与聚类的多谐波源谐波责任划分方法，其特征在于，所述步骤1具体包括：

(1.1)采集电压、电流谐波数据，进行归一化处理，并将其分为n个数据区段，每个区段包含L个数据点；

(1.2)计算每个数据区段内的互近似熵值CAE，将其与设定阈值比较，筛选出有效的谐波电压、电流数据，成新的数据集A。

3.根据权利要求2所述的基于互近似熵与聚类的多谐波源谐波责任划分方法，其特征在于，所述步骤1.2具体包括：两个不同的时间序列分别为i(t)和j(t)，对i(t)和j(t)构造N-m+1个m维矢量X_p、X_q：

X_p＝[i(p),…,i(p+m-1)],p＝1,2,…,N-m+1 (2)

X_q＝[j(q),…,j(q+m-1)],q＝1,2,…,N-m+1 (3)

使用矢量的∞—范数描述X_p与X_q之间的距离d：

d(X_p,X_q)＝||X_p-X_q||_∞

根据X_p、X_q的协方差COV可计算出相似容限r：

r＝0.2×COV(X_p,X_q) (4)

统计X_p和X_q的矢量距离小于r的个数N_p,m,r，并计算N_p,m,r与总矢量个数(N-m+1)的比值C_p,m,r：

再求两条曲线之间的互相关程度：

增加窗口长度至m+1，重复上述公式(2)～(6)的运算过程，得到Φ_m+1,r，并计算得到与m、r相关的互近似熵值：

CAE(m,r)＝Φ_m,r-Φ_m+1,r (7)。