基于量子遗传优化的非侵入式电器负荷识别方法与流程

本发明涉及非侵入式电力负荷识别技术领域，特别涉及一种基于量子遗传优化的非侵入式电器负荷识别方法。

背景技术：

随着智能电网的发展，电力负荷中居民用户用电负荷所占的比重越来越大，作为电力负荷的重要组成部分，居民用户域负荷日益引起社会的广泛关注。居民用户域的负荷在线监测是实现居民用户用电可视化的基础，它有助于用户了解家庭内不同时段各电器设备的具体能耗情况，据此来制定合理的用电计划，完善能源消费结构，促进能源有效利用，减少家庭电费开支。居民用户域负荷在线监测对推进整个社会的节能减排、缓解能源危机起着重要的作用。

传统的侵入式电力负荷监测系统需要在各个负载前加装传感器，综合成本高，施工非常复杂，对原有线路也造成破坏。而非侵入式电力负荷监测nilm系统仅需在住户进户线总开关处或者工业车间总线上安装一个传感器就能获知每个/类电器的使用情况，并且成本低，施工方便，不对用户造成用电干扰，有望发展成为新一代智能电表的核心技术，为电力用户和整个社会带来多方面的效益。

负荷辨识算法是nilm的核心内容，研究人员提出了很多不同的理论和算法。george.w.hart教授首先提出了一种通过监测电力人口处所有用电设备总的有功功率和无功功率的变化量来辨识负荷的方法。robertson等人采用小波变换技术，通过对几个未知的暂态信息进行分类的方法进行负荷辨识。由于负荷稳态特征存在局限性，利用暂态特征辨识成为研究重点。

上述方法均是通过对暂态特征量的提取、转化来实现负荷辨识，但由于切入点比较单一，考虑到实际环境下暂态特征的不确定性其准确度将受到影响，暂态、稳态特征结合的方法也受到了关注。且随着智能算法的发展，有研究者利用智能算法进行负荷辨识的研究。

技术实现要素：

本发明的目的是基于上述背景技术，提供一种基于量子遗传优化的非侵入式电器负荷识别方法，利用实测负荷电流及电压数据，获得电流有效值，通过量子遗传算法寻优，寻找最优解，最终确定电器的具体负荷类型。

为了达到上述的技术效果，本发明采取以下技术方案：

基于量子遗传优化的非侵入式电器负荷识别方法，包括以下步骤：

a.采集电器设备电流及电压数据，并通过对采集的电流及电压数据进行预处理得到有效的电流数据；

具体的，在采集电器设备电流及电压数据前会预设一个采集周期n，则将对采集周期n时间段内的电器设备的电流及电压数据进行采集，一般具体的采集周期n的时长可根据具体需求而设定，如可设定为20秒或30秒等；

通过对采集的电流及电压数据进行预处理得到有效的电流数据即为该类电器设备的单个样本数据，其中，优选同一种类电器设备需采集多份电流及电压数据并求取其平均值，从而减小误差；

同时，步骤a中对采集的电流及电压数据进行预处理即是对采集到的直接数据进行计算处理得到需要的有效的数据，比如电流有效值，功率等需求的数据，不同的数据类型其计算方式不同，且各计算处理方式均为现有技术中已有的计算方法，此处不再赘述。

其中，在得到每种电器设备的有效电流数据后再通过电压相位对齐找到m域的电流起点，根据电流起点计算每类电器m域内的电流幅值的最大值，最小值等，具体的，m域就是采样周期中的数据波形进行相位对齐获取的一个波段，且该波段内获取的数据个数具体为m个。

b.配置量子遗传算法参数；

此步骤中配置的量子遗传算法参数至少包括种群大小、迭代次数nmax、每类电器数量作为定义域的区间值、染色体长度；

其中，种群大小参数、迭代次数nmax的具体数值可根据具体需求设定；

每类电器数量作为定义域的区间值具体通过居民家庭中各类电器的数量上限决定，如居民家庭的冰箱数量上限设置为2，则冰箱的区间值即为[0,2]；

染色体长度具体是根据居民用电设备种类数来定义的，如可以定义为10或20等。

c.采用量子比特编码来初始化染色体矩阵；

步骤中采用量子比特编码进行初始化染色体矩阵时，编码初始化值设置为被初始化后的染色体矩阵具体如下：为2×n矩阵，其中，n为定义的家用电器种类数即染色体长度，此步骤c初始化是随机初始化【0,1】区间内的随机值，满足归一化原则即可，本发明中初始化为

d.计算解空间转换得到当前迭代中每类电器的数量；

具体是通过将步骤b中得到的每类电器数量作为定义域的区间值作为解空间转换的变量，进行解空间转换从而得到当前迭代中每类电器的数量，每次迭代通过解转换得到的当前迭代获得的每类电器的数量解具体的解空间转换公式如下：

其中为当前迭代解转换得到的第i类电器数量解，[mj,nj]表示第j种电器数量的区间值，表示第i条染色体的第j个量子位，根据步骤c设定的编码初始化值则具体的，i的取值范围为不小于1且不大于n，j的取值范围为不小于1且不大于m。

e.通过量子旋转门进行染色体更新及变异，计算适应度、记录最优解；

此步骤中将首先进行量子旋转门的选取以及染色体更新及变异，其中，量子旋转门选取为：其中θ为旋转角，θ取值范围为[0,2π]，对染色体进行更新变异时具体如下：

其中，为染色体中第i个基因，为通过量子旋转门更新后的基因，r(θ)为量子旋转门；

然后再进行实测数据的获取，得到各类电器设备的实测数据的电流幅值的最大值，最小值，以及对应m域内的m个点的电流数据列表，从而得到实测数据；

然后再根据步骤b中设定的迭代次数nmax，启动量子遗传算法，通过量子遗传编码，解空间转换得到的设备数量列表，分别计算适应度函数值，并获取最优值(即适应度函数求解的最小值)作为下一代的更新及变异目标，然后按照上述的步骤进行旋转门转换及染色体变异更新染色体；

其中，计算适应度函数值时所选择的适应度函数具体为：

其中qi表示实测数据中m域内的m个点的电流数据列表中第i个电流数据，pij表示第j种电器的单个样本数据的对应的m域内的m个点的电流值中的第i个电流值，nij表示第j种电器数量，表示第j种电器的权重，其中，pkj表示第j种电器的单个样本数据的对应的m域内的m个点的电流值中的第k个电流值。

f.判断是否迭代完成，若完成则进入步骤g，否则返回步骤e；

判断迭代次数是否满足迭代次数nmax，若满足则进入下一步，否则返回步骤e继续进行迭代；

g.获取最优解；

完成迭代后将得到nmax个最优个体列表，求取最小值即为最后的最优解，即最后识别的电器种类和数量。

本发明与现有技术相比，具有以下的有益效果：

本发明的基于量子遗传优化的非侵入式电器负荷识别方法，是针对之前遗传优化的负荷识别算法进行分析，将量子遗传算法应用到负荷识别系统中，与遗传算法进行对比，从识别结果上有实际性的优化，可实现全局寻优，提高计算效率，较传统遗传算法上提升了一倍的解空间，更大提升了多种设备同时运行的识别精准率，也降低了迭代次数和时间复杂度，且在设备数量比较大时，量子遗传相对于遗传算法的优化更明显。

附图说明

图1是本发明的基于量子遗传优化的非侵入式电器负荷识别方法的示意图。

具体实施方式

下面结合本发明的实施例对本发明作进一步的阐述和说明。

实施例：

如图1所示，一种基于量子遗传优化的非侵入式电器负荷识别方法，具体包括以下步骤流程：

步骤一：安装非侵入式电力监测系统装置到需要监测和负荷识别的居民用户中；

步骤二：预设一个采集周期n，对采集周期n时间段内的电器设备的电流及电压数据进行采集，一般具体的采集周期n的时长可根据具体需求而设定，如本实施例中设定为20秒。

步骤三：对采集的电流及电压数据进行预处理即是对采集到的直接数据进行计算处理得到需要的有效的数据，比如电流有效值，功率等需求的数据，并在得到每种电器设备的有效电流数据后再通过电压相位对齐找到m域的电流起点，根据电流起点计算每类电器m域内的电流幅值的最大值，最小值等。

其中，m域就是采样周期中的数据波形进行相位对齐获取的一个波段，且该波段内获取的数据个数具体为m个。

步骤四：配置量子遗传算法参数；

量子遗传算法参数至少包括种群大小、迭代次数nmax、每类电器数量作为定义域的区间值、染色体长度；

其中，种群大小参数、迭代次数nmax的具体数值可根据具体需求设定，如可设定为300、400或500等，染色体长度具体是根据居民用电设备种类数来定义的，如可以定义为10或20等。

每类电器数量作为定义域的区间值具体通过居民家庭中各类电器的数量上限决定，比如居民家用电器有冰箱，微波炉，油烟机，空调，pc/tv,风扇等，居民家庭的冰箱数量上限设置为2，冰箱数量区间值即为[0,2],同样，微波炉数量区间值设置为[0,1]，油烟机数量区间值设置为[0,1]，空调数量区间值设置为[0,4]，风扇数量区间值设置为[0,4]，pc/tv数量区间值设置为[0,3]等；

步骤五：采用量子比特编码来初始化染色体矩阵；

具体的，采用量子比特编码进行初始化染色体矩阵时，编码初始化值设置为被初始化后的染色体矩阵具体如下：为2×n矩阵，其中，n为定义的家用电器种类数即染色体长度，此步骤的初始化是随机初始化【0,1】区间内的随机值，满足归一化原则即可，本发明中初始化值为

步骤六：计算解空间转换得到当前迭代中每类电器的数量；

具体是通过将步骤四中得到的每类电器数量作为定义域的区间值作为解空间转换的变量即解空间转换的定义域，进行解空间转换从而得到当前迭代中每类电器的数量，每次迭代通过解转换得到的当前迭代获得的每类电器的数量解，具体的解空间转换公式如下：

其中为当前迭代解转换得到的第i类电器数量解，[mj,nj]表示第j种电器数量的区间值，表示第i条染色体的第j个量子位，根据步骤五设定的编码初始化值，则具体的，i的取值范围为不小于1且不大于n，j的取值范围为不小于1且不大于m。

步骤七：通过量子旋转门进行染色体更新及变异，计算适应度、记录最优解；

此步骤中将首先进行量子旋转门的选取以及染色体更新及变异，其中，量子旋转门选取为：其中θ为旋转角，θ取值范围为[0,2π]，对染色体进行更新变异时具体如下：

其中，为染色体中第i个基因，为通过量子旋转门更新后的基因，r(θ)为量子旋转门；

步骤八：获取实测数据：

进行实测数据的获取，得到各类电器设备的实测数据的电流幅值的最大值，最小值，以及对应m域内的m个点的电流数据列表，从而得到实测数据；

然后再根据步骤四中设定的迭代次数nmax，启动量子遗传算法，通过量子遗传编码，解空间转换得到的设备数量列表，分别计算适应度函数值，并获取最优值(即适应度函数求解的最小值)作为下一代的更新及变异目标，然后按照上述的步骤进行旋转门转换及染色体变异更新染色体，重复上面得到新一代的最优值；

其中，计算适应度函数值时所选择的适应度函数具体为：

其中，qi表示实测数据中m域内的m个点的电流数据列表中第i个电流数据，pij表示第j种电器的单个样本数据的对应的m域内的m个点的电流值中的第i个电流值，nij表示第j种电器的电器数量，表示第j种电器的权重，其中，pkj表示第j种电器的单个样本数据的对应的m域内的m个点的电流值中的第k个电流值。

h.判断是否迭代完成，若完成则进入步骤g，否则返回步骤e；

判断迭代次数是否满足迭代次数nmax，若满足则进入下一步，否则返回步骤e继续进行迭代；

i.获取最优解；

完成迭代后将得到nmax个最优个体列表，求取最小值即为最后的最优解，即最后识别的电器种类和数量

可以理解的是，以上实施方式仅仅是为了说明本发明的原理而采用的示例性实施方式，然而本发明并不局限于此。对于本领域内的普通技术人员而言，在不脱离本发明的精神和实质的情况下，可以做出各种变型和改进，这些变型和改进也视为本发明的保护范围。