Intelligent detection method for Biochemical Oxygen Demand based on a Self-organizing Recurrent RBF Neural Network

Abstract

Under conventional techniques, wastewater treatment has many problems such as poor production conditions, serious random interference, strong nonlinear behavior, large time-varying, and serious lagging. These problem cause difficult detection of various wastewater treatment parameter such as biochemical oxygen demand (BOD) values that are used to monitor water quality. To solve problems associated with monitoring BOD values in real-time, the present disclosure utilizes a self-organizing recurrent RBF neural network designed for intelligent detecting of BOD values. Implementations of the present disclosure build a computing model of BOD values based on the self-organizing recurrent RBF neural network to achieve real-time and more accurate detection of the BOD values (e.g., a BOD concentration). The implementations herein quickly and accurately obtain BOD concentrations and improve the quality and efficiency of wastewater treatment.

Claims

1. A method for intelligent detection of a biochemical oxygen demand (BOD) concentration, the method comprising: (1) determining input and output variables of BOD by analyzing sewage treatment process variables and selecting the input variables of a BOD computing model that include a dissolved oxygen concentration (DO), an effluent suspended solids concentration (SS), a pH value, a chemical oxygen demand (COD) for a sewage treatment process of an activated sludge system, an output value of the BOD computing model being a detected BOD concentration associated with the sewage treatment process of the activated sludge system; (2) designing the BOD computing model of the intelligent detection of the BOD concentration using a self-organizing recurrent Radial basis function (RBF) neural network; (3) training the self-organizing recurrent RBF neural network using training samples; and (4) predicting the BOD concentration associated with the sewage treatment process of the activated sludge system using the trained self-organizing recurrent RBF neural network and testing sample data.

2. The method of claim 2, wherein a structure of the self-organizing recurrent RBF neural network comprises three layers: an input layer, a hidden layer and an output layer, a connection manner of the self-organizing recurrent RBF neural network is 4-m-1 in which a number of the input layer is 4 and the hidden neurons is m, connection weights between the input layer and the hidden layer are assigned 1, the connection weights between the hidden layer and the output layer are randomly assigned values, the assignment interval is [1, 1], a number of training samples is N, and an input of the self-organizing recurrent RBF neural network is x(t)=[x.sub.1(t), x.sub.2(t), x.sub.3(t), x.sub.4(t)] at time t.

3. The method of claim 3, wherein an expectation output of the self-organizing recurrent RBF neural network is expressed as y.sub.d(t) and an actual output of the self-organizing recurrent RBF neural network is expressed as y(t), and the BOD concentration is calculated based on the self-organizing recurrent RBF neural network comprising: {circle around (1)} the input layer comprising 4 neurons representing the input variables in the input layer, and the output values of each neuron are as follows:
u.sub.i(t)=x.sub.i(t); (Equation 1) wherein u.sub.i(t) is a ith output value at time t, i=1, 2, . . . , 4, and an input vector is x(t)=[x.sub.1(t), x.sub.2(t), . . . , x.sub.4(t)]. {circle around (2)} the hidden layer comprising m neurons of the hidden layer, the outputs of hidden neurons are: $\begin{array}{cc}{}_j\left(t\right)=e^{-\frac{{.Math.h_j\left(t\right)-c_j\left(t\right).Math.}^2}{2.Math.{}_j^2\left(t\right)}},.Math.j=1,2,\mathrm{.Math.}.Math.,m,& \left(\mathrm{Equation}.Math..Math.2\right)\end{array}$ c.sub.j(t) denotes the center vector of the jth hidden neuron and c.sub.j(t)=[c.sub.1j(t), c.sub.2j(t), . . . , c.sub.4j(t)].sup.T at time t, h.sub.j(t)c.sub.j(t) is the Euclidean distance between h.sub.j(t) and c.sub.j(t), .sub.j(t) is the radius or width of the jth hidden neuron at time t, and h.sub.j(t) is input vector of the jth hidden neuron at time t described as
h.sub.j(t)=[u.sub.1(t), u.sub.2(t), . . . , u.sub.4(t), v.sub.j(t)y(t1)].sup.T, (Equation 4) wherein y(t1) is the output value of the output layer at time t1, v.sub.j(t) denotes the connection weight from output layer to the jth hidden neuron at time t, and v(t)=[v.sub.1(t), v.sub.2(t), . . . , v.sub.m(t)].sup.T, T represents transpose; {circle around (3)} the output layer comprising one node, the output is: $\begin{array}{cc}y\left(t\right)=f\left(w\left(t\right),\left(t\right)\right)=\underset{j=1}{\overset{m}{.Math.}}.Math.w_j\left(t\right){}_j\left(t\right),.Math.j=1,\mathrm{.Math.}.Math.,m,,& \left(\mathrm{Equation}.Math..Math.4\right)\end{array}$ wherein w(t)=[w.sub.1(t), w.sub.2(t), . . . , w.sub.m(t)].sup.T is the connection weights between the hidden neurons and output neuron at time t, (t)=[0.sub.1(t), 0.sub.2(t), . . . , 0.sub.m(t)].sup.T is the output vector of the hidden layer, and y(t) represents the output of recurrent RBF neural network at time t.

4. The method of claim 3, wherein an error of the self-organizing recurrent RBF neural network is calculated using: $\begin{array}{cc}E\left(t\right)=\frac{1}{N}.Math.\underset{t=1}{\overset{N}{.Math.}}.Math.{\left(y_d\left(t\right)-y\left(t\right)\right)}^2,& \left(\mathrm{Equation}.Math..Math.5\right)\end{array}$ wherein y.sub.d(t) is the expectation output of neural network and the actual output is expressed as y(t).

5. The method of claim 4, wherein the training the self-organizing recurrent RBF neural network using training samples comprising: {circle around (1)} providing the self-organizing recurrent RBF neural network, an initial number of the hidden layer neurons is m, m>2 is a positive integer, the input of the self-organizing recurrent RBF neural network is x(1), x(2), . . . , x(t), . . . , x(N), correspondingly, the expectation output is y.sub.d(1), y.sub.d(2), . . . , y.sub.d(t), . . . , y.sub.d(N), an expected error value is set to E.sub.d, E.sub.d(0, 0.01), each variable of initial center values c.sub.j(1) (-2, 2), a width value .sub.j(1) (0, 1), an initial feedback weight v.sub.j(1) (0, 1), j=1, 2, . . . , m; and an initial weight w(1) (0, 1); {circle around (2)} setting a learning step s=1; {circle around (3)} t=s, calculating the output of the self-organizing recurrent RBF neural network according to equations (1)-(4) and exploiting a fast gradient descent algorithm using parameters: $\begin{array}{cc}{c}_j\left(t+1\right)=c_j\left(t\right)-{}_c.Math.\frac{1}{{}_j^2}.Math.\left(y_d\left(t\right)-y\left(t\right)\right).Math.w_j\left(t\right)\left(t\right)\left[h_j\left(t\right)-c_j\left(t\right)\right],& \left(\mathrm{Equation}.Math..Math.6\right)\\ {}_j\left(t+1\right)={}_j\left(t\right)-{}_{}.Math.\frac{1}{{}_j^3}.Math.\left(y_d\left(t\right)-y\left(t\right)\right).Math.w_j\left(t\right)w_j\left(t\right)\left(t\right).Math.{.Math.h_j\left(t\right)-c_j\left(t\right).Math.}^2,& \left(\mathrm{Equation}.Math..Math.7\right)\\ .Math.v_j\left(t+1\right)=v_j\left(t\right)-{}_v\left(y_d\left(t\right)-y\left(t\right)\right).Math.w_j\left(t\right).Math.\left(t\right).Math.y\left(t-1\right),& \left(\mathrm{Equation}.Math..Math.8\right)\\ .Math.w_j\left(t+1\right)=w_j\left(t\right)-{}_w\left(y_d\left(t\right)-y\left(t\right)\right).Math.{}_j\left(t\right),& \left(\mathrm{Equation}.Math..Math.9\right)\end{array}$ wherein, .sub.c, .sub., .sub.v, .sub.w are the learning rate of center, width, feedback connection weight from the output layer to the hidden layer and the connection weight between the hidden layer and the output layer, respectively, .sub.c(0, 0.01], .sub.(0, 0.01], .sub.v(0, 0.02], .sub.w (0, 0.01], c.sub.j(t+1)=[c.sub.1j(t+1), c.sub.2j(t+1), . . . , c.sub.4j(t+1)] denotes the center vector of the jth hidden neuron at time t+1, .sub.j(t+1) is the radius or width of the jth hidden neuron at time t+1, v.sub.j(t+1) denotes the connection weight from output layer to the jth hidden neuron at time t+1, and w.sub.j(t+1) is the connection weights between the hidden neurons and output neuron at time t+1; {circle around (4)} t>3, calculating independent contribution: $\begin{array}{cc}{}_j\left(t\right)=\frac{q_j\left(t-1\right)+q_j\left(t\right)}{\underset{j=1}{\overset{m}{.Math.}}.Math.\left(q_j\left(t-1\right)+q_j\left(t\right)\right)},.Math.j=1,\mathrm{.Math.}.Math.,m,& \left(\mathrm{Equation}.Math..Math.10\right)\end{array}$ wherein .sub.j(t) is the independent contribution of the jth hidden neuron at time t, q.sub.j(t1) is independent contribution output of the jth hidden neuron at time t1, q.sub.j(t) is independent contribution output of the jth hidden neuron at time t, q.sub.j=[q.sub.j(t1), q.sub.j(t)] is independent contribution output vector of the jth hidden neuron, Q(t)=[q.sub.1(t), . . . q.sub.m1(t), q.sub.m(t)].sup.T is the independent contribution matrix at time t,
Q(t)=(t)(t), (Equation 11) wherein (t) is a coefficients matrix which is provided as:
(t)=D.sup.1(t)(t)B(t)z(t), (Equation 12) wherein (t)=[(t1), (t)] is output matrix of hidden layer at time t, 0(t1)=[.sub.1(t1), .sub.2(t1), . . . , .sub.m(t1)].sup.T is output vector of hidden layer at time t1, (t)=[.sub.1(t), .sub.2(t), . . . , .sub.m(t)].sup.T is output vector of hidden layer at time t, D(t), B(t) and z(t) are the covariance matrix of (t), the whitening matrix of y(t) and the whitening transformation matrix of y(t), respectively, and D(t), B(t) and z(t) are provided as: $\begin{array}{cc}D\left(t\right)=\left[\begin{array}{cc}\frac{\underset{j=1}{\overset{m}{.Math.}}.Math.{\left({}_j\left(t-1\right)-\stackrel{\_}{}\left(t-1\right)\right)}^2}{m-1}& \frac{\underset{j=1}{\overset{m}{.Math.}}.Math.\left({}_j\left(t-1\right)-\stackrel{\_}{}\left(t-1\right)\right).Math.\left({}_j\left(t\right)-\stackrel{\_}{}\left(t\right)\right)}{m-1}\\ \frac{\underset{j=1}{\overset{m}{.Math.}}.Math.\left({}_j\left(t-1\right)-\stackrel{\_}{}\left(t-1\right)\right).Math.\left({}_j\left(t\right)-\stackrel{\_}{}\left(t\right)\right)}{m-1}& \frac{\underset{j=1}{\overset{m}{.Math.}}.Math.{\left({}_j\left(t\right)-\stackrel{\_}{}\left(t\right)\right)}^2}{m-1}\end{array}\right],& \left(\mathrm{Equation}.Math..Math.13\right)\\ B\left(t\right)={}^{-1/2}\left(t\right).Math.U^T\left(t\right),& \left(\mathrm{Equation}.Math..Math.14\right)\\ z\left(t\right)={}^{-1/2}\left(t\right).Math.U^T\left(t\right).Math.y\left(t\right),& \left(\mathrm{Equation}.Math..Math.15\right)\end{array}$ wherein (t1)=.sub.1(t1)+.sub.2(t1)+. . . +.sub.m(t1))/m is the average value of elements of output vector in the hidden layer at time t1, (t)=(.sub.1(t)+.sub.2(t)+. . . +.sub.m(t))/m is the average value of elements of the output vector in the hidden layer at time t, U(t) and (t) are the eigenvector and eigenvalue matrices of y(t), y(t) is the output matrix of the self-organizing recurrent RBF neural network at time t,
y(t)=(t)(t), (Equation 16) wherein (t) is the weight matrix of hidden layer to the output layer,
(t)=[w(t1),w(t)], (Equation 17) wherein w(t1)=[w.sub.1(t1), w.sub.2(t1), . . . , w.sub.m(t1)].sup.T and w(t)=[w.sub.1(t), w.sub.2(t), . . . , w.sub.m(t)].sup.T are the output of self-organizing recurrent RBF neural network, the output vector of the hidden layer and the weight vector at time t1 and time t, respectively, {circle around (5)} t>3, calculating activity degree of the hidden neuron:
S.sub.j(t)=e.sup.h.sub.j.sup.(t)c.sub.j.sup.(t), (Equation 18) wherein S.sub.j(t) is the activity degree of the jth hidden neuron at time t, j=1, 2, . . . , m. {circle around (6)} t>3, adjusting the structure of the self-organizing recurrent RBF neural network by calculating the activity degree of the lth hidden neuron S(t) and the independent contribution of the lth hidden neuron .sub.l(t), wherein the activity degree and independent contribution of the lth hidden neuron satisfy:
S.sub.l(t)=max S(t), (Equation 19)
.sub.l(t)=max (t), (Equation 20) wherein S(t)=[S.sub.1(t), . . . , S.sub.m1(t), S.sub.m(t)] is the vector of activity degree of hidden neurons at time t, (t)=[.sub.1(t), . . . , .sub.m1(t), .sub.m(t)] is the vector of independent contribution of hidden neurons at time t, an hidden neuron is added and the number of hidden neurons is updated using M.sub.1=m+1 or the structure of self-organizing recurrent RBF neural network is not adjusted, M.sub.1=m, when the activity degree and independent contribution of the ith hidden neuron satisfy:
S.sub.i(t)=min S(t), (Equation 21)
.sub.i(t)=min (t), (Equation 22) wherein the ith hidden neuron is deleted and the number of the hidden neurons is updated using M.sub.2=M.sub.11 or the structure of self-organizing recurrent RBF neural network is not adjusted, M.sub.2=M.sub.1; {circle around (7)} increasing a learning step s: if s<N, turning to step {circle around (3)}; if s=N, turning to step {circle around (8)}; and {circle around (8)} calculating the performance of self-organizing recurrent RBF neural network according to Equation 5: if E(t)E.sub.d, then turning to step {circle around (3)}; and if E(t)<E.sub.d, stopping the training process.

6. The method of claim 5, wherein the testing sample data is the input of the self-organizing recurrent RBF neural network, and the output of the self-organizing recurrent RBF neural network is a computed value of the BOD concentration.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0049] The detailed description is described with reference to the accompanying figures. In the figures, the left-most digit(s) of a reference number identifies the figure in which the reference number first appears. The same reference numbers in different figures indicate similar or identical items.

[0050] FIG. 1 shows a structure of an intelligent detection method based on a self-organizing recurrent RBF neural network in accordance with implementations of this disclosure.

[0051] FIG. 2 shows training results of an intelligent detection method.

[0052] FIG. 3 shows a training error of an intelligent detection method.

[0053] FIG. 4 shows a predicting result of an intelligent detection method.

[0054] FIG. 5 shows a predicting error of an intelligent detection method.

[0055] FIGS. 6-17 show tables in accordance with implementations of the present disclosure. Tables 1-12 show the experimental data in this disclosure. Tables 1-4 show training samples of chemical oxygen demand (COD), DO, pH, SS concentration. Table 5 shows real output values of BOD concentrations. Table 6 shows outputs of the self-organizing recurrent RBF neural network in the training process. Tables 7-10 show the testing samples of chemical oxygen demand (COD), DO, pH, SS concentration. Table 11 shows real output values of BOD concentration. Table 12 shows the outputs of the self-organizing recurrent RBF neural network in the predicting process. Moreover, the samples are imported as the sequence from the tables. The first data is in the first row and the first column. Then, the second data is in the first row and the second column. Until all of data is imported from the first row, the data in the second row and following rows are inputted as the same way.

DETAILED DESCRIPTION

[0056] This disclosure takes suspended solids concentrations (SS), dissolved oxygen (DO), pH, chemical oxygen demand (COD) as characteristic variables for BOD, except for the pH (no unit), the unit of the above parameters is mg/L;

[0057] The experimental data comes from water quality analysis statement of a wastewater treatment plant in 2012; choosing data of SS concentrations, DO, pH and COD as experimental samples, after eliminating abnormal sample, 100 groups of data are available, and the group of 60 are used as training samples, and the remaining 40 groups are used as test samples.

[0058] This disclosure adopts the following technical scheme and implementation steps:

[0059] An intelligent detection method for the BOD concentration based on a self-organizing recurrent RBF neural network is described using the following operations.

[0060] (1) Determining the input and output variables of BOD: For sewage treatment process of activated sludge system, sewage treatment process variables are analyzed and select the input variables of BOD computing model: DO concentration, effluent SS concentration, pH, COD. The output value of computing model is the detected BOD concentration.

[0061] (2) Computing model of intelligent detection of BOD: establishing a computing model of BOD based on a self-organizing recurrent RBF neural network. The structure of recurrent RBF neural network may include three layers: input layer, hidden layer and output layer. The network is 4-m-1, named the number of input layer is 4 and hidden neurons is m. Connection weights between input layer and hidden layer are assigned one, the connection weights between hidden layer and output layer randomly assign values, the assignment interval is [1, 1]. The number of the training sample is N and the input of self-organizing recurrent RBF neural network is x(t)=[x.sub.1(t), x.sub.2(t), x.sub.3(t), x.sub.4(t)] at time t. The expectation output of neural network output is expressed as y.sub.d(t) and the actual output is expressed as y(t). Computing method of BOD can be described:

[0062] {circle around (1)} The input Layer: There are 4 neurons which represent the input variables in this layer. The output values of each neuron are as follows:

u.sub.i(t)=x.sub.i(t); (Equation 23)

[0063] wherein u.sub.i(t) is the ith output value at time t, i=1, 2, . . . , 4, and the input vector is x(t)=[x.sub.1(t), x.sub.2(t), . . . x.sub.4(t)].

[0064] {circle around (2)} The Hidden Layer: There are m neurons of hidden layer. The outputs of hidden neurons are:

[00007]$\begin{array}{cc}{}_j\left(t\right)=e^{-\frac{{.Math.h_j\left(t\right)-c_j\left(t\right).Math.}^2}{2.Math.{}_j^2\left(t\right)}},.Math.j=1,2,\mathrm{.Math.}.Math.,m,& \left(\mathrm{Equation}.Math..Math.24\right)\end{array}$

[0065] c.sub.j (t) denotes the center vector of the jth hidden neuron and c.sub.j(t)=[c.sub.1j(t), c.sub.2j(t), . . . , c.sub.4j(t)].sup.T at time t, h.sub.j(t)c.sub.j(t) is the Euclidean distance between h.sub.j(t) and c.sub.j(t), and .sub.j(t) is the radius or width of the jth hidden neuron at time t, h.sub.j(t) is input vector of the jth hidden neuron at time t described as

h.sub.j(t)=[u.sub.1(t), u.sub.2(t), . . . u.sub.4(t), v.sub.j(t)y(t1)].sup.T, (Equation 25)

[0066] y(t1) is the output value of the output layer at time t1, v.sub.j(t) denotes the connection weight from output layer to the jth hidden neuron at time t, and v(t)=[v.sub.1(t), v.sub.2(t), . . . , v.sub.m(t)].sup.T, T represents transpose;

[0067] {circle around (3)} The Output Layer: There is only one node in this layer, the output is:

[00008]$\begin{array}{cc}y\left(t\right)=f\left(w\left(t\right),\left(t\right)\right)=\underset{j=1}{\overset{m}{.Math.}}.Math.w_j\left(t\right){}_j\left(t\right),.Math.j=1,\mathrm{.Math.}.Math.,m,,& \left(\mathrm{Equation}.Math..Math.26\right)\end{array}$

[0068] wherein w(t)=[w.sub.1(t), w.sub.2(t), . . . , w.sub.m(t)].sup.T is the connection weights between the hidden neurons and output neuron at time t, (t)=[.sub.1(t), .sub.2(t), . . . , .sub.m(t)].sup.T is the output vector of the hidden layer, y(t) represents the output of recurrent RBF neural network at time t.

[0069] The error of self-organizing recurrent RBF neural network is:

[00009]$\begin{array}{cc}E\left(t\right)=\frac{1}{N}.Math.\underset{t=1}{\overset{N}{.Math.}}.Math.{\left(y_d\left(t\right)-y\left(t\right)\right)}^2,& \left(\mathrm{Equation}.Math..Math.27\right)\end{array}$

[0070] y.sub.d(t) is the expectation output of neural network and the actual output is expressed as y(t);

[0071] (3) Training the self-organizing recurrent RBF neural network;

[0072] {circle around (1)} Providing the self-organizing recurrent RBF neural network, the initial number of hidden layer neurons is m, and m>2 is a positive integer. The input of self-organizing recurrent RBF neural network is x(1), x(2), . . . , x(t), . . . , x(N), correspondingly, the expectation output is y.sub.d(1), y.sub.d(2), . . . , y.sub.d(t), . . . , y.sub.d(N), expected error value is set to E.sub.d, E.sub.d(0, 0.01). The every variable of initial centre value c.sub.j(1) (2, 2), width value .sub.j(1) (0, 1), initial feedback weight v.sub.j(1) (0, 1), j=1, 2, . . . , m; initial weight w(1) (0, 1);

[0073] {circle around (2)} Setting the learning step s=1;

[0074] {circle around (3)} t=s; According to Equations (1)-(4), calculating the output of self-organizing recurrent RBF neural network by exploiting a fast gradient descent algorithm:

[00010]$\begin{array}{cc}{c}_j\left(t+1\right)=c_j\left(t\right)-{}_c.Math.\frac{1}{{}_j^2}.Math.\left(y_d\left(t\right)-y\left(t\right)\right).Math.w_j\left(t\right)\left(t\right)\left[h_j\left(t\right)-c_j\left(t\right)\right],& \left(\mathrm{Equation}.Math..Math.28\right)\\ {}_j\left(t+1\right)={}_j\left(t\right)-{}_{}.Math.\frac{1}{{}_j^3}.Math.\left(y_d\left(t\right)-y\left(t\right)\right).Math.w_j\left(t\right)w_j\left(t\right)\left(t\right).Math.{.Math.h_j\left(t\right)-c_j\left(t\right).Math.}^2,& \left(\mathrm{Equation}.Math..Math.29\right)\\ .Math.v_j\left(t+1\right)=v_j\left(t\right)-{}_v\left(y_d\left(t\right)-y\left(t\right)\right).Math.w_j\left(t\right).Math.\left(t\right).Math.y\left(t-1\right),& \left(\mathrm{Equation}.Math..Math.30\right)\\ .Math.w_j\left(t+1\right)=w_j\left(t\right)-{}_w\left(y_d\left(t\right)-y\left(t\right)\right).Math.{}_j\left(t\right),& \left(\mathrm{Equation}.Math..Math.31\right)\end{array}$

[0075] .sub.c, .sub., .sub.v, .sub.w are the learning rate of center, width, feedback connection weight from output layer to hidden layer and the connection weight between hidden layer and output layer, respectively. In addition, .sub.c(0, 0.01], .sub.(0, 0.01], .sub.v(0, 0.02], .sub.w(0, 0.01]; c.sub.j(t1)=[c.sub.1j(t+1), c.sub.2j(t+1), c.sub.4j(t+1)] denotes the center vector of the jth hidden neuron at time t+1; .sub.j(t+1) is the radius or width of the jth hidden neuron at time t+1; v.sub.j(t+1) denotes the connection weight from output layer to the jth hidden neuron at time t+1; w.sub.j(t+1) is the connection weights between the hidden neurons and output neuron at time t+1;

[0076] {circle around (4)} t>3, calculating independent contribution:

[00011]$\begin{array}{cc}{}_j\left(t\right)=\frac{q_j\left(t-1\right)+q_j\left(t\right)}{\underset{j=1}{\overset{m}{.Math.}}.Math.\left(q_j\left(t-1\right)+q_j\left(t\right)\right)},.Math.j=1,\mathrm{.Math.}.Math.,m,& \left(\mathrm{Equation}.Math..Math.32\right)\end{array}$

[0077] wherein .sub.j(t) is the independent contribution of the jth hidden neuron at time t; q.sub.j(1) is independent contribution output of the jth hidden neuron at time t1. q.sub.j(t) is independent contribution output of the jth hidden neuron at time t; Moreover, q.sub.j[q.sub.j(t1), q.sub.j(t)] is independent contribution output vector of the jth hidden neuron; Q(t)=[q.sub.1(t), . . . q.sub.m1(t), q.sub.m(t)].sup.T is the independent contribution matrix at time t,

Q(t)=(t)(t), (Equation 33)

[0078] wherein (t) is a coefficients matrix which is provided as:

(t)=D.sup.1(t)(t)B(t)z(t), (Equation 34)

[0079] wherein (t)=[(t1), (t)] is output matrix of hidden layer at time t, (t1)=[.sub.1(t1), .sub.2(t1), . . . , .sub.m(t1)].sup.T is output vector of hidden layer at time t1, (t)=[.sub.1(t), .sub.2(t), . . . , .sub.m(t)].sup.T is output vector of hidden layer at time t; D(t), B(t) and z(t) are the covariance matrix of (t), the whitening matrix of y(t) and the whitening transformation matrix of y(t), respectively. D(t), B(t) and z(t) are provided as:

[00012]$\begin{array}{cc}D\left(t\right)=\left[\begin{array}{cc}\frac{\underset{j=1}{\overset{m}{.Math.}}.Math.{\left({}_j\left(t-1\right)-\stackrel{\_}{}\left(t-1\right)\right)}^2}{m-1}& \frac{\underset{j=1}{\overset{m}{.Math.}}.Math.\left({}_j\left(t-1\right)-\stackrel{\_}{}\left(t-1\right)\right).Math.\left({}_j\left(t\right)-\stackrel{\_}{}\left(t\right)\right)}{m-1}\\ \frac{\underset{j=1}{\overset{m}{.Math.}}.Math.\left({}_j\left(t-1\right)-\stackrel{\_}{}\left(t-1\right)\right).Math.\left({}_j\left(t\right)-\stackrel{\_}{}\left(t\right)\right)}{m-1}& \frac{\underset{j=1}{\overset{m}{.Math.}}.Math.{\left({}_j\left(t\right)-\stackrel{\_}{}\left(t\right)\right)}^2}{m-1}\end{array}\right],& \left(\mathrm{Equation}.Math..Math.35\right)\\ B\left(t\right)={}^{-1/2}\left(t\right).Math.U^T\left(t\right),& \left(\mathrm{Equation}.Math..Math.36\right)\\ z\left(t\right)={}^{-1/2}\left(t\right).Math.U^T\left(t\right).Math.y\left(t\right),& \left(\mathrm{Equation}.Math..Math.37\right)\end{array}$

[0080] wherein (t1)=(0.sub.1(t1)+0.sub.2(t1)+. . . +0.sub.m(t1))/m is the average value of elements of output vector in hidden layer at time t1, (t)=(.sub.1(t)+.sub.2(t)+ . . . +.sub.m(t))/m is the average value of elements of output vector in hidden layer at time t; (t) and (t) are the eigenvector and eigenvalue matrices of y(t); y(t) is the output matrix of self-organizing recurrent RBF neural network at time t,

y(t)=(t)(t), (Equation 38)

[0081] wherein (t) is the weight matrix of hidden layer to output layer,

(t)=[w(t1),w(t)](Equation 39)

[0082] wherein w(t1)=[w.sub.1(t1) , w.sub.2(t1), . . . , w.sub.m(t1)].sup.T and w(t)=[w.sub.1(t) , w.sub.2(t), . . . , w.sub.m(t)].sup.T are the output of self-organizing recurrent RBF neural network, the output vector of the hidden layer and the weight vector at time t1 and time t, respectively.

[0083] {circle around (5)} t>3, calculating activity degree of hidden neuron:

S.sub.j(t)=e.sup.h.sub.j.sup.(t)c.sub.j.sup.(t), (Equation 40)

[0084] wherein S.sub.j(t) is the activity degree of the jth hidden neuron at time t, j=1, 2, . . . m.

[0085] {circle around (6)} t>3, adjusting the structure of the self-organizing recurrent RBF neural network:

[0086] In the process of adjusting structure of neural network, calculating the activity degree of the lth hidden neuron S.sub./(t) and the independent contribution of the /th hidden neuron .sub./)/(t).

[0087] When the activity degree and independent contribution of the lth hidden neuron satisfy:

S.sub.l(t)=max S(t), (Equation 41)

.sub.l(t)=max (t), (Equation 42)

[0088] wherein S(t)=[S.sub.1(t), . . . , S.sub.m1(t), S.sub.m(t)] is the vector of activity degree of hidden neurons at time t, (t)=[.sub.1(t), . . . , .sub.m1(t), .sub.m(t)] is the vector of independent contribution of hidden neurons at time t; adding one hidden neuron and the number of hidden neurons is M.sub.1=m+1; Otherwise, the structure of self-organizing recurrent RBF neural network is not adjusted, M.sub.1=m;

[0089] When the activity degree and independent contribution of the ith hidden neuron satisfy:

S.sub.i(t)=min S(t), (Equation 43)

.sub.i(t)=min (t), (Equation 44)

[0090] deleting the ith hidden neuron and updating the number of hidden neurons M.sub.2=M.sub.11; otherwise the structure of self-organizing recurrent RBF neural network is not adjusted, M.sub.2=M.sub.1;

[0091] {circle around (7)} increasing one learning steps, if s<N, then turning to step {circle around (3)}; if s=N, turning to step {circle around (8)}.

[0092] {circle around (8)} according to Eq. (5), calculating the performance of self-organizing recurrent RBF neural network. If E(t)E.sub.d, then turning to step {circle around (3)}; if E(t)<E.sub.d, stopping the training process.

[0093] The training result of the intelligent detection method for BOD concentration is shown in FIG. 2. X axis shows the number of samples. Y axis shows the BOD concentration. The unit of Y axis is mg/L. The solid line presents the real values of BOD concentrations. The dotted line shows the outputs of intelligent detection method in the training process. The errors between the real values and the outputs of intelligent detection method in the training process are shown in FIG. 3. X axis shows the number of samples. Y axis shows the training error. The unit of Y axis is mg/L.

[0094] (4) BOD concentration prediction;

[0095] The testing samples used as the input of self-organizing recurrent RBF neural network, the output of neural network is the computing values of BOD concentration. The predicting result is shown in FIG. 4. X axis shows the number of testing samples. Y axis shows the BOD concentration. The unit of Y axis is mg/L. The solid line presents the real values of BOD concentration. The dotted line shows the outputs of intelligent detection method in the testing process. The errors between the real values and the outputs of intelligent detection method in the testing process are shown in FIG. 5. X axis shows the number of samples. Y axis shows the testing error. The unit of Y axis is mg/L.

CONCLUSION

[0096] Although the subject matter has been described in language specific to structural features and/or methodological acts, it is to be understood that the subject matter defined in the appended claims is not necessarily limited to the specific features or acts described above. Rather, the specific features and acts are disclosed as example forms of implementing the claims.