发明名称 一种基于鲁棒随机权神经网络的多元铁水质量软测量方法
摘要 一种基于鲁棒随机权神经网络的多元铁水质量软测量方法属于高炉冶炼自动化控制领域,特别涉及一种基于Cauchy分布加权M估计随机权神经网络(M-RVFLNs)的高炉炼铁过程多元铁水质量参数动态软测量方法。本发明运用主成分分析(PCA)方法筛选出影响高炉铁水质量的最主要参数作为模型输入变量,构造出一个具有输出自反馈结构且考虑不同时刻输入输出数据的铁水质量多元动态预测模型,可同时对表征高炉铁水质量的主要参数Si含量、P含量、S含量和铁水温度进行多元动态软测量。本发明包括以下步骤:(1)辅助变量选择与模型输入变量确定;(2)M-RVFLNs软测量模型的训练和使用。
申请公布号 CN105608492A 申请公布日期 2016.05.25
申请号 CN201610118914.7 申请日期 2016.03.02
申请人 东北大学 发明人 周平;吕友彬;王宏
分类号 G06N3/10(2006.01)I;C21B5/00(2006.01)I;C21B7/24(2006.01)I 主分类号 G06N3/10(2006.01)I
代理机构 沈阳亚泰专利商标代理有限公司 21107 代理人 许宇来
主权项 一种基于鲁棒随机权神经网络的多元铁水质量软测量方法,其特征在于包括以下步骤:(1)辅助变量选择与模型输入变量确定选择软测量的辅助变量包括:炉腹煤气量u<sub>1</sub>(m<sup>3</sup>)、热风温度u<sub>2</sub>(℃)、热风压力u<sub>3</sub>(KPa)、富氧率u<sub>4</sub>、鼓风湿度u<sub>5</sub>(RH)、喷煤量u<sub>6</sub>(m<sup>3</sup>/h);确定如下16个变量为软测量模型的输入变量:当前时刻炉腹煤气量u<sub>1</sub>(t)(m<sup>3</sup>);当前时刻热风温度u<sub>2</sub>(t)(℃);当前时刻热风压力u<sub>3</sub>(t)(KPa);当前时刻富氧率u<sub>4</sub>(t);当前时刻鼓风湿度u<sub>5</sub>(t)(RH);当前时刻设定喷煤量u<sub>6</sub>(t)(m<sup>3</sup>/h).上一时刻炉腹煤气量u<sub>1</sub>(t‑1)(m<sup>3</sup>);上一时刻热风温度u<sub>2</sub>(t‑1)(℃);上一时刻热风压力u<sub>3</sub>(t‑1)(KPa);上<sub>一</sub>时刻富氧率u<sub>4</sub>(t‑1);上一时刻鼓风湿度u<sub>5</sub>(t‑1)(RH);上一时刻设定喷煤量u<sub>6</sub>(t‑1)(m<sup>3</sup>/h);上一时刻Si含量估计值<img file="FDA0000933342570000011.GIF" wi="263" he="71" />上一时刻P含量估计值<img file="FDA0000933342570000012.GIF" wi="271" he="71" />上一时刻S含量估计值<img file="FDA0000933342570000013.GIF" wi="271" he="69" />上一时刻铁水温度估计值<img file="FDA0000933342570000014.GIF" wi="294" he="70" />(2)M‑RVFLNs软测量模型的训练和使用(A)开始:所有变量初始化;(B)若选择为M‑RVFLNs模型训练,转至(C)读取需要进行模型训练的数据集;若选择为铁水质量参数软测量,转至(I)调取已经训练完的高炉铁水质量参数M‑RVFLNs预报模型;(C)读取模型训练所需数据集:从数据库中读取或输入模型训练学习所需数据集Z={(x<sub>i</sub>,t<sub>i</sub>)|x<sub>i</sub>∈R<sup>n</sup>,t<sub>i</sub>=R<sup>m</sup>,i=1,…N},N≥L进行模型学习初始化,N是样本数目,L是隐含层节点的数目,x<sub>i</sub>为输入数据集,t<sub>i</sub>为输出数据集;n为输入数据集维度,m为输出数据集维度;(D)数据预处理:对数据进行缺失值填充,之后进行归一化处理后,作为最终预测模型的训练数据;(E)模型相关待定参数确定:M‑RVFLNs模型需要预先设定的待定参数包括:激活函数类型g,隐含层节点数目L,收敛条件E;(F)M‑RVFLNs模型初始化训练及模型参数确定:基于模型训练样本集以及预先设定的模型相关待定参数,进行模型的学习和训练;模型的训练和学习包括两个阶段,初始化阶段(initialization phase)和鲁棒学习阶段(robust learning phase);模型初始化阶段的具体步骤如下:步骤(a):从给定的数据集O={(x<sub>i</sub>,t<sub>i</sub>)|x<sub>i</sub>∈R<sup>n</sup>,t<sub>i</sub>∈R<sup>m</sup>,i=1,…}中选取训练数据集Z={(x<sub>i</sub>,t<sub>i</sub>)|x<sub>i</sub>∈R<sup>n</sup>,t<sub>i</sub>=R<sup>m</sup>,i=1,…N},N≥L,包括输入数据<img file="FDA0000933342570000021.GIF" wi="670" he="71" /><img file="FDA0000933342570000022.GIF" wi="436" he="70" />和目标数据Y={y<sub>j</sub>(t)|j=1,2,…,m},初始化隐层节点为L,设置鲁棒训练收敛条件E;步骤(b):随机选取输入权值向量a<sub>i</sub>和隐层单元偏差b<sub>i</sub>,i=1,…,L;步骤(c):计算隐含层输出矩阵H<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>H</mi><mo>=</mo><msub><mfenced open = '[' close = ']'><mtable><mtr><mtd><mrow><mi>G</mi><mrow><mo>(</mo><msub><mi>a</mi><mn>1</mn></msub><mo>,</mo><msub><mi>b</mi><mn>1</mn></msub><mo>,</mo><msub><mi>x</mi><mn>1</mn></msub><mo>)</mo></mrow></mrow></mtd><mtd><mn>...</mn></mtd><mtd><mrow><mi>G</mi><mrow><mo>(</mo><msub><mi>a</mi><mi>L</mi></msub><mo>,</mo><msub><mi>b</mi><mi>L</mi></msub><mo>,</mo><msub><mi>x</mi><mn>1</mn></msub><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mrow></mrow></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mn>...</mn></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mrow></mrow></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mrow><mi>G</mi><mrow><mo>(</mo><msub><mi>a</mi><mn>1</mn></msub><mo>,</mo><msub><mi>b</mi><mn>1</mn></msub><mo>,</mo><msub><mi>x</mi><mi>N</mi></msub><mo>)</mo></mrow></mrow></mtd><mtd><mn>...</mn></mtd><mtd><mrow><mi>G</mi><mrow><mo>(</mo><msub><mi>a</mi><mi>L</mi></msub><mo>,</mo><msub><mi>b</mi><mi>L</mi></msub><mo>,</mo><msub><mi>x</mi><mi>N</mi></msub><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mrow><mi>N</mi><mo>&times;</mo><mi>L</mi></mrow></msub></mrow>]]></math><img file="FDA0000933342570000023.GIF" wi="821" he="239" /></maths>步骤(d):按下述公式计算输出权重β:<img file="FDA0000933342570000024.GIF" wi="462" he="71" />步骤(e):计算隐含层节点数为L时的模型输出:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mover><mi>Y</mi><mo>~</mo></mover><mo>=</mo><mi>H</mi><mo>&CenterDot;</mo><mi>&beta;</mi></mrow>]]></math><img file="FDA0000933342570000025.GIF" wi="190" he="78" /></maths>步骤(f):计算建模残差r:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>r</mi><mo>=</mo><mi>T</mi><mo>-</mo><mover><mi>Y</mi><mo>~</mo></mover></mrow>]]></math><img file="FDA0000933342570000026.GIF" wi="204" he="62" /></maths>模型鲁棒学习阶段的具体步骤如下:步骤(g):计算标准化残差向量<img file="FDA0000933342570000031.GIF" wi="546" he="133" />其中median()为中值计算公式;步骤(h):将标准化残差代入到Cauchy加权函数<img file="FDA0000933342570000032.GIF" wi="426" he="118" />中得到m维数据对应的权值矩阵<img file="FDA0000933342570000033.GIF" wi="589" he="118" />进一步求得权函数方阵W<sup>(0)</sup>,其中μ=median(r<sub>i</sub>),<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>a</mi><mo>=</mo><mn>1</mn><mo>/</mo><msqrt><mrow><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msup><mrow><mo>(</mo><msub><mi>r</mi><mi>i</mi></msub><mo>-</mo><mover><msub><mi>r</mi><mi>i</mi></msub><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt></mrow>]]></math><img file="FDA0000933342570000034.GIF" wi="389" he="134" /></maths>步骤(i):根据引入M估计后输出权重的迭代公式<img file="FDA0000933342570000035.GIF" wi="551" he="69" />其中k=1,2,…为迭代次数,迭代计算求取<img file="FDA0000933342570000036.GIF" wi="391" he="71" />若每个参数估计值<img file="FDA0000933342570000037.GIF" wi="332" he="94" />都小于指定收敛条件E,则迭代计算终止,得最终输出权重<img file="FDA0000933342570000038.GIF" wi="196" he="71" />(G)建模效果评估:引入均方根误差RMSE:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>R</mi><mi>M</mi><mi>S</mi><mi>E</mi><mo>=</mo><msqrt><mrow><mfrac><mn>1</mn><mi>N</mi></mfrac><mrow><mo>(</mo><mi>H</mi><mo>(</mo><mrow><msub><mi>a</mi><mi>i</mi></msub><mo>,</mo><msub><mi>b</mi><mi>i</mi></msub><mo>,</mo><msub><mi>x</mi><mi>i</mi></msub></mrow><mo>)</mo><msup><mi>&beta;</mi><mo>*</mo></msup><mo>-</mo><mi>Y</mi><mo>)</mo></mrow><msup><mrow><mo>(</mo><mi>H</mi><mo>(</mo><mrow><msub><mi>a</mi><mi>i</mi></msub><mo>,</mo><msub><mi>b</mi><mi>i</mi></msub><mo>,</mo><msub><mi>x</mi><mi>i</mi></msub></mrow><mo>)</mo><msup><mi>&beta;</mi><mo>*</mo></msup><mo>-</mo><mi>Y</mi><mo>)</mo></mrow><mi>T</mi></msup></mrow></msqrt></mrow>]]></math><img file="FDA0000933342570000039.GIF" wi="1030" he="143" /></maths>对建模误差进行综合评价,如果建模误差符合实际工况标准,则结束本次M‑RVFLNs模型训练学习过程,转(H);若误差不符合预定标准,重新训练,转(E);(H)保存M‑RVFLNs模型:模型训练学习结束,以得到的M‑RVFLNs模型作为高炉铁水质量参数预报模型;(I)读取M‑RVFLNs模型:调出之前初始训练好的高炉铁水质量参数M‑RVFLNs软测量模型<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mover><mi>Y</mi><mo>~</mo></mover><mi>i</mi></msub><mo>=</mo><mi>H</mi><mrow><mo>(</mo><msub><mi>a</mi><mi>i</mi></msub><mo>,</mo><msub><mi>b</mi><mi>i</mi></msub><mo>,</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow><msup><mi>&beta;</mi><mo>*</mo></msup><mo>;</mo></mrow>]]></math><img file="FDA00009333425700000310.GIF" wi="431" he="79" /></maths>(J)读取模型预报的过程数据;(K)判断数据是否异常或者缺失;判断模型16个输入数据是否有数据缺失情况;若有则转(L)进行数据处理,否则转(M)进行铁水质量参数预报;(L)数据处理:若有数据缺失情况,则用前一时刻相应变量数据进行替换;(M)铁水质量参数预报:将输入变量数据归一化处理后,调用之前训练好的M‑RVFLNs模型进行铁水质量参数预报;(N)铁水质量参数预报结果显示:在预报系统人机界面上显示本次铁水质量参数预报的结果;(O)数据保存:将本次软测量的相关输入输出数据保存到相应的历史数据库,供后续系统评估、修正以及查询所用。
地址 110819 辽宁省沈阳市和平区文化路3号巷11
您可能感兴趣的专利