Enhanced Control Design of Discrete-Time Takagi-Sugeno State-Multiplicative Noisy Systems

Design conditions for existence of the H ∞ state feedback control for Takagi-Sugeno fuzzy discrete-time stochastic systems with state-multiplicative noise, stabilizing the closed-loop in such way that the quadratic performance in the mean is satisfied, are presented in the paper. Using newly introduced enhanced form of the bounded real lemma for such stochastic systems, the LMI-based procedure is provided for computation of gain matrices of the state control law, realized in the parallel distributed compensation structure. The approach is illustrated on an example, demonstrating the validity of the proposed method.


Introduction
The control design for systems with stochastic uncertainties is an area of study for several decades (see [1] and the references therein) and so far these practically and theoretically attractive fields derived numerous solutions.An approach, in which the parameter uncertainties were modeled as white noise processes in a linear system, was firstly applied for discrete-time systems [2] and continuous-time systems [3][4][5].Unfortunately, this way of doing it on linear quadratic optimal control principle and reflecting the control input as a deterministic function of system state naturally led mainly to a nonstandard Riccati equation.Representing the system uncertainty by unknown disturbance signal in the system equations,  ∞ stayed the most prominent method of dealing with the corresponding control problem of disturbance relaxation [6].Combining both representations of uncertainty, adapting  ∞ approach towards the disturbance attenuation for stochastic systems, and including those that ensure a performance bound in the  ∞ sense, new accessions have been made in diverse practical problems for finite and infinite time horizons for continuous-time [7,8] as well as for discrete-time stochastic systems [9][10][11], denoting later such systems as linear stochastic systems with multiplicative noise.The stochastic control for linear and nonlinear systems with multiplicative noises has founded its applications in many fields in control theory.Such models are encountered mainly in constrain control [12] and in gain scheduling when the scheduling parameters are corrupted with measurement noise [13].
In contrary to the linear framework, nonlinear systems are too complex to be represented by unified mathematical resources, and so a generic method has not been developed yet to design a controller valid for nonlinear stochastic systems with multiplicative noise [14].An alternative is Takagi-Sugeno (TS) fuzzy approach [15], which benefits from the advantages of suitable linear approximation.Using the TS fuzzy model, each rule utilizes the local system dynamics by a linear model and the nonlinear system is represented by a collection of fuzzy rules.Recently, TS model based fuzzy control approaches are being fast and successfully used in nonlinear control frameworks.As a result, a range of stability analysis conditions [16], as well as control design methods [17][18][19][20], have been developed for TS fuzzy systems, relying mostly on the feasibility of an associated set of linear matrix inequalities (LMI) [21,22].An important fact is that the design problem is a standard feasibility problem with several LMIs.According to the T-S fuzzy model with multiplicative noise, there exist only few approaches investigated to solve the abovementioned problems of control design and stability analysis for continuous-time [23,24], as well as for discrete-time TS models [25][26][27] and these papers only present some results on stability analysis, mean-square stabilization, and control for TS fuzzy stochastic systems with multiplicative noises.
Considering actual results in bounded real lemma forms for discrete-time stochastic systems [8,14], new design conditions, based on an enhanced form of the bounded real lemma for TS fuzzy discrete-time stochastic systems with state-multiplicative noise and  ∞ control of such systems, are derived in the paper.For some potential numerical problems, there are also posed viewpoints and clarification for discretetime TS fuzzy model [28], used in simulation.
To present this, the paper is arranged in the next sections as follows.Following the introduction in Section 1, the control design task for discrete-time TS fuzzy stochastic systems with state-multiplicative noise is presented in Section 2. The preliminary results, focused in consequence on two rudimentary bounded real lemma forms for such defined stochastic systems, are presented in Section 3 and, subsequently, Section 4 provides the controller design conditions in equivalent forms of LMIs.Section 5 illustrates the control design task by a numerical solution and Section 6 draws some conclusions.
Throughout the paper, the notations are narrowly standard in such way that x  , X  denote the transpose of the vector x and matrix X, respectively, X = X  > 0, (≥0), means that X is a symmetric positive definite (semidefinite) matrix, diag[⋅] designates a block diagonal matrix, the symbol I  indicates the th-order identity matrix, Z + is the set of all positive integers, R indicates the set of real numbers,   {⋅} remits expectation with respect to the variable , R  , R × refer to the set of all -dimensional real vectors and  ×  real matrices, respectively,  2 ⟨0, +∞) entails the space of square summable functions over ⟨0, +∞), and   is the Kronecker delta function [29].

Problem Formulation
Through the paper, the task is concerned with state feedback design to control a TS fuzzy stochastic discrete-time dynamic system, given by the set of equations where the initial condition q(0) is known, q() ∈ R  , u() ∈ R  , and y() ∈ R  are the state, input, and output vectors, respectively, d() ∈ R   is an exogenous disturbance vector, which is the vector of the premise variables, where , V are the numbers of fuzzy rules and premise variables, respectively.It is supposed in the next that all premise variables are measurable and independent of u().More details can be found, for example, in [18,21].
The problem of interest is to design in the mean square sense stable closed-loop system by the overall fuzzy state feedback controller in the parallel distributed compensation (PDF) form which uses the same set of normalized membership functions (3) and K  ∈ R × ,  = 1, 2, . . .,  is the set of feedback controller gain matrices.It is supposed that all couples (F ℎ , G ℎ ) are controllable.

Basic Preliminaries
Lemma 1 (a TS sector approximation of the function sine).
The function sine can be approximated on the interval − <  <  as where  1 ,  2 are the sector bounds and the fuzzy membership functions take the forms Proof (compare [30]).Since it can be considered that and (3) implies then (7) can be written as Thus, evidently, the real fuzzy membership functions have to take the forms Approximating for the left bounded sector, that is,  1 = 0,  2 = , then ( 12) implies (8).
It is evident that if  1 = 0 and  2 = , then (7) gives a trivial solution Remark 2. Very often the fuzzy membership functions in the above given task are defined as simple triangles [31]; that is, Moreover, the minimal sector value is often considered approximately as  1 = 10 −2 /, or  1 = (10 −2 /).
Lemma 3. If M, N are matrices of appropriate dimension and X is a symmetric positive definite matrix, then Proof.Since X = X  > 0, then respectively.It is evident that (17) implies (15).
Lemma 5 (bounded real lemma).The deterministically unforced system (1) and ( 2) is stable in the mean square sense and with the quadratic performance if there exist a positive definite matrix P ∈ R × and a positive scalar  ∈ R such that for all ℎ = 1, 2, . . ., Here and hereafter, * denotes the symmetric item in a symmetric matrix.
Proof.Let the Lyapunov function candidate be where  ∈ R is square of the  ∞ norm of the transfer function matrix for the disturbance input d and the system output y.
Taking expectation with respect to (), then it yields =   {q  ( + 1) Pq ( + 1)} − q  () Pq () Solving for the deterministically unforced system (1) and ( 2), the following is obtained: Considering (19), then the substitution of ( 5) and ( 26) in (25) gives Mathematical Problems in Engineering Defining the composite vector the next inequality has to be satisfied if Writing (30) as follows and using the Schur complement property, (31) can be rewritten in the next form: Applying twice again the Schur complement property, (32) implies (23).This concludes the proof.
Remark 6.The inequality (23) represents the standard bounded real lemma form (compare, e.g., [32, page 66], for continuous-time linear systems).If it is necessary to reduce the range of setting (defined by ), Lyapunov function candidate can be selected in the form Now the conditions that also satisfied   {V(q( + 1))} − V(q()) ≤ 0 are Note, in the same way, the structure of ( 35) is used in the following lemmas and theorems.
Lemma 7 (enhanced bounded real lemma).The deterministically unforced system (1) and ( 2) is stable in the mean square sense and with the quadratic performance if there exist positive definite matrices P, Q ∈ R × and a positive scalar  ∈ R such that for all ℎ = 1, 2, . . ., Proof.Since ∑  ℎ=1  ℎ (()) = 1, for the deterministically unforced system (1), it is implied that Then, with an arbitrary symmetric positive definite matrix and also respectively.Adding (40) as well as its transposition to (25) gives Defining the composite vector and using (19), then (41) can be written as where, for all ℎ, Using Schur complement property to construct the LMI form, (44) can be rewritten as (37).This concludes the proof.
Remark 8.The inequality (37) is an enhanced representation of bounded real lemma for a given class of TS stochastic systems.The inequality is linear with respect to the system variables and does not involve any product of the Lyapunov matrix P and the system matrices F ℎ , W ℎ , V ℎ .This offers its main possibility to be applied in systems with polytopic uncertainties.

State Feedback Control
This section addresses the problem of finding the state feedback control law (6), working on the parallel distributed compensation concept that stabilizes the system (1), ( 2) and achieves some level of attenuation.Moreover, since the LMI solvers mostly limit the maximal number of LMIs entering a solution, the number of LMIs in the following design conditions is minimized.
Applying this result directly to system (1), (2) would lead, in the case of a large number of rules, to a possibly conservative analysis, since a common symmetric positive definite matrix P is used to verify all Lyapunov inequalities.Considering the enhanced bounded real lemma form ( 36), (37), main aim of next theorem is to extend the affine TS model properties using slack matrix variables to decouple Lyapunov matrix and the system matrices in LMIs.
Note, the LMIs in Theorem 11 are linear and the equivalent Lyapunov matrix T, used to verify all inequalities, is a symmetric positive definite matrix.

Illustrative Example
To demonstrate the control properties, the trailer-like mobile robot kinematics is considered [33], described by the following set of differential equations: where () is the steering angle (rad),  1 () is the angle difference between truck and trailer (rad),  2 () is the angle of trailer (rad),  3 () is  position of rear end of trailer (m),  4 () is  position of rear end of trailer (m), V is the constant speed of the backward movement (m/sek),  is the length of truck (m), and  is the length of trailer (m).
The control purpose is to realize the backward movement of the trailer-like mobile robot along the  3 () = 0, that is, to regulate  1 () ∼  3 () by manipulating the steering angle (), where, without forward movements, lim Consequently, the variable  4 () is not necessary to take into account.
If  1 () and () are always small values, the above model can be simplified as Applying Euler discretization, the next difference equation forms are obtained (compare, e.g., [34]): where Δ (s) is the sampling period.Defining the premise variable and the sector values as then (12) implies and the solitary nonlinear difference equation can be approximated as Considering the external disturbance () and the multiplicative noise () only in the indirect measurement of  1 (), it is given that [35] where  is the constant gain of multiplicative noise and  is the constant gain of disturbance.Subsequently, the output equation is defined as Now, introducing the state vector the discrete-time TS fuzzy model takes the form where Note, the main reason of using the above approximation for  3 ( + 1) is to obtain a controllable pair (F 1 , G). Be aware that for very small value of , the stability of the system may be close to the limit.
To demonstrate the control properties, the next system parameters were considered: Solving ( 64)-(66) with respect to the LMI variables X, T, Y 1 , Y 2 , and  using Self-Dual-Minimization (SeDuMi) package of MATLAB [36], the feedback gain matrix design problem was feasible with the results ,

X
giving the feedback gain matrices In consequence, the linear state feedback control law was obtained for such system.The obtained control law insures the stable control in the mean square sense.The closed-loop subsystem matrix eigenvalue spectra are as follows: which implies also the whole closed-loop system stability.
To compare, using the same system parameters, ( 46 and, by Theorem 10, the design parameters were computed as follows: which give less relative damping of closed-loop system solutions than there were obtained in the first case.Since the conditions of Theorems 10 and 11 are all satisfied and the design tasks were feasible, it is possible to conclude that the TS fuzzy control of given stochastic system with statemultiplicative noise is asymptotically stable in the mean with the quadratic performance by applying the linear control law. The simulation presents the closed-loop system properties, where the multiplicative noise parameters were chosen as in (5), the disturbance was modeled as noise with covariance  2  = 0.2, and the initial system state vector was (in degrees) q  (0) = [0.Figure 1 shows the output response () of the controlled nonlinear system with the TS fuzzy controller (6) designed according to the conditions of Theorem 11, where the associated input signal () is presented in Figure 2.For comparison, Figures 3 and 4 show the responses of () and the control input () for the closed-loop system with the fuzzy controller designed according to the conditions of Theorem 10, respectively.

Concluding Remarks
The paper presents the control design principle for TS fuzzy discrete-time stochastic multivariable dynamic systems with state-multiplicative noise.The stability of the control scheme is established in the mean square sense and with the quadratic performance using an enhanced representation of bounded real lemma for such TS fuzzy stochastic systems.As a result, the Lyapunov matrix and the system parameter matrices are decoupled in the resulting LMIs.This provides a suitable way for determination of state control by solving these naturally affine LMI problems.Compared with the previous results, the number of assumptions is reduced since the solutions obtained according to the conditions of both presented theorems tend to receive a linear state controller.Presented applications can be considered as a task concerning the class of  ∞ stabilization control problems where the design conditions were newly formulated.This formulation poses the problem as a stabilization problem with a fuzzy state feedback controller in the parallel distributed form whose gain matrices take no special structures and allows finding a solution to the control law without restrictive assumptions and additional specifications on the design parameters.The procedures used here for discrete-time TS fuzzy stochastic systems with state-multiplicative noise can be similarly extended for the continuous-time case.

Figure 1 :Figure 2 :
Figure 1: The response of () for the system with TS fuzzy controller designed according to the conditions of Theorem 11.

Figure 3 :Figure 4 :
Figure 3: The response of () for the system with TS fuzzy controller designed according to the conditions of Theorem 10.