Stability results involving time-averaged Lyapunov function derivatives

The stability results which comprise the Direct Method of Lyapunov involve the existence of auxiliary functions (Lyapunov functions) endowed with certain definiteness properties. Although the Direct Method is very general and powerful, it has some limitations: there are dynamical systems with known stability properties for which there do not exist Lyapunov functions which satisfy the hypotheses of a Lyapunov stability theorem.In the present paper we identify a scalar switched dynamical system whose equilibrium (at the origin) has known stability properties (e.g., uniform asymptotic stability) and we prove that there does not exist a Lyapunov function which satisfies any one of the Lyapunov stability theorems (e.g., the Lyapunov theorem for uniform asymptotic stability). Using this example as motivation, we establish stability results which eliminated some of the limitations of the Direct Method alluded to. These results involve time-averaged Lyapunov function derivatives (TALFD’s). We show that these results are amenable to the analysis of the same dynamical systems for which the Direct Method fails. Furthermore, and more importantly, we prove that the stability results involving TALFD’s are less conservative than the results which comprise the Direct Method (which henceforth, we refer to as the classical Lyapunov stability results).While we confine our presentation to continuous finite-dimensional dynamical systems, the results presented herein can readily be extended to arbitrary continuous dynamical systems defined on metric spaces. Furthermore, with appropriate modifications, stability results involving TALFD’s can be generalized to discontinuous dynamical systems (DDS).     

Control vector Lyapunov functions for large-scale impulsive dynamical systems

Vector Lyapunov theory has been developed to weaken the hypothesis of standard Lyapunov theory in order to enlarge the class of Lyapunov functions that can be used for analyzing system stability. In this paper, we provide generalizations to the recent extensions of vector Lyapunov theory for continuous-time systems to address stability and control design of impulsive dynamical systems via vector Lyapunov functions. Specifically, we provide a generalized comparison principle involving hybrid comparison dynamics that are dependent on the comparison system states as well as the nonlinear impulsive dynamical system states. Furthermore, we develop stability results for impulsive dynamical systems that involve vector Lyapunov functions and hybrid comparison inequalities. Based on these results, we show that partial stability for state-dependent impulsive dynamical systems can be addressed via vector Lyapunov functions. Furthermore, we extend the recently developed notion of control vector Lyapunov functions to impulsive dynamical systems. Using control vector Lyapunov functions, we construct a universal hybrid decentralized feedback stabilizer for a decentralized affine in the control nonlinear impulsive dynamical system that possesses guaranteed gain and sector margins in each decentralized input channel. These results are then used to develop hybrid decentralized controllers for large-scale impulsive dynamical systems with robustness guarantees against full modeling and input uncertainty.     

Construction of a local and global Lyapunov function for discrete dynamical systems using radial basis functions

The basin of attraction of an asymptotically stable fixed point of the discrete dynamical system given by the iteration x_n+1=g(x_n) can be determined through sublevel sets of a Lyapunov function. In Giesl [On the determination of the basin of attraction of discrete dynamical systems. J. Difference Equ. Appl. 13(6) (2007) 523–546] a Lyapunov function is constructed by approximating the solution of a difference equation using radial basis functions. However, the resulting Lyapunov function is non-local, i.e. it has no negative discrete orbital derivative in a neighborhood of the fixed point. In this paper we modify the construction method by using the Taylor polynomial and thus obtain a Lyapunov function with negative discrete orbital derivative both locally and globally.     

Disguised toric dynamical systems

We study families of polynomial dynamical systems inspired by biochemical reaction networks. We focus on complex balanced mass-action systems, which have also been called toric. They are known or conjectured to enjoy very strong dynamical properties, such as existence and uniqueness of positive steady states, local and global stability, persistence, and permanence. We consider the class of disguised toric dynamical systems, which contains toric dynamical systems, and to which all dynamical properties mentioned above extend naturally. By means of (real) algebraic geometry we show that some reaction networks have an empty toric locus or a toric locus of Lebesgue measure zero in parameter space, while their disguised toric locus is of positive measure. We also propose some algorithms one can use to detect the disguised toric locus.     

On the stability of projected dynamical systems

A class of projected dynamical systems (PDS), whose stationary points solve the corresponding variational inequality problem (VIP), was recently studied by Dupuis and Nagurney (Ref. 1). This paper initiates the study of the stability of such PDS around their stationary points and thus gives rise to the study of the dynamical stability of VIP solutions. Examples are constructed showing that such a study can be quite distinct from the classical stability study for dynamical systems (DS). We give the definition of a regular solution to a VIP and introduce the concept of a minimal face flow induced by a PDS, which is a standard DS of a lower dimension. We then show that, at the regular solutions of the VIP, the local stability of the PDS is essentially the same as that of its minimal face flow. Hence, we reduce the problem, in this case, to one of the classical stability study of DS, a more developed discipline. In a more direct way, we then establish a series of local and global stability results of the PDS, under various conditions of monotonicity.This research was supported by the National Science Foundation under Grant DMS-9024071 under the Faculty Awards for Women Program. This support is gratefully acknowledged.     

Lyapunov Exponents versus Expansivity and Sensitivity in Cellular Automata

We establish a connection between the theory of Lyapunov exponents and the properties of expansivity and sensitivity to initial conditions for a particular class of discrete time dynamical systems; cellular automata (CA). The main contribution of this paper is the proof that all expansive cellular automata have positive Lyapunov exponents for almost all the phase space configurations. In addition, we provide an elementary proof of the non-existence of expansive CA in any dimension greater than 1. In the second part of this paper we prove that expansivity in dimension greater than 1 can be recovered by restricting the phase space to asuitablesubset of the whole space. To this extent we describe a 2-dimensional CA which is expansive over adense uncountablesubset of the whole phase space. Finally, we highlight the different behavior of expansive and sensitive CA for what concerns the speed at which perturbations propagate.     

Dimension Groups for Polynomial Odometers

Here we study a class of dynamical systems we call polynomial odometers. These are adic maps on regularly structured Bratteli diagrams and include the Pascal and Stirling adic maps as examples. We describe the dimension groups associated with these systems and use this to study spaces of invariant measures. For many, but not all, the space of invariant measures is affinely homeomorphic to the space of Borel probability measures on a closed interval in $\mathbb{R}$ , we call such polynomial odometers reasonable. We describe the possible isomorphisms between dimension groups for reasonable polynomial odometers, and use this to prove a version of a result of Giordano, Putnam and Skau for this situation. Namely, we show that there is an isomorphism between unital ordered groups associated with two reasonable polynomial odometers if and only if there is a special kind of orbit equivalence between the two.     

Simple polynomial classes of chaotic jerky dynamics

Third-order explicit autonomous differential equations, commonly called jerky dynamics, constitute a powerful approach to understand the properties of functionally very simple but nonlinear three-dimensional dynamical systems that can exhibit chaotic long-time behavior. In this paper, we investigate the dynamics that can be generated by the two simplest polynomial jerky dynamics that, up to these days, are known to show chaotic behavior in some parameter range. After deriving several analytical properties of these systems, we systematically determine the dependence of the long-time dynamical behavior on the system parameters by numerical evaluation of Lyapunov spectra. Some features of the systems that are related to the dependence on initial conditions are also addressed. The observed dynamical complexity of the two systems is discussed in connection with the existence of homoclinic orbits.     

Computational bounds on polynomial differential equations

In this paper we study from a computational perspective some properties of the solutions of polynomial ordinary differential equations.We consider elementary (in the sense of Analysis) discrete-time dynamical systems satisfying certain criteria of robustness. We show that those systems can be simulated with elementary and robust continuous-time dynamical systems which can be expanded into fully polynomial ordinary differential equations in Q[π]. This sets a computational lower bound on polynomial ODEs since the former class is large enough to include the dynamics of arbitrary Turing machines.We also apply the previous methods to show that the problem of determining whether the maximal interval of definition of an initial-value problem defined with polynomial ODEs is bounded or not is in general undecidable, even if the parameters of the system are computable and comparable and if the degree of the corresponding polynomial is at most 56.Combined with earlier results on the computability of solutions of polynomial ODEs, one can conclude that there is from a computational point of view a close connection between these systems and Turing machines.     

Stability and chaos control of regularized Prabhakar fractional dynamical systems without and with delay

In this paper, we studied the stabilization of nonlinear regularized Prabhakar fractional dynamical systems without and with time delay. We establish a Lyapunov stabiliy theorem for these systems and study the asymptotic stability of these systems without design a positive definite function V (without considering the fractional derivative of function V is negative). We design a linear feedback controller to control and stabilize the nonautonomous and autonomous chaotic regularized Prabhakar fractional dynamical systems without and with time delay. By means of the Lyapunov stability, we obtain the control parameters for these type of systems. We further present a numerical method to solve and analyze regularized Prabhakar fractional systems. Furthermore, by employing numerical simulation, we reveal chaotic attractors and asymptotic stability behaviors for four systems to illustrate the presented theorem.     

Unifying theory for stability of continuous,discontinuous, and discrete-time dynamical systems

Continuous-time dynamical systems whose motions are continuous with respect to time (called continuous dynamical systems), may be viewed as special cases of continuous-time dynamical systems whose motions are not necessarily continuous with respect to time (called discontinuous dynamical systems, or DDS). We show that the classical Lyapunov stability results for continuous dynamical systems are embedded in the authors’ stability results for DDS (given in [H. Ye, A.N. Michel, L. Hou, Stability theory for hybrid dynamical systems, IEEE Transactions on Automatic Control 43 (4) (1998) 461–474]), in the following sense: if the hypotheses for a given Lyapunov stability result for continuous dynamical systems are satisfied, then the hypotheses of the corresponding stability result for DDS are also satisfied. This shows that the stability results for DDS in [H. Ye, A.N. Michel, L. Hou, Stability theory for hybrid dynamical systems, IEEE Transactions on Automatic Control 43 (4) (1998) 461–474] are much more general than was previously known, and that the quality of the DDS results therein is consistent with that of the classical Lyapunov stability results for continuous dynamical systems.By embedding discrete-time dynamical systems into a class of DDS with equivalent stability properties, we also show that when the hypotheses of the classical Lyapunov stability results for discrete-time dynamical systems are satisfied, then the hypotheses of the corresponding DDS stability results are also satisfied. This shows that the results for DDS in [H. Ye, A.N. Michel, L. Hou Stability theory for hybrid dynamical systems, IEEE Transactions on Automatic Control 43 (4) (1998) 461–474] are much more general than previously known, having connections even with discrete-time dynamical systems!Finally, we demonstrate by the means of a specific example that the stability results for DDS are less conservative than corresponding classical Lyapunov stability results for continuous dynamical systems.     

Finite-time stabilization of nonlinear impulsive dynamical systems

Finite-time stability involves dynamical systems whose trajectories converge to a Lyapunov stable equilibrium state in finite time. For continuous-time dynamical systems finite-time convergence implies nonuniqueness of system solutions in reverse time, and hence, such systems possess non-Lipschitzian dynamics. For impulsive dynamical systems, however, it may be possible to reset the system states to an equilibrium state achieving finite-time convergence without requiring non-Lipschitzian system dynamics. In this paper, we develop sufficient conditions for finite-time stability of impulsive dynamical systems using both scalar and vector Lyapunov functions. Furthermore, we design hybrid finite-time stabilizing controllers for impulsive dynamical systems that are robust against full modelling uncertainty. Finally, we present a numerical example for finite-time stabilization of large-scale impulsive dynamical systems.     

Polynomial contractions and Lyapunov regularity

The purpose of this note is twofold: to introduce the notion of polynomial contraction for a linear nonautonomous dynamics with discrete time, and to show that it persists under sufficiently small linear and nonlinear perturbations. The notion of polynomial contraction mimics the notion of exponential contraction, but with the exponential decay replaced by a polynomial decay. We show that this behavior is exhibited by a large class of dynamics, by giving necessary conditions in terms of “polynomial” Lyapunov exponents. Finally, we establish the persistence of the asymptotic stability of a polynomial contraction under sufficiently small linear and nonlinear perturbations. We also consider the case of nonuniform polynomial contractions, for which the Lyapunov stability is not uniform.     

Polynomial diffeomorphisms of C2: V. critical points and Lyapunov exponents

There is an invariant measure μ, which is the pluri-complex version of the harmonic measure of the Julia set for polynomial maps of C.In this paper we give an integral formula for the Lyapunov exponents of a polynomial automorphism with respect to μ, analogous to the Brolin-Manning formula polynomial maps of C.Our formula relates the Lyapunov exponents to the value of a Green function at a type of critical point which we define in this paper. We show that these the critical points have a natural dynamical interpretation.     

On measure-preserving functions over ?3

This paper is devoted to (discrete) p-adic dynamical systems, an important domain of algebraic and arithmetic dynamics [31]?C[41], [5]?C[8]. In this note we study properties of measurepreserving dynamical systems in the case p = 3. This case differs crucially from the case p = 2. The latter was studied in the very detail in [43]. We state results on all compatible functions which preserve measure on the space of 3-adic integers, using previous work of A. Khrennikov and author of present paper, see [24]. To illustrate one of the obtained theorems we describe conditions for the 3-adic generalized polynomial to be measure-preserving on ?₃. The generalized polynomials with integral coefficients were studied in [17, 33] and represent an important class of T-functions. In turn, it is well known that T-functions are well-used to create secure and efficient stream ciphers, pseudorandom number generators.     

基于符号数值混合计算的混成系统Lyapunov函数构造

基于平方和松弛和有理向量恢复,提出了一种符号数值混合计算方法来构造多项式Lyapunov函数以判定非线性混成系统的稳定性,首先,为Lyapunov函数预定一个给定次数的多项式模板,则Lyapunov函数构造问题可转化为相应的带参数的多项式优化问题,然后运用平方和松弛方法求得一个近似的数值多项式Lyapunov函数,再应用高斯-牛顿精化和有理向量恢复将数值多项式转化为验证的有理多项式Lyapunov函数.     

Transformation invariance of Lyapunov exponents

Lyapunov exponents represent important quantities to characterize the properties of dynamical systems. We show that the Lyapunov exponents of two different dynamical systems that can be converted to each other by a transformation of variables are identical. Moreover, we derive sufficient conditions on the transformation for this invariance property to hold. In particular, it turns out that the transformation need not necessarily be globally invertible.     

Newton-Like Dynamics and Forward-Backward Methods for Structured Monotone Inclusions in Hilbert Spaces

In a Hilbert space setting we introduce dynamical systems, which are linked to Newton and Levenberg–Marquardt methods. They are intended to solve, by splitting methods, inclusions governed by structured monotone operators M=A+B, where A is a general maximal monotone operator, and B is monotone and locally Lipschitz continuous. Based on the Minty representation of A as a Lipschitz manifold, we show that these dynamics can be formulated as differential systems, which are relevant to the Cauchy–Lipschitz theorem, and involve separately B and the resolvents of A. In the convex subdifferential case, by using Lyapunov asymptotic analysis, we prove a descent minimizing property and weak convergence to equilibria of the trajectories. Time discretization of these dynamics gives algorithms combining Newton’s method and forward-backward methods for solving structured monotone inclusions.     

Ergodic properties for some non-expanding non-reversible systems

We study several properties of invariant measures obtained from preimages, for non-invertible maps on fractal sets which model non-reversible dynamical systems. We give two ways to describe the distribution of all preimages for endomorphisms which are not necessarily expanding on a basic set Λ. We give a topological dynamics condition which guarantees that the corresponding measures converge to a unique conformal ergodic borelian measure; this helps in estimating the unstable dimension a.e. with respect to this measure with the help of Lyapunov exponents. When there exist negative Lyapunov exponents of this limit measure, we study the conditional probabilities induced on the non-uniform local stable manifolds by the limit measure, and also its pointwise dimension on stable manifolds.     

Revised CPA method to compute Lyapunov functions for nonlinear systems

The CPA method uses linear programming to compute Continuous and Piecewise Affine Lyapunov functions for nonlinear systems with asymptotically stable equilibria. In [14] it was shown that the method always succeeds in computing a CPA Lyapunov function for such a system. The size of the domain of the computed CPA Lyapunov function is only limited by the equilibrium?s basin of attraction. However, for some systems, an arbitrary small neighborhood of the equilibrium had to be excluded from the domain a priori. This is necessary, if the equilibrium is not exponentially stable, because the existence of a CPA Lyapunov function in a neighborhood of the equilibrium is equivalent to its exponential stability as shown in [11]. However, if the equilibrium is exponentially stable, then this was an artifact of the method. In this paper we overcome this artifact by developing a revised CPA method. We show that this revised method is always able to compute a CPA Lyapunov function for a system with an exponentially stable equilibrium. The only conditions on the system are that it is C²$C^2$ and autonomous. The domain of the CPA Lyapunov function can be any a priori given compact neighborhood of the equilibrium which is contained in its basin of attraction. Whereas in a previous paper [10] we have shown these results for planar systems, in this paper we cover general n-dimensional systems.