发明名称 一种基于结构层次关系的改进相关性矩阵分析方法
摘要 本发明提出一种基于结构层次关系的改进相关性矩阵分析方法,包括构建系统的结构层次信息、梳理故障、梳理测试、建立初始相关性矩阵、判断各结构层次中扩展故障与扩展测试的关联性、建立改进相关性矩阵等几个步骤,还可以选择性的包括计算特定结构层次对象的故障检测率的步骤和计算特定结构层次对象的故障隔离率的步骤。本发明将现有二值相关性矩阵改进为三值相关性矩阵,可以剔除很多在结构上与故障不相关的测试,在故障诊断阶段减少需要检测的测试数,节省了大量的计算资源并有助于故障的快速诊断,通过改进相关性矩阵,可以实现计算任意系统层次单元的故障检测率和故障隔离率。
申请公布号 CN102243674A 申请公布日期 2011.11.16
申请号 CN201110180202.5 申请日期 2011.06.29
申请人 北京航空航天大学 发明人 石君友;王璐;侯文魁;王风武
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 北京永创新实专利事务所 11121 代理人 官汉增
主权项 1.一种基于结构层次关系的改进相关性矩阵分析方法,其特征在于:具体包括以下几个步骤:步骤一、构建系统的结构层次信息:将待诊断的工程系统按照自顶向下、连续性、长度不限和信息唯一性原则,确定系统的实际结构层次划分,构建系统结构层次信息UM;步骤二、梳理故障:扩展故障定义为包括结构层次信息的故障扩展描述形式,扩展故障的元组模型为:F=(UM FM)式中,F为扩展故障,UM是该扩展故障所属的结构层次信息,FM是该扩展故障的故障模式信息;在扩展故障定义的基础上,梳理出系统各结构层次上的故障并将其转换为扩展故障表达形式:F=UM//FM_name,name是故障模式名称;步骤三、梳理测试:扩展测试定义为包括结构信息的测试扩展描述形式,扩展测试的元组模型如下:T=(UM TM)式中,T为扩展测试,UM是测试所属结构层次信息,TM是测试功能信息;在扩展测试定义的基础上,梳理出系统各结构层次上的功能测试并将其转换为扩展测试表达形式:T=UM//TM_name′,name′是测试功能名称;步骤四、建立初始相关性矩阵:根据扩展故障和扩展测试,采用建模方法,并且要求同一结构层次下的所有扩展故障和同一层次下的所有扩展测试依次顺序书写,得到初始相关性矩阵D<sub>m×n</sub>,表达式如下:<maths num="0001"><![CDATA[<math><mrow><msub><mi>D</mi><mrow><mi>m</mi><mo>&times;</mo><mi>n</mi></mrow></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>MTTF</mi><mn>1</mn></msub></mtd><mtd><msub><mi>d</mi><mn>11</mn></msub></mtd><mtd><msub><mi>d</mi><mn>12</mn></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>d</mi><mrow><mn>1</mn><mi>n</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>MTTF</mi><mn>2</mn></msub></mtd><mtd><msub><mi>d</mi><mn>21</mn></msub></mtd><mtd><msub><mi>d</mi><mn>22</mn></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>d</mi><mrow><mn>2</mn><mi>n</mi></mrow></msub></mtd></mtr><mtr><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>MTTF</mi><mi>m</mi></msub></mtd><mtd><msub><mi>d</mi><mrow><mi>m</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>d</mi><mrow><mi>m</mi><mn>2</mn></mrow></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>d</mi><mi>mn</mi></msub></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>其中,第i行行向量为:F<sub>i</sub>=[MTTF<sub>i</sub> d<sub>i1</sub> d<sub>i2</sub>...d<sub>ij</sub>...d<sub>in</sub>]MTTF<sub>i</sub>表示第i个扩展故障F<sub>i</sub>的平均故障间隔时间;第i个扩展故障F<sub>i</sub>在各扩展测试上的反应信息是d<sub>i1</sub> d<sub>i2</sub>...d<sub>in</sub>,它表明了扩展故障F<sub>i</sub>与各个扩展测试T<sub>j</sub>(j=1,2,…,n)的相关性,d<sub>i</sub>=(d<sub>i1</sub> d<sub>i2</sub>...d<sub>ij</sub>...d<sub>in</sub>)表示扩展故障F<sub>i</sub>的相关性值,其中d<sub>ij</sub>表示第i行扩展故障和第j个扩展测试的相关性值;<img file="FDA0000072322940000021.GIF" wi="1119" he="145" />步骤五、判断各结构层次中扩展故障与扩展测试的关联性:结构层次为UM<sub>x</sub>的扩展故障与结构层次为UM<sub>y</sub>的扩展测试的关联性定义为结构层次UM<sub>x</sub>上的扩展故障和结构层次UM<sub>y</sub>上的扩展测试是否具有通信关系的描述,该通信关系的表示形式为<img file="FDA0000072322940000022.GIF" wi="253" he="60" />若<img file="FDA0000072322940000023.GIF" wi="307" he="60" />表示结构层次UM<sub>x</sub>上的扩展故障与结构层次UM<sub>y</sub>上的扩展测试相关联;若<img file="FDA0000072322940000024.GIF" wi="307" he="60" />表示结构层次为UM<sub>x</sub>的扩展故障与结构层次为UM<sub>y</sub>的扩展测试不关联;步骤5.1系统所有结构层次信息构成的结构层次信息集合为UM<sub>K</sub>,UM<sub>K</sub>={UM<sub>1</sub>,UM<sub>2</sub>,...,UM<sub>k</sub>},其中k为结构层次信息的总个数;步骤5.2建立空集合D′用于存储结构层次为UM<sub>x</sub>的所有扩展故障与结构层次为UM<sub>y</sub>的所有扩展测试的关联性;步骤5.3在结构层次信息集合{UM<sub>K</sub>}中选择第一个结构层次信息UM<sub>x</sub>,UM<sub>x</sub>∈UM<sub>K</sub>;步骤5.4确定结构层次信息为UM<sub>x</sub>的所有扩展故障,存入集合F<sub>A</sub>中;步骤5.4.1提取初始相关性矩阵的第一行扩展故障作为当前扩展故障;步骤5.4.2当前扩展故障的结构层次信息UM<sub>f</sub>;如果结构层次信息UM<sub>f</sub>等于结构层次信息UM<sub>x</sub>,则把当前扩展故障存入集合F<sub>A</sub>;否则,执行步骤5.4.3;步骤5.4.3若初始相关性矩阵中还有未遍历的扩展故障行,则选择下一行扩展故障作为当前扩展故障,返回步骤5.4.2;步骤5.5将集合F<sub>A</sub>中各扩展故障行相关性值d<sub>i</sub>按列作或运算,建立1×n的子相关性矩阵F<sub>A</sub>′,F<sub>A</sub>′=[(d<sub>11</sub>||d<sub>21</sub>||...||d<sub>i1</sub>)(d<sub>12</sub>||d<sub>22</sub>||...||d<sub>i2</sub>)...(d<sub>1j</sub>||d<sub>2j</sub>||...||d<sub>ij</sub>)...(d<sub>1n</sub>||d<sub>2n</sub>||...||d<sub>in</sub>)];步骤5.6在结构层次信息集合{UM<sub>K</sub>/UM<sub>x</sub>}中选择第一个结构层次信息UM<sub>y</sub>,UM<sub>y</sub>∈{UM<sub>K</sub>/UM<sub>x</sub>};步骤5.7将子相关性矩阵F<sub>A</sub>′中结构层次信息为UM<sub>y</sub>的所有扩展测试存入集合T<sub>A</sub>中;步骤5.7.1提取子相关性矩阵F<sub>A</sub>′的第一列扩展测试,作为当前扩展测试;步骤5.7.2当前扩展测试的结构层次信息为UM<sub>t</sub>;若结构层次信息UM<sub>t</sub>等于结构层次信息UM<sub>y</sub>,则把当前扩展测试存入集合T<sub>A</sub>;否则,执行步骤5.7.3;步骤5.7.3若子相关性矩阵F<sub>A</sub>′中还有未遍历的扩展测试列,则选择下一列扩展测试作为当前扩展测试,并返回步骤5.7.2;步骤5.8将集合T<sub>A</sub>中各测试列值d<sub>1j</sub>作或运算,记为<img file="FDA0000072322940000031.GIF" wi="666" he="60" />并将其存入集合D′中;步骤5.9在结构层次信息集合{UM<sub>K</sub>/UM<sub>x</sub>}中选取下一个结构层次信息UM<sub>y+1</sub>,UM<sub>y+1</sub>∈{UM<sub>K</sub>/UM<sub>i</sub>},令UM<sub>y</sub>=UM<sub>y+1</sub>,返回步骤5.7,直至遍历完集合{UM<sub>K</sub>/UM<sub>i</sub>}中的所有元素;否则,执行步骤5.10;步骤5.10在结构层次信息集合UM<sub>K</sub>中选取下一个结构层次信息U<sub>M+1</sub>,U<sub>M+1</sub>∈UM<sub>K</sub>,令UM<sub>x</sub>=UM<sub>x+1</sub>,返回步骤5.4,直至遍历完结构层次信息集合UM<sub>K</sub>中的所有元素,得到存有结构层次UM<sub>x</sub>上的扩展故障和结构层次UM<sub>y</sub>上的扩展测试的通信关系<img file="FDA0000072322940000032.GIF" wi="230" he="60" />的集合D′;步骤六、建立改进相关性矩阵:根据各结构层次的扩展故障与扩展测试的关联性,改进初始相关性矩阵,建立改进相关性矩阵,具体步骤如下:步骤6.1从初始相关性矩阵D<sub>m×n</sub>中,选择第一行扩展故障作为当前扩展故障F<sub>l</sub>;步骤6.2当前扩展故障F<sub>l</sub>的平均故障间隔时间MTTF保持不变,存入改进相关性矩阵;步骤6.3提取当前扩展故障F<sub>l</sub>的结构层次信息UM<sub>l</sub>;步骤6.4提取当前扩展故障F<sub>l</sub>的第一列扩展测试作为当前扩展测试T<sub>s</sub>;步骤6.5提取当前扩展测试T<sub>s</sub>的结构层次信息UM<sub>s</sub>;步骤6.6当前扩展故障F<sub>l</sub>与当前扩展测试T<sub>s</sub>的相关性值d<sub>ij</sub>改进为d′<sub>ij</sub>:步骤6.6.1如果UM<sub>l</sub>=UM<sub>s</sub>,则d<sub>ij</sub>′=d<sub>ij</sub>,将d′<sub>ij</sub>存入改进相关性矩阵;否则执行步骤6.6.2;步骤6.6.2遍历集合D′,找到通信关系<img file="FDA0000072322940000033.GIF" wi="223" he="56" />的值;步骤6.6.3如果<img file="FDA0000072322940000034.GIF" wi="302" he="56" />则d<sub>ij</sub>′=d<sub>ij</sub>,将d′<sub>ij</sub>存入改进相关性矩阵,否则,d<sub>ij</sub>′=×,将d′<sub>ij</sub>存入改进相关性矩阵;步骤6.7若当前扩展故障F<sub>l</sub>中还有未遍历的扩展测试列,则选择下一列扩展测试作为当前扩展测试T<sub>s</sub>,返回步骤6.5;步骤6.8若初始相关性矩阵中,还有未遍历的扩展故障行,则选择下一行扩展故障作为当前扩展故障F<sub>l</sub>,返回步骤6.2,直至遍历完初始相关性矩阵的全部行和列,得到{0,1,×}三值的改进相关性矩阵D′<sub>m×n</sub>:<maths num="0002"><![CDATA[<math><mrow><msubsup><mi>D</mi><mrow><mi>m</mi><mo>&times;</mo><mi>n</mi></mrow><mo>&prime;</mo></msubsup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>MTTF</mi><mn>1</mn></msub></mtd><mtd><msubsup><mi>d</mi><mn>11</mn><mo>&prime;</mo></msubsup></mtd><mtd><msubsup><mi>d</mi><mn>12</mn><mo>&prime;</mo></msubsup></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msubsup><mi>d</mi><mrow><mn>1</mn><mi>n</mi></mrow><mo>&prime;</mo></msubsup></mtd></mtr><mtr><mtd><msub><mi>MTTF</mi><mn>2</mn></msub></mtd><mtd><msubsup><mi>d</mi><mn>21</mn><mo>&prime;</mo></msubsup></mtd><mtd><msubsup><mi>d</mi><mn>22</mn><mo>&prime;</mo></msubsup></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msubsup><mi>d</mi><mrow><mn>2</mn><mi>n</mi></mrow><mo>&prime;</mo></msubsup></mtd></mtr><mtr><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>MTTF</mi><mi>m</mi></msub></mtd><mtd><msubsup><mi>d</mi><mrow><mi>m</mi><mn>1</mn></mrow><mo>&prime;</mo></msubsup></mtd><mtd><msubsup><mi>d</mi><mrow><mi>m</mi><mn>2</mn></mrow><mo>&prime;</mo></msubsup></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msubsup><mi>d</mi><mi>mn</mi><mo>&prime;</mo></msubsup></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>其中:<img file="FDA0000072322940000042.GIF" wi="1427" he="224" />
地址 100191 北京市海淀区学院路37号