American Institute of Mathematical Sciences

System specific triangulations for the construction of CPA Lyapunov functions

 1 Department of Mathematics, University of Sussex, Falmer BN1 9QH, United Kingdom 2 Faculty of Physical Sciences, University of Iceland, 107 Reykjavik, Iceland

Received  September 2020 Published  December 2020

Fund Project: The research in this paper was partly supported by the Icelandic Research Fund (Ranní s) grant number 163074-052, Complete Lyapunov functions: Efficient numerical computation

Recently, a transformation of the vertices of a regular triangulation of ${\mathbb {R}}^n$ with vertices in the lattice $\mathbb{Z}^n$ was introduced, which distributes the vertices with approximate rotational symmetry properties around the origin. We prove that the simplices of the transformed triangulation are $(h, d)$-bounded, a type of non-degeneracy particularly useful in the numerical computation of Lyapunov functions for nonlinear systems using the CPA (continuous piecewise affine) method. Additionally, we discuss and give examples of how this transformed triangulation can be used together with a Lyapunov function for a linearization to compute a Lyapunov function for a nonlinear system with the CPA method using considerably fewer simplices than when using a regular triangulation.

Citation: Peter Giesl, Sigurdur Hafstein. System specific triangulations for the construction of CPA Lyapunov functions. Discrete & Continuous Dynamical Systems - B, doi: 10.3934/dcdsb.2020378
References:

show all references

References:
Left: The triangulation ${\mathcal T}^\text{ std}_{K}$, $K = 4$, with a triangle fan at the origin. Right: The transformed approximately rotationally symmetric triangulation ${\mathcal T}_{\Phi, K}$
The vertices of the triangulations of Figure 1, right, are mapped by the linear transformation ${\bf x} \mapsto P^{-\frac{1}{2}} {\bf x}$, where $P^{-\frac{1}{2}}$ is a symmetric and positive definite matrix. This triangulation is adapted to the structure of the system with a local Lyapunov function $V( {\bf x}) = {\bf x}^\text{T}P {\bf x}$
Left: CPA Lyapunov function for system (3.1) with $\alpha = 0.5$ and $\beta = -0.3$ using a rectangular grid. Right: The rectangular grid, with a triangle fan at the origin, used for the computation. Level-sets of the Lyapunov function are drawn in red on both figures
Left: CPA Lyapunov function for system (3.1) with $\alpha = 0.5$ and $\beta = -0.3$ using a transformed grid. Right: The transformed grid, with a triangle fan at the origin, used for the computation. Level-sets of the Lyapunov function are drawn in red on both figures. Note that the triangulation is much better adapted to the shape of the level-sets than when using a rectangular grid as in Figure 3
Left: CPA Lyapunov function for system (3.1) with $\alpha = 0.5$ and $\beta = -0.4$ using a rectangular grid. Right: The rectangular grid, with a triangle fan at the origin, used for the computation. Level-sets of the Lyapunov function are drawn in red on both figures
Left: CPA Lyapunov function for system (3.1) with $\alpha = 0.5$ and $\beta = -0.4$ using a transformed grid. Right: The transformed grid, with a triangle fan at the origin, used for the computation. Level-sets of the Lyapunov function are drawn in red on both figures. Note that the triangulation is much better adapted to the shape of the level-sets than when using a rectangular grid as in Figure 5. Both the area covered by the triangle fan, in both cases with $64$ triangles, as well as the area covered overall are much larger than when using the rectangular grid, see Figure 5
 [1] Jonathan DeWitt. Local Lyapunov spectrum rigidity of nilmanifold automorphisms. Journal of Modern Dynamics, 2021, 17: 65-109. doi: 10.3934/jmd.2021003 [2] Sara Munday. On the derivative of the $\alpha$-Farey-Minkowski function. Discrete & Continuous Dynamical Systems - A, 2014, 34 (2) : 709-732. doi: 10.3934/dcds.2014.34.709 [3] Ralf Hielscher, Michael Quellmalz. Reconstructing a function on the sphere from its means along vertical slices. Inverse Problems & Imaging, 2016, 10 (3) : 711-739. doi: 10.3934/ipi.2016018 [4] Charles Fulton, David Pearson, Steven Pruess. Characterization of the spectral density function for a one-sided tridiagonal Jacobi matrix operator. Conference Publications, 2013, 2013 (special) : 247-257. doi: 10.3934/proc.2013.2013.247 [5] Dayalal Suthar, Sunil Dutt Purohit, Haile Habenom, Jagdev Singh. Class of integrals and applications of fractional kinetic equation with the generalized multi-index Bessel function. Discrete & Continuous Dynamical Systems - S, 2021  doi: 10.3934/dcdss.2021019 [6] Saima Rashid, Fahd Jarad, Zakia Hammouch. Some new bounds analogous to generalized proportional fractional integral operator with respect to another function. Discrete & Continuous Dynamical Systems - S, 2021  doi: 10.3934/dcdss.2021020 [7] Yanqin Fang, Jihui Zhang. Multiplicity of solutions for the nonlinear Schrödinger-Maxwell system. Communications on Pure & Applied Analysis, 2011, 10 (4) : 1267-1279. doi: 10.3934/cpaa.2011.10.1267 [8] Deren Han, Zehui Jia, Yongzhong Song, David Z. W. Wang. An efficient projection method for nonlinear inverse problems with sparsity constraints. Inverse Problems & Imaging, 2016, 10 (3) : 689-709. doi: 10.3934/ipi.2016017 [9] Olena Naboka. On synchronization of oscillations of two coupled Berger plates with nonlinear interior damping. Communications on Pure & Applied Analysis, 2009, 8 (6) : 1933-1956. doi: 10.3934/cpaa.2009.8.1933 [10] Jiangxing Wang. Convergence analysis of an accurate and efficient method for nonlinear Maxwell's equations. Discrete & Continuous Dynamical Systems - B, 2021, 26 (5) : 2429-2440. doi: 10.3934/dcdsb.2020185 [11] Manoel J. Dos Santos, Baowei Feng, Dilberto S. Almeida Júnior, Mauro L. Santos. Global and exponential attractors for a nonlinear porous elastic system with delay term. Discrete & Continuous Dynamical Systems - B, 2021, 26 (5) : 2805-2828. doi: 10.3934/dcdsb.2020206 [12] Thierry Cazenave, Ivan Naumkin. Local smooth solutions of the nonlinear Klein-gordon equation. Discrete & Continuous Dynamical Systems - S, 2021, 14 (5) : 1649-1672. doi: 10.3934/dcdss.2020448 [13] Scipio Cuccagna, Masaya Maeda. A survey on asymptotic stability of ground states of nonlinear Schrödinger equations II. Discrete & Continuous Dynamical Systems - S, 2021, 14 (5) : 1693-1716. doi: 10.3934/dcdss.2020450 [14] Amit Goswami, Sushila Rathore, Jagdev Singh, Devendra Kumar. Analytical study of fractional nonlinear Schrödinger equation with harmonic oscillator. Discrete & Continuous Dynamical Systems - S, 2021  doi: 10.3934/dcdss.2021021 [15] Irena PawŃow, Wojciech M. Zajączkowski. Global regular solutions to three-dimensional thermo-visco-elasticity with nonlinear temperature-dependent specific heat. Communications on Pure & Applied Analysis, 2017, 16 (4) : 1331-1372. doi: 10.3934/cpaa.2017065 [16] Xiaoming Wang. Quasi-periodic solutions for a class of second order differential equations with a nonlinear damping term. Discrete & Continuous Dynamical Systems - S, 2017, 10 (3) : 543-556. doi: 10.3934/dcdss.2017027 [17] Hong Yi, Chunlai Mu, Guangyu Xu, Pan Dai. A blow-up result for the chemotaxis system with nonlinear signal production and logistic source. Discrete & Continuous Dynamical Systems - B, 2021, 26 (5) : 2537-2559. doi: 10.3934/dcdsb.2020194

2019 Impact Factor: 1.27