发明名称 基于遗传算法优化BP神经网络的丙烯聚合生产过程最优软测量仪表及方法
摘要 一种基于遗传算法优化BP神经网络的丙烯聚合生产过程最优软测量仪表,包括丙烯聚合生产过程、现场智能仪表、控制站、存放数据的DCS数据库、基于遗传算法优化BP神经网络的最优软测量模型以及熔融指数软测量值显示仪,现场智能仪表及控制站与丙烯聚合生产过程相连,与DCS数据库相连;最优软测量模型与DCS数据库及软测量值显示仪相连。所述的基于遗传算法优化BP神经网络的最优软测量模型包括数据预处理模块、ICA独立成分分析模块、BP神经网络建模模块以及遗传算法优化BP神经网络模块。以及提供了一种用软测量仪表实现的软测量方法。本发明实现在线测量、在线参数自动优化、计算速度快、模型自动更新、抗干扰能力强、精度高。
申请公布号 CN101315557A 申请公布日期 2008.12.03
申请号 CN200810063607.9 申请日期 2008.06.25
申请人 浙江大学 发明人 刘兴高;楼巍
分类号 G05B19/418(2006.01);C08F10/06(2006.01);C08F2/00(2006.01);G06N3/12(2006.01);G06N3/02(2006.01) 主分类号 G05B19/418(2006.01)
代理机构 杭州天正专利事务所有限公司 代理人 王兵;王利强
主权项 1、一种基于遗传算法优化BP神经网络的丙烯聚合生产过程最优软测量仪表,包括丙烯聚合生产过程、用于测量易测变量的现场智能仪表、用于测量操作变量的控制站、存放数据的DCS数据库以及熔融指数软测量值显示仪,所述现场智能仪表、控制站与丙烯聚合生产过程连接,所述现场智能仪表、控制站与DCS数据库连接,其特征在于:所述软测量仪表还包括基于遗传算法优化BP神经网络的最优软测量模型,所述DCS数据库与所述基于遗传算法优化BP神经网络的最优软测量模型的输入端连接,所述基于遗传算法优化BP神经网络的最优软测量模型的输出端与熔融指数软测量值显示仪连接,所述基于遗传算法优化BP神经网络的最优软测量模型包括:数据预处理模块,用于将从DCS数据库输入的模型输入变量进行预处理,对输入变量中心化,即减去变量的平均值;再对输入变量预白化处理即变量去相关,对输入变量施加一个线性变换;通过独立成分分析方法,从经过中心化和预白化处理的线性混合数据中恢复出基本源信号;BP神经网络模型模块,用于采用BP神经网络,设BP神经网络的输出层第k个神经元的实际输出为yk,输入为netk,与此层相邻的隐含层中任一神经元j的输出为yj,则有:<math><mrow><msub><mi>net</mi><mi>k</mi></msub><mo>=</mo><munder><mi>&Sigma;</mi><mi>i</mi></munder><msub><mi>w</mi><mi>kj</mi></msub><mo>&CenterDot;</mo><msub><mi>y</mi><mi>j</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>yk=f(netk) (2)式中,wkj是神经元k与神经元j之间的连接权,f(·)为神经元的输出函数,通常取为Sigmoid函数,表示为:<math><mrow><mi>f</mi><mrow><mo>(</mo><msub><mi>net</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>=</mo><mn>1</mn><mo>/</mo><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msup><mi>e</mi><mrow><mo>-</mo><mrow><mo>(</mo><msub><mi>net</mi><mi>k</mi></msub><mo>+</mo><msub><mi>h</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>/</mo><msub><mi>&theta;</mi><mn>0</mn></msub></mrow></msup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>式中,hk为神经元k的阈值,θ0为陡度参数,用以调节Sigmoid函数的陡度;令训练样本为k,对于任一个输入模式Xp,若相应有输出层中第k个神经元的期望输出Opk,则输出层的输出方差表示为:<math><mrow><msub><mi>E</mi><mi>p</mi></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><munder><mi>&Sigma;</mi><mi>p</mi></munder><msup><mrow><mo>(</mo><msub><mi>O</mi><mi>pk</mi></msub><mo>-</mo><msub><mi>y</mi><mi>pk</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>式中,Opk代表期望输出,ypk代表实际输出;反向误差传播学习的目的是修改连接权w值,使E达到最小值;要求连接权wkj,wji应沿Ep的负梯度方向学习;所以wkj的修正量为:<math><mrow><msub><mi>&Delta;</mi><mi>p</mi></msub><msub><mi>w</mi><mi>kj</mi></msub><mo>=</mo><mo>-</mo><mi>&beta;</mi><mfrac><mrow><mo>&PartialD;</mo><msub><mi>E</mi><mi>p</mi></msub></mrow><msub><mrow><mo>&PartialD;</mo><mi>w</mi></mrow><mi>kj</mi></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>式中,β为学习速率调整因子;遗传算法优化BP神经网络模块,用于采用遗传算法优化方法对BP神经网络进行优化,具体步骤如下:①确定网络权值的编码方式,个体的位串长度,染色体采用二进制编码,网络每个链接权值都用一定长度的0/1串表示,在编码时首先假定连接权在某一预定的范围内变化,连接权的实际值和字符串表示值关系表示为:<math><mrow><msub><mi>W</mi><mi>t</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>W</mi><mi>min</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>+</mo><mfrac><mrow><mi>binrep</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mrow><msup><mn>2</mn><mi>t</mi></msup><mo>-</mo><mn>1</mn></mrow></mfrac><mo>[</mo><msub><mi>W</mi><mi>max</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>W</mi><mi>min</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>式中,binrep(t)是二进制整数,[Wmax(i,j),Wmin(i,j)]为各连接权的变化范围,将所有权值对应的二进制字符串串联在一起,得到一条染色体,代表网络的一种权值组合。②设定的种群规模n,随机产生初始种群;③设定网络的目标函数,将其转换为适应度,对个网络权值进行评价;通过网络的误差函数计算适应度函数,并认为误差大的个体适应度小,个体i的适应度函数fi表示为:fi=1/(Ei+1) (7)式中,Ei是网络的误差函数,表示为:<math><mrow><msub><mi>E</mi><mi>i</mi></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>n</mi><mn>3</mn></msub></munderover><msup><mrow><mo>(</mo><msubsup><mi>c</mi><mi>j</mi><mi>k</mi></msubsup><mo>-</mo><msub><mi>c</mi><mi>j</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>式中,m是训练集的样本个数,n3是神经元数,cj k是单元j的理想输出cj是单元j的实际输出;④选定遗传操作,设置遗传参数及自适应调整算法,交叉概率设为Pc,变异概率设为Pm,个体i被选择的概率P表示为:<math><mrow><mi>P</mi><mo>=</mo><msub><mi>f</mi><mi>i</mi></msub><mo>/</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>f</mi><mi>j</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>式中,n为群体大小,fi为个体i的适应度,fj为个体j的适应度;⑤依据适应度在遗传空间进行选择操作;⑥依据选定的交叉、变异及有关算法、参数,进行相应的操作,得到新一代种群;⑦判断是否满足性能要求,若是,结束寻优,得到一组优化的权值;否则返回步骤⑤,继续迭代寻优。
地址 310027浙江省杭州市西湖区浙大路38号