发明名称 一种基于MCSKPCA的神经网络模拟电路故障诊断方法
摘要 本发明公开了一种基于MCSKPCA的神经网络模拟电路故障诊断方法,包括以下步骤:步骤1:采集被诊断模拟电路的输出电压信号;步骤2:对采集的输出电压信号进行小波变换;步骤3:计算输出电压信号经小波变换得到的小波系数的能量特征值;步骤4:将所述的能量特征值进行MCSKPCA的特征提取和降维处理得到最优特征向量;步骤5:将所述的最优特征向量送入BP神经网络分离器中,由BP神经网络分离器输出故障诊断结果。该方法不但可以用于线性和非线性电路及其系统的诊断,同时也可以用于诊断所述电路中的硬故障和软故障。
申请公布号 CN102279358B 申请公布日期 2013.11.27
申请号 CN201110166548.X 申请日期 2011.06.20
申请人 湖南大学 发明人 何怡刚;肖迎群;方葛丰;阳辉
分类号 G01R31/316(2006.01)I;G06N3/08(2006.01)I 主分类号 G01R31/316(2006.01)I
代理机构 长沙市融智专利事务所 43114 代理人 黄美成
主权项 1.一种基于最大类别分离度核主元分析(KPCA based on Maximal Class Separability即MCSKPCA)的神经网络模拟电路故障诊断方法,其特征在于包括以下步骤:1)采集模拟电路的电信号,采集到的电信号为所述的模拟电路的输出电压信号;2)对采集的所述电信号进行小波变换;3)计算所述电信号的小波逼近系数能量特征值;4)将所述的小波逼近系数能量特征值所构成的候选特征向量数据集进行MCSKPCA的特征提取和维数降低得到最优特征向量;5)将所述的最优特征向量送入BP神经网络分离器中,由BP神经网络分离器输出故障诊断结果;步骤4中,首先进行MCSKPCA核参数优化准则的设计和参数的优化后,再利用MCSKPCA实施相应的特征提取;MCSKPCA的核参数优化准则确定为基于核的类别分离度准则为:<maths num="0001"><![CDATA[<math><mrow><mi>J</mi><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>tr</mi><mrow><mo>(</mo><msubsup><mi>S</mi><mi>B</mi><mi>&phi;</mi></msubsup><mo>)</mo></mrow></mrow><mrow><mi>tr</mi><mrow><mo>(</mo><msubsup><mi>S</mi><mi>W</mi><mi>&phi;</mi></msubsup><mo>)</mo></mrow></mrow></mfrac><mo>;</mo></mrow></math>]]></maths>而最优核参数集θ<sub>opt</sub>可由下式得到:<maths num="0002"><![CDATA[<math><mrow><msub><mi>&theta;</mi><mi>opt</mi></msub><mo>=</mo><mi>arg</mi><munder><mi>max</mi><mrow><mi>&theta;</mi><mo>&Element;</mo><mi>&Theta;</mi></mrow></munder><mo>[</mo><mfrac><mrow><mi>tr</mi><mrow><mo>(</mo><msubsup><mi>S</mi><mi>B</mi><mi>&phi;</mi></msubsup><mo>)</mo></mrow></mrow><mrow><mi>tr</mi><mrow><mo>(</mo><msubsup><mi>S</mi><mi>W</mi><mi>&phi;</mi></msubsup><mo>)</mo></mrow></mrow></mfrac><mo>]</mo><mo>;</mo></mrow></math>]]></maths>式中Θ表示核参数空间;选定高斯RBF核作为所需要的核函数即K(x<sub>i</sub>,x<sub>j</sub>)=exp(-||x<sub>i</sub>-x<sub>j</sub>||<sup>2</sup>/2σ<sup>2</sup>),经过对J(θ)的推导,得到最终的基于核的度量准则方程为:<img file="FDA00003223886900013.GIF" wi="299" he="82" />式中,tr(·)表示矩阵的迹,<img file="FDA00003223886900014.GIF" wi="142" he="84" />和<img file="FDA00003223886900015.GIF" wi="138" he="86" />分别表示核类间散布矩阵和核类内散布矩阵的迹,<img file="FDA00003223886900016.GIF" wi="60" he="78" />表示核类间散布矩阵,<img file="FDA00003223886900017.GIF" wi="73" he="79" />表示核类内散布矩阵;由于所采用的高斯RBF核只具有一个核参数即高斯宽度σ,那么,J<sub>l</sub>(θ)和θ<sub>opt</sub>就相应地变为J<sub>l</sub>(σ)和最优高斯宽度σ<sub>opt</sub>;因此,<maths num="0003"><![CDATA[<math><mrow><msub><mi>J</mi><mi>l</mi></msub><mrow><mo>(</mo><mi>&sigma;</mi><mo>)</mo></mrow><mo>=</mo><mi>tr</mi><mrow><mo>(</mo><msubsup><mi>S</mi><mi>B</mi><mi>&phi;</mi></msubsup><mo>)</mo></mrow><mo>,</mo><msub><mi>&sigma;</mi><mi>opt</mi></msub><mo>=</mo><mi>arg</mi><mi></mi><mi>max</mi><mo>[</mo><mi>tr</mi><mrow><mo>(</mo><msubsup><mi>S</mi><mi>B</mi><mi>&phi;</mi></msubsup><mo>)</mo></mrow><mo>]</mo><mo>;</mo></mrow></math>]]></maths>最后,得到MCSKPCA算法如下:给定训练样本<img file="FDA00003223886900021.GIF" wi="485" he="80" />R<sup>d</sup>为输入空间,l指样本的个数,非线性映射<img file="FDA000032238869000210.GIF" wi="483" he="70" />和核函数K(·,·),在原输入空间对训练数据进行归一化处理;根据式<img file="FDA00003223886900022.GIF" wi="464" he="94" />计算最优高斯宽度σ<sub>opt</sub>;计算核矩阵<maths num="0004"><![CDATA[<math><mrow><mover><mi>K</mi><mo>~</mo></mover><mo>=</mo><mrow><mo>(</mo><mover><mi>K</mi><mo>~</mo></mover><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><msub><mi>x</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><msup><mrow><mo>|</mo><mo>|</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mi>x</mi><mi>j</mi></msub><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><mo>/</mo><mn>2</mn><msubsup><mi>&sigma;</mi><mi>opt</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mo>)</mo></mrow><mo>,</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><mi>l</mi><mo>;</mo></mrow></math>]]></maths>执行<maths num="0005"><![CDATA[<math><mrow><mi>K</mi><mo>=</mo><mover><mi>K</mi><mo>~</mo></mover><msub><mrow><mo>-</mo><mn>1</mn></mrow><mi>l</mi></msub><mover><mi>K</mi><mo>~</mo></mover><mo>-</mo><msub><mrow><mover><mi>K</mi><mo>~</mo></mover><mn>1</mn></mrow><mi>l</mi></msub><mo>+</mo><msub><mn>1</mn><mi>l</mi></msub><msub><mrow><mover><mi>K</mi><mo>~</mo></mover><mn>1</mn></mrow><mi>l</mi></msub><mo>;</mo></mrow></math>]]></maths>根据<img file="FDA00003223886900025.GIF" wi="312" he="159" />j=1,…,p计算正交归一化特征向量v<sub>j</sub>,这里j=1,…,p,p是矩阵K分解后的特征向量的个数,不同的矩阵其特征向量的个数不同;<img file="FDA00003223886900026.GIF" wi="234" he="82" />Q=[Φ(x<sub>1</sub>),…,Φ(x<sub>l</sub>)],λ<sub>j</sub>和γ<sub>j</sub>分别为K的特征值和特征向量;给定一个测试数据x,由式(7)就可以得到对应的非线性主元集,即<maths num="0006"><![CDATA[<math><mrow><msub><mi>y</mi><mi>j</mi></msub><mo>=</mo><msubsup><mi>v</mi><mi>j</mi><mi>T</mi></msubsup><mi>&Phi;</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><msqrt><msub><mi>&lambda;</mi><mi>j</mi></msub></msqrt></mfrac><msubsup><mi>&gamma;</mi><mi>j</mi><mi>T</mi></msubsup><msup><mi>Q</mi><mi>T</mi></msup><mi>&Phi;</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow></mrow></math>]]></maths>               (7)<maths num="0007"><![CDATA[<math><mrow><mo>=</mo><mfrac><mn>1</mn><msqrt><msub><mi>&lambda;</mi><mi>j</mi></msub></msqrt></mfrac><msubsup><mi>&gamma;</mi><mi>j</mi><mi>T</mi></msubsup><mo>[</mo><mi>k</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>,</mo><mi>x</mi><mo>)</mo></mrow><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><mi>k</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>l</mi></msub><mo>,</mo><mi>x</mi><mo>)</mo></mrow><mo>]</mo><mo>,</mo><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><mi>p</mi></mrow></math>]]></maths>式中,y<sub>j</sub>为第j个抽取的主元,Φ(x)为测试数据x在特征空间<img file="FDA000032238869000211.GIF" wi="56" he="57" />中的像,k(x<sub>1</sub>,x)就是矩阵K的一个元素;p个主元就构成最终的最优特征向量。
地址 410082 湖南省长沙市岳麓区麓山南路2号