发明名称 一种动态模块化神经网络的氨氮在线软测量方法
摘要 一种动态模块化神经网络的氨氮在线软测量方法,属于检测技术及仪表研究技术领域。本发明针对当前污水处理过程中出水水质关键参数氨氮NH3-N测量周期长,不能在线检测的问题,采用了多模型的软测量方法实现了辅助变量与NH3-N之间的映射,对NH3-N进行在线软测量。与一般的在线学方法相比较具有更高的预测精度,是一种使用范围广,可靠性高,鲁棒性强的软测量方法;本发明克服了单模型建模学时间长、过程匹配不佳、外推能力差、没有知识积累等缺陷。
申请公布号 CN102662040B 申请公布日期 2014.06.18
申请号 CN201210119203.3 申请日期 2012.04.20
申请人 辽宁工程技术大学 发明人 张昭昭;郭伟;徐光宪
分类号 G01N33/18(2006.01)I;G06N3/08(2006.01)I 主分类号 G01N33/18(2006.01)I
代理机构 沈阳东大知识产权代理有限公司 21109 代理人 李运萍
主权项 1.一种动态模块化神经网络的氨氮在线软测量方法,其特征在于:包括以下步骤:步骤1:设计用于氨氮在线软测量的动态模块化神经网络拓扑结构;网络分为5层:输入层、RBF层、任务分配层、子网络层、集成输出层;输入为曝气池进水水质指标,输出为出水水质氨氮浓度;输入层:包括n个输入节点,n为输入样本的维数,k时刻的输入样本表示为x(k)=[x<sub>1</sub>(k),…,x<sub>n</sub>(k)]<sup>T</sup>;RBF层:包括l个RBF神经元,每个RBF神经元的激活函数为<img file="FDA0000492466110000011.GIF" wi="534" he="96" />在动态模块化神经网络中,激活函数Φ<sub>i</sub>(x(k))被定义为硬限幅函数:若满足公式||x(k)-c<sub>i</sub>||≤δ<sub>i</sub>,则Φ<sub>i</sub>(x(k))=1    (1)若满足公式||x(k)-c<sub>i</sub>||>δ<sub>i</sub>,则Φ<sub>i</sub>(x(k))=0式中,c<sub>i</sub>为第i个RBF神经元的数据中心,δ<sub>i</sub>为第i个RBF神经元的扩展宽度;公式(1)表示在输入样本空间中以c<sub>i</sub>为球心,以δ<sub>i</sub>为半径作一个超球面,当输入点在该超球面内时输出为1,在超球面以外时输出为0;任务分配层:该层的任务是根据RBF层中RBF神经元的输出和k时刻输入样本x(k)在输入样本空间中的位置,由模糊策略决定给k时刻输入样本x(k)分配不同的子网络进行学习;子网络层:该层由子神经网络组成,其个数与RBF层中RBF神经元个数相同,且为一一对应的关系;该层的主要任务是学习由任务分配层分配来的输入样本;集成输出层:该层的主要任务是对k时刻各子神经网络的学习结果进行集成输出,集成输出的结果Y(x(k)),公式如下:<maths num="0001"><![CDATA[<math><mrow><mi>Y</mi><mrow><mo>(</mo><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>p</mi></munderover><msub><mi>&alpha;</mi><mi>i</mi></msub><msub><mi>y</mi><mi>ik</mi></msub><mrow><mo>(</mo><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,p为被x(k)激活的子神经网络的个数;α<sub>i</sub>为被激活的第i个子神经网络的集成权值,且有<img file="FDA0000492466110000013.GIF" wi="202" he="136" />y<sub>ik</sub>为子网络层中第i个子神经网络对k时刻的输入样本x(k)的输出,且有y<sub>ik</sub>(x(k))=v<sub>i</sub>f(w<sub>i</sub>x(k)+b1<sub>i</sub>)+b2<sub>i</sub>式中,w<sub>i</sub>为第i个子神经网络输入层到隐含层的权连接矩阵;v<sub>i</sub>为第i个子神经网络隐含层到输出层的权连接矩阵;f(·)为子神经网络隐含层神经元的激活函数;b1<sub>i</sub>为第i个子神经网络隐含层神经元的偏移向量;b2<sub>i</sub>为第i个子神经网络输出层神经元的偏移向量;步骤2:输入变量的确定及归一化处理;确定输入变量为:污水出水水质氨氮预测模型输入样本包括入水氨氮、入水化学需氧量、入水悬浮物、入水PH值、反应池溶解氧浓度、反应池混合液可挥发性悬浮固体浓度,输出样本为出水氨氮预测值;由于上述自变量之间差异较大且量纲不同,采用在线递归中心化压缩方法进行归一化处理:设k时刻的输入样本和输出样本对为:Z(k)=[x<sub>1</sub>(k),x<sub>2</sub>(k),…,x<sub>6</sub>(k),x<sub>7</sub>(k)]    (3)式中,x<sub>1</sub>(k),x<sub>2</sub>(k),…,x<sub>6</sub>(k)表示输入样本在k时刻的值,作为k时刻神经网络的输入;x<sub>7</sub>(k)表示k时刻氨氮浓度的目标值,作为k时刻神经网络的期望输出值;则在线递归中心化压缩,公式为:<maths num="0002"><![CDATA[<math><mrow><msub><mi>Z</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>Z</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msub><mover><mi>Z</mi><mo>&OverBar;</mo></mover><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>&sigma;</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><msub><mover><mi>Z</mi><mo>&OverBar;</mo></mover><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>k</mi></mfrac><msub><mover><mi>Z</mi><mo>&OverBar;</mo></mover><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><mfrac><mn>1</mn><mi>k</mi></mfrac><msub><mi>Z</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msubsup><mi>&sigma;</mi><mi>j</mi><mn>2</mn></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>k</mi></mfrac><msubsup><mi>&sigma;</mi><mi>j</mi><mn>2</mn></msubsup><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><mfrac><mn>1</mn><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></mfrac><msup><mrow><mo>(</mo><msub><mi>Z</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msub><mover><mi>Z</mi><mo>&OverBar;</mo></mover><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>,</mo><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mn>7</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,Z<sub>j</sub>(k)表示k时刻输入样本的第j个分量,<img file="FDA0000492466110000024.GIF" wi="128" he="84" />表示j个分量的平均值,σ<sub>j</sub>表示j个分量之间的标准差;步骤3:对归一化后的k时刻输入样本x(k),首先,在线更新RBF层RBF神经元的个数及各RBF神经元的数据中心;其次,任务分配层根据RBF层RBF神经元的输出,由模糊策略确定并激活子网络层中不同的子神经网络学习该输入样本x(k);最后对被激活的各子神经网络的输出进行动态集成形成输出;步骤4:集成输出Y(k)即为k时刻出水氨氮的软测量结果,采集下一时刻的输入数据,返回步骤2。
地址 123000 辽宁省阜新市中华路47号