Abstract

This paper presents a convex approach for nonlinear descriptor systems with multiple delays; it allows designing delayed nonlinear controllers such that the closed-loop system holds exponential estimates for convergence. The proposal takes advantage of an equivalent convex representation of the given descriptor model together with Lyapunov-Krasovskii functionals; thus, the conditions are in the form of linear matrix inequalities, which can be efficiently solved by commercially available software. To avoid possible saturation in the actuators, conditions for bounding the control input are also given. Numerical and academic examples illustrate the performance of the proposal.

1. Introduction

In the last decades, a large number of results concerning the analysis and stabilization of systems by means of the direct Lyapunov method [1], since the publication of the book [2], linear matrix inequalities (LMIs) have become a preferred solution to many control problems [3, 4], as they can be effectively solved by means of convex optimization techniques already implemented in commercially available software [5, 6]. These ideas have been extended to the analysis of time-delay systems (TDSs) via LyapunovKrasovskii (L-K) functionals [7] or LyapunovRazumikhin (L-R) functions [8]. In this context, there are several results that provide sufficient stability conditions using LMI-based approaches for different classes of TDS, such as linear time-delay systems [913], uncertain linear time-delay systems [1416], neutral linear systems [1720], systems with uncertain time-invariant delays [21], descriptor system approach for TDS [22], linear parameter-varying (LPV) time-delay systems [23], systems with time-varying delays [2429], exponential estimates for TDS [30, 31], systems with polytopic-type uncertainties [32], singular systems [33], neural networks with time delay [34, 35], and genetic regulatory networks with probabilistic time delays [36]. Recently, in [37] convex approaches are employed to provide robust stability conditions based on quasi-polynomials.

In general, delays are undesirable phenomena, because they can destabilize or produce a poor performance in the system response. However, in recent years, it has been shown that delays can also stabilize and improve the close-loop performance of a system. Moreover, the deliberate induction of delays by means of the control law is an efficient alternative to stabilize systems [38, 39]; these types of controllers are known as delayed ones. For example, in [4042], a proportional control with an appropriate delay replaces a traditional proportional-derivative one; thus, the system response is fast and insensitive to high-frequency noise. In [43], a scheme called time-delayed feedback control (TDFC) is proposed, originating different investigations [4455].

As mentioned above, LMI-based approaches have become important in the control community; however, in the context of TDS, there is an inherent conservatism for stability and stabilization conditions even for linear setups [37]; this leaves room for improvements. Moreover, a problem little explored by the community is obtaining stability conditions for nonlinear TDS. Although, originally, LMI-based stability conditions were given for linear time-invariant (LTI) systems, these have also been used on LPV [3] and nonlinear setups via exact Takagi-Sugeno (TS) [56]. The latter case employs the sector of no linearity approach [57] which allows rewriting the original nonlinear model as a convex one by means of scalar convex functions that capture uncertainties and nonlinearities. This technique has also been applied to a class of nonlinear TDS; for instance, in [58, 59], sufficient LMI conditions are proposed; in [60], uncertain TS systems are considered; in [61], sufficient LMI conditions have been given for a class of nonlinear systems. A larger family of functionals is explored in [62]. Nonetheless, none of these previous works deal with nonlinear descriptor systems; they appear when using the EulerLagrange formalism for modeling plants [63]. In the context of convex descriptor models without delays, there are some recent works [64, 65]; time-delay nonlinear descriptor systems are a few works in the literature; for instance, in [66], LMI stability and stabilization conditions have been developed for systems with only one time-delay.

Contribution: this paper proposes an LMI methodology for analysis and stability of nonlinear descriptor systems with multiple delays, thus overcoming recent results in the literature. For example, the work [4] only considers linear in standard form systems with multiple delays, [52] only studies nonlinear systems in standard form, and [66] treats nonlinear systems in descriptor with one delay. Additionally, to avoid possible saturation in the actuators, LMI conditions for bounding the control signal are established. Numerical and academic examples illustrate that including delays in the controller can reduce noise in the control signal, which increases the useful life of the actuators.

The paper is organized as follows: the problem statement and preliminary results are shown in Section 2. LMI-based stability analysis and delayed nonlinear controller design conditions for a class of nonlinear descriptors systems with multiple delays are given in Section 3, and additionally conditions for input constraints are also given. In Section 4, the implementation and numerical validation of the previous theoretical results are provided. Concluding remarks are stated in Section 5.

2. Problem Statement and Preliminary Results

2.1. Problem Statement

Let us consider a nonlinear descriptor system under multiple delays of the following form:where is the state vector, is the input vector, , , , and are matrix functions assumed to be smooth and bounded, are time delays, and are the initial functions, where is the Banach space of real continuous functions on the intervals with the following norm:where stands for the Euclidean norm in . It is assumed that for each initial condition , there exists a unique solution of the system; moreover, . In this work, matrix is assumed to be invertible, at least in a region including the origin; this is a common assumption when studying systems (1) derived from the EulerLagrange formalism [63, 65].

In order to obtain LMI conditions, the sector nonlinearity approach [57] is employed to compute an exact convex representation of (1). This methodology begins by defining a premise vector whose entries are different nonconstant terms in , , , ; similarly, is the premise vector with nonconstant terms in . It is assumed that each entry of the vectors and is bounded in the compact set that includes the origin, that is, and . Thus, each of them can be expressed as convex sums of their bounds:whereare scalar convex functions holding the convex sum property for all , i.e., , , , and . Then, the so-called scheduling (membership) functions can be computed:where , , and indexes are chosen as a -digit binary representation of ; similarly, , , ; and the set is a -digit binary representation of . The scheduling functions also hold the convex sum property in . Finally, an equivalent convex representation of (1) is [67]where , , , , and are constants matrices; and are the number of vertices for the right and left side of (6), respectively. It is important to notice that (6) is a convex rewriting of (1); thus, all the conclusions derived from the former directly apply to the latter.

2.2. Notation and Properties

In the following, convex sums of matrices will be shortly represented by

Thus, (6) is expressed as . Additionally, an asterisk will be employed in matrix expressions to denote the transpose of the symmetric element; for in-line ones, it indicates the transpose of the terms in its left-hand side, that is, .

Usually, when deriving LMI conditions for convex descriptor models, the designer is faced to inequalities of the form ; the scheduling functions are dropped off by means of the following relaxation lemma:

Lemma 1 (see [68]). Let , , and be matrices of adequate sizes. Then,

holds if the following LMIs,are satisfied for all , .

The following results establish the exponential estimates for time-delay nonlinear systems:

Lemma 2 (see [30]). Consider system (1). If there exists a functional and positive constants , , and , such that(1)(2),then, the solutions of the system (1) satisfy the exponential estimates:

As customary, for the analysis and design of convex descriptor models, the so-called descriptor redundancy is employed [69]; in our case, the augmented vectors and , , are employed to rewrite (6) as follows:with

In what follows, the stability and stabilization conditions are derived from the augmented system (11); nevertheless, it is important to stress that the system under study has the form (1).

Let us recall previous works on the subject. The work [66] studies the stability and stabilization of a nonlinear descriptor system with only one delay, that is (1) with . For stabilization purposes, the following control law is proposed:with and ; it is a nonlinear control law with nonlinearities of both sides of the nonlinear descriptor model. In the Section 3, a generalization of this controller will be presented.

3. LMI Conditions for Descriptor Systems with Multiple Delays

In this section, the developments are based on the following LyapunovKrasovskii functional candidate:with

Note that, the functional (14) is a valid L-K functional candidate as it reduces towhich is clearly positive definite and bounded by , with and , thus fulfilling conditions (1) in Lemma 2.

For the analysis of system (1), i.e., when , we have the following result.

Theorem 1. The origin of system (1), with and an exact convex representation (6), is exponentially stable if the exist matrices , , , with , , and a scalar such that LMIs (9) hold with

where and . Additionally, the solution of (1) satisfies the exponential estimates:with and .

Proof. The time derivative of (14) along the trajectories of (11) is withwhich once the dynamics of (11) are substituted and using Leibniz’s rule yieldFrom the latter, we have that is equal toThus, in order to fulfill condition (2) in Lemma 2, can be established after some algebraic manipulations when developing vectors and , by the following inequality:Finally, in order to drop off the scheduling functions and via the relaxation scheme in Lemma 1, the proof is concluded.
Let us consider system (1). Now, the task is to design a multiple delayed PDC control law of the following form:where and are nonlinear gains to be designed via the augmented system (11); thus, (23) can be expressed as with and . The following result provides LMI conditions for the design of the control law (23). It is based on a slightly modification of the L-K functional candidate (14), that is,with

Theorem 2. The origin of system (1) with an exact convex representation (6), under the law of control (23), is exponentially stable if existing matrices , , , , , , and a scalar if the LMIs (9) hold with the following:

Then, the vertex control gains are computed as and . Moreover, and the solution satisfies the following exponential estimates:where and .

Proof. Using the augmented system (11) and its corresponding control law, the closed-loop system isSimilar to the proof of Theorem 1, consider the functional (24) and its time derivative withSubstituting the dynamics of (28) in while using Leibniz’s rule in , we haveTherefore, is equivalent toAfter some simplifications, (condition (2) in Lemma 2 holds ifholds too; however, from the previous inequality, one cannot directly obtain LMI conditions. Thus, by means of the congruence property, that is, pre- and postmultiplying by the matrix giveswith the definitions and . Now, the previous inequality can be translated into the LMI conditions in the theorem once the relaxation scheme in Lemma 1 is applied.
Recall that, by hypothesis matrix, in (1) is invertible in a region around the origin; thus, it is always possible to calculate a standard state-space form as follows:where , , and . Naturally, it is possible to obtain a convex representation of system (34), that is,with , , and , where is the number of vertex models, . Even though nonlinear systems (34) and (1) and their convex forms are equivalent, establishing exponential stability of them via LMIs may lead to different feasibility set solution. Keeping the original descriptor form (1) results in a convex representation with less vertex matrices; this, in general, yields less conservative results [64].
Thus, the following result provides stability and stabilization conditions for systems of the form (34) by means of the L-K functional (14) (for stability) and (24) (for stabilization), respectively.

Corollary 1. Stability: the origin of system (1), with and an exact convex representation (35), is exponentially stable if the exist matrices , , , and with , , and a scalar such that LMIs,hold with

Moreover, the solution satisfies the following exponential estimates:

Stabilization: the origin of system (1) with an exact convex representation (35), under the law of control , is exponentially stable if existing matrices , , , , , , and a scalar if the LMIs (36) hold with

Then, the vertex control gains are computed as and , . Then, the solution satisfies the following exponential estimates:

Proof. It follows a similar path than results in Theorems 1 and 2, respectively.
The result above employs the same L-K functional of the descriptor approach, and thus the slack matrices and are also considered into the LMI conditions. Another set of LMI conditions for establishing the exponential estimates of systems in standard form (34) and its convex representation (35) can be done via a L-K functional without slack matrices, namely,whose positiveness is inferred by , and its boundedness , with and . This is summarized in the following result.

Corollary 2. The origin of the system (1) with an exact convex representation (35) and is exponentially stable if there exists matrices , with and a scalar such that,holds for . Then, the solution satisfies the following exponential estimates:with and .

Proof. It follows a similar path than previous results.

Remark 1. Theorem 1. establishes LMI conditions for the exponential estimates for the origin of system (1), these conditions include results in [66, Theorem 1] are always included when (system (1) with only one delay). Moreover, conditions in Corollary 2 always include those for linear systems in [4, Theorem 2], to see this set .

Remark 2. The numerical complexity of the LMI problems in the above results can be approximated by , where is the number of total LMI rows and is the number of scalar decision variables [70]. For Theorem 1 we have and ; as for Theorem 2, and ; as for standard systems, Corollary 1 (stability) is and ; Corollary 1 (stabilization) is and ; while Corollary 2 is and .
Results in Theorem 2 can be directly applied for real-world setups; nevertheless, the LMIs might render controller gains whose magnitude cannot be applied in practice. To alleviate this issue as well as to avoid damages in the actuators, the following result provides conditions for bounding the control input (23); they can be combined with those of Theorem 2.

Theorem 3. Consider the delayed nonlinear controller given in (23); then, this controller satisfies that , for any , if the following inequalities hold:

Proof. Observe that with the delayed PDC (23) yieldswhich is satisfied if the following holds:On the other hand, let us consider the following inequality on functional (14):Now, combining (48) and (49), it follows thator equivalentlyThus, (51) is satisfied ifhold too. From the latter inequalities and using the Schur complement together with congruence property with , conditions (44) and (45) follow with , , and . Once again, employing the Schur complement on (52) giveswhich yields (46) after the substitution of .

4. Examples

Next, a numerical example as well as the well-known inverted pendulum on a car is employed in order to illustrate the effectiveness of the proposed results. The LMI conditions have been checked with the LMIToolbox [6] within MATLAB 2109a.

4.1. Systems in Standard Form versus Descriptor Form

The following numerical example illustrates the advantages of the descriptor structure over standard state-space representations. Firstly, it compares stability at the origin via the LMIs in Theorem 1 and those in Corollary 2 by means of their feasibility sets. Secondly, a delayed nonlinear control law is designed via Theorem 2.

Consider a nonlinear system with two delays in the descriptor form (1):where the time delays are y, and matrices are as follows:

Descriptor form: following the sector nonlinearity approach, the following nonlinear terms have been identified: , , and ; their bounds have been calculated in the region . Thus the vertex matrices are

Commonly, systems of the form (55) are analyzed in the following standard form:where

In this case, four nonlinear terms are defined as follows: , , , and ; their bounds hold within . Some of the vertex matrices are given below:

Note that for this example, we have the following:(i)With respect to its exact convex representation, the descriptor one has , i.e., 8 vertexes, while the standard has , i.e., vertexes. Additionally, the descriptor keeps a constant matrix ; thus, descriptor form requires less computational resources [64]. Indeed, in this case with , , , , , and , the computational complexity using Theorem 1 is 7.2567 with 35 LMIs, using Corollary 1 is 9.7228 with 35 LMIs, and using Corollary 2 is 4.8713 with 19 LMIs.(ii)With respect to feasibility sets for Theorem 1, and Corollaries 1 and 2, when and parameter values as , for convex representations of systems (55) and (58), respectively. In Figure 1, the regions marked with a circle correspond to the feasibility sets using Theorem 1, while the regions marked by and are the feasibility sets obtained using Corollaries 1 and 2, respectively. It can be seen that by using the descriptor form the feasibility set is larger, i.e., the descriptor approach provides more relaxed results than the standard approach.

Thus, results given in Theorem 1 improve the previously classic results found in the literature.

On the other hand, Figure 2 shows that the system response (55) does not converge to the trivial equilibrium point when , , , , , and , .

Next, a delayed nonlinear controller of the form (23) is given bywhich is employed to stabilize this system. To this end, LMI conditions in Theorem 2 with an exponential decay together with those from in Theorem 3 for render feasible solution providing the following values:

Effectively, system (55) is stabilized at the origin as shown in Figure 3. The evolution in time of the system state is shown in Figure 3(a), while the guaranteed exponential decay in the system response is depicted in Figure 3(b).

4.2. Nonlinear Controller versus Delayed Nonlinear Controller

The following example is to illustrate the advantages of the use of artificial delays in controllers when there is the presence of noise, as mentioned in the introduction and its corroboration by various results found in the literature.

Consider the system known as the car-pendulum, whose scheme is shown in Figure 4, a mathematical model is given by

For illustrative purposes, a delayed measurement in the positions of the car and pendulum are intentionally added. Also, we define , , , , and . Thus, for , , , , , , , and , the system (63) can be rewritten aswhere

Considering the region , the nonconstant terms and their bounds are , , and ; the vertex matrices are

As mentioned above, with the purpose of showing the advantages of using a controller with delayed action in the presence of noise; for this example, two controllers are used: a nonlinear controller of the following form isand a delayed nonlinear controller of the form (23) is given by

To illustrate the effectiveness of Theorem 3, the controllers (68) are conditioned to satisfy that . For the controller (67), the same condition is requested, for which Theorem 3 can be used, after some simple adjustments when considering free-delays controller. For the controller (67), the corresponding LMI conditions in Theorem 2 with an exponential decay are found feasible with the following values:

On the other hand, for the delayed controller (68), LMI conditions in Theorem 2 with an exponential decay are also feasible with the following values:

To simulate the noise present in the sensors, a random signal with variance and a step is introduced at the system input. In Figure 5, the applied control signals are plotted, and it can be seen that the control signal presents less noisy when using the delayed controller than the other.

In Figure 6, the evolution in time of the state of the system (64), with , , under the control laws (67) and (68), is shown. It can be seen that the overshoot is greater when using a controller without a delayed action.

Remark 3. Systems of the form (63) can be stabilized by a free-delay controller of the form (67). However, using delayed controllers of the form (68) or (23) to stabilize this system class may be a better option when systems have inherent noise.

5. Conclusions

In this paper, analysis and design using a convex approach for nonlinear descriptor systems with multiple delays have been presented. This analysis allows synthesizing delayed nonlinear controllers to ensure convergence of the system trajectories with a guaranteed exponential decay; moreover, conditions for bounding the control input avoid possible saturation in the actuators have been provided. It also has been shown that keeping the descriptor form increases the possibility of obtaining feasibility in the LMI conditions, unlike the use of standard forms. Also, it is observed that including deliberately delays in the controller can reduce noise in the control signal, thus avoiding mechanical wear of the actuators. As future work, an extension of the proposed results in nonlinear descriptor systems with multiple time-varying delays is in course, since it will allow the synthesis of controllers for a larger class of systems.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that there are no conflicts of interests regarding the publication of this paper.

Acknowledgments

This work was partially supported by the CONACYT scholarship 998399 (J. A. Romero-Vega).