发明名称 焦化炉压力的多目标遗传算法与RBF神经网络优化建模方法
摘要 本发明公开了一种焦化炉压力的多目标遗传算法与RBF神经网络优化建模方法。本发明通过采集过程对象的输入输出数据,结合RBF神经网络模型,使用改进的MOEA优化神经网络的网络层和参数。本发明具有较高的精确性,能很好地描述过程对象的动态特性。
申请公布号 CN105608295A 申请公布日期 2016.05.25
申请号 CN201610063623.2 申请日期 2016.01.29
申请人 杭州电子科技大学 发明人 张日东;王玉中
分类号 G06F17/50(2006.01)I;G06N3/02(2006.01)I;G06N3/12(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 浙江杭州金通专利事务所有限公司 33100 代理人 王佳健
主权项 焦化炉压力的多目标遗传算法与RBF神经网络优化建模方法,其特征在于该方法的具体步骤是:步骤1、采集过程的实时运行数据,建立过程对象RBF模型,具体步骤如下:1.1由包含输入层、输出层和隐含层的RBF神经网络结构,得到网络的映射关系即系统的输入输出模型,形式如下:<img file="FDA0000917505010000011.GIF" wi="581" he="135" />其中,x=(x<sub>1</sub>,x<sub>2</sub>,…,x<sub>n</sub>)表示n输入结点向量,y表示网络的输出变量,c<sub>i</sub>∈R<sup>n</sup>表示第i个隐含层神经元的中心向量,R<sup>n</sup>是欧氏空间,<img file="FDA0000917505010000012.GIF" wi="526" he="164" />是一个高斯函数,||x‑c<sub>i</sub>||表示x到c<sub>i</sub>的径向距离,σ<sub>i</sub>是高斯函数的基宽,1≤i≤n<sub>r</sub>,n<sub>r</sub>是隐含层的结点数,ω<sub>i</sub>表示第i个隐含层和输出层之间的连接权;步骤2、利用MOEA遗传算法优化RBF神经网络模型的参数,具体步骤是:2.1首先对神经网络模型参数进行编码,得到如下形式的第l代染色体:<img file="FDA0000917505010000013.GIF" wi="880" he="463" />其中,l=1,2,...,N,N是种群规模大小,m,n和n<sub>r</sub>是正整数,1≤m≤5,1≤n≤5,1≤n<sub>r</sub>≤60,C<sub>l</sub>中的元素满足如下条件:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>c</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>y</mi><mi>min</mi></msub><mo>+</mo><mi>r</mi><mrow><mo>(</mo><msub><mi>y</mi><mi>max</mi></msub><mo>-</mo><msub><mi>y</mi><mi>min</mi></msub><mo>)</mo></mrow></mrow></mtd><mtd><mrow><mn>1</mn><mo>&le;</mo><mi>i</mi><mo>&le;</mo><msub><mi>n</mi><mi>r</mi></msub><mn>1</mn><mo>&le;</mo><mi>j</mi><mo>&le;</mo><mi>n</mi></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>u</mi><mi>min</mi></msub><mo>+</mo><mi>r</mi><mrow><mo>(</mo><msub><mi>u</mi><mi>max</mi></msub><mo>-</mo><msub><mi>u</mi><mi>min</mi></msub><mo>)</mo></mrow></mrow></mtd><mtd><mrow><mn>1</mn><mo>&le;</mo><mi>i</mi><mo>&le;</mo><msub><mi>n</mi><mi>r</mi></msub><mn>6</mn><mo>&le;</mo><mi>j</mi><mo>&le;</mo><mn>5</mn><mo>+</mo><mi>m</mi></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000917505010000014.GIF" wi="1150" he="199" /></maths>σ<sub>i</sub>=rw<sub>max</sub> 1≤i≤n<sub>r</sub>其中,r是一个在‑0.5~1.5之间的随机数,u<sub>min</sub>,u<sub>max</sub>是系统输入的最小值和最大值,y<sub>min</sub>,y<sub>max</sub>是系统输出的最小值和最大值,w<sub>max</sub>是高斯基函数的最大宽度;2.2将采集的过程对象的数据样本分为三部分,前面1/3数据样本为训练数据样本Y<sub>1</sub>,用于计算输出层权重,中间1/3数据为第二组数据样本Y<sub>2</sub>,用于对每一代的神经网络评价,后面1/3数据为优化数据Y<sub>3</sub>,用于求取Pareto最优解,选取RBF神经网络的目标函数,形式如下:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>M</mi><mi>i</mi><mi>n</mi><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>f</mi><mn>1</mn></msub><mo>=</mo><msqrt><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msup><mrow><mo>|</mo><msub><mi>Y</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>-</mo><msub><mover><mi>Y</mi><mo>^</mo></mover><mn>1</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup></mrow></msqrt><mo>+</mo><msqrt><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mn>2</mn></msub></munderover><msup><mrow><mo>|</mo><msub><mi>Y</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>-</mo><msub><mover><mi>Y</mi><mo>^</mo></mover><mn>2</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup></mrow></msqrt></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>f</mi><mn>2</mn></msub><mo>=</mo><mrow><mo>(</mo><mi>m</mi><mo>+</mo><mi>n</mi><mo>)</mo></mrow><msub><mi>n</mi><mi>r</mi></msub></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000917505010000021.GIF" wi="1053" he="252" /></maths>其中,Min表示求最小值,f<sub>1</sub>表示Y<sub>1</sub>和Y<sub>2</sub>均方差,f<sub>2</sub>反映了输入层和隐含层结构的复杂性,<img file="FDA0000917505010000022.GIF" wi="300" he="87" />分别表示两组数据样本Y<sub>1</sub>和Y<sub>2</sub>的RBF神经网络模型预测输出,N<sub>1</sub>,N<sub>2</sub>为选取的数据样本的大小;2.3由于染色体中的基因会发生变异,p<sub>c</sub>为当前个体C<sub>l</sub>和下一个个体C<sub>l+1</sub>之间交叉概率,以交叉概率p<sub>c</sub>将被选择的染色体进行交叉操作产生下一代子染色体C′<sub>l</sub>和C′<sub>l+1</sub>;变异算子执行时,随机产生给定范围内的m,n和n<sub>r</sub>值,变异个体根据步骤2.1元素个体满足条件进行变异得到新的个体;2.4为了提高MOEA的局部搜索能力,局部算子设计为如下形式:C=αC<sub>l</sub>+(1‑α)C<sub>l'</sub>C=C<sub>l</sub>+ΔC<sub>l</sub>其中,C<sub>l</sub>是选自前λf<sub>1</sub>个体,C<sub>l'</sub>是选自前λf<sub>2</sub>个体,λ表示自然数,α∈(0,1)为随机数,当C<sub>l</sub>=C<sub>l'</sub>则ΔC<sub>l</sub>中的Δc<sub>ij</sub>=αc<sub>ij</sub>,α∈(‑1,1),局部搜索概率动态改变成如下形式:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>p</mi><mi>l</mi></msub><mo>=</mo><mn>0.02</mn><mo>+</mo><mfrac><mn>0.2</mn><mrow><mn>1</mn><mo>+</mo><mi>exp</mi><mo>&lsqb;</mo><mo>-</mo><mn>0.01</mn><mrow><mo>(</mo><mi>g</mi><mo>-</mo><mi>G</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mo>&rsqb;</mo></mrow></mfrac></mrow>]]></math><img file="FDA0000917505010000023.GIF" wi="814" he="142" /></maths>其中,G表示最大的进化代数,g表示进化的代数;2.5在个体数目大于种群规模N时,得到剪接算子,形式如下<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>Af</mi><mi>i</mi></msub><mo>=</mo><msup><mi>&rho;e</mi><mrow><mo>-</mo><mo>|</mo><mo>|</mo><mi>x</mi><mo>-</mo><msub><mi>c</mi><mi>i</mi></msub><mo>|</mo><mo>|</mo></mrow></msup><mfrac><mrow><msub><mi>&phi;</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow></mrow><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>n</mi><mi>r</mi></msub></munderover><msub><mi>&phi;</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow></mrow></mfrac><mo>,</mo><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><msub><mi>n</mi><mi>r</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000917505010000024.GIF" wi="893" he="230" /></maths>其中,Af<sub>i</sub>表式第i个隐含神经元活跃度,φ<sub>i</sub>(x)表示第i个隐含神经元的输出值,ρ&gt;1;2.6设定优良基因库的最大值为N,将优良的基因保存到基因库中,当优良基因库大于N时,快速非支配排序方法首先执行,然后将支配方法从优良基因库中移除,使所有基因满足Pareto最优解以维持优良基因的多样性和均匀性;2.7依照步骤2.2到步骤2.6中的步骤进行循环重复优化搜索,达到允许的最大进化代数式结束优化搜索计算,得到改进后的MOEA遗传算法优化后的染色体,经解码后得到优化后的RBF神经网络模型的参数。
地址 310018 浙江省杭州市下沙高教园区2号大街