Lyapunov Stability of Quasiperiodic Systems

We present some observations on the stability and reducibility of quasiperiodic systems. In a quasiperiodic system, the periodicity of parametric excitation is incommensurate with the periodicity of certain terms multiplying the state vector. We present a Lyapunov-type approach and the Lyapunov-Floquet L-F transformation to derive the stability conditions. This approach can be utilized to investigate the robustness, stability margin, and design controller for the system.


Introduction
A large class of engineering systems, such as structures subjected to quasiperiodic excitations, is described by linear ordinary differential equations with time varying coefficients.These linear systems, in general, are described as ẏ A t y, 1.1 where A t is an n × n quasiperiodic matrix and y is an n dimensional vector.In general, it is not a trivial problem to determine if 1.1 is asymptotically stable, simply stable, or unstable.The researchers have used perturbation-type techniques or numerical approaches to investigate the stability of this system 1-3 .
In this work, we address the stability of a special class of quasiperiodic systems called as periodic quasiperiodic systems where 1.1 can be written as ẏ A 0 t A 1 t y, 1.2 where A 0 t has the principal period T and A 1 t has the period T 1 .It is noted that these periods are incommensurate.These types of equations arise in parametrically excited Micro Electro Mechanical Systems MEMS 4 .It is noted that A t has a strong parametric excitation.In this paper, we present the methodology to investigate the stability of system given by 1.2 using the L-Ftransformation and Lyapunov's method.This paper is organized as follows.In Section 2, a brief mathematical background on the L-F transformation is provided.Section 3 discusses the stability conditions followed by an example.We present the L-F-transformation-type approach for quasiperiodic systems in Section 4. The discussion and conclusions are presented in Section 5.

Floquet Theory and L-F Transformation
0, that is, the system is purely time periodic, then the State Transition Matrix STM Φ t of 1.2 can be factored as 5 where the matrix Q t is real and periodic with period 2T , R is an n × n real time invariant matrix, and I is the identity matrix.Matrix Q t is known as the L-F transformation matrix 5 .
The transformation y t Q t z t produces a real-time invariant representation of purely time periodic system 1.2 with A 1 t 0 given by ż t Az t .

2.2
It is to be noted that matrix A in 2.2 is time invariant.

Construction of Lyapunov Functions
Lyapunov's direct method is widely used in the stability analysis of general dynamical systems.It makes use of a Lyapunov function V x, t .This scalar function of the state and time may be considered as some form of time-dependent generalized energy.The basic idea of the method is to utilize the time rate of energy change in V x, t for a given system to judge whether the system is stable or not.The details about Lyapunov's method and stability theorems can be found in reference 6 .For a linear system with constant coefficients, it is rather simple to find a Lyapunov function.Consider the linear system where A is a constant matrix.A quadratic form of V x may be assumed as where P is a real, symmetric, and positive definite matrix.Then

2.6
According to the Lyapunov theorem for autonomous systems, if V x is negative definite, then the system is asymptotically stable 6 .Therefore, one can write where C is a positive definite matrix.Equation 2.7 is called the Lyapunov equation.It has been shown by Kalman and Bertram 7 that if there were eigenvalues with negative real parts asymptotically stable , then for every given positive definite matrix C, there exists a unique Lyapunov matrix P. In this study, matrix C is always taken as the identity matrix.

Stability of Quasiperiodic Systems
Consider the quasiperiodic linear differential equation given by 1.2 .In order to determine the stability bounds on A 1 t , we first use the L-F transformation y t Q t z t to 1.2 .After the L-F transformation, 1.2 can be written as It is to be noted that A is a constant matrix whose eigenvalues have negative real parts.We follow the approach presented by Infante 8 to obtain stability bounds.
Theorem 3.1 see 8 .If, for some positive definite matrix B and some ε > 0, 1 is almost surely asymptotically stable in the large, where E{•} is the expectation operator and λ max is maximum real eigenvalues of a pencil [9].
Proof.Consider the quadratic Lyapunov function V z z T Bz.Then along the trajectories of 3.1 , define It is noted that the numerator and denominator in 3.3 are quadratic forms.The pencil of quadratic forms A x, x n i,k 1 a ik x i x k and B x, x n i,k 1 b ik x i x k is a matrix-valued function defined over complex numbers λ given by A x, x λ B x, x 9 .From the properties of pencils of quadratic forms 9 , we can obtain the following inequality: where λ max is defined before and λ min is the minimum real eigenvalues of a pencil.Consider z p z T Bz 1/2 ; it can be shown 10 that z p satisfies Thus, integrating and dividing 3.5 by t, z ds provides the condition for stability 10 .The solution of 3.3 can be given as It can be observed that if E{λ t } ≤ −ε for some ε > 0, V z t is bounded and that V z t → 0 as t → ∞.This is the condition imposed by inequality given by 3.4 , which proves the results.Since y t Q t z t , the stability of 3.1 implies the stability of 1.2 .It is remarked that a necessary condition for inequality 3.4 to hold is that the eigenvalues of matrix A have negative real parts.It is also possible to obtain a result that is easier to compute but not as sharp.
Corollary 3.2.If, for some positive definite matrix B and some ε > 0, 1 is almost surely asymptotically stable in the large.
Proof.The proof follows immediately from theorem by noting that

3.9
The second inequality is obtained by performing two maximizations separately.Further, using an E{•} operator yields the desired result.It is obvious that, unless the second inequality in 3.10 is an equality, the stability results obtained will not be as good as those given by the theorem.It is noted that this theorem and corollary can be extended to study stability and robustness of a linear timeperiodic system subjected to random perturbations in a straightforward fashion, and for the details, we refer the reader to reference 11 .
Example 3.3.Consider the system ẏ A t A t y, 3.11 where 12 α is a system parameter and ω 2π.The state transition matrix STM , Φ t , when the quasiperiodic term A t 0, is given as 12 Φ t e α−1 ωt cos ωt e −ωt sin ωt −e α−1 ωt sin ωt e −ωt cos ωt Q t e Rt .

3.13
Factoring the state transition matrix as shown above, the L-F transformation matrix Q t is found as

3.14
It is noted that the system is unstable for all α > 1.Using the L-F transformation z t Q t y t c.f. 3. 14 Since B 11 > 0 for α < 1 and Therefore, B is a positive definite symmetric matrix and Lyapunov stability conditions are satisfied.
Once the B matrix is constructed, the stability theorem and the corollary can be used to determine the stability conditions for the system.Simple computations yield

3.20
Setting Det D − λI 0, the eigenvalues λ of the D matrix are computed as

3.21
Application of the theorem yields  or

3.23
Using Schwarz's Inequality 13 , E{f t } 2 ≤ E{f 2 t }, and simplification yields The results obtained from condition 3.24 for α from 0 to 1 are shown in Figure 1.
In order to get the conditions for almost sure asymptotic stability from the corollary, matrices G T t BG t B −1 and A BAB −1 are calculated as

3.25
The maximum eigenvalues of matrices given by 3.25 are computed as

Mathematical Problems in Engineering
Applying the corollary

3.28
Then using Schwarz's Inequality in 3.28 , one obtains

3.31
As expected, condition 3.31 is weaker than condition 3.24 .Figure 1 displays the result obtained from 3.31 for α in the range of 0 to 1.A comparison of conditions yielding from the theorem and corollary is shown in Figure 1.

L-F Transformation Approach for Quasiperiodic System
In the previous section, we presented the theorem and corollary that provide the bounds on the quasiperiodic term so that the system described by 1.2 is stable.Alternatively, one can use L-F transformation type approach to ascertain the stability of quasiperiodic system.Unlike the theorem and corollary, this approach does not need A to have negative real parts c.f. 3.1 .Consider a quasiperiodic system given in second order form 14 by ẍ δ α cos t cos ωt x 0, 4.1 where δ, α are constants, quasiperiodicity ω ≈ p/q where p, q are integers and ẍ d 2 x/dt If eigenvalues of the FTM are inside the unit circle, then 4.3 is asymptotically stable.If the eigenvalues are on the unit circle then the system is simply stable, and if the eigenvalues are outside the unit circle, then the system is unstable.The stability or instability of the time periodic system given by 4.3 implies the stability or instability of 4.1 .It can be noted that the Floquet theory states that the STM Φ τ of 4.3 can be partitioned as where Q τ is the time periodic L-F transformation matrix and B is the constant matrix of appropriate dimensions.The eigenvalues of B are called the Floquet exponent and govern the stability of the time periodic system given by 4.3 .For computation of the L-F transformation matrix via Chebyshev polynomials, we refer the reader to reference 18 .

Conclusion
In this paper, simple and efficient computational techniques to guarantee sufficient conditions for almost sure asymptotic stability of periodic quasiperiodic systems have been presented.First, the L-F transformation has been utilized to convert the periodic part of time-periodic system to a time-invariant form.For the linear periodic-quasiperiodic system, a theorem and related corollary have been suggested using the results previously obtained by Infante 8 .In order to apply the theorem and the corollary successfully, it is observed that the eigenvalues of matrix A, which governs the stability of the system, must have negative real parts and matrix B must be positive definite.One example is presented to show the application.Another apporach pressed here is based on the Floquet-type approach, where a quasiperiodic system is approximated as a periodic system and the Floquet theory can be applied to investigate the stability.Unlike the Infante type approch, the Floquet approach does not require eigenvalues of matrix A to have negative real parts.In certain cases, Floquet type decomposition for quasiperiodic system can be used to reduce quasiperiodic system to LTI system.It is expected that these methodology would be useful in studying stability and designing controllers for a number of MEMS, where governing differential equations have time periodic quasiperiodic coefficients.The approaches presented in this paper can be extended to study stability and robustness of a linear time-periodic system subjected to random perturbations.

Figure 1 :
Figure 1: Stability results for example 1 obtained by the Theorem and Corollary.
It can be noted that 4.3 is a time periodic system with principle period π.The stability of 4.3 is governed by the Floquet theory.It is possible to find out the State Transtion Matrix STM Φ τ at the end of the principle period also called as the Floquet Transition Matrix FTM Φ τ π numerically 14, 15 or analytically using Picard iterations 16, 17 .For the details on the analytical computation of STM using Picard iteration approximation, we refer to reference 16 .It is noted that the Picard iteration approach yields an approximate closed form symbolic expression of the STM for time periodic system.