CN105868790A - 用电负载类型识别器 - Google Patents

Abstract

一种用电负载类型识别器，包括信息采集模块、信息处理模块、通信模块。所述装置同时采用电器的启动电流特征、电器的基波电压电流相位差和负载电流频谱特征作为用电负载类型的识别特征，特征信息丰富；采用包括BP神经网络分类器和贝叶斯分类器的组合分类器进行识别分类，兼顾BP神经网络分类器和贝叶斯分类器的特点进行综合识别，识别准确率高；提供的基波电压电流相位差和启动电流特征、负载电流频谱特征获取方法简单、可靠。所述用电负载类型识别器可以用在学生集体宿舍、办公场所、集体市场等一些需要进行用电电器管理的集体公共场所，也可以用于需要进行用电负载类型识别与统计的其他需要进行用电设备管理的场合。

Description

用电负载类型识别器

技术领域

本发明涉及一种设备识别及分类装置，尤其是涉及一种用电负载类型识别器。

背景技术

目前，主流的用电负载性质识别方法包括基于负载功率综合系数算法的用电负载识别方法、基于电磁感应的用电负载识别方法、基于神经网络算法的用电负载识别方法、基于周期性离散变换算法的用电负载识别方法等。各种方法均能够在一定程度是实现用电负载性质的识别，但由于特征性质单一，识别手段单一，普遍存在泛化能力不够及不能完全准确识别的问题。

发明内容

本发明的目的在于，针对现在已有技术的缺陷，提供一种能够实现高效识别的用电负载类型识别器。所述用电负载类型识别器包括信息采集模块、信息处理模块、通信模块。

所述信息采集模块用于采集用电负载的负载电流并转换成电流数字信号；所述电流数字信号被送至信息处理模块；所述信息处理模块依据输入的电流数字信号，采用组合分类器进行用电负载类型识别；所述通信模块用于发送信息处理模块的用电负载类型识别结果至上位机。

所述组合分类器的输入特征包括用电负载的启动电流特征、用电负载的负载电流频谱特征和用电负载的基波电压电流相位差；所述组合分类器包括BP神经网络分类器和贝叶斯分类器；所述启动电流特征包括启动冲激电流、启动平均电流、启动电流冲量。

所述信息采集模块包括电流传感器、前置放大器、滤波器、A/D转换器；所述信息处理模块的核心为DSP，或者为ARM，或者为单片机，或者为FPGA。

所述A/D转换器可以采用信息处理模块的核心中包括的A/D转换器。

所述信息采集模块、信息处理模块、通信模块的全部或者部分功能集成在一片SoC上。

所述通信模块还接收上位机的相关工作指令；所述通信模块与上位机之间的通信方式包括无线通信方式与有线通信方式；所述无线通信方式包括ZigBee、蓝牙、WiFi、433MHz数传方式；所述有线通信方式包括485总线、CAN总线、互联网、电力载波方式。

所述负载电流频谱特征通过以下方法获得：

步骤一、获取用电负载的稳态电流信号，并将其转换为对应的稳态电流数字信号；

步骤二、对稳态电流数字信号进行傅立叶变换，得到负载电流频谱特性；

步骤三、将负载电流频谱特性中的n次谐波信号相对幅值作为负载电流频谱特征，其中，n＝1，2，…，M；所述M表示谐波最高次数且M大于等于3。

所述组合分类器中，BP神经网络分类器为主分类器，贝叶斯分类器为辅助分类器。

所述组合分类器进行用电负载类型识别的方法是：当主分类器成功实现用电负载类型识别时，主分类器的用电负载类型识别结果为组合分类器的识别结果；当主分类器未能实现用电负载类型识别，且主分类器的识别结果为2种或者2种以上用电负载类型，将主分类器输出的2种或者2种以上用电负载类型识别结果中，辅助分类器输出中概率最高的用电负载类型作为组合分类器的用电负载类型识别结果；当主分类器未能实现用电负载类型识别，且主分类器的识别结果中未能给出识别的用电负载类型时，将辅助分类器输出中概率最高的用电负载类型作为组合分类器的用电负载类型识别结果。

所述启动电流特征由信息处理模块通过以下方法获得：

步骤1、用电负载启动前，开始对用电负载的负载电流连续采样并对负载电流大小进行判断；当负载电流有效值大于ε时，判定用电负载开始启动并转向步骤2；所述ε为大于0的数值；

步骤2、对用电负载的负载电流进行连续采样，以工频周期为单位计算负载电流有效值并保存；计算最近N个工频周期的负载电流有效值的平均值；当最近N个工频周期之内的每个工频周期的负载电流有效值与该N个工频周期的负载电流有效值的平均值相比较，波动幅度均小于设定的相对误差范围E时，判定用电负载进入稳定状态，转向步骤3；所述N的取值范围为50-500；所述E的取值范围为2％-20％；

步骤3、将最近N个工频周期之内的负载电流有效值的平均值作为用电负载稳态电流；

将用电负载开始启动时刻至最近N个工频周期起始时刻之间的时间作为启动过程时间；计算用电负载开始启动后L个工频周期之内的用电负载电流有效值的平均值与用电负载稳态电流之间的比值，将该比值作为用电负载的启动冲激电流；计算用电负载的启动过程时间之内的用电负载电流有效值的平均值与用电负载稳态电流之间的比值，将该比值作为用电负载的启动平均电流；计算用电负载的启动平均电流与启动过程时间之间的乘积，将该乘积作为用电负载的启动电流冲量；所述L的取值范围为1-5。

所述组合分类器的输入特征还包括用电负载稳态电流。

所述用电负载的基波电压电流相位差通过以下方法获得：

步骤①、待用电负载进入稳定状态后，同步获取用电负载的稳态电压信号、稳态电流信号，并将其转换为相应的稳态电压数字信号、稳态电流数字信号；

步骤②、对稳态电压数字信号、稳态电流数字信号分别进行数字滤波，提取出基波电压信号、基波电流信号；

步骤③、分析计算基波电压信号与基波电流信号之间的相位差，将基波电压信号与基波电流信号之间的相位差作为用电负载的基波电压电流相位差。

本发明的有益效果是：同时采用用电负载的启动电流特征、用电负载的负载电流频谱特征以及用电负载的基波电压电流相位差作为所述用电负载类型识别器的识别特征，特征信息丰富；采用包括BP神经网络分类器和贝叶斯分类器的组合分类器进行识别分类，兼顾BP神经网络分类器和贝叶斯分类器的特点进行综合识别，泛化能力与识别准确率高；提供的包括启动冲激电流、启动平均电流、启动电流冲量在内的启动电流特征获取方法，以及负载电流频谱特征获取方法简单、可靠。

附图说明

图1为本发明用电负载类型识别器实施例的结构示意图；

图2为白炽灯台灯的启动过程电流波形；

图3为电阻炉等电阻性负载的启动过程电流波形；

图4为单相电机类负载的启动过程电流波形；

图5为计算机及开关电源类负载的启动过程电流波形；

图6为用电负载类型识别器进行用电负载类型识别的流程图。

具体实施方式

以下结合附图对本发明作进一步说明。

图1为本发明用电负载类型识别器实施例的结构示意图，包括信息采集模块101、信息处理模块102、通信模块103。

信息采集模块102用于采集用电负载的负载电压、负载电流并将负载电压、负载电流转换成电压数字信号、电流数字信号，电压数字信号、电流数字信号被送至信息处理模块102。信息采集模块中包括电压传感器、电流传感器、前置放大器、滤波器、A/D转换器等组成部分，分别完成负载电压、负载电流信号的传感、放大、滤波与模数转换功能。当负载电流范围较大时，可以选择具有程控功能的前置放大器，或者是在A/D转换器前再增加一个独立的程控放大器，对范围较大的负载电流实行分段控制放大，使输入至A/D转换器的电压信号范围保持在合理的区间，保证转换精度。滤波器用于滤除高频分量，避免频谱混叠。

信息处理模块102依据输入的电压数字信号、电流数字信号，采用包括BP神经网络分类器和贝叶斯分类器的组合分类器实现用电负载类型识别。组合分类器的输入特征包括用电负载的启动电流特征、用电负载的负载电流频谱特征和用电负载的基波电压电流相位差。信息处理模块102的核心为DSP、ARM、单片机，或者为FPGA。当信息处理模块的核心中包括有A/D转换器且该A/D转换器满足要求时，信息采集模块101中的A/D转换器可以采用信息处理模块102的核心中包括的A/D转换器。

通信模块103用于实现与上位机之间的通信，将识别结果发送至上位机。通信模块102与上位机之间的通信方式包括无线通信方式与有线通信方式，可以采用的无线通信方式包括ZigBee、蓝牙、WiFi、433MHz数传等方式，可以采用的有线通信方式包括485总线、CAN总线、互联网、电力载波等方式。通信模块103还可以接收上位机的相关工作指令，完成指定的工作任务。上位机可以是管理部门的服务器，也可以是各种工作站，或者是各种移动终端。

信息采集模块101、信息处理模块102、通信模块103的全部或者部分功能可以集成在一片SoC上，减小识别器体积，方便安装。

不同的用电负载设备具有不同的启动电流特征。如图2所示为白炽灯台灯的启动过程电流波形。白炽灯是将灯丝通电加热到白炽状态，利用热辐射发出可见光的电光源。白炽灯的灯丝通常用耐高温的金属钨制造，但金属钨的电阻随温度变化大，以R_t表示钨丝在t℃时的电阻，以R₀表示钨丝在0℃时的电阻，则两者有下述的关系

R_t＝R₀(1+0.0045t)

例如，设白炽灯的灯丝(钨丝)在正常工作时的温度为2000℃，一只“220V 100W”

的白炽灯的灯丝在2000℃正常工作时的电阻为

$R_t=\frac{U^2}{P}=\frac{220\×220}{100}=484\mathrm{\Ω}$

其在不通电时0℃的电阻为

$R_0=\frac{R_t}{1+0.0045t}=\frac{484}{1+0.0045\×2000}=48.4\mathrm{\Ω}$

其在不通电时20℃的电阻为

R₂₀＝R₀(1+0.0045t)＝52.8Ω

即白炽灯在启动通电的瞬间电流超过其额定电流的9倍，且最大启动电流发生在启动时刻。随着白炽灯钨丝温度的升高，白炽灯的负载电流按照指数规律减小，然后进入稳定状态。

设用电负载稳态电流有效值为I_W，且定义用电负载电流有效值进入用电负载稳态电流有效值的一个设定的相对误差范围之内并稳定在这个相对误差范围之内，则用电负载进入稳定状态。相对误差范围可以设定为10％，也可以设定为2％、5％、15％、20％等2％-20％之间的值。图2中，设定的相对误差范围为10％，当白炽灯的负载电流按照指数规律减小到其I_W的10％误差范围时，如图2中的时刻T_S，启动过程结束。白炽灯的启动过程时间为T_S。I_W为有效值。

选择启动冲激电流I_G、启动平均电流I_D、启动电流冲量Q_I作为用电负载的启动电流特征；启动冲激电流I_G、启动平均电流I_D均为标么值。具体定义是：启动冲激电流I_G为用电负载启动开始后T₂时间之内的用电负载电流平均值与用电负载稳态电流I_W的比值；启动平均电流I_D为用电负载启动时间T_S之内的用电负载电流平均值与用电负载稳态电流I_W的比值；启动电流冲量Q_I为启动平均电流I_D与启动过程时间T_S的乘积，量纲为ms。用电负载电流、用电负载稳态电流均为有效值。T₂的取值范围为20-100ms，或者是1-5个工频周期；例如，T₂取值40ms，即2个工频周期。启动冲激电流I_G反映的是用电负载启动后短时间内的电流冲激大小。在部分用电负载的启动过程中，当有用电负载的实际启动过程时间T_S小于设定的T₂时，令用电负载的启动过程时间T_S等于T₂。启动平均电流I_D反映的是用电负载启动过程中的电流整体大小。启动电流冲量Q_I反映的是用电负载启动的整体强度。

图2中，白炽灯的启动冲激电流I_G为T₀(白炽灯启动时刻，电流为I₀)至T₂(设定的时刻，电流为I₂)之间白炽灯的电流平均值与白炽灯的稳态电流I_W的比值。启动平均电流I_D为T₀(白炽灯启动时刻)至T_S(白炽灯启动过程结束时间)之间白炽灯的电流平均值与白炽灯的稳态电流I_W的比值。启动电流冲量Q_I为白炽灯启动平均电流I_D与启动过程时间T_S的乘积。

如图3所示为电阻炉等电阻性负载的启动过程电流波形。电阻炉等电阻性负载通常采用镍铬、铁铬铝等电热合金丝，其共同特点是电阻温度修正系数小，电阻值稳定。以牌号为Cr20Ni80的镍铬电热丝为例，其在1000℃时的电阻修正系数为1.014，即1000℃时相对于20℃时，牌号为Cr20Ni80的镍铬电热丝电阻只增加1.4％。电阻炉等电阻性负载在通电启动时即进入稳定状态，电阻炉等电阻性负载的实际启动过程时间T_S＝0，因此，令电阻炉等电阻性负载的实际启动过程时间T_S＝T₂；例如，当T₂设定为40ms时，则此时的启动过程时间T_S也为40ms。由于电阻性负载T₀时刻电流I₀、T₂时刻电流I₂与电阻性负载的稳态电流I_W相等，因此，电阻性负载的启动冲激电流I_G＝1，启动平均电流I_D＝1。

如图4所示为单相电机类负载的启动过程电流波形。单相电机类负载既具有电感性负载特性，又具有反电动势负载特性。启动时刻，由于电感的作用，启动时刻的启动电流I₀为0；随后电流迅速上升，在电机反电动势未建立之前，达到电流峰值I_M；此后，电机转速增加，电机负载电流逐步减小，直到进入稳定状态。图4中，单相电机类负载的启动冲激电流I_G为T₀(单相电机类负载启动时刻，电流为I₀)至T₂(设定的时刻，电流为I₂)之间单相电机类负载的电流平均值与稳态电流I_W的比值。启动平均电流I_D为T₀(单相电机类负载启动时刻)至T_S(单相电机类负载启动过程结束时间)之间单相电机类负载的电流平均值与稳态电流I_W的比值。启动电流冲量Q_I为单相电机类负载启动平均电流I_D与启动过程时间T_S的乘积。

如图5所示为计算机及开关电源类负载的启动过程电流波形。计算机及开关电源类负载因为对电容充电的影响，在启动瞬间会产生一个很大的浪涌电流，其峰值可达到稳态电流有效值I_W的几倍至十几倍，时间为1至2个工频周期。由于计算机及开关电源类负载的启动时间短，其启动过程时间T_S有可能小于设定的T₂；当其启动过程时间T_S小于设定的T₂时，令T_S等于T₂。图5中，计算机及开关电源类负载的启动冲激电流I_G为T₀(计算机及开关电源类负载启动时刻，电流为I₀)至T₂(设定的时刻，电流为I₂)之间计算机及开关电源类负载的电流平均值与稳态电流I_W的比值。启动平均电流I_D为T₀(计算机及开关电源类负载启动时刻)至T_S(计算机及开关电源类负载启动过程结束时间)之间计算机及开关电源类负载的电流平均值与稳态电流I_W的比值。启动电流冲量Q_I为计算机及开关电源类负载启动平均电流I_D与启动过程时间T_S的乘积。

获取用电负载的启动电流特征的方法是：

用电负载启动前，负载电流值为0(未开机)或者很小(处于待机状态)时，信息处理模块102即开始对负载电流进行连续采样；当采样得到的负载电流值有效值开始大于0或者是开始大于用电负载的待机电流时，即判断出用电负载已经启动，记录该时刻为T₀。用一个较小的非负阈值ε来区分用电负载启动前后的负载电流值，当ε取值特别小时，例如，ε取值1mA时，所述识别器不考虑待机情况，即认为待机也是用电负载的启动状态；当ε取值较小但大于用电负载的待机电流时，例如，ε取值20mA时，所述识别器会将用电负载的待机状态认为是未启动状态，但同时也会的部分功率特别小的用电负载造成漏识别。

信息处理模块102对负载电流进行连续采样，且以工频周期为单位计算负载电流有效值并保存；当用电负载已经启动，且连续采样达到N个工频周期后，采样的同时连续计算最近N个工频周期的负载电流有效值的平均值I_V；信息处理模块102对最近N个工频周期之内每个工频周期的负载电流有效值与该N个工频周期的负载电流有效值的平均值进行比较，误差(或波动)幅度均小于设定的相对误差范围E时，判定用电负载进入稳定状态，该最近N个工频周期的起始时刻为启动过程的结束时刻，记录该时刻为T₁(如图2-图5所示)。

将最近N个工频周期之内的负载电流有效值的平均值作为电器负载稳态电流I_W；将电器开始启动时刻T₀至最近N个工频周期起始时刻T₁之间的时间作为启动过程时间T_S。计算T₀至设定的T₂之间(即电器开始启动后1-5个工频周期之内)的负载电流平均值与稳态电流I_W的比值，将该比值作为电器的启动冲激电流I_G。计算T₀至T_S之间的负载电流平均值与稳态电流I_W的比值，将该比值作为电器的启动平均电流I_D。计算电器的启动平均电流I_D与启动过程时间T_S的乘积，将该乘积作为电器的启动电流冲量Q_I。

由于预先不知道用电负载稳态电流有效值I_W，因此，将N个工频周期，即一段持续时间T_P之内波动范围小于设定的相对误差范围E时的负载电流有效值的平均值作为用电负载稳态电流有效值I_W。由于普通用电负载的启动过程较快，所以，T_P的取值范围为1-10s，典型取值是2s，相应的工频周期数量N的取值范围为50-500，N的典型取值是100。所述相对误差范围E的取值范围为2％-20％，E的典型取值是10％。

组合分类器的输入特征还包括用电负载的负载电流频谱特征。用电负载的负载电流频谱特征由信息处理模块102控制信息采集模块101，通过以下步骤获得：

步骤一、待用电负载进入稳定状态后，获取用电负载的稳态电流信号，并将其转换为对应的稳态电流数字信号。

步骤二、对稳态电流数字信号进行傅立叶变换，得到负载电流频谱特性。为保证傅立叶变换的顺利进行，在前述获取用电负载的稳态电流信号，并将其转换为对应的稳态电流数字信号的过程中，A/D转换器的精度和速度需要满足傅立叶变换的要求，采样频率可以设定为10kHz，或者是其他数值；信息处理模块102对采集到的稳态电流数字信号进行FFT运算，计算其频谱。

步骤三、将负载电流频谱特性中的n次谐波信号相对幅值作为负载电流频谱特征，其中，n＝1，2，…，M；在组成组合分类器的输入特征向量时，n次谐波信号相对幅值在输入特征向量中按照1，2，…，M的顺序依次排列。由于负载电流频谱特性主要由奇次谐波组成，除少数用电负载设备外，偶次谐波分量几乎为0，因此，也可以将负载电流频谱特性中谐波次数为n的次奇次谐波信号相对幅值作为负载电流频谱特征，其中，n＝1，3，…，M。所述谐波信号相对幅值为谐波信号幅值与用电负载稳态电流有效值I_W的比值。n＝1时的1次谐波为工频基波。所述M表示谐波最高次数，一般情况下，M大于等于3。

组合分类器的输入特征还包括用电负载的基波电压电流相位差。基波电压电流相位差可以对电阻性、电容性、电感性负载进行区分，还可以对一般的电感性负载和大电感性负载进行区分。用电负载的基波电压电流相位差由信息处理模块102控制信息采集模块101，通过以下步骤获得：

步骤①、待用电负载进入稳定状态后，获取用电负载的稳态电压信号、稳态电流信号，并将其转换为对应的稳态电压数字信号、稳态电流数字信号。

步骤②、对稳态电压数字信号、稳态电流数字信号分别进行数字滤波，提取出基波电压信号、基波电流信号。

步骤③、分析计算基波电压信号与基波电流信号之间的相位差，将基波电压信号与基波电流信号之间的相位差作为用电负载的基波电压电流相位差。

步骤②中对稳态电压数字信号、稳态电流数字信号分别进行数字滤波，其数字滤波算法可以选择卡尔曼滤波法、小波变换法、维纳滤波法、自适应滤波等数字滤波器算法。

组合分类器中，BP神经网络分类器为主分类器，贝叶斯分类器为辅助分类器。组合分类器的输入特征包括前述的启动电流特征和负载电流频谱特征，组合分类器的输入特征同时作为BP神经网络分类器的输入特征和贝叶斯分类器的输入特征。

如图6所示为用电负载类型识别器进行用电负载类型识别的流程图，用电负载类型识别器进行用电负载类型识别的方法是：

步骤A、等待用电负载启动；

步骤B、采集用电负载启动电流数据并保存，直至用电负载启动过程结束；

步骤C、分析采集的用电负载启动电流数据，获取用电负载的启动电流特征；

步骤D、采集用电负载稳态工作时的电压、电流数据并保存；

步骤E、分析采集的用电负载稳态工作时的电压、电流数据，获取用电负载的负载电流频谱特征、基波电压电流相位差；

步骤F、将启动电流特征、负载电流频谱特征、基波电压电流相位差作为组合分类器的输入特征；组合分类器进行用电负载类型识别；

步骤G、输出用电负载类型识别结果。

所述组合分类器进行用电负载类型识别的方法是：当主分类器成功实现用电负载类型识别，即主分类器输出的识别结果为唯一的用电负载类型，即识别结果中唯一的用电负载类型为是时，将主分类器识别的用电负载类型作为组合分类器的用电负载类型识别结果；当主分类器未能实现用电负载类型识别，且主分类器的识别结果为2种或者2种以上用电负载类型，即识别结果中有2种或者2种以上用电负载类型为是时，将主分类器输出的2种或者2种以上用电负载类型识别结果中，辅助分类器输出中概率最高的用电负载类型作为组合分类器的用电负载类型识别结果；当主分类器未能实现用电负载类型识别，且主分类器的识别结果中未能给出识别的用电负载类型，即识别结果中没有用电负载类型为是时，将辅助分类器输出中概率最高的用电负载类型作为组合分类器的用电负载类型识别结果。

以一个简单的实施例1为例，来说明组合分类器进行用电负载类型识别的方法。设有一个组合分类器，其输入特征为x＝{I_G，I_D，Q_I，A₁，A₂，A₃，A₄，A₅，}，其中，I_G是启动冲激电流；I_D是启动平均电流；Q_I是启动电流冲量；A₁、A₂、A₃、A₄、A₅为负载电流频谱特性中的1-5次谐波信号相对幅值，为用电负载的基波电压电流相位差。组合分类器的输出是{B₁，B₂，B₃，B₄}，B₁、B₂、B₃、B₄分别代表组合分类器对白炽灯、电阻炉、吹风机、计算机的识别结果输出，识别结果B₁、B₂、B₃、B₄的取值均为二值分类标记。主分类器的输入特征也是x＝{I_G，I_D，Q_I，A₁，A₂，A₃，A₄，A₅，}，其输出是{F₁，F₂，F₃，F₄}，F₁、F₂、F₃、F₄分别代表主分类器对白炽灯、电阻炉、吹风机、计算机的识别结果输出，识别结果F₁、F₂、F₃、F₄的取值也均为二值分类标记。辅助分类器的输入特征同样为x＝{I_G，I_D，Q_I，A₁，A₂，A₃，A₄，A₅，}，其输出是{P(y₁|x)，P(y₂|x)，P(y₃|x)，P(y₄|x)}，P(y₁|x)、P(y₂|x)、P(y₃|x)、P(y₄|x)为辅助分类器输出的后验概率，P(y₁|x)、P(y₂|x)、P(y₃|x)、P(y₄|x)之间的相互大小表明辅助分类器的当前输入特征表示所识别的用电负载属于白炽灯、电阻炉、吹风机、计算机的可能性大小。

在实施例1中，B₁、B₂、B₃、B₄的分类标记和F₁、F₂、F₃、F₄的分类标记均取1、0。分类标记为1时，相应的用电负载类型与当前输入特征匹配，为确认的识别结果，或者说相应的用电负载类型识别结果为是；分类标记为0时，相应的用电负载类型与输入特征不匹配，未能成为确认的识别结果，或者说相应的用电负载类型识别结果为否。

在实施例1中，设某次的主分类器的识别结果分类标记为F₁F₂F₃F₄＝0100，则认为主分类器成功实现用电负载类型识别，因此，不考虑辅助分类器的识别结果，直接令B₁B₂B₃B₄＝0100，即组合分类器的识别结果是：被识别的用电负载为电阻炉。

在实施例1中，设某次的主分类器的识别结果分类标记为F₁F₂F₃F₄＝1010，则认为主分类器未能实现用电负载类型识别，且主分类器的识别结果为2种或者2种以上用电负载类型；再设此时辅助分类器的识别结果满足P(y₁|x)<P(y₃|x)，则令B₁B₂B₃B₄＝0010，即组合分类器的识别结果是：被识别的用电负载为吹风机。

在实施例1中，设某次的主分类器的识别结果分类标记为F₁F₂F₃F₄＝0000，则认为主分类器未能实现用电负载类型识别，且主分类器的识别结果中未能给出识别的用电负载类型；再设此时辅助分类器的识别结果满足P(y₁|x)>P(y₂|x)且P(y₁|x)>P(y₃|x)且P(y₁|x)>P(y₄|x)，则令B₁B₂B₃B₄＝1000，即组合分类器的识别结果是：被识别的用电负载为白炽灯。

组合分类器、主分类器的识别结果分类标记也可以采用其他的方案，例如，分别用分类标记1、-1，或者是0、1，或者是-1、1，以及其他方案来表示相应用电负载识别结果为是、否。组合分类器与主分类器的分类标记方案可以相同，也可以不相同。

所述组合分类器的输入特征中，还可以包括用电负载稳态电流有效值I_W。例如，有2种不同的用电负载，电烙铁和电阻炉需要识别，电烙铁、电阻炉都是纯电阻负载，且都具有电阻温度修正系数小，电阻值稳定的共同特点。因此，单纯依靠前述的启动电流特征和负载电流频谱特征、用电负载的基波电压电流相位差特征无法将他们进行区分。输入特征中增加用电负载稳态电流有效值I_W后，电烙铁功率小，用电负载稳态电流有效值I_W小；电阻炉功率大，用电负载稳态电流有效值I_W大，特征不同，组合分类器可以进行并完成识别。

辅助分类器为贝叶斯分类器。可以选择NBC分类器(朴素贝叶斯分类器)、TAN分类器(树扩展朴素贝叶斯分类器)、BAN分类器(增强的贝叶斯分类器)等三种贝叶斯分类器之中的一种作为辅助分类器。

实施例2选择NBC分类器作为辅助分类器。朴素贝叶斯分类的定义如下：

⑴设x＝{a₁,a₂,…,a_m}为一个待分类项，而每个a为x的一个特征属性；

⑵有类别集合C＝{y₁,y₂,…,y_n}；

⑶计算P(y₁|x),P(y₂|x),…,P(y_n|x)；

⑷如果P(y_k|x)＝max{P(y₁|x),P(y₂|x),…,P(y_n|x)}，则x∈y_k。

计算第⑶步中的各个条件概率的具体方法是：

①找到一个已知分类的待分类项集合作为训练样本集；

②统计得到各类别下各个特征属性的条件概率估计；

P(a₁|y₁),P(a₂|y₁),…,P(a_m|y₁)；

P(a₁|y₂),P(a₂|y₂),…,P(a_m|y₂)；

…；

P(a₁|y_n),P(a₂|y_n),…,P(a_m|y_n)。

③根据贝叶斯定理，有：

$P\left(y_i\right|x)=\frac{P\left(x\right|y_i)P\left(y_i\right)}{P\left(x\right)}---\left(1\right)$

因为分母对于所有类别为常数，因此我们只要将分子最大化即可；又因为在朴素贝叶斯中各特征属性是条件独立的，所以有：

$P\left(x\right|y_i)P\left(y_i\right)=P\left(a_1\right|y_i)P\left(a_2\right|y_i)...P\left(a_m\right|y_i)P\left(y_i\right)=P\left(y_i\right)\underset{j=1}{\overset{m}{\&Pi;}}P\left(a_j\right|y_i)$

实施例2中，组合分类器的输入特征是{I_G，I_D，Q_I，A₁，A₃，I_W}，其中，I_G是启动冲激电流；I_D是启动平均电流；Q_I是启动电流冲量；A₁、A₃为负载电流频谱特性中的1、3次奇次谐波信号相对幅值；为用电负载的基波电压电流相位差，单位为度，且基波电压超前于基波电流时，I_W为用电负载稳态电流有效值，单位是安培。要求识别的用电负载类别是白炽灯、电阻炉、电风扇、计算机、电烙铁。令朴素贝叶斯分类器的特征属性组合x＝{a₀，a₀,a₃,a₄,a₅,a₆,a₇}中的元素与组合分类器的输入特征集合中的元素按序{I_G，I_D，Q_I，A₁，A₃，I_W}一一对应；朴素贝叶斯分类器的输出类别集合C＝{y₁,y₂,y₃,y₄,y₅}则分别与用电负载类别白炽灯、电阻炉、电风扇、计算机、电烙铁一一对应。

训练NBC分类器的过程包括：

1、对特征属性进行分段划分，进行离散化处理。实施例2中，采取的特征属性离散化方法是：

a₁：{a₁<3.9,3.9≤a₁≤6.5,a₁>6.5}；

a₂：{a₂<1.2,1.2≤a₂≤2.6,a₂>2.6}；

a₃：{a₃<150,150≤a₃≤600,a₃>600}；

a₄：{a₄<0.7,0.7≤a₄≤0.9,a₄>0.9}；

a₅：{a₅<0.02,0.02≤a₅≤0.05,a₅>0.05}；

a₆：{a₆<-6,-6≤a₆≤18,a₆>18}；

a₇：{a₇<0.45,a₇≥0.45}。

2、对每类用电负载类型均采集多组样本作为训练样本，同时计算每类用电负载类型样本在所有用电负载类型样本中所占有的比例，即分别计算P(y₁)、P(y₂)、P(y₃)、P(y₄)、P(y₅)。当每类用电负载均采集相同的样本数量时，例如，每类用电负载均采集超过100组的样本，其中每类用电负载随机选择100组样本作为训练样本，其他则作为测试样本，总的训练样本为500组，且有

P(y₁)＝P(y₂)＝P(y₃)＝P(y₄)＝P(y₅)＝0.2。

3、计算训练样本每个类别条件下各个特征属性分段的频率(比例)，统计得到各类别下各个特征属性的条件概率估计，即分别统计计算

P(a₁<3.9|y₁)、P(3.9≤a₁≤6.5|y₁)、P(a₁>6.5|y₁)；

P(a₁<3.9|y₂)、P(3.9≤a₁≤6.5|y₂)、P(a₁>6.5|y₂)；

…；

P(a₁<3.9|y₅)、P(3.9≤a₁≤6.5|y₅)、P(a₁>6.5|y₅)；

P(a₂<1.2|y₁)、P(1.2≤a₂≤2.6|y₁)、P(a₂>2.6|y₁)；

P(a₂<1.2|y₂)、P(1.2≤a₂≤2.6|y₂)、P(a₂>2.6|y₂)；

…；

P(a₂<1.2|y₅)、P(1.2≤a₂≤2.6|y₅)、P(a₂>2.6|y₅)；

P(a₃<150|y₁)、P(150≤a₃≤600|y₁)、P(a₃>600|y₁)；

P(a₃<150|y₂)、P(150≤a₃≤600|y₂)、P(a₃>600|y₂)；

…；

P(a₃<150|y₅)、P(150≤a₃≤600|y₅)、P(a₃>600|y₅)；

P(a₄<0.7|y₁)、P(0.7≤a₄≤0.9|y₁)、P(a₄>0.9|y₁)；

P(a₄<0.7|y₂)、P(0.7≤a₄≤0.9|y₂)、P(a₄>0.9|y₂)；

…；

P(a₄<0.7|y₅)、P(0.7≤a₄≤0.9|y₅)、P(a₄>0.9|y₅)；

P(a₅<0.02|y₁)、P(0.02≤a₅≤0.05|y₁)、P(a₅>0.05|y₁)；

P(a₅<0.02|y₂)、P(0.02≤a₅≤0.05|y₂)、P(a₅>0.05|y₂)；

P(a₅<0.02|y₅)、P(0.02≤a₅≤0.05|y₅)、P(a₅>0.05|y₅)；

P(a₆<-6|y₁)、P(-6≤a₆≤18|y₁)、P(a₆>18|y₁)；

P(a₆<-6|y₂)、P(-6≤a₆≤18|y₂)、P(a₆>18|y₂)；

…；

P(a₆<-6|y₅)、P(-6≤a₆≤18|y₅)、P(a₆>18|y₅)；

P(a₇<0.45|y₁)、P(a₇≥0.45|y₁)；

P(a₇<0.45|y₂)、P(a₇≥0.45|y₂)；

…；

P(a₇<0.45|y₅)、P(a₇≥0.45|y₅)。

经过上述的步骤1、步骤2、步骤3，NBC分类器训练完成。其中，步骤1对特征属性进行分段划分由人工确定，对每一个输入特征进行分段离散化时，分段的数量为2段或者2段以上，例如，实施例2中，特征a₁-a₆都分为3段，特征a₇分为2段。每一个特征具体分为多少段，分段阈值的选择可以根据训练后的贝叶斯分类器对测试样本测试后的结果进行调整。步骤2、步骤3由信息处理模块102或者是计算机计算完成。

本发明中采用贝叶斯分类器进行分类的方法是：

1、将组合分类器的输入特征作为贝叶斯分类器的输入特征。在实施例2中，将组合分类器的输入特征集合{I_G，I_D，Q_I，A₁，A₃，I_W}作为贝叶斯分类器的输入特征x，且有x＝{a₁,a₂,a₃,a₄,a₅,a₆,a₇}。

2、根据训练得到的各类别下各个特征属性的条件概率估计，分别确定各输入特征属性的分段所在并确定其对每类用电负载类别的概率P(a₁|y₁)～P(a_m|y_n)，其中，用电负载类别集合为C＝{y₁,y₂,…,y_n}。实施例2中，用电负载类别集合C＝{y₁,y₂,y₃,y₄,y₅}对应代表的用电负载类别是白炽灯、电阻炉、电风扇、计算机、电烙铁，确定P(a₁|y₁)～P(a₇|y₅)的方法是采用训练NBC分类器过程中得到的各个特征属性的条件概率估计。

3、按照式

$P\left(y_i\right|x)=\frac{P\left(x\right|y_i)P\left(y_i\right)}{P\left(x\right)}$

计算每种用电负载类别的后验概率。因为分母P(x)对于所有用电负载类别为常数，令P(x)＝1替代实际的P(x)值，不影响每种用电负载类别后验概率之间的相互大小比较，此时有

$P\left(y_i\right|x)=P\left(x\right|y_i)P\left(y_i\right)=P\left(y_i\right)\underset{j=1}{\overset{m}{\&Pi;}}P\left(a_j\right|y_i)$

实施例2中，有

$P\left(y_1\right|x)=P\left(x\right|y_1)P\left(y_1\right)=P\left(y_1\right)\underset{j=1}{\overset{7}{\&Pi;}}P\left(a_j\right|y_1);$

$P\left(y_2\right|x)=P\left(x\right|y_2)P\left(y_2\right)=P\left(y_2\right)\underset{j=1}{\overset{7}{\&Pi;}}P\left(a_j\right|y_2);$

$P\left(y_3\right|x)=P\left(x\right|y_3)P\left(y_3\right)=P\left(y_3\right)\underset{j=1}{\overset{7}{\&Pi;}}P\left(a_j\right|y_3);$

$P\left(y_4\right|x)=P\left(x\right|y_4)P\left(y_4\right)=P\left(y_4\right)\underset{j=1}{\overset{7}{\&Pi;}}P\left(a_j\right|y_4);$

$P\left(y_5\right|x)=P\left(x\right|y_5)P\left(y_5\right)=P\left(y_5\right)\underset{j=1}{\overset{7}{\&Pi;}}P\left(a_j\right|y_5).$

采用测试样本对训练好的贝叶斯分类器进行测试，根据测试结果决定是否调整对输入特征的离散化方法(即调整分段数量与阈值)，重新训练贝叶斯分类器。

主分类器为BP神经网络分类器，选择3层BP神经网络分类器作为主分类器。将BP神经网络分类器输入特征向量中元素的数量，即输入特征的数量作为输入层的节点数，例如，实施例1中的9个，或者是实施例2中的7个。将需要识别的用电负载类型的数量作为输出层节点数，例如，实施例1中，输出层节点为4个，分别输出识别白炽灯、电阻炉、吹风机、计算机的结果；实施例2中，输出层节点为5个，分别输出识别白炽灯、电阻炉、电风扇、计算机、电烙铁的结果。中间隐层的节点数量根据经验来取，例如，实施例1、实施例2中隐层的节点数量可以在6-18的范围内选取。对每类用电负载类型均采集多组样本，例如，均采集200组样本；随机选取其中的若干组，例如150组样本作为训练样本，剩余的作为测试样本，对BP神经网络分类器进行训练与测试。多输入、多输出的3层BP神经网络分类器由于多输出之间的耦合作用，有可能在训练或者测试时不能对样本进行完全识别；即使是在训练或者测试时能够对样本进行完全识别，受到泛化能力的制约，对新输入的某一特征属性进行识别时，主分类器有可能输出的识别结果为唯一的用电负载类型，或者识别结果为2种或者2种以上用电负载类型，或者未能给出识别的用电负载类型。

主分类器还可以选择多个单节点输出的3层BP神经网络分类器共同组成，每个单节点输出的3层BP神经网络分类器对应识别一种用电负载类型，例如，实施例1中可以采用4个单节点输出的3层BP神经网络分类器分别识别白炽灯、电阻炉、吹风机、计算机；实施例2中可以采用5个单节点输出的3层BP神经网络分类器分别识别白炽灯、电阻炉、电风扇、计算机、电烙铁。主分类器选择多个单节点输出的3层BP神经网络分类器共同组成时，所有单节点输出的3层BP神经网络分类器的输入层节点数为主分类器输入特征向量中元素的数量；中间隐层的节点数量根据经验来取，各单节点输出的3层BP神经网络分类器中间隐层的节点数量可以相同，也可以不同，按照各自的需要选择。与非单节点输出的神经网络一样，需要对每类用电负载类型均采集多组样本，例如，均采集200组样本；随机选取其中的若干组，例如150组样本作为训练样本，剩余的作为测试样本，对每个单节点输出的BP神经网络分类器进行训练与测试。主分类器选择多个单节点输出的3层BP神经网络分类器共同组成时，每个单节点输出的3层BP神经网络分类器只需要完成一种用电负载类型的识别，每个网络的训练相对简单。由于此时主分类器由多个单节点输出的3层BP神经网络分类器组成，各单节点输出的3层BP神经网络分类器之间相互独立，因此，对某一特征属性进行识别时，主分类器有可能输出的识别结果为唯一的用电负载类型，或者识别结果为2种或者2种以上用电负载类型，或者未能给出识别的用电负载类型。

BP神经网络分类器的训练方法可以采用梯度下降法，也可以采用粒子群优化、遗传算法等优化方法。样本采集采用前述的获取用电负载的启动电流特征的方法和获取用电负载的负载电流频谱特征以及获取用电负载的基波电压电流相位差特征的方法。

Claims (10)

1.一种用电负载类型识别器，其特征在于，包括信息采集模块、信息处理模块、通信模块；

所述信息采集模块用于采集用电负载的负载电流并转换成电流数字信号；所述电流数字信号被送至信息处理模块；

所述信息处理模块依据输入的电流数字信号，采用组合分类器进行用电负载类型识别；

所述通信模块用于发送信息处理模块的用电负载类型识别结果至上位机；

所述组合分类器的输入特征包括用电负载的启动电流特征、用电负载的负载电流频谱特征和用电负载的基波电压电流相位差；

所述组合分类器包括BP神经网络分类器和贝叶斯分类器；

所述启动电流特征包括启动冲激电流、启动平均电流、启动电流冲量。

2.如权利要求1所述的用电负载类型识别器，其特征在于，所述信息采集模块包括电流传感器、前置放大器、滤波器、A/D转换器；所述信息处理模块的核心为DSP，或者为ARM，或者为单片机，或者为FPGA。

3.如权利要求1所述的用电负载类型识别器，其特征在于，所述信息采集模块、信息处理模块、通信模块的全部或者部分功能集成在一片SoC上。

4.如权利要求1所述的用电负载类型识别器，其特征在于，所述通信模块还接收上位机的相关工作指令；所述通信模块与上位机之间的通信方式包括无线通信方式与有线通信方式；所述无线通信方式包括ZigBee、蓝牙、WiFi、433MHz数传方式；所述有线通信方式包括485总线、CAN总线、互联网、电力载波方式。

5.如权利要求1-4中任一项所述的用电负载类型识别器，其特征在于，所述组合分类器中，BP神经网络分类器为主分类器，贝叶斯分类器为辅助分类器。

6.如权利要求5所述的用电负载类型识别器，其特征在于，所述组合分类器进行用电负载类型识别的方法是：当主分类器成功实现用电负载类型识别时，主分类器的用电负载类型识别结果为组合分类器的识别结果；当主分类器未能实现用电负载类型识别，且主分类器的识别结果为2种或者2种以上用电负载类型，将主分类器输出的2种或者2种以上用电负载类型识别结果中，辅助分类器输出中概率最高的用电负载类型作为组合分类器的用电负载类型识别结果；当主分类器未能实现用电负载类型识别，且主分类器的识别结果中未能给出识别的用电负载类型时，将辅助分类器输出中概率最高的用电负载类型作为组合分类器的用电负载类型识别结果。

7.如权利要求5所述的用电负载类型识别器，其特征在于，所述负载电流频谱特征通过以下方法获得：

步骤一、获取用电负载的稳态电流信号，并将其转换为对应的稳态电流数字信号；

步骤二、对稳态电流数字信号进行傅立叶变换，得到负载电流频谱特性；

步骤三、将负载电流频谱特性中的n次谐波信号相对幅值作为负载电流频谱特征，其中，n＝1，2，…，M；所述M表示谐波最高次数且M大于等于3。

8.如权利要求5所述的用电负载类型识别器，其特征在于，所述启动电流特征由信息处理模块通过以下方法获得：

步骤1、用电负载启动前，开始对用电负载的负载电流连续采样并对负载电流大小进行判断；当负载电流有效值大于ε时，判定用电负载开始启动并转向步骤2；所述ε为大于0的数值；

步骤2、对用电负载的负载电流进行连续采样，以工频周期为单位计算负载电流有效值并保存；计算最近N个工频周期的用电负载电流有效值的平均值；当最近N个工频周期之内的每个工频周期的负载电流有效值与该N个工频周期的负载电流有效值的平均值相比较，波动幅度均小于设定的相对误差范围E时，判定用电负载进入稳定状态，转向步骤3；所述N的取值范围为50-500；所述E的取值范围为2％-20％；

步骤3、将最近N个工频周期之内的负载电流有效值的平均值作为用电负载稳态电流；将用电负载开始启动时刻至最近N个工频周期起始时刻之间的时间作为启动过程时间；计算用电负载开始启动后L个工频周期之内的用电负载电流有效值的平均值与用电负载稳态电流之间的比值，将该比值作为用电负载的启动冲激电流；计算用电负载的启动过程时间之内的用电负载电流有效值的平均值与用电负载稳态电流之间的比值，将该比值作为用电负载的启动平均电流；计算用电负载的启动平均电流与启动过程时间之间的乘积，将该乘积作为用电负载的启动电流冲量；所述L的取值范围为1-5。

9.如权利要求8所述的用电负载类型识别器，其特征在于，所述组合分类器的输入特征还包括用电负载稳态电流。

10.如权利要求5所述的用电负载类型识别器，其特征在于，所述用电负载的基波电压电流相位差通过以下方法获得：

步骤①、待用电负载进入稳定状态后，同步获取用电负载的稳态电压信号、稳态电流信号，并将其转换为相应的稳态电压数字信号、稳态电流数字信号；

步骤②、对稳态电压数字信号、稳态电流数字信号分别进行数字滤波，提取出基波电压信号、基波电流信号；

步骤③、分析计算基波电压信号与基波电流信号之间的相位差，将基波电压信号与基波电流信号之间的相位差作为用电负载的基波电压电流相位差。