本发明涉及一种电器设备识别及分类方法，尤其是涉及一种学生宿舍用电器类型分类方法。
背景技术：
目前，主流的电器性质或者电器类型识别方法包括基于负载功率综合系数算法的电器识别方法、基于电磁感应的电器识别方法、基于神经网络算法的电器识别方法、基于周期性离散变换算法的电器识别方法等。各种方法均能够在一定程度是实现电器性质的识别，但由于特征性质单一，识别手段单一，普遍存在泛化能力不够及不能完全准确识别的问题。
技术实现要素：
本发明的目的在于，针对现在已有技术的缺陷，提供一种能够实现高效识别的学生宿舍用电器类型分类方法。所述学生宿舍用电器类型分类方法采用包括支持向量机分类器和贝叶斯分类器的组合分类器进行电器类型识别；所述组合分类器的输入特征包括电器的启动电流特征、电器的负载电流频谱特征和电器的基波电压电流相位差。
所述组合分类器中，支持向量机分类器为主分类器，贝叶斯分类器为辅助分类器。所述组合分类器进行电器类型识别的方法是：当主分类器成功实现电器类型识别时，主分类器的电器类型识别结果为组合分类器的识别结果；当主分类器未能实现电器类型识别，且主分类器的识别结果为2种或者2种以上电器类型，将主分类器输出的2种或者2种以上电器类型识别结果中，辅助分类器输出中概率最高的电器类型作为组合分类器的电器类型识别结果；当主分类器未能实现电器类型识别，且主分类器的识别结果中未能给出识别的电器类型时，将辅助分类器输出中概率最高的电器类型作为组合分类器的电器类型识别结果。
所述负载电流频谱特征通过以下方法获得：
步骤一、获取电器的稳态电流信号，并将其转换为对应的稳态电流数字信号；
步骤二、对稳态电流数字信号进行傅立叶变换，得到负载电流频谱特性；
步骤三、将负载电流频谱特性中谐波次数为n次的奇次谐波信号相对幅值作为负载电流频谱特征，其中，n＝1，3，…，m；所述m表示谐波最高次数且m大于等于3。所述谐波信号相对幅值为谐波信号幅值与电器负载稳态电流有效值的比值。
所述电器的基波电压电流相位差通过以下方法获得：
步骤①、待电器进入稳定状态后，同步获取电器的稳态电压信号、稳态电流信号，并将其转换为相应的稳态电压数字信号、稳态电流数字信号；
步骤②、对稳态电压数字信号、稳态电流数字信号分别进行数字滤波，提取出基波电压信号、基波电流信号；
步骤③、分析计算基波电压信号与基波电流信号之间的相位差，将基波电压信号与基波电流信号之间的相位差作为电器的基波电压电流相位差。
所述待电器进入稳定状态后的判断，根据最近n个工频周期之内的每个工频周期的负载电流有效值来进行，具体方法是：对电器的负载电流进行连续采样，以工频周期为单位计算负载电流有效值并保存；计算最近n个工频周期的负载电流有效值的平均值；当最近n个工频周期之内的每个工频周期的负载电流有效值与该n个工频周期的负载电流有效值的平均值相比较，波动幅度均小于设定的相对误差范围e时，判定电器进入稳定状态。所述n的取值范围为50-500，e的取值范围为2％-20％。
所述启动电流特征包括启动冲激电流、启动平均电流、启动电流冲量，通过以下方法获得：
步骤1、电器启动前，开始对电器的负载电流连续采样并对负载电流大小进行判断；当负载电流有效值大于ε时，判定电器开始启动并转向步骤2；所述ε为大于0的数值；
步骤2、待判定电器进入稳定状态后，转向步骤3；
步骤3、将最近n个工频周期之内的负载电流有效值的平均值作为电器负载稳态电流；将电器开始启动时刻至最近n个工频周期起始时刻之间的时间作为启动过程时间；计算电器开始启动后l个工频周期之内的电器负载电流有效值的平均值与电器负载稳态电流之间的比值，将该比值作为电器的启动冲激电流；计算电器的启动过程时间之内的电器负载电流有效值的平均值与电器负载稳态电流之间的比值，将该比值作为电器的启动平均电流；计算电器的启动平均电流与启动过程时间之间的乘积，将该乘积作为电器的启动电流冲量；所述l的取值范围为1-5。
所述组合分类器的输入特征还包括电器负载稳态电流。
所述学生宿舍用电器类型分类方法由包括信息采集模块、信息处理模块、通信模块的装置实现。所述信息采集模块用于采集电器的负载电流信息和负载电压信息并送至信息处理模块；所述信息处理模块依据输入的信息进行电器类型识别；所述通信模块用于发送信息处理模块的电器类型识别结果至上位机。
所述信息采集模块包括电流传感器、前置放大器、滤波器、a/d转换器；所述信息处理模块的核心为dsp，或者为arm，或者为单片机，或者为fpga。
所述a/d转换器可以采用信息处理模块的核心中包括的a/d转换器。
所述信息采集模块、信息处理模块、通信模块的全部或者部分功能集成在一片soc(systemonchip，片上系统)上。
所述通信模块还接收上位机的相关工作指令；所述通信模块与上位机之间的通信方式包括无线通信方式与有线通信方式；所述无线通信方式包括zigbee、蓝牙、wifi、433mhz数传方式；所述有线通信方式包括485总线、can总线、互联网、电力载波方式。
本发明的有益效果是：采用包括支持向量机分类器和贝叶斯分类器的组合分类器进行识别分类，兼顾支持向量机分类器和贝叶斯分类器的特点进行综合识别，泛化能力与识别准确率高；同时采用电器的启动电流特征、电器的负载电流频谱特征以及电器的基波电压电流相位差作为所述学生宿舍用电器类型分类方法的识别特征，特征信息丰富；提供的包括启动冲激电流、启动平均电流、启动电流冲量在内的启动电流特征获取方法，以及负载电流频谱特征获取方法简单、可靠。
附图说明
图1为实现学生宿舍用电器类型分类方法的装置实施例结构示意图；
图2为白炽灯台灯的启动过程电流波形；
图3为电阻炉等电阻性负载的启动过程电流波形；
图4为单相电机类负载的启动过程电流波形；
图5为计算机及开关电源类负载的启动过程电流波形；
图6为学生宿舍用电器类型分类方法流程图。
具体实施方式
以下结合附图对本发明作进一步说明。
图1为实现学生宿舍用电器类型分类方法的装置实施例结构示意图，包括信息采集模块101、信息处理模块102、通信模块103。
信息采集模块102用于采集电器的负载电压、负载电流并将负载电压、负载电流转换成电压数字信号、电流数字信号，电压数字信号、电流数字信号被送至信息处理模块102。信息采集模块中包括电压传感器、电流传感器、前置放大器、滤波器、a/d转换器等组成部分，分别完成负载电压、负载电流信号的传感、放大、滤波与模数转换功能。当负载电流范围较大时，可以选择具有程控功能的前置放大器，或者是在a/d转换器前再增加一个独立的程控放大器，对范围较大的负载电流实行分段控制放大，使输入至a/d转换器的电压信号范围保持在合理的区间，保证转换精度。滤波器用于滤除高频分量，避免频谱混叠。
信息处理模块102依据输入的电压数字信号、电流数字信号，采用包括支持向量机分类器和贝叶斯分类器的组合分类器实现电器类型识别。组合分类器的输入特征包括电器的启动电流特征、电器的负载电流频谱特征和电器的基波电压电流相位差。信息处理模块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所示为白炽灯台灯的启动过程电流波形。白炽灯是将灯丝通电加热到白炽状态，利用热辐射发出可见光的电光源。白炽灯的灯丝通常用耐高温的金属钨制造，但金属钨的电阻随温度变化大，以rt表示钨丝在t℃时的电阻，以r0表示钨丝在0℃时的电阻，则两者有下述的关系
rt＝r0(1+0.0045t)
例如，设白炽灯的灯丝(钨丝)在正常工作时的温度为2000℃，一只“220v100w”的白炽灯的灯丝在2000℃正常工作时的电阻为
其在不通电时0℃的电阻为
其在不通电时20℃的电阻为
r20＝r0(1+0.0045t)＝52.8ω
即白炽灯在启动通电的瞬间电流超过其额定电流的9倍，且最大启动电流发生在启动时刻。随着白炽灯钨丝温度的升高，白炽灯的负载电流按照指数规律减小，然后进入稳定状态。
设电器稳态电流有效值为iw，且定义电器电流有效值进入电器稳态电流有效值的一个设定的相对误差范围之内并稳定在这个相对误差范围之内，则电器进入稳定状态。相对误差范围可以设定为10％，也可以设定为2％、5％、15％、20％等2％-20％之间的值。图2中，设定的相对误差范围为10％，当白炽灯的负载电流按照指数规律减小到其iw的10％误差范围时，如图2中的时刻ts，启动过程结束。白炽灯的启动过程时间为ts。iw为有效值。
选择启动冲激电流ig、启动平均电流id、启动电流冲量qi作为电器的启动电流特征；启动冲激电流ig、启动平均电流id均为标么值。具体定义是：启动冲激电流ig为电器启动开始后t2时间之内的电器负载电流平均值与电器负载稳态电流iw的比值；启动平均电流id为电器启动时间ts之内的电器负载电流平均值与电器负载稳态电流iw的比值；启动电流冲量qi为启动平均电流id与启动过程时间ts的乘积，量纲为ms。电器负载电流、电器负载稳态电流均为有效值。t2的取值范围为20-100ms，或者是1-5个工频周期；例如，t2取值40ms，即2个工频周期。启动冲激电流ig反映的是电器负载启动后短时间内的电流冲激大小。在部分电器的启动过程中，当有电器的实际启动过程时间ts小于设定的t2时，令电器的启动过程时间ts等于t2。启动平均电流id反映的是电器负载启动过程中的电流整体大小。启动电流冲量qi反映的是电器负载启动的整体强度。
图2中，白炽灯的启动冲激电流ig为t0(白炽灯启动时刻，电流为i0)至t2(设定的时刻，电流为i2)之间白炽灯的电流平均值与白炽灯的稳态电流iw的比值。启动平均电流id为t0(白炽灯启动时刻)至ts(白炽灯启动过程结束时间)之间白炽灯的电流平均值与白炽灯的稳态电流iw的比值。启动电流冲量qi为白炽灯启动平均电流id与启动过程时间ts的乘积。
如图3所示为电阻炉等电阻性负载的启动过程电流波形。电阻炉等电阻性负载通常采用镍铬、铁铬铝等电热合金丝，其共同特点是电阻温度修正系数小，电阻值稳定。以牌号为cr20ni80的镍铬电热丝为例，其在1000℃时的电阻修正系数为1.014，即1000℃时相对于20℃时，牌号为cr20ni80的镍铬电热丝电阻只增加1.4％。电阻炉等电阻性负载在通电启动时即进入稳定状态，电阻炉等电阻性负载的实际启动过程时间ts＝0，因此，令电阻炉等电阻性负载的实际启动过程时间ts＝t2；例如，当t2设定为40ms时，则此时的启动过程时间ts也为40ms。由于电阻性负载t0时刻电流i0、t2时刻电流i2与电阻性负载的稳态电流iw相等，因此，电阻性负载的启动冲激电流ig＝1，启动平均电流id＝1。
如图4所示为单相电机类负载的启动过程电流波形。单相电机类负载既具有电感性负载特性，又具有反电动势负载特性。启动时刻，由于电感的作用，启动时刻的启动电流i0为0；随后电流迅速上升，在电机反电动势未建立之前，达到电流峰值im；此后，电机转速增加，电机负载电流逐步减小，直到进入稳定状态。图4中，单相电机类负载的启动冲激电流ig为t0(单相电机类负载启动时刻，电流为i0)至t2(设定的时刻，电流为i2)之间单相电机类负载的电流平均值与稳态电流iw的比值。启动平均电流id为t0(单相电机类负载启动时刻)至ts(单相电机类负载启动过程结束时间)之间单相电机类负载的电流平均值与稳态电流iw的比值。启动电流冲量qi为单相电机类负载启动平均电流id与启动过程时间ts的乘积。
如图5所示为计算机及开关电源类负载的启动过程电流波形。计算机及开关电源类负载因为对电容充电的影响，在启动瞬间会产生一个很大的浪涌电流，其峰值可达到稳态电流有效值iw的几倍至十几倍，时间为1至2个工频周期。由于计算机及开关电源类负载的启动时间短，其启动过程时间ts有可能小于设定的t2；当其启动过程时间ts小于设定的t2时，令ts等于t2。图5中，计算机及开关电源类负载的启动冲激电流ig为t0(计算机及开关电源类负载启动时刻，电流为i0)至t2(设定的时刻，电流为i2)之间计算机及开关电源类负载的电流平均值与稳态电流iw的比值。启动平均电流id为t0(计算机及开关电源类负载启动时刻)至ts(计算机及开关电源类负载启动过程结束时间)之间计算机及开关电源类负载的电流平均值与稳态电流iw的比值。启动电流冲量qi为计算机及开关电源类负载启动平均电流id与启动过程时间ts的乘积。
获取电器的启动电流特征的方法是：
电器启动前，负载电流值为0(未开机)或者很小(处于待机状态)时，信息处理模块102即开始对负载电流进行连续采样；当采样得到的负载电流值有效值开始大于0或者是开始大于电器的待机电流时，即判断出电器已经启动，记录该时刻为t0。用一个较小的非负阈值ε来区分电器启动前后的负载电流值，当ε取值特别小时，例如，ε取值1ma时，所述识别装置不考虑待机情况，即认为待机也是电器的启动状态；当ε取值较小但大于电器的待机电流时，例如，ε取值20ma时，所述识别装置会将电器的待机状态认为是未启动状态，但同时也会的部分功率特别小的电器造成漏识别。
信息处理模块102对负载电流进行连续采样，且以工频周期为单位计算负载电流有效值并保存；当电器已经启动，且连续采样达到n个工频周期后，采样的同时连续计算最近n个工频周期的负载电流有效值的平均值iv；信息处理模块102对最近n个工频周期之内每个工频周期的负载电流有效值与该n个工频周期的负载电流有效值的平均值进行比较，误差(或波动)幅度均小于设定的相对误差范围e时，判定电器进入稳定状态，该最近n个工频周期的起始时刻为启动过程的结束时刻，记录该时刻为t1(如图2-图5所示)。
将最近n个工频周期之内的负载电流有效值的平均值作为电器负载稳态电流iw；将电器开始启动时刻t0至最近n个工频周期起始时刻t1之间的时间作为启动过程时间ts。计算t0至设定的t2之间(即电器开始启动后1-5个工频周期之内)的负载电流平均值与稳态电流iw的比值，将该比值作为电器的启动冲激电流ig。计算t0至ts之间的负载电流平均值与稳态电流iw的比值，将该比值作为电器的启动平均电流id。计算电器的启动平均电流id与启动过程时间ts的乘积，将该乘积作为电器的启动电流冲量qi。
由于预先不知道电器稳态电流有效值iw，因此，将n个工频周期，即一段持续时间tp之内波动范围小于设定的相对误差范围e时的负载电流有效值的平均值作为电器稳态电流有效值iw。由于普通电器的启动过程较快，所以，tp的取值范围为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。n＝1时的1次谐波为工频基波。所述谐波信号相对幅值为谐波信号幅值与电器稳态电流有效值iw的比值。所述m表示谐波最高次数，一般情况下，m大于等于3。
组合分类器的输入特征还包括电器的基波电压电流相位差。基波电压电流相位差可以对电阻性、电容性、电感性负载进行区分，还可以对一般的电感性负载和大电感性负载进行区分。电器的基波电压电流相位差由信息处理模块102控制信息采集模块101，通过以下步骤获得：
步骤①、待电器进入稳定状态后，同步获取电器的稳态电压信号、稳态电流信号，并将其转换为对应的稳态电压数字信号、稳态电流数字信号。
步骤②、对稳态电压数字信号、稳态电流数字信号分别进行数字滤波，提取出基波电压信号、基波电流信号。
步骤③、分析计算基波电压信号与基波电流信号之间的相位差，将基波电压信号与基波电流信号之间的相位差作为电器的基波电压电流相位差。
步骤②中对稳态电压数字信号、稳态电流数字信号分别进行数字滤波，其数字滤波算法可以选择卡尔曼滤波法、小波变换法、维纳滤波法、自适应滤波等数字滤波器算法。
组合分类器中，支持向量机分类器为主分类器，贝叶斯分类器为辅助分类器。组合分类器的输入特征包括前述的启动电流特征和负载电流频谱特征，组合分类器的输入特征同时作为支持向量机分类器的输入特征和贝叶斯分类器的输入特征。
如图6所示为学生宿舍用电器类型分类方法流程图，方法流程包括：
步骤a、等待电器启动；
步骤b、采集电器启动电流数据并保存，直至电器启动过程结束；
步骤c、分析采集的电器启动电流数据，获取电器的启动电流特征；
步骤d、采集电器稳态工作时的电压、电流数据并保存；
步骤e、分析采集的电器稳态工作时的电压、电流数据，获取电器的负载电流频谱特征、基波电压电流相位差；
步骤f、将启动电流特征、负载电流频谱特征、基波电压电流相位差作为组合分类器的输入特征；组合分类器进行电器类型识别；
步骤g、输出电器类型识别结果。
所述组合分类器进行电器类型识别的方法是：当主分类器成功实现电器类型识别，即主分类器输出的识别结果为唯一的电器类型，即识别结果中唯一的电器类型为是时，将主分类器识别的电器类型作为组合分类器的电器类型识别结果；当主分类器未能实现电器类型识别，且主分类器的识别结果为2种或者2种以上电器类型，即识别结果中有2种或者2种以上电器类型为是时，将主分类器输出的2种或者2种以上电器类型识别结果中，辅助分类器输出中概率最高的电器类型作为组合分类器的电器类型识别结果；当主分类器未能实现电器类型识别，且主分类器的识别结果中未能给出识别的电器类型，即识别结果中没有电器类型为是时，将辅助分类器输出中概率最高的电器类型作为组合分类器的电器类型识别结果。
以一个简单的实施例1为例，来说明组合分类器进行电器类型识别的方法。设有一个组合分类器，其输入特征为
在实施例1中，b1、b2、b3、b4的分类标记和f1、f2、f3、f4的分类标记均取1、0。分类标记为1时，相应的电器类型与当前输入特征匹配，为确认的识别结果，或者说相应的电器类型识别结果为是；分类标记为0时，相应的电器类型与输入特征不匹配，未能成为确认的识别结果，或者说相应的电器类型识别结果为否。
在实施例1中，设某次的主分类器的识别分类标记为f1f2f3f4＝0100，则认为主分类器成功实现电器类型识别，因此，不考虑辅助分类器的识别结果，直接令b1b2b3b4＝0100，即组合分类器的识别结果是：被识别的电器为电阻炉。
在实施例1中，设某次的主分类器的识别分类标记为f1f2f3f4＝1010，则认为主分类器未能实现电器类型识别，且主分类器的识别结果为2种或者2种以上电器类型；再设此时辅助分类器的识别结果满足p(y1|x)<p(y3|x)，则令b1b2b3b4＝0010，即组合分类器的识别结果是：被识别的电器为吹风机。
在实施例1中，设某次的主分类器的识别分类标记为f1f2f3f4＝0000，则认为主分类器未能实现电器类型识别，且主分类器的识别结果中未能给出识别的电器类型；再设此时辅助分类器的识别结果满足p(y1|x)>p(y2|x)且p(y1|x)>p(y3|x)且p(y1|x)>p(y4|x)，则令b1b2b3b4＝1000，即组合分类器的识别结果是：被识别的电器为白炽灯。
组合分类器、主分类器的识别结果分类标记也可以采用其他的方案，例如，分别用分类标记1、-1，或者是0、1，或者是-1、1，以及其他方案来表示相应电器识别结果为是、否。组合分类器与主分类器的分类标记方案可以相同，也可以不相同。
所述组合分类器的输入特征中，还可以包括电器稳态电流有效值iw。例如，有2种不同的电器，电烙铁和电阻炉需要识别，电烙铁、电阻炉都是纯电阻负载，且都具有电阻温度修正系数小，电阻值稳定的共同特点。因此，单纯依靠前述的启动电流特征和负载电流频谱特征、电器的基波电压电流相位差特征无法将他们进行区分。输入特征中增加电器稳态电流有效值iw后，电烙铁功率小，电器稳态电流有效值iw小；电阻炉功率大，电器稳态电流有效值iw大，特征不同，组合分类器可以进行并完成识别。
辅助分类器为贝叶斯分类器。可以选择nbc分类器(朴素贝叶斯分类器)、tan分类器(树扩展朴素贝叶斯分类器)、ban分类器(增强的贝叶斯分类器)等三种贝叶斯分类器之中的一种作为辅助分类器。
实施例2选择nbc分类器作为辅助分类器。朴素贝叶斯分类的定义如下：
⑴设x＝{a1，a2，…，am}为一个待分类项，而每个a为x的一个特征属性；
⑵有类别集合c＝{y1，y2，…，yn}；
⑶计算p(y1|x)，p(y2|x)，…，p(yn|x)；
⑷如果p(yk|k)＝max{p(y1|x)，p(y2|x)，…，p(yn|x)}，则x∈yk。
计算第⑶步中的各个条件概率的具体方法是：
①找到一个已知分类的待分类项集合作为训练样本集；
②统计得到各类别下各个特征属性的条件概率估计；
p(a1|y1)，p(a2|y1)，…，p(am|y1)；
p(a1|y2)，p(a2|y2)，…，p(am|y2)；
…；
p(a1|yn)，p(a2|yn)，…，p(am|yn)。
③根据贝叶斯定理，有：
因为分母对于所有类别为常数，因此我们只要将分子最大化即可；又因为在朴素贝叶斯中各特征属性是条件独立的，所以有：
实施例2中，组合分类器的输入特征是
训练nbc分类器的过程包括：
1、对特征属性进行分段划分，进行离散化处理。实施例2中，采取的特征属性离散化方法是：
a1：{a1＜3.5，3.5≤a1≤7，a1＞7}；
a2：{a2＜1.2，1.2≤a2≤1.8，a2＞1.8}；
a3：{a3＜120，120≤a3≤560，a3＞560}；
a4：{a4＜0.7，0.7≤a4≤0.9，a4＞0.9}；
a5：{a5＜0.02，0.02≤a5≤0.05，a5＞0.05}；
a6：{a6＜-5，-5≤a6≤12，a6＞12}。
2、对每类电器类型均采集多组样本作为训练样本，同时计算每类电器类型样本在所有电器类型样本中所占有的比例，即分别计算p(y1)、p(y2)、p(y3)、p(y4)、p(y5)。当每类电器均采集相同的样本数量时，例如，每类电器均采集超过100组的样本，其中每类电器随机选择100组样本作为训练样本，其他则作为测试样本，总的训练样本为500组，且有
p(y1)＝p(y2)＝p(y3)＝p(y4)＝p(y5)＝0.2。
3、计算训练样本每个类别条件下各个特征属性分段的频率(比例)，统计得到各类别下各个特征属性的条件概率估计，即分别统计计算
p(a1＜3.5|y1)、p(3.5≤a1≤7|y1)、p(a1＞7|y1)；
p(a1＜3.5|y2)、p(3.5≤a1≤7|y2)、p(a1＞7|y2)；
…；
p(a1＜3.5|y5)、p(3.5≤a1≤7|y5)、p(a1＞7|y5)；
p(a2＜1.2|y1)、p(1.2≤a2≤1.8|y1)、p(a2＞1.8|y1)；
p(a2＜1.2|y2)、p(1.2≤a2≤1.8|y2)、p(a2＞1.8|y2)；
…；
p(a2＜1.2|y5)、p(1.2≤a2≤1.8|y5)、p(a2＜1.8|y5)；
p(a3＜120|y1)、p(120≤a3≤560|y1)、p(a3＞560|y1)；
p(a2＜120|y2)、p(120≤a3≤560|y2)、p(a3＞560|y2)；
…；
p(a3＜120|y5)、p(120≤a3≤560|y5)、p(a3＞560|y5)；
p(a4＜0.7|y1)、p(0.7≤a4≤0.9|y1)、p(a4＞0.9|y1)；
p(a4＜0.7|y2)、p(0.7≤a4≤0.9|y2)、p(a4＞0.9|y2)；
…；
p(a4＜0.7|y5)、p(0.7≤a4≤0.9|y5)、p(a4＞009|y5)；
p(a5＜0.02|y1)、p(0.02≤a5≤0.05|y1)、p(a5＞0.05|y1)；
p(a5＜0.02|y2)、p(0.02≤a5≤0.05|y2)、p(a5＞0.05|y2)；
p(a5＜0.02|y5)、p(0.02≤a5≤0.05|y5)、p(a5＞0.05|y5)；
p(a6＜-5|y1)、p(-5≤a6≤12|y1)、p(a6≥12|y1)；
p(a6＜-5|y2)、p(-5≤a6≤12|y2)、p(a6≥12|y2)；
…；
p(a6＜-5|y5)、p(-5≤a6≤12|y5)、p(a6≥12|y5)。
经过上述的步骤1、步骤2、步骤3，nbc分类器训练完成。其中，步骤1对特征属性进行分段划分由人工确定，对每一个输入特征进行分段离散化时，分段的数量为2段或者2段以上，例如，实施例2中，特征a1-a6都分为3段。每一个特征具体分为多少段，分段阈值的选择可以根据训练后的贝叶斯分类器对测试样本测试后的结果进行调整。步骤2、步骤3由信息处理模块102或者是计算机计算完成。
本发明中采用贝叶斯分类器进行分类的方法是：
1、将组合分类器的输入特征作为贝叶斯分类器的输入特征。在实施例2中，将组合分类器的输入特征集合
2、根据训练得到的各类别下各个特征属性的条件概率估计，分别确定各输入特征属性的分段所在并确定其对每类电器类别的概率p(a1|y1)～p(am|yn)，其中，电器类别集合为c＝{y1，y2，…，yn}。实施例2中，电器类别集合c＝{y1，y2，y3，y4，y5}对应代表的电器类别是白炽灯、电阻炉、电风扇、计算机、吹风机，确定p(a1|y1)～p(a6|y5)的方法是采用训练nbc分类器过程中得到的各个特征属性的条件概率估计。
3、按照式
计算每种电器类别的后验概率。因为分母p(x)对于所有电器类别为常数，令p(x)＝1替代实际的p(x)值，不影响每种电器类别后验概率之间的相互大小比较，此时有
实施例2中，有
采用测试样本对训练好的贝叶斯分类器进行测试，根据测试结果决定是否调整对输入特征的离散化方法(即调整分段数量与阈值)，重新训练贝叶斯分类器。
主分类器为支持向量机分类器，或称为svm分类器。svm分类器特别适合解决二值分类情况，因此，主分类器采用多个二类输出svm分类器组成，每个二类输出svm分类器对应识别一种电器类型，例如，实施例1中可以采用4个二类输出svm分类器分别识别白炽灯、电阻炉、吹风机、计算机，实施例2中可以采用5个二类输出svm分类器分别识别白炽灯、电阻炉、电风扇、计算机、吹风机。主分类器选择多个二类输出svm分类器共同组成时，所有二类输出svm分类器的输入特征均为主分类器的输入特征。
训练各二类输出svm分类器时，对每类电器类型均采集多组样本，随机抽取部分作为训练样本，其余的作为测试样本。样本采集采用前述的获取电器的启动电流特征的方法和获取电器的负载电流频谱特征以及获取电器的基波电压电流相位差特征的方法。所有的训练样本均作为每个二类输出svm分类器的训练样本。例如，在实施例2中，可以分别对白炽灯、电阻炉、电风扇、计算机、吹风机等负载均采集超过100组样本，随机抽取其中每种100组、共500组样本组成训练样本，剩余的样本组成测试样本；当然，某种负载或者所有负载采集的样本数量达不到100组样本，svm分类器也能够取得较好的分类效果。
主分类器所选用的二类输出svm分类器选择径向基rbf核函数，并采用粒子群算法(pso)对各二类输出svm分类器的惩罚参数c和核参数y进行优化选择。
每个二类输出svm分类器只需要完成一种电器类型的分类识别，svm分类器的训练相对简单。主分类器由多个二类输出svm分类器组成，各二类输出svm分类器之间相互独立，因此，对某一特征属性进行识别时，主分类器有可能输出的识别结果为唯一的电器类型，或者识别结果为2种或者2种以上电器类型，或者未能给出分类识别结果的电器类型。