一种引入工况信息的电路故障诊断扩展决策树构建方法

摘要

本发明提出一种引入工况信息的电路故障诊断扩展决策树构建方法，属于故障诊断技术领域。包括以下步骤：步骤一，确定被测电路的测试点、工况组成和电路状态种类；步骤二，根据电路历史运行数据或者仿真数据，获取被测电路在各电路状态和各电路工况下的电路测试点数据集；步骤三，构建扩展决策树；步骤四，应用决策树，诊断电路。本发明通过采集各测试点信息和工况种类信息，构建一个既包括测试信息、也包括工况类型的特征矩阵，再根据基尼指数划分属性，构建扩展决策树，实现多工况综合故障诊断。与传统方法中每个工况建立一个决策树相比较，本发明方法只需构建一个决策树，大大减少了诊断资源，有效提高了诊断效率。

说明书

技术领域

本发明涉及一种引入工况信息的电路故障诊断扩展决策树构建方法，属于故障诊断技术领域。

背景技术

电路故障诊断，主要是通过测试点的设计和布置来监测系统的运行状态，并在诊断出异常后进行分析处理。工况是指电路的工作状态。在实际的生活中一旦负载、输入等发生改变，那么电路工况也会发生变化。故障诊断对于电路的维修和后期的保障有很重要的作用。传统的决策树故障诊断方法只考虑单一工作状况下的情况，却很少考虑多种工况下的情况。在实际应用中，电路普遍存在着多工况频繁切换的过程，电路的诊断模型应随着多工况的切换而变化。

发明内容

本发明为了提高多工况电路决策树故障诊断效率，提出一种引入工况信息的电路故障诊断扩展决策树构建方法，主要是通过采集各测试点信息和工况种类信息，构建一个既包括测试信息、也包括工况类型的特征矩阵，再根据基尼指数划分属性，构建扩展决策树，实现多工况综合故障诊断。与传统方法中每个工况建立一个决策树相比较，本发明方法只需构建一个决策树，大大减少了诊断资源，有效提高了诊断效率。

本发明提出的一种引入工况信息的电路故障诊断扩展决策树构建方法，具体包括以下几个步骤：

步骤一，确定被测电路的测试点、工况组成和电路状态种类。

设被测电路有n个测试点，有m个工况，有p种故障状态，加上健康状态，则被测电路共有(p+1)种电路状态。n、m、p均为正整数。

步骤二，根据电路历史运行数据或者仿真数据，获取被测电路在各电路状态和各电路工况下的电路测试点数据集。

将每个工况用唯一的m位的二进制序列编码，二进制序列中只有一位是1，其他位都是0。设获得测试数据集D'如下：

D'＝{(x₁₁,y₁),(x₁₂,y₁),…,(x_1n,y₁),(x_1(n+1),y₁),…,(x_1(n+m),y₁),(x₂₁,y₂),(x₂₂,y₂),…,

(x_2n,y₂),(x_2(n+1),y₂),…,(x_2(n+m),y₂),…,(x_d(n+m),y_d)}

其中，d为样本个数，x_in表示第i个样本中测试点T_n所测得的数据值，y_i表示采集第i个样本时的电路状态，x_i(n+1),…,x_i(n+m)表示采集第i个样本时的电路工况，为二进制序列。

步骤三，以数据集D'为训练样本构建扩展决策树；

利用分类回归树法，以数据集D'中的n个测试点T₁,T₂,…T_n和m位工况编码A₁,A₂,…,A_m为属性，以电路状态S₀,S₁……S_p为所属类别，构建扩展决策树；

步骤四，采集被测电路测试点数据，将测试点数据和工况编码依据扩展决策树，诊断电路状态。

本发明与现有技术相比，具有以下明显优势：

(1)本发明只需构建一个决策树既可完成故障诊断，节省了运算资源；

(2)本发明通过单个决策树对电路状态进行分类，提高了诊断效率；

(3)本发明能够对多工况电路系统进行故障诊断。

附图说明

图1是本发明实施例中的电路两级负反馈放大电路图；

图2是本发明实施例中的电路模拟输入不同幅值作为不同工况的电压图；

图3是本发明的引入工况信息的电路故障诊断扩展决策树构建方法的整体步骤示意图；

图4是本发明实施例数据经过CART法得到的扩展决策树。

具体实施方式

为了便于本领域普通技术人员理解和实施本发明，下面结合附图对本发明作进一步的详细描述。

本发明实施例以图1所示的多级负反馈放大电路为研究对象，来说明本发明方法的具体实现步骤，图1所示电路由两个三极管和其附属外围器件组成，耦合方式采用电容耦合。

其中，Q1三极管及其附属器件组成一级放大，Q2三极管及其附属器件组成二级放大，之间用电容C2耦合。输入信号为正弦波Vsin，频率为1kHz。电容C6和电阻R8组成反馈回路，提升电路的稳定性。在该电路中设置了8个测试点，通过模拟探针监测所连接点处的电压。

本发明实施例通过设置不同输入电压幅值(频率不变)模拟工况的变化，如图2所示，拟定为3种不同的输入，10mv(1kHz)、100mv(1kHz)、500mv(1kHz)，对应三种不同的工况。设置8个测试点采集电压信息，对电路中相关器件进行故障模拟注入。电路的状态如表1所示，包括9种。

表1电路的状态列表

Class4C6短路Class5R1-30％Class6R7-30％Class7R8-30％Class8R8开路

下面结合上面所述实施例，来说明本发明提出的一种引入工况信息的电路故障诊断扩展决策树构建方法的如图3所示的各步骤。

步骤一，确定电路的测试点组成、工况组成和状态组成。

根据电路的设计资料，确定对象产品的测试点组成，工况组成，以及正常和故障状态组成，并分别记录到电路测试点组成表、电路工况组成表和电路状态组成表中。设被测电路有n个测试点T₁,T₂……T_n，有m个工况O₁,O₂……O_m，有p种故障状态，加上健康状态，则被测电路共有(p+1)种电路状态S₀,S₁……S_p。n、m、p均为正整数。

针对本发明实施例，如表2至表4所示。表2中左侧列为测试点的编号，右侧列为测试点位置的描述；表3左侧列为工况的编号，右侧列为工况内容的描述；表4中左侧列为电路工作状态的编号，右侧列为电路健康或故障状态的描述。

表2电路测试点组成

测试点描述T₁C1与Q1之间的电路T₂R3与Q1之间的电路T₃R4与R5之间的电路T₄C2与Q2之间的电路T₅C4与R11之间的电路T₆R10与Q2之间的电路T₇R9与Q2之间的电路T₈C6与Q1之间的电路

表3电路工况组成

工况描述O₁输入电压幅值为10mv，频率为1kHz的正弦波O₂输入电压幅值为100mv，频率为1kHz的正弦波O₃输入电压幅值为500mv，频率为1kHz的正弦波

表4电路状态组成

状态描述S₀正常S₁R4短路S₂R5开路S₃C3短路S₄C6短路S₅R1-30％S₆R7-30％S₇R8-30％S₈R8开路

本发明实施例的被测电路有8个测试点；有3个工况；有8种故障状态，加上健康状态，则该电路系统共有9种电路状态。

步骤二，根据电路历史运行数据及电路仿真补充数据，确定对象产品在各电路状态和各电路工况下的电路测试点数据集。

设从n个测试点采集的一个样本i表示为(x_i1,x_i2,...,x_in)，采集该样本时对应的电路状态标记为y_i，y_i是(p+1)种电路状态中的一种，则得到的数据集如下：

数据集D＝{(x₁₁,y₁),(x₁₂,y₁),…,(x_1n,y₁),(x₂₁,y₂),(x₂₂,y₂),…,(x_2n,y₂),…,(x_dn,y_d)}；

d表示样本总数。

本步骤包括如下子步骤2.1～2.3。

步骤2.1，根据电路的历史运行数据，确定出不同状态、不同工况下的测试点数据。对于每个状态与工况组合情况下，对每个测试点结果，都分别选取最大值和最小值情况作为样本数据，记录在测试点样本数据表中。

步骤2.2，对于历史数据中缺少的状态工况组合数据，通过电路仿真补充获取。同样，对于每个状态与工况组合情况下，对每个测试点结果，都分别选取最大值和最小值情况作为样本数据，记录在测试点样本数据表中。

具体步骤如下：

(1)采用商用仿真工具软件，建立电路的状态工况组合情况的仿真模型；

(2)在电路模型中的测试点位置上放置探针；

(3)开始仿真，获得测试点数据，对每个测试点，都分别选取最大值和最小值情况作为样本数据。将数据记录到表5之中。

(4)所有需要仿真的状态工况组合，仿真完毕后，结束仿真。

如表5所示，为本发明实施例获得的部分测试点样本数据表，表中的第一行为测试点编号，每一列(除了最后一列)为对应第一行测试点的数据，最后一列为电路的状态标签。

表5测试点样本数据表

T₁T₂T₃T₄T₅T₆T₇T₈状态5.6151697.8228894.8053823.9966526.7920935.97181010.3474630.082870S₂1.9280461.1892511.1085180.0000006.7916345.97140310.2821690.017585S₃3.09686623.6623742.3015972.3015977.1609896.32577815.3491884.954675S₂5.35333415.5134824.5323053.3903825.1877854.36140023.6257307.414710S₇5.6983568.4018944.8856243.9929895.6916464.85934722.8086408.862589S₄5.6151697.8229444.8053813.9966535.5572824.74290013.9976520.082830S₇4.2234563.7543223.3364283.1452985.5264864.85163215.2485610.775422S₀

步骤2.3，对工况进行编码，并加入到样本数据中。

(1)将工况种类以二进制序列来表示，并记录在工况编码表中。

每种工况对应一个唯一的m位的二进制序列编码，二进制序列中只有一位是1，其他位都是0。例如，设第i个工况O_i的二进制序列为

本发明实施例中共有3种工况，将工况O_i编码为并记录到表6中，表6的第一列表示工况，第二列表示工况编码。

表6工况编码表

工况工况编码O₁001O₂010O₃100

(2)将工况编码加入到样本数据中，本发明实施例的部分测试样本数据如表7所示。

表7将工况编码加入到测试数据后的样本数据

表7是在表5的基础上，增加了3列工况信息。由此组成新的数据集D'。

D'＝{(x₁₁,y₁),(x₁₂,y₁),…,(x_1n,y₁),(x_1(n+1),y₁),…,(x_1(n+m),y₁),(x₂₁,y₂),(x₂₂,y2),…,

(x_2n,y₂),(x_2(n+1),y₂),…,(x_2(n+m),y₂),…,(x_d(n+m),y_d)}

其中，相对于数据集D，D'中增加了工况的数据信息，对于样本i的数据增加了(x_i(n+1),y_i),...,(x_i(n+m),y_i)。

步骤三，构建扩展决策树。

以数据集D'为训练样本利用标准CART(Classification And Regression Tree，分类回归树)法构建决策树。以n个测试点数据T₁,T₂,…T_n和m位工况编码A₁,A₂,…,A_m为属性值，共(n+m)个属性，以电路状态S₀,S₁……S_p为所属类别。

本发明实施例中训练数据集D'中共有204个训练样本；以测试点数据(每个测点表示一个属性，因为有8个测点，所以有8个属性)和工况编码(每一列编码作为一个属性，因为每个工况编码有3位数，所以有3个属性)为属性值，共有11个属性(T₁,T₂,T₃,T₄,T₅,T₆,T₇,T₈,A₁,A₂,A₃)。以电路状态为分类类别，共有8种故障状态和健康状态共9类。利用标准CART法构建决策树。

数据集D'的基尼值Gini(D')为：

其中，p_k代表样本集D'中第k类样本所占的比例；p_k'代表样本集D'中第k'类样本所占的比例；k＝1,2,…,|γ|，|γ|代表样本集D'的样本种类个数。

属性T_i(i＝1,2,…,n)和A_j(j＝1,2,…,m)的基尼指数为：

其中，V表示属性T_i或A_j的划分种类个数，对应本发明实施例中V为9；D^'v表示属性第v个分类的样本集，|D^'v|和|D'|分别表示集合D^'v和D'的样本个数。

将所计算的各属性的基尼指数记录到根节点或叶节点各属性基尼指数表中。本发明实施例如表8所示，表中的第一行表示属性编号，第二行对应第一行属性的根节点或叶的节点基尼指数。

表8根节点或叶节点各属性基尼指数表

属性T₁T₂T₃T₄T₅T₆T₇T₈A₁A₂A₃基尼指数0.070.020.050.0450.0770.120.030.250.0240.050.053

在候选属性集合中，选择基尼指数最小的对应属性作为最优划分属性。这里选择T₂作为根节点划分属性。从节点依据属性值将数据集划分为两支，对每支数据集进行重复属性划分工作，依次将节点属性基尼指数记录到对应的表格之中，重复以上操作，直到将数据集划分到数据集中的样本都表示一种电路状态为止，决策树构建完成。在本发明实施例中，构建成的决策树如图4所示。图4中，根节点T₂所对应的属性值为3.824932，3.824932的由来是：将采集到的所有样本中的T₂属性值按从小到大排列，分别取相邻两个属性值的中间值，每个中间值都对应一个基尼值，这里选取3.824932作为T₂属性值分类的分类点是因为——3.824932对应的基尼值在所有中间值中最小。其他属性的分类点值类似，都根据CART算法和样本数据来确定。属性值大小代表测试点采集到的电压值。

步骤四，应用决策树，诊断电路。

采集被测电路测试点数据和工况种类编码信息，将所采集的数据记录到表9之中，表9的第一行表示各属性编号，第二、三行表示各属性值。

表9信息采集得到的各属性值

属性T₁T₂T₃T₄T₅T₆T₇T₈A₁A₂A₃属性值4.6487213.5978135.2597884.8294234.5147834.49765322.4513820.841123010属性值4.2365483.6548715.3478563.1145244.8532155.648529723.5485840.758463100

最后，依据表9中的第二行属性值数据，诊断出电路在100mv的工况下电路发生故障“S₁”；依据表9中的第三行属性值数据，诊断出电路在10mv的工况下电路健康。