发明名称 用于铝电解槽况故障诊断的优化权重相对主元分析方法
摘要 一种用于铝电解槽况故障诊断的优化权重相对主元分析方法,其特征在于:一,随机采集n组铝电解槽况数据组成原始测量样本集;二,得到标准化后样本矩阵X<sup>*</sup>;三,随机产生相对转换矩阵Λ;四,得到相对化样本矩阵X<sup>R</sup>;五,对X<sup>R</sup>进行主元分析,计算SPE检验的误报率;六,利用遗传算法优化相对转换矩阵Λ,得到最优的相对转换矩阵Λ<sup>*</sup>;七,利用最优的相对转换矩阵Λ<sup>*</sup>,得到最优的相对化样本矩阵X<sup>RZ</sup>;八,实现铝电解槽况故障诊断。本发明利用遗传算法优化主元权重,将“均匀”分布的量突显出来,以便更好地提取出具有代表性的主元,从而提高铝电解槽况故障诊断的精确度。
申请公布号 CN103952724B 申请公布日期 2016.04.20
申请号 CN201410165336.3 申请日期 2014.04.22
申请人 重庆科技学院 发明人 易军;黄迪;李太福;周伟;张元涛;姚立忠;田应甫
分类号 C25C3/20(2006.01)I;G06F17/16(2006.01)I 主分类号 C25C3/20(2006.01)I
代理机构 重庆为信知识产权代理事务所(普通合伙) 50216 代理人 余锦曦
主权项 一种用于铝电解槽况故障诊断的优化权重相对主元分析方法,其特征在于:包括下列步骤:步骤一,随机采集n组铝电解槽况数据组成原始测量样本集X=[x<sub>1</sub>,x<sub>2</sub>,…,x<sub>n</sub>]∈R<sup>m×n</sup>,每个样本含有m个独立的铝电解槽况参数采样值;步骤二,对原始测量样本集X进行标准化处理,得到标准化后样本矩阵X<sup>*</sup>;步骤三,在[0,50)范围内随机产生相对转换矩阵Λ,所述相对转换矩阵Λ为对角矩阵:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>&Lambda;</mi><mo>=</mo><mfenced open = "(" close = ")"><mtable><mtr><mtd><msub><mi>&lambda;</mi><mn>1</mn></msub></mtd><mtd><mn>0</mn></mtd><mtd><mo>...</mo></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msub><mi>&lambda;</mi><mn>2</mn></msub></mtd><mtd><mo>...</mo></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mo>...</mo></mtd><mtd><msub><mi>&lambda;</mi><mi>m</mi></msub></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000869066830000011.GIF" wi="486" he="303" /></maths>步骤四,对标准化后样本矩阵X<sup>*</sup>进行相对化转换,得到相对化样本矩阵X<sup>R</sup>,所述相对化转换按下式进行:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msup><mi>X</mi><mi>R</mi></msup><mo>=</mo><mi>&Lambda;</mi><mo>&times;</mo><msup><mi>X</mi><mo>*</mo></msup><mo>=</mo><mfenced open = "(" close = ")"><mtable><mtr><mtd><msub><mi>&lambda;</mi><mn>1</mn></msub></mtd><mtd><mn>0</mn></mtd><mtd><mo>...</mo></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msub><mi>&lambda;</mi><mn>2</mn></msub></mtd><mtd><mo>...</mo></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mo>...</mo></mtd><mtd><msub><mi>&lambda;</mi><mi>m</mi></msub></mtd></mtr></mtable></mfenced><mo>&times;</mo><mfenced open = "(" close = ")"><mtable><mtr><mtd><msubsup><mi>x</mi><mn>11</mn><mo>*</mo></msubsup></mtd><mtd><msubsup><mi>x</mi><mrow><mi>2</mi><mn>1</mn></mrow><mo>*</mo></msubsup></mtd><mtd><mo>...</mo></mtd><mtd><msubsup><mi>x</mi><mrow><mi>n</mi><mn>1</mn></mrow><mo>*</mo></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>x</mi><mrow><mn>1</mn><mi>2</mi></mrow><mo>*</mo></msubsup></mtd><mtd><msubsup><mi>x</mi><mi>22</mi><mo>*</mo></msubsup></mtd><mtd><mo>...</mo></mtd><mtd><msubsup><mi>x</mi><mrow><mi>n</mi><mn>2</mn></mrow><mo>*</mo></msubsup></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><msubsup><mi>x</mi><mrow><mn>1</mn><mi>m</mi></mrow><mo>*</mo></msubsup></mtd><mtd><msubsup><mi>x</mi><mrow><mn>2</mn><mi>m</mi></mrow><mo>*</mo></msubsup></mtd><mtd><mo>...</mo></mtd><mtd><msubsup><mi>x</mi><mrow><mi>n</mi><mi>m</mi></mrow><mo>*</mo></msubsup></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open = "(" close = ")"><mtable><mtr><mtd><msubsup><mi>x</mi><mn>11</mn><mi>R</mi></msubsup></mtd><mtd><msubsup><mi>x</mi><mrow><mi>2</mi><mn>1</mn></mrow><mi>R</mi></msubsup></mtd><mtd><mo>...</mo></mtd><mtd><msubsup><mi>x</mi><mrow><mi>n</mi><mn>1</mn></mrow><mi>R</mi></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>x</mi><mrow><mn>1</mn><mi>2</mi></mrow><mi>R</mi></msubsup></mtd><mtd><msubsup><mi>x</mi><mi>22</mi><mi>R</mi></msubsup></mtd><mtd><mo>...</mo></mtd><mtd><msubsup><mi>x</mi><mrow><mi>n</mi><mn>2</mn></mrow><mi>R</mi></msubsup></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><msubsup><mi>x</mi><mrow><mn>1</mn><mi>m</mi></mrow><mi>R</mi></msubsup></mtd><mtd><msubsup><mi>x</mi><mrow><mn>2</mn><mi>m</mi></mrow><mi>R</mi></msubsup></mtd><mtd><mo>...</mo></mtd><mtd><msubsup><mi>x</mi><mrow><mi>n</mi><mi>m</mi></mrow><mi>R</mi></msubsup></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000869066830000012.GIF" wi="1691" he="311" /></maths>步骤五,对X<sup>R</sup>进行主元分析,计算SPE检验中超过控制限的点的数量q,计算误报率C:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>C</mi><mo>=</mo><mfrac><mi>q</mi><mi>n</mi></mfrac><mo>&times;</mo><mn>100</mn><mi>%</mi></mrow>]]></math><img file="FDA0000869066830000013.GIF" wi="294" he="117" /></maths>步骤六,利用遗传算法优化相对转换矩阵Λ,得到最优的相对转换矩阵Λ<sup>*</sup>;选取误报率为遗传算法中的适应度值,按步骤五中的误报率C计算式计算适应度值;遗传算法的终止条件有两个:第一个终止条件是:适应度值在连续两次迭代中均低于6%;第二个终止条件是:连续迭代次数达到K次;在迭代过程中,当满足任一终止条件时,停止迭代,结束相对转换矩阵的优化,此时对应的最优解即为最优的相对转换矩阵Λ<sup>*</sup>;步骤七,利用最优的相对转换矩阵Λ<sup>*</sup>,对标准化后样本矩阵X<sup>*</sup>进行相对化转换,得到最优的相对化样本矩阵X<sup>RZ</sup>:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msup><mi>X</mi><mrow><mi>R</mi><mi>Z</mi></mrow></msup><mo>=</mo><msup><mi>&Lambda;</mi><mo>*</mo></msup><mo>&times;</mo><msup><mi>X</mi><mo>*</mo></msup><mo>=</mo><mfenced open = "(" close = ")"><mtable><mtr><mtd><msubsup><mi>&lambda;</mi><mn>1</mn><mo>*</mo></msubsup></mtd><mtd><mn>0</mn></mtd><mtd><mo>...</mo></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msubsup><mi>&lambda;</mi><mn>2</mn><mo>*</mo></msubsup></mtd><mtd><mo>...</mo></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mo>...</mo></mtd><mtd><msubsup><mi>&lambda;</mi><mi>m</mi><mo>*</mo></msubsup></mtd></mtr></mtable></mfenced><mo>&times;</mo><mfenced open = "(" close = ")"><mtable><mtr><mtd><msubsup><mi>x</mi><mn>11</mn><mo>*</mo></msubsup></mtd><mtd><msubsup><mi>x</mi><mrow><mi>2</mi><mn>1</mn></mrow><mo>*</mo></msubsup></mtd><mtd><mo>...</mo></mtd><mtd><msubsup><mi>x</mi><mrow><mi>n</mi><mn>1</mn></mrow><mo>*</mo></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>x</mi><mrow><mn>1</mn><mi>2</mi></mrow><mo>*</mo></msubsup></mtd><mtd><msubsup><mi>x</mi><mi>22</mi><mo>*</mo></msubsup></mtd><mtd><mo>...</mo></mtd><mtd><msubsup><mi>x</mi><mrow><mi>n</mi><mn>2</mn></mrow><mo>*</mo></msubsup></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><msubsup><mi>x</mi><mrow><mn>1</mn><mi>m</mi></mrow><mo>*</mo></msubsup></mtd><mtd><msubsup><mi>x</mi><mrow><mn>2</mn><mi>m</mi></mrow><mo>*</mo></msubsup></mtd><mtd><mo>...</mo></mtd><mtd><msubsup><mi>x</mi><mrow><mi>n</mi><mi>m</mi></mrow><mo>*</mo></msubsup></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open = "(" close = ")"><mtable><mtr><mtd><msubsup><mi>x</mi><mn>11</mn><mrow><mi>R</mi><mi>Z</mi></mrow></msubsup></mtd><mtd><msubsup><mi>x</mi><mrow><mi>2</mi><mn>1</mn></mrow><mi>R</mi></msubsup></mtd><mtd><mo>...</mo></mtd><mtd><msubsup><mi>x</mi><mrow><mi>n</mi><mn>1</mn></mrow><mrow><mi>R</mi><mi>Z</mi></mrow></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>x</mi><mrow><mn>1</mn><mi>2</mi></mrow><mrow><mi>R</mi><mi>Z</mi></mrow></msubsup></mtd><mtd><msubsup><mi>x</mi><mi>22</mi><mi>R</mi></msubsup></mtd><mtd><mo>...</mo></mtd><mtd><msubsup><mi>x</mi><mrow><mi>n</mi><mn>2</mn></mrow><mrow><mi>R</mi><mi>Z</mi></mrow></msubsup></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><msubsup><mi>x</mi><mrow><mn>1</mn><mi>m</mi></mrow><mrow><mi>R</mi><mi>Z</mi></mrow></msubsup></mtd><mtd><msubsup><mi>x</mi><mrow><mn>2</mn><mi>m</mi></mrow><mrow><mi>R</mi><mi>Z</mi></mrow></msubsup></mtd><mtd><mo>...</mo></mtd><mtd><msubsup><mi>x</mi><mrow><mi>n</mi><mi>m</mi></mrow><mrow><mi>R</mi><mi>Z</mi></mrow></msubsup></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000869066830000021.GIF" wi="1747" he="311" /></maths>步骤八,对最优的相对化样本矩阵X<sup>RZ</sup>进行主元分析,实现铝电解槽况故障诊断。
地址 401331 重庆市沙坪坝区大学城东路20号