Stability and Limit Oscillations of a Control Event-Based Sampling Criterion

Abstract

This paper investigates the presence of limit oscillations in an adaptive sampling system. The basic sampling criterion operates in the sense that each next sampling occurs when the absolute difference of the signal amplitude with respect to its currently sampled signal equalizes a prescribed threshold amplitude. The sampling criterion is extended involving a prescribed set of amplitudes. The limit oscillations might be interpreted through the equivalence of the adaptive sampling and hold device with a nonlinear one consisting of a relay with multiple hysteresis whose parameterization is, in general, dependent on the initial conditions of the dynamic system. The performed study is performed on the time domain.

1. Introduction

Nonperiodic sampling theory opens a set of new technical possibilities compared with the classical sampling with periodic sampling period [1–17]. Those possibilities are as follows: (1)to adapt the sampled signals to get better performances [1–3, 6, 8, 9, 16]: for instance, if that signal varies rapidly, then the sampling period is made smaller and vice versa. In general, some constraints in terms of bandwidth, stability, and technical requirements related to circuitry or computing should be respected so that the sampling rate has to belong to some appropriate admissibility domain;(2)transmission errors from data to results in algebraic problems like controllability, and observability, might be reduced by a judicious selection of the sampling instants when choosing a fixed sampling period. The reason is that the condition number of the matrix defining the problem depends on such a choice and one can convert a one-parameter optimization problem (a fixed sampling period) into a multiple one (the whole set of distinct sampling periods). In particular, the smaller the condition number of the coefficient matrix is, the smaller are the relative transmission errors from the data to the results depending on each particular problem dealt with [13, 14, 16]. The technique might be used by its “ad-hoc” implementation in a great variety of problems like biology measurements, economics, control theory and engineering, [16], statistics, random sampling [18–22];(3)to improve the adaptation transients in recursive identification or adaptive control of both classical or hybrid systems by combining the estimation algorithm with the signal adaptation, [2, 3, 9, 16]. Related adaptive sampling techniques can be used in the context of expert systems to improve the performances under supervisory rules (see, e.g., [17] and references therein).

Nonperiodic sampling being updated under certain adaptive sampling laws can often be interpreted as event-driven [23, 24], since, although sampling occurs through time, most of sampling rules involve signal comparison rules related to their immediate previous sampled values or involve certain performance tests. There are a set of background interesting papers, available in the literature, in which sampling is considered either state-dependent, random, or based in stochastic considerations, in general, and used in a number of applications. See, for instance, [19–22, 25–27] and references therein. The constant difference of amplitudes sampling criterion consists of keeping constant the absolute increment of the signal being sampled inbetween each two consecutive sampling instants. The sampling criterion together with its associate sampling and zero-order-hold device is equivalent to a separate nonlinearity which is fully equivalent to a multiple relay with hysteresis (i.e., a multiple bang-bang device with hysteresis). See [1, 15, 18, 28–32] and some references there in. In particular, the sampling criterion based on constant difference of amplitudes was generalized in [1] to the use of several threshold amplitudes the initial sampling criterion proposal of [30] based on a single constant difference of amplitudes. This equivalence motivates that the discretized system exhibits some properties being commonly associated with certain nonlinear systems, like for instance, the potential existence of limit oscillations. A close nonlinear model was proposed in [21] for feedback-based stabilization by triggering the plant output samples through the crossings, with hysteresis, of the signal through its quantization levels. In [22], a close problem related to saturating quantized measurements is focused on. It is well known that limit cycles are highly unsuitable in applications where the objective is to get a zero asymptotic tracking errors. However, they are pursued as an objective for the design of oscillators is some applications as in the design of tank circuits for tuning a suited frequency in radio or TV. There are unified sampling formulations available in the background literature including the presence of sampling constraints in [33, 34] and references there in, and work is also in progress to extend results to the presence of internal delays [35].

This paper characterizes and formalizes mathematically in the time domain the above sampling criterion by extending some previous background results in [1, 18], where the study of oscillations was only approximate and made in a first-harmonic approximation in the frequency domain by using the describing function approach, while the stability properties were not investigated. The results are obtained in the time domain rather than in the frequency one. This allows not to necessarily assume in the problem statement that the linear dynamics exhibits low-pass filtering properties so as to justify the use of a first-harmonic approximation method, as it was done in [1, 15, 18] which was an important limitation in those papers. In this way, such an assumption is no longer needed in the subsequent study which is performed with an exact analytic treatment rather than involving an approximate one. Also, a set of difference amplitudes, rather than a constant fixed one, are allowed in a generalized version of that sampling criterion in order to generalize the problem and to improve its potential applicability. If there is just a single amplitude available to be used as adaptive sampling threshold, then the sampling criterion is referred to as constant amplitude difference sampling criterion (CADSC). If several amplitudes are used, then the sampling criterion is referred to as sampling-dependent amplitude difference sampling criterion (SDADSC). Note that the model obtained in [4] for nonperiodically sampled systems is basically a linear time-varying difference equation. This model is useful to describe discretized systems under varying sampling periods. The time-varying coefficients of the discrete equation depend on both the sequence of sampling periods and the continuous-time parameters. Thus, it may be applied also to the criterion of constant difference of amplitudes. However, some properties like, for instance, the ability of generating limit oscillations are not easily discovered from an earlier inspection when using such a time- varying linear equation. The analysis method for stability and limit oscillations consists basically of the following steps:(a)describe the linear uncontrolled continuous-timer system by an ordinary differential equation of nth order submitted to a piecewise constant control input which varies at a set of sampling instants with, in general, time-varying sampling periods. The “ad hoc” control device for this purpose is referred to as a sampling and hold device. The solution of such a differential system is referred to as the “output” of the system;(b)discretize the equivalent differential system of nth order at generic sampling instants. Since the input is piecewise constant with discontinuities at such time instants, the solution of the differential equation for any given initial conditions coincides with that of the discretized system at sampling instants. The feedback law for a regulator with unity feedback is introduced so that the piecewise constant feedback control takes the minus values of the output at sampling instants;(c)define the generic sampling instants as those generated by the event-driven law of constant absolute difference of amplitudes of the feedback error inbetween each two consecutive sampling instants. This is generalized for a set of prescribed amplitudes in a more general sampling criterion. The amplitude, or the set of amplitudes, parameterize the solution together with the parameters of the continuous-time differential equation. It is seen that the zero-order and hold device together with the sampling criterion is equivalent to a relay with a multiple hysteresis. This suggests that limit cycles of the solution can potentially exist;(d)limit cycles are found by investigating double points of the solution in the time domain.

The dynamic system studied in this paper is complex in the sense that a continuous-time dynamic system is controlled by a feedback law consisting of an adaptive sampling criterion which is based on the use of a set of threshold amplitudes to calculate the sequence of sampling instants. For a second-order case study given in Section 5, it is shown that the zero-order hold used for discretization plus the adaptive sampling criterion itself are jointly equivalent to a relay device with multiple hysteresis. The whole feedback type is hybrid since it consist of a continuous-time system under nonlinear feedback and, in this sense, the whole system is a complex dynamic system. The equivalent multiple relay with hysteresis nonlinearity in the feedback-loop allows to interpret the presence of sustained limit oscillations as an asymptotic solution of the state-space trajectory of the closed-loop system.

2. Some Preliminary Framework and Basic Results

Notation 2. 𝐑 is the set of real numbers, 𝐑₀₊∶={𝐑∋𝑧≥0} and 𝐑₊∶={𝐑∋𝑧>0}≡𝐑₀₊⧵{0}:(i)𝐍 the set of natural numbers, 𝐍₀=𝐍∪{0} and 𝑘∶={1,2,…,𝑘}⊂𝐍 is the set of natural numbers ranging from 1 to 𝑘;(ii)PC(𝐑₀₊,𝐑) is the set of piecewise continuous functions on 𝐑₀₊;(iii)PC^(𝑛−1)([0,𝑇_per];𝐑) is the set of real almost everywhere piecewise (𝑛−1)th continuous-time differentiable functions on the definition domain[0,𝑇_per];(iv)𝐼_𝑛 is the 𝑛th order identity matrix;(v)the disjunction logic rule (spelled “or”) and the conjunction logic rule (spelled “and”) are denoted by the symbols ∨ and ∧, respectively;(vi)the ℓ₂ (or spectral) vector norm of 𝑧∈𝐑^𝑞 is defined as ‖𝑧‖₂=√𝑧^𝑇𝑧 (with the superscript “𝑇” standing for transposition. The ℓ₂-vector norm coincides with the Froebenius or Euclidean vector norm; (vii)for a real matrix 𝑀∈𝐑^𝑝×𝑞, its ℓ₂-induced matrix norm is ‖𝑀‖₂∶=max‖𝑀𝑧‖²‖𝑧‖₂=max||𝜆_𝑖∶0<‖𝑧‖₂≤1=max√𝑧^𝑇𝑧∶‖𝑧‖₂=1𝑀^𝑇𝑀||^1/2∶𝜆_𝑖∈𝜎𝑀^𝑇𝑀;∀𝑖∈𝑛_𝜎,(2.1)where 𝜎(𝑀^𝑇𝑀) is the spectrum of the square matrix 𝑀^𝑇𝑀 consisting of 1≤𝑛_𝜎≤𝑞 distinct real eigenvalues 𝜆_𝑖; 𝑖∈𝑛_𝜎. The above positive real maximum defining the spectral ‖𝑀‖₂ will be denoted by 𝜆_max(𝑀^𝑇𝑀). If 𝑞=𝑝, then ‖𝑀‖²²=‖𝑀^𝑇𝑀‖₂=𝜆_max(𝑀^𝑇𝑀)=|𝜆_max(𝑀)|²;(viii)𝑓∈𝐶_𝑇(𝐑^𝑝×[𝑡_𝑘,𝑡_𝑘+1);𝐑^𝑞) is a testing real vector function 𝑓∶𝐑^𝑝×[𝑡_𝑘,𝑡_𝑘+1)→𝐑^𝑞 within a testing class 𝐶_𝑇 being of the form 𝑓(𝑥_𝜏,𝜏), where “𝑠” stands for cartesian product of sets, with 𝑥_𝜏 being a real 𝑝-dimensional strip on [𝑡_𝑘,𝑡_𝑘+1) where 𝑡_𝑘 and 𝑡_𝑘+1 are two consecutive sampling instants from some sampling criterion SC. Thus, 𝑓is a piecewise real vector function from 𝐑^𝑝to 𝐑^𝑞 on [𝑡_𝑘,𝑡_𝑘+1) valued at some argument vector function 𝑥∶𝐑^𝑝×[𝑡_𝑘,𝑡_𝑘+1)→𝐑^𝑞.Consider the ordinary linear time-invariant differential equation: 𝐴(𝐷)𝑦(𝑡)=𝐵(𝐷)𝑢(𝑡),𝐷^𝑖𝑦(0)=𝑦^(𝑖)(0)∈𝐑,(2.2) under a piecewise continuous control input 𝑢∈PC(𝐑₀₊,𝐑) to be specified later on, where 𝐑₀₊∶={𝐑∋𝑧≥0} and the polynomials 𝐴(𝐷) and 𝐵(𝐷) of real coefficients in the time-derivative operator 𝐷∶=𝑑/𝑑𝑡 (subject to 𝐷⁰=1 and 𝐷^𝑖=𝐷𝐷^𝑖−1; for all 𝑖∈𝐍), which is formally equivalent to the Laplace transform argument “𝑠,” are 𝐴(𝐷)=^𝑛_𝑖=0𝑎_𝑖𝐷^𝑛−𝑖,𝐵(𝐷)=^𝑚_𝑖=0𝑏_𝑖𝐷^𝑚−𝑖,(2.3) where 𝑎₀≠0, 𝑏₀≠0 and 𝑛∶=deg(𝐴(𝐷))≥𝑚∶=deg(𝐵(𝐷)) so that the transfer function 𝐺(𝑠)=𝐵(𝑠)/𝐴(𝑠) is realizable, where the Laplace argument “𝑠” is formally equivalent to the time derivative operator “𝐷=𝑑/𝑑𝑡.” It is assumed with no loss in generality that the polynomial 𝐴(𝐷) is monic, that is, 𝑎₀=1, since any other nonzero value 𝑎₀≠1 can also lead to the differential equation (2.2) after normalization by 𝑎₀ of all the remaining coefficients of 𝐴(𝐷)and 𝐵(𝐷). It is also assumed that any potential zero cancellations in those polynomials, if any, are stable. This guarantees that the state-space realization is either minimal (i.e., no such cancellations exist) or, otherwise, any existing uncontrollable/unobservable mode is stable so that it does not contribute to the asymptotic solution as time tends to infinity. It is well known that the differential equation (2.2) can be described by a 𝑛th order differential system of first-order differential equations. Through the paper, it will be assumed that the differential system (2.2)-(2.3) will be controlled by a unity feedback control using, in general, nonperiodically samples 𝑦(𝑡_𝑖); 𝑖∈𝐍₀ of the solution 𝑦(𝑡). The unity feedback control law is from a zero-order sampling and hold device by 𝑢(𝑡)=𝑒𝑡_𝑘∶=𝑦𝑡_𝑘−𝑟𝑡_𝑘;∀𝑡∈𝑡_𝑘,𝑡_𝑘+1,(2.4) where 𝑟∈PC(𝐑₀₊,𝐑) and 𝑒∈PC(𝐑₀₊,𝐑) are, the reference function and the error feedback respectively, and SSIP∶=∶=𝑡_𝑘𝑇∈_𝑘𝐑₀₊∶𝑡₀∈𝐑₀₊,∞≥𝑡≥𝑡_𝑘+1>𝑡_𝑘,∀𝑘∈ID⊂𝐍₀,∈𝐑₊∶𝑇_𝑘∶=𝑡_𝑘+1−𝑡_𝑘≤𝑇≤∞,∀𝑘∈ID⊂𝐍₀,(2.5) are, the totally ordered set of sampling instants of indicator set ID⊂𝐍₀and the associated set of sampling periods with the same indicator set, respectively. Note that 𝑢∈PC(𝐑₀₊,𝐑) and it is, in particular, piecewise constant. The following simple sampling process consistency result holds directly from (2.5).

Lemma 2.1. Assume by convention, and with no loss in generality, that the first sampling instant 𝑡₀=0 and that 𝑘∈𝐍∩ID⇒(𝑘−1)∈ID. Then, 𝑡<∧∞∃⟺𝑡_𝑘(Card(ID)=Card(SI)=Card(SP)=𝑘+1<∞)∶=max_{𝑗∈𝐼𝐷}𝑡_𝑗≤𝑡<∞∧𝑇=limsup_𝑗→∞𝑇_𝑗=𝑇_𝑘=∞,𝑡=∞⟺Card(ID)=Card(SI)=Card(SP)=ℵ₀∧¬∃𝑡_𝑘∶=max_{𝑗∈𝐼𝐷}𝑡_𝑗<∞∧𝑇<∞.(2.6)

The convention 𝑡₀=0 does not imply loss in generality and it is adopted to simplify the exposition. The convention 𝑘∈𝐍∩ID⇒(𝑘−1)∈ID means that no natural number is missed inbetween any two consecutive ones in the enumeration of the members of SI and SP. The first part of Lemma 2.1 related to 𝑡<∞ means that the sampling process stops in finite time so that there is a maximum and last finite sampling instant and a last unbounded sampling period (therefore, the sequence of sampling periods is unbounded with infinite superior limit), and also that the number of sampling instants and periods is finite. The part of Lemma 2.1 for 𝑡=∞ states that the sampling process never ends so that there are infinitely many sampling instants and periods belonging to their respective numerable sets. Therefore, the cardinal of those sets is denoted by ℵ₀ related to infinite cardinals of numerable sets while the ∞ symbol is usually applied to cardinals of nonnumerable sets of infinitely many elements.

2.1. General Sampling Criterion and a Particular Sampling Criterion of Interest

A general sampling criterion SC is defined as an iterative procedure for some given testing function of the error and/or some of its time derivatives on a next tentative sampling period:𝑡𝑡_𝑘+1∈SI=SI(SC)isgeneratedbythesamplingcriterionSC;∀𝑘∈𝐍₀if_𝑘+1∶=Argmin𝐑₀₊∋𝑡>𝑡_𝑘∶𝑓∈𝐶_𝑇𝐷sa^𝑖𝑒ti×sfie𝑡s_𝑘,S𝑡C;⊂𝑡_𝑘𝐑^𝐽×𝑡_𝑘,𝑡,some𝑖∈𝐽⊆𝑛−1∪{0};𝐑∈SI∈SI,(2.7) where 𝑡₀∈SI. Therefore, given a set of sampling instants 𝑡_𝑗∈SI for all 𝑗∈𝑘∪{0}, then 𝑡_𝑘+1∈SI or 𝑇_𝑘=∞ if the sampling criterion ends such that Card(SI)=𝑘+1<∞. Sampling criteria through testing functions have been obtained in [1–3, 9, 14]. Some of them generate sampling periods in-between consecutive sampling instants as being inversely proportional to the time derivative of the sampled function, or to a combination of consecutive time derivatives, between a maximum admissibility interval (chosen from engineering requirements a such a stability or suited bandwidth). Other types of sampling criteria are chosen through integral criteria over the current sampling period of a quadratic, or some higher even power, of the error time integral between the sampled function and its previous sampled value. A very important one is the so-called criterion of constant difference of amplitudes, firstly proposed in [30], and then generalized formally in [1], and intuitively focused on in [18], to the use of a set of amplitudes which are thresholds of the variation of the sampled signal for each next sampling process. The whole element consisting of the sampling and hold device plus the CADSC (or the more general SDADSC) is, equivalently, modelled with a multiple-hysteresis relay, [1, 15, 18]. This curious nonlinearity in the control law allows an easy interpretation about why sustained oscillations can appear even when the main forward dynamics is linear. If a tracker is being designed, then the use of multiple amplitudes as signal sampling thresholds allows to decrease the amplitude of eventual sustained oscillations and then to improve the tracking servo from a control engineering point of view. Those sampling criteria have the important property, that they are able to generate sustained oscillations of great interest in oscillator design but unsuitable in tracking control problems since a permanent error between the tracked reference and the governed output signals always exists (see also [16, 17, 23]). In particular, some more general sampling criteria are obtained in [17] which include as particular cases many of those ones existing by that date in the background literature. Some of the results in this paper apply to generic sampling criteria (2.7) irrespective of each particular SC. Other specific results are mainly concerned with a particular sampling criterion, the so-called, SDADSC [1], which is defined implicitly as follows:𝑡_𝑘+1=Argmin𝐑₀₊∋𝑡>𝑡_𝑘∶||𝑒(𝑡)−𝑒𝑡_𝑘||=𝛿_𝑘∈𝐑₊,𝑡_𝑘∈SI∈SI,(2.8) for some given sampling set ST𝛿∶={𝛿_𝑘∈𝐑₊∶0<𝛿≤𝛿_𝑘≤𝛿<∞,forall𝑘∈ID} of amplitude thresholds of the SC. Note from (2.2)–(2.4) that the solution of (2.2) is unique for each set of initial conditions from Picard-Lindelöff theorem for existence and uniqueness of systems of differential equations from continuity and complete induction arguments as follows. Provided that a unique solution exists on [𝑡₀,𝑡_𝑘) for given initial conditions 𝐷^𝑖𝑦(0)=𝑦^(𝑖)(0)∈𝐑, a continuous and time-differentiable solution also exists and it is unique on [𝑡₀,𝑡_𝑘], and since the input is piecewise constant on [𝑡_𝑘,𝑡_𝑘+1), it is also continuous and time-differentiable on [𝑡_𝑘,𝑡_𝑘+1)forall𝑡_𝑘∈SI. Furthermore, the solution is everywhere continuously time-differentiable if 𝑛−𝑚≥2. This follows from the uniqueness of the solutions of ordinary differential equations (ODE) for each given set of initial conditions. The following consistency lemma follows. If 𝛿_𝑘=𝛿∈𝐑₊forall𝑘∈𝐍₀in (2.8), then the sampling criterion becomes, in particular, the CADSC [1].

Lemma 2.2. 𝑡₀∈SI⇒𝑡_𝑘∈SIforall𝑘∈𝐍, via the sampling rule (2.7), irrespective of the sampling set of amplitudes ST𝛿.

Proof. Proceed by complete induction by assuming that 𝑡_𝑗∈SI for all 𝑗∈𝑘 so that from (2.8): 𝑡_𝑗+1=⟺Arg𝑡𝑡>_𝑗+1𝑡_𝑗∶||𝑒(𝑡)−𝑒𝑡_𝑗||=𝛿_𝑗∈𝐑₊,𝑡_𝑗∈SI∈SI,∀𝑗∈𝑘−1∪{0}⊂ID=𝛿_𝑗=Arg𝑡>𝑡_𝑗∶||𝑒(𝑡)−𝑒𝑡_𝑗||∧𝑡∈𝐑₊,∀𝑡_𝑗≤𝑡_𝑘−1∈SI∈SI;∀𝑗∈𝑘−1∪{0}⊂ID_𝑘+1⟺𝑡_𝑗+1=Arg𝑡>𝑡_𝑘∶||𝑒(𝑡)−𝑒𝑡_𝑘||=𝛿_𝑘∈𝐑₊,𝑡_𝑘∈SI∈SI=Arg𝑡>𝑡_𝑗∶||𝑒(𝑡)−𝑒𝑡_𝑗||=𝛿_𝑗∈𝐑₊,∀𝑡_𝑗≤𝑡_𝑘∈SI∈SI;∀𝑗∈𝑘∪{0}⊂ID.(2.9)

The sampling criterion (2.8) and its particular version for constant amplitude is of major theoretical interest because the study of the dynamics it generates combines properties of discrete-time systems with some properties of nonlinear systems since, in particular, limit cycles appear in the solution of (2.2)–(2.4) for both the SDADSC, [1, 18]. If 𝛿_𝑘=𝛿 is a positive real constant then the sampling criterion is the CADSC, [1, 15–18]. It has been also used in some practical applications, in particular, for tuning PID controllers, [1]. Since the generation of each next sampling period is given by an implicit function in such sampling criteria, the whole control scheme might be considered in the framework of event-driven processes. It is known that for any nonsingular real matrix𝐓∈𝐑^𝑛×𝑛, a time-differentiable real-state vector function 𝑥∶𝐑₀₊→𝐑^𝑛 satisfying 𝑥(𝑡)=𝑇(𝑦(𝑡),𝐷𝑦(𝑡),…,𝐷^𝑛−1𝑦(𝑡))^𝑇 may be defined so that (2.2)–(2.4) is equivalent to the 𝑛th order dynamic feedback system:𝑥̇(𝑡)=𝐴𝑥(𝑡)+𝑏𝑟𝑡_𝑘−𝑦𝑡_𝑘,𝑦(𝑡)=𝑐^𝑇𝑥(𝑡),∀𝑡∈𝑡_𝑘,𝑡_𝑘+1,(2.10) where 𝐴∈𝐑^𝑛×𝑛 and 𝑏,𝑐∈𝐑^𝑛 are the matrix of dynamics and the control and output vectors of the continuous-time system, which depend on the coefficients of the polynomials 𝐴(𝐷) and 𝐵(𝐷) and on the entries to the matrix 𝐓, subject to initial conditions 𝑥(0)=𝐓(𝑦(0),𝐷𝑦(0),…,𝐷^𝑛−1𝑦(0))^𝑇 at 𝑡=0. Equation (2.10) holds for any set of sampling instants SI independent of the particular sampling criterion (2.7). The solution of the first equation in (2.10) within [𝑡_𝑘,𝑡_𝑘+1) yields directly, again irrespective of SC, the following discrete-time system: 𝑥𝑡_𝑘+1==Φ_𝑘𝑇_𝑘𝑥𝑡_𝑘+Γ𝑇_𝑘𝑢𝑡_𝑘=Ψ𝑇_𝑘𝑥𝑡_𝑘+Γ𝑇_𝑘𝑟𝑡_𝑘,_𝑖=0Ψ𝑇_𝑖𝑥(0)+^𝑘_𝑘𝑖=0_𝑗=𝑖+1Ψ𝑇_𝑗Γ𝑇_𝑖𝑟𝑡_𝑖,(2.11a)𝑦𝑡_𝑘+1=𝑐^𝑇𝑥𝑡_𝑘+1,(2.11b)whereΨ𝑇_𝑘∶=Φ𝑇_𝑘−Γ𝑇_𝑘𝑐^𝑇=𝑒^𝐴𝑇𝑘𝐼^𝑛−^𝑇𝑘₀𝑒^{𝐴(𝑇𝑘−𝜏)}𝑑𝜏𝑏𝑐^𝑇(2.12a)Φ𝑇_𝑘∶=𝑒^𝐴𝑇𝑘,Γ𝑇_𝑘∶=^𝑇𝑘₀𝑒^{𝐴(𝑇𝑘−𝜏)}𝑑𝜏𝑏,(2.12b)where Φ(𝑇_𝑘) and Ψ(𝑇_𝑘) are the open-loop (i.e., control-free) and closed-loop (i.e., controlled) matrices of dynamics, respectively, and Γ(𝑇_𝑘) and 𝑐 are the control and output vectors, respectively.

2.2. Basic Stability Results

The global BIBO (bounded-input bounded-output) stability of the controlled closed-loop system is discussed provided that the uncontrolled transfer function: 𝐺(𝑠)∶=𝑐^𝑇(𝑠𝐼_𝑛−𝐴)⁻¹𝑏 is stable and it possesses a sufficiently small static gain related to the admissible variation domain of the time-varying sampling periods. In the regulation case (i.e., the case of identically zero reference signal 𝑟(𝑡)), the closed-loop system is globally asymptotically Lyapunov stable. Note that the static gain of |𝐺(𝑠)|, |𝐺(0)|=|𝑐^𝑇𝐴⁻¹𝑏|, varies linearly with |𝑏^𝑇𝑐| since det(𝐴)≠0 if 𝐴 is a stability matrix. Note that the assumption of smallness of the static gain of the open-loop transfer function is always achievable via incorporation of an amplifier of sufficiently small gain 𝐾 to the forward loop provided that such a condition is not directly satisfied by the given transfer function so that |𝐾𝑐^𝑇𝑏| is as sufficiently small as requested. It is proven in the next result that if the maximum allowable time-varying sampling period 𝑇_max increases, then the allowed |𝐺(0)| being compatible with stability decreases correspondingly. If the minimum allowable sampling period 𝑇_min increases then such a gain may increase while keeping the stability. The inequality useful for stability in Theorem 2.3 below is |𝑏^𝑇𝑐|≤(1−𝜀−𝑒^{−𝑟𝑇min})𝑟/𝑇_max(1−𝑒^{−𝑟𝑇max}), with 𝜀<1−𝑒^{−𝑟𝑇min} and max(Re𝜆_𝑖(𝐴)<−𝑟<0;𝜆_𝑖(𝐴)∈𝜎(𝐴)forall𝑖∈𝑛_𝜎𝐴). The stability abscissa of the system matrix is also relevant in the sense that the gain is allowed to increase as such an abscissa increases. The subsequent result is concerned with such considerations.

Theorem 2.3. Assume that there is an admissibility bounded interval [𝑇_min,𝑇_max] such that 𝑇_𝑘∈[𝑇_min,𝑇_max] forall𝑇_𝑘∈SP for some given sampling criterion SC. Assume also that A is a stability matrix (i.e., 𝐺(𝑠)∶=𝑐^𝑇(𝑠𝐼_𝑛−𝐴)⁻¹𝑏 is a stable transfer function). Then, if |𝑏^𝑇𝑐| is sufficiently small according to an explicit trade-off related to the size of [𝑇_min,𝑇_max] and the stability abscissa of the matrix A, then the closed-loop system is BIBO stable. Furthermore, it is globally asymptotically Lyapunov stable in the regulation case without any extra assumptions on the uncontrolled transfer function.

Proof. Direct calculations with (2.12a) and (2.12b) yield: ‖‖Ψ𝑇_𝑘‖‖₂≤‖‖𝑒^𝐴𝑇𝑘‖‖₂‖‖‖𝐼_𝑛−^𝑇𝑘₀𝑒^−𝐴𝜏𝑑𝜏𝑏𝑐^𝑇‖‖‖₂≤𝑒^{−𝑟𝑇𝑘}+1−𝑒^{−𝑟𝑇𝑘}𝑟𝑇^𝑘≤𝑒_{−𝑟𝑇min}||𝑏^𝑇𝑐||+1−𝑒^{−𝑟𝑇max}𝑟𝑇^max||𝑏^𝑇𝑐||≤1−𝜀<1(2.13) provided that |𝑏^𝑇𝑐|≤(1−𝜀−𝑒^{−𝑟𝑇min})𝑟/𝑇_max(1−𝑒^{−𝑟𝑇max}), for some prefixed real constant 𝜀∈(0,1) satisfying 0<𝜀<1−𝑒^{−𝑟𝑇min}, max(Re𝜆_𝑖(𝐴)<−𝑟<0;𝜆_𝑖(𝐴)∈𝜎(𝐴)forall𝑖∈𝑛_𝜎𝐴) (for some 1≤𝑛_𝜎𝐴≤𝑛 being the number of distinct eigenvalues of 𝐴), and 𝑇_𝑘∈⌊𝑇_min,𝑇_max⌋forall𝑇_𝑘∈SP and any given sampling criterion SC. Proceeding recursively and taking ℓ₂-vector and matrix norms in (2.11a) and (2.11b), one gets since ‖‖Γ𝑇‖‖Ψ𝑇_𝑘‖‖₂≤𝜌∶=1−𝜀<1,_𝑘‖‖₂≤1−𝑒^{−𝑟𝑇𝑘}𝑟||𝑏^𝑇𝑐||≤1−𝜀−𝑒^{−𝑟𝑇min}𝑇_max,(2.14)‖‖𝑥𝑡_𝑘+1‖‖₂≤‖‖‖‖^𝑘_𝑖=0Ψ𝑇_𝑖‖‖‖‖₂‖𝑥(0)‖₂+^𝑘_𝑖₌₀‖‖‖‖^𝑘_𝑗_=𝑖+1Ψ𝑇_𝑗‖‖‖‖₂≤𝜀_𝑘‖‖Γ𝑇_𝑖‖‖₂||𝑟𝑡_𝑖||‖𝑥(0)‖₂+1−𝜀−𝑒^{−𝑟𝑇min}𝑇_max^𝑘_𝑖=0𝜀^𝑘−𝑖max_{𝑖∈𝑘∪{0}}≤𝜀_𝑘||𝑟𝑡_𝑖||‖𝑥(0)‖₂+1−𝜀−𝑒^{−𝑟𝑇min}(1−𝜀)𝑇_max≤𝐾_𝑥max_𝑖∈𝐍0||𝑟𝑡_𝑖||‖𝑥(0)‖₂+𝐾_𝑟0𝐾_𝑟<∞,∀𝑡_𝑘+1∈𝑆𝐼,(2.15) where 𝐾_𝑟0∶=(1−𝜀−𝑒^{−𝑟𝑇min})/(1−𝜀)𝑇_max, 𝐾_𝑟∶=max_𝑖∈𝐍0(‖𝑟(𝑡_𝑖)‖₂), and then 0≤limsup_𝑘→∞‖‖𝑥𝑡_𝑘+1‖‖₂≤𝐾_𝑟0max_𝑖∈𝐍0||𝑟𝑡_𝑖||≤𝐾_𝑟0𝐾_𝑟<∞,(2.16) since 𝜀^𝑘→0 as 𝑘→∞. Equations (2.11a) and (2.11b)-(2.12a) and (2.12b) are replaced within the inter-sample time intervals with: 𝑥𝑡_𝑘+𝜏=Ψ(𝜏)𝑥𝑡_𝑘+Γ(𝜏)𝑟𝑡_𝑘,∀𝜏∈⌊0,𝑇_𝑘,(2.17)Ψ(𝜏)∶=Φ(𝜏)−Γ(𝜏)𝑐^𝑇=𝑒^𝐴𝜏𝐼^𝑛−^𝜏0𝑒^{𝐴(𝜏−𝜏′)}𝑑𝜏^Φ(𝜏)∶=𝑒_𝐴𝜏𝑏𝑐^𝑇,∀𝜏∈⌊0,𝑇^𝑘,,Γ(𝜏)∶=^𝜏0𝑒^{𝐴(𝜏−𝜏′)}𝑑𝜏^𝑏,∀𝜏∈⌊0,𝑇^𝑘,∀𝑇_𝑘∈SP,(2.18) so that, one gets from (2.17)-(2.18) by using (2.15)-(2.16): ‖‖𝑥𝑡_𝑘+𝜏‖‖₂≤≤𝐾‖Ψ_(𝜏)‖₂‖‖𝑥𝑡_𝑘‖‖₂+‖Γ(𝜏)‖₂||𝑟𝑡_𝑘||𝜀^𝑘‖𝑥(0)‖₂+𝐾_𝑟0+𝑟⁻¹||𝑏^𝑇𝑐||𝐾_𝑟<∞,∀𝜏∈⌊0,𝑇_𝑘,∀𝑇_𝑘∈SP,(2.19)0≤limsup_𝑘→∞‖‖𝑥𝑡_𝑘+𝜏‖‖₂≤𝐾_𝑟0+𝑟⁻¹||𝑏^𝑇𝑐||𝐾_𝑟<∞,∀𝜏∈⌊0,𝑇_𝑘,∀𝑇_𝑘∈SP.(2.20) Thus, the dynamic system (2.11a) and (2.11b)-(2.12a) and (2.12b) is bounded-input-bounded-output (BIBO) stable for any uniformly bounded reference 𝑟(𝑡). If the reference is identically zero (regulation), then 𝐾_𝑟=0 so that, one gets from (2.16)–(2.20), that lim_𝑘→∞𝑥(𝑡_𝑘+𝜏)=0forall𝜏∈[0,𝑇_𝑘)forall𝑇_𝑘∈SP, irrespective of the initial conditions so that the dynamic system is globally asymptotically Lyapunov stable. The proof is complete.

It is obvious that Ψ(𝑇_𝑘) is a convergent matrix (i.e., a stability matrix in the discrete sense then all its eigenvalues have modulus less than unity) under the conditions of Theorem 2.3. Note that, otherwise, the state at sampling instants would be at least critically stable and would diverge for certain bounded inputs which could be fixed by construction. Since the system is globally Lyapunov stable then it also exhibits ultimate boundedness in the usual Lyapunov sense as direct conclusion from Theorem 2.3. A more general ultimate boundedness results is now derived without invoking a sufficiently small static gain of the uncontrolled system for a string of consecutive products of the matrix Ψ(𝑇_𝑘) being convergent.

Theorem 2.4. The following properties hold:(i)assume that for a given SC and each 𝑘∈ID, there exists 1≤𝑖=𝑖(𝑘)≤𝑖<∞ such that ‖Ψ(𝑡_𝑘,𝑡_𝑘+𝑖)‖₂≤𝜀₁<1, where Ψ(𝑡_𝑘,𝑡_𝑘+𝑖)∶=∏^{𝑘+𝑖−1𝑗=𝑘}[Ψ(𝑇_𝑗)]. Thus, the system is BIBO-stable for any bounded initial state and possesses the ultimate boundedness property for any bounded reference sequence 𝑟(𝑡_𝑘)forall𝑡_𝑘∈SI;(ii)if ‖Ψ(𝑡_𝑘,𝑡_𝑘+ℓ)‖₂≤𝐾(𝑘,ℓ)≤𝐾_Ψ<∞ (which holds in particular if ‖Ψ(𝑇_𝑘)‖₂≤1forall𝑡_𝑘∈SP) and ‖∑^𝑘𝑖=0∏^{𝑘𝑗=𝑖+1}[Ψ(𝑇_𝑗)]Γ(𝑇_𝑖)𝑟(𝑡_𝑖)‖₂<∞forall𝑘∈𝐍, then the system is BIBO-stable for any bounded initial state.

Proof. For any finite 𝑘∈ID, define the nonnegative scalar function 𝑣(‖𝑥(𝑡_𝑘)‖₂)∶=𝑥^𝑇(𝑡_𝑘)𝑥(𝑡_𝑘). Then, for some finite 1≤𝑖=𝑖(𝑘)≤𝑖: 𝑣‖‖𝑥𝑡_𝑘+𝑖‖‖₂−𝑣‖‖𝑥𝑡_𝑘‖‖₂=≤⎡⎢⎢⎣𝑥^𝑇𝑡_{𝑘+𝑖𝑘+𝑖−1}𝑥𝑡_𝑘+𝑖−𝑥^𝑇𝑡_𝑘𝑥𝑡_𝑘_ℓ=𝑖⎛⎜⎜⎝‖‖‖‖^{𝑘+𝑖−1}_𝑗_=𝑘+1Ψ𝑇_𝑗‖‖‖‖₂‖‖Γ𝑇_ℓ‖‖₂||𝑟𝑇_ℓ||⎞⎟⎟⎠⎤⎥⎥⎦²+⎛⎜⎜⎝2‖‖‖‖^{𝑘+𝑖−1}_𝑗=𝑘Ψ𝑇_𝑗‖‖‖‖_{2𝑘+𝑖−1ℓ=𝑖}⎛⎜⎜⎝‖‖‖‖^{𝑘+𝑖−1}_𝑗_=𝑘+1Ψ𝑇_𝑗‖‖‖‖₂−1−𝜀₂₁‖‖Γ𝑇_ℓ‖‖₂||𝑟𝑇_ℓ||⎞⎟⎟⎠‖‖‖‖^{𝑘+𝑖−1}_𝑗=𝑘Ψ𝑇_𝑗‖‖‖‖²²‖‖𝑥𝑡_𝑘‖‖₂×‖‖𝑥𝑡_𝑘⎞⎟⎟⎠=𝑥_𝑇‖‖₂𝑣‖‖𝑥𝑡_𝑘+𝑖‖‖₂−𝑣‖‖𝑥𝑡^𝑘‖‖₂≤𝐾_2𝑘𝑟𝑡_𝑘+𝑖𝑥𝑡_𝑘+𝑖−𝑥^𝑇𝑡_𝑘𝑥𝑡_𝑘−1−𝑡_𝑘𝜀,₂₁𝑡_𝑘+𝑖+2𝐾_𝑘Ψ𝑡_𝑘,𝑡_𝑘+𝑖𝐾_𝑘𝑟𝑡_𝑘,𝑡_𝑘+𝑖𝐾^2𝑘Ψ𝑡_𝑘,𝑡_𝑘+𝑖‖‖𝑥𝑡_𝑘‖‖₂‖‖𝑥𝑡_𝑘‖‖₂≤𝐾^2𝑘𝑟+2𝐾_𝑘Ψ𝐾_𝑘𝑟−1−𝜀²¹𝐾^2𝑘Ψ𝑡_𝑘,𝑡_𝑘+𝑖‖‖𝑥𝑡_𝑘‖‖₂‖‖𝑥𝑡^𝑘‖‖₂,(2.21)where 0≤𝐾_𝑘Ψ𝑡_𝑘,𝑡_𝑘+𝑖∶=‖‖‖‖^{𝑘+𝑖−1}_𝑗=𝑘Ψ𝑇_𝑗‖‖‖‖₂0≤𝐾_𝑟𝑘≤𝐾_Ψ<∞,𝑡_𝑘,𝑡_𝑘+𝑖∶=^{𝑘+𝑖−1}_ℓ=𝑖⎛⎜⎜⎝‖‖‖‖^{𝑘+𝑖−1}_𝑗_=𝑘+1Ψ𝑇_𝑗‖‖‖‖₂‖‖Γ𝑇_ℓ‖‖₂||𝑟𝑇_ℓ||⎞⎟⎟⎠≤𝐾_𝑟<∞,(2.22) with 𝐾_Ψ∶=max_{𝑘≤𝑗≤𝑘+𝑖}(𝐾_𝑘Ψ(𝑡_𝑘,𝑡_𝑘+𝑖)∶𝑘∈ID) and 𝐾_𝑟∶=max_{𝑘≤𝑗≤𝑘+𝑖}(𝐾_𝑟𝑘(𝑡_𝑘,𝑡_𝑘+𝑖)∶𝑘∈ID) being independent of 𝑘. Now, proceed by contradiction by assuming that {‖𝑥(𝑡_ℓ)‖₂}_ℓ∈IDis unbounded. Thus, there exists a subsequence {‖𝑥(𝑡_𝑘)‖₂}_{𝑘∈ID0⊆ID}which diverges for some numerable subset ID0 of ID so that 𝑣(‖𝑥(𝑡_𝑘+𝑖)‖₂)>𝑣(‖𝑥(𝑡_𝑘)‖₂) and lim_𝑘→∞𝑣(‖𝑥(𝑡_𝑘)‖₂)=∞. Choose 𝑘∈ID0 being arbitrarily large but finite, 𝑘+𝑖(𝑘)∈ID0 so that ‖𝑥(𝑡_𝑘)‖₂≥max((2𝐾_𝑘Ψ𝐾_𝑘𝑟+𝜀₂)/(1−𝜀²¹)𝐾^2𝑘Ψ(𝑡_𝑘,𝑡_𝑘+𝑖),𝑀_𝑘) with 𝑀_𝑘 being an arbitrarily large positive real number depending on 𝑘 and 𝜀₂∈𝐑₊. This is always possible since {‖𝑥(𝑡_𝑘)‖₂}_𝑘∈ID0diverges and 0≤𝜀₁<1. From (2.21), one gets: 0<≤𝑣‖‖𝑥𝑡_𝑘+𝑖‖‖₂−𝑣‖‖𝑥𝑡_𝑘‖‖₂𝐾^2𝑘𝑟+2𝐾_𝑘Ψ𝐾_𝑘𝑟−1−𝜀²¹𝐾^2𝑘Ψ𝑡_𝑘,𝑡_𝑘+𝑖‖‖𝑥𝑡_𝑘‖‖₂‖‖𝑥𝑡^𝑘‖‖₂≤𝐾^2𝑘𝑟−𝜀₃𝑀_𝑘≤0,(2.23) if 𝑀_𝑘≥𝐾^2𝑘𝑟/𝜀₃ which leads to a contradiction. As a result, {‖𝑥(𝑡_𝑘)‖₂}_𝑘∈IDis bounded from above if there exists 1≤𝑖=𝑖(𝑘)≤𝑖<∞ such that ‖Ψ(𝑡_𝑘,𝑡_𝑘+𝑖)‖₂≤𝜀₁<1, and the reference sequence is uniformly bounded. Using a similar technique as in Theorem 2.3, it may be proved that the state is bounded within any intersample period on [𝑡_𝑘,∞), forall𝑘∈ID and finite and any sampling criterion SC (i.e., for any 𝑡_𝑘∈SI generated from any SC). Therefore, the system possesses the property of ultimate boundedness. On the other hand, the recursive equations (2.15) lead to a bounded solution sequence on any interval [0,𝑡_𝑘) of finite measure if the reference sequence is uniformly bounded since all the matrices of parameters of the discrete-time dynamic system are bounded for finite time irrespective of further considerations. Therefore, ultimate boundedness implies BIBO stability and global stability of the unforced discrete-time system. The proof of Property (i) is complete. Property (ii) follows directly by taking upper bounds via the use of norms in (2.12b).

The following result parallel to Theorem 2.4 is concerned with instability:

Theorem 2.5. Assume that for a given SC and each 𝑘∈ID, there exists ∞>𝑖=𝑖(𝑘)≥𝑗∈𝐍such that ‖Ψ(𝑡_𝑘,𝑡_𝑘+𝑖)‖₂≥𝜀₁>1. Thus, the discrete-time (2.17)-(2.18) system is unstable.

Proof. Take the set of sampling instants 𝑡_𝑗∈SI for the given SC and zero reference input. Now, take initial conditions 𝑥(𝑡_𝑘) at a finite 𝑡_𝑘∈SI which are a nonzero eigenvector of Ψ(𝑡_𝑘,𝑡_𝑘+𝑗ℓ) so that lim_ℓ→∞inf‖Ψ(𝑡_𝑘,𝑡_𝑘+𝑗ℓ)𝑥(𝑡_𝑘)‖₂≥lim_ℓ→∞inf𝜀^ℓ1‖𝑥(𝑡_𝑘)‖₂=∞. Then, the system is unstable.

Note that the stability condition in terms of the modulus of eigenvalues being less than unity is equivalent in terms of positive definiteness of the matrix of dynamics to:𝜀²⁰𝐼_𝑛≤Ψ^𝑇𝑡_𝑘,𝑡_𝑘+𝑖Ψ𝑡_𝑘,𝑡_𝑘+𝑖≤𝜀²¹𝐼_𝑛,(2.24) where 0≤𝜀₀≤𝜀₁<1, which could be alternative used in both the statement and proof of Theorem 2.4.

3. Oscillations and Periodic Oscillations

Concerning the discrete-time system (2.17)-(2.18), whose expression at sampling instants are (2.11a) and (2.11b)-(2.12a) and (2.12b) for any sampling criterion SC, the following definitions for weak and strong oscillatory solutions will apply.

Definition 3.1. The discrete-time system (2.11a) and (2.11b)-(2.12a) and (2.12b) has a weak oscillatory output solution for a given sampling criterion SC and some initial condition𝑥(0)∈𝐑^𝑛 if for any given 𝑡∈𝐑₀₊, such that 𝑦(𝑡)≠0, there exist finite real numbers 𝛼(𝑡)≥𝜀_𝛼 and 𝛽(𝑡)≥𝜀_𝛽, being in general dependent on 𝑡, for some 𝜀_𝛼𝜀_𝛽∈𝐑₊, such that sign(𝛿_𝑦(𝑡,𝑡+𝛼(𝑡))𝛿_𝑦(𝑡_𝑘,𝑡+𝛼(𝑡)+𝛽(𝑡)))≤0, where 𝛿_𝑦(𝑡,𝑡^)∶=𝑦(𝑡^)−𝑦(𝑡).

Definition 3.2. The discrete-time system (2.11a) and (2.11b)-(2.12a) and (2.12b) has a strong oscillatory output solution for some initial condition 𝑥(0)∈𝐑^𝑛 if for any given 𝑡∈𝐑₀₊, such that 𝑦(𝑡)≠0, sign(𝛿_𝑦(𝑡,𝑡+𝛼(𝑡))𝛿_𝑦(𝑡_𝑘,𝑡+𝛼(𝑡)+𝛽(𝑡)))<0 and 𝑦(𝑡+𝛼(𝑡)) and 𝑦(𝑡+𝛼(𝑡)+𝛽(𝑡)) are not both zero.

Definition 3.3. The discrete-time system (2.17)-(2.18) has a periodic weak oscillatory output solution of oscillation period 𝑇_per∈𝐑₊, for some initial condition 𝑥(0)∈𝐑^𝑛, if sign(𝛿_𝑦(𝑡,𝑡+𝑇_per/2)𝛿_𝑦(𝑡,𝑡+𝑇_per))≤0, 𝑦(𝑡,𝑡+𝑇_per)=𝑦(𝑡), forall𝑡∈𝐑₊.

Definition 3.4. The discrete-time system (2.17)-(2.18) has a periodic strong oscillatory output solution of oscillation period 𝑇_per for some initial condition 𝑥(0)∈𝐑^𝑛 if it has a periodic weak oscillatory output for such a period and, furthermore, sign(𝛿_𝑦(𝑡,𝑡+𝑇_per/2)𝛿_𝑦(𝑡,𝑡+𝑇_per))<0 if 𝑦(𝑡,𝑡+𝑇_per)=𝑦(𝑡)=0, forall𝑡∈𝐑₊.

Note that a solution may be oscillatory (Definitions 3.1-3.2) without being periodic (Definitions 3.3-3.4) when there are changes in the sign of the incremental output along intervals of finite duration. A weak oscillation compared to a strong oscillation allows positive or negative increments of the output at finite intervals always of the same sign. The above Definitions 3.1-3.4 might also be refereed to, in general, to nonsymmetric oscillations related to their deviations from zero. Note that trivial solutions, that is, those being identically zero are not periodic solutions according to the given definitions. Note also that periodic solutions can possess an oscillation period which is not the sum of any fixed set of consecutive sampling periods even for such a set obeying a rule implying some repetitive sequence of periods. It turns out that the concepts of oscillation and periodic oscillation may be extended to any of the components of the state vector. The next result establishes clear implications among Definitions 3.1-3.4.

Theorem 3.5. If an output solution is strongly oscillatory, then it is also weakly oscillatory.
If an output solution is strongly periodic oscillatory then it is also weakly periodic oscillatory.
If an output solution is weakly (strongly) periodic oscillatory, then it is also weakly (strongly) oscillatory.

Note that oscillations are not always detectable for any given sampling criterion at arbitrary sampling instants since hidden oscillations can exist which cannot be detected at sampling instants. However, sufficient conditions for existence of oscillations can be formulated at sampling instants as stated in the subsequent results, whose proofs are direct conclusions of Definitions 3.1-3.2.

Theorem 3.6. The discrete-time system (2.11a) and (2.11b)-(2.12a) and (2.12b) exhibits a weak oscillatory output at sampling instants for a given sampling criterion SC and some initial condition𝑥(0)∈𝐑^𝑛 if for any 𝑡_𝑘∈SI, such that 𝑦(𝑡_𝑘)≠0, there exist finite natural numbers 𝑘₁(𝑘) and 𝑘₂(𝑘), being in general dependent on 𝑘∈ID, such that sign(𝛿_𝑦(𝑡_𝑘,𝑡_{𝑘+𝑘1(𝑘)})𝛿_𝑦(𝑡_𝑘,𝑡_{𝑘+𝑘1(𝑘)+𝑘2(𝑘)}))≤0, where 𝛿_𝑦(𝑡_𝑘,𝑡_𝑗)∶=𝑦(𝑡_𝑗)−𝑦(𝑡_𝑘).

Theorem 3.7. The discrete-time system (2.11a) and (2.11b)-(2.12a) and (2.12b) has a strong oscillatory output at sampling instants for some initial condition 𝑥(0)∈𝐑^𝑛 if sign(𝛿_𝑦(𝑡_𝑘,𝑡_{𝑘+𝑘1(𝑘)})𝛿_𝑦(𝑡_𝑘,𝑡_{𝑘+𝑘1(𝑘)+𝑘2(𝑘)}))<0 but 𝑦(𝑡_{𝑘+𝑘1(𝑘)}) and 𝑦(𝑡_{𝑘+𝑘1(𝑘)+𝑘2(𝑘)}) are not both zero.

Remark 3.8. The existence of weak and strong oscillations under the sufficient conditions of Theorems 3.6 and 3.7, respectively, may be investigated explicitly by the use of the state evolution over a finite number of consecutive sampling instants through (2.11b) together with the output expression at sampling instants in the second formula of (2.11a).

Note that the detection of periodic oscillations involving sampling instants only is not feasible even in terms of sufficient-type conditions since the period of such oscillations is not necessarily the exact sum of a consecutive number of limit sampling periods. See, for instance, [1, 15, 16], for SDADSC and CADSC, respectively. The following result states that stable uncontrolled systems which are closed-loop stable under unity feedback, fulfil the conditions of Theorem 2.4 and which do not have stable equilibrium points exhibit oscillatory responses.

Theorem 3.9. Assume that the closed-loop discrete-time system has no stable equilibrium point, while the uncontrolled system is stable under the conditions of Theorem 2.4, and the sampling criterion also fulfils the conditions of Theorem 2.4. Then, any solution of the discrete-time closed-loop system is at least weakly oscillatory and bounded.

Proof. Since any state solution is bounded for bounded initial conditions and do not converge to a constant equilibrium point, it follows that all the state components verify the incremental changes of sign of Definition 4.1, since no one can either converge to a constant or to be unbounded.

A direct related result which follows from Theorem 2.5 is now stated by simple inspection without a formal proof.

Theorem 3.10. Assume that the closed-loop discrete-time system has no stable equilibrium point while the discrete-time system is unstable under the conditions of Theorem 2.5 for some sampling criterion SC. Then, no solution of the discrete-time closed-loop system can be bounded, while it can be weakly oscillatory and unbounded.

Limit cycles are asymptotic isolated limit periodic oscillations in certain nonlinear systems which are usually independent of the initial conditions (as, e.g., the well-known Van der Pol equation). Since some nonlinear systems can also posses oscillations which depend on initial conditions, as for instance, the also well-known Duffing equation modelling certain nonlinear strings with combined linear and cubic effects, no difference is made at the moment between both situations. More precisely, a limit cycle on a plane or a 𝑛th-dimensional manifold is a closed trajectory having the property that at least one another trajectory spirals into either as time tends to infinity (stable limit cycles or a self-sustained periodic oscillations) or as time tends to minus infinity (unstable limit cycles). It turns out that a limit cycle exist in the dynamic system of Section 2 only if for a given sampling criterion SC:lim_𝑘→∞𝐷^𝑖𝑦𝑡_𝑘+𝜏=𝐷^𝑖𝑦^∗(𝜏),∀𝑖∈𝑛−1∪{0},(3.1) for some periodic function 𝑦^∗∈PC^(𝑛−1)([0,𝑇_per];𝐑) of period 𝑇_per>0 such that:(1)all its time derivatives until order (𝑛−1) exist and are almost everywhere continuous except at the sampling instants;(2)𝐷^𝑛𝑦^∗(𝑡) exists everywhere on its definition domain, but it is not required to be continuous in-between sampling instants, so that it is not required for the limit cycle to satisfy 𝑦^∗∈PC^(𝑛)([0,𝑇_per];𝐑); (3)𝐷^𝑛𝑦^∗(0⁺)=𝐷^(𝑛)𝑦^∗(𝑇⁺_per) forall 𝜏∈[0,𝑇_per) such that 𝑡_𝑘+𝜏∉SI and for (𝑡_𝑘+𝜏)⁺=𝑡^+𝑘+𝑖 if 𝑡_𝑘+𝑖∈SI for some 𝑖∈𝐍.

4. Limit Oscillations under Sampling Criteria

Note from (2.11a) and (2.11b)-(2.12a) and (2.12b) and (2.17)-(2.18) that for any SC:𝑥𝑡_𝑘+ℓ+𝜏==ΨΨ((𝜏𝜏))Ψ𝑇_{𝑘𝑘+ℓ−1}𝑥𝑡_𝑘+Γ𝑇_𝑘𝑟𝑡_𝑘+Γ(𝜏)𝑟𝑡_𝑘+ℓ+Γ(𝜏)𝑟𝑡_{𝑖=𝑘𝑘+ℓ}Ψ𝑇_𝑖𝑥𝑡_𝑘^𝑘+ℓ−1_{𝑖=𝑘𝑘+ℓ−1𝑗=𝑖+1}Ψ𝑇_𝑗Γ𝑇_𝑖𝑟𝑡_𝑖=𝑐_𝑇𝑦𝑡_𝑘+ℓ+𝜏=𝑐^𝑇𝑥𝑡_𝑘+𝜏Ψ(𝜏)^𝑘+ℓ−1_𝑖=𝑘Ψ𝑇_𝑖𝑥𝑡_𝑘+^𝑘+ℓ−1_{𝑖=𝑘𝑘+ℓ−1𝑗=𝑖+1}Ψ𝑇_𝑗Γ𝑇_𝑖𝑟𝑡_𝑖+Γ(𝜏)𝑟𝑡_𝑘+ℓ,(4.1) for all 𝑡_𝑘∈𝑆𝐼, for all 𝜏∈⌊0,𝑇_𝑘+ℓ), for all 𝑇_𝑘∈SP, and for all 𝑘∈ID, for all ℓ∈𝐍, subject to the parameterizations (2.18) becoming (2.12a) and (2.12b) at sampling instants. If a limit oscillation exists then, one gets for (4.1):∃=lim_{SI∋𝑡𝑘+ℓ→∞}𝑥lim𝑡_𝑘+ℓ+_{SP∋𝑇𝑘+𝑖,ID∋𝑘→∞}𝜏𝑥^∗^𝑘+ℓ−1_𝑖=𝑘=lim_{SI∋𝑡𝑘→∞,𝑇𝑖∈SP}𝑇_𝑘+𝑖+𝜏Ψ(𝜏)^𝑘+ℓ−1_𝑖=𝑘=lim_{SP∋𝑇𝑘+𝑖,ID∋𝑘→∞}Ψ𝑇_𝑖𝑥𝑡_𝑘+^𝑘+ℓ−1_{𝑖=𝑘𝑘+ℓ−1𝑗=𝑖+1}Ψ𝑇_𝑗Γ𝑇_𝑖𝑟𝑡_𝑖+Γ(𝜏)𝑟𝑡_𝑘+ℓΨ(𝜏)^𝑘+ℓ−1_𝑖=𝑘Ψ𝑇_𝑖𝑥^∗𝑡_𝑘+^𝑘+ℓ−1_{𝑖=𝑘𝑘+ℓ−1𝑗=𝑖+1}Ψ𝑇_𝑗Γ𝑇_𝑖𝑟𝑡_𝑖+Γ(𝜏)𝑟𝑡_𝑘+ℓ(4.,2) for all ℓ∈𝑝 and some finite 𝑝∈𝐍₀, for some 𝜏∈𝐑₀₊. Thus,𝑥^∗𝑡_𝑖+^𝑝_𝑖₌₁=Ψ(𝜏)𝑇_𝑖+𝜏_{𝑖+𝑝−1}=𝑥^∗𝑡_𝑖_𝑗=𝑖Ψ+_{𝑖+𝑝−1}𝑇_𝑗𝑥^∗𝑡_𝑖_{𝑗=1𝑖+𝑝−1ℓ=𝑗+1}Ψ𝑇_ℓΓ𝑇𝑗𝑟𝑡_𝑗Γ(𝜏)𝑟𝑡_𝑖+𝑝,(4.3) where the following limits have to exist for 𝑇_𝑘∈SP, 𝑡_𝑘∈SI, 𝑘∈IDforall𝑖∈𝑝, that is, for the sampling periods and sampling instants with a certain repeated string sequence, where 𝑡₀ in arbitrary starting limit reference sampling instant:∃∃limlim_{ID∋𝑘→∞ID∋𝑘→∞}𝑇_𝑘+𝑖∶=𝑇_𝑖=𝑇_𝑝+𝑖,∃lim_{ID∋𝑘→∞}𝑡_𝑘+𝑖∶=𝑡_𝑖=𝑡_𝑝+𝑖,𝑟𝑡_𝑘+𝑖∶=𝑟𝑡_𝑖=𝑟𝑡_𝑝+𝑖.(4.4) Then, the period of the limit oscillation is 𝑇_per∶=∑^𝑝𝑖=1𝑇_𝑖+𝜏, some real 𝜏∈[0,𝑇₁). A similar limiting equation using (4.2) into the output equation: 𝑦(𝑡_𝑘+ℓ+𝜏)=𝑐^𝑇𝑥(𝑡_𝑘+𝜏) describes the limit oscillation in the output as 𝑡_𝑘→∞. The following four lemmas related to the necessary condition of the existence of a limit cycle independent of a particular SC follow from (4.3)-(4.4) and simple topological considerations about uniqueness of the state- trajectory solution:

Lemma 4.1. If (4.3), subject to (4.4), holds, then ∃lim_𝑡→∞𝑥(𝑡+𝜏)=𝑥^∗(𝜏)=𝑥^∗(𝜏+𝑇_per); for all 𝜏∈[0,𝑇_per) and then a limit oscillation of the state-trajectory solution exists.

Proof. Note that lim_𝑘→∞𝑥(𝑡_𝑘+𝑖)=𝑥^∗(𝑡_𝑖)=𝑥^∗(𝑡_𝑖+𝑇_per) and lim_𝑘→∞𝑥(𝑡_𝑘+𝜏)=𝑥^∗(𝜏)=𝑥^∗(𝜏+𝑇_per)for 𝑇_per∶=∑^𝑝𝑖=1𝑇_𝑖+𝜏, for some parameterizing 𝜏∈[0,𝑇₁), forall𝑖∈𝑝 that is, at a discrete set of (𝑝+1) limit sampling instants as time tends to infinity, some 𝑝∈𝐍₀ and this sequence of identities is repeated with period 𝑇_per. The statetrajectory inbetween consecutive samples is prescribed according to the values of the limit reference and the state trajectory components cannot intersect at any time so that the periodic limit identity holds in continuous-time as time tends to infinity, and the result is proven.

Lemma 4.2. Assume that distinct double points 𝑥^∗(𝑡_𝑖) (𝑖∈𝑝) exist satisfying (4.3), subject to (4.4) for some 𝑝∈𝐍₀, or equivalently, 𝐼^𝑛−Ψ(𝜏)^{𝑖+𝑝−1}_𝑗=𝑖Ψ𝑇_𝑗𝑥^∗𝑡_𝑖=^{𝑖+𝑝−1}_{𝑗=1𝑖+𝑝−1ℓ=𝑗+1}Ψ𝑇_ℓΓ𝑇𝑗𝑟𝑡_𝑗,(4.5) for all𝑖∈𝑝 and some 𝜏∈[0,𝑇₁). Assume also that if 𝑝=1 then 𝑥^∗(𝑡_𝑖) satisfying the above identity is not an equilibrium point. Then, the existing limit oscillation may be tested by any of the double points, in particular, by the limit double point 𝑥^∗(𝑡₁) satisfying: 𝐼^𝑛−Ψ(𝜏)^𝑝_𝑗=1Ψ𝑇_𝑗𝑥^∗𝑡₁=^𝑝_𝑝𝑗=1_ℓ=𝑗+1Ψ𝑇_ℓΓ𝑇𝑗𝑟𝑡_𝑗.(4.6) If the reference sequence is identically zero, then a limit oscillation exists verifying double points 𝑥∈Ker𝐼^𝑛−Ψ(𝜏)∏^𝑝𝑗=1Ψ𝑇_𝑗P_eqifKer𝐼^𝑛−Ψ(𝜏)∏^𝑝𝑗=1Ψ𝑇_𝑗P_eq≠{0}.(4.7)

Proof. It follows from (4.3)-(4.4) and Lemma 4.1 since the state-trajectory solution is unique for any initial conditions, sampling periods and reference sequence and a periodic limit oscillation exist. Since the limit double points are distinct, they are not equilibrium points since the state-trajectory solution is unique if 𝑝>1. If 𝑝=1, the double point is not an equilibrium one as a requirement of the lemma statement.

Lemma 4.3. If Lemmas 4.1-4.2 hold for a given set of 𝑝∈𝐍₀ limit sampling periods 𝑇_𝑖forall𝑖∈𝑝 and some real𝜏∈[0,𝑇₁), then there is no other limit oscillation for the same sets of limit sampling periods and limit reference sequence neighboring the one with oscillation period 𝑇_per∶=∑^𝑝𝑖=1𝑇_𝑖+𝜏.

Proof. If 𝜏→𝜏+Δ𝜏 then 𝑇_per→𝑇_per+Δ𝜏 provided identical limit sampling periods 𝑇_𝑖forall𝑖∈𝑝. Since all state-trajectories are distinct, any two closed trajectories cannot be everywhere identical. Thus, two trajectories with identical initial conditions should bifurcate to different subtrajectories to complete both distinct closed paths at points inside the common parts of both trajectories. This contradicts the fact that state-trajectory solutions are unique.

Lemma 4.4. All the closed state trajectory solutions verifying Lemmas 4.1–4.3 are either stable or any unstable one, if any, is surrounded by two stable ones, namely, point-wise strictly bounded from above and below by two distinct stable closed state-trajectory solutions. Furthermore, any two closed stable trajectories cannot be arbitrarily close to each other.

Proof. The two matrices Ψ(𝜏)(∏^𝑝𝑗=1[Ψ(𝑇_𝑗)]) have to posses at least to complex conjugate eigenvalues at the unit circumference for both tuples (𝑝_ℓ,𝜏_ℓ,𝑇_𝑗ℓ;𝑗∈𝑝₁); ℓ=1,2 associated with the limit closed state-trajectory solutions. Otherwise, the system would be either BIBO stable or unstable from Theorems 2.4-2.5. Then, since all the eigenvalues are within the closed unity circle, the system is BIBO stable from Theorem 2.4 so that any state-trajectory solution can be unbounded. Thus, all existing limit oscillations are bounded for all time and then either stable or surrounded by two stable ones. On the other hand, if any two stable trajectories are arbitrary and close to each other then it would be destroyed by any arbitrarily small disturbance so they would not be stable.

Lemma 4.4 dictates that potential limit cycles of the solutions are separated to each other so that there is no accumulation closed attractor of the state-space trajectories. The interpretation of the implications of Lemma 4.4 for a linear dynamics of dimension 𝑛=2 is direct. For 𝑛>2, it is possible to interpret the lemma consequences in a plane corresponding to a 2nd-dimensional system for two of the state components in the same above way while, for the remaining components, we can consider the surrounding trajectories being equal to that one under consideration. The whole surrounding closed trajectories are still distinct from the study for the second-order subsystem.

5. Limit Oscillations for the Constant and Sampling-Dependent Amplitude Difference Sampling Criteria

5.1. The CADSC

The presence of limit oscillations is now discussed for the CADSC from the study of oscillations for sampling criteria in the above section. The conditions of stable limit oscillation for the CADSC are from (4.4)–(4.6) for a double point 𝑥^∗(𝑡₁) to exist satisfying:𝐼^𝑛−Ψ(𝜏)^𝑝_𝑗=1Ψ𝑇_𝑗𝑥^∗𝑡₁=^𝑝_𝑝_𝑗=1_ℓ=𝑗+1Ψ𝑇_ℓΓ𝑇𝑗𝑟𝑡_𝑗,(5.1) (𝑖∈𝑝) for some 𝑝∈𝐍₀ and some real 𝜏∈[0,𝑇₁); and, furthermore,=𝑦^∗𝑐_𝑇𝑡₁,𝑦^∗𝑡₂,…,𝑦^∗𝑡_𝑝,𝑦^∗𝑡_𝑝+𝜏𝐼^𝑛,Ψ𝑇₁,…,^𝑝_𝑗=1Ψ𝑇_𝑗,Ψ(𝜏)^𝑝_𝑗=1Ψ=ℓ𝛿,(ℓ+1)𝛿,…,ℓ+𝑝₁𝑇_𝑗𝑥^∗𝑡₁ℓ+𝑝₁−1𝛿,ℓ+𝑝₁−2𝛿,…−𝑝₂−1𝛿,ℓ+𝑝₁−𝑝₂𝛿,…,(ℓ+𝑝−1)𝛿,(ℓ+𝑝−1)𝛿+𝜏(5.2) for some 𝑝₁=𝑝₁(ℓ),𝑝₂=𝑝₂(ℓ)≤𝑝∈𝐍₀ and some finite ℓ∈𝐍₀. Since the limit oscillation may be tested starting at any sampling point, it turns out that any limit oscillation verifying (5.2) for some ℓ∈𝐍 is also verified 1≤𝑗≤ℓ by redefining the integers 𝑝₁(𝑗); 𝑝₂(𝑗)≤𝑝. A brief intuitive explanation of (5.2) follows. Take any positive value of the limit oscillation for a certain ℓ∈𝐍 such that (5.1)-(5.2) hold together for some 𝑝₁(ℓ),𝑝₂(ℓ)≤𝑝∈𝐍₀ for some 𝑝∈𝐍₀ and some real 𝜏∈[0,𝑇₁). Then, a limit oscillation with the system output satisfying (5.2) starting with no loss in generality by a positive value, continuing to increase 𝑝₁ times by step-by-step positive increments 𝛿 at each sampling instant, then decreasing 𝑝₂ times, then increasing again to complete the closed trajectory. Note that for any ℓ∈𝐍₀ satisfying the given constraints, 𝑝=(𝑝−𝑝₂(ℓ))+𝑝₂(ℓ) and |𝑝−2𝑝₂(ℓ)|≤2 since the number of negative increments is some 𝑝₂(ℓ)≤𝑝∈𝐍₀, and the number of positive increments is then 𝑝−2𝑝₂(ℓ), while the absolute difference of amplitude increments is of at most two.

5.2. The SDADSC

The study of oscillations can be directly generalized to the SDADSC as follows in the subsequent technical result whose proofs is obvious from (5.1)-(5.2).

Theorem 5.1. (1) Assume that a SDADSC is defined to generate the set of sampling instants SI with a potential set of amplitudes ST𝛿∶={𝛿_𝑘∈𝐑₊∶0<𝛿≤𝛿_𝑘≤𝛿<∞,forall𝑘∈ID} so that ID has a finite cardinal cardID=𝛾₀≥1 (If 𝛾₀=1, one has the particular CADSC);
(2) consider any strictly ordered finite sequence of 𝛾≥𝛾₀ amplitudes with possible repetitionsST𝛿=ST𝛿(𝛾)∶={𝛿₁,𝛿₂,…,𝛿_𝛾∶𝛿_𝑖=𝛿_𝑗∈ST𝐷,∀𝑖∈𝛾,some𝑗∈𝛾₀};
Let a finite real number M be defined as 𝑀=𝑀(𝑗)∶=∑^ℓ𝑖=0𝛿(𝑡_𝑖), where ℓ∈𝐍₀ is some finite positive integer defined according to 𝛿(𝑡_𝑖)=𝛿_𝑘∈ST𝛿forall𝑡_𝑖∈SI and some chosen 𝑘=𝑘(𝑡_𝑖)∈𝛾₀. Also, define accordingly a set of real numbers 𝑀₁=𝑀+𝛿₁ and 𝑀_𝑖+1=𝑀_𝑖+𝛿_ìforall𝑖∈𝛾is defined from the given set ST𝛿 of amplitudes. If 𝛾>𝛾₀, then ST𝛿 contains (𝛾−𝛾₀) repeated elements.
Thus, if (5.1) is defined with 𝑝=𝑗+𝛾+1 and (5.2) holds with its right-hand side being replaced with the tuple [𝑀,𝑀₁,…,𝑀_𝛾+1+𝜏], then a limit oscillation exists which satisfies the extended version of (5.1) under the above replacements.

Note from Theorem 5.1 that by appropriate choice of the limiting sequence ST𝛿 of amplitudes, the amplitudes of sustained oscillations might be reduced compared to the use of a single amplitude.

Example 5.2. First, consider the linear dynamic system of transfer function: 𝐺(𝑠)=𝑌(𝑠)𝑈(𝑠)=𝐾𝑠𝑠+²1,(5.3) where 𝑌(𝑠) and 𝑈(𝑠) are the Laplace transforms in the Laplace argument “𝑠” of the output 𝑦(𝑡) of a linear time-invariant dynamic system and its time-differentiable control input 𝑢(𝑡) under zero initial conditions. Under linear unit control feedback 𝑢(𝑡)=−𝑦(𝑡), the closed-loop differential equation becomes 𝑦^(𝑡)+𝐾(𝑦^(𝑡)+𝑦(𝑡))=0 subject to any initial conditions 𝑦(0)=𝑦₀ and 𝑦^(0)=𝑦₀₁. Note the following features in the context of limit sustained oscillations:
(1) this system is globally asymptotically stable (then, its solution 𝑦(𝑡) is oscillation-free) for any 𝐾>0 and might describe a wide set of real processes, for instance, a mechanical system subject to damping and stiffness or the control of the angular position of a satellite with respect to its axis under a derivative or tachometric control. It can also describe mathematically a linear electric circuit with two energy storing devices specified by capacitors and/or inductors with at least one dissipative device, that is, a resistor which can be in practice either a separated dissipative device or dissipative effects of the inductors/capacitors;
(2) on the other hand, note that, in order to design an electronic oscillator, that is, an electronic system whose asymptotic solution is periodic irrespective of the initial conditions, a nonlinear effect should be included in the system. In this context, note that the solution of the above damped second-order differential equation converges asymptotically to zero for any initial conditions and then it is not periodic so that it cannot be used in that way for the design of oscillators. Note, furthermore, that a typical and well-used class of electric oscillators in applications consists of those being typically synthesized with a saturation function 𝑓(𝑢(𝑡))=sat_{𝑢𝑀0,𝑢𝑀}(𝑢(𝑡)) of a certain amplifier linear gain𝐾=𝑢_𝑀/𝑢_𝑀0 in the linear mode, of saturation threshold 𝑢_𝑀0and saturated value 𝑢_𝑀, that is, sat_{𝑢𝑀0,𝑢𝑀}(𝑢(𝑡))=𝐾𝑢(𝑡)if||𝑢(𝑡)||≤𝑢^𝑀0𝑢_𝑀,sign(𝑢(𝑡))otherwise,(5.4) together with a linear electric network being of at least third order. The oscillation condition at a frequency 𝜔₀ is that the first-harmonic of the closed-loop response of the frequency domain satisfies 1+𝐾𝐺(𝐣𝜔₀)=0, obtained under the replacement 𝑠→𝐣𝜔₀ with 𝐣=√−1, provided that such an equation has a real solution 𝜔₀>0. The amplitude 𝑦 of the such a first harmonic 𝑦(𝑡)=𝑦sin(𝜔₀𝑡+𝜑) is approximately calculated from the companion complex identity 𝐺(𝐣𝜔₀)=𝐶_sat(𝑢_𝑀0,𝑢_𝑀), where 𝐶_sat(𝑢_𝑀0,𝑢_𝑀) is the critical locus (i.e., the minus describing function −(sat𝑢(𝑡))/𝑦(𝑡), an extended concept of the frequency response for certain separable or analytical nonlinearities, [36]) which is real for the case of saturations parameterized by the pair (𝑢_𝑀0,𝑢_𝑀) so that the gain in the linear mode of the saturation is 𝐾=𝑢_𝑀/𝑢_𝑀0. Note that 1+𝐶_sat(𝑢_𝑀0,𝑢_𝑀)𝐺(𝐣𝜔₀)=0 replaces intuitively the condition 1+𝐾𝐺(𝐣𝜔₀)=0 of complex conjugate modes for the replacement (−1/𝐾)→𝐶_sat(𝑢_𝑀0,𝑢_𝑀). The precision of the computation of the locus 𝐶_sat(𝑢_𝑀0,𝑢_𝑀), and then the precision of the calculated amplitude of the oscillation first-harmonic, depends on the type of describing function calculated for the saturation. See [36] and references therein for a number of useful describing functions/critical locus for different nonlinearities through different, but mutually close, useful definitions of describing function. The temporal asymptotic solution y(𝑡) tends to the limit cycle of first-harmonic 𝑦(𝑡)=𝑦sin(𝜔₀𝑡+𝜑) for any initial conditions;
(3) it is well known that electronic oscillators with basic saturated amplifiers of gain 𝐾 (in their linear mode) require also linear network of at least third- order to be synthesized. See, for instance, [36]. This is because the impulse response hodograph 𝐺(𝐣𝜔) (being the Fourier transform, if it exists, of the impulse response of the dynamic system) of the linear feed-forward part of first- and second-order jointly stable and inversely stable systems (i.e., both poles and zeros are in Re𝑠<0) are always in the third and fourth quadrants of the complex plane. As a result, they cannot cut the critical locus of a saturation nonlinearity for some frequency since such a critical locus is always allocated in the negative real semi axis;
(4) it is now described, in the context of the current problem at hand, how sustained oscillations can be obtained from the above described CADSC and SDADSC criteria by using just second-order systems of transfer functions 𝐺(𝑠)=𝐾(𝑠+1)/𝑠² in the feedforward loop. This implies that the order of the auxiliary linear network to synthesize the oscillator can be diminished related to the typical design using electronic circuitry whose basic amplifier in saturation mode needs the use of an auxiliary network of at least third-order. Then, consider again the feedback differential equation referred to above but under discrete control at, in general, nonperiodic sampling for CADSC and SDADSC: 𝑦”(𝑡)+𝐾(𝑦′(𝑡)+1)=𝐾𝑢(𝑡),𝑢(𝑡)=−𝑦𝑡_𝑖,∀𝑡∈𝑡_𝑖,𝑡_𝑖+1,(5.5) where {𝑡_𝑖} is the real sequence of sampling instants, and {𝑇_𝑖=𝑡_𝑖+1−𝑡_𝑖}is the real sequence of sampling periods under the CADSC or the SDADSC sampling criteria forall𝑖∈𝐍₀=𝐍∪{0}. The above feedback system is a regulator since the control signal 𝑢(𝑡)=𝑒(𝑡)=−𝑦(𝑡_𝑖)forall𝑡∈[𝑡_𝑖,𝑡_𝑖+1) is generated by a zero-order hold under an identically zero external reference signal. If the control is identically zero then the resulting linear feedback is globally asymptotically Lyapunov stable to the origin which is the sole stable equilibrium point. However, the use of the sampling criteria translates into the presence of limit cycles that is, asymptotic oscillations being the limits of the solution trajectories in the phase plane. The following values are taken: 𝐾=1, 𝐴=0.2154 for the single threshold case as the sampling amplitude in the CADSC criterion and the set of amplitudes 𝑆_amp≡{𝐴₁,𝐴₂,𝐴₃,𝐴₄}={0.12,0.1677,0.2154,0.2631} for the SDADSC criterion injected in this order to implement the sampling criterion. Since the control is a regulator, r≡0 for all time, then the sampling criterion becomes ||𝑒(𝑡)−𝑒𝑡_𝑖||=||𝑦(𝑡)−𝑦𝑡_𝑖||=𝐴,∀𝑡∈𝑡_𝑖,𝑡_𝑖+1,(5.6) for the CADSC, and ||𝑒(𝑡)−𝑒𝑡_𝑖||=||𝑦(𝑡)−𝑦𝑡_𝑖||=𝐴𝑡_𝑖=𝐴_[𝑖/4],∀𝑡∈𝑡_𝑖,𝑡_𝑖+1,(5.7) where [𝑖/4]=IntegerPart(𝑖/4); forall𝑖∈𝐍₀=𝐍∪{0} for the SDADSC. Note that the sampling criteria (5.6) and (5.7) can be interpreted as a separated nonlinearity of the dynamic systems (5.5) consisting of a multi-relay with hysteresis displayed in Figure 1 where 𝑚(𝑡_𝑖)=−𝑢(𝑡_𝑖). In such a way, the sampling and hold device with the sampling criteria is equivalent to such a nonlinearity which could be potentially be generated in a completely different way by using relays with hysteresis. In the case of the SDADSC, the amplitudes are taken to vary consecutively in the defined order in the set of amplitudes of the sampling criterion. The asymptotic phase plot for the CADSC criterion and SDADSC criterion are, a limit cycle of fundamental amplitude and frequency 0.39 and 0.3079 cycles/sec., and another limit cycle of amplitude 0.31 and frequency 0.3095 cycles/sec., respectively. Both limit cycles to which the phase portraits of the trajectory solutions asymptotically converge are shown in Figure 2. A direct interpretation of why the asymptotic solution is a stable limit cyclic, so that the solution is bounded and the whole system is (nonasymptotically) stable relies on the equivalence of the tandem sampling criteria CADSC and SDADSC through a companion zero-order hold to two variants of the multiple hysteretic relay nonlinearity of Figure 1. Note that it is well known that nonlinear systems under certain conditions can generate limit cycles. The sequences of constant asymptotic sampling periods reached in both cases are also listed in the figure. Note that the fundamental amplitude of the second limit cycle corresponding to the SDADSC is reduced more that 20% with respect to the first one while the fundamental frequencies differ only in about 0.50%, the second one being very slightly larger than the first one. It has been also observed under exhaustive inspection of related examples by modifying their parameterizations that the duration of the transient time interval towards the limit cycle solution is slightly shorter under the first criterion compared to the second one.
If the sampling criterion is modified to the constant amplitude-based sampling criterion |𝑦(𝑡)+𝑦(𝑡_𝑖)|=𝐴forall𝑡∈[𝑡_𝑖,𝑡_𝑖+1) (i.e., modified CADSC) and to the multithreshold sampling criterion obtained via its right-hand-side replacement by the same multithreshold sequence as above, that is, 𝑆_amp≡{𝐴₁,𝐴₂,𝐴₃,𝐴₄}={0.12,0.1677,0.2154,0.2631} (i.e., modified SDADSC) then one gets the results of Figure 3 below. Note the complex geometry of the asymptotic oscillations of the standard criteria displayed in Figure 2 compared to the more smooth shaped ones of the modified ones displayed in Figure 3. It can be pointed out that other potential sampling criteria, not been subject to constant or varying (within a prescribed set) differences of amplitude can lead to asymptotically stable solutions provided that the admissibility domain for the sampling intervals defined by such sampling criteria is constrained to the stability domain of a constant sampling provided that the continuous part of the dynamic system is globally asymptotically stable. See, for instance, Theorem 2.3.

Figure 1 

Multirelay with hysteresis nonlinear characteristics of the sampling criterion.

Figure 2 

Phase plane plot for the solution versus its first time derivative exhibiting closed limit trajectories.

Figure 3 

Phase plane plot of the modified sampling criteria for the solution versus its first time derivative exhibiting limit cycles.

Remark 5.3. Note that Example 5.2 is based on a transfer function description of the linear part. Thus, the above mathematical results on the limit asymptotic solutions are applicable to any minimal state-space realization, since in this case the dimension of the linear system coincides with the order of the transfer function (i.e., its number of poles). In the case of nonminimal realizations (then being either noncontrollable or nonobservable or both), the above discussed results still hold if the cancelled modes are strictly stable since their contribution to the state-space trajectories and their timederivatives the relevant order vanish asymptotically as time tends to infinity.

6. Conclusions

This paper has been devoted to investigate the solutions and, in particular, their stability and instability properties as well as the possible presence of sustained oscillations in discrete linear dynamic systems under sampling laws which generate time-varying sampling periods in general. Two sampling criteria have been specially emphasized, namely, (a) the so-called constant amplitude difference sampling criterion (CADSC), under which the signal of interest is sampled at each time that it reaches a prescribed threshold variation which is the positive real constant defining the sampling criterion; (b) the more general sampling criterion is referred to as a sampling-dependent amplitude difference sampling criterion (SDADSC) which involves a set of at least two distinct of such amplitudes. Both sampling criteria possess the property that, together with their associate sampling and zero-order hold device, are characterized as a relay with multiple hysteresis. Such a nonlinear model is expected to potentially generate potentially sustained limit oscillations of the solution. The analysis has been fully performed in the time domain so that, contrarily to the case of the use of frequency-domain analysis methods, no specific assumption is needed about low-pass filtering constraints of the linear auxiliary network in order to perform the analysis of the first-harmonic of the existing sustained oscillations. It is noticed that, the proposed analysis, no separation of the first-order harmonic of the whole oscillation has to be taken in mind.

Acknowledgments

The authors thank the Spanish Ministry of Education and the Basque Government for their support of this work through Grants DPI2009-07197, IT378-10, and SAIOTEK S-PE09UN12. They also thank the editor of the special issue of “Modelling and Control of Complex Dynamic Systems: Applied Mathematics Aspects” and the anonymous referees and editor for their useful comments which helped the authors to improve the former versions of the paper.