xinxishang.rar_信息熵_故障分类_熵 matlab

clc; clear; tic %%% 高压电器 %%% 信息增益计算方法 %%% 测试手段 % A1 油色谱分析 % A2 各绕组直流电阻 % A3 介质损耗 % A4 空载损耗及电流 % A5 绕组对地绝缘电阻 % A6 铁心接地电流 % A7 铁心层间电阻 %%% 故障类型 % C1 线圈层间、匝间短路 % C2 线圈接地 % C3 线圈与铁心接地 % C4 铁心层间短路 % C5 铁心多点接地 % C6 铁心接缝不良、硅钢片填充不足 % C7 裸金属过热 % C8 火花放电 % C9 变压器油流受阻过热 % s 样本总数 s=35 a11=5 %故障类型为线圈层间、匝间短路的样本数 a21=6 %故障类型为线圈接地的样本数 a22=2 %故障类型为线圈与铁心接地的样本数 a23=2 %故障类型为铁心层间短路的样本数 a24=6 %故障类型为铁心多点接地的样本数 a25=2 %故障类型为铁心接缝不良、硅钢片填充不足的样本数 a31=8 %故障类型为裸金属过热的样本数 a41=2 %故障类型为火花放电的样本数 a51=2 %故障类型为变压器油流受阻过热的样本数 %计算系统初始熵 b1=log2(a11/s) b2=log2(a21/s) b3=log2(a22/s) b4=log2(a23/s) b5=log2(a24/s) b6=log2(a25/s) b7=log2(a31/s) b8=log2(a41/s) b9=log2(a51/s) Es=-(a11/s*b1+a21/s*b2+a22/s*b3+a23/s*b4+a24/s*b5+a25/s*b6+a31/s*b7+a41/s*b8+a51/s*b9) %系统初始熵值 %计算油色谱分析的系统平均熵 b11=log2(a11/a11) e1=-a11/a11*b11 a2=a21+a22+a23+a24+a25 %油色谱分析:过热的样本数 b21=-log2(a21/a2) b22=-log2(a22/a2) b23=-log2(a23/a2) b24=-log2(a24/a2) b25=-log2(a25/a2) e2=(a21/a2)*b21+(a22/a2)*b22+(a23/a2)*b23+(a24/a2)*b24+(a25/a2)*b25 b31=log2(a31/a31) e3=-a31/a31*b31 b41=log2(a41/a41) e4=-a41/a41*b41 b51=log2(a51/a51) e5=a51/a51*b51 H1=a21/s*e1+a2/s*e2+a31/s*e3+a41/s*e4+a51/s*e5 %油色谱分析的系统平均熵值 GainA1=Es-H1 %油色谱分析的信息增益值 a1=5 %故障类型为线圈层间、匝间短路的样本数 a21=6 %故障类型为线圈接地的样本数 a22=2 %故障类型为线圈与铁心接地的样本数 a23=2 %故障类型为铁心层间短路的样本数 a24=6 %故障类型为铁心多点接地的样本数 a25=2 %故障类型为铁心接缝不良、硅钢片填充不足的样本数 a26=8 %故障类型为裸金属过热的样本数 a27=2 %故障类型为火花放电的样本数 a28=2 %故障类型为变压器油流受阻过热的样本数 a2=30 %绕组直流电阻平衡的样本数 e1=-a1/a1*log2(a1/a1) e2=-a21/a2*log2(a21/a2)-a22/a2*log2(a22/a2)-a23/a2*log2(a23/a2)-a24/a2*log2(a24/a2)-a25/a2*log2(a25/a2)-a26/a2*log2(a26/a2)-a27/a2*log2(a27/a2)-a28/a2*log2(a28/a2) H2=a1/s*e1+a2/s*e2 %A2的系统平均熵值 GainA2=Es-H2 %A2的信息增益值 a11=5 %故障类型为线圈层间、匝间短路的样本数 a21=6 %故障类型为线圈接地的样本数 a22=2 %故障类型为线圈与铁心接地的样本数 a12=2 %故障类型为铁心层间短路的样本数 a13=6 %故障类型为铁心多点接地的样本数 a14=2 %故障类型为铁心接缝不良、硅钢片填充不足的样本数 a15=8 %故障类型为裸金属过热的样本数 a23=2 %故障类型为火花放电的样本数 a24=2 %故障类型为变压器油流受阻过热的样本数 a1=a11+a12+a13+a14+a15 %介质损耗正常的样本数 a2=a21+a22+a23+a24 %介质损耗超标的样本数 e1=-a11/a1*log2(a11/a1)-a12/a1*log2(a12/a1)-a13/a1*log2(a13/a1)-a14/a1*log2(a14/a1)-a15/a1*log2(a15/a1) e2=-a21/a2*log2(a21/a2)-a22/a2*log2(a22/a2)-a23/a2*log2(a23/a2)-a24/a2*log2(a24/a2) H3=a1/s*e1+a2/s*e2 %A3的系统平均熵值 GainA3=Es-H3 %A3的信息增益值 a11=5 %故障类型为线圈层间、匝间短路的样本数 a12=6 %故障类型为线圈接地的样本数 a13=2 %故障类型为线圈与铁心接地的样本数 a14=2 %故障类型为铁心层间短路的样本数 a15=6 %故障类型为铁心多点接地的样本数 a16=2 %故障类型为铁心接缝不良、硅钢片填充不足的样本数 a21=8 %故障类型为裸金属过热的样本数 a22=2 %故障类型为火花放电的样本数 a23=2 %故障类型为变压器油流受阻过热的样本数 a1=a11+a12+a13+a14+a15+a16 %空载损耗及电流正常的样本数 a2=a21+a22+a23 %空载损耗及电流大的样本数 e1=-a11/a1*log2(a11/a1)-a12/a1*log2(a12/a1)-a13/a1*log2(a13/a1)-a14/a1*log2(a14/a1)-a15/a1*log2(a15/a1)-a16/a1*log2(a16/a1) e2=-a21/a2*log2(a21/a2)-a22/a2*log2(a22/a2)-a23/a2*log2(a23/a2) H4=a1/s*e1+a2/s*e2 %A4的系统平均熵值 GainA4=Es-H4 %A4的信息增益值 a11=5 %故障类型为线圈层间、匝间短路的样本数 a21=6 %故障类型为线圈接地的样本数 a22=2 %故障类型为线圈与铁心接地的样本数 a12=2 %故障类型为铁心层间短路的样本数 a13=6 %故障类型为铁心多点接地的样本数 a14=2 %故障类型为铁心接缝不良、硅钢片填充不足的样本数 a15=8 %故障类型为裸金属过热的样本数 a16=2 %故障类型为火花放电的样本数 a17=2 %故障类型为变压器油流受阻过热的样本数 a1=a11+a12+a13+a14+a15+a16+a17 %绕组对地绝缘电阻正常的样本数 a2=a21+a22 %绕组对地绝缘电阻下降的样本数 e1=-a11/a1*log2(a11/a1)-a12/a1*log2(a12/a1)-a13/a1*log2(a13/a1)-a14/a1*log2(a14/a1)-a15/a1*log2(a15/a1)-a16/a1*log2(a16/a1)-a17/a1*log2(a17/a1) e2=-a21/a2*log2(a21/a2)-a22/a2*log2(a22/a2) H5=a1/s*e1+a2/s*e2 %A5的系统平均熵值 GainA5=Es-H5 %A5的信息增益值 a11=5 %故障类型为线圈层间、匝间短路的样本数 a12=6 %故障类型为线圈接地的样本数 a21=2 %故障类型为线圈与铁心接地的样本数 a13=2 %故障类型为铁心层间短路的样本数 a22=6 %故障类型为铁心多点接地的样本数 a14=2 %故障类型为铁心接缝不良、硅钢片填充不足的样本数 a15=8 %故障类型为裸金属过热的样本数 a16=2 %故障类型为火花放电的样本数 a17=2 %故障类型为变压器油流受阻过热的样本数 a1=a11+a12+a13+a14+a15+a16+a17 %A6无的样本数 a2=a21+a22 %A6有的样本数 e1=-a11/a1*log2(a11/a1)-a12/a1*log2(a12/a1)-a13/a1*log2(a13/a1)-a14/a1*log2(a14/a1)-a15/a1*log2(a15/a1)-a16/a1*log2(a16/a1)-a17/a1*log2(a17/a1) e2=-a21/a2*log2(a21/a2)-a22/a2*log2(a22/a2) H6=a1/s*e1+a2/s*e2 %A6的系统平均熵值 GainA6=Es-H6 %A6的信息增益值 a11=5 %故障类型为线圈层间、匝间短路的样本数 a12=6 %故障类型为线圈接地的样本数 a13=2 %故障类型为线圈与铁心接地的样本数 a21=2 %故障类型为铁心层间短路的样本数 a22=6 %故障类型为铁心多点接地的样本数 a14=2 %故障类型为铁心接缝不良、硅钢片填充不足的样本数 a15=8 %故障类型为裸金属过热的样本数 a16=2 %故障类型为火花放电的样本数 a17=2 %故障类型为变压器油流受阻过热的样本数 a1=a11+a12+a13+a14+a15+a16+a17 %A7正常的样本数 a2=a21+a22 %A6下降的样本数 e1=-a11/a1*log2(a11/a1)-a12/a1*log2(a12/a1)-a13/a1*log2(a13/a1)-a14/a1*log2(a14/a1)-a15/a1*log2(a15/a1)-a16/a1*log2(a16/a1)-a17/a1*log2(a17/a1) e2=-a21/a2*log2(a21/a2)-a22/a2*log2(a22/a2) H7=a1/s*e1+a2/s*e2 %A7的系统平均熵值 GainA7=Es-H7 %A7的信息增益值 Es H=[H1 H2 H3 H4 H5 H6 H7] A=[GainA1 GainA2 GainA3 GainA4 GainA5 GainA6 GainA7] b=max(A)