Abstract

A nonlinear Lur’e-type plant with a sector bound nonlinearity is considered. The plant is stabilized by a discrete-time feedback signal with a nonperiodic uncertain sampling. The sampling control function is nonlinear and also obeys some sectoral constraints at discrete (sampling) times. The linear matrix inequality (LMI) conditions for the stability of the closed-loop system are obtained.

1. Introduction

Absolute stability theory of nonlinear systems with sectoral constraints goes back to works of A. I. Lur’e (see [1], some historical reviews can be found in [2, 3]). In [4] and subsequent papers the Lur’e problem was reduced to feasibility of a special system of Linear Matrix Inequalities (LMIs). Later, the advantage of the LMI approach to different problems of applied mathematics was comprehensively discussed in monograph [5] that launched a broad development of a specific computer software for exploring LMIs. In [68] the Lur’e theory was extended to multiple nonlinearities.

We will mention two basic concepts put forward by V. A. Yakubovich in 1960s–1970s. The first one is -procedure [2, 9, 10] that is especially useful when we have to deal with several nonlinearities. The second one is integral-quadratic constraint (IQC), the concept that was initially introduced by V. A. Yakubovich in connection with the study of pulse-width modulated systems (a special class of sampled-data systems) [11]. Regrettably, the last paper was never translated and is almost unknown for a non-Russian reader. Notice that the impact of the early works of V. A. Yakubovich on the modern IQC theory was recognized in the review part of a widely known paper [12]. Sectoral constraints and IQCs proved to be instrumental for stability analysis of various classes of nonlinear control systems (see, e.g., [2]).

The third type of constrains, that can be named discrete-time constraints, was put forward by A. Kh. Gelig for nonlinear sampled-data systems [13]. Unlike the usual sectoral constraints, discrete-time constraints are valid not at all times, by only at some discrete-time instants lying in the sampling interval. The exact position of these instants depends on the type of a pulse modulation. In the further development of this approach it was proposed to exploit a Lyapunov–Krasovskii functional [14], but later it was found that for the stability analysis it is more convenient to combine discrete-time constraints with IQCs. Namely, in [15] it was proposed to employ the IQC based on Wirtinger integral inequality [16]. Since that time, the Wirtinger-based IQCs were used in a great number of publications for various types of sampled-data systems [1726]. In particular, the problem of stabilization of a linear plant by means of a pulse-modulated signal was considered [21, 22]. The statements of the above works were formulated in terms of frequency-domain inequalities, but later some of them were restated in terms of LMIs [2326].

The main idea of the Gelig’s approach is a substitution of the initial train of pulses for a sequence of the average values of these pulses, with a supposition that these averages satisfy some instant constraints. The errors of such a substitution are estimated with the help of IQCs. Unlike other averaging theorems, the results of Gelig were not asymptotical, but could be used for an estimation of the sampling frequency from below. For sufficiently high sampling frequencies the Gelig-type stability conditions reduce to the conventional absolute stability criteria (the circle criterion, the Popov criterion, and some others).

The problem of stabilization of a continuous-time plant by a sampled-data signal attracted much attention last decade; see a review paper [27] where the existing modern approaches are outlined. We can distinguish two main competing methods. The input delay approach was contributed by E. Fridman [28, 29] who considered a sampled-data signal as a special case of delayed signal. If , , are sampling times, then a sampled signal can be reformulated as a delayed signal with , , then Lyapunov–Krasovskii functionals can be applied.

An alternative is the IQCs approach that is more closely related to the absolute stability theory [3033]. In [30] it was firstly proved that specially chosen IQCs give the same results as those previously obtained by the input delay method. In [31] it was demonstrated that by extending the IQC approach these results can be refined. Notice that the estimate for a -gain used in [30] can be considered as a reformulated Wirtinger inequality. In [32, 33] the stability problem was reduced to feasibility of an infinite number of LMIs with coefficients depending on time . It was shown that under certain assumptions they can be checked only at a finite number of points. Though being more laborious, this approach leads to improvements of the previous results.

The most publications on sampled-data stabilization treat the case when the plant is linear and the discrete-time control implements the zero-order hold strategy. As for nonlinear plants (with single or multiple nonlinearities), their stabilization problem by a linear zero-order hold sampled-data control was considered in [3437]. The technique used in these papers was based on the method of input delays and on constructing special Lyapunov–Krasovskii functionals. The paper [38] considered a multiple nonlinearity Lur’e system stabilized by a multirate discrete-data control. The consideration was based on the integral estimate from [39]; the result was formulated in the form of a frequency-domain inequality.

This paper aims to demonstrate how the classical absolute stability technique (Yakubovich’s -procedure, Lur’e sectoral constraints, IQCs, and Gelig’s discrete-time constraints) can be applied to stability analysis of a nonlinear Lur’e-type plant under a sampled-data nonlinear stabilizing signal. The results thus obtained are parsimonious and easily verifiable, but they are compatible with those obtained by more sophisticated mathematical methods.

The paper is organized as follows. Firstly we describe a model that consists of a Lur’e-type nonlinear plant under a sampled-data nonlinear control. The sampling is supposed to be nonuniform, with a known upper bound of dwell-times. Further, we demonstrate how the sampled-data stabilization problem for a Lur’e system can be treated in the frame of the Gelig–Yakubovich approach to the absolute stability theory. The key role is played by the Yakubovich’s -procedure that in fact provides an alternative technique to using Lyapunov–Krasovski functionals. It is proved that for a sufficiently high sampling frequency the stability criterion obtained is reduced to the circle criterion for absolute stability of continuous-time systems. Thus for high sampling frequencies the conservatism of the considered method is the same as that of the classical absolute stability criteria. Finally the main result is illustrated by an application to a simple first-order problem and to a sampled-data control of a mathematical pendulum.

2. System Model

Here we address a sampling-data counterpart of the circle stability criterion (see, e.g., [2, 40, 41]).

Consider a nonlinear Lur’e-type plant with a sectoral bound uncertaintyHere the nonlinearity describes an intrinsic nonlinear feedback; it is continuous and obeys the sectoral boundfor all real , . In other words, satisfies a quadratic constraint [2]for all . Here , are scalars, and , , are constant matrices of sizes , , , respectively. Obviously, system (1) has a zero equilibrium , whose stability can be investigated with the help of the circle criterion of absolute stability [2]. However, we will be interested in the case when the zero equilibrium is unstable, so a sampled-data external feedback is used for its stabilization.

Let plant (1) be governed by a sampled-data external signal. Assume that we have a strictly increasing sequence of sampling times with the lengths of the sampling intervals (dwell-times) estimated aswhere , are some positive constants. The sequence is uncertain, and it does not need to be periodic, and only estimates (4) matter. The ratio can be considered as an instant sampling frequency. Let plant (1) be controlled by a zero-order hold signal :Here , are vectors of sizes , , respectively, and is a nonlinear function (a modulation characteristic) that satisfies a sectoral constraintfor all real and some scalars , , hencefor all . Thus the quadratic constraint (9) holds not for all times , but only for discrete times (see [13]).

The zero-order hold signal defined by (6), (7) can be considered as a special case of a square (rectangular) pulse where , , and (with ) are the amplitude, the width and the instant frequency, respectively. The specific of the zero-order hold is that ; thus all the three impulse parameters, the amplitude, the width, and the frequency, are modulated. The average of the pulse signal considered on the th sampling interval is Hence , ; i.e., the error of a substitution of a pulse for its average is equal to zero.

Formula (7) reads that the amplitude modulation is of the first kind [20, 42], so discrete constraints can be imposed at the points . (In the theory of hybrid systems this type of modulation is termed as self-triggered control [43, 44].) Notice that for more elaborate types of pulse modulation the averaging method considers discrete constraints at some intermediate points , . Then the value can be interpreted as a signal with a deviating argument that can be not only delayed, but also advanced: with , .

Here we will demonstrate how the approach developed in Theorem 3.3 [20] can be reformulated for this special case. We will use not the frequency-domain inequalities (as in [20]), but LMIs. The result will be augmented by an additional IQC taken from [31].

3. The Main Statement

Let us make an additional assumption on the nonlinearity . Suppose that there exists a scalar such that the function is bounded for all , . The following theorem presents LMI conditions for the zero asymptotic of the closed-loop system (5), (6), (7).

Theorem 1. Consider a nonlinear system (5), (6), (7) with a nonlinearity satisfying (3) and with a sampled-data control satisfying (4) and (9). Assume that there exist a symmetric matrix and scalars , , such that the following set of matrix inequalities is feasible:withHere , denotes matrix transpose and asterisks stand for the matrix blocks symmetric with respect to the main diagonal. Inequalities (13), (14) are understood in the sense of positive and negative definiteness of quadratic forms. Then any solution of (5), (6), (7) is asymptotically zero: as and as .

Notice that if we abandon the above assumption on the function we can assert only that any solution is square integrable that is a weaker property than the zero asymptotic.

If the control gains vector is given and fixed, inequalities (14), (15) present LMIs with respect to variables , , , , . However, if is considered as a design parameter and needs to be chosen, then we get a problem with nonlinear constraints.

4. Proof of the Main Statement

Let us introduce an auxiliary functionfor . As well as , the function is discontinuous with jumps at the points , .

The proof follows the mathematical technique conventional for the absolute stability theory (see, e.g., [2]). The three types of constraints will be used. Firstly, we will use a quadratic (Lur’e-type) constraint (3) that can be rewritten asfor all . Secondly, consider a discrete-time (Gelig-type) constraint (9). With the help of (16) it can be rewritten asfor , . Recall that .

Thirdly, let us take two integral-quadratic (Yakubovich type) constraints. The first one is based on the Wirtinger integral inequality (see [15, 17, 19, 20]) for all . From (4) the last inequality implieswith .

The second IQC can be extracted from [31]. Since , , and , we get for . Thus

Following the -procedure (see [2, 9]), let us define a quadratic formwhere , , are some nonnegative parameters. From (17), (18), (20), (22) it follows thatfor any solution of (5), (6), (7) and any (see also [11]).

Let us define a quadratic Lyapunov function , where is a symmetric matrix satisfying (14). Denote . Since linear matrix inequality (14) can be rewritten in terms of quadratic forms asfor all vectors and scalars , , , where is some (sufficiently small) positive number and is the Euclidean vector norm. Along the solutions of (5), (6), (7) inequality (26) impliesfor all , . Notice the vector function is continuous in ; hence is also continuous for all . Integrating (27) and using (24) we getfor all . Since , (28) impliesFrom (29) and (4) it follows that that implies as .

Let us prove that is asymptotically zero. From (28) we get for . Since , we conclude that the sequence is bounded for . The sequence vanishes as , thus it is bounded, and the function is also bounded for . Introduce the matrix . Then the first equation (5) can be rewritten aswhere , and the function is bounded for . Integrating (31), we getfor . Hencefor . Because is bounded, estimate (33) implies that is bounded for . Thus the right-hand side of (31) is also bounded, and so the function is uniformly continuous. Applying Barbalat’s lemma (see, e.g., [8]) we conclude that as .

5. Some Remarks to Theorem 1

Remark 2. Let us multiply formulas (15) elementwise by and change the variables Using the Schur complement [5] inequality (14) can be rewritten aswith where Assume that the sampling frequency is sufficiently high; i.e. Then (35) can be reduced to . The last inequality ensures the fulfillment of the circle criterion of absolute stability for a continuous-time system with quadratic constraints (17) and for all .

Remark 3. Let the conditions of Theorem 1 be satisfied. As it was shown above, this implies inequality (26). Let us set , , in (26), where , are some numbers such thatThen from (26) we obtain withSince , we conclude that the matrix defined by (42) is Hurwitz stable for any numbers , satisfying (41). This gives necessary conditions for the fulfillment of Theorem 1.

Remark 4. Let us discuss a relation of Theorem 1 of this paper to Theorem 3.2 [38]. The difference of the problem setting in Theorem 3.2 with that considered in this paper is threefold. Firstly, Theorem 3.2 studies not zero asymptotic, but a stronger property of exponential stability. Secondly, Theorem 3.2 considers a more general case of multiple nonlinearities and a multirate control. Thirdly, Theorem 3.2 is formulated not as an LMI, but as a frequency-domain inequality. However, a reformulation of Theorem 1 to this more general case presents no problem. (Frequency-domain counterparts of the LMI problem given here can be found in Chapter 3 [20].) The more interesting is to compare integral estimates used in both works. The proof of Theorem 3.2 is based on an inequality from Lemma 1 [39] that can be written aswhere . The advantage of estimate (43) is that it is valid for any function , . However, for a special type of delay , , the estimate based on (43) can be refined. Integrating (43) with this special delay we obtainObviously, (44) is more conservative that the Wirtinger inequality (20) with the multiplier in the right-hand side.

6. Example: First-Order System

Consider a simplest first-order modelwhere is a scalar function and the nonlinearity satisfies the sectoral bounds (8) with some scalars , , .

Let us apply Theorem 1. Equation (45) can be rewritten in the form (5) with , , , , . Let us take . Then inequality (14) takes the formwhere the asterisks stand for the symmetric entries with respect to the main diagonal. Take any positive number for and Then , and inequality (46) is reduced towhere . By applying Sylvester’s criterion we conclude that (48) is satisfied provided that . The latter inequality can be rewritten asLet us estimate the conservatism of estimate (49). Introduce notation . Integrating (45) we come to a discrete-time mapMap (50) is contracting if that is guaranteed if .

The case when the function is linear was considered previously in [29, 30]. Assume , then we can apply Theorem 1 with , , then we come to the estimate . This inequality is consistent with the result obtained in [30] (with the help of IQCs), but Proposition 1 of [29] gives the more accurate estimate (with the help of Lyapunov–Krasovskii functionals).

Notice that unlike the linear case, the model considered in Theorem 1 is much more general. Observe also that from the practical standpoint the systems’ stability is not the only issue that should be taken into account. If a sampling frequency is close to the boundary of the stability region, the decay rate of solutions may decrease dramatically. On the other hand, (50) reduces to for the linear system considered. Thus if the sampling period is chosen , the zero equilibrium is attained in one iteration.

7. Numerical Example: Mathematical Pendulum

Following [34], consider an equation of a computer controlled pendulumfor . Let the sampling times satisfy the bound , . The parameters of (53) are , , . The function is bounded and satisfies a sectoral constraint (2) with , (see [34]). Assume that the function obeys (8) with , (i.e., its slope can deviate as .

Define , . Then we have where , . With the help of Remark 3 let us obtain upper bounds for feasible values of feedback gains , . Suppose that numbers , satisfy (41). Then the matrix defined by (42) must be Hurwitz stable. Since The conditions for Hurwitz stability of the matrix are

Let us apply MATLAB software with YALMIP package for interface and SeDuMi solver for semidefinite programming [45, 46]. Feasible values of , were found by a manual search within region (56), then , were fixed while was minimized. It was discovered that the conditions of Theorem 1 are fulfilled for with Thus the lower bound for the sampling frequency is Hz.

Let us compare the above result with Example 4.1 of [34], where the linear discrete-time control of (53) was treated (in our notation , so ). We are primarily interested in the case of a nonuniform sampling which was considered in Theorem 1 [34]. From that theorem it was found that the pendulum system is stable with and , . It is seen that even for this special case the lower bound for the sampling frequency provided in [34] is Hz, which is greater than ours. However, for the case of a periodic sampling and a linear discrete-time control the estimate given in [34] is better than ours: Hz. Notice that the simulation for the case of a periodic sampling and a linear discrete-time control gave the minimal sampling rate Hz (see Table 1 [34]).

8. Conclusion

The paper discusses an application of the absolute stability theory to a sampled-data stabilization of a nonlinear Lur’e-type system. The design of the stabilizing feedback is reduced to optimization problem for some system of matrix inequalities. The mathematical considerations are based on the Gelig–Yakubovich approach to the stability of sampled-data systems, including -procedure and specific integral-quadratic constraints. When the sampling frequency is sufficiently high, the main statement of this paper reduces to the circle criterion of absolute stability with two nonlinearities. Illustrative examples demonstrate sufficiently good agreement with the previously known results on linear sampled-data control.

Data Availability

All data generated or analysed during this study are included in the present article.

Conflicts of Interest

The author declares that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

Alexander N. Churilov was supported by the Russian Foundation for Basic Research, Grant 17-01-00102-a.