This paper deals with the stability and stabilization problems for a class of discrete-time nonlinear systems. The systems are composed of a linear constant part perturbated by an additive nonlinear function which satisfies a quadratic constraint. A new approach to design a static output feedback controller is proposed. A sufficient condition, formulated as an LMI optimization convex problem, is developed. In fact, the approach is based on a family of LMI parameterized by a scalar, offering an additional degree of freedom.
The problem of performance taking into account an H∞ criterion is also investigated. Numerical examples are provided to illustrate the effectiveness of the proposed conditions.
1. Introduction
Modeling a real process is generally a complex and difficult task. Even if in numerous cases, a linear model can capture the main dynamical characteristics of a process. In some situations, it is necessary to take into account model uncertainties in order to design an efficient control law.
There exists an extensive literature dealing with this problem which is in fact the main problem in robust control design [1–8].
Among the numerous solutions allowing taking into account model uncertainties, a way which has been frequently investigated in literature consists in adding to the linear part of the model a nonlinear one which captures model uncertainties and frequently referred in literature as nonlinear systems with separated nonlinearity.
Nonlinear systems with separated nonlinearity are a class of nonlinear systems composed of a linear constant part to which another nonlinear function part is added. This function depends on both time and state and satisfies a quadratic constraint [9–16]. This class of nonlinear system can be considered as a generalized model for linear systems with parametric uncertainties where uncertainties can be norm bounded [2, 17] or polytopic [9, 10, 18, 19].
Many papers have investigated robust stability, analysis, and synthesis using essentially Lyapunov theory which has proved to be efficient in this context. Recent proposed approaches [9–13] are based on convex optimization problems involving linear matrix inequality (LMI) where the objective is to maximize the bounds on the nonlinearity that systems can tolerate without unstabilities. In particular, sufficient conditions where developed in the context of static state feedback and static or dynamic output feedback controllers [9–13, 15].
Even if the static output feedback stabilization (SOF) problem is considered as NP-hard [20] and still one of the most important open questions in the control theory, it concentrates the efforts of many researchers. SOF gains, which stabilize the system, are not easy to find due to the nonconvexity of the SOF formulation. In some papers the design of SOF controllers for a class of discrete-time nonlinear systems is proposed [11, 15, 17]. In this paper, we propose a new approach for robust static output feedback stabilization of a class of discrete-time nonlinear systems using LMI techniques. In fact, our approach is based on the introduction of a relaxation scheme to the SOF problem similar to [15, 21–23]. Our major objective is to maximize the admissible bounds on the nonlinearity guaranteeing the stability of system, with a prescribed degree μ. The main contribution is the possibility of decoupling the Lyapunov matrix to the SOF gain leading to less restrictive conditions.
The problem of performance is also treated in the context of H∞ settings. A new H∞ norm characterization is proposed for this class of nonlinear discrete time systems in terms of LMI formulation.
The paper is organized as follows. Section 2 presents robust stability condition for the class of nonlinear discrete time systems. Then, we develop our main results for robust stabilization by SOF. In Section 3, the problem of robust H∞ synthesis via SOF is presented. Section 4 presents numerical examples for robust stabilization and H∞ synthesis to illustrate the potential of the proposed conditions.
Notation 1.
For conciseness the following notations are used: sym(A)=A+AT,diag(A,B)=[A00B], [AB•C]=[ABBTC], and P=PT>0 is a symmetric and positive definite.
2. Preliminaries
In this section, we consider nonlinear discrete-time system with the following state-space representation:
x(k+1)=Ax(k)+f(k,x),
where x(k)∈ℜn is the state vector of the system. A∈ℜn×n is a constant matrix and f(k,x) a nonlinear function in both arguments k and x satisfying f(k,0)=0. This means that the origin is an equilibrium point of the system.
The nonlinear function f is bounded by the following quadratic constraints:
fT(k,x)f(k,x)≤α2xTMTMx,
where α>0 is the bounding parameter of the nonlinear function f and M is a constant matrix of appropriate dimensions.
The parameter α can be defined as a degree of robustness, because its maximization leads to an increase of robustness against uncertain perturbations. Note that constraint (2.2) is equivalent to
[xf]T[-α2MTM00I][xf]≤0.
Remark 2.1.
The nonlinear function f(k,x), which satisfies the quadratic constraint (2.2), can be considered as parameter uncertainty [17].
In the sequel, we will use the following definition to present the concept of robust stability of the system (2.1), (2.2).
Definition 2.2.
System (2.1) is robustly stable with degree α>0 if the equilibrium x=0 is globally asymptotically stable for all f(k,x) satisfying constraint (2.2).
In this section, we develop a method for studying robust stability of system (2.1). Before, we introduce some instrumentals tools which will be used in the proof of characterization of stability of system (2.1).
Lemma 2.3 (S-procedure lemma [1]).
Let Ω0(x) and Ω1(x) be two arbitrary quadratic forms over ℜn, then Ω0(x)<0 for all x∈ℜn-{0} satisfying Ω1(x)≤0 if and only if there exist τ≥0:
Ω0(x)-τΩ1(x)<0,∀x∈ℜn-{0}.
Proof.
See [1].
Lemma 2.4 (Projection lemma [1]).
Given a symmetric matrix ψ∈ℜn×n, and two matrices P, Q of column dimensions n, there exists X such that the following LMI holds:
ψ+sym(PTXTQ)<0,
if and only if the projection inequalities with respect to X are satisfied:
𝒩Pψ𝒩pT<0,𝒩QTψ𝒩Q<0,
where 𝒩P and 𝒩Q denote arbitrary bases of the nullspaces of P and Q, respectively.
Proof.
See [1].
Lemma 2.5.
Let Φ a symmetric matrix and N, J be matrices of appropriate dimensions. The following statements are equivalent:
Φ<0andΦ+NJT+JNT<0,
there exists a matrix G such that
[ΦJ+NGJT+GTNT-G-GT]<0.
Proof.
The proof is obtained remarking that (2.7) can be developed as follows:
[ΦJ+NGJT+GTNT-G-GT]=[ΦJJT0]+sym{[0I]GT[NT-I]}<0,
and by applying lemma 2.
2.1. Stability Characterization
We first introduce the following theorem which gives a robust stability conditions for system (2.1). In fact, it is described by a convex optimization problem where we try to maximize the nonlinear bounding parameter without loss of the stability of system.
Theorem 2.6.
System (2.1) is robustly stable with degree α>0 if there exist a positive definite symmetric matrix Q and a positive scalar β=1/α such that:
{minimizeβ[-Q0QATQMT•-II0••-Q0•••-βI]<0.
Proof.
See [13].
Remark 2.7.
The stability condition given by Theorem 2.6 is equivalent to the one introduced in [11] and [13].
A new condition for robust stability is proposed in the following theorem.
Theorem 2.8.
System (2.1) is stable with degree α>0 if there exist a positive definite symmetric matrix Q , a matrix G of appropriate dimensions, and a positive scalar β=1/α, such that the following optimization problem:
{minimizeβ[-Q0μQ0Q•-II00••-Q0(A-μI)G•••-βIMG••••-G-GT]<0,
is feasible for any prescribed scalar μ∈]-11[.
Proof.
Inequality (2.9) can be expressed as follows:
[-Q0QATQMT•-II0••-Q0•••-βI]=[-Q0μQ0•-II0••-Q0•••-βI]︸ϕ+sym([00A-μIM]︸N[Q000]︸J)<0.
It is not difficult to proof that ϕ<0 for any μ∈]-11[. Now writing:
N=[00A-μIM],
J=[Q000]T,
and by lemma 2, there exists a matrix G of appropriate dimensions such that inequality is satisfied.
Remark 2.9.
The two optimization problems (2.9) and (2.10) are equivalent. In the case of stability analysis or state feedback control synthesis, no improvement is obtained by problem (2.10). The main advantage of problem (2.10) will appear when dealing with static output feedback. In that case, we will see that it theoretically improves the obtained results.
2.2. Static Output Feedback Control
In this section, we investigate the static output feedback stabilization problem for nonlinear discrete systems.
We consider the nonlinear discrete-time system described as follows:
x(k+1)=Ax(k)+f(k,x,u)+Bu(k),y(k)=Cx(k),
where u(k)∈ℜm is the control input,y(k)∈ℜp is the measured output, and B∈ℜn×m and C∈ℜp×n are constant matrices. We assume that the pair (A,B) is stabilizable and C is full rank matrices. Also f(k,x,u) is a nonlinear function which satisfies the following quadratic constraints:
fT(k,x,u)f(k,x,u)≤α2(xTFTFx+uTHTHu),
where α>0 is the bounding parameter of the function f and F and H are constant matrices of appropriate dimensions.
The objective is to find a static output feedback control law such as
u(k)=Ky(k),
where K∈ℜm×p.
The closed loop system is given by the following state space representation:
x(k+1)=(A+BKC)x(k)+f(k,x,u),
and function f satisfies:
fT(k,x,u)f(k,x,u)≤α2xT(FTF+(HKC)THKC)x.
Note that in this case, the constraint (2.16) is equivalent to
[xf]T[-α2(FTF+(HKC)THKC)00I][xf]≤0.
To establish a robust stabilization theory for system (2.12) with (2.13) by SOF, we give in the following theorem, an optimization problem which allows to stabilize the linear part of system (2.15) with (2.17) and at the same time to maximize the value of parameter α.
Theorem 2.10.
System (2.12) is asymptotically stable by static output feedback with degree α>0 if there exist a positive definite symmetric matrix Q, a matrix R∈m×p, and a positive scalar β=1/α such that the following optimization problem is solvable:
{minimizeβ[-Q0(AQ+BRC)TQFT(HRC)T•-II00••-Q00•••-βI0••••-βI]<0,
where
Q=V[Q100Q2]VT.
The static output feedback gain is given by
K=RUC0Q1-1C0-1UT,
with U∈ℜp×p and V∈ℜn×n are unitary matrices, and C0∈ℜp×p matrix which are obtained by using the singular value decomposition of the matrix C :
C=U[C00]VT.
Proof.
According to the Theorem 2.6, system (2.15) is robustly stable if there exist a positive definite symmetric matrix Q and a positive scalar β=1/α such that the following optimization problem is solvable:
{minimizeβ[-Q0QÃTQMT•-II0••-Q0•••-βI]<0.
Now we defin
Ã=A+BKC,
M=(FHKC),
Q is replaced by (2.19) with: Q1∈ℜp×p, Q2∈ℜ(n-p)×(n-p).
Equation (2.22) becomes:
{minimizeβ[-Q0Q(A+BKC)TQFTQCTKTHT•-II00••-Q00•••-βI0••••-βI]<0.
Unfortunately, (2.23) is not convex in K and Q and cannot be solved by the LMI tools. We can introduce some transformations to simplify the KCQ term of the inequality (2.23) by using (2.19) and (2.21) as follows:
KCQ=KU[C00]VTV[Q100Q2]VT=KUC0Q1C0-1U-1︸RU[C00]VT︸,C
where R∈ℜm×p.
With this transformation we obtain the optimization problem given in Theorem 2.10.
Remark 2.11.
The optimization problem given by Theorem 2.10 presents a sufficient condition for the robust stabilization by SOF of discrete-time nonlinear system (2.12). In fact to solve the BMI problem, we impose a diagonal structure to Lyapunov matrix Q as in [24] and we obtain a LMI convex problem.
2.3. Main Results
In this section, we introduce a new approach for robust stabilization by SOF of nonlinear discrete time system (2.12). The following results present solutions to static output feedback problem in which an improved sufficient condition is presented. In fact, this approach is derived from Theorem 2.8.
Theorem 2.12.
The system (2.12) is robustly stable by static output feedback with degree α>0, for an arbitrary prescribed number μ in ]-11[ if there exist a positive definite symmetric matrix Q, matrices R∈ℜm×p, G∈ℜn×n, and a positive scalar β=1/α such that the following optimization problem is solvable:
{minimizeβ[-Q0μQ00Q•-II000••-Q00AG+BRC-μG•••-βI0FG••••-βIHRC•••••-G-GT]<0,
where
G=V[G10G2G3]VT.
The static output feedback gain is given by
K=RUC0G1-1C0-1UT,
with U,V, and C0 are given in (2.21).
Proof.
According to Theorem 2.8, the closed loop system (2.15) is robustly stable if there exist a positive definite symmetric matrix Q, a matrix G of appropriate dimensions, and a positive scalar β=1/α such that
{minimizeβ[-Q0μQ0Q•-II00••-Q0(Ã-μI)G•••-βIMG••••-G-GT]<0,
where the following hold:
Ã=A+BKC,
M=(FHKC),
G replaced by (2.26) where G1∈ℜp×p, G2∈ℜ(n-p)×p, G3∈ℜ(n-p)×(n-p),
with U,V, and C0 are given by (2.21).
Then, inequality (2.28) becomes
[-Q0μQ00Q•-II000••-Q00AG+BKCG-μG•••-βI0FG••••-βIHKCG•••••-G-GT]<0.
Unfortunately, (2.29) is not convex and cannot be solved by the LMI tools.
For this reason, we introduce some transformations to simplify the KCG term of the inequality (2.29) by using (2.26) and (2.21) as follows:
KCG=KU[C00]VTV[G10G2G3]VT=KUC0G1C0-1U-1︸RU[C00]VT︸C,
where R∈ℜm×p.
With this transformation, we obtain the optimization problem given by Theorem 2.12.
The following lemma gives a connection of the results of Theorem 2.10 with the one of Theorem 2.12.
Lemma 2.13.
If the SOF stabilization problem is solvable by Theorem 2.10, then it is solvable by Theorem 2.12.
Proof.
If we consider the optimization problem (2.25) with G=Q andμ=0, we obtain the optimization problem (2.18) with Q satisfying (2.19). Therefore, if (2.25) is feasible, then (2.19) is feasible too.
Stability is the minimum requirement and in practice a performance level has to be guaranteed. Performance objectives can be achieved via H∞ norm optimization. In this section, we study the H∞ control problem for the following nonlinear discrete-time system:
x(k+1)=Ax(k)+Bw(k)+f(k,x),z(k)=Cx(k)+Dw(k),
where w(k)∈Rq is the exogenous disturbance and z(k)∈Rr is the controlled output. A,B,C, and D are known constant matrices of appropriate dimensions. f(k,x) is a nonlinear function satisfying the following quadratic constraints:
[xwf]T[-α2MTM0000000I][xwf]≤0.
We note Twz the transfer matrix from input w to output z when f(k,x)=0 is expressed by
Twz=C(sI-A)-1B+D.
Theorem 3.1.
System (3.1) is robustly stable with ∥Twz∥∞<γ for a prescribed constant value γ>0, if there exist a positive definite symmetric matrix Q and positive scalars τ>0 and β=1/α such that the following optimization problem is feasible:
{minβ[-Q00QATQMTQCT•-τγ2I0τBT0τDT••-II00•••-Q00••••-βI0•••••-τI]<0.
Proof.
Define the Lyapunov function:
V(x)=xT(k)Px(k),
where P=PT>0.
By using the dissipative theory, we show that
V(k+1)-V(k)+zT(k)z(k)-γ2wT(k)w(k)<0,
where γ>0 is a prescribed scalar such that
∥Twz∥∞<γ∞<γ,
where γ∞ is the corresponding bounding norm bound when f=0.
Evaluating (3.6) leads to
[xTwTfT][ATPA-P+CTCATPB+CTDATPBTPA+DTCBTPB+DTD-γ2IBTPPAPBP][xTwTfT]T<0.
Now applying the S-procedure Lemma to (3.8) with (3.2), there exits τ>0 such that
[ATP̃A-P̃+τ-1CTC+α2MTMATP̃B+τ-1CTDATP̃BTP̃A+τ-1DTCBTP̃B+τ-1DTD-τ-1γ2IBTP̃P̃AP̃BP̃-I]<0,
where P̃=P/τ.
By Schur complement, (3.9) is equivalent to
[-P00ATMTCT•-τ-1γ2I0BT0DT••-II00•••-P-100••••-βI0•••••-τI]<0,
where β=α-2.
Multiplying by diag(Q,τ,I,I,I) with Q=P̃-1 the both sides of (3.10), we obtain the optimization problem (3.4).
Now, we introduce the following theorem which can be seen as an alternate characterization of upper bounds of the H∞ norm of.
Theorem 3.2.
System (3.1) is robustly stable with ∥Twz∥∞<γ for a prescribed constant value γ>0, if there exist a positive definite symmetric matrix Q, a matrix G of appropriate dimensions, and positive scalars τ>0 and β=1/α such that for any prescribed scalar μ in ]-11[ , the following optimization problem is feasible:
{minβ[-Q00μQ00Q•-τγ2I0τBT0τDT0••-II000•••-Q00(A-μI)G••••-βI0MG•••••-τICG••••••-G-GT]<0.
Proof.
Inequality (3.4) can be written as follows:
[-Q00QATQMTQCT•-τγ2I0τBT0τDT••-II00•••-Q00••••-βI0•••••-τI]=[-Q00μQ00•-τγ2I0τBT0τDT••-II00•••-Q00••••-βI0•••••-τI]+sym([000A-μIMC][Q00000])<0.
By lemma 2, denoting
ϕ=[-Q00μQ00•-τγ2I0τBT0τDT••-II00•••-Q00••••-βI0•••••-τI]<0,for anyμ∈]-11[N=[000A-μIMC],J=[Q00000],
there exists a matrix G of appropriate dimensions such that (3.11) holds where γ,β, and τ are scalars previously defined.
3.1. Static Output Feedback H∞ Synthesis
In this section, we consider the static output feedback stabilization problem for the following nonlinear discrete system (3.14):
x(k+1)=Ax(k)+Buu(k)+Bww(k)+f(k,x,u),z(k)=Czx(k)+Dzuu(k)+Dzww(k),y(k)=Cyx(k),
where u(k)∈ℜm is the control input, w(k)∈ℜq is the exogenous disturbance, z(k)∈ℜr is the controlled output, and y(k)∈ℜp is the measured output. Also, Bu,Bw,Cz,Dzu,Dzw and Cy are known constant matrices of appropriate dimensions.
The system closed by SOF is written as
x(k+1)=Aclx(k)+Bclw(k)+f(k,x,u),z(k)=Cclx(k)+Dclw(k),
where the following hold:
Acl=A+BuKCy,Bcl=Bw,Ccl=Cz+DzuKCy,Dcl=Dzw.
The objective of this section is to design static output feedback H∞ controllers for nonlinear discrete time system (3.14).
Theorem 3.3.
System (3.14) is robustly static output feedback stabilizable with ∥Twz∥∞<γ for a prescribed constant value γ>0 if there exist a positive definite symmetric matrix Q, a matrix R of appropriate dimensions, and positive scalars τ>0 and β=1/α such that the following optimization problem is feasible:
{minβ[-Q00QAT+(BuRCy)TQFT(HRCy)TQCzT+(DzuRCy)T•-τγ2I0τBwT00τDzwT••-II000•••-Q000••••-βI00•••••-βI0••••••-τI]<0,
where
Q=V[Q100Q2]VT.
The static output feedback gain is then given by
K=RUC0Q1-1C0-1UT,
where U,V, and C0 are given in (2.21).
Proof.
According to Theorem 3.1, system (3.14) is robustly stable with ∥Twz∥∞<γ for a prescribed constant value γ>0 if there exist a positive definite symmetric matrix Q and positive scalars τ>0 and β=1/α such that the following optimization problem is feasible:
{minβ[-Q00QAclTQMTQCclT•-τγ2I0τBclT0τDclT••-II00•••-Q00••••-βI0•••••-τI]<0.
Now:
Acl,Bcl,Ccl and Dcl are replaced by (3.16),
M=(FHKCy),
Q is replaced by (3.18) where : Q1∈ℜp×p, Q2∈ℜ(n-p)×(n-p).
After same direct developments the result follows.
3.2. An Improved Approach of Static Output Feedback Synthesis for H∞ Robust Control
In this paragraph, an H∞ robust control for nonlinear systems (3.14) improving the previous approach is proposed.
Theorem 3.4.
System (3.14) is robustly SOF stabilisable with ∥Twz∥∞<γ for a prescribed constant value γ>0 if there exist a positive definite symmetric matrix Q, matrices G and R of appropriate dimensions, and positive scalars τ>0, and β=1/α such that the following optimization problem is feasible for any prescribed scalar μ in ]-11[:
[-Q00μQ000Q•-τγ2I0τBwT00τDzwT0••-II0000•••-Q000AG-μG+BuRCy••••-βI00FG•••••-βI0HRCy••••••-τICzG+DzuRCy•••••••-G-GT]<0,
where
G=V[G10G2G3]VT.
The static output feedback gain is given by
K=RUC0G1-1C0-1UT,
with R∈ℜm×p, U,V,and C0 are given in (2.21).
Proof.
According to the Theorem 3.3, the system (3.14) is robustly stable with ∥Twz∥∞<γ for a prescribed degree γ>0 if there exist a positive definite symmetric matrix Q, a matrix G of appropriate dimensions, and positive scalars τ>0 and β=1/α, for any prescribed scalar μ in ]-11[, such that the following optimization problem is feasible:
{minβ[-Q00μQ00Q•-τγ2I0τBclT0τDclT0••-II000•••-Q00(Acl-μI)G••••-βI0MG•••••-τICclG••••••-G-GT]<0,
where
Acl,Bcl,Ccl, and Dcl are replaced by (3.16),
M=(FHKCy),
G replaced by (3.22), where G1∈ℜp×p, G2∈ℜ(n-p)×p.
The optimization problem (3.24) is expressed as follows:
{minβ<0[-Q00μQ000Q•-τγ2I0τBwT00τDzwT0••-II0000•••-Q000AG+BuKCyG-μG••••-βI00FG•••••-βI0HKCyG••••••-τICzG+DzuKCy•••••••-G-GT]<0.
For this reason, we express differently the KCyG term of the inequality (3.25), by using (3.22) and (2.21), in the same way as in (2.30). Consequently, we obtain
KCyG=RCy.
With this transformation, we obtain the optimization problem given in Theorem 3.4.
4. Numerical Examples
We present in this section two numerical examples to illustrate the proposed theory for SOF synthesis.
Our approach Theorem 2.12 is compared to the methods presented in [11].
Example 4.1 (see [5]).
We consider the nonlinear discrete-time system of inverted pendulum (see Figure 1).The system can be described by
x(k+1)=Ax(k)+B(u(k)+αsin y(t)),y(k)=Cx(k),
where x(k)=[x1(k)x2(k)]T, x1(k) and x2(k) being, respectively, the angular displacement and velocity; u(k) represents the field current of the DC motor; the matrices A∈ℜ2×2,B∈ℜ2×1, and C∈ℜ1×2are given by:
A=[10.095200.9048];B=[0.0480.952];C=[10].
The nonlinear function f is given by the following expression:
f(k,x)=αBsin y(k)=αBsinx1(k).
This function satisfies
fT(k,x)f(k,x)=α2BTBsin2x1(k)=0.9α2sin2x1(k)≤α2xT[0.95000][0.95000]x.
In this case, we take F=[0.95000];H=0.
The nonlinear system is unstable. The matrix A is unstable which means that the linear part of systems is unstable. Therefore, we apply the approach given by Theorems 2.10 and 2.12 to stabilize system by SOF. We summarize in Table 1 the obtained results.Figures 2 and 3 show the output trajectory obtained respectively by applying Theorems 2.10 and 2.12 to system (3.26) with α=αmax.We can see a drastic improvement obtained by Theorem 2.12.
Numerical Evaluation for Example 4.1.
Approach
αmax
K
Theorem 2.10 condition
0.8254·10-4
-1.9987
Theorem 2.12 condition for μ=0.7
0.06688
-0.3101
Inverted pendulum scheme.
Output y=x1.
Output y=x1.
Example 4.2.
We consider now the nonlinear discrete-time system:
x(k+1)=[-0.07250.19571.59310.03010.04040.60840.3764-0.16350.9024]x(k)+[1.21230.28950.61740.06510.93790.5110]u(k)+f(k,x(k),u(k)),y(k)=[0.9313-0.7534-0.0335120.5]x(k),
with the nonlinear function f satisfying the constraint (2.17), where F=I3 and H=I2. The results obtained byapproaches of Theorem 2.12 and in [11] are given in Table 2.The new approach (Theorem 2.12) improves the results obtained by [11].
Numerical Evaluation for Example 4.2.
Approach
αmax
K
[11]
0.1954
[-0.1329-0.1268-0.4309-0.0405]
New approach, (2.25), Theorem 2.12 for μ=0
0.3084
[-0.0019-0.2766-1.15420.1067]
New approach, (2.25), Theorem 2.12 for μ=-0.35
0.3329
[0.0357-0.3141-0.98120.0924]
Example 4.3.
We consider the nonlinear discrete-time system (3.14) with
A=[0.81890.08630.09000.08130.25241.00330.03130.2004-0.05450.01020.7901-0.2580-0.1918-0.10340.16020.8604],Bu=[0.00450.00440.10010.01000.0003-0.0136-0.00510.0926],Bw=[0.0953000.0145000.086200-0.001100],Cz=[10-1000000000],Dzu=[001001],Cy=[10000010],Dyw=0,
with the nonlinear function f satisfying the constraint (3.2), where F=I4 and H=I4. In Table 3, we present numerical result for H∞ performance via SOF by applying Theorem 3.3 and for different values of γ which satisfy (3.7).In Table 4, we present numerical result for H∞ performance via SOF by applying Theorem 3.4 and for different values of γ which satisfy (3.7).
Numerical Evaluation for Example4.3 by Theorem 3.3.
Approach
αmax
K
Approach, (3.17), Theorem 3.3 for γ∞=1.9918,γ=2.2361
0.0032
[-2.0416-0.17960.6545-1.0362]
Approach, (3.17), Theorem 3.3 for γ∞=1.9918,γ=5
0.0288
[-2.1323-0.08081.1624-1.2381]
Numerical Evaluation for Example4.3 by Theorem 3.4.
Approach
αmax
K
Approach, (3.21), Theorem 3.4 for γ∞=1.9918,γ=2.2361,μ=0.22
0.3643
[-0.2993-0.1328-0.0685-0.0743]
Approach, (3.21), Theorem 3.4 for γ∞=1.9918,γ=5,μ=0.2400
0.4087
[-0.3376-0.1185-0.0660-0.0737]
5. Conclusion
In this paper, the stabilization problem by static output feedback (SOF) for a particular class of nonlinear discrete time systems is investigated. A new sufficient condition is elaborated by using Lyapunov theory and formulated by LMI constraints. We obtain a convex optimization problem for maximizing the bound of the nonlinearity preserving the stability of the systems.
Finally, the proposed controller design method was extended to incorporate H∞ synthesis. An optimization problem, which is linear both in the admissible nonlinearity bound and the disturbance attenuation, is developed. Numerical comparisons with existing methods in literature illustrate the improvement obtained by our approaches. All of them can also be extended to the dynamic output feedback.
Other classes of nonlinear discrete time nonlinear systems exist and present some interesting characteristic from a practical point of view. For example, the ones where the linear part is affected by polytopic uncertainties. It would be interesting to extend the results of this paper to those classes. This will be exploited in a near future.
BoydS.El GhaouiL.FeronE.BalakrishnanV.199415Philadelphia, Pa, USASIAMxii+193SIAM Studies in Applied MathematicsMR1284712ZBL0816.93004GarciaG.BernussouJ.ArzelierD.Robust stabilization of discrete-time linear systems with norm-bounded time-varying uncertainty1994225327339MR1274907ZBL0820.93059HaddadW. M.BernsteinD. S.Parameter-dependent Lyapunov functions and the Popov criterion in robust analysis and synthesis1995403536543MR131926110.1109/9.376077ZBL0821.93063FridmanE.Robust sampled-data control of linear singulary perturbed systemsProceedings of the 44th IEEE Conference on decision and Control, and the European Control ConferenceDecember 2005Seville, Spain4324432910.1080/00207170110120026MahmoudM. S.XieL.Positive real analysis and synthesis of uncertain discrete time systems200047340340610.1109/81.841924BahnasawiA. A.Al-FuhaidA. S.MahmoudM. S.Linear feedback approach to the stabilisation of uncertain discrete systems198913614752MR976147ZBL0664.93072GalkowskiK.RogersE.XuS.LamJ.OwensD. H.LMIs-a fundamental tool in analysis and controller design for discrete linear repetitive processes2002496768778MR191128210.1109/TCSI.2002.1010032XuS.LamJ.New positive realness conditions for uncertain discrete descriptor systems: analysis and synthesis200451918971905MR209035010.1109/TCSI.2004.834524StipanovićD. M.ŠiljakD. D.Robust stability and stabilization of discrete-time non-linear systems: the LMI approach2001749873879MR184712610.1080/00207170010038712ZBL1022.93034ŠiljakD. D.StipanovićD. M.Robust stabilization of nonlinear systems: the LMI approach200065461493MR2034962ZBL0968.93075HoD. W. C.LuG.Robust stabilization for a class of discrete-time non-linear systems via output feedback: the unified LMI approach2003762105115MR196584710.1080/0020717031000067367LuG.HoD. W. C.Generalized quadratic stability for continuous-time singular systems with nonlinear perturbation2006515818823MR223260610.1109/TAC.2006.875017ZuoZ.WangJ.HuangL.Robust stabilization for non-linear discrete-time systems2004774384388MR205293310.1080/00207170410001663543StankovićS. S.StipanovićD. M.ŠiljakD. D.Decentralized dynamic output feedback for robust stabilization of a class of nonlinear interconnected systems2007435861867MR2306733ZBL1117.93057BediouiN.bedioui_neila_enit@yahoo.frSalhiS.salhis@lycos.comKsouriM.mekki.ksouri@insat.rnu.tnAn improved LMI approach for robust static output feedback stabilization of nonlinear systemsProceedings of the 5th IEEE International Multi-Conference on Systems, Signals and Devices (SSD '08)July 2008Amman, Jordan1610.1109/SSD.2008.4632859AbbaszadehM.MarquezH. J.Robust static ouput feedback dtabilization of discrete time nonlinear uncertain systems with H∞ performanceProceedings of the 16th Meditteranean Conference on Control and Automation, CongressJune 2008Ajaccio, France22623110.1109/MED.2008.4602121CauëtS.RambaultL.BachelierO.Parameter-dependent Lyapunov stability analysis: application to an induction motor1Proceedings of the American Control Conference (ACC '01)June 2001Arlington, Va, USA16016510.1109/ACC.2001.945532AbbaszadehM.MarquezH. J.Robust H∞ observer design for a class of nonlinear uncertain systems via convex optimizationProceedings of the American Control Conference (ACC '07)July 2007New York, NY, USA1699170410.1109/ACC.2007.4282280KauS.-W.LiuY.-S.HongL.LeeC.-H.FangC.-H.LeeL.A new LMI condition for robust stability of discrete-time uncertain systems2005541211951203MR2175634ZBL1129.93482TokerO.ÖzbayH.On the NP-hardness of solving bilinear matrix inequalities and simultaneous stabilization with static output feedbackProceedings of the American Control Conference (ACC '95)June 1995Seattle, Wash, USA25252526de OliveiraM. C.GeromelJ. C.BernussouJ.Extended H2 and H∞ norm characterizations and controller parametrizations for discrete-time systems2002759666679MR191643210.1080/00207170210140212ZBL1029.93020MehdiD.BoukasE. K.BachelierO.Static output feedback design for uncertain linear discrete time systems2004211113MR203700810.1093/imamci/21.1.1ZBL1049.93069de OliveiraM. C.GeromelJ. C.BernussouJ.LMI optimization approach to multiobjective controller design for discrete-time systems4Proceedings of the 38th IEEE Conference on Decision and Control (CDC '99)December 1999Phoenix, Ariz, USA3611361610.1109/CDC.1999.827913CrusiusC. A. R.TrofinoA.Sufficient LMI conditions for output feedback control problems199944510531057MR169055510.1109/9.763227ZBL0956.93028