发明名称 一种高炉炼铁专家系统决策合成方法
摘要 本发明公开了一种高炉炼铁专家系统决策合成方法。以多种类型高炉专家系统的单独决策为输入,构造了四种决策合成算子模型,采用线性和非线性回归技术对决策合成算子模型进行参数辨识,以高炉炉况特征向量及各类型高炉专家系统的历史表现为选择依据,把最优化决策合成算子选择问题转化为一个分类器构造问题,采用非线性多分类支持向量机模型进行分类器建模,通过快速分类器训练算法得到最优分类器,进而利用该分类器对决策合成算子模型参数进行递推式的二次校正,由最优分类器动态选择出适合于当前炉况的最优决策合成算子形成高炉操作优化决策。对高炉冶炼过程的智能控制与决策的合成具有普遍的通用性,能提高高炉专家系统决策准确度与系统鲁棒性。
申请公布号 CN100583136C 申请公布日期 2010.01.20
申请号 CN200710164605.4 申请日期 2007.12.21
申请人 浙江大学 发明人 吴铁军;李艳君;于玲;朱学其;曾建明;杜方;崔承刚;孙丽丽;周微
分类号 G06N5/00(2006.01)I 主分类号 G06N5/00(2006.01)I
代理机构 杭州求是专利事务所有限公司 代理人 张法高
主权项 1.一种高炉炼铁专家系统决策合成方法,其特征在于包括如下步骤:1)以多种类型高炉专家系统的单独决策为输入,构造四种应用于决策值为连续实数的决策合成算子模型,即朴素选择决策合成算子、简单平均决策合成算子、加权平均决策合成算子和非线性回归决策合成算子;2)采用线性回归和非线性回归技术分别辨识出上述四种决策合成算子的模型参数;3)将最优决策合成算子选择问题转化为多分类器问题;4)以多分类支持向量机为分类器模型,采用快速SMO算法为训练方法求解出最优分类器;5)利用最优分类器对决策合成算子的模型参数进行递推式的二次校正,得到最终的决策合成算子模型;6)依次利用最优分类器和最终决策合成算子进行决策合成,获得高炉操作优化决策;所述的朴素选择决策合成算子:选择其中一个专家系统决策为合成后的决策,假设有M个高炉专家系统参与决策合成,它们单独给出的决策值分别为d<sub>1</sub>,d<sub>2</sub>,K,d<sub>M</sub>,合成后的决策记为h,那么共有M种具体的决策合成算子,记为S<sub>1</sub>{h=d<sub>1</sub>},S<sub>2</sub>{h=d<sub>2</sub>},K,S<sub>M</sub>{h=d<sub>M</sub>},其中S<sub>i</sub>{h=d<sub>i</sub>},i=1,2,K M,表示决策合成算子S<sub>i</sub>按h=d<sub>i</sub>的方式产生合成后的决策<img file="C2007101646050002C1.GIF" wi="168" he="88" />所述的简单平均决策合成算子:以多种类型高炉专家系统决策的平均值为合成后的决策,假设有M个高炉专家系统参与决策合成,它们单独给出的决策值分别为d<sub>1</sub>,d<sub>2</sub>,K,d<sub>M</sub>,合成后的决策记为h,记该决策合算子为<maths num="0001"><![CDATA[<math><mrow><msub><mi>S</mi><mrow><mi>M</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>{</mo><mi>h</mi><mo>=</mo><mfrac><mn>1</mn><mi>M</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>d</mi><mi>i</mi></msub><mo>}</mo><mo>,</mo></mrow></math>]]></maths>表示决策合成算子S<sub>M+1</sub>按<maths num="0002"><![CDATA[<math><mrow><mi>h</mi><mo>=</mo><mfrac><mn>1</mn><mi>M</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>d</mi><mi>i</mi></msub></mrow></math>]]></maths>的方式产生合成后的决策;所述的加权平均决策合成算子:以多种类型高炉专家系统决策的加权平均值为合成后的决策,假设有M个高炉专家系统参与决策合成,M个高炉专家系统单独给出的决策值分别为d<sub>1</sub>,d<sub>2</sub>,K,d<sub>M</sub>,合成后的决策记为h,记该决策合成算子为:<maths num="0003"><![CDATA[<math><mrow><msub><mi>S</mi><mrow><mi>M</mi><mo>+</mo><mn>2</mn></mrow></msub><mo>{</mo><mi>h</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>&lambda;</mi><mi>i</mi></msub><msub><mi>d</mi><mi>i</mi></msub><mo>}</mo><mo>,</mo></mrow></math>]]></maths>表示决策合成算子S<sub>M+2</sub>按<maths num="0004"><![CDATA[<math><mrow><mi>h</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>&lambda;</mi><mi>i</mi></msub><msub><mi>d</mi><mi>i</mi></msub></mrow></math>]]></maths>的方式产生合成后的决策,其中λ<sub>i</sub>>0,i=1,2,K M是各专家系统决策的权重分配,且满足<maths num="0005"><![CDATA[<math><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>&lambda;</mi><mi>i</mi></msub><mo>=</mo><mn>1</mn><mo>;</mo></mrow></math>]]></maths>所述的非线性回归决策合成算子:以多种类型高炉专家系统决策为回归项构造非线性回归模型,并以回归模型输出值作为合成后的决策,假设有M个高炉专家系统参与决策合成,它们单独给出的决策值分别为d<sub>1</sub>,d<sub>2</sub>,K,d<sub>M</sub>,合成后的决策记为h,记该决策合成算子为S<sub>M+3</sub>{h=f(d<sub>1</sub>,d<sub>2</sub>,K,d<sub>M</sub>)},其中f是非线性函数,采用非线性回归支持向量机模型来表示f,记该决策合成算子为<maths num="0006"><![CDATA[<math><mrow><msub><mi>S</mi><mrow><mi>M</mi><mo>+</mo><mn>3</mn></mrow></msub><mo>{</mo><mi>h</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>V</mi></munderover><msub><mi>&beta;</mi><mi>i</mi></msub><msub><mi>K</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>d</mi><mn>1</mn></msub><mo>,</mo><msub><mi>d</mi><mn>2</mn></msub><mo>,</mo><mi>K</mi><mo>,</mo><msub><mi>d</mi><mi>M</mi></msub><mo>)</mo></mrow><mo>}</mo><mo>,</mo></mrow></math>]]></maths>表示决策合成算子S<sub>M+3</sub>按<maths num="0007"><![CDATA[<math><mrow><mi>h</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>V</mi></munderover><msub><mi>&beta;</mi><mi>i</mi></msub><msub><mi>K</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>d</mi><mn>1</mn></msub><mo>,</mo><msub><mi>d</mi><mn>2</mn></msub><mo>,</mo><mi>K</mi><mo>,</mo><msub><mi>d</mi><mi>M</mi></msub><mo>)</mo></mrow></mrow></math>]]></maths>的方式产生合成后的决策,其中V是支持向量的个数,β<sub>i</sub>是核函数系数,K<sub>i</sub>是与第i个支持向量对应的非线性核函数;所述的将最优决策合成算子选择问题转化为多分类器问题:将步骤2)中辨识后的四种决策合成算子模型分别应用于训练样本Z={Z<sub>1</sub>,Z<sub>2</sub>,K Z<sub>N</sub>},对于每组样本Z<sub>i</sub>={d<sub>i1</sub>,d<sub>i2</sub>,K,d<sub>iM</sub>,y<sub>i</sub>},得出各种决策合成算子产生的决策合成值{h<sub>i1</sub>,h<sub>i2</sub>,K,h<sub>iM</sub>,h<sub>i,M+1</sub>,h<sub>i,M+2</sub>,h<sub>i,M+3</sub>},以真实决策值y<sub>i</sub>为参考标准,可确认出针对样本Z<sub>i</sub>的最优决策合成算子<img file="C2007101646050003C6.GIF" wi="85" he="69" />其数学表达式为<img file="C2007101646050003C7.GIF" wi="675" he="165" />采用交叉验证方法选择出高炉炉况特征向量为输入向量,记训练样本Z={Z<sub>1</sub>,Z<sub>2</sub>,K Z<sub>N</sub>}对应的炉况特征向量数据为X={X<sub>1</sub>,X<sub>2</sub>,K X<sub>N</sub>},X<sub>i</sub>=[x<sub>i1</sub>,x<sub>i2</sub>,K,x<sub>iQ</sub>],其中Q是高炉特征向量的维数<img file="C2007101646050003C8.GIF" wi="127" he="198" />据式<maths num="0008"><![CDATA[<math><mrow><msub><mi>o</mi><mi>i</mi></msub><mo>=</mo><munder><mrow><mi>arg</mi><mi>min</mi></mrow><mrow><mi>m</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mi>K</mi><mo>,</mo><mi>M</mi><mo>+</mo><mn>3</mn></mrow></munder><mrow><mo>(</mo><mo>|</mo><msub><mi>h</mi><mi>im</mi></msub><mo>-</mo><msub><mi>y</mi><mi>i</mi></msub><mo>|</mo><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>以[x<sub>i1</sub>,x<sub>i2</sub>,K,x<sub>iQ</sub>]为输入向量,以o<sub>i</sub>为正确的类别标签,形成用于分类问题的训练样本C={C<sub>1</sub>,C<sub>2</sub>,K C<sub>N</sub>},C<sub>i</sub>={x<sub>i1</sub>,x<sub>i2</sub>,K x<sub>iQ</sub>,o<sub>i</sub>},至此,将最优决策合成算子选择问题转化为多分类器问题;所述的利用最优分类器对决策合成算子的模型参数进行递推式的二次校正:将步骤4)的最优分类器,F:X→{S<sub>1</sub>,S<sub>2</sub>,K,S<sub>M</sub>,S<sub>M+1</sub>,S<sub>M+2</sub>,S<sub>M+3</sub>},应用于炉况数据X={X<sub>1</sub>,X<sub>2</sub>,K X<sub>N</sub>},根据映射F将X划分了M+3个子集,记为X<sup>1</sup>,X<sup>2</sup>,K,X<sup>M</sup>,X<sup>M+1</sup>,X<sup>M+2</sup>,X<sup>M+3</sup>,由X={X<sub>1</sub>,X<sub>2</sub>,K X<sub>N</sub>}与决策合成算子模型原始训练样本Z={Z<sub>1</sub>,Z<sub>2</sub>,K Z<sub>N</sub>}的一一对应关系,F亦将Z划分成M+3个子集,记为Z<sup>1</sup>,Z<sup>2</sup>,K,Z<sup>M</sup>,Z<sup>M+1</sup>,Z<sup>M+2</sup>,Z<sup>M+3</sup>,与加权平均决策合成算子S<sub>M+2</sub>对应的样本子集是Z<sup>M+2</sup>,由于Z<sup>M+2</sup>对S<sub>M+2</sub>模型参数辨识更具指导意义,因此增强样本子集Z<sup>M+2</sup>的信度,对S<sub>M+2</sub>模型进行二次校正,以S<sub>M+2</sub>为初始模型,以Z<sup>M+2</sup>为需要强调的更新样本,采用带遗忘因子的递推最小二乘辨识算法完成对S<sub>M+2</sub>的二次校正,修正后的模型记为S′<sub>M+2</sub>,同理,与非线性回归决策合成算子S<sub>M+3</sub>对应的样本子集是Z<sup>M+3</sup>,对S<sub>M+3</sub>模型进行二次校正,以S<sub>M+3</sub>为初始模型,以Z<sup>M+3</sup>为需要更新样本,采用带遗忘因子的递推支非线性回归持向量机训练算法完成对S<sub>M+3</sub>的二次校正,修正后的模型记为S′<sub>M+3</sub>;所述的依次利用最优分类器和最终决策合成算子形成高炉操作优化决策:以当前炉况数据特征向量X′=[x′<sub>1</sub>,x′<sub>2</sub>,K,x′<sub>Q</sub>]为最优分类器的输入向量,其中Q是炉况数据特征向量的维数,利用式F:X→{S<sub>1</sub>,S<sub>2</sub>,K,S<sub>M</sub>,S<sub>M+1</sub>,S<sub>M+2</sub>,S<sub>M+3</sub>}最优分类选出最适合于当前炉况的决策合成算子S<sub>o</sub>∈{S<sub>1</sub>,S<sub>2</sub>,K,S<sub>M</sub>,S<sub>M+1</sub>,S′<sub>M+2</sub>,S′<sub>M+3</sub>},被选出的最终决策合成算子S<sub>o</sub>以各类型高炉专家系统单独决策值[d′<sub>1</sub>,d′<sub>2</sub>,K,d′<sub>M</sub>]为模型的输入,产生合成后的高炉操作优化决策h′=S<sub>o</sub>(d′<sub>1</sub>,d′<sub>2</sub>,K,d′<sub>M</sub>)。
地址 310027浙江省杭州市浙大路38号