CN104899463B

CN104899463B - 高炉铁水硅含量四分类趋势预报模型的建立方法及应用

Info

Publication number: CN104899463B
Application number: CN201510340996.5A
Authority: CN
Inventors: 蒋朝辉; 尹菊萍; 桂卫华; 阳春华; 谢永芳
Original assignee: Central South University
Current assignee: Central South University
Priority date: 2015-06-18
Filing date: 2015-06-18
Publication date: 2017-11-21
Anticipated expiration: 2035-06-18
Also published as: CN104899463A

Abstract

本发明涉及一种高炉铁水硅含量四分类趋势预报模型的建立方法及应用，属于自动化检测技术领域。所述模型的建立方法包括采集历史数据；进行数据预处理；确定模型输入变量和输出变量；依据模糊均值聚类的方法对铁水硅含量样本有效的聚类，获取模型输出变量四类趋势变化区间的划分标准；利用极限学习机建立四分类趋势预报模型。利用本发明所建立的模型不仅可以预报硅含量变化的趋势，同时可以得到趋势变化的幅度大小，即可以预报得到下一炉铁水硅含量是大幅上升、小幅上升、大幅下降、小幅下降的四分类趋势变化情况。这对高炉操作者提早判断炉况，并采取小调、早调等措施，避免炉况发生急剧变化有很好的指导。

Description

高炉铁水硅含量四分类趋势预报模型的建立方法及应用

技术领域

本发明涉及一种高炉冶炼过程中铁水硅含量趋势预报模型的建立方法及应用，属于自动化检测技术领域。

背景技术

铁水硅含量是表征高炉炼铁过程炉温及其变化趋势的关键信息，也是反应铁水质量、能耗等指标的重要物理量。但铁水硅含量及其变化趋势无法直接在线检测，导致对炉况调控不及时或者盲目，致使炉温大幅波动、炉况不顺，若炉温过低，则铁水物理热不足，炉缸热储备不够，不仅铁水质量差，且容易造成严重的炉缸冻结事故；若炉温过高，则炉内煤气流太过旺盛，导致悬料、崩料等事故，严重影响高炉的焦比、排放。因此研究如何实时预报铁水硅含量及其变化趋势，对稳定高炉热状态、减少炉况的波动、提高生铁质量和降低焦比等具有重要意义。

目前，对高炉铁水硅含量的预测可分为两大方向，一是利用检测得到的大量数据建立基于数据驱动的数值预测模型，也称为黑箱模型。常用的黑箱模型有：时间序列预测模型、神经网络预测模型、混沌模型、SVM预测模型、非线性动力学理论预测模型等。二是通过将高炉铁水硅含量的数值转变为上升和下降的趋势变化，进而建立铁水硅含量的二分类趋势预报。目前，关于铁水硅含量趋势预报方面的文献还比较少，而实际应用中，高炉工厂一般根据经验判断炉温的变化趋势，进而进行操作调节。

基于数据驱动的预测模型以铁水硅含量为高炉炉温的衡量标准，研究高炉冶炼过程参数与硅含量的函数关系，是一种仅依赖于系统输入输出的黑箱模型方法，对高炉铁水硅含量的预报有一定的效果，但是仍存在自身局限性，如时间序列模型更适应于炉况平稳的情形；神经网络能够获得比时间序列更好的预测效果，但是在实际应用中学习时间长且容易产生拟合现象；SVM的预测模型具有很好的泛化性能，在大多数情况下能获得优于其他方法的预测效果。

然而，上述的数据驱动模型主要集中在铁水硅含量的数值预测上，很少关注趋势的变化。尽管数值预测给出了铁水硅含量直接的结果，但是数值预测结果有时候会对铁水硅含量趋势的变化无法做出诊断或者做出错误的诊断。例如，在连续的两炉铁水中硅含量数值的变化很小，甚至比评价是否成功命中的临界值都小，这种情况下没有必要去构建数值预测模型，因为单一的预测模型完全依据目标命中率来评判，而此时的命中率虽然显示成功命中，却无法表针硅含量的这种微小变化，从这一点也可以得到，预测高炉炉缸热状态，也就是预测高炉铁水硅含量的趋势变化至关重要。

目前，对铁水硅含量的趋势变化的研究主要集中在二分类预报问题(即铁水硅含量上升和下降)，但是二分类的预报只能得到趋势变化的方向，无法得到这种方向变化的幅度。

中国专利申请公布号CN 103160629 A，申请公布日2013.06.19，公开了一种预报高炉趋热的方法，其综合考虑了整个炼铁工艺中影响铁水硅含量的因素。通过建立三个预测模型，分别为：高炉铁水中的Si含量的RBF神经网络计算模型、高炉料速计算模型、理论燃烧温度计算模型，然后利用这三个模型对高炉趋热结果的判断，最终得到高炉是否趋热。但是该方法仅仅可以预测高炉是否趋热，对于其他变化趋势无法得到判断，比如高炉是否趋冷，或者趋热趋冷的程度等等，都无法得到准确的信息。

中国专利申请公布号CN103160626 A，申请公布日2013.06.19，公开了一种判断高炉炉缸过凉的方法，该方法与CN 103160629 A类似，不同点在于该方法揭示了一种判断高炉炉缸过凉的方法。分别从物理传热、化学传热、历史数据等影响高炉向凉发展的角度出发，建立了三个预测模型，分别为：高炉铁水中的Si含量的RBF神经网络计算模型、高炉料速计算模型、理论燃烧温度计算模型。综合三个模型所得结果判断高炉炉缸是否过凉。但是该方法也是只能判断一种炉温变化情况，无法获取炉温是否趋热、以及趋冷趋热程度的变化信息。另外，三个模型得到的数值对铁水硅含量影响的权重判断无法准确划分，这对预报结果的命中率有一定影响。

中国专利申请公布号CN10211383A，申请公布日2008.7.2，公开了一种高炉铁水硅含量的特征分析预报方法，该方法建立了一种高炉铁水硅含量的预报模型，采取改进的动态独立成分分析方法对输入变量的样本数据进行特征提取，提取了测量对象的高阶统计信息，通过小样本建模设计的基于遗传算法优化的最小二乘支持向量机方法建立铁水硅含量预报回归模型。但是该方法只能预测得到下一炉铁水硅含量的数值，无法准确判断炉温的变化，更无法获取炉温趋势的变化程度，可以反馈给高炉操作者的信息少，可靠性不高。

中国专利申请公布号CN102031319A，申请公布日20011.04.27，公开了一种高炉铁水硅含量的预报方法，该方法将炉况分为小波动和大波动两种情况分别建模，分别对炉况稳定和波动都进行了预报，当炉况波动较小时主要依靠均线系统，在炉况波动大时自动加入风口理论燃烧温度以及前次铁水含硫量，达到预测铁水硅含量的目的。该方法虽然综合考虑了炉况稳定和异常两种情况，但是在实际应用中，无法提前预判何时出现异常炉况，进而也就无法及时更换预测模型，对结果的命中率有一定影响。

综上所述，现有的各种预测炉温的方法，或无法全面预报，或只能预报数值，而对于炉温趋势预报的方法较少。

发明内容

本发明的目的是提出一种高炉铁水硅含量四分类趋势预报模型的建立方法及应用，根据与高炉铁水硅含量有密切关系的可控以及不可控变量的历史数据，利用极限学习机(extreme learning machine，ELM)建立四分类趋势预报模型，优选并采用状态转移算法优化模型参数，进而利用所建模型来预报未来一段时间硅含量的变化趋势。利用本发明所建立的模型不仅可以预报硅含量变化的趋势，同时可以得到趋势变化的幅度大小，即可以预报得到下一炉铁水硅含量是大幅上升、小幅上升、大幅下降、小幅下降的四分类趋势情况。这对高炉操作者提早判断炉况，并采取小调、早调等措施，避免炉况发生急剧变化有很好的指导。

本发明目的是通过以下技术方案实现的，一种高炉铁水硅含量四分类趋势预报模型的建立方法，包括：采集历史数据；进行数据预处理，以确定模型输入变量(即铁水硅含量影响因子)和输出变量；优选结合时滞分析和\或专家经验分析确定模型输入变量；依据模糊均值聚类的方法对铁水硅含量样本有效的聚类，达到对数据集进行正确的分类，获取模型输出变量——铁水硅含量的四类趋势变化区间的划分标准，即大幅上升、小幅上升、大幅下降、小幅下降的评判指标；利用极限学习机建立四分类趋势预报模型即ELM分类器；优选采用状态转移算法优化所述预报模型参数。

具体地，所述高炉铁水硅含量四分类趋势预报模型的建立方法，包括以下步骤：

S1、采集历史数据

所述采集历史数据为采集整个炼铁过程中影响铁水硅含量变化的变量及铁水硅含量数据。

由于高炉内部的复杂物理化学反应，其间接影响铁水硅含量变化的变量有很多，包括上部的布料方式、原料性质，下部的控制参数如风量、风温等。常用的变量数据有富氧率、透气性指数、标准风速、富氧流量、冷风流量、喷煤量、铁水成分、煤气成分、顶压、全压差、热风压力、实际风速、冷风压力、理论燃烧温度、热风温度、鼓风动能、富氧流量、富氧压力、炉腹煤气量等等。

S2、数据预处理

所述数据预处理即通过对现场炼铁工艺了解，对可采集到的历史数据进行分析，获得时间点对应的时间序列样本、剔除异常数据、确定模型输入变量(即与铁水硅含量相关性强的变量)和输出变量、归一化处理及差分处理。具体包括以下步骤：

S21、获得时间点对应的时间序列样本；

采集历史数据时，不同变量有不同的采样周期(或间隔)，因此在数据预处理之前需要折中考虑。以一个变量的时间间隔为参照(例如以采样周期最长的变量的时间间隔，本发明的一种具体实施方式为1h)，依据时间点对应对所有变量进行匹配，最终得到所需时间序列样本。比如：铁水成分(如：Si，Ti，S，P，Mn)的采集周期大概为40min，煤气成分(如：H₂,CO,CO₂)的采集周期为1s，喷煤量的采集周期为30min，而一些控制变量(喷煤量、热风温度、富氧量、风量等)和状态变量(透气性指数、炉顶温度)的采集周期为1h。在建立本发明所述模型时，需要相同时间点对应的各个变量值作为模型的一个输入样本。本发明的一种实施方式中以最长的采集周期1h为参照，通过时间点对应对所有变量进行匹配，最终提取时间序列样本。

S22、剔除异常数据；

在采样过程中，受到高温高压等环境影响或者高炉休风、减风等非正常状况，数据存在异常值。这些异常数据在一定程度上会改变数据的变化趋势，影响模型建立的准确性，因此需要对异常数据进行剔除。

本发明通过3σ准则进行异常值(和缺失值)处理，剔除异常数据，即若残差|V_i|＞3σ，则该数据为异常数据，该样本为异常样本，予以剔除。重复使用该方法剔除异常数据，直至没有异常数据为止。

其中：

x_i为所处理的数据，N为样本个数。为x_i的残差。

S23、确定模型输入变量和输出变量；

对所有采集到的影响铁水硅含量变化的变量与铁水硅含量之间进行相关性分析，以与铁水硅含量的变化相关性强的变量，作为本发明所述模型的输入变量；以铁水硅含量作为输出变量。由于高炉内部的复杂物理化学反应，其间接影响铁水硅含量的因素有很多，包括上部的布料方式、原料性质，下部的控制参数如风量、风温等。本发明在建立所述模型时，将影响铁水硅含量的变化强的变量作为模型输入变量，将铁水硅含量作为模型输出变量。铁水中硅含量的变化与固体原料，包含铁矿石、烧结矿和焦炭；待加热的气态物质，包含空气和一些辅助燃料；以及下部的风量、风温等参数的变化有密切的关系。表1列出了20个待选变量。过多的输入变量会增加模型的复杂度，而过少的输入变量又会降低模型精度。因此，在建立模型时有必要删除待选变量中与铁水硅含量相关性不强的变量。本发明在一种具体实施方式中，通过SPSS(Statistical Product and Service Solutions，即“统计产品与服务解决方案”软件)，计算了所有待选变量与铁水硅含量的相关系数，结果如表1所示。

值得说明的是，虽然透气性指数与硅含量的相关性系数不高，但在实际高炉调控中，它是现场操作者判断铁水硅含量变化的重要依据之一，因此，本发明在建模时也将其作为模型输入变量。

同时，由于高炉自身特点所决定，高炉的工艺参数对高炉炉温[Si]具有时滞影响，即存在一定的滞后性，及存在滞后时间。所谓的滞后时间指的是某一参数的变化将在n小时或n分钟之后才会导致高炉炉温[Si]发生变化。那么所谓的滞后时间也可以这样理解，某一时刻的高炉炉温[Si]受n小时或n分钟以前的工艺参数的影响，对应着那一段时间的工艺参数，相关性是最大的，也是最密切的。因此，要确定相关变量的滞后时间，需要借助于高炉现场专家人员的经验及相关性分析的手段。具体地，本发明通过SPSS进行时滞分析，得到上一炉铁水硅含量与本炉铁水硅含量有很强的相关性，因此，将上一炉铁水硅含量也作为待选变量进行选择。本发明在建模时也将时滞(即上一炉铁水硅含量)作为模型输入变量。同时，专家经验具有重要参考价值，也是值得考虑的重要因素。优选地，本发明同时结合时滞分析和/或专家经验确定模型输入变量。综上，在本发明的一种具体实施方式中，结合相关性分析、时滞分析以及专家经验得到一种优选的方案是选取富氧率、透气性指数、鼓风动能、全压差、热风温度、喷煤量、利用系数、冷风流量以及上一炉铁水硅含量共9个变量作为模型的输入变量。

表1输入变量与硅含量相关性系数

S24、归一化处理；

即将各输入变量分别归一化处理；优选地，将输出变量(即铁水硅含量)也进行归一化处理。因选取的各输入变量量纲不同，对模型的收敛速度和复杂度有很大影响，在建模前需要对其分别进行归一化处理。通过数值变换来消除变量间的量纲影响。方法如下：

其中x_i，分别为第i个变量归一化前、后取值，max(x_i)，min(x_i)分别为第i个变量归一化前的最大值、最小值。

S25、差分处理；

即分别对各输入变量和输出变量进行差分处理(例如进行前向差分)得到所有变量的变化率；这样处理可以更好的与趋势变化相吻合。前向差分公式如下：

Δx_i＝x_i-x_i-1,i＝1,2,…,n (Ⅲ)

经差分处理分别得到所有输入变量的变化率，组成输入样本；得到所有输出变量(即铁水硅含量)的变化率，组成输出样本。

S3、数据分类

所述数据分类即将所述输出样本依据模糊均值聚类算法划分为四类区间；具体包括以下步骤：

S31、以零为边界将所述输出样本划分为两大类，大于零的样本划分为第一类，简称为正类样本，小于零的样本划分为第二类，简称为负类样本。本发明将等于零的样本划分为正类样本。

S32、在所述正类样本中获取两个聚类中心c₁、c₂，在所述负类样本中获取两个聚类c′₁、c′₂，以这四个聚类中心得到两个边界值：

S33、根据步骤S32所述边界值，将所述正类样本和负类样本划分为四类区间：这四类区间依次代表四类变化趋势：大幅下降、小幅下降、小幅上升、大幅上升。

所述S32在正类样本中获取两个聚类中心c₁、c₂，在负类样本中获取两个聚类c′₁、c′₂，具体包括以下步骤：

S321、初始化算法参数；具体方法如下：

初始化聚类个数c、模糊指数m、误差ε、算法迭代次数T，初始化聚类中心V＝[v₁,v₂,...,v_c]。在[0,1]范围内随机初始化隶属度矩阵U，使u_ij满足公式(Ⅳ)约束条件

其中c是聚类个数，n为样本数，x_i为聚类样本，v_i聚类中心向量，u_ij为第i个样本属于第j个中心的隶属度，m为模糊指数，用来控制聚类的模糊程度。

S322、根据式(VII)计算K(x_i,v_j)，式(VIII)计算参数ρ；具体方法如下：

构建模糊C均值聚类算法的目标函数如下：

式(Ⅴ)中，

||φ(x_i)-φ(v_j)||＝K(x_i,x_i)+K(v_j,v_j)-2K(x_i,v_j) (Ⅵ)

K(x_i,v_j)＝exp{-(||x_i-v_j||²)/ρ²} (VII)

依据S321和S322将所述正类样本和所述负类样本初步划分在不同类别中，然而，这种初步划分准确性不高。因此，通过如下公式更新聚类中心和隶属度，为步骤S324的进一步划分提供条件。

S323、依据式(VI)、式(VII)更新聚类中心V和隶属度矩阵U；具体方法如下：

分别对J_φ(U,V)关于U,V求偏导，得到新的聚类中心和隶属度，公式如下：

S324、依据算法目标函数进行算法迭代；即依据S323更新一次聚类中心和隶属度，带入式(V)便得到一个新的目标函数值。

迭代终止的条件如下：

即更新前后目标函数值的差值小于给定误差或者达到给定的最大迭代次数时算法结束，否则继续更新。

算法迭代结束即可以在所述正类样本中获得两个聚类中心c₁、c₂；在所述负类样本中获取两个聚类c′₁、c′₂。

S4、建立预报模型

即依据所述步骤S3对所述输出样本的四类区间划分标准，将所述输出样本进行二进制编码；然后结合所述输入样本利用极限学习机建立预报模型。所述建立预报模型具体包括以下步骤：

S41、编码；

即依据上述步骤S3对对输出样本的四类区间划分标准，将所有输出样本进行二进制编码。所述输出样本的四类区间划分标准为：这四类区间依次代表四类变化趋势：大幅下降、小幅下降、小幅上升、大幅上升。

也就是说，对所有样本(x_j,t_j)的输出t_j进行二进制编码。x_j指输入变量；t_j指输出变量(即铁水硅含量变化率)。对c类分类问题，需要[log₂c]个字节编码c个不同的类。本发明中的四分类问题，输出编码如下表1所示：

表1四分类编码

其中：f2代表趋势变化方向，该列中“-1”代表下降，“1”代表上升；f1代表趋势变化幅度，该列中“-1”代表小幅变化，“1”代表大幅变化。因此得到编码矩阵：

通过如上编码，将所有输出样本t_j(即铁水硅含量)按编码的不同依次划分在各自的类中。例如其中每一个输出值t_j都由两个编码字节组成。

S42、利用极限学习机(ELM)建立预报模型；所述建立好的预报模型也称ELM分类器；具体过程包括：经过上述编码，对于N个不同样本(x_j,t_j)，其中x_j＝[x_j1,x_j2,…,x_jn]^T∈Rⁿ,t_j＝[t_j1,t_j2,…,t_jm]^T∈R^m，m＝log₂c，t_j＝{-1,1}，具有K个隐含层神经元数目，并且激励函数为g(x)的ELM模型可以表示为：

式中ω_i＝[ω_1i,ω_2i,…,ω_ni]是连接输入神经元和第i个隐含层神经元的ELM模型输入权值；b_i是第i个隐元偏差；β_i＝[β_i1,β_i2,…,β_im]^T是连接第i个隐含层神经元和输出神经元的ELM模型输出权值；ω_i·x_j表示ω_i和x_j的内积，激励函数g(x)选用sigmoid函数，即：

最小化：

使得：

其中：h(x_i)＝[g(α₁x_i+b₁),…，g(α_Kx_i+B_K)],ξ_i＝[ξ_i,1,…,ξ_i,m]^T是样本x_i对应输出的训练误差，依据KKT条件(库恩塔克条件，即Karush-Kuhn-Tucker Conditions)，训练ELM等价于解决如下的对偶优化问题：引入拉格朗日乘子α_i(i＝1,2,…n)：

其中，每个拉格朗日乘子α_i对应第i个训练样本。β_j连接隐含层节点到第j个输出节点的向量，β＝[β₁,…,β_m]，相应的KKT的优化条件如下：

将(XVII)、(XVIII)式带入(XIX)式，得到如下方程：

其中：

将(XX)带入(XVII)式，得到隐元偏置：

其中：

通过式(XXII)计算得到的隐层权值β，将结果带入式(XII)，得到(测试样本)对应的输出。

ELM的网络结构图1所示。

优选地，本发明所述高炉铁水硅含量四分类趋势预报模型的建立方法还包括解码的步骤；所述解码是指对建立好的预报模型(即ELM分类器)的输出进行解码，即对步骤S42得到的(测试样本)输出进行解码。

优选地，所述建立预报模型还包括训练和测试的步骤；

所述测试样本就是用来检测模型准确性的样本。也就是说，先用一大批历史数据(称为训练样本，包含输入变量和输出变量，都是已知的)对上述建立好的预报模型进行训练，训练好之后就需要用另一批数据(称为测试样本，也包含输入变量和输出变量)对所训练好的所述模型进行测试。测试的时候只需将输入变量带入所述模型，该模型会得到相应的输出，然后用该模型得到的输出与实际的输出比较，进而得到结果的准确性，如果准确性高，说明所建模型好，可以用于现场对高炉铁水硅含量进行预测，如果不高，则继续修正模型参数。

以下以测试样本为例，对解码步骤进行描述：

如图2所示，测试样本指图2中的②对应的那部分数据，测试样本的输出即由本发明预报模型(上述步骤S42所建立的ELM分类器)计算得到的铁水硅含量变化趋势。

给定测试样本x，依据上述建立好的预报模型(ELM分类器)，得到m＝log₂c维的输出向量：

f(x)＝(f_m(x),…,f₁(x)) (XXIV)

这些输出需要解码才能得到样本x最终所属的类。在此，运用基于损失函数的解码方法对ELM分类器的输出进行解码。在这种解码方法中，如要样本x被标签为i类，样本(x,i)的总损失函数值最小。样本(x,i)的总损失函数定义为：

M(i)指的是编码矩阵c*m的第i个行向量，表示为：M(i)＝(M(i,m),…,M(i,1))。

基于上述损失函数(方程(XXV))，训练样本x最终的类输出为：

式(XXVI)也称决策函数。

依据上述描述得到如下ELM分类过程：

输入：训练集X_train；输出：训练集Y_train；

输入：测试集X_test；输出：测试集Y_test；

Step1：编码过程，依据编码矩阵对训练集类别Y_train进行编码；

Step2：训练过程，得到隐含层输出矩阵H以及连接隐含层节点和输出层节点的权值β，即完成图2中的①部分；

Step3：测试过程，输入测试集，即图2中的②部分，基于上述参数，获得ELM的输出，即图2中的③部分；

Step4：解码过程，将上述输出依据式(XXV)的损失函数以及式(XXVI)的决策函数得到测试集所属的类别。

建立所述预报模型(即ELM分类器)时，由于模型的输入权值和隐元偏差理论上可以随机赋值得到，这无疑是ELM的一个较大的缺陷，使得预报结果出现了很大的不稳定性，对高炉炼铁过程中铁水硅含量的趋势预报带来了很大的影响。因此，如何设置影响算法性能的关键参数，对预测模型的学习精度和泛化能力的好坏起着决定性作用。优选地，本发明用状态转移算法(STA)优化得到输入权值和隐元偏差，不仅可以提高模型稳定性，而且使得模型的预测结果更可靠。具体流程如下：

输入：随机产生ELM的输入权值ω_K*n和隐元偏差β_K*m；

输出：最优ELM输入权值和隐元偏差；

Step1：初始化搜索力度SE，问题维数Dim，搜索范围，迭代次数，极限学习机中的输入权值和隐元偏差以及隐含层节点数K；

Step2:引入STA算法，将初始化得到的输入权值和隐元偏差作为STA的输入样本，均方根误差MSE作为STA算法的目标函数；

Step3:如果算法满足收敛精度或是达到了最大的迭代次数，转step6，否则转step4；

Step4：对STA算法中的样本执行：

Step41扩张转变(expend transformation)

x_k+1＝x_k+γR_cx_k (XXVII)

其中，γ是一个正常数，称为扩张因子；R_c∈R^n*n是一个随机对角矩阵，其中每一个元素都服从高斯分布。扩张转变的功能是将x_k的范围扩张到[-∞,+∞]，即在整个空间搜索。

Step42旋转转变(rotation transformation)

其中，x_k∈Rⁿ，α代表正常数，称为旋转因子。R_r∈R^n*n是一个随机矩阵，服从[-1,1]的均匀分布，||·||₂代表二范数。旋转转变可以在超平面进行搜索。

Step43压缩转变(axesion transformation)

x_k+1＝x_k+δR_ax_k (XXIX)

其中，δ是一个正常数，称为压缩因子；R_a∈R_n*n是一个随机对角矩阵，其中每一个元素都服从高斯分布，并且只有一个随机位置是非零值。压缩转变目的是搜索中心，增强单维搜索。

Step5：判断算法是否满足运行的终止条件，一般终止条件为最大的迭代次数以及算法的精度，若是满足，转step6，否则g＝g+1，转step4；

Step6：输出STA算法的最优值，即ELM的输入权值和隐元偏差。

本发明还包括按上述方法所建立的预报模型在高炉铁水硅含量四分类趋势预报方面的应用。

进一步，本发明还提供一种高炉铁水硅含量四分类趋势预报方法，包括按上述方法建立高炉铁水硅含量四分类趋势预报模型，选取一组变量数据作为输入变量，输入到所述预报模型，得到高炉铁水硅含量四分类趋势预报结果；优选地，选取的输入到所述预报模型的输入变量与建立所述预报模型时采用的输入变量相同；所述变量相同是指变量名称相同，并非指变量数据相同；例如，建立所述预报模型时选取的变量为富氧率、透气型指数、鼓风动能、全压差、热风温度、喷煤、利用系数、冷风流量以及上一炉铁水硅含量共9个变量，那么在利用所述模型进行实际预报时也选取这9个变量，将其具体数据输入所述模进行预报；优选地，按上述方法对所述预报模型的输出进行解码，得到该输出所属的类别，即可得到该输出的四分类趋势。

本发明所述高炉铁水硅含量四分类趋势是指大幅上升、小幅上升、大幅下降、小幅下降。

本发明的关键点

1、本发明通过分析高炉内部的复杂物理化学反应以及高炉炼铁原理，获取了影响铁水硅含量的各种变量在一段时间内的数据，这些变量包括上部的布料方式、原料性质，下部的风温、风量等变量。通过时间点对应选取所需样本，依据拉伊达准则进行异常值处理，通过SPSS分析采集到的各变量与铁水硅含量之间的相关性系数以及时滞影响，得到与铁水硅含量有强相关性的变量作为所建模型的输入变量。

2、本发明采用模糊均值聚类算法确定区间划分依据，首先对历史数据进行聚类划分，形成若干个子样本集，其次对模型输出样本集进行区间划分，最终确定了模型输出变化趋势的四类区间的划分依据。

3.本发明采用二进制代码对上述四类进行编码，得到编码矩阵。用编码区分不同的类，更容易建立模型，也更容易识别输出类别。在对模型输出进行解码时，运用汉明距离来衡量模型输出与编码矩阵之间的相似程度，进而确定模型输出所属的类。

4、采用极限学习机对样本集进行建模，将经过异常值处理后的样本划分为训练集和测试集。为了避免极限学习机在参数选取上的随机性，进而造成结果的不稳定，本发明通过状态转移算法优化极限学习机的输入权值和隐元偏差。通过训练集对模型进行训练，获得铁水硅含量趋势变化幅度的分类器。再通过测试集进行铁水硅含量四分类趋势变化的准确性验证。

本发明的效果

本发明基于采样数据以及炼铁工艺的特点，以每炉次出铁时间为周期，时间点对应的方式选取数据样本，采用SPSS选取模型输入变量，并用模糊均值聚类的方法确定了模型输出变量的四分类趋势变化的控制边界。通过STA-ELM建模方法实现铁水硅含量四分类趋势预测的目的，本发明充分利用高炉可检测到的影响铁水硅含量的数据，使得铁水硅含量的趋势预测结果命中率更高，对高炉操作现场有更好的指导作用。

附图说明

图1为本发明ELM网络结构图。

图2为本发明STA-ELM趋势预测流程图。

图3为本发明模糊均值聚类划分区间图。

图4为本发明STA-ELM模型优化过程示意图。

图5为本发明实施例2样本实际铁水硅含量变化趋势图。

图6为本发明实施例2样本预测铁水硅含量变化趋势图。

图5和图6中：纵轴“1”“2”“3”“4”分别代表大幅下降、小幅下降、大幅上升、小幅上升；横轴代表炉次。

具体实施方式

以下实施例用于说明本发明，但不用来限制本发明的范围。

实施例1

本实施例在某钢厂2650m³高炉进行试验测试。

一种高炉铁水硅含量四分类趋势预报模型的建立方法，具体包括如下步骤：

1)采集历史数据。通过自动化系统的组态软件实现炼铁过程控制及数据采集。所述的自动化系统包括高炉本体、上料系统、热风炉系统、喷煤系统组成。其中来自高炉本体的数据主要由：炉顶压力、热风压力、炉顶温度等相关数据。来自热风炉系统的数据主要有：高炉煤气量、送风、炉顶温度、烟道温度等送风相关数据。来自喷煤系统的数据主要有：喷吹压力、喷吹流量等喷煤有关的数据。

2)数据预处理。将上述采集到的数据进行数据预处理。具体包括：

①获得时间点对应的时间序列样本，即通过时间点对应选取样本。以最长的采集周期1h为参照，通过时间点对应对所有采集到的变量进行时间点匹配，提取时间序列样本，本实施例最终提取931个时间序列样本。

②利用3σ准则删除异常数据。通过异常值处理，本实施例共剔除31个异常数据。

③确定模型输入变量和输出变量。本实施例通过SPSS选取的模型输入变量为：富氧率、透气型指数、鼓风动能、全压差、热风温度、喷煤、利用系数、冷风流量以及上一炉铁水硅含量共9个变量作为模型的输入变量。

④归一化处理。

⑤差分处理。本实施例通过前向差分对所有变量进行差分，得到变化率。所有输入变量的变化率，组成输入样本；所有输出变量的变化率，组成输出样本。

3)模糊均值聚类得到趋势变化区间。通过模糊均值聚类对输出样本(铁水硅含量)进行聚类划分，最终确定四个趋势变化区间为：这四个区间分别对应趋势变化的：大幅下降、小幅下降、

小幅上升、大幅上升。

4)建立预报模型

①编码：对数据预处理之后的输出样本通过编码矩阵进行二进制编码。

②建立预报模型。利用状态转移算法优化极限学习机的输入权值和隐元偏差，建立优化极限学习机的预报模型。将数据预处理后的900组样本划分为训练样本和测试样本，选取前700组作为训练样本，用来训练模型参数，后200组作为测试样本，用来验证模型准确性。通过不断修正模型参数，得到最优的预报模型。

实施例2

本实施例涉及一种利用实施例1所建立高炉铁水硅含量四分类趋势预报模型对高炉铁水硅含量四分类趋势预报方法，具体为选取一组变量数据作为输入变量，输入到所述预报模型，对预报模型的输出结果进行解码，得到最终的硅含量变化趋势(高炉铁水硅含量四分类趋势预报结果)。

具体而言，选取了从2013年1月9号21:00到2013年2月27号10:00总共1166组数据利用实施例1所述数据处理方法进行处理。其中200组测试样本对应的实际铁水硅含量变化趋势如图5所示，由图5可知，变化趋势大多都落在小幅上升和小幅下降的区间内，只有小部分样本落在大幅上升和大幅下降的区间，这也说明炉况比较稳定，符合现场实际情况。利用实施例1所建立的高炉铁水硅含量四分类趋势预报模型进行预报，预报结果如图6所示。对比图5和图6，可见大多炉次趋势预测都准确，表明本实施例预报结果与实际值相吻合，符合实际要求，取得了满意的预报效果。

虽然，上文中已经用一般性说明及具体实施方案对本发明作了详尽的描述，但在本发明基础上，可以对之作一些修改或改进，这对本领域技术人员而言是显而易见的。因此，在不偏离本发明精神的基础上所做的这些修改或改进，均属于本发明要求保护的范围。

Claims

1.一种高炉铁水硅含量四分类趋势预报模型的建立方法，包括以下步骤：

S1、采集历史数据；即采集整个炼铁过程中影响铁水硅含量变化的变量及铁水硅含量数据；

S2、数据预处理；包括获得时间点对应的时间序列样本、剔除异常数据、确定模型输入变量和输出变量、归一化处理及差分处理；经差分处理分别得到所有输入变量的变化率，组成输入样本；得到所有输出变量的变化率，组成输出样本；

S3、数据分类；即将所述输出样本依据模糊均值聚类算法划分为四类区间；

S4、建立预报模型；即依据所述步骤S3对所述输出样本的四类区间划分标准，将所述输出样本进行二进制编码；然后结合所述输入样本利用极限学习机建立预报模型；

所述利用极限学习机建立预报模型的过程包括：将所述输出样本进行二进制编码后，对于N个不同样本(x_j,t_j)，其中x_j＝[x_j1,x_j2,…,x_jn]^T∈Rⁿ,t_j＝[t_j1,t_j2,…,t_jm]^T∈R^m，m＝log₂c，t_j＝{-1,1}，具有K个隐含层神经元数目，并且激励函数为g(x)的ELM模型表示为：

<mrow> <msub> <mi>f</mi> <mi>K</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>K</mi> </munderover> <msub> <mi>&beta;</mi> <mi>i</mi> </msub> <mi>g</mi> <mrow> <mo>(</mo> <msub> <mi>&omega;</mi> <mi>i</mi> </msub> <mo>&CenterDot;</mo> <msub> <mi>x</mi> <mi>j</mi> </msub> <mo>+</mo> <msub> <mi>b</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mi>h</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mi>&beta;</mi> <mo>,</mo> <mi>j</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mi>N</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mi>X</mi> <mi>I</mi> <mi>I</mi> <mo>)</mo> </mrow> </mrow>

<mrow> <mi>g</mi> <mrow> <mo>(</mo> <msub> <mi>&omega;</mi> <mi>i</mi> </msub> <mo>&CenterDot;</mo> <msub> <mi>x</mi> <mi>j</mi> </msub> <mo>+</mo> <msub> <mi>b</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mn>1</mn> <mo>+</mo> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mo>(</mo> <mrow> <msub> <mi>&omega;</mi> <mi>i</mi> </msub> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <msub> <mi>x</mi> <mi>j</mi> </msub> <mo>+</mo> <msub> <mi>b</mi> <mi>i</mi> </msub> </mrow> <mo>)</mo> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mi>X</mi> <mi>I</mi> <mi>I</mi> <mi>I</mi> <mo>)</mo> </mrow> </mrow>

最小化：

使得：

其中：h(x_i)＝[g(α₁x_i+b₁),…,g(α_Kx_i+B_K)],ξ_i＝[ξ_i,1,…,ξ_i,m]^T是样本x_i对应输出的训练误差，依据KKT条件，训练ELM等价于解决如下的对偶优化问题：引入拉格朗日乘子α_i(i＝1,2,…n)：

<mrow> <msub> <mi>L</mi> <msub> <mi>D</mi> <mrow> <mi>E</mi> <mi>L</mi> <mi>M</mi> </mrow> </msub> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mo>|</mo> <mo>|</mo> <mi>&beta;</mi> <mo>|</mo> <msup> <mo>|</mo> <mn>2</mn> </msup> <mo>+</mo> <mi>C</mi> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <mo>|</mo> <mo>|</mo> <mi>&xi;</mi> <mo>|</mo> <msup> <mo>|</mo> <mn>2</mn> </msup> <mo>-</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <munderover> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>m</mi> </munderover> <msub> <mi>&alpha;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>h</mi> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> <msub> <mi>&beta;</mi> <mi>j</mi> </msub> <mo>-</mo> <msub> <mi>t</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>&xi;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mi>X</mi> <mi>V</mi> <mi>I</mi> <mo>)</mo> </mrow> </mrow>

其中，每个拉格朗日乘子α_i对应第i个训练样本；β_j连接隐含层节点到第j个输出节点的向量，β＝[β₁,…,β_m]，相应的KKT的优化条件如下：

<mrow> <mfrac> <mrow> <mo>&part;</mo> <msub> <mi>L</mi> <msub> <mi>D</mi> <mrow> <mi>E</mi> <mi>L</mi> <mi>M</mi> </mrow> </msub> </msub> </mrow> <mrow> <mo>&part;</mo> <msub> <mi>&beta;</mi> <mi>j</mi> </msub> </mrow> </mfrac> <mo>=</mo> <mn>0</mn> <mo>&DoubleRightArrow;</mo> <msub> <mi>&beta;</mi> <mi>j</mi> </msub> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>&alpha;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mi>k</mi> <msup> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mi>T</mi> </msup> <mo>&DoubleRightArrow;</mo> <mi>&beta;</mi> <mo>=</mo> <msup> <mi>H</mi> <mi>T</mi> </msup> <mi>&alpha;</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mi>X</mi> <mi>V</mi> <mi>I</mi> <mi>I</mi> <mo>)</mo> </mrow> </mrow>

<mrow> <mfrac> <mrow> <mo>&part;</mo> <msub> <mi>L</mi> <msub> <mi>D</mi> <mrow> <mi>E</mi> <mi>L</mi> <mi>M</mi> </mrow> </msub> </msub> </mrow> <mrow> <mo>&part;</mo> <msub> <mi>&xi;</mi> <mi>i</mi> </msub> </mrow> </mfrac> <mo>=</mo> <mn>0</mn> <mo>&DoubleRightArrow;</mo> <msub> <mi>&alpha;</mi> <mi>i</mi> </msub> <mo>=</mo> <msub> <mi>C&xi;</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>i</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mi>N</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mi>X</mi> <mi>V</mi> <mi>I</mi> <mi>I</mi> <mi>I</mi> <mo>)</mo> </mrow> </mrow>

<mrow> <mfrac> <mrow> <mo>&part;</mo> <msub> <mi>L</mi> <msub> <mi>D</mi> <mrow> <mi>E</mi> <mi>L</mi> <mi>M</mi> </mrow> </msub> </msub> </mrow> <mrow> <mo>&part;</mo> <msub> <mi>&alpha;</mi> <mi>i</mi> </msub> </mrow> </mfrac> <mo>=</mo> <mn>0</mn> <mo>&DoubleRightArrow;</mo> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mi>&beta;</mi> <mo>-</mo> <msubsup> <mi>t</mi> <mi>i</mi> <mi>T</mi> </msubsup> <mo>+</mo> <msubsup> <mi>&xi;</mi> <mi>i</mi> <mi>T</mi> </msubsup> <mo>=</mo> <mn>0</mn> <mo>,</mo> <mi>i</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mi>N</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mi>X</mi> <mi>I</mi> <mi>X</mi> <mo>)</mo> </mrow> </mrow> 1

将(XVII)、(XVIII)式带入(XIX)式，得到如下方程：

<mrow> <mo>(</mo> <mfrac> <mi>I</mi> <mi>C</mi> </mfrac> <mo>+</mo> <msup> <mi>HH</mi> <mi>T</mi> </msup> <mo>)</mo> <mi>&alpha;</mi> <mo>=</mo> <mi>T</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mo>(</mo> <mi>X</mi> <mi>X</mi> <mo>)</mo> </mrow>

其中：

将(XX)带入(XVII)式，得到隐元偏置：

其中：

<mrow> <mi>H</mi> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>N</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <msub> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mi>g</mi> <mrow> <mo>(</mo> <msub> <mi>&omega;</mi> <mn>1</mn> </msub> <mo>&CenterDot;</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>b</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mrow> <mi>g</mi> <mrow> <mo>(</mo> <msub> <mi>&omega;</mi> <mi>K</mi> </msub> <mo>&CenterDot;</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>b</mi> <mi>K</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>g</mi> <mrow> <mo>(</mo> <msub> <mi>&omega;</mi> <mn>1</mn> </msub> <mo>&CenterDot;</mo> <msub> <mi>x</mi> <mn>2</mn> </msub> <mo>+</mo> <msub> <mi>b</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mrow> <mi>g</mi> <mrow> <mo>(</mo> <msub> <mi>&omega;</mi> <mi>K</mi> </msub> <mo>&CenterDot;</mo> <msub> <mi>x</mi> <mn>2</mn> </msub> <mo>+</mo> <msub> <mi>b</mi> <mi>K</mi> </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> <mrow></mrow> </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>&omega;</mi> <mn>1</mn> </msub> <mo>&CenterDot;</mo> <msub> <mi>x</mi> <mi>N</mi> </msub> <mo>+</mo> <msub> <mi>b</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mrow> <mi>g</mi> <mrow> <mo>(</mo> <msub> <mi>&omega;</mi> <mi>K</mi> </msub> <mo>&CenterDot;</mo> <msub> <mi>x</mi> <mi>N</mi> </msub> <mo>+</mo> <msub> <mi>b</mi> <mi>K</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mrow> <mi>N</mi> <mo>&times;</mo> <mi>K</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mi>X</mi> <mi>X</mi> <mi>I</mi> <mi>I</mi> <mi>I</mi> <mo>)</mo> </mrow> </mrow>

通过式(XXII)计算得到的隐层权值β，将结果带入式(XII)，得到样本对应的输出。

2.根据权利要求1所述的建立方法，其特征在于，所述获得时间点对应的时间序列样本是指以一个变量的时间间隔为参照依据时间点对应对所有变量进行匹配，最终得到所需时间序列样本。

3.根据权利要求1所述的建立方法，其特征在于，所述确定模型输入变量的方法为通过相关性分析，以与铁水硅含量的变化相关性强的变量作为所述模型输入变量。

4.根据权利要求3所述的建立方法，其特征在于，结合时滞分析和/或专家经验确定模型输入变量。

5.根据权利要求1-4任一项所述的建立方法，其特征在于，所述S3数据分类具体包括以下步骤：

S31、以零为边界将所述输出样本划分为两大类，大于零和等于零的样本划分为正类样本，小于零的样本划分为负类样本；

S32、在所述正类样本中获取两个聚类中心c₁、c₂，在所述负类样本中获取两个聚类c₁′、c₂′，以这四个聚类中心得到两个边界值：

6.根据权利要求5所述的建立方法，其特征在于，所述S32在正类样本中获取两个聚类中心c₁、c₂，在负类样本中获取两个聚类c₁′、c₂′，具体包括以下步骤：

S321、初始化算法参数；具体方法如下：

初始化聚类个数c、模糊指数m、误差ε、算法迭代次数T，初始化聚类中心V＝[v₁,v₂,...,v_c]；在[0,1]范围内随机初始化隶属度矩阵U，使u_ij满足公式(Ⅳ)约束条件

<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mn>0</mn> <mo><</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msub> <mi>u</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo><</mo> <mi>n</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <munderover> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>c</mi> </munderover> <msub> <mi>u</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>=</mo> <mn>1</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>0</mn> <mo>&le;</mo> <msub> <mi>u</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>&le;</mo> <mn>1</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> <mn>1</mn> <mo>&le;</mo> <mi>i</mi> <mo>&le;</mo> <mi>n</mi> <mo>,</mo> <mn>1</mn> <mo>&le;</mo> <mi>j</mi> <mo>&le;</mo> <mi>c</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mi>I</mi> <mi>V</mi> <mo>)</mo> </mrow> </mrow>

其中c是聚类个数，n为样本数，x_i为聚类样本，u_ij为第i个样本属于第j个中心的隶属度，m为模糊指数，用来控制聚类的模糊程度；

S322、根据式(VII)计算K(x_i,v_j)，式(XIII)计算参数ρ；具体方法如下，构建模糊C均值聚类算法的目标函数如下：

<mrow> <msub> <mi>J</mi> <mi>&phi;</mi> </msub> <mrow> <mo>(</mo> <mi>U</mi> <mo>,</mo> <mi>V</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <munderover> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>c</mi> </munderover> <msubsup> <mi>u</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>m</mi> </msubsup> <mo>|</mo> <mo>|</mo> <mi>&phi;</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mi>&phi;</mi> <mrow> <mo>(</mo> <msub> <mi>v</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>|</mo> <msup> <mo>|</mo> <mn>2</mn> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mi>V</mi> <mo>)</mo> </mrow> </mrow>

式(Ⅴ)中，

||φ(x_i)-φ(v_j)||＝K(x_i,x_i)+K(v_j,v_j)-2K(x_i,v_j) (Ⅵ)

K(x_i,v_j)＝exp{-(||x_i-v_j||²)/ρ²} (VII)

<mrow> <mi>&rho;</mi> <mo>=</mo> <msqrt> <mrow> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <mfrac> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>-</mo> <mover> <mi>x</mi> <mo>&OverBar;</mo> </mover> <mo>|</mo> <msup> <mo>|</mo> <mn>2</mn> </msup> </mrow> <mi>n</mi> </mfrac> </mrow> </msqrt> <mo>,</mo> <mover> <mi>x</mi> <mo>&OverBar;</mo> </mover> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <mfrac> <msub> <mi>x</mi> <mi>i</mi> </msub> <mi>n</mi> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mi>V</mi> <mi>I</mi> <mi>I</mi> <mi>I</mi> <mo>)</mo> </mrow> </mrow>

<mrow> <msubsup> <mi>v</mi> <mi>j</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msubsup> <mo>=</mo> <mfrac> <mrow> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msup> <mrow> <mo>(</mo> <msubsup> <mi>u</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> </msubsup> <mo>)</mo> </mrow> <mi>m</mi> </msup> <mi>K</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>,</mo> <msubsup> <mi>v</mi> <mi>j</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> </msubsup> <mo>)</mo> </mrow> <msub> <mi>x</mi> <mi>i</mi> </msub> </mrow> <mrow> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msup> <mrow> <mo>(</mo> <msubsup> <mi>u</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> </msubsup> <mo>)</mo> </mrow> <mi>m</mi> </msup> <mi>K</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>,</mo> <msubsup> <mi>v</mi> <mi>j</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> </msubsup> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mi>I</mi> <mi>X</mi> <mo>)</mo> </mrow> </mrow>

<mrow> <msubsup> <mi>u</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mrow> <mo>(</mo> <mi>k</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msubsup> <mo>=</mo> <mfrac> <msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>K</mi> <mo>(</mo> <mrow> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>,</mo> <msubsup> <mi>v</mi> <mi>j</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msubsup> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> <mo>/</mo> <mrow> <mo>(</mo> <mi>m</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </msup> <mrow> <munderover> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>c</mi> </munderover> <msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>K</mi> <mo>(</mo> <mrow> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>,</mo> <msubsup> <mi>v</mi> <mi>j</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msubsup> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> <mo>/</mo> <mrow> <mo>(</mo> <mi>m</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </msup> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mi>X</mi> <mo>)</mo> </mrow> </mrow>

S324、依据算法目标函数进行算法迭代；即依据S323更新一次聚类中心和隶属度，带入式(V)得到一个新的目标函数值；

迭代终止的条件如下：

<mrow> <mo>|</mo> <mo>|</mo> <msubsup> <mi>J</mi> <mi>&phi;</mi> <mi>k</mi> </msubsup> <mo>-</mo> <msubsup> <mi>J</mi> <mi>&phi;</mi> <mrow> <mi>k</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>|</mo> <mo>|</mo> <mo>&le;</mo> <mi>&epsiv;</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mi>X</mi> <mi>I</mi> <mo>)</mo> </mrow> </mrow>

即更新前后目标函数值的差值小于给定误差或者达到给定的最大迭代次数时算法结束，否则继续更新；算法迭代结束即在所述正类样本中获得两个聚类中心c₁、c₂；在所述负类样本中获取两个聚类c₁′、c₂′。

7.根据权利要求1-4、6任一项所述的建立方法，其特征在于，所述依据步骤S3对所述输出样本的四类区间划分标准，将所述输出样本进行二进制编码即对所有样本(x_j,t_j)的输出t_j进行二进制编码；x_j指输入变量；t_j指输出变量；编码矩阵如下：

8.根据权利要求1所述的建立方法，其特征在于，还包括对所述预报模型的输出进行解码的步骤；解码过程如下：给定样本x，依据所述预报模型，得到m＝log₂c维的输出向量：

f(x)＝(f_m(x),…,f₁(x)) (XXIV)

样本(x,i)的总损失函数定义为：

<mrow> <msub> <mi>d</mi> <mi>H</mi> </msub> <mrow> <mo>(</mo> <mi>M</mi> <mo>(</mo> <mi>i</mi> <mo>)</mo> <mo>,</mo> <mi>y</mi> <mo>(</mo> <mi>x</mi> <mo>)</mo> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>q</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>m</mi> </munderover> <mrow> <mo>(</mo> <mfrac> <mrow> <mn>1</mn> <mo>-</mo> <mi>s</mi> <mi>i</mi> <mi>g</mi> <mi>n</mi> <mrow> <mo>(</mo> <mi>M</mi> <mo>(</mo> <mrow> <mi>i</mi> <mo>,</mo> <mi>q</mi> </mrow> <mo>)</mo> <msub> <mi>f</mi> <mi>q</mi> </msub> <mo>(</mo> <mi>x</mi> <mo>)</mo> <mo>)</mo> </mrow> </mrow> <mn>2</mn> </mfrac> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mi>X</mi> <mi>X</mi> <mi>V</mi> <mo>)</mo> </mrow> </mrow>

M(i)指编码矩阵c*m的第i个行向量，表示为：M(i)＝(M(i,m),…,M(i,1))；

基于上述损失函数，样本x最终的类输出为：

9.根据权利要求1-4、6、8任一项所述的建立方法，其特征在于，用状态转移算法对所述预报模型的输入权值和隐元偏差进行优化，得到最优的输入权值和隐元偏差；所述优化具体流程如下：

输入：随机产生ELM的输入权值ω_K*n和隐元偏差β_K*m；

输出：最优ELM输入权值和隐元偏差；

Step4：对STA算法中的样本执行：

Step41扩张转变

x_k+1＝x_k+γR_cx_k (XXVII)

其中，γ是一个正常数，称为扩张因子；R_c∈R^n*n是一个随机对角矩阵，其中每一个元素都服从高斯分布；

Step42旋转转变

<mrow> <msub> <mi>x</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <msub> <mi>x</mi> <mi>k</mi> </msub> <mo>+</mo> <mi>&alpha;</mi> <mfrac> <mn>1</mn> <mrow> <mi>n</mi> <mo>|</mo> <mo>|</mo> <msub> <mi>x</mi> <mi>k</mi> </msub> <mo>|</mo> <msub> <mo>|</mo> <mn>2</mn> </msub> </mrow> </mfrac> <msub> <mi>R</mi> <mi>r</mi> </msub> <msub> <mi>x</mi> <mi>k</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mi>X</mi> <mi>X</mi> <mi>V</mi> <mi>I</mi> <mi>I</mi> <mi>I</mi> <mo>)</mo> </mrow> </mrow>

其中，x_k∈Rⁿ，α代表正常数，称为旋转因子；R_r∈R^n*n是一个随机矩阵，服从[-1,1]的均匀分布，||·||₂代表二范数；

Step43压缩转变

x_k+1＝x_k+δR_ax_k (XXIX)

其中，δ是一个正常数，称为压缩因子；R_a∈R^n*n是一个随机对角矩阵，其中每一个元素都服从高斯分布，并且只有一个随机位置是非零值；

Step5：判断算法是否满足运行的终止条件，终止条件为最大的迭代次数以及算法的精度，若是满足，转step6，否则g＝g+1，转step4；

Step6：输出STA算法的最优值，即ELM的输入权值和隐元偏差。

10.一种高炉铁水硅含量四分类趋势预报方法，包括按权利要求1-9任一项所述方法建立高炉铁水硅含量四分类趋势预报模型，选取一组变量数据作为输入变量，输入到所述预报模型，得到高炉铁水硅含量四分类趋势预报结果；选取的输入到所述预报模型的输入变量与建立所述预报模型时采用的输入变量相同。