发明名称 基于非高斯成分提取和支持向量描述的过程监控方法
摘要 本发明公开了一种基于非高斯成分提取和支持向量描述的过程监控方法,该方法包括:读入训练数据和待诊断数据,数据预处理,建立主元分析模型,微粒群优化算法,非高斯投影计算,支持向量数据描述,残差分析,主元估计,故障检测及模型更新。通过本发明提出的方法,能自动地从工业过程运行数据中提取非高斯成份,从而避免了传统统计过程监控方法假设数据服从正态分布的不足;并且,基于微粒群优化算法的非高斯投影算法,保证了提取的独立成分非高斯性的最大化,避免了独立成分分析方法易陷入局部最优解的问题。与传统的统计过程监控方法相比,本发明提出的方法能够及时发现异常,有效减小误报率,取得良好监控效果。
申请公布号 CN101403923A 申请公布日期 2009.04.08
申请号 CN200810122086.X 申请日期 2008.10.31
申请人 浙江大学 发明人 许仙珍;谢磊;王树青
分类号 G05B23/02(2006.01)I 主分类号 G05B23/02(2006.01)I
代理机构 杭州求是专利事务所有限公司 代理人 张法高
主权项 1、一种基于非高斯成分提取和支持向量描述的过程监控方法,其特征在于包括以下步骤:1)读取生产过程正常运行时关键变量的数据作为训练样本TX;2)对训练样本TX进行预处理,使得各变量的均值为0,方差为1,得到输入矩阵X∈R<sup>N×n</sup>,步骤为:(1)计算均值:<maths num="0001"><![CDATA[<math><mrow><mover><mi>TX</mi><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>T</mi><msub><mi>X</mi><mi>i</mi></msub></mrow></math>]]></maths>(2)计算方差:<maths num="0002"><![CDATA[<math><mrow><msubsup><mi>&sigma;</mi><mi>x</mi><mn>2</mn></msubsup><mo>=</mo><mfrac><mn>1</mn><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msup><mrow><mo>(</mo><msub><mi>TX</mi><mi>i</mi></msub><mo>-</mo><mover><mi>TX</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mn>2</mn></msup></mrow></math>]]></maths>(3)白化处理:<maths num="0003"><![CDATA[<math><mrow><mi>X</mi><mo>=</mo><mfrac><mrow><mi>TX</mi><mo>-</mo><mover><mi>TX</mi><mo>&OverBar;</mo></mover></mrow><msubsup><mi>&sigma;</mi><mi>x</mi><mn>2</mn></msubsup></mfrac></mrow></math>]]></maths>其中,TX为训练样本,N为训练样本数,n为变量数;3)建立主元分析模型;4)基于微粒群优化算法的非高斯投影计算,提取数据中的非高斯成分;5)基于支持向量数据描述,构建非高斯信号的统计变量和控制限;求取非高斯信号分布的超球面,求解以下二次规划问题:<maths num="0004"><![CDATA[<math><mrow><munder><mi>max</mi><mi>&alpha;</mi></munder><munder><mi>&Sigma;</mi><mi>i</mi></munder><msub><mi>&alpha;</mi><mi>i</mi></msub><mo>&lt;</mo><msub><mi>x</mi><mi>i</mi></msub><mo>&CenterDot;</mo><msub><mi>x</mi><mi>i</mi></msub><mo>></mo><mo>-</mo><munder><mi>&Sigma;</mi><mi>i</mi></munder><munder><mi>&Sigma;</mi><mi>j</mi></munder><msub><mi>&alpha;</mi><mi>i</mi></msub><msub><mi>&alpha;</mi><mi>j</mi></msub><mo>&lt;</mo><msub><mi>x</mi><mi>i</mi></msub><mo>&CenterDot;</mo><msub><mi>x</mi><mi>j</mi></msub><mo>></mo></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><munder><mi>&Sigma;</mi><mi>i</mi></munder><msub><mi>&alpha;</mi><mi>i</mi></msub><mo>=</mo><mn>1</mn><mo>,</mo><msub><mi>&alpha;</mi><mi>i</mi></msub><mo>&Element;</mo><mo>[</mo><mn>0</mn><mo>,</mo><mi>C</mi><mo>]</mo><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><mi>N</mi><mo>.</mo></mrow></math>]]></maths>得到超球面的中心<maths num="0006"><![CDATA[<math><mrow><mi>a</mi><mo>=</mo><munder><mi>&Sigma;</mi><mi>i</mi></munder><msub><mi>&alpha;</mi><mi>i</mi></msub><msub><mi>x</mi><mi>i</mi></msub></mrow></math>]]></maths>和半径:<maths num="0007"><![CDATA[<math><mrow><msup><mi>R</mi><mn>2</mn></msup><mo>=</mo><mo>&lt;</mo><msub><mi>x</mi><mi>k</mi></msub><mo>&CenterDot;</mo><msub><mi>x</mi><mi>k</mi></msub><mo>></mo><mo>-</mo><mn>2</mn><munder><mi>&Sigma;</mi><mi>i</mi></munder><msub><mi>&alpha;</mi><mi>i</mi></msub><mo>&lt;</mo><msub><mi>x</mi><mi>k</mi></msub><mo>&CenterDot;</mo><msub><mi>x</mi><mi>i</mi></msub><mo>></mo><mo>+</mo><munder><mi>&Sigma;</mi><mi>i</mi></munder><munder><mi>&Sigma;</mi><mi>j</mi></munder><msub><mi>&alpha;</mi><mi>i</mi></msub><msub><mi>&alpha;</mi><mi>j</mi></msub><mo>&lt;</mo><msub><mi>x</mi><mi>i</mi></msub><mo>&CenterDot;</mo><msub><mi>x</mi><mi>j</mi></msub><mo>></mo><mo>,</mo></mrow></math>]]></maths>x<sub>i</sub>、x<sub>j</sub>为非高斯成分的样本点,x<sub>k</sub>为超球面边界上的支持向量;6)主元估计:构建主元高斯信号的T<sup>2</sup>统计量,计算控制限;当检验水平为α,控制限按下式计算:<maths num="0008"><![CDATA[<math><mrow><msubsup><mi>T</mi><mi>&alpha;</mi><mn>2</mn></msubsup><mo>=</mo><mfrac><mrow><mi>l</mi><mrow><mo>(</mo><msup><mi>N</mi><mn>2</mn></msup><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow><mrow><mi>N</mi><mo>-</mo><mi>l</mi></mrow></mfrac><msub><mi>F</mi><mrow><mi>l</mi><mo>,</mo><mi>N</mi><mo>-</mo><mi>l</mi><mo>;</mo><mi>&alpha;</mi></mrow></msub><mo>,</mo></mrow></math>]]></maths>l为主元个数。7)残差分析:构建残差高斯信号Q统计量,计算控制限;对于任一输入残差e<sub>i</sub>,Q统计量为:<maths num="0009"><![CDATA[<math><mrow><mi>Q</mi><mo>=</mo><msup><mrow><mo>|</mo><mo>|</mo><msub><mi>e</mi><mi>i</mi></msub><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><mo>=</mo><msubsup><mi>e</mi><mi>i</mi><mi>T</mi></msubsup><msub><mi>e</mi><mi>i</mi></msub><mo>;</mo></mrow></math>]]></maths>当检验水平为α,控制限按下式计算:<maths num="0010"><![CDATA[<math><mrow><msubsup><mi>&delta;</mi><mi>&alpha;</mi><mn>2</mn></msubsup><mo>=</mo><msubsup><mi>g&chi;</mi><mrow><mi>h</mi><mo>;</mo><mi>&alpha;</mi></mrow><mn>2</mn></msubsup><mo>,</mo></mrow></math>]]></maths>其中<maths num="0011"><![CDATA[<math><mrow><mi>g</mi><mo>=</mo><mfrac><msup><mi>&rho;</mi><mn>2</mn></msup><mrow><mn>2</mn><mi>&mu;</mi></mrow></mfrac><mo>,</mo></mrow></math>]]></maths><maths num="0012"><![CDATA[<math><mrow><mi>h</mi><mo>=</mo><mfrac><msup><mrow><mn>2</mn><mi>&mu;</mi></mrow><mn>2</mn></msup><msup><mi>&rho;</mi><mn>2</mn></msup></mfrac><mo>,</mo></mrow></math>]]></maths>ρ和μ分别为Q统计量的方差和均值。8)读取生产过程中最新的变量数据作为待诊断数据VX;9)故障检测;10)定期将过程状态正常的点添加到训练集TX中,重复2)~7)的训练过程,以便及时更新支持向量描述、残差分析和主元统计等模型。
地址 310027浙江省杭州市浙大路38号