基于自组织T-S模糊神经网络污泥沉降指数软测量技术

　　申请日2013.11.12

　　公开(公告)日2014.02.26

　　IPC分类号G06N3/02; G01N15/04; G06N3/08

　　摘要

　　基于自组织T-S模糊神经网络的污泥沉降指数软测量方法既属于控制领域，又属于污水处理领域。污泥沉降指数SVI的准确预测是污水处理过程正常运行的保证，本发明首先以规则层的输出量，即规则层的空间激活强度作为判定模糊规则是否增加的依据;其次，在生成新的模糊规则的基础上，以隶属函数层输出量作为判定模糊集是否增加的依据;最后，利用梯度下降算法调整模型的权值参数和高斯函数的中心值和宽度，获得一种自组织T-S模糊递归神经网络，并基于SOTSFEN建立了SVI的在线软测量模型，实现了SVI的实时检测，为预防污泥膨胀提供了一种有效方法。

　　权利要求书

　　1.一种SVI软测量方法，其特征在于，包括以下步骤：

　　(1)数据预处理及辅助变量精选;

　　样本集数据用零均值标准化方法进行归一化处理，通过主元分析PCA进行辅助变量的精选，最终确定混合液悬浮物浓度MLSS、酸碱度pH、曝气池水温T、曝气池氨NH4作为模型的输入变量;

　　(2)建立SVI的软测量的递归模糊神经网络模型，输入量为MLSS、pH、T、 NH4，模型的输出量为SVI;递归模糊神经网络拓扑结构：输入层即第一层、隶属函数层即第二层、规则层即第三层、参数层即第四层、输出层即第五层、反馈层;

　　神经网络结构：输入层有4个输入节点，每个输入层节点连接m个隶属函数层节点，规则层有m个节点，m表示模糊规则数，在网络结构训练过程中生成的模糊规则数确定m的值，参数层节点个数、反馈层节点数与规则层的节点数相等，输出层有1个节点;x=[x1,x2,x3,x4]表示神经网络的输入，yd表示神经网络的期望输出;设第k组样本数据为x(k)=[x1(k),x2(k),x3(k),x4(k)];第k组样本数据输入时：

　　输入层第i个节点的输出表示为：

　　 $o_{i}^{(1)} (k) = x_{i} (k), i = 1,2,3,4$

　　(1)

　　其中，表示输入层的第i个节点在输入第k组样本时的输出;

　　隶属函数层的节点总数为：4m，每个输入层的节点均连接m个隶属函数层的节点，隶属函数层的输出为：

　　 $o_{ij}^{(2)} (k) = \exp (- \frac{{(o_{i}^{(1)} (k) - c_{ij} (k))}^{2}}{{(σ_{ij} (k))}^{2}}), j = 1,2 . . . m$

　　(2)

　　其中，隶属函数采用高斯函数，cij和σij分别代表输入层第i个节点对应的隶属函数层第j个节点的高斯函数的中心值和宽度，每一个高斯函数均为一个模糊集，在结构调整阶段对高斯函数赋初值;表示输入层的第i个节点对应的隶属函数层的第j个节点在输入第k组样本时的输出;

　　规则层的节点数为m，输入层的每个节点对应的隶属函数层的第j个节点均连接到规则层的第j个节点;在规则层引入反馈连接，在反馈层加入内部变量; 反馈层节点包含两种类型的节点：承接节点，对规则层的输出量和权值进行加权求和运算，计算内部变量;反馈节点，采用sigmoid函数作为隶属函数，计算反馈层的输出量;规则层每个节点均连接所有的反馈层承接节点;反馈层的承接节点与反馈节点一一对应，节点数相等，反馈节点与规则层节点一一对应，节点数相等，并随着模糊规则数的变化而变化，始终等于模糊规则数;规则层第j个节点的输出为：

　　 $h_{q} = Σ_{j = 1}^{m} o_{j}^{(3)} (k - 1) ω_{jq}, j = 1,2, . . ., m; q = 1,2, . . ., m$

　　(3)

　　 $f_{q} = \frac{1}{1 + \exp ({- h}_{q})}$

　　(4)

　　 $o_{j}^{(3)} (k) = f_{q} Π_{i = 1}^{4} o_{ij}^{(2)} (k)$

　　(5)

　　其中，ωjq是规则层的第j个节点与反馈层的第q个承接节点的连接权值，初始赋值为0到1间的随机数;表示规则层第j个节点在第k-1组样本的输出量，hq表示反馈层第q个承接节点的内部变量;fq表示反馈层第q个反馈节点的输出量;规则层的输出即模糊规则的激活强度，其中表示空间激活强度，fq表示时间激活强度;

　　规则层的节点与参数层的节点是一一对应的，参数层有m个节点，该层的输出表示为：

　　 $W_{j} (k) = Σ_{i = 1}^{4} a_{ij} x_{i} (k)$

　　(6)

　　 $o_{j}^{(4)} (k) = o_{j}^{(3)} (k) W_{j} (k)$

　　(7)

　　其中，aij表示线性参数，初值赋值为0至1的随机数，aj=[a1j,a2j,a3j,a4j];Wj(k) 表示参数层第j个节点在输入第k组样本时的线性参数加权求和的值，表示参数层的第j个节点在第k组样本的输出;

　　网络模型为多输入单输出，输出层有1个节点，所有的参数层节点连接到输出节点;网络输出表示为：

　　 $y (k) = \frac{Σ_{j = 1}^{m} o_{j}^{(4)} (k)}{Σ_{j = 1}^{m} o_{j}^{(3)} (k)}$

　　(8)

　　其中，y(k)表示第k样本的网络输出;

　　(3)模糊神经网络首先通过结构自组织调整确定网络结构;网络的结构自组织调整：x(k)=[x1(k),x2(k),x3(k),x4(k)]表示模型当前输入第k组样本，从第一组样本数据开始输入，到全部数据输入完成，每输入一组样本数据均采用以下步骤判断是否增加新的模糊规则，在增加模糊规则的基础上判断是否增加新的模糊集;

　　①网络模型初始结构中模糊规则数为0，第一组数据输入，模糊规则数变为 1，并增加新的模糊集;模糊集的初始化，即高斯函数中心值和宽度值的初始化表示为：

　　c(1)=x(1)=[x1(1),x2(1),x3(1),x4(1),] (9)

　　σ(1)=[σ11,σ21,σ31,σ41]=[0.5,0.5,0.5,0.5] (10)

　　其中c(1)表示生成的第一组模糊集的中心值，即隶属函数的中心值;σ(1)表示生成的第一个模糊集的宽度值，即隶属函数的宽度值;

　　②输入数据依次输入，每输入一组数据，均判断是否增加新的模糊规则，判断公式表示为：

　　(11)

　　(12)

　　(13)

　　其中表示规则层第j个节点的空间激活强度;J表示当取最大值时j的值;N表示当前模糊规则数;为预先设定的阈值，取值为0.24;

　　如果满足式(13)，则增加一条新的模糊规则，N'=N+1，执行③;

　　如果不满足式(13)，重复②，输入下一组数据;

　　③在增加一条新的模糊规则的基础上，判断是否增加新的模糊集，判断公式表示为：

　　 $I = \arg \max_{1 < = j < = h} (o_{ij}^{(2)} (k)) - - - (14)$

　　 $o_{iI}^{(2)} (k) > I_{th}$

　　(15)

　　其中I表示当取最大值时j的值，h表示当前模型的模糊集数，h=N;Ith为预先设定的阈值，取值为0.92;

　　如果满足式(15)，则增加一条新的模糊集，h'=h+1，h'=N';对新增的模糊集赋初值即对隶属函数层的新增的高斯函数中心值和宽度值赋初值;初始化表示为：

　　cN+1=x(k)

　　(16)

　　 $c^{+} \arg \min_{1 < = p < = N} | | x (k) - c_{p} | | - - - (17)$

　　 $σ_{i, N + 1} = r | x_{i} (k) - c_{i}^{+} |, i = 1,2,3,4 - - - (18)$

　　其中，N表示当前模型已有的模糊规则数，p=1,2,...,N，cN+1表示新增隶属函数的中心值的初始值，r表示重叠系数，取值为0.6;cp表示为当前模型中第p个高斯函数的中心值;c+表示当x(k)与cp空间距离最小时cp的值;σi,N+1表示新增隶属函数的宽度初始值，输入量x(k)与c+的差值的绝对值与重叠系数的乘积;增加新的模糊集，即在每个输入层节点对应的隶属函数层增加一个新的节点，规则层、反馈层、参数层相应的各增加对应节点;

　　如果不满足(15)，则不增加新的模糊集h'=h，并且N''=N'-1，即将新增加的模糊规则删除;

　　④继续调整神经网络的结构，输入数据依次输入，重复②③，所有输入数据输入完成后，神经网络结构调整训练完成;

　　(4)确定了网络模型的结构，然后对网络参数就行调整，用校正后的数据来训练神经网络，其中共150组样本数据，训练样本数据90组，测试样本数据60 组，训练步数为1000步，每一步的训练过程中，90组训练样本数据全部输入，每组训练样本数据输入均采用梯度下降算法对网络参数进行调整，其中调整的参数包括：参数层的线性参数a，规则层到反馈层的连接权值ω，隶属函数层高斯函数的中心值c和宽度值σ;

　　定义训练的目标函数即系统误差定义为：

　　 $E (k) = \frac{1}{2} {(y (k) - y_{d} (k))}^{2}$

　　(19)

　　其中y(k)为网络的实际输出量，yd(k)为网络的期望输出量，E(k)表示系统误差;

　　梯度下降算法：采用目标函数对每一层的输出量求偏导数，计算每层输出的误差传播项;通过复合函数求偏导数，计算目标函数对参数值的偏导数，偏导数即为各参数的调整量;

　　参数层的误差传播项δ(4)，即目标函数对参数层输出的偏导数：

　　 $δ_{j}^{(4)} (k) = - \frac{\partial E (k)}{\partial y (k)} \frac{\partial y (k)}{\partial o_{j}^{(4)} (k)} - - - (20)$

　　对线性参数aj的调整量为：

　　 $\frac{\partial E (k)}{{\partial a}_{j} (k)} = \frac{\partial E (k)}{\partial o_{j}^{(4)} (k)} \frac{{\partial o}_{j}^{(4)} (k)}{\partial a_{j} (k)} = - δ_{j}^{(4)} (k) \frac{o_{j}^{(3)} (k)}{Σ_{j = 1}^{m} o_{j}^{(3)} (k)}$

　　 $(21)$ $a_{j} (k + 1) = a_{j} (k) - η \frac{\partial E (k)}{{\partial a}_{j} (k)} - - - (22)$

　　规则层误差传播项如下：

　　 $δ_{j}^{(3)} (k) = - \frac{\partial E (k)}{{\partial o}_{j}^{(3)} (k)} = δ^{(4)} (k) \frac{{\partial o}^{(4)} (k)}{{\partial o}_{j}^{(3)} (k)} - - - (23)$

　　递归层到规则层连接权值ωjq的调整规则如下：

　　 $\frac{\partial E (k)}{{\partial ω}_{jq} (k)} = \frac{\partial E (k)}{{\partial o}_{j}^{(3)} (k)} \frac{{\partial o}_{j}^{(3)} (k)}{{\partial ω}_{jq} (k)} = - δ_{j}^{(3)} (k) \frac{{\partial o}_{j}^{(3)} (k)}{{\partial ω}_{jq} (k)} - - - (24)$

　　 $ω_{jq} (k + 1) = ω_{jq} (k) - η \frac{\partial E}{{\partial ω}_{jq}} (k) - - - (25)$

　　隶属函数层的误差传播项如下：

　　 $δ_{ij}^{(2)} (k) = - \frac{\partial E (k)}{{\partial o}_{ij}^{(2)} (k)} = δ_{j}^{(3)} (k) \frac{{\partial o}_{j}^{(3)} (k)}{{\partial o}_{ij}^{(2)} (k)} - - - (26)$

　　隶属函数层高斯函数的中心值cij的调整规则如下：

　　 $\frac{\partial E (k)}{{\partial c}_{ij} (k)} = \frac{\partial E (k)}{{\partial o}_{j}^{(2)} (k)} \frac{{\partial o}_{j}^{(2)} (k)}{{\partial c}_{ij} (k)} (k) = - δ_{j}^{(2)} (k) \frac{{\partial o}_{j}^{(2)} (k)}{{\partial c}_{ij} (k)} - - - (27)$

　　 $c_{ij} (k + 1) = c_{ij} (k) - η \frac{\partial E (k)}{{\partial c}_{ij} (k)}$

　　(28)

　　隶属函数层高斯函数的宽度σij的调整规则如下：

　　 $\frac{\partial E (k)}{{\partial σ}_{ij} (k)} = \frac{\partial E (k)}{{\partial o}_{j}^{(3)} (k)} \frac{{\partial o}_{j}^{(3)} (k)}{{\partial σ}_{ij} (k)} = - δ_{j}^{(3)} (k) \frac{{\partial o}_{j}^{(3)} (k)}{{\partial σ}_{ij} (k)} - - - (29)$

　　 $σ_{ij} (k + 1) = σ_{ij} (k) - η \frac{\partial E (k)}{{\partial σ}_{ij} (k)}$

　　(30)

　　其中η为参数的学习率，取值为0.15;

　　(5)对测试样本进行预测，将测试样本数据作为训练好的神经网络的输入，神经网络的输出即为SVI的预测结果。

　　说明书

　　基于自组织T-S模糊神经网络的污泥沉降指数软测量方法

　　技术领域

　　本发明利用自组织T-S模糊递归网络建立污泥容积指数SVI的软测量模型，实现对污泥沉降指标SVI的实时预测。污泥沉降指数SVI的准确预测是污水处理过程正常运行的保证，本发明既属于控制领域，又属于污水处理领域。

　　背景技术

　　污水处理是我国政府水资源综合利用的重要举措，也是我国可持续发展战略的重要组成部分。目前，全国各城市、县基本上建立了城镇污水处理厂，污水处理能力与美国等国家发达国家相当。但污水处理运行状况不容乐观，其中污泥膨胀问题严重制约着污水处理的发展。污泥膨胀一旦发生，丝状菌大量繁殖，污泥沉降性能变差，固液分离困难，导致出水水质超标，污泥溢出流失，甚至可能引发泡沫的产生，造成污水处理系统崩溃。因此，本发明基于自组织T-S模糊递归神经网络SOTSRFNN的SVI的软测量研究具有广泛的应用前景。

　　污泥容积指数SVI是污泥沉降性能的重要评价指标之一。目前，针对SVI 的检测方法主要有两类：①人工检测法，利用量筒定时采样检测，计算SVI值，但该方法耗时且误差大，难以满足污水处理日益复杂的现实要求;②自动检测法，但该方法存在设备造价高，寿命短，稳定性差等缺点，而且受现场环境和人工操作的影响，检测精度得不到保障。软测量技术利用系统变化及参量之间的关系，建立输入输出之间的模型，通过易测水质变量估计SVI值。具有投资少、时间短、反应迅速，易于保养和维护等优点。因此，研究SOTSRFNN的软测量方法对解决SVI实时测量问题具有重要的现实意义。

　　本发明提出了一种SVI的在线软测量方法：首先，以规则层的空间激活强度，即规则层的输出作为判定模糊规则是否增加的依据;其次，在生成新的模糊规则的基础上，以隶属函数层输出量作为判定模糊集是否增加的依据;最后，利用梯度下降算法调整模型的权值参数和高斯函数的中心值和宽度，获得一种自组织T-S模糊递归神经网络，并基于SOTSRFNN建立了SVI的在线软测量模型，实现了SVI的实时检测，为预防污泥膨胀提供了一种有效方法。

　　发明内容

　　本发明针对SVI在线测量困难的问题，分析污泥膨胀的形成原因，总结与 SVI密切相关的易测水质参量，利用主元分析法PCA确定了模型的输入量;并提出了一种改进的模糊递归神经网络，基于结构自组织算法，设计了 SOTSRFNN，建立了SVI的在线软测量模型;最后，利用建立的模型进行SVI 的软测量，实现SVI的在线测量;

　　本发明采用了如下的技术方案及实现步骤：

　　1一种SVI软测量方法，其特征在于，包括以下步骤：

　　(1)数据预处理及辅助变量精选;

　　样本集数据用零均值标准化方法进行归一化处理，通过主元分析PCA进行辅助变量的精选，最终确定混合液悬浮物浓度MLSS、酸碱度pH、曝气池水温T、曝气池氨NH4作为模型的输入变量。

　　(2)建立SVI的软测量的递归模糊神经网络模型，输入量为MLSS、pH、T、 NH4，模型的输出量为SVI。递归模糊神经网络拓扑结构：输入层即第一层、隶属函数层即第二层、规则层即第三层、参数层即第四层、输出层即第五层、反馈层。

　　神经网络结构：输入层有4个输入节点，每个输入层节点连接m个隶属函数层节点，规则层有m个节点，m表示模糊规则数，在网络结构训练过程中生成的模糊规则数确定m的值，参数层节点个数、反馈层节点数与规则层的节点数相等，输出层有1个节点。x=[x1,x2,x3,x4]表示神经网络的输入，yd表示神经网络的期望输出。设第k组样本数据为x(k)=[x1(k),x2(k),x3(k),x4(k)]。第k组样本数据输入时：

　　输入层第i个节点的输出表示为：

　　 $o_{i}^{(1)} (k) = x_{i} (k), i = 1,2,3,4 - - - (1)$

　　其中，表示输入层的第i个节点在输入第k组样本时的输出;

　　隶属函数层的节点总数为：4m，每个输入层的节点均连接m个隶属函数层的节点，隶属函数层的输出为：

　　 $o_{ij}^{(2)} (k) = \exp (- \frac{{(o_{i}^{(1)} (k) - c_{ij} (k))}^{2}}{{(σ_{ij} (k))}^{2}}), j = 1,2 . . . m - - - (2)$

　　其中，隶属函数采用高斯函数，cij和σij分别代表输入层第i个节点对应的隶属函数层第j个节点的高斯函数的中心值和宽度，每一个高斯函数均为一个模糊集，在结构调整阶段对高斯函数赋初值。表示输入层的第i个节点对应的隶属函数层的第j个节点在输入第k组样本时的输出;

　　规则层的节点数为m，输入层的每个节点对应的隶属函数层的第j个节点均连接到规则层的第j个节点。在规则层引入反馈连接，在反馈层加入内部变量。反馈层节点包含两种类型的节点：承接节点，对规则层的输出量和权值进行加权求和运算，计算内部变量;反馈节点，采用sigmoid函数作为隶属函数，计算反馈层的输出量。规则层每个节点均连接所有的反馈层承接节点。反馈层的承接节点与反馈节点一一对应，节点数相等，反馈节点与规则层节点一一对应，节点数相等，并随着模糊规则数的变化而变化，始终等于模糊规则数。规则层第j个节点的输出为：

　　 $h_{q} = Σ_{j = 1}^{m} o_{j}^{(3)} (k - 1) ω_{jq}, j = 1,2, . . ., m; q = 1,2, . . ., m - - - (3)$

　　 $f_{q} = \frac{1}{1 + \exp ({- h}_{q})} - - - (4)$

　　 $o_{j}^{(3)} (k) = f_{q} Π_{i = 1}^{4} o_{ij}^{(2)} (k)$

　　(5)

　　其中，ωjq是规则层的第j个节点与反馈层的第q个承接节点的连接权值，初始赋值为0到1间的随机数。表示规则层第j个节点在第k-1组样本的输出量， hq表示反馈层第q个承接节点的内部变量。fq表示反馈层第q个反馈节点的输出量。规则层的输出即模糊规则的激活强度，其中表示空间激活强度，fq表示时间激活强度。

　　规则层的节点与参数层的节点是一一对应的，参数层有m个节点，该层的输出表示为：

　　 $W_{j} (k) = Σ_{i = 1}^{4} a_{ij} x_{i} (k) - - - (6)$

　　 $o_{j}^{(4)} (k) = o_{j}^{(3)} (k) W_{j} (k)$

　　(7)

　　其中，aij表示线性参数，初值赋值为0至1的随机数，aj=[a1j,a2j,a3j,a4j]。Wj(k)表示参数层第j个节点在输入第k组样本时的线性参数加权求和的值，表示参数层的第j个节点在第k组样本的输出。

　　网络模型为多输入单输出，输出层有1个节点，所有的参数层节点连接到输出节点。网络输出表示为：

　　 $y (k) = \frac{Σ_{j = 1}^{m} o_{j}^{(4)} (k)}{Σ_{j = 1}^{m} o_{j}^{(3)} (k)} - - - (8)$

　　其中，y(k)表示第k样本的网络输出。

　　(3)模糊神经网络首先通过结构自组织调整确定网络结构。网络的结构自组织调整：x(k)=[x1(k),x2(k),x3(k),x4(k)]表示模型当前输入第k组样本，从第一组样本数据开始输入，到全部数据输入完成，每输入一组样本数据均采用以下步骤判断是否增加新的模糊规则，在增加模糊规则的基础上判断是否增加新的模糊集。

　　①网络模型初始结构中模糊规则数为0，第一组数据输入，模糊规则数变为 1，并增加新的模糊集。模糊集的初始化，即高斯函数中心值和宽度值的初始化表示为：

　　c(1)=x(1)=[x1(1),x2(1),x3(1),x4(1),] (9)

　　σ(1)=[σ11,σ21,σ31,σ41]=[0.5,0.5,0.5,0.5] (10) 其中c(1)表示生成的第一组模糊集的中心值，即隶属函数的中心值。σ(1)表示生成的第一个模糊集的宽度值，即隶属函数的宽度值。

　　②输入数据依次输入，每输入一组数据，均判断是否增加新的模糊规则，判断公式表示为：

　　其中表示规则层第j个节点的空间激活强度。J表示当取最大值时j的值。N 表示当前模糊规则数。为预先设定的阈值，取值为0.24。

　　如果满足式(13)，则增加一条新的模糊规则，N'=N+1，执行③。

　　如果不满足式(13)，重复②，输入下一组数据;

　　③在增加一条新的模糊规则的基础上，判断是否增加新的模糊集，判断公式表示为：

　　 $I = \arg \max_{1 < = j < = h} (o_{ij}^{(2)} (k)) - - - (14)$

　　 $o_{iI}^{(2)} (k) > I_{th}$

　　(15)

　　其中I表示当取最大值时j的值，h表示当前模型的模糊集数，h=N。Ith为预先设定的阈值，取值为0.92。

　　如果满足式(15)，则增加一条新的模糊集，h'=h+1，h'=N'。对新增的模糊集赋初值即对隶属函数层的新增的高斯函数中心值和宽度值赋初值。初始化表示为：

　　cN+1=x(k) (16)

　　 $c^{+} = \arg \min_{1 < = p < = N} | | x (k) - c_{p} | | - - - (17)$

　　 $σ_{i, N + 1} = r | x_{i} (k) - c_{i}^{+} |, i = 1,2,3,4 - - - (18)$

　　其中，N表示当前模型已有的模糊规则数，p=1,2,...,N，cN+1表示新增隶属函数的中心值的初始值，r表示重叠系数，取值为0.6。cp表示为当前模型中第p个高斯函数的中心值。c+表示当x(k)与cp空间距离最小时cp的值。σi,N+1表示新增隶属函数的宽度初始值，输入量x(k)与c+的差值的绝对值与重叠系数的乘积。增加新的模糊集，即在每个输入层节点对应的隶属函数层增加一个新的节点，规则层、反馈层、参数层相应的各增加对应节点。

　　如果不满足(15)，则不增加新的模糊集h'=h，并且N''=N'-1，即将新增加的模糊规则删除。

　　④继续调整神经网络的结构，输入数据依次输入，重复②③，所有输入数据输入完成后，神经网络结构调整训练完成。

　　定义训练的目标函数即系统误差定义为：

　　 $E (k) = \frac{1}{2} {(y (k) - y_{d} (k))}^{2} - - - (19)$

　　其中y(k)为网络的实际输出量，yd(k)为网络的期望输出量，E(k)表示系统误差。

　　梯度下降算法：采用目标函数对每一层的输出量求偏导数，计算每层输出的误差传播项。通过复合函数求偏导数，计算目标函数对参数值的偏导数，偏导数即为各参数的调整量。

　　参数层的误差传播项δ(4)，即目标函数对参数层输出的偏导数：

　　 $δ_{j}^{(4)} (k) = - \frac{\partial E (k)}{\partial y (k)} \frac{\partial y (k)}{\partial o_{j}^{(4)} (k)} - - - (20)$

　　对线性参数aj的调整量为：

　　 $a_{j} (k + 1) {= a}_{j} (k) - η \frac{\partial E (k)}{{\partial a}_{j} (k)} - - - (22)$

　　规则层误差传播项如下：

　　 $δ_{j}^{(3)} (k) = - \frac{\partial E (k)}{{\partial o}_{j}^{(3)} (k)} = δ^{(4)} (k) \frac{{\partial o}^{(4)} (k)}{{\partial o}_{j}^{(3)} (k)} - - - (23)$

　　递归层到规则层连接权值ωjq的调整规则如下：

　　 $ω_{jq} (k + 1) = ω_{jq} (k) - η \frac{\partial E}{{\partial ω}_{jq}} (k) - - - (25)$

　　隶属函数层的误差传播项如下：

　　 $δ_{ij}^{(2)} (k) = - \frac{\partial E (k)}{{\partial o}_{ij}^{(2)} (k)} = δ_{j}^{(3)} (k) \frac{{\partial o}_{j}^{(3)} (k)}{{\partial o}_{ij}^{(2)} (k)} - - - (26)$

　　隶属函数层高斯函数的中心值cij的调整规则如下：

　　 $c_{ij} (k + 1) = c_{ij} (k) - η \frac{\partial E (k)}{{\partial c}_{ij} (k)} - - - (28)$

　　隶属函数层高斯函数的宽度σij的调整规则如下：

　　 $σ_{ij} (k + 1) = σ_{ij} (k) - η \frac{\partial E (k)}{{\partial σ}_{ij} (k)} - - - (30)$

　　其中η为参数的学习率，取值为0.15。

　　(5)对测试样本进行预测，将测试样本数据作为训练好的神经网络的输入，神经网络的输出即为SVI的预测结果。

　　本发明的创造性主要体现在：

　　(1).本发明针对当前污泥容积指数SVI难以在线检测的问题，提出了基于一种自组织T-S模糊递归网络的软测量方法，实现了辅助变量和SVI之间的映射关系及SVI的在线检测;该模型具有精度高、稳定性好、实时性强等特点，为预防污泥膨胀提供了一种有效的检测方法，促进了污水处理自动化水平的提高;

　　(2).本发明针对污泥膨胀的特性，对模糊递归神经网络的反馈层加以改进。基于规则产生准则，自动的生成模糊规则，通过有效的模糊集生成算法，动态调整网络结构。解决了网络结构难以确定的问题，在保证模型精度的同时，有效的简化了网络结构。采用梯度下降参数学习算法，提高了网络的学习能力。