一种基于工业大数据挖掘的状态预测方法与流程

文档序号:12669584阅读:来源:国知局

技术特征:

1.一种基于工业大数据挖掘的状态预测方法,其特征在于,包括以下步骤:

步骤一,数据采集:将反映系统历史运行状态的样本作为训练集其中xi是系统状态变量,即模型的输入,ti是关注的预测指标,即模型的输出;

步骤二,OS-ELM模型:采用步骤一的训练样本建立若干个OS-ELM模型,并计算得到若干个预测值;

步骤三,EOS-ELM模型:对OS-ELM模型的的预测结果取平均值,得到EOS-ELM模型预测结果。

2.如权利要求1所述的一种基于工业大数据挖掘的状态预测方法,其特征在于,所述OS-ELM模型建立过程包括:

初始化之前,先确定网络初始参数:网络有L个隐含节点,首先确定隐含节点类型,隐含节点类型包括为RBF或additive隐含节点;

初始化阶段,从训练样本中选取部分样本N0≥L进行初始化,包括以下步骤:

步骤1:随机对输入参数赋值

其中,对RBF隐含节点,参数为中心点ai和影响因子bi;对additive隐含节点,参数为输入权重ai和偏差bi

步骤2:计算初始的隐含层输出矩阵H0

<mrow> <msub> <mi>H</mi> <mn>0</mn> </msub> <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> <msub> <mi>N</mi> <mn>0</mn> </msub> </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> <msub> <mi>N</mi> <mn>0</mn> </msub> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mrow> <msub> <mi>N</mi> <mn>0</mn> </msub> <mo>&times;</mo> <mi>L</mi> </mrow> </msub> </mrow>

其中,利用RBF隐含节点时,G(ai,bi,xj)=g(bi||xj-ai||),bi∈R+,当利用additive隐含节点时,G(ai,bi,xj)=g(ai·xj+bi),bi∈R;

步骤3:估计初始的输出权重β(0)

权重求解问题转化为最小化||H0β-t0||;

由ELM算法的求解可知,可知最小化||H0β-t0||的求解结果为其中

步骤4:令k=0,k表示加入网络的数据块的数量;

连续学习阶段,包括以下步骤:

k+1数据块的观测值为:

其中,Nk+1是k+1数据块的观测值的数量;

步骤5:计算局部的隐含层输出矩阵Hk+1

<mrow> <msub> <mi>H</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <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> <mrow> <mi>i</mi> <mo>=</mo> <mrow> <mo>(</mo> <msubsup> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mi>k</mi> </msubsup> <msub> <mi>N</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mn>1</mn> </mrow> </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> <mrow> <mi>i</mi> <mo>=</mo> <mrow> <mo>(</mo> <msubsup> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mi>k</mi> </msubsup> <msub> <mi>N</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mn>1</mn> </mrow> </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> <mrow> <mi>i</mi> <mo>=</mo> <mrow> <mo>(</mo> <msubsup> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msubsup> <msub> <mi>N</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mn>1</mn> </mrow> </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> <mrow> <mi>i</mi> <mo>=</mo> <mrow> <mo>(</mo> <msubsup> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msubsup> <msub> <mi>N</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mrow> <msub> <mi>N</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>&times;</mo> <mi>L</mi> </mrow> </msub> </mrow>

步骤6:设置参数

<mrow> <msub> <mi>T</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <msubsup> <mrow> <mo>&lsqb;</mo> <mrow> <msub> <mi>t</mi> <mrow> <mi>i</mi> <mo>=</mo> <mrow> <mo>(</mo> <msubsup> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mi>k</mi> </msubsup> <msub> <mi>N</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>,</mo> <mn>...</mn> <mo>,</mo> <msub> <mi>t</mi> <mrow> <msubsup> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msubsup> <msub> <mi>N</mi> <mi>j</mi> </msub> </mrow> </msub> </mrow> <mo>&rsqb;</mo> </mrow> <mrow> <msub> <mi>N</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>&times;</mo> <mi>m</mi> </mrow> <mi>T</mi> </msubsup> </mrow>

步骤7:计算输出权重β(k+1)

<mrow> <msub> <mi>P</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <msub> <mi>P</mi> <mi>k</mi> </msub> <mo>-</mo> <msub> <mi>P</mi> <mi>k</mi> </msub> <msubsup> <mi>H</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> <mi>T</mi> </msubsup> <msup> <mrow> <mo>(</mo> <mi>I</mi> <mo>+</mo> <msub> <mi>H</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <msub> <mi>P</mi> <mi>k</mi> </msub> <msubsup> <mi>H</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> <mi>T</mi> </msubsup> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msub> <mi>H</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <msub> <mi>P</mi> <mi>k</mi> </msub> </mrow>

<mrow> <msup> <mi>&beta;</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msup> <mo>=</mo> <msup> <mi>&beta;</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> </msup> <mo>+</mo> <msub> <mi>P</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <msubsup> <mi>H</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> <mi>T</mi> </msubsup> <mrow> <mo>(</mo> <msub> <mi>T</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>H</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <msup> <mi>&beta;</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> </msup> <mo>)</mo> </mrow> </mrow>

步骤8:令k=k+1,返回步骤5;

当所有训练数据都参与训练时,循环结束,计算预测输出,即OS-ELM预测值;

步骤1-8重复J次,J为OS-ELM模型数量。

3.如权利要求2所述的一种基于工业大数据挖掘的状态预测方法,其特征在于,步骤二中,在所述连续学习阶段后还包括模型评估阶段,该阶段具体为:

利用连续学习阶段产生的参数a、b、权重β以及测试的输入数据,得到预测输出,即EOS-ELM预测值。

4.如权利要求2所述的一种基于工业大数据挖掘的状态预测方法,其特征在于,对所述OS-ELM模型的的预测结果取平均值具体为:对所述变量xi,每个OS-ELM模型的输出为fj(xi),j=1,...,J,则EOS-ELM的预测输出为

5.如权利要求1所述的一种基于工业大数据挖掘的状态预测方法,其特征在于:在所述EOS-ELM算法前,使用变量选择与特征提取方法对步骤二得到的若干个预测值确定与输出相关的变量,以降低数据维数。

当前第2页1 2 3 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1