一种模拟电路软故障诊断的模糊推理方法

摘要

本发明涉及一种模拟电路软故障的诊断方法，属于模拟电路故障诊断技术领域。本方法首先构造软故障诊断的模糊规则，当软故障特征参数的在线监测值选中软故障诊断的模糊规则后，将被选中模糊规则前项的归一化置信度推理到后项，得到对模糊规则后项的归一化置信度。然后，对后项的置信度进行加权处理，得到对模拟电路软故障集合中每个软故障的置信度，依照置信度最大准则判断是何故障发生。本方法可以诊断已知软故障，还可以检测未知软故障，适用于较为复杂的模拟电路软故障诊断环境。根据本发明方法编制的程序(编译环境LabVIEW，C++等)可以在监控计算机上运行，并联合传感器、数据采集器等硬件组成在线监测系统，进行实时的模拟电路软故障的检测与诊断。

权利要求

1.一种模拟电路软故障的诊断方法，其特征在于该方法包括以下各步骤：（1）设定模拟电路的软故障集合为Θ={F1,…,Fj,…,FN,FUN}，其中Fj代表软故障集合Θ中的第j个已知软故障，j=1,2,…,N，N为已知软故障的个数，FUN为软故障集合Θ中的未知软故障；

（2）设模拟电路软故障特征参数集合E={e1,…,ei,…,eM}，其中ei代表特征参数集合E中的第i个软故障特征参数，i=1,2,…,M，M为软故障特征参数的个数；

（3）通过对模拟电路软故障的在线检测，确定上述软故障集合Θ中每个已知软故障发生时，软故障特征参数ei的取值区间ESi,j，取值区间集合为ESi={ESi,1,…,ESi,j,…,ESi,N}，表示在第j个已知软故障发生时第i个软故障特征参数的取值区间，分别为取值区间ESi,j的左、右端点值，上标L和R分别代表取值区间ESi,j的左、右端点；

（4）根据上述取值区间集合ESi，得到每个软故障特征参数ei的最大取值区间为其中 表示N个软故障特征参数取值区间ESi,1,…,ESi,j,…,ESi,N的左端点值中的最小值， 表示N个软故障特征参数取值区间

ESi,1,…,ESi,j,…,ESi,N的右端点中值的最大值，

（5）根据每个软故障特征参数ei的取值区间集合ESi和最大取值区间EIi，构造软故障诊断的模糊规则，具体过程如下：（5－1）对每个软故障特征参数ei的最大取值区间EIi进行三角均分，得到模糊语言项集 其中下标p代表模糊语言项的个数，共有Ji个模糊语言项，对于模糊语言项集Ui中的第2到第Ji-1个模糊语言项，即p=2,…,Ji-1，相应的模糊语言项为一个等腰三角型隶属度函数，其中 和 分别代表该等腰三角形的左端点、中点和右端点的取值，对于模糊语言项集Ui中的第1个和第Ji个模糊语言项，即p=1,Ji，相应的模糊语言项 是一个直角三角形隶属度函数，三角形的左端点和中点取值相同，即 与上述每个模糊语言项Ai,p对应的模糊区间记为 该区间的左端点、右端点的取值即为模糊语言项Ai,p的左端点、右端点的取值；

（5－2）判断第j个已知软故障的第i个软故障特征参数的取值区间ESi,j与模糊语言项集Ui中各个模糊语言项 之间的关系，若ESi,j∩IAi,p≠0，则第j个已知软故障的第i个软故障特征参数的取值区间ESi,j落入与模糊区间IAi,p对应的模糊语言项Ai,p中，若ESi,j∩IAi,p=0，则第j个已知软故障的第i个特征参数的取值区间ESi,j未落入与模糊区间IAi,p对应的模糊语言项Ai,p中；

（5－3）根据步骤（5－2）中的判断结果，得到三类共 条软故障诊断的模糊规则，第一类为单个已知软故障的模糊规则，共NR1条，第二类为已知软故障子集的模糊规则，共MNR条，第三类为未知软故障FUN的模糊规则，共UNR条，得到三类模糊规则的具体过程如下：

（5－3－1）得到第一类单个已知软故障的模糊规则：

若[且且…且]，则[模拟电路的软故障=Fj] （A）即若第j个已知软故障Fj的每个软故障特征参数ei的取值区间ES1,j,ES2,j,…,ESM,j分别依次落入模糊语言项A1,p1,A2,p2,…,AM,pM，则模拟电路出现已知软故障Fj，式（A）中的“[且且…且]”为模糊规则的“前项”，以组合形式A1,p1×A2,p2×…×AM,pM表示，“[模拟电路的软故障=Fj]”为模糊规则的“后项”，该后项表示模拟电路出现已知软故障Fj，若模糊语言项A1,p1,A2,p2,…,AM,pM分别为模糊语言项集U1,U2,…,UM中的n1,n2,…,nM个模糊语言项，则将n1,n2,…,nM个模糊语言项进行组合，得到NRj=n1n2…nM条软故障Fj的模糊规则，对于软故障集合Θ中的所有已知软故障，共计得到 条第一类单个已知软故障的模糊规则，其中包括前项和后项各不相同的NR1条模糊规则和前项相同、后项不同的NR2条模糊规则，即NR＝NR1+NR2，其中n1≤J1,n2≤J2,…,nM≤JM；

（5－3－2）得到第二类已知软故障子集 的模糊规则：

上述步骤（5－3－1）的NR2条已知单个软故障的模糊规则中，当出现模糊规则的前项相同，记为[且且…且]，但后项不同，即后项对应不同的单个已知软故障时，将相应模糊规则后项中不同的单个已知软故障合并成已知软故障子集 并将相应模糊规则合并，得到已知软故障子集的模糊规则为：

若[且且…且]，则[模拟电路的软故障

其中“[且且…且]”为模糊规则的“前项”，以组合形式A1,q1×A2,q2×…×AM,qM表示，“[模拟电路的软故障 ”为模糊规则的“后项”，该后项表示模拟电路出现的软故障是软故障子集 中的一

个；

当模糊语言项A1,q1,A2,q2,…,AM,qM分别为模糊语言项集U1,U2,…,UM中的m1,m2,…,mM个模糊语言项时，将得到MNR条已知软故障子集 的模糊规则，并有MNR≤m1m2…mM，其中m1≤n1,m2≤n2,…,mM≤nM；

（5－3－3）得到第三类未知软故障FUN的模糊规则：

若[且且…且]，则[模拟电路的软故障=FUN] （C）其中“[且且…且]”为模糊规则的“前项”，以组合形式A1,r1×A2,r2×…×AM,rM表示，模糊语言项A1,r1,A2,r2,…,AM,rM分别属于模糊语言项集U1,U2,…,UM，“[模拟电路的软故障=FUN]”为模糊规则的“后项”，该后项表示模拟电路出现已知软故障FUN，可生成UNR条未知软故障FUN的模糊规则，UNR＝TNR－MNR－NR1；

（6）当模拟电路在线运行时，对M个软故障特征参数进行观测，分别得到M个监测值xi，i=1,2,…,M，计算监测值xi隶属于模糊语言项集 中每个模糊语言项 的隶属度：

其中，xi为监测值，每个监测值xi分别会对Ui中编号连续的两个模糊语言项Ai,p和,Ai,p+1,p∈{1,2,…,Ji-1}的隶属度μi,p(xi)和μi,p+1(xi)大于零，对其他模糊语言项的隶属度等于零，则称xi选中了模糊语言项Ai,p和Ai,p+1，xi对模糊语言项Ai,p和Ai,p+1归一化后的置信度分别为：每个软故障特征参数的监测值xi选中两个模糊语言项Ai,p和Ai,p+1，对于M个软故障特征参数的监测值为x1,…,xi,…,xM，选中M对模糊语言项{A1,p,A1,p+1…,Ai,p,Ai,p+1,…M,AM,p,AM,p+1}，将M对模糊语言项进行组合得到共JNR=2 个软故障诊断的模糊规则的前项，则称M个软故障特征参数的监测值x1,…,xi,…,xM选中与JNR个前项所对应的软故障诊断的模糊规则，将其中第t条被选中的软故障诊断的模糊规则记为：

1,t 2,t M,t

若[且且…且]，则[模拟电路的软故障=Rt] （G）

1,t 2,t M,t

其中t=1,2,…,JNR，“[且且…且]”为被选中的软故障诊

1,t 2,t M,t 1,t

断的模糊规则的“前项”，以组合形式A ×A ×…×A 表示，并有A ∈{A1,p,A1,p+1}，

2,t M,t

A ∈{A2,p,A2,p+1}，…，A ∈{AM,p,AM,p+1},“[模拟电路的软故障=Rt]“为被选中的软故障诊断的模糊规则的“后项”，该后项中的Rt为被选中的单个已知软故障或已知软故障子集或未知软故障；

（7）根据式（E）和式（F）的被选中语言项置信度的计算方法，得到第t条被选中的软故障诊断的模糊规则前项的置信度为

1,t 2,t M,t 1,t 2,t M,t

m(A ×A ×…×A )=m(A )m(A )…m(A ) （H）则由式（H）推理出第t条被选中的软故障诊断的模糊规则后项的置信度为：

1,t 2,t M,t

ρ(Rt)=m(A ×A ×…×A ) （I）

并有 亦即该置信度是归一化的；

（8）对每条被选中的软故障诊断的模糊规则后项的置信度ρ(Rt)进行加权平均推理，得到软故障集合Θ中每个软故障的置信度ρ′(Fj)，j=1,2,…,N,UN，为：并有 即对每个软故障的置信度都是归一化，其中|Rt|代表Rt中软故障的个数，则式（J）表示把赋予整个Rt的置信度平均地分配给关于软故障集合Θ中的每个软故障，对于不包含在Rt中的单个软故障，其置信度为零；

（9）根据上述式（J）所获取的置信度进行故障决策，具体决策准则是：当模拟电路出现软故障特征参数监测值x1,…,xi,…,xM时，由式（J）获取的置信度值最大的软故障即为模拟电路此时发生的软故障。

说明书

一种模拟电路软故障诊断的模糊推理方法

技术领域

[0001] 本发明涉及一种模拟电路软故障诊断的模糊推理方法，属于模拟电路故障诊断技术领域。

背景技术

[0002] 在模拟电路系统的故障诊断中，由于故障的随机性和模糊性，使得故障与其特征之间的关系复杂，并存在不确定性。通常某一故障可以由多种特征或征兆反应，同时某一特
征也可能由多种故障引起的，所以难以借助单一的特征来诊断故障，也难以用精确的数学
模型来描述故障和特征之间多对多的映射关系。而模糊推理系统可以用“IF前项，THEN后
项”的产生式模糊规则表示和建模含有不确定性信息的诊断知识。该类模糊规则的前项中
包含多种故障特征，各特征模糊语言项之间的不同组合可对应后项的不同故障模式。实际
上，模糊推理系统利用模糊规则来形象地描述故障和特征之间的映射关系，并通过模糊推
理来模拟人类的推理及决策过程。所以模糊推理方法在实际的故障诊断中得到了广泛应
用。

[0003] 但是，由于软故障诊断的复杂性，在利用模糊推理方法进行软故障诊断时，还存在以下问题需要进一步研究：第一：模拟电路系统的软故障常常是由于电路中各部件电气参
数偏离允许的容差范围所造成的，与短路或断路硬故障相比，软故障特征参数的取值不再
是单点，可以用连续变化的区间来表示。在此情况下，如何将区间型软故障特征信息与软
故障之间的关系转化为产生式模糊规则规则，进而建立软故障模糊规则库是值得研究的问
题。第二：以往的模糊推理系统只能推理出已知故障，但是由于模拟电路软故障的多样性，
也需要模糊推理系统除了能够诊断已知软故障，还要能检测出未知软故障，亦即它的软故
障特征参数取值明显异于已知软故障的取值。第三：对于原有的故障诊断模糊推理系统，
在故障特征参数的在线监测值输入系统后，经推理后输出结果为各个已知故障发生的置信
度。此时，若故障特征参数的在线监测值对应未知故障，则系统输出结果未能准确反映该故
障，无法准确做出进一步的判断。所以针对以上问题，需要设计一个新的模糊推理系统，满
足模拟电路软故障诊断的需要。

发明内容

[0004] 本发明的目的是提出一种模拟电路软故障诊断的模糊推理方法，利用软故障特征参数的取值区间得到软故障诊断的模糊规则，运用归一化映射原理从模糊规则中推理出诊
断结果。该方法不仅可以诊断已知软故障，还可以检测未知软故障。

[0005] 本发明提出的一种模拟电路软故障诊断的模糊推理方法，包括以下各步骤：

[0006] (1)设定模拟电路的软故障集合为Θ＝{F1，…，Fj，…，FN，FUN}，其中Fj代表软故障集合Θ中的第j个已知软故障，j＝1，2，…，N，N为已知软故障的个数，FUN为软故障
集合Θ中的未知软故障；

[0007] (2)设模拟电路软故障特征参数集合E＝{e1，…，ei，…，eM}，其中ei代表特征参数集合E中的第i个软故障特征参数，i＝1，2，…，M，M为软故障特征参数的个数；

[0008] (3)通过对模拟电路软故障的在线检测，确定上述软故障集合Θ中每个已知软故障发生时，软故障特征参数ei的取值区间ESi，j，取值区间集合为ESi＝{ESi，1，…，ESi，
j，…，ESi，N}， 表示在第j个已知软故障发生时第i个软故障特征参数的取
值区间， 分别为取值区间ESi，j的左、右端点值，上标L和R分别代表取值区间ESi，j
的左、右端点；

[0009] (4)根据上述取值区间集合ESi，得到每个软故障特征参数ei的最大取值区间为其中 表示N个软故障特征参数取值区间ESi，1，…，ESi，j，…，ESi，N的左端
点值中的最小值， 表示N个软故障特征参数取值区间ESi，
1，…，ESi，j，…，ESi，N的右端点中值的最大值，

[0010] (5)根据每个软故障特征参数ei的取值区间集合ESi和最大取值区间EIi，构造软故障诊断的模糊规则，具体过程如下：

[0011] (5-1)对每个软故障特征参数ei的最大取值区间EIi进行三角均分，得到模糊语言项集Ui＝{Ai，1，…，Ai，p，…，Ai，Ji}，其中下标p代表模糊语言项的个数，共有Ji个模糊语
言项，对于模糊语言项集Ui中的第2到第Ji-1个模糊语言项，即p＝2，…，Ji-1，相应的模
糊语言项 为一个等腰三角型隶属度函数，其中 分别代表
该等腰三角形的左端点、中点和右端点的取值，对于模糊语言项集Ui中的第1个和第Ji个
模糊语言项，即p＝1，Ji，相应的模糊语言项 是一个直角三角形隶属度函
数，三角形的左端点和中点取值相同，即 与上述每个模糊语言项Ai，p对应的模糊区
间记为 该区间的左端点、右端点的取值即为模糊语言项Ai，p的左端点、右端
点的取值；

[0012] (5-2)判断第j个已知软故障的第i个软故障特征参数的取值区间ESi，j与模糊语言项集Ui中各个模糊语言项Ai，1，…，Ai，p，…，Ai，Ji之间的关系，若ESi，j∩IAi，p≠0，则第
j个已知软故障的第i个软故障特征参数的取值区间ESi，j落入与模糊区间IAi，p对应的模
糊语言项Ai，p中，若ESi，j∩IAi，p＝0，则第j个已知软故障的第i个特征参数的取值区间
ESi，j未落入与模糊区间IAi，p对应的模糊语言项Ai，p中；

[0013] (5-3)根据步骤(5-2)中的判断结果，得到三类共 条软故障诊断的模糊规则，第一类为单个已知软故障的模糊规则，共NR1条，第二类为已知软故障子集
的模糊规则，共MNR条，第三类为未知软故障FUN的模糊规则，共UNR
条，得到三类模糊规则的具体过程如下：

[0014] (5-3-1)得到第一类单个已知软故障的模糊规则：

[0015] 若[且且…且]，则[模拟电路的软故障＝Fj] (1)即若第j个已知软故障Fj的每个软故障特征参数ei的取值区间ES1，j，ES2，
j，…，ESM，j分别依次落入模糊语言项A1，p1，A2，p2，…，AM，pM，则模拟电路出现已知软故障Fj，
式(1)中的“[且且…且]”为模糊规则的“前项”，以组
合形式A1，p1×A2，p2×…×AM，pM表示，“[模拟电路的软故障＝Fj]”为模糊规则的“后项”，该
后项表示模拟电路出现已知软故障Fj，若模糊语言项A1，p1，A2，p2，…，AM，pM分别为模糊语言
项集U1，U2，…，UM中的n1，n2，…，nM(n1≤J1，n2≤J2，…，nM≤JM)个模糊语言项，则将
n1，n2，…，nM个模糊语言项进行组合，得到NRj＝n1n2…nM条软故障Fj的模糊规则，对于软
故障集合Θ中的所有已知软故障，共计得到 条第一类单个已知软故障的模糊规
则，其中包括前项和后项各不相同的NR1条模糊规则和前项相同、后项不同的NR2条模糊规
则，即NR＝NR1+NR2；

[0016] (5-3-2)得到第二类已知软故障子集 的模糊规则：

[0017] 上述步骤(5-3-1)的NR2条已知单个软故障的模糊规则中，当出现模糊规则的前项相同，记为[且且…且]，但后项不同，即后项
对应不同的单个已知软故障时，将相应模糊规则后项中不同的单个已知软故障合并成
已知软故障子集 并将相应模糊规则合并，得到已知软故障子集
的模糊规则为：

[0018] 若[且且…且]，

[0019] 则

[0020] 其中“[且且…且]”为模糊规则的“前项”，以组合形式A1，q1×A2，q2×…×AM，qM表示， 为模
糊规则的“后项”，该后项表示模拟电路出现的软故障是软故障子集
中的一个；

[0021] 当模糊语言项A1，q1，A2，q2，…，AM，qM分别为模糊语言项集U1，U2，…，UM中的m1，m2，…，mM(m1≤n1，m2＜n2，…，mM≤nM)个模糊语言项时，将得到MNR条已知软故障子集
的模糊规则，并有MNR≤m1m2…mM；

[0022] (5-3-3)得到第三类未知软故障FUN的模糊规则：

[0023] 若[且且…且]，则[模拟电路的软故障＝FUN](3)其中“[且且…且]”为模糊规则的“前项”，以组合
形式A1，r1×A2，r2×…×AM，rM表示，模糊语言项A1，r1，A2，r2，…，AM，rM分别属于模糊语言项集
U1，U2，…，UM，“[模拟电路的软故障＝FUN]”为模糊规则的“后项”，该后项表示模拟电路出
现已知软故障FUN，可生成UNR条未知软故障FUN的模糊规则，UNR＝TNR-MNR-NR1；

[0024] (6)当模拟电路在线运行时，对M个软故障特征参数进行观测，分别得到M个监测值xi，i＝1，2，…，M，计算监测值xi隶属于模糊语言项集Ui＝{Ai，1，…，Ai，p，…，Ai，Ji}
中每个模糊语言项 的隶属度：

[0025]

[0026] 其中，xi为监测值，每个监测值xi分别会对Ui中编号连续的两个模糊语言项Ai，p和，Ai，p+1，(p∈{1，2，…，Ji-1})的隶属度μi，p(xi)和μi，p+1(xi)大于零，对其他模糊语言
项的隶属度等于零，则称xi选中了模糊语言项Ai，p和Ai，p+1，xi对模糊语言项Ai，p和Ai，p+1归
一化后的置信度分别为：

[0027]

[0028]

[0029] 每个软故障特征参数的监测值xi选中两个模糊语言项Ai，p和Ai，p+1，对于M个软故障特征参数的监测值为x1，…，xi，…，xM，选中M对模糊语言项{A1，p，A1，p+1…，Ai，p，Ai，
p+1，…，AM，p，AM，p+1}，将M对模糊语言项进行组合得到共JNR＝2M个软故障诊断的模糊规则
的前项，则称M个软故障特征参数的监测值x1，…，xi，…，xM选中与JNR个前项所对应的
软故障诊断的模糊规则，将其中第t(t＝1，2，…，JNR)条被选中的软故障诊断的模糊规则
记为：

[0030] 若[且且…且]，则[模拟电路的软故障＝Rt](7)其中“[且且…且]”为被选中的软故障诊断的模糊
规则的“前项”，以组合形式A1，t×A2，t×…×AM，t表示，并有A1，t∈{A1，p，A1，p+1}，A2，t∈{A2，
p，A2，p+1}，…，AM，t∈{AM，p，AM，p+1}，“[模拟电路的软故障＝Rt]“为被选中的软故障诊断的
模糊规则的“后项”，该后项中的Rt为被选中的单个已知软故障或已知软故障子集或未知软
故障；

[0031] (7)根据式(5)和式(6)的被选中语言项置信度的计算方法，得到第t条被选中的软故障诊断的模糊规则前项的置信度为

[0032] m(A1，t×A2，t×…×AM，t)＝m(A1，t)m(A2，t)…m(AM，t) (8)

[0033] 则由式(8)推理出第t条被选中的软故障诊断的模糊规则后项的置信度为：

[0034] ρ(Rt)＝m(A1，t×A2，t×…×AM，t) (9)

[0035] 并有 亦即该置信度是归一化的；

[0036] (8)对每条被选中的软故障诊断的模糊规则后项的置信度ρ(Rt)进行加权平均推理，得到软故障集合Θ中每个软故障的置信度ρ’(Fj)，j＝1，2，…，N，UN，为：

[0037]

[0038] 并有 即对每个软故障的置信度都是归一化，其中|Rt|代表Rt中软故障的个数，则式(10)表示把赋予整个Rt的置信度平均地分配给关于软故障集合Θ中的
每个软故障，对于不包含在Rt中的单个软故障，其置信度为零；

[0039] (9)根据上述式(10)所获取的置信度进行故障决策，具体决策准则是：当模拟电路出现软故障特征参数监测值x1，…，xi，…，xM时，由式(10)获取的置信度值最大的软故
障即为模拟电路此时发生的软故障。

[0040] 本发明涉及一种模拟电路软故障诊断的模糊推理方法，属于模拟电路故障诊断技术领域。基于已知软故障的故障特征参数的取值区间，首先构造软故障诊断的模糊规则，这
些模糊规则的后项可以指向单个已知软故障，也可以指向已知软故障子集，或未知软故障。
当软故障特征参数的在线监测值选中软故障诊断的模糊规则后，可将这些被选中模糊规则
前项的归一化置信度推理到后项，得到对模糊规则后项的归一化置信度。然后，对后项的置
信度进行加权处理，得到对模拟电路软故障集合中每个软故障的置信度，依照置信度最大
准则判断是何故障发生。利用所设计的模糊推理方法，不仅可以诊断已知软故障，还可以检
测未知软故障，适用于较为复杂的模拟电路软故障诊断环境。根据本发明方法编制的程序
(编译环境LabVIEW，C++等)可以在监控计算机上运行，并联合传感器、数据采集器等硬件
组成在线监测系统，进行实时的模拟电路软故障的检测与诊断。

附图说明

[0041] 图1是本发明方法的流程框图。

[0042] 图2是本发明具体实施方式中软故障特征参数1在3种已知软故障下的取值区间及模糊语言项。

[0043] 图3是本发明具体实施方式中软故障特征参数2在3种已知软故障下的取值区间及模糊语言项。

[0044] 图4是本发明具体实施方式中软故障特征参数3在3种已知软故障下的取值区间及模糊语言项。

[0045] 图5是本发明实施例中ZPW-2000无绝缘轨道电路结构图。

[0046] 图6是本发明实施例中软故障特征参数1在9种已知软故障下的取值区间。

[0047] 图7是本发明实施例中软故障特征参数2在9种已知软故障下的取值区间。

[0048] 图8是本发明实施例中软故障特征参数3在9种已知软故障下的取值区间。

[0049] 图9是本发明实施例中软故障特征参数4在9种已知软故障下的取值区间。

[0050] 图10是本发明实施例中软故障特征参数5在9种已知软故障下的取值区间。

具体实施方式

[0051] 本发明提出的一种模拟电路软故障的诊断的模糊推理方法，其流程框图如图1所示，包括以下各步骤：

[0052] (1)设定模拟电路的软故障集合为Θ＝{F1，…，Fj，…，FN，FUN}，其中Fj代表软故障集合Θ中的第j个已知软故障，j＝1，2，…，N，N为已知软故障的个数，FUN为软故障
集合Θ中的未知软故障；

[0053] (2)设模拟电路软故障特征参数集合E＝{e1，…，ei，…，eM}，其中ei代表特征参数集合E中的第i个软故障特征参数，i＝1，2，…，M，M为软故障特征参数的个数；

[0054] (3)通过对模拟电路软故障的在线检测，确定上述软故障集合Θ中每个已知软故障发生时，软故障特征参数ei的取值区间ESi，j，取值区间集合为ESi＝{ESi，1，…，ESi，
j，…，ESi，N}， 表示在第j个已知软故障发生时第i个软故障特征参数的取
值区间， 分别为取值区间ESi，j的左、右端点值，上标L和R分别代表取值区间ESi，j
的左、右端点；

[0055] (4)根据上述取值区间集合ESi，得到每个软故障特征参数ei的最大取值区间为其中 表示N个软故障特征参数取值区间ESi，1，…，ESi，j，…，ESi，N的左端
点值中的最小值， 表示N个软故障特征参数取值区间ESi，
1，…，ESi，j，…，ESi，N的右端点中值的最大值，

[0056] 为了加深对各步骤的理解，这里举例说明。设某模拟电路的软故障集合为Θ＝{F1，F2，F3，FUN}，即j＝1，2，3，N＝3，它们共同的软故障特征参数集合为E＝
{e1，e2，e3}，即i＝1，2，3，M＝3。每个软故障特征参数ei在每种软故障Fj下的取
值区间分别如图2、图3和图4所示，并根据步骤(3)得到取值区间集合ES1＝{ES1，
1，ES1，2，ES1，3}，ES2＝{ES2，1，ES2，2，ES2，3}，ES3＝{ES3，1，ES3，2，ES3，3}，各取值区间的端
点值在横轴上标出。根据步骤(4)确定三个软故障特征参数的最大取值区间分别为

[0057] (5)根据每个软故障特征参数ei的取值区间集合ESi和最大取值区间EIi，构造软故障诊断的模糊规则，具体过程如下：

[0058] (5-1)对每个软故障特征参数ei的最大取值区间EIi进行三角均分，得到模糊语言项集Ui＝{Ai，1，…，Ai，p，…，Ai，Ji}，其中下标p代表模糊语言项的个数，共有Ji个模糊语
言项，对于模糊语言项集Ui中的第2到第Ji-1个模糊语言项，即p＝2，…，Ji-1，相应的模
糊语言项 为一个等腰三角型隶属度函数，其中 分别代表
该等腰三角形的左端点、中点和右端点的取值，对于模糊语言项集Ui中的第1个和第Ji个
模糊语言项，即p＝1，Ji，相应的模糊语言项 是一个直角三角形隶属度函
数，三角形的左端点和中点取值相同，即 与上述每个模糊语言项Ai，p对应的模糊区
间记为 该区间的左端点、右端点的取值即为模糊语言项Ai，p的左端点、右端
点的取值；

[0059] 例如，在图2至图4所示例子中，根据步骤(5-1)对每个软故障特征参数ei的最大取值区间EIi进行三角均分，分别得到三个软故障特征参数的模糊语言项集U1＝{A1，1，A1，2，
A1，3，A1，4，A1，5}，U2＝{A2，1，A2，2，A2，3，A2，4，A2，5}，U3＝{A3，1，A3，2，A3，3，A3，4，A3，5}，每个模糊语言项集中，模糊语言项的个数分别为J1＝J2J3＝5。当p＝2，…，4时，
i＝1，2，3，为一等腰三角形隶属度函数，其中 和 分别代表该等腰三角形的左端
点、中点和右端点的取值；当p＝1，5时，相应的模糊语言项 是一个直角
三角形隶属度函数，三角形的左端点和中点取值相同，即 与上述每个模糊语言项
Ai，p对应的模糊区间记为 该区间的左端点、右端点的取值即为模糊语言项
Ai，p的左端点、右端点的取值；

[0060] (5-2)判断第j个已知软故障的第i个软故障特征参数的取值区间ESi，j与模糊语言项集Ui中各个模糊语言项Ai，1，…，Ai，p，…，Ai，Ji之间的关系，若ESi，j∩IAi，p≠0，则第
j个已知软故障的第i个软故障特征参数的取值区间ESi，j落入与模糊区间IAi，p对应的模
糊语言项Ai，p中，若ESi，j∩IAi，p＝0，则第j个已知软故障的第i个特征参数的取值区间
ESi，j未落入与模糊区间IAi，p对应的模糊语言项Ai，p中；

[0061] 例如，在图2中，判断第一个已知软故障F1的第一个特征特征ES1，1与模糊语言项集U1中模糊语言项A1，1，A1，2，A1，3，A1，4，A1，5之间的关系，因为ES1，1∩IA1，1≠0，所以ES1，1落
入与模糊区间IA1，1对应的模糊语言项A1，1中，ES1，1∩IA1，2≠0，所以ES1，1落入与模糊区间
IA1，2对应的模糊语言项A1，2中，此外，ES1，1与模糊区间IA1，3，IA1，4，IA1，5的交集都为零，所以
ES1，1未落入与IA1，3，IA1，4，IA1，5对应的模糊语言项A1，3，A1，4，A1，5中。

[0062] (5-3)根据步骤(5-2)中的判断结果，得到三类共 条软故障诊断的模糊规则，第一类为单个已知软故障的模糊规则，共NR1条，第二类为已知软故障子集
的模糊规则，共MNR条，第三类为未知软故障FUN的模糊规则，共UNR
条，得到三类模糊规则的具体过程如下：

[0063] 在图2、图3和图4所示例子中，含有三个软故障特征参数，它们的模糊语言项集分别为U1＝{A1，1，A1，2，A1，3，A1，4，A1，5}，U2＝{A2，1，A2，2，A2，3，A2，4，A2，5}，U3＝{A3，1，A3，2，A3，3，A3，
4，A3，5}，每个Ui，i＝1，2，3，中含有的模糊语言项的个数分别为J1＝J2＝J3＝5，则可以得
到共TNR＝J1J2J3＝125条软故障诊断的模糊规则。

[0064] (5-3-1)得到第一类单个已知软故障的模糊规则：

[0065] 若[且且…且]，则[模拟电路的软故障＝Fj](1)即若第j个已知软故障Fj的每个软故障特征参数ei的取值区间ES1，j，ES2，j，…，ESM，
j分别依次落入模糊语言项A1，p1，A2，p2，…，AM，pM，则模拟电路出现已知软故障Fj，式(1)中
的“[且且…且]”为模糊规则的“前项”，以组合形式
A1，p1×A2，p2×…×AM，pM表示，“[模拟电路的软故障＝Fj]”为模糊规则的“后项”，该后项表
示模拟电路出现已知软故障Fj，若模糊语言项A1，p1，A2，p2，…，AM，pM分别为模糊语言项集U1，
U2，…，UM中的n1，n2，…，nM(n1≤J1，n2≤J2，…，nM≤JM)个模糊语言项，则将n1，n2，…，
nM个模糊语言项进行组合，得到NRj＝n1n2…nM条软故障Fj的模糊规则，对于软故障集合
Θ中的所有已知软故障，共计得到 条第一类单个已知软故障的模糊规则，其中
包括前项和后项各不相同的NR1条模糊规则和前项相同、后项不同的NR2条模糊规则，即NR
＝NR1+NR2；

[0066] 在图2、图3和图4所示例子中，例如第2个已知软故障F2的三个软故障特征参数的取值区间ES1，2，ES2，2，ES3，2分别依次落入模糊语言项A1，2，A2，1，A3，4，则模拟电路出现已知
故障F2，此时已知故障F2的软故障诊断的模糊规则为：

[0067] 若[且且]，则[模拟电路的软故障＝F2]已知软故障F2的软故障特征参数取值区间ES1，2，ES2，2，ES3，2分别落入模糊语言项集U1的模糊语
言项A1，2，A1，3，A1，4中(n1＝3)，落入U2模糊语言项A2，1，A2，2，A2，3中(n2＝3)，落入U3的模
糊语言项A3，2，A3，2，A3，4中(n3＝3)，则将这三组模糊语言项进行组合，得到NR2＝3·3·3
＝9条单个软故障F2的模糊规则。同理，可以分别获得软故障F1和F3的NR1＝2·3·2＝
12条和NR3＝4·3·3＝36条软故障诊断的模糊规则，共计得到NR＝12+27+36＝75条
第一类单个已知软故障的模糊规则，其中包括前项和后项各不相同的NR1＝65条模糊规则
和前项相同、后项不同的NR2＝10条模糊规则，即NR＝NR1+NR2；

[0068] (5-3-2)得到第二类已知软故障子集 的模糊规则：

[0069] 上述步骤(5-3-1)的NR2条已知单个软故障的模糊规则中，当出现模糊规则的前项相同，记为[且且…且]，但后项不同，即后项
对应不同的单个已知软故障时，将相应模糊规则后项中不同的单个已知软故障合并成
已知软故障子集 并将相应模糊规则合并，得到已知软故障子集
的模糊规则为：

[0070] 若[且且…且]，

[0071] 则

[0072] 其中“[且且…且]”为模糊规则的“前项”，以组合形式A1，q1×A2，q2×…×AM，qM表示， 为模
糊规则的“后项”，该后项表示模拟电路出现的软故障是软故障子集
中的一个；

[0073] 当模糊语言项A1，q1，A2，q2，…，AM，qM分别为模糊语言项集U1，U2，…，UM中的m1，m2，…，mM(m1≤n1，m2≤n2，…，mM≤nM)个模糊语言项时，将得到MNR条已知软故障子集
的模糊规则，并有MNR≤m1m2…mM；

[0074] 在图2、图3和图4所示例子中，得到的NR2＝10条前项相同、后项不同的第一类单个软故障的模糊规则为：

[0075] a.若[且且]，则[模拟电路的软故障＝F1]

[0076] b.若[且且]，则[模拟电路的软故障＝F2]

[0077] c.若[且且]，则[模拟电路的软故障＝F2]

[0078] d.若[且且]，则[模拟电路的软故障＝F3]

[0079] e.若[且且]，则[模拟电路的软故障＝F2]

[0080] f.若[且且]，则[模拟电路的软故障＝F3]

[0081] g.若[且且]，则[模拟电路的软故障＝F2]

[0082] h.若[且且]，则[模拟电路的软故障＝F3]

[0083] i.若[且且]，则[模拟电路的软故障＝F2]

[0084] j.若[且且]，则[模拟电路的软故障＝F3]

[0085] 模糊规则a和b的前项相同，后项分别为F1和F2，根据步骤(5-3-2)可将两者合并为一条软故障子集{F1，F2}的模糊规则为：

[0086] 若 [< e1 ＝ A1，2> 且 且 ]， 则同理，分别将模糊规则c和d、e和f、g和h、i
和j进行合并，得到以下四条软故障子集{F2，F3}的模糊规则为：

[0087] 若 [< e1 ＝ A1，2> 且 且 ]， 则

[0088] 若 [< e1 ＝ A1，2> 且 且 ]， 则

[0089] 若 [< e1 ＝ A1，3> 且 且 ]， 则

[0090] 若 [< e1 ＝ A1，3> 且 且 ]， 则所以，共计得到MNR＝5条第二类已知软故障
子集的模糊规则，分别涉及模糊语言项集U1，U2，U3中的m1＝m2＝m3＝2个模糊语言项，并
有MNR≤m1＝m2＝m3＝8。

[0091] (5-3-3)得到第三类未知软故障FUN的模糊规则：

[0092] 若[且且…且]，则[模拟电路的软故障＝FUN](3)其中“[且且…且]”为模糊规则的“前项”，以组合
形式A1，r1×A2，r2×…×AM，rM表示，模糊语言项A1，r1，A2，r2，…，AM，rM分别属于模糊语言项集
U1，U2，…，UM，“[模拟电路的软故障＝FUN]”为模糊规则的“后项”，该后项表示模拟电路出
现已知软故障FUN，可生成UNR条未知软故障FUN的模糊规则，UNR＝TNR-MNR-NR1；

[0093] 在图2、图3和图4所示例子中，得到UNR＝TNR-MNR-NR1＝125-5-65＝55条第三类未知故障FUN的模糊规则，例如模糊规则：

[0094] 若[且且]，则[模拟电路的软故障＝FUN]；

[0095] (6)当模拟电路在线运行时，对M个软故障特征参数进行观测，分别得到M个监测值xi，i＝1，2，…，M，计算监测值xi隶属于模糊语言项集Ui＝{Ai，1，…，Ai，p，…，Ai，Ji}
中每个模糊语言项 的隶属度：

[0096]

[0097] 其中，xi为监测值，每个监测值xi分别会对Ui中编号连续的两个模糊语言项Ai，p和，Ai，p+1，(p∈{1，2，…，Ji-1})的隶属度μi，p(xi)和μi，p+1(xi)大于零，对其他模糊语言
项的隶属度等于零，则称xi选中了模糊语言项Ai，p和Ai，p+1，xi对模糊语言项Ai，p和Ai，p+1归
一化后的置信度分别为：

[0098]

[0099]

[0100] 每个软故障特征参数的监测值xi选中两个模糊语言项Ai，p和Ai，p+1，对于M个软故障特征参数的监测值为x1，…，xi，…，xM，选中M对模糊语言项{A1，p，A1，p+1…，Ai，p，Ai，
p+1，…，AM，p，AM，p+1}，将M对模糊语言项进行组合得到共JNR＝2M个软故障诊断的模糊规则
的前项，则称M个软故障特征参数的监测值x1，…，xi，…，xM选中与JNR个前项所对应的
软故障诊断的模糊规则，将其中第t(t＝1，2，…，JNR)条被选中的软故障诊断的模糊规则
记为：

[0101] 若[且且…且]，则[模拟电路的软故障＝Rt](7)其中“[且且…且]”为被选中的软故障诊断的模糊
规则的“前项”，以组合形式A1，t×A2，t×…×AM，t表示，并有A1，t∈{A1，p，A1，p+1}，A2，t∈{A2，
p，A2，p+1}，…，AM，t∈{AM，p，AM，p+1}，“[模拟电路的软故障＝Rt]“为被选中的软故障诊断的
模糊规则的“后项”，该后项中的Rt为被选中的单个已知软故障或已知软故障子集或未知软
故障；

[0102] (7)根据式(5)和式(6)的被选中语言项置信度的计算方法，得到第t条被选中的软故障诊断的模糊规则前项的置信度为

[0103] m(A1，t×A2，t×…×AM，t)＝m(A1，t)m(A2，t)…m(AM，t) (8)

[0104] 则由式(8)推理出第t条被选中的软故障诊断的模糊规则后项的置信度为：

[0105] ρ(Rt)＝m(A1，t×A2，t×…×AM，t) (9)

[0106] 并有 亦即该置信度是归一化的；

[0107] (8)对每条被选中的软故障诊断的模糊规则后项的置信度ρ(Rt)进行加权平均推理，得到软故障集合Θ中每个软故障的置信度ρ’(Fj)，j＝1，2，…，N，UN，为：

[0108]

[0109] 并有 即对每个软故障的置信度都是归一化，其中|Rt|代表Rt中软故障的个数，则式(10)表示把赋予整个Rt的置信度平均地分配给关于软故障集合Θ中的
每个软故障，对于不包含在Rt中的单个软故障，其置信度为零；

[0110] (9)根据上述式(10)所获取的置信度进行故障决策，具体决策准则是：当模拟电路出现软故障特征参数监测值x1，…，xi，…，xM时，由式(10)获取的置信度值最大的软故
障即为模拟电路此时发生的软故障。

[0111] 以下结合附图，详细介绍本发明方法的实施例：

[0112] 本发明方法的流程框图如图1所示，核心部分是：首先确定模拟电路的软故障集合和软故障特征参数；然后在基于每个软故障特征参数在各个已知软故障下的取值区间，
获取软故障诊断的模糊规则；获取软故障特征参数的在线监测值后，确定其选中的软故障
诊断的模糊规则，可将这些被选中模糊规则前项的归一化置信度推理到后项，得到对模糊
规则后项的归一化置信度；然后，对后项的置信度进行加权处理，得到对模拟电路软故障集
合中每个软故障的置信度，依照置信度最大准则判断是何故障发生。

[0113] 以下结合图5中ZPW-2000A无绝缘轨道电路系统软故障诊断的最佳实施例，详细介绍本发明方法的各个步骤。

[0114] 1、ZPW-2000A无绝缘轨道电路系统的功能与结构

[0115] ZPW-2000A无绝缘轨道电路的结构如图5所示。它是由主轨道电路(550m～1900m)和调谐区小轨道电路(29m)组成的电路网络，其工作过程如下：发送器发出正弦信
号经钢轨传送到接收器，当钢轨无车占用时，与发送器连接的继电器吸合，信号灯显示为绿
色。当列车通过时，传送信号被车轮短路，接收器电压降低，继电器落下，信号灯随即显示红
色表示“占用”状态。

[0116] 该系统的主要部件包括：发送器、SPT电缆、站防雷与电缆模拟网络、匹配变压器、电气绝缘节(由空心线圈、调谐单元及29m钢轨组成)、补偿电容、接收器和塞钉、引(导)
接线等连接元件等。它们的主要功能如下：

[0117] (1)发送器发出不同载频(1700+n×300Hz，n＝0～3)传送18种低频调制信号(10.3+n×1.1Hz，n＝0～17)，用于地车通信及运行控制。

[0118] (2)接收器用于接收主轨道电路信号，并在检查所属调谐区小轨道电路状态(XGJ、XGJH)条件下，动作本轨道电路的轨道继电器(GJ)；

[0119] (3)SPT电缆即铁路内屏蔽数字信号电缆用于信号的传输；

[0120] (4)电缆模拟网络与站防雷，前者是为了调整区间轨道电路传输特性，补偿实际SPT电缆，以便于轨道电路在列车不同运行方向时的电路调整，保证传输电路工作的稳定
性。后者是实现对传输电缆引入室内雷电冲击的防护，以保护模拟网络及室内发送、接收设
备；

[0121] (5)匹配变压器实现轨道电路和传输电缆的匹配连接；

[0122] (6)相邻轨道区段采用不同的载频频率，电气绝缘节根据谐振原理，限制某种载频的移频信号只能在本区段传送，而不能向相邻区段传送，防止传送信号的混淆，从而实现相
邻区段信号的电气绝缘。采用这种“无绝缘”方式，代替故障率较高的机械绝缘接头，在长
轨区段安装不用锯轨，这样可以提高轨道电路的可靠性，改善钢轨线路的运营质量。

[0123] (7)由于钢轨对信号呈现较高的感抗值，使轨道电路的传输衰耗较大。所以采取分段加补偿电容的方法减弱电感的影响，使轨道电路趋于阻性，增加了轨道电路的传输距离，
保证了轨道电路入口端的信号与干扰比。

[0124] 从电路分析的角度来看，除发送器和接收器为数字电子设备之外，其他部件都可等效成由电阻、电感和电容等模拟元件组成的模拟电路系统。由于钢轨振动和撞击、环境温
度、湿度变化等原因，这些部件易出现参数偏移的软故障。任意故障都会造成轨道电路功能
的丧失或部分失效，这将直接导致接收器输入电压超出容差范围，引起继电器的误动，造成
铁路占用信号的错误。并且，这些软故障都会引起发送端、接收端、匹配变压器两端、电气绝
缘节两端等多处节点电压、电流及某些部件温度等可检测量的变化，可以将这些量作为故
障特征信息加以分析处理。

[0125] 2、轨道电路Simulink仿真模型上的软故障设置及软故障特征参数的选取

[0126] 利用Matlab提供的Simulink仿真工具，建立ZPW-2000A无绝缘轨道电路的计算机仿真模型，该模型可以模拟载频1700Hz下，主轨长度为1200m的轨道电路，其中利用发送
端和接收端的电缆模拟网络模拟长度为10km的发送端和接收端SPT电缆。在此模型上模
拟表1中的10种软故障，采集该模型上的5个软故障特征参数，如表2所示

[0127] 表1 软故障设置列表

[0128]故障编号 软故障模式
1 正常状态(以接收端电压不超过容差范围为准)
2 发送端近端调谐单元电阻原件阻值增加2倍到5倍
3 发送端空心线圈电阻阻值增加20倍到50倍
4 发送端电缆模拟网络阻值增加20倍到50倍
5 发送端变压器电阻原件阻值增加100倍到200倍
6 接收端近端调谐单元电阻原件阻值增加2倍到5倍
7 接收端空心线圈电阻阻值增加20倍到50倍
8 接收端电缆模拟网络阻值增加20倍到50倍
9 接收端变压器电阻原件阻值增加100倍到200倍
10 发送端远端调谐单元电阻原件阻值增加2倍到5倍

[0129] 表2 软故障特征参数列表

[0130]标号 软故障特征参数信息采集位置与特征量
e1 发送端站防雷出口处交流电压(有效值)单位：伏特
e2 发送端站防雷出口处交流电流(有效值)单位：安培
e3 发送端电缆模拟网络出口处交流电压(有效值)单位：伏特
e4 接收端电缆模拟网络入口处交流电压(有效值)单位：伏特
e5 接收端站防雷出口处交流电压(有效值)单位：伏特

[0131] 这里将软故障1至软故障9设为已知软故障，将软故障10设定为未知软故障，以便后续对本发明方法进行测试。由步骤(1)-(2)可知，本实施例中所处理的故障集合为Θ
＝{F1，F2，F3，F4，F5，F6，F7，F8，F9，F10}，亦即j＝1，2，…9，JN＝9，其中F10＝FUN。它们共
同的软故障特征参数为ei i＝1，2，…，5，M＝5。根据步骤(3)-(4)可确定5个软故障特
征参数在9种已知软故障下的取值区间分别如表3、表4、表5、表6和表7所示，分别对应如
图6、图7、图8、图9和图10所示。

[0132] 表3软故障特征参数e1在9个已知软故障下的取值区间

[0133]

[0134] 表4软故障特征参数e2在9个已知软故障下的取值区间

[0135]

[0136] 表5软故障特征参数e3在9个已知软故障下的取值区间

[0137]

[0138] 表6软故障特征参数e4在9个已知软故障下的取值区间

[0139]

[0140] 表7软故障特征参数e5在9个已知软故障下的取值区间

[0141]

[0142] 3、得到三类软故障诊断的模糊规则

[0143] 利用本发明步骤(5)构造关于实施例中9个已知软故障、5个软故障特征参数的三类软故障诊断模糊规则。根据步骤(5-1)中的三角均分模糊区间划分方法，可以构造5个
故障特征参数的模糊语言集为

[0144] U1＝{A1，1，A1，2，A1，3，A1，4，A1，5，A1，6，A1，7，A1，8，A1，9，A1，10A1，11，A1，12}，J1＝12[0145] U2＝{A2，1，A2，2，A2，3，A2，4，A2，5，A2，6，A2，7，A2，8，A2，9，A2，10 A2，11，A2，12，A2，13}，J2＝13[0146] U3＝{A3，1，A3，2，A3，3，A3，4，A3，5，A3，6，A3，7，A3，8，A3，9，A3，10}，J3＝10

[0147] U4＝{A4，1，A4，2，A4，3，A4，4，A4，5，A4，6，A4，7}，J4＝7

[0148] U5＝{A5，1，A5，2，A5，3，A5，4，A5，5，A5，6，A5，7，A5，8}，J5＝8

[0149] 根据步骤(5-2)和(5-3-1)，得到第一类9个软故障的模糊规则的个数分别为NR1＝108、NR2＝240、NR3＝144、NR4＝2400、NR5＝540、NR6＝216、NR7＝72、NR8＝147、NR9
＝205。根据步骤(5-3-2)将前项相同但后项不同的模糊规则进行合并，得到第二类15个
软故障子集的模糊规则，这些模糊规则后项中的软故障子集(软故障用其标号表示)和这
些模糊规则的个数如表8所示。需要指出的是，各个软故障特征参数的语言项个数的选择
标准是：所生成的所有模糊规则中，至少各有一条模糊规则的后项是单个软故障，当模糊规
则的后项为软故障子集时，该软故障子集对应的模糊规则数应该最小。

[0150] 表8 所有软故障诊断的模糊规则的后项故障集类型及模糊规则个数的分配

[0151]}8，3
8 56 ，2{ 3

}9
7 04 ，8{ 92

}
8 9，8， }6，
6 81 3{ 5 1{ 4
}
5 045 }8，2{ 51 8，7，6{ 61
}9
004 ，3，2 }8，7 6
4 2 { 1 { 1
}9
50 }9，2 ，8，2
3 1 { 6 { 1
2 222 }3，2{ 3 }8，2{ 3
1 401 }9，3{ 11 }8，6{ 8
}9，8
NU 1243 26 ，3，2
F 8 9 1 { 1
数 数 数
项 则 项 则 项 则
后 规 后 规 后 规

[0152] 4、在轨道电路在线运行时，在某工况下，得到5种软故障特征参数的在线监测值，确定被监测值选中的模糊规则，并经模糊推理判断是何故障发生。

[0153] 以故障F8发生为例，轨道电路在线运行状态下，对5种软故障特征参数进行观测得到监测值，根据步骤(6)可以计算该组监测值所选中的模糊语言项及其归一化后的置信
度取值，如表9所示。

[0154] 表9 5种软故障特征参数的在线监测值及被选中语言项的归一化置信度

[0155]

[0156] 根据步骤(6)可以得到在该组监测值所选中的软故障诊断的模糊规则的前项如表10所示(规则前项用相应模糊语言项的组合形式表示)：

[0157] 表10 被选中的软故障诊断的模糊规则的前项及置信度

[0158]

[0159]序号 被选中的前项 置信度 序号 被选中的前项 置信度
3 8 3 1 2 4 0 1 0 8 3 4 3 3 9
841. 600. 860. 300. 250. 200. 420. 100. 810. 000. 800. 000. 600. 000. 200.
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
A6.1 A5.1 A6.1 A5.1 A6.1 A5.1 A6.1 A5.1 A6.1 A5.1 A6.1 A5.1 A6.1 A5.1 A6.1
×A5，1 ×A5，1 ×A6，1 ×A6，1 ×A5，1 ×A5，1 ×A6，1 ×A6，1 ×A5，1 ×A5，1 ×A6，1 ×A6，1 ×A5，1 ×A5，1 ×A6，1
×A6，1 ×A6，1 ×A6，1 ×A6，1 ×A7，1 ×A7，1 ×A7，1 ×A7，1 ×A6，1 ×A6，1 ×A6，1 ×A6，1 ×A7，1 ×A7，1 ×A7，1
×A3，1 ×A3，1 ×A3，1 ×A3，1 ×A3，1 ×A3，1 ×A3，1 ×A3，1 ×A4，1 ×A4，1 ×A4，1 ×A4，1 ×A4，1 ×A4，1 ×A4，1
×A 6，1 ×A 6，1 ×A 6，1 ×A 6，1 ×A 6，1 ×A 6，1 ×A 6，1 ×A 6，1 ×A 6，1 ×A 6，1 ×A 6，1 ×A 6，1 ×A 6，1 ×A 6，1 ×A 6，1

71 81 91 02 12 22 32 42 52 62 72 82 92 03 13

33 13 60 06 79 64 06 12 44 61 95 70 12 60 65
82.0 10.0 31.0 00.0 90.0 00.0 40.0 00.0 30.0 00.0 10.0 00.0 10.0 00.0 00.0

A6，1 A5，1 A6，1 A5，1 A6，1 A5，1 A6，1 A5，1 A6，1 A5，1 A6，1 A5，1 A6，1 A5.1 A6.1
×A5，1 ×A5，1 ×A6，1 ×A6，1 ×A5，1 ×A5，1 ×A6，1 ×A6，1 ×A5，1 ×A5，1 ×A6，1 ×A6，1 ×A5，1 ×A5，1 ×A6，1
×A6，1 ×A6，1 ×A6，1 ×A6，1 ×A7，1 ×A7，1 ×A7，1 ×A7，1 ×A6，1 ×A6，1 ×A6，1 ×A6，1 ×A7，1 ×A7，1 ×A7，1
×A3，1 ×A3，1 ×A3，1 ×A3，1 ×A3，1 ×A3，1 ×A3，1 ×A3，1 ×A4，1 ×A4，1 ×A4，1 ×A4，1 ×A4，1 ×A4，1 ×A4，1
×A 5，1 ×A 5，1 ×A 5，1 ×A 5，1 ×A 5，1 ×A 5，1 ×A 5，1 ×A 5，1 ×A 5，1 ×A 5，1 ×A 5，1 ×A 5，1 ×A 5，1 ×A 5，1 ×A 5，1
0 1 2 3 4 5
1 2 3 4 5 6 7 8 9 1 1 1 1 1 1
1
000.
0
A5.1
×A6，1
×A7，1
×A4，1
×A 6，1

23

30
00.0

A5.1
×A6，1
×A7，1
×A4，1
×A 5，1
6
1

[0160] 表10中实际上给出的是被选中的JNR＝32条模糊规则的前项，根据步骤(7)可得它们各自对应的后项及其置信度如表11所示：

[0161] 表11 被选中模糊规则的后项及置信度

[0162]序号 选中规则的后项 置信度 序号 被选中的后项 置信度
1 {F7，F8} 0.2833 17 {F7，F8} 0.1483
2 {F8，F9} 0.0131 18 {F8，F9} 0.0068
3 {F8} 0.1306 19 {F8} 0.0683
4 {F8，F9} 0.0060 20 {F8，F9} 0.0031
5 {F7，F8} 0.0997 21 {F7，F8} 0.0522
6 {F8，F9} 0.0046 22 {F8，F9} 0.0024
7 {F8} 0.0460 23 {F8} 0.0240
8 {F8，F9} 0.0021 24 {F8，F9} 0.0011
9 {F7，F8} 0.0344 25 {F7，F8} 0.0180
10 {F8，F9} 0.0016 26 {F8，F9} 0.0008
11 {F8} 0.0159 27 {F8} 0.0083
12 {F8，F9} 0.0007 28 {F8，F9} 0.0004
13 {F7，F8} 0.0121 29 {F7，F8} 0.0063
14 {F8，F9} 0.0006 30 {F8，F9} 0.0003
15 {F8} 0.0056 31 {F8} 0.0029
16 {F8，F9} 0.0003 32 {F8，F9} 0.0001

[0163]

[0164] 根据步骤(8)，对被选中模糊规则后项的置信度进行加权处理获得单个软故障的置信度，如表12所示：

[0165] 表12 模糊推理结果

[0166]NU
F 0
0220
F 9 .0
805
F 8 6.0

2
7 723.
F 0
F 5 0
F4 0

F 3 0
F2 0

F1 0

′
ρ

[0167] 根据步骤(9)中的故障决策准则可判断是故障F8发生，与实际设置的已知软故障模式一致。

[0168] 对于事先设定的未软知故障F10＝FUN，，其在5个软故障特征参数下的取值区间如表13所示：

[0169] 表13 未知故障F10在5种故障特征参数下的取值

[0170]]538
1.73
1，07
SE 01，5 15.631[

]
7484
.831
，068
01，4 0.73
SE 1[
]38
44.4
31，7
SE01，3 835.331[

]1
961.
531，
3274
S01，2 .431
E [
]73
34.6
31，1
SE01，1 988.531[

[0171] 在以上表13的5个取值区间中分别选取200个监测样本，对它们分别加0.5％～7％不等的扰动，用这些加扰动后的监测样本测试本发明所给出的模糊推理系统，检测出未
知故障F10发生的机率为98％。若实际中出现该种情况，即可判断是未知软故障发生，可以
及时排查轨道电路确定是何故障发生。

[0172] 对于9种已知软故障，在对应各软故障的5个软故障特征参数的取值区间中，随机选取200组监测样本，对它们分别加0.5％～7％不等的扰动，用加扰动后的监测样本进行
测试，软故障平均确诊率为93％，达到了一般模拟电路故障诊断系统的确诊率要求。

[0173] 需要说明的是，选择0.5％～7％的扰动量是根据实际软故障特征采集装置的观测误差确定的，加过扰动后的监测样本会落在各自故障特征取值区间之外，但是大多数此
类情况下，本发明所提出的软故障诊断方法还能正确诊断故障。