发明名称 最优模糊网络的工业熔融指数软测量仪表及方法
摘要 本发明公开了一种最优模糊网络的工业熔融指数软测量仪表及方法。该方法通过引入支持向量机对原有的模糊神经网络进行优化,解决了模糊神经网络构建过程中参数难设定的问题。在本发明中,现场智能仪表、控制站与DCS数据库连接,软测量值显示仪包括最优模糊网络的工业熔融指数软测量模型,DCS数据库与软测量模型的输入端连接,所述最优模糊网络的工业熔融指数软测量模型的输出端与熔融指数软测量值显示仪连接。最后,本发明具有在线测量、计算速度快、抗噪声能力强、推广性能好的特点。
申请公布号 CN103675005B 申请公布日期 2015.09.30
申请号 CN201310432289.X 申请日期 2013.09.22
申请人 浙江大学 发明人 刘兴高;张明明;李见会
分类号 G01N25/04(2006.01)I;G05B13/04(2006.01)I 主分类号 G01N25/04(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 周烽
主权项 一种最优模糊网络的工业熔融指数软测量仪表,包括用于测量易测变量的现场智能仪表、用于测量操作变量的控制站、存放数据的DCS数据库以及熔融指数软测量值显示仪,所述现场智能仪表、控制站与DCS数据库连接,其特征在于:所述软测量仪表还包括最优模糊网络的工业熔融指数软测量模型,所述DCS数据库与所述最优模糊网络的工业熔融指数软测量模型的输入端连接,所述最优模糊网络的工业熔融指数软测量模型的输出端与熔融指数软测量值显示仪连接,所述最优模糊网络的工业熔融指数软测量模型包括:数据预处理模块,用于将从DCS数据库输入的模型训练样本进行预处理,使得训练样本的均值为0,方差为1,该处理采用以下算式过程来完成:计算均值:<maths num="0001" id="cmaths0001"><math><![CDATA[<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><msub><mi>TX</mi><mi>i</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000729114440000011.GIF" wi="1430" he="148" /></maths>计算方差:<maths num="0002" id="cmaths0002"><math><![CDATA[<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><mrow><mo>(</mo><msub><mi>TX</mi><mi>i</mi></msub><mo>-</mo><mover><mi>TX</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000729114440000012.GIF" wi="1405" he="142" /></maths>标准化:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>X</mi><mo>=</mo><mfrac><mrow><mi>TX</mi><mo>-</mo><mover><mi>TX</mi><mo>&OverBar;</mo></mover></mrow><msub><mi>&sigma;</mi><mi>x</mi></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000729114440000013.GIF" wi="1467" he="162" /></maths>其中,TX<sub>i</sub>为第i个训练样本,N为训练样本数,<img file="FDA0000729114440000014.GIF" wi="85" he="70" />为训练样本的均值,X为标准化后的训练样本;σ<sub>x</sub>表示训练样本的标准差,σ<sup>2</sup><sub>x</sub>表示训练样本的方差;模糊神经网络模块,对从数据预处理模块传过来的输入变量,进行模糊推理和建立模糊规则;对从数据预处理模块传过来的经过预处理过的训练样本X进行模糊分类,得到模糊规则库中每个模糊聚类的中心和宽度;设第p个标准化后的训练样本X<sub>p</sub>=[X<sub>p1</sub>,...,X<sub>pn</sub>],其中n是输入变量的个数;设模糊神经网络有R个模糊规则,为了求得每个模糊规则对于训练样本X<sub>p</sub>的每个输入变量X<sub>pj</sub>,j=1,…,n,下面的模糊化方程将求出其对第i个模糊规则的隶属度:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>M</mi><mi>ij</mi></msub><mo>=</mo><mi>exp</mi><mo>{</mo><mo>-</mo><mfrac><msup><mrow><mo>(</mo><msub><mi>X</mi><mi>pj</mi></msub><mo>-</mo><msub><mi>m</mi><mi>ij</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><msubsup><mi>&sigma;</mi><mi>ij</mi><mn>2</mn></msubsup></mfrac><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000729114440000015.GIF" wi="1680" he="213" /></maths>其中m<sub>ij</sub>和σ<sub>ij</sub>分别表示第i个模糊规则的第j个高斯成员函数的中心和宽度,由模糊聚类求得;设标准化后的训练样本X<sub>p</sub>对模糊规则i的适应度为μ<sup>(i)</sup>(X<sub>p</sub>),则μ<sup>(i)</sup>(X<sub>p</sub>)的大小可由下式决定:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msup><mi>&mu;</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msup><mrow><mo>(</mo><msub><mi>X</mi><mi>p</mi></msub><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Pi;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>M</mi><mi>ij</mi></msub><mrow><mo>(</mo><msub><mi>X</mi><mi>p</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mo>{</mo><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mfrac><msup><mrow><mo>(</mo><msub><mi>X</mi><mi>pj</mi></msub><mo>-</mo><msub><mi>m</mi><mi>ij</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><msubsup><mi>&sigma;</mi><mi>ij</mi><mn>2</mn></msubsup></mfrac><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000729114440000021.GIF" wi="1655" he="205" /></maths>求得输入训练样本对于每个规则的适应度之后,模糊神经网络对模糊规则输出进行推导以得到最后的解析解;在常用的模糊神经网络结构中,每个模糊规则推导的过程都可以表示为如下:首先求得训练样本中所有输入变量的线性乘积和,然后用此线性乘积和与规则的适用度μ<sup>i</sup>(X<sub>p</sub>)相乘,得到最终的每条模糊规则的输出;模糊规则i的推导输出可以表示如下:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msup><mi>f</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msup><mo>=</mo><msup><mi>&mu;</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msup><mrow><mo>(</mo><msub><mi>X</mi><mi>p</mi></msub><mo>)</mo></mrow><mo>&times;</mo><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>a</mi><mi>ij</mi></msub><mo>&times;</mo><msub><mi>X</mi><mi>pj</mi></msub><mo>+</mo><msub><mi>a</mi><mrow><mi>i</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000729114440000022.GIF" wi="1647" he="168" /></maths><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mover><mi>y</mi><mo>^</mo></mover><mi>p</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>R</mi></munderover><msup><mi>f</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msup><mo>+</mo><mi>b</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>R</mi></munderover><mo>[</mo><msup><mi>&mu;</mi><mi>i</mi></msup><mrow><mo>(</mo><msub><mi>X</mi><mi>p</mi></msub><mo>)</mo></mrow><mo>&times;</mo><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>a</mi><mi>ij</mi></msub><mo>&times;</mo><msub><mi>X</mi><mi>pj</mi></msub><mo>+</mo><msub><mi>a</mi><mrow><mi>i</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow><mo>]</mo><mo>+</mo><mi>b</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000729114440000023.GIF" wi="1673" he="170" /></maths>式中,f<sup>(i)</sup>为第i条模糊规则的输出,<img file="FDA0000729114440000024.GIF" wi="62" he="75" />是模糊神经网络模型对第p个训练样本的预测输出,a<sub>ij</sub>,j=1,…,n是第i条模糊规则中第j个变量的线性系数,a<sub>i0</sub>是第i条模糊规则中输入变量线性乘积和的常数项,b是输出偏置量;支持向量机优化模块,在式(7)中,输入变量线性乘积和中的参数的确定是模糊神经网络使用中用到的一个主要问题,这里我们采用把原有的模糊规则推导输出形式转换为支持向量机优化问题,再使用支持向量机进行线性优化,转换过程如下:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mfenced open='' close=''><mtable><mtr><mtd><msub><mover><mi>y</mi><mo>^</mo></mover><mi>p</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>R</mi></munderover><msup><mi>f</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msup><mo>+</mo><mi>b</mi></mtd></mtr><mtr><mtd><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>R</mi></munderover><mo>[</mo><msup><mi>&mu;</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msup><mrow><mo>(</mo><msub><mi>X</mi><mi>p</mi></msub><mo>)</mo></mrow><mo>&times;</mo><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>a</mi><mi>ij</mi></msub><mo>&times;</mo><msub><mi>X</mi><mi>pj</mi></msub><mo>+</mo><msub><mi>a</mi><mrow><mi>i</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow><mo>]</mo><mo>+</mo><mi>b</mi></mtd></mtr><mtr><mtd><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>R</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mi>n</mi></munderover><msub><mi>a</mi><mi>ij</mi></msub><mo>&times;</mo><msup><mi>&mu;</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msup><mrow><mo>(</mo><msub><mi>X</mi><mi>p</mi></msub><mo>)</mo></mrow><mo>&times;</mo><msub><mi>X</mi><mi>pj</mi></msub><mo>+</mo><mi>b</mi></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000729114440000025.GIF" wi="1664" he="463" /></maths>其中X<sub>p0</sub>为常数项且恒等于1;令<img file="FDA0000729114440000026.GIF" wi="1683" he="104" />其中,<img file="FDA0000729114440000027.GIF" wi="161" he="100" />表示原训练样本的转化形式,即把原来的训练样本转换为如上式形式,作为支持向量机的训练样本:<img file="FDA0000729114440000028.GIF" wi="1675" he="106" />其中y<sub>1</sub>,…,y<sub>N</sub>是训练样本的目标输出,取S作为新的输入训练样本集合,那么原有问题可以转化为如下的支持向量机对偶优化问题:<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><mi>R</mi><mrow><mo>(</mo><mi>&omega;</mi><mo>,</mo><mi>b</mi><mo>)</mo></mrow><mo>=</mo><mi>&gamma;</mi><mfrac><mn>1</mn><mi>N</mi></mfrac><msubsup><mi>&Sigma;</mi><mrow><mi>p</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></msubsup><msub><mi>L</mi><mi>&epsiv;</mi></msub><mrow><mo>(</mo><msub><mi>y</mi><mi>p</mi></msub><mo>,</mo><mi>f</mi><mrow><mo>(</mo><msub><mi>X</mi><mi>p</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mi>&omega;</mi><mi>T</mi></msup><mi>&omega;</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000729114440000031.GIF" wi="1672" he="131" /></maths>其中y<sub>p</sub>是输入训练样本X<sub>p</sub>的目标输出,ω是支持向量机超平面的法向量,f(X<sub>p</sub>)是对应于X<sub>p</sub>的模型输出,γ是支持向量机的惩罚因子,R(ω,b)是优化问题的目标函数,N是训练样本数,L<sub>ε</sub>(y<sub>p</sub>,f(X<sub>p</sub>))表达式如下:<img file="FDA0000729114440000032.GIF" wi="1679" he="221" />其中ε是优化问题的误差容限,接下来使用支持向量机求得模糊神经网络的模糊规则最优推导线性参数和对偶优化问题的预报输出:<maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><msub><mi>a</mi><mi>ij</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mrow><mo>(</mo><msubsup><mi>&alpha;</mi><mi>k</mi><mo>*</mo></msubsup><mo>-</mo><msub><mi>&alpha;</mi><mi>k</mi></msub><mo>)</mo></mrow><msup><mi>&mu;</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msup><msub><mi>X</mi><mi>kj</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>&Element;</mo><mi>SV</mi></mrow><mi>N</mi></munderover><mrow><mo>(</mo><msubsup><mi>&alpha;</mi><mi>k</mi><mo>*</mo></msubsup><mo>-</mo><msub><mi>&alpha;</mi><mi>k</mi></msub><mo>)</mo></mrow><msup><mi>&mu;</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msup><msub><mi>X</mi><mi>kj</mi></msub><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>R</mi><mo>;</mo><mi>j</mi><mo>=</mo><mn>0</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000729114440000033.GIF" wi="1663" he="147" /></maths><img file="FDA0000729114440000034.GIF" wi="1657" he="144" />其中α<sub>k</sub>,<img file="FDA0000729114440000035.GIF" wi="79" he="80" />(k=1,…,N)分别是y<sub>p</sub>‑f(X<sub>p</sub>)大于0和小于0时对应的拉格朗日乘子,<img file="FDA0000729114440000036.GIF" wi="62" he="81" />即为对应于第p个标准化后的训练样本X<sub>p</sub>的MI预报值;所述最优模糊网络的工业熔融指数软测量模型还包括:模型更新模块,用于模型的在线更新,定期将离线化验数据输入到训练集中,更新模糊神经网络模型。
地址 310058 浙江省杭州市西湖区余杭塘路866号