发明名称 基于核主元分析贡献图的非线性过程故障辨识方法
摘要 一种基于核主元分析贡献图的非线性故障辨识方法,包括采集数据、提取特征空间的基、提取非线性主元、故障检测和故障辨识五个步骤,本发明提出了一种新的特征空间基的提取方法,通过基的提取来减少了数据映射到特征空间后的线性冗余,以及在训练样本数量很大的时候减少KPCA的计算量。本发明方法采用贡献图进行故障辨识,对过程的每个变量计算其贡献图和控制限,通过贡献图和控制限的关系判断出各变量失控对故障的发生所承担的责任,此辨识方法克服了输入空间和特征空间不能自由转换所带来的故障辨识的困难。
申请公布号 CN101169623A 申请公布日期 2008.04.30
申请号 CN200710158468.3 申请日期 2007.11.22
申请人 东北大学 发明人 张颖伟;秦泗钊;王滢
分类号 G05B13/04(2006.01) 主分类号 G05B13/04(2006.01)
代理机构 沈阳东大专利代理有限公司 代理人 梁焱
主权项 1.一种基于核主元分析贡献图的非线性故障辨识方法,其特征在于包括以下步骤:步骤一、采集数据采集过程相关变量的数据,并对其进行中心化和归一化,即规范化处理;步骤二、提取特征空间的基通过选取一个特征空间的基来提取特征,根据特征因子依次判断规范化后的数据与基中数据的线性关系,并去除线性关系较强的数据,其余数据加入基中,最后得到基Ω;步骤三、提取非线性主元通过非线性映射把输入空间的观测数据x映射到特征空间F,接着在特征空间F执行线性主元分析,非线性主元计算公式如下:<maths num="0001"><![CDATA[<math><mrow><msub><mi>t</mi><mi>k</mi></msub><mo>=</mo><mo>&lang;</mo><msub><mi>v</mi><mi>k</mi></msub><mo>,</mo><mover><mi>&Phi;</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>&rang;</mo><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mi>&alpha;</mi><mi>i</mi><mi>k</mi></msubsup><mo>&lang;</mo><mover><mi>&Phi;</mi><mo>~</mo></mover><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>,</mo><mover><mi>&Phi;</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>&rang;</mo><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mi>&alpha;</mi><mi>i</mi><mi>k</mi></msubsup><mover><mi>k</mi><mo>~</mo></mover><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><mi>x</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>24</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中t<sub>k</sub>-非线性主元t的第k个元素;k=1,2,...,p,其中p是非线性主元要提取的个数;v<sub>k</sub>-训练数据协方差矩阵的第k特征向量;<img file="S2007101584683C00012.GIF" wi="98" he="56" />-Φ为非线性映射函数,<img file="S2007101584683C00013.GIF" wi="98" he="57" />为中心化处理后的数据;α<sub>i</sub><sup>k</sup>-中心化核矩阵的第k个特征向量的第i个元素;x-要进行非线性主元提取的数据;x<sub>i</sub>-训练数据中的第i个数据,i=1,2,…,N;<img file="S2007101584683C00014.GIF" wi="143" he="62" />-中心化核向量<img file="S2007101584683C00015.GIF" wi="28" he="49" />的第i个值,其中<maths num="0002"><![CDATA[<math><mrow><mover><mi>k</mi><mo>~</mo></mover><mo>=</mo><mi>k</mi><mo>-</mo><msub><mn>1</mn><mi>t</mi></msub><mi>K</mi><mo>-</mo><msub><mrow><mi>k</mi><mn>1</mn></mrow><mi>N</mi></msub><mo>+</mo><msub><mn>1</mn><mi>t</mi></msub><msub><mrow><mi>K</mi><mn>1</mn></mrow><mi>N</mi></msub><mo>,</mo></mrow></math>]]></maths>k=[k(x,x<sub>1</sub>),...,k(x,x<sub>N</sub>)];步骤四、故障检测利用Hotelling的T<sup>2</sup>统计和平方预测误差SPE统计进行故障检测,通过控制图即可判断出是否发生故障,当T<sup>2</sup>统计和SPE统计超出各自的控制限时,认为有故障发生,反之,整个过程正常;步骤五、故障辨识当步骤四检测到故障发生,则采用贡献图进行故障辨识,对过程的每个变量计算其贡献图和控制限,通过贡献图和控制限的关系即可判断出各变量失控对故障的发生所承担的责任;当变量的贡献图小于其控制限时,则此变量对故障的发生要承担一定的责任;根据差值的大小判断责任的大小。
地址 110004辽宁省沈阳市和平区文化路3号巷11号