Abstract

This paper deals with the self-tuning regulator for large-scale stochastic nonlinear systems, which are composed of several interconnected nonlinear monovariable subsystems. Each interconnected subsystem is described by discrete Hammerstein model with unknown and time-varying parameters. This self-tuning control is developed on the basis of the minimum variance approach and is combined by a recursive algorithm in the estimation step. The parametric estimation step is performed on the basis of the prediction error method and the least-squares techniques. Simulation results of the proposed self-tuning regulator for two interconnected nonlinear hydraulic systems show the reliability and effectiveness of the developed method.

1. Introduction

In the past few years, large-scale systems, which can be efficiently applied to many practical situations, such as power systems, transportation systems, industrial processes, and communication networks, have attracted much attention [125]. There exists extensive literature concerning control for large-scale systems. In this case, several results deal with some schemes for control of these systems on the basis of various approaches, such as minimum variance approach [7], generalized minimum variance strategy [8], decentralization [1, 2, 12], sliding mode control technique [23], artificial neuronal networks, and fuzzy logical techniques [911].

Self-tuning regulators of stochastic systems with time-varying parameters are a problem of theoretical and practical importance. Since the synthesis of the first self-tuning regulator by Åström and Wittenmark in 1973, great numbers of theoretical and practical results are published in the literature in order to adjust automatically the regulator parameters online in response to change in the process and the environment. These results concerned more particularly the linear stochastic systems [7, 8, 2628], the dual-rate sampled-data systems [2931], and so forth, when various types of adaptive controllers are developed. In the field of adaptive control, Kamoun presented optimal self-tuning regulators for large-scale linear systems, which are described by discrete input-output mathematical models, based on minimum variance and generalized minimum variance strategies [7, 8]; Zhang et al. developed a self-tuning control scheme based on multi-innovation stochastic gradient parameter estimation for discrete-time systems [28]; Ding et al. proposed an adaptive control algorithm for Hammerstein nonlinear dual-rate systems [29]; and Ding and Chen developed least-squares and gradient based adaptive control algorithms for dual-rate linear systems [30, 31].

The minimum variance regulator was proposed by Åström and Wittenmark [26], which allows solving the regulation problem for a small size dynamic system described by an input-output mathematical model with known parameters. Afterwards, this approach has been extended by Kamoun for the control of large-scale linear systems, which are composed of several Single-Input Single-Output (SISO) interconnected subsystems and described by linear input-output mathematical models [7, 8]. In the same way, Elloumi and Kamoun have considered the class of large-scale nonlinear systems, which are comprised of several SISO interconnected nonlinear subsystems, and developed an implicit self-tuning regulator for this class of dynamic systems [21].

Our objective is to develop a self-tuning regulator with explicit scheme which can be applied to the class of large-scale nonlinear systems with unknown time-varying parameters in order to resolve the regulation problem for this class of systems. We particularly focus on the dynamic large-scale nonlinear systems which are constituted by several SISO interconnected nonlinear systems and described by stochastic Hammerstein mathematical models with unknown time-varying parameters. The developed algorithm is combined by a recursive parametric estimation algorithm for determining the system parameters.

This paper is organized into four sections. The second section deals with the parametric estimation problem of the large-scale nonlinear systems operating in a stochastic environment. We focus here on the dynamic large-scale systems consisting of several interconnected nonlinear monovariable subsystems, which can be described by the class of Hammerstein mathematical models, with unknown and time-varying parameters. Section 3, which is the principal part of the paper, detailed the self-tuning control problem of the nonlinear interconnected subsystems based on the minimum variance approach. This control problem is formulated while being based on a control strategy minimizing a certain quadratic criterion. Finally, a simulation example is provided to illustrate the effectiveness and the efficiency of the proposed self-tuning regulator scheme in Section 4. Conclusions are drawn in Section 5.

2. Description and Parametric Estimation

We consider a large-scale nonlinear system composed of SISO interconnected nonlinear subsystems . We assume that these nonlinear subsystems operate in a stochastic environment with unknown time-varying parameters. Thereby, the general structure of this considered system can be illustrated by Figure 1.

The dynamic linear part of Hammerstein structure, as shown in Figure 1, is described by the following equation:where and are, respectively, the input and the output of the dynamic linear part at the discrete-time , , , and represent the inputs of the nonlinear statics parts, and denote the inputs from the other interconnected dynamics systems , , , is an independent random variable and is generated as Gaussian distribution with zero average and constant variance , is an intrinsic delay of the interconnected system , and represent the delays of the interactions, which are related, respectively, to the inputs and the outputs of the other interconnected nonlinear systems , , , , and are polynomials with unknown and time-varying parameters, and is a polynomial with unknown but constant parameters, which are defined bywith , where , , , , and are the orders of the polynomials , , , , and , respectively.

The nonlinear statics parts of the Hammerstein structure are represented by the following nonlinear functions:where , , and are nonlinear functions.

In the literature, several presentations of these nonlinear functions are used like the polynomial functions that are frequently used in the description of these nonlinear elements. In this case, we propose to present these nonlinear functions by the following polynomials:where , , and represent approximated errors that can be assumed as noise acting on the system output. Their variances values depend on the chosen value of the nonlinearity degree .

In order to facilitate the formulation of the parametric estimation problem for large-scale nonlinear systems, we retain the following assumptions.

Assumption 1. The choice of the nonlinearity degree is made in such a way that the approximated errors , , and can be neglected.

Assumption 2. The polynomials , , , , and of the mathematical model of the linear dynamic part have the same order .

Taking into account these assumptions and using (1) and (4), we can express the system output as follows:with , .

The stochastic nonlinear system output , which is described by the Hammerstein mathematical model (5), can be rewritten in the following developed form:

We remark that model (6) contains some redundant parameters, such as the parameters , , and , , which are nonlinear over the parameters , , and , , , , , respectively. In that event, the formulation of the parametric estimation scheme presents some implementation difficulties. To overcome these difficulties, we proposed one of the following two situations [22]:(1)To have knowledge about one of the parameters , , and , , . For example, when the parameters , , and are well known at the discrete-time , we can estimate the parameters , , and , , , . Thereafter, we can easily determine the parameters , and , , based on the knowledge of , , and .(2)To have knowledge about one of the parameters , , and , , . For example, in the case where the parameters , , and are well known at the discrete-time , we can estimate the parameters , , and , . Afterwards, based on the knowledge of these estimated parameters, we can easily determine the parameters , , and .

To overcome this problem, we assume, in what follows, that , , and are constant and known, in such a way that , . In this case, the system output can be expressed byor equivalently in the matrix formwhere and represent, respectively, the parameter and the observations vectors, which are defined bywith

In the identification area of nonlinear systems, a large amount of research exploring different approaches has been published in the literature [22, 3234]. Ding et al. developed an auxiliary model-based least-squares identification algorithm for nonlinear Hammerstein systems [32]. Ding and Chen presented an iterative and a recursive least-squares identification method for Hammerstein nonlinear ARMAX systems [33]. Recently, Ding et al. studied the parameter estimation problem of a class of nonlinear systems and proposed a RLS algorithm based on the model decomposition [34]. For the class of large-scale nonlinear systems, Elloumi and Kamoun developed an iterative estimator for identification of Hammerstein large-scale systems, based on prediction error method and least-squares techniques [22]. In this case, the estimate of the parameter vector , as defined by (9), is ensured starting from the following recursive extended least-squares (RELS) algorithm [22]:where is the forgetting factor (), is the vector of the approximated observations , and is the estimated parameters vector, which are described by

The practical implementation of the recursive extended least-squares (RELS) algorithm requires its initialization. Thus, we must give the values of the initial conditions of the estimated parameters vector and the adaptation gain matrix . The convergence analysis of the RELS algorithm, which permits estimating the parameters of this class of systems, was proven by Elloumi and Kamoun [22].

3. Minimum Variance Self-Tuning Regulation

In this part, we propose to formulate the regulation problem of interconnected nonlinear systems, which can be described by Hammerstein stochastic mathematical model (1). In this case, we present an explicit self-tuning algorithm for the control of the subsystem . This algorithm is obtained by combining a least-squares estimator of the explicit model. The purpose of the developed controller is to regulate the state of the overall system and to maintain stability in spite of the presence of nonlinear interconnections, environment variations, and ignorance and time variations of the parameters.

To develop an optimal self-tuning regulator for this class of large-scale systems, we retain these assumptions.

Assumption 1. The delays , , and are well known. In addition to that, we suppose also that and to facilitate the implementation of the control law .

Assumption 2. The parameters , , and are well known at each discrete-time , in such a way that , .

Assumption 3. The sequences and are independent.

Let us consider a stochastic large-scale system, which consists of SISO interconnected nonlinear subsystems , with unknown time-varying parameters. This considered system can be described by the proposed Hammerstein model (1). Based on the cited assumptions, we can write the system output in the following form:with , , , and .

Or, equivalently, at the discrete-time ,

The problem arising here consists in the development of a self-tuning regulator with explicit scheme that can be applied to the considered stochastic system, which is described by the discrete model (15). The formulation of this control problem can be ensured starting from the minimization of the following quadratic criterion:where represents the expected value and is the reference output of the considered system .

By carrying out an Euclidian division of until the order , we can write (15) as follows:where the polynomials and are defined byThese two polynomials are solutions of the equation

Taking into account (14), the system output , as defined by (17), can be rewritten in the following form:where the polynomial is described as follows:with and .

We can show easily that the quadratic criteria , as given by (16), can be put in the following form:where the polynomials and are defined as follows:

The optimal controller making it possible to minimize the quadratic criterion, as given by (16), is obtained by solving the following equation:Equation (25) can be written in the following polynomial form:

It should be noted that (25) admits real and/or complex solutions for each interconnected subsystem , . In this instance, we regard only the real solutions, in order to have a feasible controller . Moreover, if we get more than a real solution, the choice of the optimal control law , to minimize the effects of the noise acting on each output of the interconnected nonlinear system, can be conducted based on the following relation:where is the real solutions obtained starting from solving (25) and represents the optimal controller for the considered system, allowing the system output to follow the desired constant reference signal with a minimum variance output.

The implementation of the optimal control law requires the knowledge of the different parameters involved in the polynomials , , , , and , . In this context, we need to add a parametric estimation step in the implementation.

In general, the minimum variance self-tuning regulator with explicit scheme can proceed by considering the following three steps.

Step 1. Estimate the different parameters involved in the Hammerstein mathematical model (1), based on the RELS algorithm.

Step 2. Determine the parameters of the polynomials and by resolution of the polynomial (19) and based on the estimated polynomials and . We can also calculate the parameters of , , as given by (21), from the knowledge of the estimated parameters values and the values of the parameters involved in the polynomials and . Similarly, from the knowledge of the estimated parameters values , , and the values of the parameters involved in the polynomials and , we can calculate the parameters of the polynomial , as defined by (23). In addition to that, we can compute the parameters of the polynomial , as described by (24), from the knowledge of the estimated parameters values , , and the values of the parameters involved in the polynomials and .

Step 3. Compute the control law starting from the resolution of polynomial (25) or (26).

It should be noted that the developed self-tuning regulator uses real-time estimation of the parameters of the interconnected nonlinear system , . In addition to that, the control law is obtained by solving (25) and based on relation (25), which must be made at each iteration step.

Let us note that the stability of the interconnected nonlinear system is ensured in closed loop if the controller when tends to infinity. This occurs to show the conditions which allow obtaining a bounded control signal. These requirements can be achieved on the basis of an extended version of Kharitonov’s method. This will be studied in a future research.

4. Simulation Results

The interconnected hydraulic system that we will deal with in this section is installed in petroleum society and depicted in Figure 2, in order to test the performance and the efficiency of the developed self-tuning regulator. This process is composed of two separators (V120 and V121), four motor-pumps shipments (P300, P305, P310, and P315), and four valves controls (LCV120, LCV121, LCV130, and LCV413). The docket of these separators is to deal with the fluid reservoir as it arrives from other separation stages, to segregate the three phases, water, oil, and gas, with stabilization of oil level in each one.

Figure 3 represents the interaction structure of the considered two interconnected subsystems and .

This considered system , , can be described by the following mathematical model:where the polynomials , , , and are defined byThe nonlinear statics parts of the considered system are represented by the following polynomials:

We will apply in this simulation study the control method, developed in Section 3, to each interconnected nonlinear system , which is described by the mathematical model (28). Thus, we can express the interconnected system output , , by the following developed forms:In this example of our numerical simulation, the noise sequences and are assumed to be independent and corresponding to a Gaussian distribution with zero mean and variances and , respectively.

The output of the interconnected nonlinear subsystem can be rewritten as follows:where the parameters vector and the observations vector are given as follows:with , , , and .

Similarly, the output of the interconnected nonlinear subsystem can be rewritten in the following form:where the parameters vector and the observations vector are described by with , , , and .

In the first step of the control method, we propose to estimate the parameters involved in the parameter vectors (34) and (37) using the RELS algorithm (12). The second step is reserved to the resolution of the polynomial equation:The resolution of (39) permits determining the parameters of the polynomials:withAs for the third step, it concerns resolving the following equation:wherewithwith .

Based on discriminant , the control law is given, such that(i)if , then ;(ii)if , then ;(iii)if is strictly negative, then the polynomial (42) does not admit solution. In this case, we propose that .

It can be remarked that the boundedness of these obtained controllers depends on the parameter . However, if this parameter is equal to zero at discrete-time , then the control law diverges. To overcome this problem, we must append a test on the estimated parameter value . For example, if , we suppose that in order to complete the computation of the controller involved in the iterative control algorithm.

The data for the practical implementation of the developed self-tuning regulator is given below:(1)The parameters values intervening in the mathematical model (31) are chosen as follows: , , , , , , , , and .(2)The parameters values of the mathematical model (32) are selected as follows: , , , , , , , , and .(3)The initial conditions of the estimated parameters vector and the adaptation gain matrix are chosen in such a way that and , where is the identity matrix, .(4)The forgetting factors and are chosen in such a way that and .(5)The number of measurements is selected as .(6)The two control laws and must minimize the following quadratic criteria:with and .

Some results of this simulation example of the considered system , , are given. Thus, Figures 4 and 5 represent the evolution curves of the prediction error, the parametric distance, and their overall variances for each interconnected nonlinear subsystem.

Table 1 presents the statistical average values of the prediction errors and , the parametric distances and , and their overall variances , , , and , where the statistical average and the variance of the different sizes illustrated in Table 1 can be made starting from the calculation of the following expressions:

It can be remarked that the parametric distances (statistical average) have some quite low values. We can remark also that the shapes of prediction errors have some relatively low values too. Consequently, the evolution curves of the different elements, as delineated in Figures 4 and 5, indicate the good quality of the estimate, which is obtained using the RELS algorithm (12). This estimate quality ensures exact quality of control, and it is justified by Figures 6 and 7, which present the evolution curves of the control laws and , the system outputs and , and their variances , , , and .

Table 2 presents the statistical average values and of the control law and the system output , respectively, and their variances and for each interconnected system , . These values were obtained using (46). Note that the evolution curves of the variances and converge towards minimum values. Therefore, exact regulation is achieved with the developed self-tuning controller despite the parameters variations of the interconnected nonlinear system, the presence of interactions signals, and disturbances interim on each interconnected system output.

5. Conclusions

The self-tuning control problem of large-scale nonlinear discrete-time systems composed of several SISO interconnected nonlinear stochastic subsystems has been proposed in this paper. Each interconnected nonlinear subsystem is described by Hammerstein mathematical model with unknown and time-varying parameters. A recursive algorithm for self-tuning regulation of large-scale interconnected stochastic systems has been developed on the basis of the minimum variance approach with explicit scheme. The presented adaptive control algorithm requires three steps at each control cycle. The general algorithm proposed here is combined with a real-time estimation step and it is carried out to solve the regulation problem for this class of interconnected discrete-time systems by reducing the disturbances acting on the output. The estimation step has been achieved through the adjustable model method and the least-squares techniques. An illustrative simulation example of hydraulic system was treated in order to validate our developed theoretical results.

Competing Interests

The authors declare that they have no competing interests.