发明名称 |
基于动态径向基神经网络的溶解氧的控制方法 |
摘要 |
本发明公开了一种基于动态径向基神经网络的溶解氧的控制方法,采用以下步骤:确定控制对象;设计用于污水处理过程中溶解氧DO控制器的动态RBF神经网络拓扑结构;对样本数据进行校正;用校正后的部分数据训练神经网络,利用训练好的RBF神经网络对溶解氧(DO)进行控制,期望DO浓度与实际输出DO浓度的误差及误差变化率作为RBF神经网络的输入,RBF神经网络的输出即为变频器的输入,变频器通过调节电动机的转速达到控制鼓风机的目的,最终控制曝气量,整个控制系统的输出为实际DO浓度值,提高控制器的控制效果,能够快速、准确地使溶解氧达到期望要求;解决了当前基于开关控制和PID控制自适应能力较差的问题。 |
申请公布号 |
CN101576734B |
申请公布日期 |
2010.12.08 |
申请号 |
CN200910086991.9 |
申请日期 |
2009.06.12 |
申请人 |
北京工业大学 |
发明人 |
乔俊飞;韩红桂;郭迎春 |
分类号 |
G05B13/00(2006.01)I;C02F3/12(2006.01)I |
主分类号 |
G05B13/00(2006.01)I |
代理机构 |
北京思海天达知识产权代理有限公司 11203 |
代理人 |
张慧 |
主权项 |
1.一种基于动态径向基神经网络的溶解氧的控制方法,其特征在于,包括以下步骤:(1)确定控制对象;针对序批式间歇活性污泥系统中溶解氧进行控制,以曝气量为控制量,溶解氧浓度为被控量;(2)设计用于污水处理过程中溶解氧DO控制器的动态RBF神经网络拓扑结构;网络分为三层:输入层、隐含层、输出层;初始化RBF神经网络:确定神经网络p-n-1的连接方式,即输入层神经元为p个,隐含层神经元为n个,输出层神经元为1个;对神经网络的权值进行随机赋值;RBF神经网络的输入表示为x<sub>1</sub>,x<sub>2</sub>,…,x<sub>p</sub>,神经网络的期望输出表示为y<sub>d</sub>,实际输出表示为y;设第k时刻RBF神经网络输入为x<sub>1</sub>(k),x<sub>2</sub>(k),…,x<sub>p</sub>(k),RBF神经网络各层的计算功能是:输入层,该层由p个神经元组成:<img file="FSB00000108139200011.GIF" wi="530" he="69" /><img file="FSB00000108139200012.GIF" wi="620" he="70" />i=1,2,…,p;(1)In<sub>i</sub><sup>(1)</sup>(k),Out<sub>i</sub><sup>(1)</sup>(k)分别表示输入层的输入和输出;隐含层,该层由n个神经元组成:<img file="FSB00000108139200013.GIF" wi="744" he="87" /><img file="FSB00000108139200014.GIF" wi="608" he="75" />i=1,2,…,p;j=1,2,…,n; (2)In<sub>j</sub><sup>(2)</sup>(k),Out<sub>j</sub><sup>(2)</sup>(k)表示隐含层的输入和输出,x(k)=[x<sub>1</sub>(k),x<sub>2</sub>(k),…,x<sub>p</sub>(k)]<sup>T</sup>表示第k时刻的输入值,c<sub>j</sub>表示隐含层第j个神经元的中心值,||x(k)-c<sub>j</sub>||表示c<sub>j</sub>与x(k)之间的欧式距离,<img file="FSB00000108139200015.GIF" wi="100" he="59" />为高斯函数,其形式为:<img file="FSB00000108139200016.GIF" wi="1459" he="160" />c<sub>j</sub>表示隐含层第j个神经元中心值,δ<sub>j</sub>表示隐含层第j个神经元的中心宽度;输出层,该层只有1个神经元:<img file="FSB00000108139200017.GIF" wi="864" he="125" />j=1,2,…,n (4)Out<sup>(3)</sup>(k)表示输出层的输出,w<sub>j</sub><sup>3</sup>(k)表示第j个隐含层神经元与输出层神经元之间的联结权值,y(k)为神经网络的实际输出;定义误差函数为<img file="FSB00000108139200018.GIF" wi="1274" he="131" />T为(y(k)-y<sub>d</sub>(k))的转置,m为训练样本总数,训练神经网络的目的是使得式(5)定义的误差函数达到期望值E<sub>d</sub>;(3)对样本数据进行校正;设t个数据样本x(1),x(2),…,x(t),均值为<img file="FSB00000108139200021.GIF" wi="70" he="60" />每一个样本的偏差为<img file="FSB00000108139200022.GIF" wi="382" he="69" />q=1,2,…,t,按照Bessel公式计算出标准偏差:<img file="FSB00000108139200023.GIF" wi="1183" he="159" />若某一个样本x(q)的偏差满足:|D(q)|≥3σ,q=1,2,…,t;(7)则认为样本x(q)是异常数据,应予以剔除,得到校正后的数据,该数据作为神经网络的训练样本;(4)用校正后的部分数据训练神经网络,在训练过程中利用模型输出敏感度分析对神经网络的冗余隐神经元进行修剪,利用最临近插值法对隐含层神经元进行增加;以简化神经网络结构,提高神经网络的性能;具体为:①训练给定初始神经网络达到预先设定训练步骤m;记下训练过程中每个隐含层神经元的输出权值w<sub>1</sub><sup>3</sup>,w<sub>2</sub><sup>3</sup>,…,w<sub>n</sub><sup>3</sup>,找出最大值和最小值;通过计算参数w<sub>1</sub><sup>3</sup>,w<sub>2</sub><sup>3</sup>,…,w<sub>n</sub><sup>3</sup>对于网络输出y的灵敏度;②第j个隐神经元输出在训练m步中的最大值是b<sub>j</sub>,最小值是a<sub>j</sub>,则假定w<sub>j</sub><sup>3</sup>以指定频率<img file="FSB00000108139200024.GIF" wi="61" he="54" />在[a<sub>j</sub>,b<sub>j</sub>]内振荡,即:<img file="FSB00000108139200025.GIF" wi="1284" he="81" />其中,j=1,2,…,n,s<sup>(j)</sup>=2π/m,<img file="FSB00000108139200026.GIF" wi="61" he="54" />指定频率,m是总训练步骤;结合式(4),神经网络的输出y表示为:<img file="FSB00000108139200027.GIF" wi="841" he="124" />j=1,2,…,n (9)③计算各个神经元的灵敏度:如果要计算第j个隐神经元的灵敏度,则该神经元对应的傅立叶系数为:<img file="FSB00000108139200028.GIF" wi="657" he="140" /><img file="FSB00000108139200029.GIF" wi="866" he="175" />其中,j=1,2,…,n,s的取值范围是[-π,π];由RBF神经网络的各隐含层神经元的输出之间没有相互作用,傅立叶振幅值主要集中在基频上,采用下式计算第j个隐神经元的灵敏度;<img file="DEST_PATH_FSB00000198866000011.GIF" wi="1163" he="384" />则j个隐神经元的总灵敏度之和为:<img file="DEST_PATH_FSB00000198866000012.GIF" wi="1062" he="190" />第j个输入参数的灵敏度为:S<sub>j</sub>=ST<sub>j</sub>/Sum_ST (13)④进行结构调整,若第j个隐含层神经元的显著性S<sub>j</sub>小于ε<sub>1</sub>,则删除该神经元;若第j个隐含层神经元的显著性S<sub>j</sub>大于ε<sub>2</sub>,利用最临近插值法在该神经元附近插入新的神经元,其中ε<sub>1</sub>和ε<sub>2</sub>预先设定;插入的新神经元初始条件:<img file="DEST_PATH_FSB00000198866000013.GIF" wi="278" he="138" />v<sub>new</sub>=v<sub>j</sub>,c<sub>new</sub>=c<sub>j</sub> (14)其中,<img file="DEST_PATH_FSB00000198866000014.GIF" wi="116" he="68" />v<sub>new</sub>,c<sub>new</sub>分别表示插入新神经元的联结权值,中心值,中心宽度;⑤继续训练神经网络,然后每隔一定步数重复步骤②-④,直到新神经网络中所有隐含层神经元显著性S<sub>j</sub>处于ε<sub>1</sub>和ε<sub>2</sub>之间停止结构调整;⑥调整校正网络中所有的参数:<img file="DEST_PATH_FSB00000198866000015.GIF" wi="645" he="143" /><img file="DEST_PATH_FSB00000198866000016.GIF" wi="627" he="143" /><img file="DEST_PATH_FSB00000198866000017.GIF" wi="692" he="143" />(15)其中,i=1,2,…,n,η<sub>1</sub>,η<sub>2</sub>,η<sub>3</sub>是正常数,它们取值范围是[0,1];⑦继续对样本数据进行训练,直到误差E小于期望误差E<sub>d</sub>;(5)利用训练好的RBF神经网络对溶解氧DO进行控制,控制系统的输出为实际DO浓度。 |
地址 |
100124 北京市平乐园100号 |