发明名称 基于特征子空间优化相对矩阵的铝电解槽况诊断方法
摘要 一种基于特征子空间优化相对矩阵的铝电解槽况诊断方法,其特征在于:一,采集原始测量样本集,对该原始测量样本集进行预处理后投影到核空间;二,对中心化矩阵K进行相对主元分析,建立铝电解槽况诊断模型,并对铝电解槽况进行诊断;三,通过细菌觅食算法在搜索区域内寻得最优相对转换矩阵;四,按步骤二所述方法,利用最优相对转换矩阵建立铝电解槽况诊断模型,实现对铝电解槽况的准确诊断。本发明充分考虑铝电解槽况非线性的特性,通过核函数,将非线性参数投影到高维线性特征空间,在核空间内借助细菌觅食法优化相对转换矩阵,进行相对主元分析,大大提高了铝电解槽况故障诊断的精确度。
申请公布号 CN104499001A 申请公布日期 2015.04.08
申请号 CN201510027742.8 申请日期 2015.01.20
申请人 重庆科技学院 发明人 易军;黄迪;李太福;周伟;张元涛;姚立忠;田应甫
分类号 C25C3/20(2006.01)I 主分类号 C25C3/20(2006.01)I
代理机构 重庆为信知识产权代理事务所(普通合伙) 50216 代理人 余锦曦
主权项 一种基于特征子空间优化相对矩阵的铝电解槽况诊断方法,其特征在于:包括下列步骤:步骤一,采集原始测量样本集,对该原始测量样本集进行预处理后投影到核空间,包括:第一步:采集n组铝电解槽况数据组成原始测量样本集<img file="FDA0000658823820000011.GIF" wi="691" he="113" />每个样本含有m个独立的铝电解槽况参数采样值;第二步:对原始测量样本集X<sup>0</sup>进行标准化处理,得到标准化后样本矩阵X;第三步:利用核函数,将标准化后样本矩阵X投影到高维特征空间后得到矩阵K<sup>0</sup>;第四步:对矩阵K<sup>0</sup>进行中心化处理,得到中心化矩阵K,中心化处理按下式进行:K=K<sup>0</sup>‑I<sub>n</sub>K<sup>0</sup>‑K<sup>0</sup>I<sub>n</sub>+I<sub>n</sub>K<sup>0</sup>I<sub>n</sub>其中,<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>I</mi><mi>n</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>n</mi></mfrac><msub><mfenced open='(' close=')'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mi>K</mi></mtd><mtd><mn>1</mn></mtd></mtr><mtr><mtd><mi>M</mi></mtd><mtd><mi>O</mi></mtd><mtd><mi>M</mi></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mi>L</mi></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mrow><mi>n</mi><mo>&times;</mo><mi>n</mi></mrow></msub></mrow>]]></math><img file="FDA0000658823820000012.GIF" wi="565" he="307" /></maths>步骤二,对中心化矩阵K进行相对主元分析,建立铝电解槽况诊断模型,并对铝电解槽况进行诊断,包括:第一步:在[0,50)范围内随机产生相对转换矩阵Λ,所述相对转换矩阵Λ为对角矩阵:<maths num="0002" id="cmaths0002"><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><mi>&Lambda;</mi></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msub><mi>&lambda;</mi><mn>2</mn></msub></mtd><mtd><mi>&Lambda;</mi></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mi>M</mi></mtd><mtd><mi>M</mi></mtd><mtd><mi>M</mi></mtd><mtd><mi>M</mi></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mi>&Lambda;</mi></mtd><mtd><msub><mi>&lambda;</mi><mi>m</mi></msub></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000658823820000013.GIF" wi="492" he="309" /></maths>第二步:对中心化矩阵K进行相对化转换,得到相对化样本矩阵K<sup>R</sup>=K·Λ,即:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msup><mi>K</mi><mi>R</mi></msup><mo>=</mo><mfenced open='(' close=')'><mtable><mtr><mtd><msub><mi>&lambda;</mi><mn>1</mn></msub><mi>k</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>,</mo><msub><mi>x</mi><mn>1</mn></msub><mo>)</mo></mrow></mtd><mtd><mi>K</mi></mtd><mtd><msub><mi>&lambda;</mi><mn>1</mn></msub><mi>k</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>,</mo><msub><mi>x</mi><mi>m</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>M</mi></mtd><mtd><mi>O</mi></mtd><mtd><mi>M</mi></mtd></mtr><mtr><mtd><msub><mi>&lambda;</mi><mi>m</mi></msub><mi>k</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>m</mi></msub><mo>,</mo><msub><mi>x</mi><mn>1</mn></msub><mo>)</mo></mrow></mtd><mtd><mi>L</mi></mtd><mtd><msub><mi>&lambda;</mi><mi>m</mi></msub><mi>k</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>m</mi></msub><mo>,</mo><msub><mi>x</mi><mi>m</mi></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open='(' close=')'><mtable><mtr><mtd><msup><mi>k</mi><mo>*</mo></msup><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>,</mo><msub><mi>x</mi><mn>1</mn></msub><mo>)</mo></mrow></mtd><mtd><mi>K</mi></mtd><mtd><msup><mi>k</mi><mo>*</mo></msup><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>,</mo><msub><mi>x</mi><mi>m</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>M</mi></mtd><mtd><mi>O</mi></mtd><mtd><mi>M</mi></mtd></mtr><mtr><mtd><msup><mi>k</mi><mo>*</mo></msup><mrow><mo>(</mo><msub><mi>x</mi><mi>m</mi></msub><mo>,</mo><msub><mi>x</mi><mn>1</mn></msub><mo>)</mo></mrow></mtd><mtd><mi>L</mi></mtd><mtd><msup><mi>k</mi><mo>*</mo></msup><mrow><mo>(</mo><msub><mi>x</mi><mi>m</mi></msub><mo>,</mo><msub><mi>x</mi><mi>m</mi></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000658823820000021.GIF" wi="1334" he="228" /></maths>第三步:对K<sup>R</sup>进行主元分析,计算<img file="FDA0000658823820000022.GIF" wi="56" he="70" />检验统计量及对应控制限<img file="FDA0000658823820000023.GIF" wi="99" he="77" />对铝电解槽况进行诊断;<img file="FDA0000658823820000024.GIF" wi="60" he="73" />检验统计量按下式计算:<img file="FDA0000658823820000025.GIF" wi="424" he="185" />其中,SPE和SPE<sub>0</sub>分别表示SPE检验统计量和对应的控制限;T<sup>2</sup>和<img file="FDA0000658823820000026.GIF" wi="78" he="92" />分别表示T<sup>2</sup>检验统计量和对应的控制限;若<img file="FDA0000658823820000027.GIF" wi="222" he="76" />则认为铝电解槽况正常;若<img file="FDA0000658823820000028.GIF" wi="218" he="72" />则认为铝电解槽况异常;第四步:统计铝电解槽况诊断情况,记录诊断错误的样本点个数q,按下式计算漏检率C:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>C</mi><mo>=</mo><mfrac><mi>q</mi><mi>n</mi></mfrac><mo>&times;</mo><mn>100</mn><mo>%</mo></mrow>]]></math><img file="FDA0000658823820000029.GIF" wi="301" he="119" /></maths>步骤三,通过细菌觅食算法在搜索区域内寻得最优相对转换矩阵,具体包括以下步骤:第一步:初始化细菌觅食算法相关参数:细菌群体大小BIOsize=50,趋向次数Nc=100,趋向行为执行中前进次数NS=4,繁殖次数Nre=6,驱散次数Ned=4,执行驱散行为的概率Ped=0.25;第二步:将诊断模型中的漏检率C作为评价函数,利用细菌觅食算法在给定范围内优化相对转换矩阵,当漏检率到达最低时,得到的相对转换矩阵即为最优;步骤四,按步骤二所述方法,利用最优相对转换矩阵建立铝电解槽况诊断模型,实现对铝电解槽况的准确诊断。
地址 401331 重庆市沙坪坝区大学城东路20号