发明名称 基于故障类型分类能力评价矩阵的多分类器融合故障诊断方法
摘要 本发明公开了基于故障类型分类能力评价矩阵的多分类器融合故障诊断方法。本发明基于分类器度量级输出形式,提出了以分类器输出结果熵值评价分类器多故障类型分类能力的方法,计算得到评价矩阵,并构建了基于模糊综合评判方式的多分类器融合基本模型,进行决策级融合,得出最终的诊断结论。
申请公布号 CN104484678A 申请公布日期 2015.04.01
申请号 CN201510007001.3 申请日期 2015.01.07
申请人 青岛理工大学 发明人 文妍;谭继文;战卫侠;战红;孙显彬
分类号 G06K9/62(2006.01)I 主分类号 G06K9/62(2006.01)I
代理机构 北京中北知识产权代理有限公司 11253 代理人 段秋玲
主权项 基于故障类型分类能力评价矩阵的多分类器融合故障诊断方法,包括以下步骤:设故障类型有m种,则组成模式空间D可以记为D=C<sub>1</sub>∪C<sub>2</sub>∪...∪C<sub>m</sub>其中C<sub>i</sub>,<img file="FDA0000652900470000019.GIF" wi="409" he="75" />成为一个类,且要求C<sub>1</sub>∩C<sub>2</sub>∩...∩C<sub>m</sub>=Φ;若用J个分类器e<sub>j</sub>(j=1,2,...,J)对来自模式空间D的样本x进行分类,分类器e<sub>j</sub>的输出可以记为y<sub>j</sub>=e<sub>j</sub>(x)分类器e<sub>j</sub>的输出形式为<img file="FDA0000652900470000011.GIF" wi="486" he="89" />该向量给每一个类别标签分配了一个数值,用来度量样本x属于该类的程度,即该种故障类型的发生概率;这里要求<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mn>0</mn><mo>&le;</mo><msubsup><mi>y</mi><mi>j</mi><mi>i</mi></msubsup><mo>&le;</mo><mn>1</mn><mo>,</mo></mrow>]]></math><img file="FDA0000652900470000012.GIF" wi="236" he="90" /></maths>且<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msubsup><mi>y</mi><mi>j</mi><mi>i</mi></msubsup><mo>=</mo><mn>1</mn><mo>;</mo></mrow>]]></math><img file="FDA0000652900470000013.GIF" wi="206" he="134" /></maths>融合步骤:(1)为了衡量分类器分类能力,需要一定数目N样本进行测试;在测试数据中随机抽取N个样本,要求m种故障类型的样本数目比为1:1:...:1;(2)对于故障类型C<sub>i</sub>,<img file="FDA0000652900470000014.GIF" wi="444" he="75" />按照下式分别计算J个分类器输出的熵<img file="FDA0000652900470000015.GIF" wi="111" he="90" />求其均值<img file="FDA0000652900470000016.GIF" wi="94" he="93" />作为第j分类器故障类型C<sub>i</sub>输出的熵值;<img file="FDA0000652900470000017.GIF" wi="835" he="253" />上式中常数K的取值与故障类型的数目m有关;因为根据信息熵的定义,若分类器输出<img file="FDA0000652900470000018.GIF" wi="482" he="129" />该输出认为样本m种故障的概率相同,也就是说该判断无效,此时熵具有最大值;为了分析方便,使熵值在[0,1]范围内,因此取<img file="FDA0000652900470000021.GIF" wi="225" he="126" />如果分类器对某样本的判断类型与期望类型不同时,意味着诊断错误,也可以理解为该判断无效,因此这种情况下,直接令熵值为最大值1;(3)偏差度<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>d</mi><msub><mi>y</mi><mi>j</mi></msub><msub><mi>C</mi><mi>i</mi></msub></msubsup><mo>=</mo><mn>1</mn><mo>-</mo><msubsup><mover><mi>H</mi><mo>&OverBar;</mo></mover><msub><mi>y</mi><mi>j</mi></msub><msub><mi>C</mi><mi>i</mi></msub></msubsup><mo>,</mo><mrow><mo>(</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>J</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000652900470000022.GIF" wi="789" he="93" /></maths>对于故障类型C<sub>i</sub>,第J个分类器的权重为<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mi>&omega;</mi><msub><mi>y</mi><mi>j</mi></msub><msub><mi>C</mi><mi>i</mi></msub></msubsup><mo>=</mo><mfrac><msubsup><mi>d</mi><msub><mi>y</mi><mi>j</mi></msub><msub><mi>C</mi><mi>i</mi></msub></msubsup><mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>J</mi></munderover><msubsup><mi>d</mi><msub><mi>y</mi><mi>j</mi></msub><msub><mi>C</mi><mi>i</mi></msub></msubsup></mrow></mfrac><mo>,</mo><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>m</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000652900470000023.GIF" wi="753" he="235" /></maths>(4)得到权重矩阵A为<img file="FDA0000652900470000024.GIF" wi="644" he="331" />矩阵A中的元素<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msubsup><mi>&omega;</mi><msub><mi>y</mi><mi>j</mi></msub><msub><mi>C</mi><mi>i</mi></msub></msubsup><mo>,</mo><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>m</mi><mo>;</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>J</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000652900470000025.GIF" wi="948" he="94" /></maths>代表分类器y<sub>j</sub>对故障类型C<sub>i</sub>,<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mo>&ForAll;</mo><mi>i</mi><mo>&Element;</mo><mi>&Lambda;</mi><mo>=</mo><mo>{</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>m</mi><mo>}</mo></mrow>]]></math><img file="FDA0000652900470000026.GIF" wi="417" he="69" /></maths>的权重;(4)对于某一样本X,J个分类器的输出组成待评判矩阵R<img file="FDA0000652900470000027.GIF" wi="555" he="325" />(5)因此多分类器融合的结果为B<sub>m×m</sub>=A·R取该矩阵的对角线元素组成向量B<sub>out</sub>,根据最大隶属原则,确定故障类别。
地址 266000 山东省青岛市市北区抚顺路11号