一种基于网络拓扑结构的故障诊断测点选择方法
一种基于网络拓扑结构的故障诊断测点选择方法
【技术领域】
[0001] 本发明属于电子电路系统的可测性设计领域，应用二分复杂网络模型，设及一种 基于网络拓扑结构的故障诊断测点选择方法。
【背景技术】
[0002] 模拟电路系统故障测试与诊断中测点选择是其中一个重要研究课题。随着电路系 统复杂度与集成度的不断增加，模拟电路系统测点选择也越来越复杂。
[0003] 目前，复杂网络理论正是研究具有大规模，高复杂程度网络的一口交叉性和综合 性学科。但是到目前为止，复杂网络理论和电子电路系统的故障诊断、测点选择所结合，得 到的研究成果还较少，主要研究成果仅限于电子电路系统的小世界性。并且，大都采用传统 的网络建模方法，只能研究电路系统的单一网络模型:元器件网络模型或者测点网络模型， 运种方法丢失了原始电路系统的许多有用信息。而二分网络模型比传统的单顶点网络模型 在研究模拟电路系统上更有优势，不仅能够弥补单顶点网络模型不能够区分电路中拓扑同 构的缺陷，也能够得到除元器件网络W外的测点网络。
[0004] 在复杂网络中，网络节点中屯、性是用来衡量网络节点重要性的一个指标。本发明 立足拓扑结构分析网络节点中屯、性，主要有7种，分别是:节点度中屯、性判断方法、介数中屯、 性判断方法、凝聚度中屯、性判断方法、网络流中屯、性判断方法、随机行走中屯、性判断方法、 子图中屯、性判断方法W及特征向量中屯、性判断方法。方法概略为首先计算测点网络的各个 节点中屯、性，然后应用模糊数学中的一些模糊决策方法(模糊意见集中决策、模糊二元对比 决策等），得到网络中各个测点重要性排序，最后讨论网络中测点重要性排序与测点选择的 关系。
【发明内容】
[0005] 本发明的目的在于克服现有技术的不足，提供一种基于网络拓扑结构的故障诊断 测点选择方法，使城域网与接入网得到有效融合，同时具有可扩展、高灵活W及低成本的性 能。
[0006] 为实现上述发明目的，本发明基于网络拓扑结构的故障诊断测点选择方法，其特 征在于，包括W下步骤：
[0007] (1 )、计算网络拓扑中各个备选测点的中屯、性系数；
[000引（2)、将中屯、性系数归一化到0-1，得到归一化的中屯、性系数矩阵；
[0009] (3)、根据中屯、性系数系数矩阵，利用Borda数法和平均法计算各个备选测点的 Borda数和隶属函数值，得到备选测点的两种重要性排序表；
[0010] (4)、根据备选测点的两种重要性排序表，运行基于节点中屯、性的测点选择算法得 到最优测点集Snpt，选择测点数更少的Snpt作为最终选择结果。
[ocm]其中，所述的中屯、性系数包括:节点度、介数、凝聚度、网络流、随机行走、子图和特 征向量；
[0012] 所述的节点度中屯、性系数的计算方法为：
[0013] 设备选测点i的度是k( i)，则有：
[0014] Cd(i)=k(i)；
[0015] 所述的介数中屯、性系数的计算方法为：
[0016]
[0017] 其中，s<t表示Ws为源节点，t为目标节点，ns康示源节点S和目标节点帥存在所 有最短路径的路径数，η表示网络中的节点数，gst,i表示经过备选测点i的节点对(s，t)的最 短路径数目；
[0018] 所述的凝聚度中屯、性系数的计算方法为：
[0019]
[0020] 其中，dst表示备选测点i到网络中其他所有节点的最短距离；
[0021] 所述的网络流中屯、性系数的计算方法为：
[0022]
[0023] 其中，gst表示网络中节点对(s，t)之间的所有路径数，但不包含回路；
[0024] 所述的随机行走中屯、性系数的计算方法为：
[002引（a)、计算Laplace矩阵L，L = D-A，其中，D表示节点度组成的对角矩阵，A表示目标 网络的邻接矩阵；
[0026] 化）、除去L的最后一行、最后一列，得到矩阵L%
[0027] (C)、计算矩阵L嘴逆矩阵L^，再对逆矩阵增添一个行向量和一个列向量，且 添加的向量分量均为0，得到矩阵T;
[00%] (d)、得到随机行走中屯、性系数的计算公式：
[0029]
[0030] 其中，j表示目标网络测点编号，Αυ表示目标网络邻接矩阵的第i行第j列元素， Tls、Tlt、Ls、Tパ的下标表示其所在矩阵T的行列位置，Ist,康示从源节点巧l旧标节点t的随 机游走过程中经过节点i的次数；
[0031] 所述的子图中屯、性系数的计算方法为：
[0032]
[0033] 化 a) = (Ak)ii
[0034] 其中，yk(i)表示W备选节点i为起点经过k个连边回到备选节点i的路径数目， (Ak)ii表示目标网络的邻接矩阵A的k次幕的第i个对角线元素；
[0035] 所述的特征向量中屯、性系数的计算方法为：
[0036] 设网络拓扑中具有η个节点，且该目标网络的邻接矩阵为A的特征值为λι，? = 1， 2，…，η，λ?对应的特征向量为a = ei，A的最大特征值为λ，则节点i的特征向量中屯、性系数的 计算公式为：
[0037]
[003引其中，aij为邻接矩阵A的第i行第巧阮素。
[0039] 进一步的，所述步骤(3)中，计算各个备选测点的Borda数的方法为：
[0040] (S3.1)、根据中屯、性系数矩阵中每列中屯、性数值大小进行节点排序，得到备选测 点的中屯、性系数优先关系矩阵M(mab)nX7,其中a，b分别为矩阵M(mab)nX7的行列编号；
[0041] (S3.2)、根据M(mab )nx7计算各个备选测点的波达数：
[0042] 设网络拓扑中某个测点ueU，U表示网络拓扑测点集合，Bi'（u)表示中屯、性系数优 先关系矩阵M(mab)nX满i '列中排在U之后的元素个数；
[0043] 若U在第i'列中位于第j'个，则有Bi'（u)=n-j'，称
为测点U的 Borda 数；
[0044] 其中，η表示中屯、性系数矩阵的行数，即备选测点数，最后将U中的各个元素按照 Borda数大小降序排列，得到U的一个综合排序；
[0045] 计算各个备选测点的隶属函数值的方法为：
[0046] (T3.1)、根据中屯、性系数矩阵，对每个节点进行中屯、性求和并排序，得到测点的模 糊优先关系矩阵R= (rij)nXn;
[0047] (Τ3.2)、利用平均法计算公式，得到各个测点的重要性排序
[004引设论域11=山，叫一，站}，4臣巧巧是一模糊集，3=(。^)。乂。为模糊优先关系矩 阵，模糊优先关系排序决策中的平均法如下：
[0049]
[0050] 如V,)即为各个备选测点的隶属函数值。
[0051] 更进一步的，所述的步骤(4)中，基于节点中屯、性的测点选择算法获取最优测点集 Snpt的具体步骤为：
[0052] (4.1)、将最优测点集Sopt置空;将所有备选测点放置于备选测点集S。进入(4.2);
[0053] (4.2)、从备选测点集Sc选取节点重要性集合{cv，v=l，2，…p}中评估值为fflfx柄} 对应的测点nv放入Sopt中，其中，cv为备选测点的中屯、性指标，体现为Borda数和隶属函数值， P为备选测点的总个数；
[0054] (4.3)、从Sdp冲随机移出一个nv，将最优测点集标记为5；,,，进入步骤(4.4);
[0055] (4.4)、将添加1^化。。*与移除1^的馬,进行对比，看移除1^勺馬,是否能隔离出新的 故障，如果不能，则进入步骤(4.3);如果能，则判断移除nv的Χ?,是否可W分隔出所有的模 拟电路系统故障，如果可W，进入(4.5)，否则跳转到步骤(4.3)，直到移除1^的^可^分隔 出所有的模拟电路系统故障，再进入步骤(4.5);
[0化6] (4.5)、将集合5；,中的备选测点按照按评估值大小从高到低排序，然后移除一个 具有最低评估值min k', ！的备选测点nv，再进入步骤(4.6);
[0057] (4.6)、检测剩下的测点集是否可W分隔所有故障，若可W，则删除Χ?,中对应移除 的备选测点ην，再跳转至步骤(4.5);若不能，则将该移除的备选测点ην放置于最优测点集 Sdp冲，同时将5；,-而置为馬在跳转至步骤(4.5)，直到所有的备选测点检查完毕，从而得 到最优测点集Snpt。
[005引本发明的发明目的是运样实现的：
[0059] 本发明基于网络拓扑结构的故障诊断测点选择方法，借用二分网络拓扑模型对模 拟电路的研究优势，计算复杂网络的节点中屯、性，并且采用Borda数法和平均法对节点的重 要性进行分 别的综合排序，提供一种基于节点中屯、性的测点选择算法得到最优测点集，最 后选择两种模糊决策方法中测点数目更少的方法作为最终的选择结果，具有很好的实践价 值，更够节省电路诊断的时间人力成本。
[0060] 同时，本发明基于网络拓扑结构的故障诊断测点选择方法还具有W下有益效果：
[0061] (1)、通过二分网络拓扑结构作为基础，可区分电路同构的情况，也能得到除元器 件网络W外的测点网络；
[0062] (2)、采用屯种网络节点中屯、性判断方法和两种模糊决策方法，充分考虑了不同评 估角度的重要性测度差别，得到的最终排序较为准确地反映真实情况；
[0063] (3)、基于节点中屯、性的测点选择算法是包含法和排除法的综合，能够确保得到的 最优测点集不仅能检测出所有的故障还尽可能包含最少的测点，能够很好地提高电路诊断 的效率。
【附图说明】
[0064] 图1带通滤波器电路的电路原理图；
[0065] 图2基于节点中屯、性系数的测点选择算法流程图；
[0066] 图3是最优测点集的获取流程图。
【具体实施方式】
[0067] 下面结合附图对本发明的【具体实施方式】进行描述，W便本领域的技术人员更好地 理解本发明。需要特别提醒注意的是，在W下的描述中，当已知功能和设计的详细描述也许 会淡化本发明的主要内容时，运些描述在运里将被忽略。
[006引实施例
[0069] 图带通滤波器电路为例。测点nO是电源测点，在下面的测点选择过程中，主要 研究测点nl~nil。下面应用Borda数法与平均法对带通滤波器电路进行测点网络的测点重 要性进行研究。
[0070] 图2基于节点中屯、性系数的测点选择算法流程图。
[0071] -种基于网络拓扑结构的故障诊断测点选择方法，如图2所示，包括W下步骤：
[0072] S1、计算网络拓扑中各个备选测点的中屯、性系数；
[0073] 中屯、性系数包括:节点度、介数、凝聚度、网络流、随机行走、子图和特征向量屯种 中屯、性系数；
[0074] 其中，节点度中屯、性系数的计算方法为：
[0075] 节点度中屯、性系数判断方法是根据网络中节点的度来判断节点的重要程度，因 此，设备选测点i的度是k( i)，则有：
[0076] Cd(i)=k(i)；
[0077] 介数中屯、性系数的计算方法为：
[0078] 介数中屯、性判断方法是根据网络中节点的介数中屯、性来衡量节点的重要程度，若 备选测点i的介数中屯、性较大，则表示该节点虽然具有较少的连边，但对整个网络的连通性 具有十分重要的作用。其计算公式为：
[0079]
[0080] 其中，s<t表示Ws为源节点，t为目标节点，ns康示源节点S和目标节点t中存在所 有最短路径的路径数，η表示网络中的节点数，gst,i表示经过备选测点i的节点对(s，t)的最 短路径数目；
[0081] 凝聚度中屯、性系数的计算方法为：
[0082] 凝聚度中屯、性判断方法是指根据网络中的节点与其他节点的流通程度，它表示备 选测点i到其他节点最短路径长度的之和。其计算公式为：
[0083]
[0084] 其中，dst表示备选测点i到网络中其他所有节点的最短距离；
[0085] 网络流中屯、性系数的计算方法为：
[0086] 网络流中屯、性判断方法是根据网络的几何中屯、来评估网络中的节点的重要性，其 计算公式为：
[0087]
[0088] 其中，gst表示网络中节点对(s，t)之间的所有路径数，但不包含回路；
[0089] 随机行走中屯、性系数的计算方法为：
[0090] 随机行走中屯、性判断方法是根据随机游走的思想来评价网络中节点的重要性程 度。下面给出由Newman提出的随机行走算法进行具体计算：
[0091] (a)、计算Laplace矩阵L，L = D-A，其中，D表示节点度组成的对角矩阵，A表示目标 网络的邻接矩阵；
[0092] 化）、除去L的最后一行、最后一列，得到矩阵L%
[009引（C)、计算矩阵L嘴逆矩阵L^，再对逆矩阵增添一个行向量和一个列向量，且 添加的向量分量均为0，得到矩阵T;
[0094] (d)、得到随机行走中屯、性系数的计算公式：
[0095]
[0096] 其中，j表示目标网络测点编号，Αυ表示目标网络邻接矩阵的第i行第j列元素， Tls、Tlt、Ls、Tパ的下标表示其所在矩阵T的行列位置，Ist,康示从源节点巧l旧标节点t的随 机游走过程中经过节点i的次数；
[0097] 子图中屯、性系数的计算方法为：
[0098] 子图中屯、性判断方法是指由网络的各种子图中某个节点的参与水平来判断该节 点的在网络中的重要性，其计算公式如下：
[0099]
[0100] 化(i) = (Ak)ii
[0101] 其中，yk(i)表示W备选节点i为起点经过k个连边回到备选节点i的路径数目， (Ak)ii表示目标网络的邻接矩阵A的k次幕的第i个对角线元素；
[0102] 特征向量中屯、性系数的计算方法为：
[0103] 特征向量中屯、性判断方法是根据目标网络的邻接矩阵A的特征值来衡量一个节点 的重要程度，它表示其他节点对该节点的影响。
[0104] 设网络拓扑中具有η个节点，且该目标网络的邻接矩阵为A的特征值为λι，? = 1， 2，··· ,η,λ?对应的特征向量为a = ei，
[0105] A的最大特征值为λ，则节点i的特征向量中屯、性系数的计算公式为：
[0106]
[0107]其中，ai功邻接矩阵A的第i行第巧阮素。
[010引综上，根据上述屯种节点中屯、性系数的计算，可W进一步得到屯种拓扑节点中屯、 性算法的性质和实际应用背景。
[0109] 表1是巧巾基于拓扑节点中屯、性算法的性质与实际应用背景；
[0110]
[0111] 表1
[0112] S2、将中屯、性系数归一化到0-1，得到归一化的中屯、性系数矩阵；
[0113] S3、根据中屯、性系数系数矩阵，利用Borda数法和平均法计算各个备选测点的 Borda数和隶属函数值，得到备选测点的两种重要性排序表；
[0114] 计算各个备选测点的Borda数的方法为：
[0115] S3.1、根据中屯、性系数矩阵中每列中屯、性数值大小进行节点排序，得到备选测点 的中屯、性系数优先关系矩阵M(mab)nX7,其中a，b分别为矩阵M(mab)nX7的行列编号；
[0116] S3.2、计算各个备选测点的波达数：
[0117] 设网络拓扑中某个测点uEU，U表示网络拓扑测点集合，Bi'（u)表示中屯、性系数优 先关系矩阵M(mab)nX满i '列中排在U之后的元素个数；
[011引若U在第i'列中位于第j'个，贝贿Bi'（u)=n-j'，称
为测点u的Borda 数；
[0119] 其中，η表示中屯、性系数矩阵的行数，即备选测点数，最后将U中的各个元素按照 Borda数大小降序排列，得到U的一个综合排序；
[0120] 计算各个备选测点的隶属函数值的方法为：
[0121] T3.1、根据中屯、性系数矩阵，对每个节点进行中屯、性求和并排序，得到测点的模糊 优先关系矩阵R=(ru)nXn;
[0122] T3.2、利用平均法计算公式，得到各个测点的重要性排序
[0123] 设论域U= {xi，X2,…，Xn}，d e /^(心')是一模糊集，R= (rij)nxn为模糊优先关系矩 阵，模糊优先关系排序决策中的平均法如下：
[0124]
[0125] ^(.τ)即为各个备选测点的隶属函数值。
[01%]在本实施例中，采用Borda数法与平均法得到排序结果如表2所示。
[0127]表2是带通滤波器电路的测点重要性排序；
[012 引
[0129]
[0130] 表2
[0131] 由表2可知，采用平均法与Borda数法得到的测点重要性排序结果比较相似，但也 存在差异性，主要是因为两种方法的评价角度不同。
[0132] S4、根据备选测点的两种重要性排序表，运行基于节点中屯、性的测点选择算法得 到最优测点集Sopt，选择测点数更少的Sopt作为最终选择结果。
[0133] 如图3所示，其具体步骤为：
[0134] S4.1、将最优测点集Snpt置空;将所有备选测点放置于备选测点集S。进入S4.2;
[0135] S4.2、从备选测点集S。选取节点重要性集合レv，v=l，2，…p}中评估 值为mfxK} 对应的测点nv放入Sopt中，其中，cv为备选测点的中屯、性指标，体现为Borda数和隶属函数值， P为备选测点的总个数；
[0136] S4.3、从Sopt中随机移出一个nv，将最优测点集标记为忠口，进入步骤S4.4;
[0137] S4.4、将添加 nv的Sopt与移除1^的&，进行对比，看移除nv的是否能隔离出新的 故障，如果不能，则进入步骤S4.3;如果能，则判断移除nv的是否可W分隔出所有的模拟 电路系统故障，如果可W，进入S4.5,否则跳转到步骤S4.3,直到移除nv的写W可W分隔出所 有的模拟电路系统故障，再进入步骤S4.5;
[0138] S4.5、将集合Χ?,.中的备选测点按照按评估值大小从高到低排序，然后移除一个具 有最低评估值的备选测点ην，再进入步骤S4.6;
[0139] S4.6、检测剩下的测点集是否可W分隔所有故障，若可W，则删除式Ρ,中对应移除 的备选测点nv，再跳转至步骤S4.5;若不能，则将该移除的备选测点nv放置于最优测点集Sopt 中，同时将5；,-而置为馬《在跳转至步骤S4.5，直到所有的备选测点检查完毕，从而得到最 优测点集Snpt。
[0140] 在本实施例中，根据表2可知，当测点排序靠前时，表明该测点相对于排序靠后的 ii点被选入最优测点集Snpt的优先级较高；相反，排序靠后的表明被选入最优测点集Snpt的 优先级较低。
[0141] 实例
[0142] 在本实施例中，根据上述方法，计算图1所示的带通滤波器电路的两种测点排序方 法，可W得到该带通滤波器电路的测点选择表，其中，基于节点中屯、性系数的测点选择结果 如表3所示：
[0143] 表3是基于节点中屯、性系数的测点选择结果。
[0144]
[0145] 表3
[0146] 由结果可W看到，本发明提出的基于节点中屯、性的测点选择算法是有效的，且运 行结果显示平均法的测点选择更佳。于是最终的选择结果为测点nl，n5，n8，n9，nl0，nll。
[0147] 尽管上面对本发明说明性的【具体实施方式】进行了描述，W便于本技术领域的技术 人员理解本发明，但应该清楚，本发明不限于【具体实施方式】的范围，对本技术领域的普通技 术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，运些 变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。
【主权项】
1. 一种基于网络拓扑结构的故障诊断测点选择方法，其特征在于，包括以下步骤： (1 )、计算网络拓扑中各个备选测点的中心性系数； (2) 、将中心性系数归一化到0-1，得到归一化的中心性系数矩阵； (3) 、根据中心性系数系数矩阵，利用Borda数法和平均法计算各个备选测点的Borda数 和隶属函数值，得到备选测点的两种重要性排序表； (4) 、根据备选测点的两种重要性排序表，运行基于节点中心性的测点选择算法得到最 优测点集Scipt，选择测点数更少的S cipt作为最终选择结果。2. 根据权利要求1所述的基于网络拓扑结构的故障诊断测点选择方法，其特征在于，所 述的中心性系数包括:节点度、介数、凝聚度、网络流、随机行走、子图和特征向量； 所述的节点度中心性系数的计算方法为： 设备选测点i的度是k(i)，则有： Cd(i)=k(i)； 所述的介数中心性系数的计算方法为：其中，s<t表示以s为源节点，t为目标节点，nst表示源节点s和目标节点t中存在所有最 短路径的路径数，η表示网络中的节点数，gst>1表示经过备选测点i的节点对(s，t)的最短路 径数目； 所述的凝聚度中心性系数的计算方法为：其中，Clst表示备选测点i到网络中其他所有节点的最短距离； 所述的网络流中心性系数的计算方法为：其中，gst表示网络中节点对(s，t)之间的所有路径数，但不包含回路。 所述的随机行走中心性系数的计算方法为： (a) 、计算Laplace矩阵L，L = D_A，其中，D表示节点度组成的对角矩阵，A表示目标网络 的邻接矩阵； (b) 、除去L的最后一行、最后一列，得到矩阵ΙΛ (c) 、计算矩阵C的逆矩阵广1，再对逆矩阵广1增添一个行向量和一个列向量，且添加的 向量分量均为〇,得到矩阵T; (d) 、得到随机行走中心性系数的计算公式：其中，j表示目标网络测点编号，^表示目标网络邻接矩阵的第i行第j列元素，Tls、Tlt、 TJS、Lt的下标表示其所在矩阵T的行列位置，Ist>1表示从源节点S到目标节点t的随机游走过 程中经过节点i的次数； 所述的子图中心性系数的计算方法为：yk(i) = (Ak) ? 其中，yk(i)表示以备选节点i为起点经过k个连边回到备选节点i的路径数目，以"^表 示目标网络的邻接矩阵A的k次幂的第i个对角线元素； 所述的特征向量中心性系数的计算方法为： 设网络拓扑中具有η个节点，且该目标网络的邻接矩阵为A的特征值为M，i = l，2,…，n， M对应的特征向量为a = ei，A的最大特征值为λ，则节点i的特征向量中心性系数的计算公式 为：其中，m沩邻接矩阵A的第i行第j列元素。3. 根据权利要求1所述的基于网络拓扑结构的故障诊断测点选择方法，其特征在于，所 述的步骤(3)中，计算各个备选测点的Borda数的方法为： (53.1) 、根据中心性系数矩阵中每列中心性数值大小进行节点排序，得到备选测点的 中心性系数优先先关系矩阵M(mab) nx7,其中a，b分别为矩阵M(mab)nx7的行列编号； (53.2) 、根据M(mab) nX7计算各个备选测点的波达数： 设网络拓扑中某个测点u e U，U表示网络拓扑测点集合，Bi，（ u)表示中心性系数优先关 系矩阵M(Hiab)nx7第i '列中排在u之后的元素个数； 若u在第i '列中位于第j '个，则有Bi1 (u) = n-j '，称1为测点U的Borda数； 其中，η表示中心性系数矩阵的行数，即备选测点数，最后将U中的各个元素按照Borda 数大小降序排列，得到U的一个综合排序； 计算各个备选测点的隶属函数值的方法为： (Τ3.1)、根据中心性系数矩阵，对每个节点进行中心性求和并排序，得到测点的模糊优 先关系矩阵R= (rij)nxn; (T3.2)、利用平均法计算公式，得到各个测点的重要性排序 设论域U= {xi，X2，···，Xn}，么GF(U)是一模模糊集，R=(rij)nXn为模糊优先关系矩阵，模 糊优先关系排序决策中的平均7去如下：A(X1)即为各个备选测点的隶属函数值。4. 根据权利要求1所述的基于网络拓扑结构的故障诊断测点选择方法，其特征在于，所 述的步骤(3)中，基于节点中心性的测点选择算法获取最优测点集S cipt的具体步骤为： (4.1) 、将最优测点集Scipt置空;将所有备选测点放置于备选测点集S。进入(4.2); (4.2) 、从备选测点集S。选取节点重要性集合{Cv，V = 1，2，…p}中评估值为X k丨对应 的测点nv放入Scipt中，其中，Cv为备选测点的中心性指标，体现为Borda数和隶属函数值，p为 备选测点的总个数； (4.3) 、从Scip冲随机移出一个nv，将最优测点集标记为进入步骤(4.4); (4.4) 、将添加nv的Scipt与移除nv的进行比，看移除nv的是否能隔离出新的故障， 如果不能，则进入步骤(4.3);如果能，则判断移除nv的5；?:是否可以分隔出所有的模拟电路 系统故障，如果可以，进入(4.5)，否则跳转到步骤(4.3)，直到移除ην的仏可以分隔出所有 的模拟电路系统故障，再进入步骤(4.5); (4.5) 、将集合^?中的备选测点按照按评估值大小从高到低排序，然后移除一个具有 最低评估值minkJ的备选测点η ν，再进入步骤(4.6); (4.6) 、检测剩下的测点集是否可以分隔所有故障，若可以，则删除中对应移除的备 选测点ην，再跳转至步骤(4.5);若不能，则将该移除的备选测点ην放置于最优测点集S cip冲， 同时将*5；,-置为^在跳转至步骤(4.5)，直到所有的备选测点检查完毕，从而得到最优 测点集S cipt。
【专利摘要】本发明公开了一种基于网络拓扑结构的故障诊断测点选择方法，借用二分网络拓扑模型对模拟电路的研究优势，计算复杂网络的节点中心性，并且采用Borda数法和平均法对节点的重要性进行分别的综合排序，提供一种基于节点中心性的测点选择算法得到最优测点集，最后选择两种模糊决策方法中测点数目更少的方法作为最终的选择结果，具有很好的实践价值，更够节省电路诊断的时间人力成本。
【IPC分类】H04L12/24
【公开号】CN105490840
【申请号】CN201510834935
【发明人】刘震, 付家敏, 杨成林, 田书林, 龙伊雯
【申请人】电子科技大学
【公开日】2016年4月13日
【申请日】2015年11月26日
【技术领域】
[0001] 本发明属于电子电路系统的可测性设计领域，应用二分复杂网络模型，设及一种 基于网络拓扑结构的故障诊断测点选择方法。
【背景技术】
[0002] 模拟电路系统故障测试与诊断中测点选择是其中一个重要研究课题。随着电路系 统复杂度与集成度的不断增加，模拟电路系统测点选择也越来越复杂。
[0003] 目前，复杂网络理论正是研究具有大规模，高复杂程度网络的一口交叉性和综合 性学科。但是到目前为止，复杂网络理论和电子电路系统的故障诊断、测点选择所结合，得 到的研究成果还较少，主要研究成果仅限于电子电路系统的小世界性。并且，大都采用传统 的网络建模方法，只能研究电路系统的单一网络模型:元器件网络模型或者测点网络模型， 运种方法丢失了原始电路系统的许多有用信息。而二分网络模型比传统的单顶点网络模型 在研究模拟电路系统上更有优势，不仅能够弥补单顶点网络模型不能够区分电路中拓扑同 构的缺陷，也能够得到除元器件网络W外的测点网络。
[0004] 在复杂网络中，网络节点中屯、性是用来衡量网络节点重要性的一个指标。本发明 立足拓扑结构分析网络节点中屯、性，主要有7种，分别是:节点度中屯、性判断方法、介数中屯、 性判断方法、凝聚度中屯、性判断方法、网络流中屯、性判断方法、随机行走中屯、性判断方法、 子图中屯、性判断方法W及特征向量中屯、性判断方法。方法概略为首先计算测点网络的各个 节点中屯、性，然后应用模糊数学中的一些模糊决策方法(模糊意见集中决策、模糊二元对比 决策等），得到网络中各个测点重要性排序，最后讨论网络中测点重要性排序与测点选择的 关系。
【发明内容】
[0005] 本发明的目的在于克服现有技术的不足，提供一种基于网络拓扑结构的故障诊断 测点选择方法，使城域网与接入网得到有效融合，同时具有可扩展、高灵活W及低成本的性 能。
[0006] 为实现上述发明目的，本发明基于网络拓扑结构的故障诊断测点选择方法，其特 征在于，包括W下步骤：
[0007] (1 )、计算网络拓扑中各个备选测点的中屯、性系数；
[000引（2)、将中屯、性系数归一化到0-1，得到归一化的中屯、性系数矩阵；
[0009] (3)、根据中屯、性系数系数矩阵，利用Borda数法和平均法计算各个备选测点的 Borda数和隶属函数值，得到备选测点的两种重要性排序表；
[0010] (4)、根据备选测点的两种重要性排序表，运行基于节点中屯、性的测点选择算法得 到最优测点集Snpt，选择测点数更少的Snpt作为最终选择结果。
[ocm]其中，所述的中屯、性系数包括:节点度、介数、凝聚度、网络流、随机行走、子图和特 征向量；
[0012] 所述的节点度中屯、性系数的计算方法为：
[0013] 设备选测点i的度是k( i)，则有：
[0014] Cd(i)=k(i)；
[0015] 所述的介数中屯、性系数的计算方法为：
[0016]
[0017] 其中，s<t表示Ws为源节点，t为目标节点，ns康示源节点S和目标节点帥存在所 有最短路径的路径数，η表示网络中的节点数，gst,i表示经过备选测点i的节点对(s，t)的最 短路径数目；
[0018] 所述的凝聚度中屯、性系数的计算方法为：
[0019]
[0020] 其中，dst表示备选测点i到网络中其他所有节点的最短距离；
[0021] 所述的网络流中屯、性系数的计算方法为：
[0022]
[0023] 其中，gst表示网络中节点对(s，t)之间的所有路径数，但不包含回路；
[0024] 所述的随机行走中屯、性系数的计算方法为：
[002引（a)、计算Laplace矩阵L，L = D-A，其中，D表示节点度组成的对角矩阵，A表示目标 网络的邻接矩阵；
[0026] 化）、除去L的最后一行、最后一列，得到矩阵L%
[0027] (C)、计算矩阵L嘴逆矩阵L^，再对逆矩阵增添一个行向量和一个列向量，且 添加的向量分量均为0，得到矩阵T;
[00%] (d)、得到随机行走中屯、性系数的计算公式：
[0029]
[0030] 其中，j表示目标网络测点编号，Αυ表示目标网络邻接矩阵的第i行第j列元素， Tls、Tlt、Ls、Tパ的下标表示其所在矩阵T的行列位置，Ist,康示从源节点巧l旧标节点t的随 机游走过程中经过节点i的次数；
[0031] 所述的子图中屯、性系数的计算方法为：
[0032]
[0033] 化 a) = (Ak)ii
[0034] 其中，yk(i)表示W备选节点i为起点经过k个连边回到备选节点i的路径数目， (Ak)ii表示目标网络的邻接矩阵A的k次幕的第i个对角线元素；
[0035] 所述的特征向量中屯、性系数的计算方法为：
[0036] 设网络拓扑中具有η个节点，且该目标网络的邻接矩阵为A的特征值为λι，? = 1， 2，…，η，λ?对应的特征向量为a = ei，A的最大特征值为λ，则节点i的特征向量中屯、性系数的 计算公式为：
[0037]
[003引其中，aij为邻接矩阵A的第i行第巧阮素。
[0039] 进一步的，所述步骤(3)中，计算各个备选测点的Borda数的方法为：
[0040] (S3.1)、根据中屯、性系数矩阵中每列中屯、性数值大小进行节点排序，得到备选测 点的中屯、性系数优先关系矩阵M(mab)nX7,其中a，b分别为矩阵M(mab)nX7的行列编号；
[0041] (S3.2)、根据M(mab )nx7计算各个备选测点的波达数：
[0042] 设网络拓扑中某个测点ueU，U表示网络拓扑测点集合，Bi'（u)表示中屯、性系数优 先关系矩阵M(mab)nX满i '列中排在U之后的元素个数；
[0043] 若U在第i'列中位于第j'个，则有Bi'（u)=n-j'，称
为测点U的 Borda 数；
[0044] 其中，η表示中屯、性系数矩阵的行数，即备选测点数，最后将U中的各个元素按照 Borda数大小降序排列，得到U的一个综合排序；
[0045] 计算各个备选测点的隶属函数值的方法为：
[0046] (T3.1)、根据中屯、性系数矩阵，对每个节点进行中屯、性求和并排序，得到测点的模 糊优先关系矩阵R= (rij)nXn;
[0047] (Τ3.2)、利用平均法计算公式，得到各个测点的重要性排序
[004引设论域11=山，叫一，站}，4臣巧巧是一模糊集，3=(。^)。乂。为模糊优先关系矩 阵，模糊优先关系排序决策中的平均法如下：
[0049]
[0050] 如V,)即为各个备选测点的隶属函数值。
[0051] 更进一步的，所述的步骤(4)中，基于节点中屯、性的测点选择算法获取最优测点集 Snpt的具体步骤为：
[0052] (4.1)、将最优测点集Sopt置空;将所有备选测点放置于备选测点集S。进入(4.2);
[0053] (4.2)、从备选测点集Sc选取节点重要性集合{cv，v=l，2，…p}中评估值为fflfx柄} 对应的测点nv放入Sopt中，其中，cv为备选测点的中屯、性指标，体现为Borda数和隶属函数值， P为备选测点的总个数；
[0054] (4.3)、从Sdp冲随机移出一个nv，将最优测点集标记为5；,,，进入步骤(4.4);
[0055] (4.4)、将添加1^化。。*与移除1^的馬,进行对比，看移除1^勺馬,是否能隔离出新的 故障，如果不能，则进入步骤(4.3);如果能，则判断移除nv的Χ?,是否可W分隔出所有的模 拟电路系统故障，如果可W，进入(4.5)，否则跳转到步骤(4.3)，直到移除1^的^可^分隔 出所有的模拟电路系统故障，再进入步骤(4.5);
[0化6] (4.5)、将集合5；,中的备选测点按照按评估值大小从高到低排序，然后移除一个 具有最低评估值min k', ！的备选测点nv，再进入步骤(4.6);
[0057] (4.6)、检测剩下的测点集是否可W分隔所有故障，若可W，则删除Χ?,中对应移除 的备选测点ην，再跳转至步骤(4.5);若不能，则将该移除的备选测点ην放置于最优测点集 Sdp冲，同时将5；,-而置为馬在跳转至步骤(4.5)，直到所有的备选测点检查完毕，从而得 到最优测点集Snpt。
[005引本发明的发明目的是运样实现的：
[0059] 本发明基于网络拓扑结构的故障诊断测点选择方法，借用二分网络拓扑模型对模 拟电路的研究优势，计算复杂网络的节点中屯、性，并且采用Borda数法和平均法对节点的重 要性进行分 别的综合排序，提供一种基于节点中屯、性的测点选择算法得到最优测点集，最 后选择两种模糊决策方法中测点数目更少的方法作为最终的选择结果，具有很好的实践价 值，更够节省电路诊断的时间人力成本。
[0060] 同时，本发明基于网络拓扑结构的故障诊断测点选择方法还具有W下有益效果：
[0061] (1)、通过二分网络拓扑结构作为基础，可区分电路同构的情况，也能得到除元器 件网络W外的测点网络；
[0062] (2)、采用屯种网络节点中屯、性判断方法和两种模糊决策方法，充分考虑了不同评 估角度的重要性测度差别，得到的最终排序较为准确地反映真实情况；
[0063] (3)、基于节点中屯、性的测点选择算法是包含法和排除法的综合，能够确保得到的 最优测点集不仅能检测出所有的故障还尽可能包含最少的测点，能够很好地提高电路诊断 的效率。
【附图说明】
[0064] 图1带通滤波器电路的电路原理图；
[0065] 图2基于节点中屯、性系数的测点选择算法流程图；
[0066] 图3是最优测点集的获取流程图。
【具体实施方式】
[0067] 下面结合附图对本发明的【具体实施方式】进行描述，W便本领域的技术人员更好地 理解本发明。需要特别提醒注意的是，在W下的描述中，当已知功能和设计的详细描述也许 会淡化本发明的主要内容时，运些描述在运里将被忽略。
[006引实施例
[0069] 图带通滤波器电路为例。测点nO是电源测点，在下面的测点选择过程中，主要 研究测点nl~nil。下面应用Borda数法与平均法对带通滤波器电路进行测点网络的测点重 要性进行研究。
[0070] 图2基于节点中屯、性系数的测点选择算法流程图。
[0071] -种基于网络拓扑结构的故障诊断测点选择方法，如图2所示，包括W下步骤：
[0072] S1、计算网络拓扑中各个备选测点的中屯、性系数；
[0073] 中屯、性系数包括:节点度、介数、凝聚度、网络流、随机行走、子图和特征向量屯种 中屯、性系数；
[0074] 其中，节点度中屯、性系数的计算方法为：
[0075] 节点度中屯、性系数判断方法是根据网络中节点的度来判断节点的重要程度，因 此，设备选测点i的度是k( i)，则有：
[0076] Cd(i)=k(i)；
[0077] 介数中屯、性系数的计算方法为：
[0078] 介数中屯、性判断方法是根据网络中节点的介数中屯、性来衡量节点的重要程度，若 备选测点i的介数中屯、性较大，则表示该节点虽然具有较少的连边，但对整个网络的连通性 具有十分重要的作用。其计算公式为：
[0079]
[0080] 其中，s<t表示Ws为源节点，t为目标节点，ns康示源节点S和目标节点t中存在所 有最短路径的路径数，η表示网络中的节点数，gst,i表示经过备选测点i的节点对(s，t)的最 短路径数目；
[0081] 凝聚度中屯、性系数的计算方法为：
[0082] 凝聚度中屯、性判断方法是指根据网络中的节点与其他节点的流通程度，它表示备 选测点i到其他节点最短路径长度的之和。其计算公式为：
[0083]
[0084] 其中，dst表示备选测点i到网络中其他所有节点的最短距离；
[0085] 网络流中屯、性系数的计算方法为：
[0086] 网络流中屯、性判断方法是根据网络的几何中屯、来评估网络中的节点的重要性，其 计算公式为：
[0087]
[0088] 其中，gst表示网络中节点对(s，t)之间的所有路径数，但不包含回路；
[0089] 随机行走中屯、性系数的计算方法为：
[0090] 随机行走中屯、性判断方法是根据随机游走的思想来评价网络中节点的重要性程 度。下面给出由Newman提出的随机行走算法进行具体计算：
[0091] (a)、计算Laplace矩阵L，L = D-A，其中，D表示节点度组成的对角矩阵，A表示目标 网络的邻接矩阵；
[0092] 化）、除去L的最后一行、最后一列，得到矩阵L%
[009引（C)、计算矩阵L嘴逆矩阵L^，再对逆矩阵增添一个行向量和一个列向量，且 添加的向量分量均为0，得到矩阵T;
[0094] (d)、得到随机行走中屯、性系数的计算公式：
[0095]
[0096] 其中，j表示目标网络测点编号，Αυ表示目标网络邻接矩阵的第i行第j列元素， Tls、Tlt、Ls、Tパ的下标表示其所在矩阵T的行列位置，Ist,康示从源节点巧l旧标节点t的随 机游走过程中经过节点i的次数；
[0097] 子图中屯、性系数的计算方法为：
[0098] 子图中屯、性判断方法是指由网络的各种子图中某个节点的参与水平来判断该节 点的在网络中的重要性，其计算公式如下：
[0099]
[0100] 化(i) = (Ak)ii
[0101] 其中，yk(i)表示W备选节点i为起点经过k个连边回到备选节点i的路径数目， (Ak)ii表示目标网络的邻接矩阵A的k次幕的第i个对角线元素；
[0102] 特征向量中屯、性系数的计算方法为：
[0103] 特征向量中屯、性判断方法是根据目标网络的邻接矩阵A的特征值来衡量一个节点 的重要程度，它表示其他节点对该节点的影响。
[0104] 设网络拓扑中具有η个节点，且该目标网络的邻接矩阵为A的特征值为λι，? = 1， 2，··· ,η,λ?对应的特征向量为a = ei，
[0105] A的最大特征值为λ，则节点i的特征向量中屯、性系数的计算公式为：
[0106]
[0107]其中，ai功邻接矩阵A的第i行第巧阮素。
[010引综上，根据上述屯种节点中屯、性系数的计算，可W进一步得到屯种拓扑节点中屯、 性算法的性质和实际应用背景。
[0109] 表1是巧巾基于拓扑节点中屯、性算法的性质与实际应用背景；
[0110]
[0111] 表1
[0112] S2、将中屯、性系数归一化到0-1，得到归一化的中屯、性系数矩阵；
[0113] S3、根据中屯、性系数系数矩阵，利用Borda数法和平均法计算各个备选测点的 Borda数和隶属函数值，得到备选测点的两种重要性排序表；
[0114] 计算各个备选测点的Borda数的方法为：
[0115] S3.1、根据中屯、性系数矩阵中每列中屯、性数值大小进行节点排序，得到备选测点 的中屯、性系数优先关系矩阵M(mab)nX7,其中a，b分别为矩阵M(mab)nX7的行列编号；
[0116] S3.2、计算各个备选测点的波达数：
[0117] 设网络拓扑中某个测点uEU，U表示网络拓扑测点集合，Bi'（u)表示中屯、性系数优 先关系矩阵M(mab)nX满i '列中排在U之后的元素个数；
[011引若U在第i'列中位于第j'个，贝贿Bi'（u)=n-j'，称
为测点u的Borda 数；
[0119] 其中，η表示中屯、性系数矩阵的行数，即备选测点数，最后将U中的各个元素按照 Borda数大小降序排列，得到U的一个综合排序；
[0120] 计算各个备选测点的隶属函数值的方法为：
[0121] T3.1、根据中屯、性系数矩阵，对每个节点进行中屯、性求和并排序，得到测点的模糊 优先关系矩阵R=(ru)nXn;
[0122] T3.2、利用平均法计算公式，得到各个测点的重要性排序
[0123] 设论域U= {xi，X2,…，Xn}，d e /^(心')是一模糊集，R= (rij)nxn为模糊优先关系矩 阵，模糊优先关系排序决策中的平均法如下：
[0124]
[0125] ^(.τ)即为各个备选测点的隶属函数值。
[01%]在本实施例中，采用Borda数法与平均法得到排序结果如表2所示。
[0127]表2是带通滤波器电路的测点重要性排序；
[012 引
[0129]
[0130] 表2
[0131] 由表2可知，采用平均法与Borda数法得到的测点重要性排序结果比较相似，但也 存在差异性，主要是因为两种方法的评价角度不同。
[0132] S4、根据备选测点的两种重要性排序表，运行基于节点中屯、性的测点选择算法得 到最优测点集Sopt，选择测点数更少的Sopt作为最终选择结果。
[0133] 如图3所示，其具体步骤为：
[0134] S4.1、将最优测点集Snpt置空;将所有备选测点放置于备选测点集S。进入S4.2;
[0135] S4.2、从备选测点集S。选取节点重要性集合レv，v=l，2，…p}中评估 值为mfxK} 对应的测点nv放入Sopt中，其中，cv为备选测点的中屯、性指标，体现为Borda数和隶属函数值， P为备选测点的总个数；
[0136] S4.3、从Sopt中随机移出一个nv，将最优测点集标记为忠口，进入步骤S4.4;
[0137] S4.4、将添加 nv的Sopt与移除1^的&，进行对比，看移除nv的是否能隔离出新的 故障，如果不能，则进入步骤S4.3;如果能，则判断移除nv的是否可W分隔出所有的模拟 电路系统故障，如果可W，进入S4.5,否则跳转到步骤S4.3,直到移除nv的写W可W分隔出所 有的模拟电路系统故障，再进入步骤S4.5;
[0138] S4.5、将集合Χ?,.中的备选测点按照按评估值大小从高到低排序，然后移除一个具 有最低评估值的备选测点ην，再进入步骤S4.6;
[0139] S4.6、检测剩下的测点集是否可W分隔所有故障，若可W，则删除式Ρ,中对应移除 的备选测点nv，再跳转至步骤S4.5;若不能，则将该移除的备选测点nv放置于最优测点集Sopt 中，同时将5；,-而置为馬《在跳转至步骤S4.5，直到所有的备选测点检查完毕，从而得到最 优测点集Snpt。
[0140] 在本实施例中，根据表2可知，当测点排序靠前时，表明该测点相对于排序靠后的 ii点被选入最优测点集Snpt的优先级较高；相反，排序靠后的表明被选入最优测点集Snpt的 优先级较低。
[0141] 实例
[0142] 在本实施例中，根据上述方法，计算图1所示的带通滤波器电路的两种测点排序方 法，可W得到该带通滤波器电路的测点选择表，其中，基于节点中屯、性系数的测点选择结果 如表3所示：
[0143] 表3是基于节点中屯、性系数的测点选择结果。
[0144]
[0145] 表3
[0146] 由结果可W看到，本发明提出的基于节点中屯、性的测点选择算法是有效的，且运 行结果显示平均法的测点选择更佳。于是最终的选择结果为测点nl，n5，n8，n9，nl0，nll。
[0147] 尽管上面对本发明说明性的【具体实施方式】进行了描述，W便于本技术领域的技术 人员理解本发明，但应该清楚，本发明不限于【具体实施方式】的范围，对本技术领域的普通技 术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，运些 变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。
【主权项】
1. 一种基于网络拓扑结构的故障诊断测点选择方法，其特征在于，包括以下步骤： (1 )、计算网络拓扑中各个备选测点的中心性系数； (2) 、将中心性系数归一化到0-1，得到归一化的中心性系数矩阵； (3) 、根据中心性系数系数矩阵，利用Borda数法和平均法计算各个备选测点的Borda数 和隶属函数值，得到备选测点的两种重要性排序表； (4) 、根据备选测点的两种重要性排序表，运行基于节点中心性的测点选择算法得到最 优测点集Scipt，选择测点数更少的S cipt作为最终选择结果。2. 根据权利要求1所述的基于网络拓扑结构的故障诊断测点选择方法，其特征在于，所 述的中心性系数包括:节点度、介数、凝聚度、网络流、随机行走、子图和特征向量； 所述的节点度中心性系数的计算方法为： 设备选测点i的度是k(i)，则有： Cd(i)=k(i)； 所述的介数中心性系数的计算方法为：其中，s<t表示以s为源节点，t为目标节点，nst表示源节点s和目标节点t中存在所有最 短路径的路径数，η表示网络中的节点数，gst>1表示经过备选测点i的节点对(s，t)的最短路 径数目； 所述的凝聚度中心性系数的计算方法为：其中，Clst表示备选测点i到网络中其他所有节点的最短距离； 所述的网络流中心性系数的计算方法为：其中，gst表示网络中节点对(s，t)之间的所有路径数，但不包含回路。 所述的随机行走中心性系数的计算方法为： (a) 、计算Laplace矩阵L，L = D_A，其中，D表示节点度组成的对角矩阵，A表示目标网络 的邻接矩阵； (b) 、除去L的最后一行、最后一列，得到矩阵ΙΛ (c) 、计算矩阵C的逆矩阵广1，再对逆矩阵广1增添一个行向量和一个列向量，且添加的 向量分量均为〇,得到矩阵T; (d) 、得到随机行走中心性系数的计算公式：其中，j表示目标网络测点编号，^表示目标网络邻接矩阵的第i行第j列元素，Tls、Tlt、 TJS、Lt的下标表示其所在矩阵T的行列位置，Ist>1表示从源节点S到目标节点t的随机游走过 程中经过节点i的次数； 所述的子图中心性系数的计算方法为：yk(i) = (Ak) ? 其中，yk(i)表示以备选节点i为起点经过k个连边回到备选节点i的路径数目，以"^表 示目标网络的邻接矩阵A的k次幂的第i个对角线元素； 所述的特征向量中心性系数的计算方法为： 设网络拓扑中具有η个节点，且该目标网络的邻接矩阵为A的特征值为M，i = l，2,…，n， M对应的特征向量为a = ei，A的最大特征值为λ，则节点i的特征向量中心性系数的计算公式 为：其中，m沩邻接矩阵A的第i行第j列元素。3. 根据权利要求1所述的基于网络拓扑结构的故障诊断测点选择方法，其特征在于，所 述的步骤(3)中，计算各个备选测点的Borda数的方法为： (53.1) 、根据中心性系数矩阵中每列中心性数值大小进行节点排序，得到备选测点的 中心性系数优先先关系矩阵M(mab) nx7,其中a，b分别为矩阵M(mab)nx7的行列编号； (53.2) 、根据M(mab) nX7计算各个备选测点的波达数： 设网络拓扑中某个测点u e U，U表示网络拓扑测点集合，Bi，（ u)表示中心性系数优先关 系矩阵M(Hiab)nx7第i '列中排在u之后的元素个数； 若u在第i '列中位于第j '个，则有Bi1 (u) = n-j '，称1为测点U的Borda数； 其中，η表示中心性系数矩阵的行数，即备选测点数，最后将U中的各个元素按照Borda 数大小降序排列，得到U的一个综合排序； 计算各个备选测点的隶属函数值的方法为： (Τ3.1)、根据中心性系数矩阵，对每个节点进行中心性求和并排序，得到测点的模糊优 先关系矩阵R= (rij)nxn; (T3.2)、利用平均法计算公式，得到各个测点的重要性排序 设论域U= {xi，X2，···，Xn}，么GF(U)是一模模糊集，R=(rij)nXn为模糊优先关系矩阵，模 糊优先关系排序决策中的平均7去如下：A(X1)即为各个备选测点的隶属函数值。4. 根据权利要求1所述的基于网络拓扑结构的故障诊断测点选择方法，其特征在于，所 述的步骤(3)中，基于节点中心性的测点选择算法获取最优测点集S cipt的具体步骤为： (4.1) 、将最优测点集Scipt置空;将所有备选测点放置于备选测点集S。进入(4.2); (4.2) 、从备选测点集S。选取节点重要性集合{Cv，V = 1，2，…p}中评估值为X k丨对应 的测点nv放入Scipt中，其中，Cv为备选测点的中心性指标，体现为Borda数和隶属函数值，p为 备选测点的总个数； (4.3) 、从Scip冲随机移出一个nv，将最优测点集标记为进入步骤(4.4); (4.4) 、将添加nv的Scipt与移除nv的进行比，看移除nv的是否能隔离出新的故障， 如果不能，则进入步骤(4.3);如果能，则判断移除nv的5；?:是否可以分隔出所有的模拟电路 系统故障，如果可以，进入(4.5)，否则跳转到步骤(4.3)，直到移除ην的仏可以分隔出所有 的模拟电路系统故障，再进入步骤(4.5); (4.5) 、将集合^?中的备选测点按照按评估值大小从高到低排序，然后移除一个具有 最低评估值minkJ的备选测点η ν，再进入步骤(4.6); (4.6) 、检测剩下的测点集是否可以分隔所有故障，若可以，则删除中对应移除的备 选测点ην，再跳转至步骤(4.5);若不能，则将该移除的备选测点ην放置于最优测点集S cip冲， 同时将*5；,-置为^在跳转至步骤(4.5)，直到所有的备选测点检查完毕，从而得到最优 测点集S cipt。
【专利摘要】本发明公开了一种基于网络拓扑结构的故障诊断测点选择方法，借用二分网络拓扑模型对模拟电路的研究优势，计算复杂网络的节点中心性，并且采用Borda数法和平均法对节点的重要性进行分别的综合排序，提供一种基于节点中心性的测点选择算法得到最优测点集，最后选择两种模糊决策方法中测点数目更少的方法作为最终的选择结果，具有很好的实践价值，更够节省电路诊断的时间人力成本。
【IPC分类】H04L12/24
【公开号】CN105490840
【申请号】CN201510834935
【发明人】刘震, 付家敏, 杨成林, 田书林, 龙伊雯
【申请人】电子科技大学
【公开日】2016年4月13日
【申请日】2015年11月26日
最新回复(0)