电器负载类型智能识别方法与装置的制造方法

电器负载类型智能识别方法与装置的制造方法
【专利摘要】一种电器负载类型智能识别装置，包括信息采集模块、信息处理模块、通信模块。所述装置同时采用电器负载的启动电流特征、电器负载的基波电压电流相位差和负载电流频谱特征作为电器负载的识别特征，特征信息丰富；采用包括支持向量机分类器和贝叶斯分类器的组合分类器进行识别分类，兼顾两种分类器的特点进行综合识别，识别准确率高；提供的基波电压电流相位差、启动电流特征、负载电流频谱特征获取方法简单、可靠。所述装置可以用在学生集体宿舍、大型集贸市场等一些需要进行电器负载管理的集体公共场所，也可以用于需要进行电器负载类型与统计的其他需要进行用电设备管理的场合。
【专利说明】
电器负载类型智能识别方法与装置
技术领域
[0001] 本发明涉及一种设备识别及分类装置，尤其是涉及一种电器负载类型智能识别方 法与装置。
【背景技术】
[0002] 目前，主流的电器负载性质或者电器类型识别方法包括基于负载功率综合系数算 法的电器负载识别方法、基于电磁感应的电器负载识别方法、基于神经网络算法的电器负 载识别方法、基于周期性离散变换算法的电器负载识别方法等。各种方法均能够在一定程 度是实现电器负载性质的识别，但由于特征性质单一，识别手段单一，普遍存在泛化能力不 够及不能完全准确识别的问题。

【发明内容】

[0003] 本发明的目的在于，针对现在已有技术的缺陷，提供一种能够实现高效识别的电 器负载类型智能识别方法与装置。所述装置包括信息采集模块、信息处理模块、通信模块。
[0004] 所述信息采集模块用于采集电器负载的负载电流并转换成电流数字信号;所述电 流数字信号被送至信息处理模块;所述信息处理模块依据输入的电流数字信号，采用组合 分类器电器负载类型智能识别;所述通信模块用于发送信息处理模块的电器负载类型识别 结果至上位机。
[0005] 所述组合分类器的输入特征包括电器负载的启动电流特征、电器负载的负载电流 频谱特征和电器负载的基波电压电流相位差;所述组合分类器包括支持向量机分类器和贝 叶斯分类器;所述启动电流特征包括启动过程时间、启动电流最大值、启动电流最大值时 间。
[0006] 所述信息采集模块包括电流传感器、前置放大器、滤波器、A/D转换器;所述信息处 理模块的核心为DSP，或者为ARM，或者为单片机，或者为FPGA。
[0007] 所述A/D转换器可以采用信息处理模块的核心中包括的A/D转换器。
[0008] 所述信息采集模块、信息处理模块、通信模块的全部或者部分功能集成在一片SoC 上。
[0009] 所述通信模块还接收上位机的相关工作指令;所述通信模块与上位机之间的通信 方式包括无线通信方式与有线通信方式；所述无线通信方式包括ZigBee、蓝牙、WiFi、 433MHz数传方式;所述有线通信方式包括485总线、CAN总线、互联网、电力载波方式。
[0010] 所述负载电流频谱特征通过以下方法获得：
[0011]步骤一、获取电器负载的稳态电流信号，并将其转换为对应的稳态电流数字信号；
[0012] 步骤二、对稳态电流数字信号进行傅立叶变换，得到负载电流频谱特性；
[0013] 步骤三、将负载电流频谱特性中谐波次数为n次的奇次谐波信号相对幅值作为负 载电流频谱特征，n = 1，3，…，M;所述M表示谐波最高次数且M大于等于3。
[0014]所述组合分类器中，支持向量机分类器为主分类器，贝叶斯分类器为辅助分类器。
[0015]所述组合分类器进行电器负载类型智能识别的方法是：当主分类器成功实现电器 负载类型识别时，主分类器的电器负载类型识别结果为组合分类器的识别结果；当主分类 器未能实现电器负载类型识别，且主分类器的识别结果为2种或者2种以上电器负载类型， 将主分类器输出的2种或者2种以上电器负载类型识别结果中，辅助分类器输出中概率最高 的电器负载类型作为组合分类器的电器负载类型识别结果；当主分类器未能实现电器负载 类型识别，且主分类器的识别结果中未能给出识别的电器负载类型时，将辅助分类器输出 中概率最高的电器负载类型作为组合分类器的电器负载类型识别结果。
[0016]所述启动电流特征由信息处理模块通过以下方法获得：
[0017] 步骤1、电器负载启动前，开始对电器负载的负载电流连续采样并对负载电流大小 进行判断；当负载电流有效值大于e时，判定电器负载开始启动并转向步骤2;所述e为大于0 的数值；
[0018] 步骤2、对电器负载的负载电流进行连续采样，以工频周期为单位计算负载电流有 效值并保存;计算最近N个工频周期的负载电流有效值的平均值；当最近N个工频周期之内 的每个工频周期的负载电流有效值与该N个工频周期的负载电流有效值的平均值相比较， 波动幅度均小于设定的相对误差范围E时，判定电器负载进入稳定状态，转向步骤3;所述N 的取值范围为50-500;所述E的取值范围为2%-20% ;
[0019]步骤3、将最近N个工频周期之内的负载电流有效值的平均值作为电器负载稳态电 流有效值;将电器负载开始启动时刻至最近N个工频周期起始时刻之间的时间作为启动过 程时间;将电器负载开始启动时刻至启动过程时间之内负载电流有效值最大的工频周期之 间的时间作为启动电流最大值时间;将启动电流最大值时间所在工频周期的负载电流有效 值与电器负载稳态电流有效值之间的比值作为启动电流最大值。
[0020] 所述组合分类器的输入特征还包括电器负载稳态电流有效值。
[0021] 所述电器负载的基波电压电流相位差通过以下方法获得：
[0022] 步骤①、待电器负载进入稳定状态后，同步获取电器负载的稳态电压信号、稳态电 流信号，并将其转换为相应的稳态电压数字信号、稳态电流数字信号；
[0023] 步骤②、对稳态电压数字信号、稳态电流数字信号分别进行数字滤波，提取出基波 电压信号、基波电流信号；
[0024] 步骤③、分析计算基波电压信号与基波电流信号之间的相位差，将基波电压信号 与基波电流信号之间的相位差作为电器负载的基波电压电流相位差。
[0025] 本发明的有益效果是：同时采用电器负载的启动电流特征、电器负载的负载电流 频谱特征以及电器负载的基波电压电流相位差作为所述电器负载类型智能识别装置的识 别特征，特征信息丰富；采用包括支持向量机分类器和贝叶斯分类器的组合分类器进行识 别分类，兼顾支持向量机分类器和贝叶斯分类器的特点进行综合识别，泛化能力与识别准 确率高;提供的包括启动过程时间、启动电流最大值、启动电流最大值时间在内的启动电流 特征获取方法，以及负载电流频谱特征获取方法简单、可靠。
【附图说明】
[0026] 图1为本发明电器负载识别装置实施例的结构示意图；
[0027] 图2为白炽灯台灯的启动过程电流波形；
[0028] 图3为电阻炉等电阻性负载的启动过程电流波形；
[0029] 图4为单相电机类负载的启动过程电流波形；
[0030] 图5为计算机及开关电源类负载的启动过程电流波形；
[0031 ]图6为电器负载类型智能识别方法流程图。
【具体实施方式】
[0032] 以下结合附图对本发明作进一步说明。
[0033] 图1为本发明电器负载识别装置实施例的结构示意图，包括信息采集模块101、信 息处理模块102、通信模块103。
[0034] 信息采集模块102用于采集电器负载的负载电压、负载电流并将负载电压、负载电 流转换成电压数字信号、电流数字信号，电压数字信号、电流数字信号被送至信息处理模块 102。信息采集模块中包括电压传感器、电流传感器、前置放大器、滤波器、A/D转换器等组成 部分，分别完成负载电压、负载电流信号的传感、放大、滤波与模数转换功能。当负载电流范 围较大时，可以选择具有程控功能的前置放大器，或者是在A/D转换器前再增加一个独立的 程控放大器，对范围较大的负载电流实行分段控制放大，使输入至A/D转换器的电压信号范 围保持在合理的区间，保证转换精度。滤波器用于滤除高频分量，避免频谱混叠。
[0035]信息处理模块102依据输入的电压数字信号、电流数字信号，采用包括支持向量机 分类器和贝叶斯分类器的组合分类器实现电器负载类型识别。组合分类器的输入特征包括 电器负载的启动电流特征、电器负载的负载电流频谱特征和电器负载的基波电压电流相位 差。信息处理模块10 2的核心为DSP、ARM、单片机，或者为FPGA。当信息处理模块的核心中包 括有A/D转换器且该A/D转换器满足要求时，信息采集模块101中的A/D转换器可以采用信息 处理模块102的核心中包括的A/D转换器。
[0036]通信模块103用于实现与上位机之间的通信，将识别结果发送至上位机。通信模块 102与上位机之间的通信方式包括无线通信方式与有线通信方式，可以采用的无线通信方 式包括ZigBee、蓝牙、WiFi、433MHz数传等方式，可以采用的有线通信方式包括485总线、CAN 总线、互联网、电力载波等方式。通信模块103还可以接收上位机的相关工作指令，完成指定 的工作任务。上位机可以是管理部门的服务器，也可以是各种工作站，或者是各种移动终 端。
[0037]信息采集模块101、信息处理模块102、通信模块103的全部或者部分功能可以集成 在一片SoC上，减小装置体积，方便安装。
[0038] 不同的电器负载设备具有不同的启动电流特征。如图2所示为白炽灯台灯的启动 过程电流波形。白炽灯是将灯丝通电加热到白炽状态，利用热辐射发出可见光的电光源。白 炽灯的灯丝通常用耐高温的金属钨制造，但金属钨的电阻随温度变化大，以Rt表示钨丝在t °(：时的电阻，以Ro表示钨丝在0°C时的电阻，则两者有下述的关系
[0039] Rt = R〇(l+0.0045t)
[0040] 例如，设白炽灯的灯丝(钨丝)在正常工作时的温度为2000°C，一只"220V 100W"的 白炽灯的灯丝在2000 °C正常工作时的电阻为
[0042] 其在不通电时0°C的电阻为
[0044] 其在不通电时20 °C的电阻为
[0045] R2〇 = Ro( 1+0.0045t) =52.8 Q
[0046] 即白炽灯在启动通电的瞬间电流超过其额定电流的9倍，且最大启动电流发生在 启动时刻。随着白炽灯钨丝温度的升高，白炽灯的负载电流按照指数规律减小，然后进入稳 定状态。
[0047] 设电器负载稳态电流有效值为IW，且定义电器负载电流有效值进入电器负载稳态 电流有效值的一个设定的相对误差范围之内并稳定在这个相对误差范围之内，则电器负载 进入稳定状态。相对误差范围可以设定为10%，也可以设定为2%、5%、15%、20%等2%-20%之间的值。图2中，设定的相对误差范围为10%，当白炽灯的负载电流按照指数规律减 小到其Iw的10 %误差范围时，如图2中的时刻Ts，启动过程结束。白炽灯的启动过程时间为 Ts。Iw为有效值。
[0048] 选择启动过程时间、启动电流最大值I*、启动电流最大值时间作为电器负载的启 动电流特征；启动电流最大值为标么值，即启动电流最大值I*为启动电流的最大有效值Im 与电器负载稳态电流有效值Iw的比值。
[0049] 图2中，白炽灯的启动过程时间为Ts;启动电流最大值I*为Im/Iw，其值约在9-10之 间；启动电流最大值时间为Tm，Tm=0。
[0050] 如图3所示为电阻炉等电阻性负载的启动过程电流波形。电阻炉等电阻性负载通 常采用镍铬、铁铬铝等电热合金丝，其共同特点是电阻温度修正系数小，电阻值稳定。以牌 号为Cr20Ni80的镍铬电热丝为例，其在1000°C时的电阻修正系数为1.014，即1000°(：时相对 于20 °C时，牌号为Cr20Ni80的镍铬电热丝电阻只增加1.4 %。因此，电阻炉等电阻性负载在 通电启动时即进入稳定状态，电阻炉等电阻性负载的启动过程时间Ts = 0;启动电流最大值 I* = 1;启动电流最大值时间Tm=0。
[0051] 如图4所示为单相电机类负载的启动过程电流波形。单相电机类负载既具有电感 性负载特性，又具有反电动势负载特性。启动时刻，由于电感的作用，启动时刻的启动电流 为〇;随后电流迅速上升，在电机反电动势未建立之前，达到电流峰值Im;此后，电机转速增 加，电机负载电流逐步减小，直到进入稳定状态。图4中，单相电机类负载的启动过程时间为 Ts;启动电流最大值I*为Im/Iw;启动电流最大值时间为Tm。
[0052] 如图5所示为计算机及开关电源类负载的启动过程电流波形。计算机及开关电源 类负载因为对电容充电的影响，在启动瞬间会产生一个很大的浪涌电流，其峰值可达到稳 态电流有效值Iw的几倍至十几倍，时间为1至2个工频周期。图5中，计算机及开关电源类负 载的启动过程时间为Ts，约1至2个工频周期；启动电流最大值I*为Im/Iw;启动电流最大值时 间为Tm=0。
[0053]获取电器负载的启动电流特征的方法是：
[0054]电器负载启动前，负载电流值为0(未开机)或者很小(处于待机状态)时，信息处理 模块102即开始对负载电流进行连续采样；当采样得到的负载电流值有效值开始大于0或者 是开始大于电器负载的待机电流时，即判断出电器负载已经启动，记录该时刻为To。用一个 较小的非负阈值e来区分电器负载启动前后的负载电流值，当e取值特别小时，例如，e取值 1mA时，所述识别装置不考虑待机情况，即认为待机也是电器负载的启动状态；当￡取值较小 但大于电器负载的待机电流时，例如，e取值20mA时，所述识别装置会将电器负载的待机状 态认为是未启动状态，但同时也会的部分功率特别小的电器负载造成漏识别。
[0055] 信息处理模块102对负载电流进行连续采样，且以工频周期为单位计算负载电流 有效值并保存；当电器负载已经启动，且连续采样达到N个工频周期后，采样的同时连续计 算最近N个工频周期的负载电流有效值的平均值Iv;信息处理模块102对最近N个工频周期 之内每个工频周期的负载电流有效值与该N个工频周期的负载电流有效值的平均值进行比 较，误差(或波动)幅度均小于设定的相对误差范围E时，判定电器负载进入稳定状态，该最 近N个工频周期的起始时刻为启动过程的结束时刻，记录该时刻为h。
[0056] 将最近N个工频周期之内的负载电流有效值的平均值作为电器负载稳态电流有效 值If;将电器负载开始启动时刻To至最近N个工频周期起始时刻h之间的时间作为启动过程 时间Ts;将~至!^之内负载电流有效值最大的工频周期所在时刻记录为T 2,将^至!^之间的 时间作为启动电流最大值时间Tm ;将1~2所在工频周期的负载电流有效值与电器负载稳态电 流有效值之间的比值作为启动电流最大值I*。
[0057]由于预先不知道电器负载稳态电流有效值Iw，因此，将N个工频周期，即一段持续 时间TP之内波动范围小于设定的相对误差范围E时的负载电流有效值的平均值作为电器负 载稳态电流有效值Iw。由于普通电器负载的启动过程较快，所以，T P的取值范围为1 -10s，典 型取值是2s，相应的工频周期数量N的取值范围为50-500，N的典型取值是100。所述相对误 差范围E的取值范围为2%-20%，E的典型取值是10%。
[0058] 组合分类器的输入特征还包括电器负载的负载电流频谱特征。电器负载的负载电 流频谱特征由信息处理模块102控制信息采集模块101，通过以下步骤获得：
[0059] 步骤一、待电器负载进入稳定状态后，获取电器负载的稳态电流信号，并将其转换 为对应的稳态电流数字信号。
[0060] 步骤二、对稳态电流数字信号进行傅立叶变换，得到负载电流频谱特性。为保证傅 立叶变换的顺利进行，在前述获取电器负载的稳态电流信号，并将其转换为对应的稳态电 流数字信号的过程中，A/D转换器的精度和速度需要满足傅立叶变换的要求，采样频率可以 设定为10kHz，或者是其他数值;信息处理模块102对采集到的稳态电流数字信号进行FFT运 算，计算其频谱。
[0061 ]步骤三、将负载电流频谱特性中的n次谐波信号相对幅值作为负载电流频谱特征， 其中，n= 1，2，…，M;在组成组合分类器的输入特征向量时，n次谐波信号相对幅值在输入特 征向量中按照1，2,…，M的顺序依次排列。由于负载电流频谱特性主要由奇次谐波组成，除 少数电器负载设备外，偶次谐波分量几乎为0,因此，也可以将负载电流频谱特性中谐波次 数为n次的奇次谐波信号相对幅值依序作为负载电流频谱特征，其中，n=l，3,…，时 的1次谐波为工频基波。所述谐波信号相对幅值为谐波信号幅值与电器负载稳态电流有效 值Iw的比值。所述M表示谐波最高次数，一般情况下，M大于等于3。
[0062]组合分类器的输入特征还包括电器负载的基波电压电流相位差。基波电压电流相 位差可以对电阻性、电容性、电感性负载进行区分，还可以对一般的电感性负载和大电感性 负载进行区分。电器负载的基波电压电流相位差由信息处理模块102控制信息采集模块 101，通过以下步骤获得：
[0063] 步骤①、待电器负载进入稳定状态后，同步获取电器负载的稳态电压信号、稳态电 流信号，并将其转换为对应的稳态电压数字信号、稳态电流数字信号。
[0064] 步骤②、对稳态电压数字信号、稳态电流数字信号分别进行数字滤波，提取出基波 电压信号、基波电流信号。
[0065] 步骤③、分析计算基波电压信号与基波电流信号之间的相位差，将基波电压信号 与基波电流信号之间的相位差作为电器负载的基波电压电流相位差。
[0066] 步骤②中对稳态电压数字信号、稳态电流数字信号分别进行数字滤波，其数字滤 波算法可以选择卡尔曼滤波法、小波变换法、维纳滤波法、自适应滤波等数字滤波器算法。
[0067] 组合分类器中，支持向量机分类器为主分类器，贝叶斯分类器为辅助分类器。组合 分类器的输入特征包括前述的启动电流特征和负载电流频谱特征，组合分类器的输入特征 同时作为支持向量机分类器的输入特征和贝叶斯分类器的输入特征。
[0068] 如图6所示为电器负载类型智能识别方法流程图，具体步骤包括：
[0069]步骤A、等待电器负载启动；
[0070] 步骤B、采集电器负载启动电流数据并保存，直至电器负载启动过程结束；
[0071] 步骤C、分析采集的电器负载启动电流数据，获取电器负载的启动电流特征；
[0072] 步骤D、采集电器负载稳态工作时的电压、电流数据并保存；
[0073] 步骤E、分析采集的电器负载稳态工作时的电压、电流数据，获取电器负载的负载 电流频谱特征、基波电压电流相位差；
[0074] 步骤F、将启动电流特征、负载电流频谱特征、基波电压电流相位差作为组合分类 器的输入特征;组合分类器进行电器负载类型识别；
[0075]步骤G、输出电器负载类型识别结果。
[0076] 所述组合分类器进行电器负载类型识别的方法是：当主分类器成功实现电器负载 类型识别，即主分类器输出的识别结果为唯一的电器负载类型，即识别结果中唯一的电器 负载类型为是时，将主分类器识别的电器负载类型作为组合分类器的电器负载类型识别结 果；当主分类器未能实现电器负载类型识别，且主分类器的识别结果为2种或者2种以上电 器负载类型，即识别结果中有2种或者2种以上电器负载类型为是时，将主分类器输出的2种 或者2种以上电器负载类型识别结果中，辅助分类器输出中概率最高的电器负载类型作为 组合分类器的电器负载类型识别结果；当主分类器未能实现电器负载类型识别，且主分类 器的识别结果中未能给出识别的电器负载类型，即识别结果中没有电器负载类型为是时， 将辅助分类器输出中概率最高的电器负载类型作为组合分类器的电器负载类型识别结果。
[0077] 以一个简单的实施例1为例，来说明组合分类器进行电器负载类型识别的方法。设 有一个组合分类器，其输入特征为％ = {2V .1*，7^:，南，爲> 南，▲；，為.，_於}，其中，Ts 是启动过程时间，单位是ms; I*是启动电流最大值;Tm是启动电流最大值时间，单位是ms ;Ai、 知、知^4、心为负载电流频谱特性中的1-5次谐波信号相对幅值，@为电器负载的基波电压电 流相位差。组合分类器的输出是{Bl，B2，B3，B4}，Bl、B2、B3、B4分别代表组合分类器对白炽灯、 电阻炉、吹风机、计算机的识别结果输出，识别结果8 1、82、8334的取值均为二值分类标记。 主分类器的输入特征也是x = U's，P，7\卜也，七，木，南，J5>，，其输出是{FhFs， 分别代表主分类器对白炽灯、电阻炉、吹风机、计算机的识别结果输出， 识别结果、F4的取值也均为二值分类标记。辅助分类器的输入特征同样为 x = {rs, /*, Ik,為，士，為，j4,為，免}，其输出是{P(yi|x)，P(y2|x)，P(y 3|x)，P(y4 x)}，P(yi | x)、P(y21 x)、P(y31 x)、P(y41 x)为辅助分类器输出的后验概率，P(yi | x)、P(y21 x)、P (y31 x)、P(y41 x)之间的相互大小表明辅助分类器的当前输入特征表示所识别的电器负载属 于白炽灯、电阻炉、吹风机、计算机的可能性大小。
[0078] 在实施例1中，Bl、B2、B3、B4的分类标记和Fl、F2、F3、F4的分类标记均取1、0。分类标记 为1时，相应的电器负载类型与当前输入特征匹配，为确认的识别结果，或者说相应的电器 负载类型识别结果为是;分类标记为0时，相应的电器负载类型与输入特征不匹配，未能成 为确认的识别结果，或者说相应的电器负载类型识别结果为否。
[0079] 在实施例1中，设某次的主分类器的识别分类标记为FJ2F3F4 = 0100，则认为主分 类器成功实现电器负载类型识别，因此，不考虑辅助分类器的识别结果，直接令 0100,即组合分类器的识别结果是:被识别的电器负载为电阻炉。
[0080] 在实施例1中，设某次的主分类器的识别分类标记为FW2F3F4 = 1010，则认为主分 类器未能实现电器负载类型识别，且主分类器的识别结果为2种或者2种以上电器负载类 型;再设此时辅助分类器的识别结果满足P( yi I x)〈P(y31 x)，则令BiB.BFOOlO，即组合分 类器的识别结果是:被识别的电器负载为吹风机。
[0081] 在实施例1中，设某次的主分类器的识别分类标记为?正疋#4 = 0000,则认为主分 类器未能实现电器负载类型识别，且主分类器的识别结果中未能给出识别的电器负载类 型;再设此时辅助分类器的识别结果满足P( yi I x)>P(y21 x)且P(yi I x)>P(y31 x)且P(yi I x)>P (y4|x)，则令關出出4=1000,即组合分类器的识别结果是:被识别的电器负载为白炽灯。
[0082] 组合分类器、主分类器的识别结果分类标记也可以采用其他的方案，例如，分别用 分类标记1、_1，或者是〇、1，或者是_1、1，以及其他方案来表示相应电器负载识别结果为是、 否。组合分类器与主分类器的分类标记方案可以相同，也可以不相同。
[0083] 所述组合分类器的输入特征中，还可以包括电器负载稳态电流有效值Iw。例如，有 2种不同的电器负载，电烙铁和电阻炉需要识别，电烙铁、电阻炉都是纯电阻负载，且都具有 电阻温度修正系数小，电阻值稳定的共同特点。因此，单纯依靠前述的启动电流特征和负载 电流频谱特征、电器负载的基波电压电流相位差特征无法将他们进行区分。输入特征中增 加电器负载稳态电流有效值Iw后，电烙铁功率小，电器负载稳态电流有效值Iw小；电阻炉功 率大，电器负载稳态电流有效值Iw大，特征不同，组合分类器可以进行并完成识别。
[0084]辅助分类器为贝叶斯分类器。可以选择NBC分类器(朴素贝叶斯分类器）、TAN分类 器(树扩展朴素贝叶斯分类器）、BAN分类器(增强的贝叶斯分类器)等三种贝叶斯分类器之 中的一种作为辅助分类器。
[0085]实施例2选择NBC分类器作为辅助分类器。朴素贝叶斯分类的定义如下：
[0086] ⑴设x = {ai，a2，…，am}为一个待分类项，而每个a为x的一个特征属性；
[0087] (2)有类别集合〇={71，72，."，711};
[0088] (3)计算P(yi | x)，P(y21 x)，…，P(yn | x);
[0089] (4)如果P(yk | x) =max{P(yi | x)，P(y21 x)，…，P(yn | x)}，贝iJxG yko
[0090] 计算第⑶步中的各个条件概率的具体方法是：
[0091 ]①找到一个已知分类的待分类项集合作为训练样本集；
[0092] ②统计得到各类别下各个特征属性的条件概率估计；
[0093] P(ai |yi) ,P(a2 |yi) ,??? ,P(am|yi)；
[0094] P(ai | y2)，P(a21 y2)，…，P(am| y2);
[0095] ???；
[0096] P(ai | yn)，P(a21 yn)，…，P(am | yn) 〇 [0097]③根据贝叶斯定理，有：
(1)
[0099]因为分母对于所有类别为常数，因此我们只要将分子最大化即可；又因为在朴素 贝叶斯中各特征属性是条件独立的，所以有：
[0101] 实施例2中，组合分类器的输入特征是{ Fs.，尸7 rM，為.，為，W ,其中，Ts是启 动过程时间，单位是ms; I*是启动电流最大值;Tm是启动电流最大值时间，单位是ms ;Ai、A3为 负载电流频谱特性中的1、3次奇次谐波信号相对幅值#为电器负载的基波电压电流相位 差，单位为度，且基波电压超前于基波电流时#> 〇 <?要求识别的电器负载类别是白炽灯、电 阻炉、电风扇、计算机、吹风机。令朴素贝叶斯分类器的特征属性组合x= {ai，a2，a3，a4，as， a6}中的元素与组合分类器的输入特征集合中的元素按序{ rs, /、rY1，义，木，p }-一 对应;朴素贝叶斯分类器的输出类别集合0={71，72,73, 74,75}则分别与电器负载类别白炽 灯、电阻炉、电风扇、计算机、吹风机一一对应。
[0102] 训练NBC分类器的过程包括：
[0103] 1、对特征属性进行分段划分，进行离散化处理。实施例2中，采取的特征属性离散 化方法是：
[0104] ai： {ai<50,50^ai^100 0,ai>1000}；
[0105] a2：{a2<7,7^a2^11,a2>ll}；
[0106] a3: {a3〈20,20彡a3彡300，a3>300};
[0107] a4： {a4<0.7,0.7^a4^0.9 ,a4>0.9}；
[0108] a5: {a5〈0 ? 02，0 ? 02彡a5彡0 ? 05，a5>0 ? 05};
[0109] as： {ae<-5, -5^as^ 12,ae> 12} 〇
[0110] 2、对每类电器负载类型均采集多组样本作为训练样本，同时计算每类电器负载类 型样本在所有电器负载类型样本中所占有的比例，即分别计算P( yi)、P(y2)、P(y3)、P(y4)、P (y5)。当每类电器负载均采集相同的样本数量时，例如，每类电器负载均采集超过100组的 样本，其中每类电器负载随机选择100组样本作为训练样本，其他则作为测试样本，总的训 练样本为500组，且有
[0111] p(yi)=p(y2)=p(y3)=p( y4)=p(y5)=0.2〇
[0112] 3、计算训练样本每个类别条件下各个特征属性分段的频率(比例），统计得到各类 别下各个特征属性的条件概率估计，即分别统计计算
[0113] P(ai〈50|yi)、P(50彡ai彡1000|yi)、P(ai>1000|yi);
[0114] P(ai〈50|y2)、P(50彡ai彡1000|y2)、P(ai>1000|y 2);
[0115] ???；
[0116] P(ai〈50|y5)、P(50彡ai彡1000|y5)、P(ai>1000|y 5);
[0117] P(a2〈7|yi)、P(7彡 a2彡11 lyO'PUAlllyO;
[0118] P(a2〈7|y2)、P(7彡a2彡11|y 2)、P(a2>ll|y2);
[0119] ???；
[0120] P(a2〈7|y5)、P(7彡a2彡11|y 5)、P(a2>ll|y5);
[0121] P(a3〈20 | yi)、P(20彡a3<300 | yi)、P(a3>300 | yi);
[0122] P(a3〈20 | y2)、P(20彡a3<300 | y2)、P(a3>300 | y2);
[0123] ???；
[0124] P(a3〈20 | y5)、P(20彡a3<300 | y5)、P(a3>300 | y5);
[0125] P(a4〈0.7|yi)、P(0.7<a4<0.9|yi)、P(a4>0.9|yi);
[0126] P(a4〈0.7|y2)、P(0.7<a4<0.9|y2)、P(a4>0.9|y 2);
[0127] ???；
[0128] P(a4〈0.7|y5)、P(0.7<a4<0.9|y5)、P(a4>0.9|y5);
[0129] P(a5〈0.02|yi)、P(0.02<a5<0.05|yi)、P(a5>0.05|yi);
[0130] P(a5〈0.02|y2)、P(0.02<a5<0.05|y2)、P(a5>0.05|y 2);
[0131] P(a5〈0.02|y5)、P(0.02<a5<0.05|y5)、P(a5>0.05|y5);
[0132] P(a6〈_5 | y〇、P(_5彡a6彡 12 | y〇、P(a6彡 12 | ;
[0133] P(a6〈-5|y2)、P(_5彡a6彡12|y 2)、P(a6彡12|y2);
[0134] ???；
[0135] P(a6〈_5|y5)、P(_5彡a6彡12|y 5)、P(a6彡12|y5) 〇
[0136] 经过上述的步骤1、步骤2、步骤3，NBC分类器训练完成。其中，步骤1对特征属性进 行分段划分由人工确定，对每一个输入特征进行分段离散化时，分段的数量为2段或者2段 以上，例如，实施例2中，特征ai_a6都分为3段。每一个特征具体分为多少段，分段阈值的选择 可以根据训练后的贝叶斯分类器对测试样本测试后的结果进行调整。步骤2、步骤3由信息 处理模块102或者是计算机计算完成。
[0137] 本发明中采用贝叶斯分类器进行分类的方法是：
[0138] 1、将组合分类器的输入特征作为贝叶斯分类器的输入特征。在实施例2中，将组合 分类器的输入特征集合{r s，尸，Fm，為，南，_}作为贝叶斯分类器的输入特征X，且有x --(ell ,￡L2 , &3 , &4 , cL5 , cL6 } 〇
[0139] 2、根据训练得到的各类别下各个特征属性的条件概率估计，分别确定各输入特征 属性的分段所在并确定其对每类电器负载类别的概率P( ai|yi)~P(am|yn)，其中，电器负载 类别集合为C= {yi，y2，"_，yn}。实施例2中，电器负载类别集合C= {71，72,73,74,75}对应代表 的电器负载类别是白炽灯、电阻炉、电风扇、计算机、吹风机，确定P(ai | yi)~P(a61 ys)的方 法是采用训练NBC分类器过程中得到的各个特征属性的条件概率估计。
[0140] 3、按照式
[0142]计算每种电器负载类别的后验概率。因为分母P(x)对于所有电器负载类别为常 数，令P(x) = l替代实际的P(x)值，不影响每种电器负载类别后验概率之间的相互大小比 较，此时有
[0150] 采用测试样本对训练好的贝叶斯分类器进行测试，根据测试结果决定是否调整对 输入特征的离散化方法(即调整分段数量与阈值），重新训练贝叶斯分类器。
[0151] 主分类器为支持向量机分类器，或称为SVM分类器。SVM分类器特别适合解决二值 分类情况，因此，主分类器采用多个二类输出SVM分类器组成，每个二类输出SVM分类器对应 识别一种电器负载类型，例如，实施例1中可以采用4个二类输出SVM分类器分别识别白炽 灯、电阻炉、吹风机、计算机，实施例2中可以采用5个二类输出SVM分类器分别识别白炽灯、 电阻炉、电风扇、计算机、吹风机。主分类器选择多个二类输出SVM分类器共同组成时，所有 二类输出SVM分类器的输入特征均为主分类器的输入特征。
[0152] 训练各二类输出SVM分类器时，对每类电器负载类型均采集多组样本，随机抽取部 分作为训练样本，其余的作为测试样本。样本采集采用前述的获取电器负载的启动电流特 征的方法和获取电器负载的负载电流频谱特征以及获取电器负载的基波电压电流相位差 特征的方法。所有的训练样本均作为每个二类输出SVM分类器的训练样本。例如，在实施例2 中，可以分别对白炽灯、电阻炉、电风扇、计算机、吹风机等负载均采集超过100组样本，随机 抽取其中每种100组、共500组样本组成训练样本，剩余的样本组成测试样本；当然，某种负 载或者所有负载采集的样本数量达不到100组样本，SVM分类器也能够取得较好的分类效 果。
[0153]主分类器所选用的二类输出SVM分类器选择径向基RBF核函数，并采用粒子群算法 (PS0)对各二类输出SVM分类器的惩罚参数C和核参数Y进行优化选择。
[0154]每个二类输出SVM分类器只需要完成一种电器负载类型的识别，SVM分类器的训练 相对简单。主分类器由多个二类输出SVM分类器组成，各二类输出SVM分类器之间相互独立， 因此，对某一特征属性进行识别时，主分类器有可能输出的识别结果为唯一的电器负载类 型，或者识别结果为2种或者2种以上电器负载类型，或者未能给出识别的电器负载类型。
【主权项】
1. 一种电器负载类型智能识别装置，其特征在于，包括信息采集模块、信息处理模块、 通信模块； 所述信息采集模块用于采集电器负载的负载电流并转换成电流数字信号;所述电流数 字信号被送至信息处理模块； 所述信息处理模块依据输入的电流数字信号，采用组合分类器进行电器负载类型智能 识别;所述通信模块用于发送信息处理模块的电器负载类型识别结果至上位机； 所述组合分类器的输入特征包括电器负载的启动电流特征、电器负载的负载电流频谱 特征和电器负载的基波电压电流相位差； 所述组合分类器包括支持向量机分类器和贝叶斯分类器； 所述启动电流特征包括启动过程时间、启动电流最大值、启动电流最大值时间。2. 如权利要求1所述的电器负载类型智能识别装置，其特征在于，所述信息采集模块包 括电流传感器、前置放大器、滤波器、A/D转换器;所述信息处理模块的核心为DSP，或者为 ARM，或者为单片机，或者为FPGA。3. 如权利要求2所述的电器负载类型智能识别装置，其特征在于，所述A/D转换器采用 信息处理模块的核心中包括的A/D转换器。4. 如权利要求1所述的电器负载类型智能识别装置，其特征在于，所述通信模块还接收 上位机的相关工作指令;所述通信模块与上位机之间的通信方式包括无线通信方式与有线 通信方式;所述无线通信方式包括ZigBee、蓝牙、WiFi、433MHz数传方式;所述有线通信方式 包括485总线、CAN总线、互联网、电力载波方式。5. 如权利要求1-4中任一项所述的电器负载类型智能识别装置，其特征在于，所述组合 分类器中，支持向量机分类器为主分类器，贝叶斯分类器为辅助分类器。6. 如权利要求5所述的电器负载类型智能识别装置，其特征在于，所述组合分类器进行 电器负载类型智能识别的方法是：当主分类器成功实现电器负载类型识别时，主分类器的 电器负载类型识别结果为组合分类器的识别结果；当主分类器未能实现电器负载类型识 另Ij，且主分类器的识别结果为2种或者2种以上电器负载类型，将主分类器输出的2种或者2 种以上电器负载类型识别结果中，辅助分类器输出中概率最高的电器负载类型作为组合分 类器的电器负载类型识别结果；当主分类器未能实现电器负载类型识别，且主分类器的识 别结果中未能给出识别的电器负载类型时，将辅助分类器输出中概率最高的电器负载类型 作为组合分类器的电器负载类型识别结果。7. 如权利要求5所述的电器负载类型智能识别装置，其特征在于，所述负载电流频谱特 征通过以下方法获得： 步骤一、获取电器负载的稳态电流信号，并将其转换为对应的稳态电流数字信号； 步骤二、对稳态电流数字信号进行傅立叶变换，得到负载电流频谱特性； 步骤三、将负载电流频谱特性中谐波次数为η次的奇次谐波信号相对幅值作为负载电 流频谱特征，其中，η = 1，3，···，Μ;所述M表示谐波最高次数且M大于等于3。8. 如权利要求5所述的电器负载类型智能识别装置，其特征在于，所述启动电流特征由 信息处理模块通过以下方法获得： 步骤1、电器负载启动前，开始对电器负载的负载电流连续采样并对负载电流大小进行 判断；当负载电流有效值大于ε时，判定电器负载开始启动并转向步骤2;所述ε为大于〇的数 值； 步骤2、对电器负载的负载电流进行连续采样，以工频周期为单位计算负载电流有效值 并保存;计算最近N个工频周期的负载电流有效值的平均值；当最近N个工频周期之内的每 个工频周期的负载电流有效值与该N个工频周期的负载电流有效值的平均值相比较，波动 幅度均小于设定的相对误差范围E时，判定电器负载进入稳定状态，转向步骤3;所述N的取 值范围为50-500;所述E的取值范围为2%-20% ; 步骤3、将最近N个工频周期之内的负载电流有效值的平均值作为电器负载稳态电流有 效值;将电器负载开始启动时刻至最近N个工频周期起始时刻之间的时间作为启动过程时 间；将电器负载开始启动时刻至启动过程时间之内负载电流有效值最大的工频周期之间的 时间作为启动电流最大值时间;将启动电流最大值时间所在工频周期的负载电流有效值与 电器负载稳态电流有效值之间的比值作为启动电流最大值。9. 如权利要求8所述的电器负载类型智能识别装置，其特征在于，所述组合分类器的输 入特征还包括电器负载稳态电流有效值。10. 如权利要求5所述的电器负载类型智能识别装置，其特征在于，所述电器负载的基 波电压电流相位差通过以下方法获得： 步骤①、待电器负载进入稳定状态后，同步获取电器负载的稳态电压信号、稳态电流信 号，并将其转换为相应的稳态电压数字信号、稳态电流数字信号； 步骤②、对稳态电压数字信号、稳态电流数字信号分别进行数字滤波，提取出基波电压 信号、基波电流信号； 步骤③、分析计算基波电压信号与基波电流信号之间的相位差，将基波电压信号与基 波电流信号之间的相位差作为电器负载的基波电压电流相位差。
【文档编号】G06K9/00GK105913005SQ201610213352
【公开日】2016年8月31日
【申请日】2016年4月8日
【发明人】肖伸平, 练红海, 凌云, 陈刚
【申请人】湖南工业大学