发明名称 一种基于动态前馈神经网络的生化需氧量在线软测量方法
摘要 本发明一种基于动态前馈神经网络的生化需氧量在线软测量方法,通过设计用于污水曝气池生化需氧量软测量的动态前馈神经网络拓扑结构,确定动态前馈神经网络的输入样本,并对输入样本进行在线归一化处理;采用标准化后的数据训练神经网络,计算神经网络中与隐节点相连的所有权连接值在每次训练过程中的变化情况,判断该隐节点的活跃性,对活跃度较大的隐节点进行分裂;通过计算隐节点在训练过程中绝对输出的变化情况,判断该隐节点学信息的能力,删除没有学能力的隐节点;对神经网络参数进行调整;动态前馈神经网络训练过程结束后,确定曝气池出水生化需氧量BOD;本发明具有实时性好、稳定性好、精度高、神经网络的泛化能力强的优点。
申请公布号 CN102879541A 申请公布日期 2013.01.16
申请号 CN201210271656.8 申请日期 2012.07.31
申请人 辽宁工程技术大学 发明人 张昭昭;郭伟;张美金
分类号 G01N33/18(2006.01)I;G06N3/02(2006.01)I 主分类号 G01N33/18(2006.01)I
代理机构 沈阳东大专利代理有限公司 21109 代理人 李运萍
主权项 1.一种基于动态前馈神经网络的生化需氧量在线软测量方法,其特征在于:包括以下步骤:步骤1、设计用于污水曝气池生化需氧量软测量的动态前馈神经网络拓扑结构;其中,网络分为三层:输入层、隐含层和输出层;所述的输入层为污水曝气池进水水质指标,所述的输出层为出水生化需氧量BOD;步骤1-1、初始化神经网络结构:确定神经网络结构为l-1-1的连接方式,即输入层神经元为l个,在网络开始训练的时刻,隐含层的神经元为1个,输出层神经元1个;步骤1-2、对神经网络的权值进行随机赋值;神经网络的输入表示为x=(x<sub>1</sub>,x<sub>2</sub>,…,x<sub>l</sub>)<sup>T</sup>,其中(x<sub>1</sub>,x<sub>2</sub>,…,x<sub>l</sub>)<sup>T</sup>为(x<sub>1</sub>,x<sub>2</sub>,…,x<sub>l</sub>)的转置,神经网络的期望输出表示为y<sub>d</sub>;k时刻的训练样本为x(k)=(x<sub>1</sub>(k),x<sub>2</sub>(k),…,x<sub>l</sub>(k))<sup>T</sup>;假设k时刻动态动态前馈神经网络的隐含层中有m个隐节点,则用k时刻的训练样本训练神经网络时,动态前馈神经网络的输出描述为:<maths num="0001"><![CDATA[<math><mrow><mi>y</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>w</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mi>f</mi><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>l</mi></munderover><msub><mi>v</mi><mi>ji</mi></msub><msub><mi>x</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,m表示k时刻当前动态前馈神经网络中隐含层中的神经元数;x<sub>i</sub>(k)表示k时刻第i个输入节点的输入,其等于k时刻中x(k)=(x<sub>1</sub>(k),x<sub>2</sub>(k),…,x<sub>l</sub>(k))<sup>T</sup>的第i个分量;v<sub>ji</sub>表示k时刻第j个隐节点与第i个输入节点之间的权连接值;w<sub>j</sub>(k)表示k时刻第j个隐节点与输出节点之间的权连接值;<img file="FDA00001964372900012.GIF" wi="273" he="120" />表示k时刻动态前馈神经网络中隐含层中第j个隐节点的输出,且有<maths num="0002"><![CDATA[<math><mrow><msub><mi>h</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mi>f</mi><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>l</mi></munderover><msub><mi>v</mi><mi>ji</mi></msub><msub><mi>x</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>所述的隐节点输出f(·)的函数表达式为:<maths num="0003"><![CDATA[<math><mrow><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>x</mi></mrow></msup></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>步骤1-3、定义误差函数为:<maths num="0004"><![CDATA[<math><mrow><mi>e</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mi>M</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mrow><mo>(</mo><mi>y</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>y</mi><mi>d</mi></msub><msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mi>T</mi></msup><mi>y</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>y</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,M为训练样本总数;y(k)和y<sub>d</sub>(k)分别表示k时刻神经网络的实际输出和期望输出;步骤2、确定动态前馈神经网络的输入样本,并对输入样本进行在线归一化处理;步骤2-1、确定动态前馈神经网络的输入样本,包括曝气池进水溶解氧浓度DO、化学需氧量COD、悬浮物SS浓度和曝气池PH值;输出为曝气池出水生化需氧量BOD;步骤2-2、对动态前馈神经网络的输入样本进行在线归一化处理;设k时刻的输入样本和输出样本对为:z(k)=(x<sub>1</sub>(k),x<sub>2</sub>(k),x<sub>3</sub>(k),x<sub>4</sub>(k),x<sub>5</sub>(k))(4)其中,x<sub>1</sub>(k),x<sub>2</sub>(k),x<sub>3</sub>(k),x<sub>4</sub>(k)表示输入样本在k时刻的值,作为k时刻神经网络的输入;x<sub>5</sub>(k)表示k时刻表示曝气池出水生化需氧量的目标值,作为k时刻神经网络的期望输出;在线递归法对输入输出数据进行标准化方法如下:<maths num="0005"><![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><msubsup><mi>&delta;</mi><mi>j</mi><mn>2</mn></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></mfrac></mrow></math>]]></maths>j=[1,5],k=2,3,4,…                            (5)其中,z<sub>j</sub>(k)表示k时刻z(k)中的第j个分量;<img file="FDA00001964372900022.GIF" wi="106" he="54" />表示k时刻第j个分量的平均值;δ<sub>j</sub>(k)表示k时刻第j个分量之间的标准差;Z<sub>j</sub>(k)表示标准化之后的k时刻第j个分量;通过在线递归计算<img file="FDA00001964372900023.GIF" wi="105" he="55" />和<img file="FDA00001964372900024.GIF" wi="136" he="63" /><maths num="0006"><![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><mn>1</mn><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths><maths num="0007"><![CDATA[<math><mrow><msub><mover><mi>z</mi><mo>&OverBar;</mo></mover><mi>j</mi></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0008"><![CDATA[<math><mrow><msubsup><mi>&delta;</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>&delta;</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><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><mn>1</mn><mo>)</mo></mrow><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>δ<sub>j</sub>(1)=0(7)步骤3、采用标准化后的数据训练神经网络,计算神经网络中与隐节点相连的所有权连接值在每次训练过程中的变化情况,判断该隐节点的活跃性,对活跃度较大的隐节点进行分裂;通过计算隐节点在训练过程中绝对输出的变化情况,判断该隐节点学习信息的能力,删除没有学习能力的隐节点;对神经网络参数进行调整;步骤4、动态前馈神经网络训练过程结束后,确定曝气池出水生化需氧量BOD;将测试样本数据按照公式(5)、(6)(7)进行标准化处理,标准化处理后的测试样本数据作为训练好的神经网络的输入样本,神经网络的输出按照公式(8)进行反标准化处理:<maths num="0009"><![CDATA[<math><mrow><msub><mi>z</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>Z</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msubsup><mi>&delta;</mi><mi>j</mi><mn>2</mn></msubsup><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><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths>反标准化后的值即为出水BOD的软测量结果。
地址 110043 辽宁省阜新市中华路47号