Hindawi Publishing Corporation Mathematical Problems in Engineering Volume 2013, Article ID 679039,10 pages http://dx.doi.org/10.1155/2013/679039

Research Article

Robust Adaptive Control for Nonlinear Discrete-Time Systems by Using Multiple Models

Xiao-Li Li, De-Xin Liu, Jiang-Yun Li, and Da-Wei Ding

School of Automation and Electrical Engineering and the Key Laboratory of Advanced Control of Iron and Steel Process (Ministry of Education), University of Science and Technology Beijing, Beijing 100083, China

Correspondence should be addressed to Xiao-Li Li; lixiaoli@hotmail.com

Received 23 August 2013; Accepted 28 September 2013

Academic Editor: Zhiguang Feng

Copyright © 2013 Xiao-Li Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Back propagation (BP) neural network is used to approximate the dynamic character of nonlinear discrete-time system. Considering the unmodeling dynamics of the system, the weights of neural network are updated by using a dead-zone algorithm and a robust adaptive controller based on the BP neural network is proposed. For the situation that jumping change parameters exist, multiple neural networks with multiple weights are built to cover the uncertainty of parameters, and multiple controllers based on these models are set up. At every sample time, a performance index function based on the identification error will be used to choose the optimal model and the corresponding controller. Different kinds of combinations of fixed model and adaptive model will be used for robust multiple models adaptive control (MMAC). The proof of stability and convergence of MMAC are given, and the significant efficacy of the proposed methods is tested by simulation.

1. Introduction

Due to the strong ability of approximation, neural network has been widely used in the identification of nonlinear system. It is also a very useful tool for prediction, pattern recognition, and control [1]. The network structure comprises the interconnected group of nodes and the weight. There are many kinds of neural networks such as back propagation (BP), radial basis function (RBF), cerebellar model articulation controller (CMAC). As the most effective learning algorithm for feedforward networks [2], BP neural network has been the research focus for manyyears [3-6].

Adaptive control of nonlinear systems using neural network has been an active research area for over two decades [7-9]. The controller will be set up by adjusting the weights of the neural network [10,11]. But adaptive control using neural network still has the same shortcomings as conventional adaptive control; it is extensively studied in time-invariant system with unknown parameters or time-variant system with slow drifting parameters [12,13]. While the system has abrupt changes in parameters, the algorithm cannot find the exact identification model and will respond slowly to system parameter variations. To solve this kind of problem, MMAC has been a very useful tool in recent years.

Since MMAC was presented in 1970s, it has attracted a lot of attention of experts [14-17]. MMAC is an effective approach to solve problems such as time variations and uncertainties. It has the ability to improve the transient responses and the control performance. According to the dynamic character of controlled plant, multiple models are set up to cover the uncertainty of parameter. Much research has been done on continuous-time and discrete-time linear systems [18, 19]. For nonlinear system, only a few results have been given. In recent years, the MMAC based on neural network has been considered by some researchers [20, 21]. But in these papers, the nonlinear system has been modeled by the combination of linear model (the main part) and neural network model (the unmodeled dynamics). The multiple models are still multiple linear models with different parameters, and neural network is used only to compensate for the modeling error of linear model. In this case, the nonlinear system should not be very complex, and too big modeling error between the system and linear model is forbidden. The parameter and structure uncertainty of a relatively complex nonlinear system cannot be modeled by this method. This kind of MMAC with neural network still follows the main ideas of linear MMAC.

In this paper, a kind of robust MMAC is proposed for nonlinear system. Multiple BP neural networks with different weights will be used to cover the uncertainty of the parameters of the system. A performance index function based on the identification errors will be used to choose the best model and the corresponding controller. Considering the unmodeling error of neural network, a dead-zone recursive algorithm will be used, and the proof of robust property and stability of the MMAC are given. Different combinations of adaptive models and fixed models will be used for MMAC, and the effectiveness of the proposed method has been tested in simulations.

2. Robust Adaptive Control Using Neural Network

The single-input/single-output nonlinear discrete-time system can be represented as follows:

y(k+l) = f0 (■) + g0 (■)u(k-d+1), (1)

where f0 and g0 are infinitely differentiable functions of

y(k-n + 1),..., y (k) ,u(k - d - m + 1),... ,u(k - d),

where y is the output, u is the input, m < n, d is the relative degree of the system, and g0 is bounded away from zero. The arguments of f0 and g0 are real variables.

Due to the existence of noncausal problem, normally state transformation should be made first [11], and a causal system as follows can be given:

y(k + d) = fd_1 [x (k)] + gd_1 [x (k)] u (k). (3)

Assumptions 1 and 2 in [11] about g0(x),..., gd_1(x) and minimum phase assumption should still be satisfied.

As Assumption 3 in [11], there exist the weights w, v of neural network; the functions fd_1 [x(k), w] and gd_1 [x(k), v] can approximate the functions fd_1 and gd_1 with any accuracy e.

Plant (3) can be modeled by the neural network.

Consider

y(k + d) = fd_1 [x (k), w] + gd_1 [x (k), v] u (k). (4)

The functions fd_1 [■, ■] and gd_1 [■, ■] depend on the structure of the neural network and the number of neurons. For example, if fd-1[^r] and gd-1 [■, ■] are three-layer neural networks with p and q hidden neurons, respectively, then they can be expressed as

p j m+n+d_1 \

fd_1 [x (k),w] = £wtH ( T wtjxi(k) +B¡J ■,

q / m+n+d_1 \

gd_1 [x (k),v] = ^¡H ( T Vtjxi (k) +VtJ' where H is a hyperbolic tangent function.

Let w(k) and v(k) denote the estimates of w and v at time k. Rewrite (3) and (4) as follows:

y(k+\) = fd_1 [x (k-d+1)]

+ dd_1 [x (k-d+l)]u(k-d+l), y(k+\) = fd_1 [x (k-d+\), w]

+ gd_1 [x (k-d+1), v]u(k-d + 1). We have the estimated plant output as:

/ (k+1) = fd_1 [x(k-d+1),w(k)]

+ gd_1 [x(k-d+ 1),v(k)] u(k-d + 1). Define e*(k + 1) as

e* (k+1) = y* (k+1)- y(k+1).

If the neural network could approximate the nonlinear system with zero error; that is, e = 0, the following weight 0 = [ W ] updating rule can be used:

0 (k+1) = 0 (k)--^-e* (k+1) J (k-d+1), (9)

where J (k-d+1)

dy(k+1)*

dfd_1 [x(k-d+1), w]

ddd-i [x(k-d+1),v]

, w(k)

u(k - d + 1)

and r(k) is the reference command. One has

r(k) = 1+ J' (k-d+1) J (k-d+1). (11)

Due to the existence of unmodeling dynamics e > 0, the design of robust adaptive controller should be considered. A dead-zone algorithm will be used instead of (9) for updating the weights. Therefor,

D(e) =

0, ee +1

if \e\ < d0, if e > d0, if e < -dn,

where d0 will be the function of e which can be seen in proof procedure of the Theorem 1.

The output of the dead-zone function is used in the following updating rule:

0 (k+1) = 0 (k)--^-D(e* (k+1)) J (k-d+1). (13) r(k)

Define the parameter error as

0(k) = 0 (k)- 0 and give the control u(k) as follows:

-fd_1 [X (k), W (k)]+r(k) U( ) g^ [x (k), k (k)]

We have Theorem 1 for the feedback control system.

Theorem 1 (see [11]). Suppose \r(k)\ < d1 for all k > 0. Given any constant p > 0 and any small constant d0 > 0, there exist positive constants p1 = p1(p,d1), p2 = p2(p,d1), e* = e*(p,d0,d1), and S* = S*(p,d0,d1) such that if Assumptions 1 and 3 are satisfied on S D Bp1, with e < e*, Assumption 2 is

satisfied on Bp2, \x(0)\ < p, and \0(O)\ < S < S*, then

(1) the x(k) and u(k) are bounded for all k,

(2) \0(fc)\ will be monotonically nonincreasing, and 0(k + 1) - 0(k) will converge to zero,

(3) XZ1([D(e*(k+1))]2/(1+ j'(k-d+1)j(k-d+1))) < 0,

(4) the tracking error between the plant output and the reference command will converge to a ball of radius d0 centered at the origin.

3. Robust Multiple Model Adaptive Control

The conventional adaptive control systems are usually based on a fixed or slowly adaptive model. It cannot react quickly to abrupt changes and will result in large transient errors before convergence. For this kind of problem, MMAC algorithm is presented as a useful tool. The rationale for using MMAC is to ensure that there is at least one model with parameters sufficiently close to those of the unknown plant. By the switching rule, the control strategy is to determine the best model for the current environment at every instant and activate the corresponding controller. The structure of the multiple model adaptive control is shown in Figure 1.

3.1. Architecture of the Control System. Multiple adaptive models can be regarded as an extension of conventional indirect adaptive control. The objective is to make the control error ec = r-y tend to zero, where r is the desired output. The control system contains N identification models, denoted by 7®, I e {1,2,..., N] according to (7), operating in parallel. Consider the following:

{l) : (k+1) = [x (k-d+ 1), w® (k)]

+ gf^ [x (k-d+1), v® (k)] (16) xum (k-d+1),

ff_, [x (k), w® (k)]

W? (k) HI

m+n+d_\

W® (k)X: -ë? (k)

Figure 1: Structure of multiple model adaptive control.

^ [x (k), v® (k)]

q I m+n+d_l \

= YV? (k)H( 1 V® (k)Xj - vf (k)\.

The identification error between the output y of I® and y of the plant is denoted as e= y*^ - y. Corresponding to each I® is a parameterized controller C® which achieves the control objective for 7®. The output of C® is denoted by w®.

One has

-fd»! [x (k), W® (k)]+T(k)

[x (k), V® (k)]

where I e {1,2,..., N], and r(k) is the reference command.

At every instant, one of the models I® is selected by a switching rule, and the corresponding control input is used to control the plant.

Given prior knowledge of the different possible environments, the control problem is to determine suitable rules for switching and tuning these parameters to yield the best performance for the given objective while assuring stability.

3.2. Choice of Multiple Models. The following three different combinations have been considered [22].

3.2.1. N Adaptive Models. N adaptive models with different initial parameter values can be viewed as an extension of conventional indirect adaptive control. When the parameters of the plant change abruptly, the change can be detected by identification error. Then, the parameters of models are reset to initial values, and the model with the smallest error is selected. Hence, we can construct multiple adaptive models with different initial parameters which adjust dynamically in any instant. The method was considered in detail in [8,19].

3.2.2. N-1 Fixed Models and One Adaptive Model. The previous method reveals that massive calculation may be produced because each adaptive model needs to adjust dynamically. Hence, if the models are fixed, the same strategy can be used in stationary and time-varying environments. However, fixed models can represent exactly only a finite number of environments. Thus, N-1 parallel fixed model and one adaptive model are combined, and the efficiency can be improved by the multiple fixed models and the accuracy can be increased by the adaptive model.

Theorem 2. Suppose \r(k)\ < d1 for all k > 0. Given any constant p > 0, for all the model (adaptive model or fixed

model), if \%(0)\ < p, \0(°(0)\ < 5 < 5*, I = 1,2,..., N, and the conditions in Theorem 1 are satisfied, then when index switching (19) is used,

(1) all the signals in the system are bounded,

(2) the tracking error between r(k) and y(k + d)

lim \r(k) - у (k + d)\ <

3.2.3. N-2 Fixed Models, One Free Running Adaptive Model, and One Reinitialized Adaptive Model. It is commonly accepted that the convergence time of an adaptive model will be large for large initial parametric errors. Hence, in the configuration described above, a large number of fixed models may be needed to keep the transient response under control until the adaptive model has converged. If the fixed model, which is the closest to the given plant, is assumed to be known, faster convergence can be obtained by initiating a newadaptivemodel from thelocationofthe former. Thesame objective can be achieved on-line by starting adaptation from the location of each different fixed model that is successively chosen by the switching scheme.

The reinitialized adaptive model Ir included is introduced, and its parameters are determined as follows: if a fixed model I{1), I e {1,2,..., N - 2}, is activated by the switching rule at any instant k, then the parameters of Ir are reinitialized to the value of I®. Thereafter, this adaptive model will be left to adapt until the next reinitialization.

3.3. Choice of the Switching Rule. A natural way to decide when and to which controller one should switch is to determine performance cost indexes for each controller C®, I e {1,2,..., N}, and switch to the one with the minimum index at every instant. However, since only one control input can be used at any instant, the performance of any candidate controller can be evaluated only after it has been used. On the other hand, the performance of all the identification models can be evaluated in parallel at every instant. Hence, the indexes must be based on the performance of the models rather than the controllers, that is, using identification errors e® rather than the control error ec. From an adaptive control point of view, this rationale extends the principle of certainty equivalence from tuning to switching.

Considering the unmodeling error of neural network and robustness of the adaptive controller, the specific performance index proposed has the form

[D(e* (k+1))]2

+ J' (k-d+1) J (k-d+1)'

where e*(k+1) = y* (k + 1) - y(k + 1).

The switching scheme consists of monitoring the performance indexes }m(k) at every instant. After every switching, the controller corresponding to the model with the minimum index is chosen (switched) to control the plant.

3.4. Proof of the Multiple Models Stability

3.4.1. N Adaptive Model. At time k, Ith adaptive model will be selected, I e {1,2,..., N}. The control input

r (k) = fd_! [x (k), w® (k)] + gd_x [x (k), v® (k)] u (k).

From Theorem 1, if |x(0)| < p, |©®(0)| < S < S* ,withcontrol input (21), we can have that x(k) and u(k) are bounded, for all k, and

lim |e*® (k)\= lim \y*w (к) - у (k)\ <

lim (k)\ = lim 10® (k+1)- ©w (k)\ 0.

k —> <ю I \ k —ю ' '

Then, the control error at time к is given by

|/® (k + d)-r (k)\ < If,-! [x (k), W® (k + d- 1)]

[x (k), W® (k)]\ + \{gd_! [x (к), V® (k + d-1)]

-0d-\ [x (к), V® (k)]}u(k)\

< к |0® (k-d+1)- 0® (k)\

0 as к

\r (k)-y(k + d)\= \r (к) - /® (k + d)

+ /® (k + d)-y(k + d)\

< \r(k)-f ® (k + d)

+ 1/® (k + d)-y(k + d)\

< \y*(l) (k + d)-y(k + d)\

< d0 as k —> >x>.

Let one of N controller be chosen at random, for any instant of time k and any model chosen; then (24) holds.

So we have all the signals in the system bounded, and

(1) lirnfc^m\e*m(k)\ < d0,

(2) lirnfc^m\r(k) - y(k + d)\ <dQ.

3.4.2. N-1 Fixed Model and One Adaptive Model. Consider the following index function:

% (k) = 1:

K (e* (k+1))f

k=i 1 + (J® (k-d+ 1))J® (k-d+1)'

At every instant k, model 7®,

I = arg min Jm (k),

1<1<N m

will be selected.

Proof. From Theorem 1, we have that, for adaptive model,

lim im (k)

= lim 1

[D(a) (e* (k+1))]'

1 + (J(a) (k-d+ 1))'j(a) (k-d + 1)

For the fixed models I, lim^ TO/®(fc) is either bounded or <<, if the performance index J®(k) tends to <; there exist a time t, J^ (k) < ¡m(k), k>t, which implies that the adaptive model will be selected finally.

If limk^mjm)(k) is bounded, then

[D(l\e*(k+1))]2

lim 1— .

1 + (J® (k-d+ 1))'J® (k-d+1)

[D(l\e*(k+1))f

lim ----

1 + (J© (k-d+ 1))'J® (k-d+1)

< œ, (28)

If the fixed model is selected, the relationship of p, e, S in Theorem 1 is satisfied; the proof procedure will be similar to that of multiple adaptive controller. We also have the following:

(1) all the signals in the system are bounded,

(2) lirnfc^m\e*(l\k)\ < dQ,

(3) limk^œ\r(k)-y(k + d)\<d0.

3.4.3. N-2 Fixed Model, One Free Running Adaptive Model, and One Reinitialized Adaptive Model. The introduction of the reinitialized adaptive model will not affect the stability of the whole system, and the proof of the stability will be similar to the case of N-1 fixed model and one adaptive model.

4. Simulation

4.1. The Problem. PH neutralization is a very important procedure in the chemical industry. Usually, we use the logarithmic behavior to present pH characteristic; the existing nonlinearity always makes the identification and control of pH neutralization more difficult. A strong acid flows into a tank and is thoroughly mixed with a strong base whose inward rate of flow is controlled in such a way to produce a neutral outward flow from the tank. Because the acid and the base are strong, they are completely dissociated, and also the dissociation of the water can be disregarded [23]. The equation describing this model is

V% = F{k) (a-y {k)) -u(k)(h + y W)> (3°)

where y(k) = [H+] - [pH-] is the distance from neutrality,

V = volume of the tank,

F(k) = rate of flow of the acid,

a = concentration of the acid,

u(k) = rate of flow of the base,

b = concentration of the base.

Note that pH value can be determined from the y(k), pH(fc) by the following nonlinear transformation:

pH (k) = -logi

y(k) + ^(y(k))2 + 4KW

where Kw = water equilibrium constant - 1014.

We suppose that b and a are fixed and known, that F(k) can be measured online, and that u(k) can be given assigned values within certain limits.

An approximate discrete-time model can be developed, incorporating measurement and input actuator errors, as follows:

yk+1 -yk +y[F (k) {a - yk) -uk {b + yk)]. (32)

The following values were adopted for the various quantities of interest:

0.1 <F(k)< 0.125 l/min, 0<u(k)< 0.2 l/min,

a=10_3 mol/l, V=21, (33)

b=10_3 mol/l, T=2 min. Then our goal is to control the plant as follows:

y(k+1) = y(k)[1- 0.5F (k)] + 0.5 x 10_3F (k)

- 0.5 [10_3 + y(k)]u(k),

where F(x) is the variable parameter. And the single-input/ single-output nonlinear discrete-time system represented as (1) can be modeled by

/ (k+1) = f[y(k), w (k)] +g[y(k), k (k)] u (k), (35)

where f and g are the output of neural network.

IW- ........

50 100

200 250

u 0.06

50 100

(a) System output pH(k) (b) Control input u(k)

Figure 2: One adaptive mode with fixed parameter.

200 250 300

The neural networks f and g are 3 layered with 4 neurons in each hidden layer. Based on the error between the plant output y(k + 1) and the model output y* (k+ 1), the network parameters w(k) and v(k) are updated to w(k + 1) and v(k + 1) using the standard back propagation algorithm.

4.2. Experiments

Adaptive Model Only. In the process of parameter identification, let F(k) = 0.1, d0 = 0.001, y(1) = 0, y(2) = 0, the reference command r(k) = 0 (according to (31), pH(fc) = 7), and k = 1,2,..., 300. The initial weights w and v are given a random number in the range [-1,1]. The identification results are shown in Figure 2.

After 300 sample times, the weights will converge to the following values:

w = [0.865,0.799, -0.823,1.478,0.232, 0.484,

- 0.251, 0.794,1.052, -0.742, 0.591, -0.019]T, v = [-0.483, -0.665, -0.288, -0.767, -0.519, -0.187, 0.506, -0.236, -0.268, 0.965,1.201, 0.891]T

Now, we consider a worse case of the plant where abruptly changing parameters appear:

F(k) =

0.1, 1<k<100, 0.113, 100 < k < 200, 0.125, 200 <k< 300.

Given the initial weights as the convergent weights of w and v in (36).

As the parameters change at k = 100 and k = 200, the overshoot of the system is big and the settling time is long. The nonlinear system cannot track the reference trajectory in time (Figure 3). When MMAC is used, the following simulation results can be obtained.

4.2.1. Three Adaptive Models. Three adaptive models I(-al\ I(a2), and I("3 are established. According to the three different values of F(k) obtained (k = 1,100,300), each group of weights can be got using the same method as that in (36). One has

( a1) ( a1)

w = w, v = w

w(a2) = [-0.423, 0.482, -0.054,1.023,0.487, 0.248,

0.072,0.825, 0.230, -0.214, -0.678,0.032]T,

v(a2) = [-1.064,-0.103,0.563,0.578,-0.439,-0.487,

0.258, -0.557, -0.820, 0.385,0.842, 0.831]T,

w(a3) = [-0.584,0.341,0.376,0.733,0.394,0.476,

1.067,1.089, -0.144, -0.577,0.051, -0.058]T,

v(a3) = [0.811,0.076,-0.532,0.399,0.412,-0.268,

0.897, -0.384,0.933, -0.131, -0.313, -0.758]T

The multiple models based on neural networks are chosen as in (16). Figures 4(a) and 4(b) present the responses of the plant. Switching sequence of controllers is shown in Figure 4(c). Obviously, this method can track the reference trajectory fast and improve the transient response. According to the index function, the system can choose an approximate model to identify the unknown plant. Once the parameters change, the weights and the index functions of neural networkmodels will be initialized and the system will choose the optimal model again to conduct identification. In this way, the overshoot of the system can be decreased and the reference trajectory can be tracked fast at the same time.

4.2.2. Three Fixed Models and One Adaptive Model. In this case, three fixed models I^1, I(f2\ and I^3 are used to

(a) System output pH(k)

0.14 0.12 0.1 0.08 0.06 0.04 0.02

200 250 300

50 100

(b) Control input u(k)

Figure 3: One adaptive mode with variable parameter.

250 300

If" V

(a) System output pH( k) 4

8 2.5 o

200 250 300

0.14 0.12 0.1 0.08 0.06 0.04 0.02 0

(b) Control input u( k)

200 250 300

(c) Switching scheme Figure 4: Three adaptive models.

200 250 300

^ 2 C z

(a) System output pH(fc) 5

(b) Control input u( k)

4.5 4 3.5 3 2.5 2 1.5 1

(c) Switching scheme Figure 5: Three fixed models and one adaptive model.

improve the transient response, and one adaptive model I(a4) is used to guarantee the stability. The initial weights of the three fixed models are the same as those of the three adaptive models, and the weights of one adaptive model I(a4) are equal to those of I(f1).

In the process of parameter identification, this method could improve the transient response compared with the conventional adaptive control (Figures 5(a) and 5(b)). Switching sequence of controllers is shown in Figure 5(c). Once the parameter changes abruptly at 100 or 200, the controller will switch to the nearest fixed model I(f2\ I(f3) to reduce the error. When the adaptive model gradually converges to the true value, the system will switch to the adaptive model I(a4\ Multiple fixed models play a transitional role in the process of identification. This method can reduce massive calculation compared with the case of three adaptive models, but it produces a larger overshoot compared with Figures 4(a) and 4(b).

4.2.3. Three Fixed Models, One Free Running Adaptive Model, and One Reinitialized Adaptive Model. In this case, we

establish three fixed models I(^1'), I(f2\ and I^33 with different initial weights; I(a4> is the free adaptive model and I(55 is the reinitialized adaptive model. The reinitialized adaptive model can achieve the initial weights by choosing a set of fixed models based on the past performance of the plant. If at any instant one of them is determined to be the best, the reinitialized adaptive model can be adapted from this model.

From the simulation, we can see that this method can improve the control quality dramatically (Figures 6(a) and 4(b)). Switching sequence of controllers is shown in Figure 6(c). Compared with the other algorithm proposed before, this method show perfect performance in reducing the overshoot and tracking the reference trajectory, and computation time is reduced greatly.

5. Conclusion

In this paper, multiple models are used to establish robust multiple models adaptive controller for a class of nonlinear discrete-time systems by using neural networks. Three kinds of combinations of adaptive model and fixed model are used

K 7 di

50 100 150 200

(a) System output pH( k)

250 300

(b) Control input u( k)

50 100 150 200 250 300

(c) Switching scheme

Figure 6: Three fixed models, one free running adaptive model and one reinitialized adaptive model.

to make the multiple model set, and a switching law is suitably defined to make the decision of the best model. The principal contribution of this paper is the proof of stability of robust MMAC by using neural networks. Multiple neural network models with different weights represent different dynamical characters ofthe plant when it operates in different environments, which can be described by a mount of input and output data. So the design of the model set can also be regarded as a kind of data driven problem [24, 25]. How to divide the region of data into suitable numbers of subregions which can be represented by multiple neural network models will decide the accuracy of MMAC. A moving or dynamically optimal model set will be an important problem that needs to be solved in the future.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was supported by the Fundamental Research Funds for the Central Universities under Grant FRF-TP-12-005B, the Program for New Century Excellent Talents in Universities under Grant NCET-11-0578, and the National Natural Science Foundation of China under Grant 61074055.

References

[1] D. Wang and J. Huang, "Adaptive neural network control for a class of uncertain nonlinear systems in pure-feedback form," Automatica, vol. 38, no. 8, pp. 1365-1372, 2002.

[2] D. E. Rumelhart, G. E. Hinton, and R. J. Williams, "Learning representations by back-propagating errors," Nature, vol. 323, no. 6088, pp. 533-536,1986.

[3] X. J. Jing and L. Cheng, "An optimal PID control algorithm for training feedforward neural networks," IEEE Transactions on Industrial Electronics, vol. 60, no. 6, pp. 2273-2283, 2013.

^ 3 c 3

[4] B. William and L. Stanislao, "Robot navigation control based on monocular images: an image processing algorithm for obstacle avoidance decisions," Mathematical Problems in Engineering, vol. 2012, Article ID 240476,14 pages, 2012.

[5] M. Arvaneh, C. Guan, K. K. Ang, and C. Quek, "Optimizing spatial filters by minimizing within-class dissimilarities in electroencephalogram-based brain—computer interface," IEEE Transactions on Neural Networks and Learning Systems, vol. 24, no. 4, pp. 610-619, 2013.

[6] H. A. Talebi and K. Khorasani, "A neural network-based multiplicative actuator fault detection and isolation of nonlinear systems," IEEE Transactions on Control Systems Technology, vol. 21, no. 3, pp. 842-851, 2012.

[7] K. J. Astrom and B. Wittenmark, Adaptive Control, Addison-Wesley, Longman Publishing, Boston, Mass, USA, 1995.

[8] K. S. Narendra and A. Annaswamy, Stable Adaptive Systems, Prentice-Hall, Englewood Cliffs, NJ, USA, 1989.

[9] P. A. Ioannou and J. Sun, Robust Adaptive Control, Prentice-Hall, Upper Saddle River, NJ, USA, 1996.

[10] K. S. Narendra and K. Parthasarathy, "Identification and control of dynamical systems using neural networks," IEEE Transactions on Neural Networks, vol. 1, no. 1, pp. 4-27,1990.

[11] F. C. Chen and H. K. Khalil, "Adaptive control of a class of nonlinear discrete-time systems using neural networks," IEEE Transactions on Automatic Control, vol. 40, no. 5, pp. 791-801, 1995.

[12] F. C. Chen and H. K. Khalil, "Adaptive control of nonlinear systems using neural networks," International Journal of Control, vol. 55, no. 6, pp. 1299-1317, 1992.

[13] G. Kreisselmeier and B. D. O. Anderson, "Robust model reference adaptive control," IEEE Transactions on Automatic Control, vol. 31, no. 2, pp. 127-133,1986.

[14] L. Giovanini, "Robust adaptive control using multiple models, switching and tuning," IET Control Theory & Applications, vol. 5, no. 18, pp. 2168-2178, 2011.

[15] W. Wang and X. L. Li, Multiple Model Adaptive Control, Science Press, Beijing, China, 2001.

[16] D. T. Magill, "Optimal adaptive estimation of sampled stochastic processes," IEEE Transactions on Automatic Control, vol. 10, pp. 434-439,1965.

[17] D. G. Lainiotis, "Partitioning: a unifying framework for adaptive systems. I. Estimation," vol. 64, no. 8, pp. 1126-1143,1976.

[18] M. Athans, D. Castanon, K.-P. Dunn et al., "The stochastic control of the F-8C aircraft using a multiple model adaptive control (MMAC) method—part I: equilibrium flight," IEEE Transactions on Automatic Control, vol. 22, no. 5, pp. 768-780, 1977.

[19] K. S. Narendra and J. Balakrishnan, "Adaptive control using multiple models," IEEE Transactions on Automatic Control, vol.

42, no. 2, pp. 171-187,1997.

[20] L. Chen and K. S. Narendra, "Nonlinear adaptive control using neural networks and multiple models," Automatica, vol. 37, no. 8, pp. 1245-1255, 2001.

[21] Y. Fu and T. Chai, "Nonlinear multivariable adaptive control using multiple models and neural networks," Automatica, vol.

43, no. 6, pp. 1101-1110, 2007.

[22] K. S. Narendra and C. Xiang, "Adaptive control of discrete-time systems using multiple models," IEEE Transactions on Automatic Control, vol. 45, no. 9, pp. 1669-1686, 2000.

[23] G. C. Goodwin and K. S. Sin, Adaptive Filtering Prediction and Control, Prentice-Hall, Englewood Cliffs, NJ, USA, 1984.

[24] S. Yin, S. X. Ding, A. Haghani, H. Y. Hao, and P. Zhang, "A comparison study of basic data-driven fault diagnosis and process monitoring methods on the benchmark Tennessee Eastman process," Journal of Process Control, vol. 22, no. 9, pp. 1567-1581, 2012.

[25] S. Yin, H. Luo, and S. X. Ding, "Real-time implementation of fault-tolerant control systems with performance optimization," IEEE Transactions on Industrial Electronics, no. 99, p. 1, 2013.

Copyright of Mathematical Problems in Engineering is the property of Hindawi Publishing Corporation and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use.