The Hill equation is a second order differential equation of the form
where are real numbers and T is the minimal period of. Without loss of generality, it is assumed through this document that the periodic function is of zero average, i.e.
Equation (1) can always be rewritten as a two dimensional first order Linear Periodic (LP) system
As Lyapunov proved in  , the asymptotic stability of an LP system is equivalent to the asymptotic stability of a Linear Time Invariant (LTI) system under a periodic change of coordinates. Unfortunately, in order to obtain the LTI system, it is required to explicitly have the solution of the original LP system, and with some few exceptions, this is in general impossible.
Stability analysis of periodic systems is a topic of high interest in e.g. engineering. In particular, due to the existence of parametric resonance    , it becomes crucial to be able to design controllers that can overcome such effects. Some examples of physical systems modeled by a Hill equation and the effects of parametric resonance can be found in  -  . Moreover, the understanding and control of parametric resonance for electric applications is currently receiving much attention  -  . Then, techniques to eliminate (or promote) parametric resonance have immediate impact in such applications. As shown below, Arnold Tongues are a graphical way to visualize the parameter set of a Hill equation prone to such effect. Thus, developing new and fast techniques to obtain Arnold Tongues and to modify them is very relevant for the above context, and in fact those are the contributions of this paper.
In this paper, the first result shows a new method for computing Arnold Tongues which is valid for a wide range of values of parameters. It also offers infor- mation on whether the given parameters correspond to a stable or unstable solution of (1). One advantage of the presented numerical method is that it is faster compared to numerical integration. Our second contribution is the design of a vibrational controller that stabilizes an operating point of a Hill equation. The basic idea of our second contribution is presented on “On Vibrational Control and Vibrational Stabilization”  . The main difference with  is that Meerkov considers LTI systems with negative trace and introduces a periodic variation of some parameters. In this paper, in contrast, we consider an unstable linear periodic systems with no restriction on the trace of the system matrix. So, in a way, we extend the results of  .
The organization of the paper is as follows. In section 2, some preliminaries as Floquet theory, the stability of linear periodic systems and some results of linear algebra are reviewed. Next, in section 3, (1) is transformed to a Linear Shift Invariant approxi- mation of it; then, the proposed numerical method to compute the Arnold Tongues is presented and a computational time comparison is provided. In section 4, the main result of this document is presented, namely, a novel scheme to stabilize a linear periodic system is provided and afterwards, the result is illustrated with an example. Finally, in section 5 some concluding remarks as well as some discussion are provided.
The following notation is used through the paper: represents the natural numbers, given a square matrix, its determinant is denoted as, the set of its eigenvalues counting multiplicities is denoted as. The polyno- mials and represent the characteristic and minimal polynomial of respectively.
This section reviews the Floquet Theorem, which provides a factorization of the state transition matrix of a Linear Periodic System, this factorization allows to determine the stability of the LP system (3) from the algebraic localization of the eigenvalues of the Monodromy Matrix. Also, a result concerning the determinant of some special sum of matrices is given.
Consider (3), its solution is given by
where is the state transition matrix and is given by the solution of   
The periodicity of the system leads to the periodicity of the state transition matrix
The transition matrix over one period is defined as the Monodromy matrix M
Remark 2.1. is independent on   .
One of the basic tools used for the stability analysis of periodic systems is based on Floquet theory.
Theorem 2.1 (Floquet  ). Consider the homogeneous system given by (3), then there exists a periodic invertible matrix and a constant matrix such that the state transition matrix of the system can be written as
Evaluating (8) over one period T leads to
That is, the Monodromy matrix M is similar to. Assume, then, where and, then
In (10), the factor and are bounded, therefore will be bounded if and only if is bounded; i.e., boundedness of the solution of (3) depends only on the eigenvalues of the Monodromy Matrix which from (9) are the same as the eigenvalues of, there the importance of the calculation of M.
Theorem 2.2 (Lyapunov  ). The periodic system (3) is:
1) Exponentially stable if.
2) Stable if and if such that then is a simple root of.
3) Unstable if such that, or if and such that and it is a multiple root of
The following results are useful in the analysis of the rest of the paper. Let represent the identity matrix of dimension.
1In (*) we have used the fact:, and in (**) we have used: if, then.
Theorem 2.3. Let R denote an matrix, S an matrix, T an matrix, and U an matrix. If R and T are non singular then
Lemma 2.1 (  ) For all, and for all:
The previous Lemma states that, and therefore,.
Corollary 2.1. Consider a set of matrices:, where, all matrices of dimensions, then
The same conclusion of Lemma 2.1 applies, i.e.
3. From a Linear Periodic System to a Discrete Shift Invariant System
Consider a system given by (1). There are, in principle, two ways to determine its stability. One using numerical integration of the solution and the other is to approximate (1) by a discrete time system and solve it. Then, the solution of the latter is an approximation of. In this paper the last method is used. The advantage of such approach is that it is always possible to find an analytical solution of the approximated discrete problem.
Lemma 3.1 (  ). Any p-norm of the approximation error, under finite difference methods, of functions defined on a finite interval is bounded.
Remark 3.1. Observe that knowing the solution of (3) in the interval allows to know the solution everywhere.
Performing a discretization based on the centered finite-difference method  to (1), the following approximation is obtained
where, and with, is the sampling period.
・ Centered finite-difference has an error of approximation of of order.
・ In the sampling process, the continuous-time period T has to be divided in an integer number of samples per period, K; i.e., , with K a positive integer. This restriction guarantees that the obtained discrete-time system is K-periodic. If this constrain is not satisfied, the resulting discrete-time system is quasi-periodic and the Floquet Theory is not longer valid  .
As decreases, (16) approximates better to the behavior of the original Equation (1). In state space (16) is rewritten as
In (17) is K-periodic. During one period K, the index k takes the values, so it follows that
but in the next period
and so on. Then, using the Lifting Technique  the next Shift Invariant system is obtained.
Remark 3.3. Roughly speaking, the Lifting Technique augments the dimension of the state K times. Then the augmented system evaluated every K samples turns to be shift-invariant.
In (20), so the dimension of the Shift Invariant System increases as the sampling period decreases. Equation (20) is equivalent to (17), so as decreases, the solution of (20) gives a better approximation of the solution of (1), but it is important to notice that the stability problem of (20) can be solved analytically. The stability analysis of (17) is now simpler: is a block diagonal Matrix and because of Corollary 2.1 each block of has the same eigenvalues. Moreover, the last block of results to be the Discrete Monodromy Matrix associated to (17), then the stability of the Shift Invariant System (20) is reduced to the analysis of any of the diagonal blocks of. The first contribution of this paper is then stated as follows.
Theorem 3.1. Consider an homogeneous periodic differential equation
where and T is the minimum period, and its sampled approximation
where with K an integer. Then
1) (23) can always be solved analytically via Lifting Technique, which leads to a Shift Invariant System
2) The stability of (23)2 is given by the eigenvalues of any of the blocks of the constant matrix .
2Recall that the shift invariant system is
(a) Exponentially stable if.
(b) Stable if and if such that then is a simple root of.
(c) Unstable if such that, or if and it is a multiple root of.
Remark 3.4. The solution of the Shift Invariant System (23) and its stability analysis under small sampling periods, are approximately equivalent to their counterparts in the continuous time system (22).
Example 3.1. As an example let the periodic function in (22) by given by. In Figure 1 the Arnold Tongues computed with the method proposed in this document are shown. The gray zones represent unstable solutions of (22) while white zones represent stable solutions. Only this graph is presented since its difference to the one obtained via numerical integration is indistinguishable.
Notation: The first instability region is labeled as 0, and the following accordingly
with natural numbers, the Arnold Tongues starts at and we refer to the k-th
Arnold Tongue  .
The same stability diagrams were also computed using numerical integration, both algorithms on a computer with a Intel Core 2 Duo processor at 2.6 GHz and 1 Gb of RAM, having in both cases 360 samples per period, and the same increment for the parameters and, the computational time comparison is shown in Table 1.
Table 1 shows that the computational time of the proposed method is much lower (approximately 20 times faster) than the time needed for computing the same chart using numerical integration.
Figure 1. Arnold Tongues computed using the method proposed in this document. Observe in Table 1 that the computation based on the lifting technique is considerably faster than the classical numerical integration method.
Table 1. Time comparison of computation methods of the stability diagram of Figure 1. The parameters are and.
4. Vibrational Stabilization
3, where we denote by n the number of the Arnold Tongue to which the pair belongs.
In this section the second contribution of this paper is presented. Suppose that the pair is fixed, where, such that
has unbounded solutions. Let us introduce a vibrational control scheme which consists on modifying by, where is -periodic3, i.e.,
. This results in the forced system
The goal is to design a controller such that the system
is stable for. The following result shows how to design and.
Theorem 4.1. Consider two linear periodic systems, of the same dimensions, and of the same structure
Denote by and the corresponding monodromy matrices and suppose that both systems are unstable, i.e., and, under the ad- ditional condition that. Then there exists a constant such that for the system:
Note that, and therefore (32) it is stable.
Remark 4.1. If belongs to the n-th Arnold Tongue for (28) and belongs to the Arnold Tongue for (30), then the condition holds.
Proof. The proof is performed in discrete time. The Monodromy Matrix associated to (32) is:
which is a polynomial and therefore a continuous function of. It follows that
Therefore, from the hypothesis and from the con- tinuity of, we have that there exists a constant for which.
・ Theorem 4.1 does not only guarantees the existence of a stabilizing constant, but provided with Theorem 3.1, it is possible to explicitly compute such constant.
・ In the proof it is possible to choose and the same stability result holds. By choosing a stable point more or less at the middle of the corresponding stable interval is obtained.
Figure 2 shows schematically how the vibrational stabilization method performs.
Theorem 4.1 might be also written for Hill systems as follows.
Corollary 4.1. Assume that the pair, where, is given such that the solution of the system
is unstable. Then there exists a constant and a -periodic function
such that for the same pair, the system
has a stable solution.
Proof. Recall that denotes the number of the Arnold Tongue to which the pair belongs for the stability chart of (35). From what has been mentioned in this paper, if for it happens that, the effect is that the right contiguous Tongue is going to be modified  .
Remark 4.3. Let the Fourier series of be written as. Then the boundaries of the k-th Arnold Tongue are tangent at and if and
only if, and are transversal if and only if  . Such results means that we can design a controller that modifies a specific Tongue.
It only rests to find. Let
Figure 2. A 3 dimensional view of the vibrational stabilization method. The control scheme introduced changes the shape of the stability chart in such a way that for some, the pair belongs to a stable zone.
therefore (36) is rewritten as
It is clear that (38) is in the same format given in Theorem 4.1. Then, it is guaranteed the existence of such that the solution of (36) is stable.
Consider the system (Kapitza Pendulum   ) shown in Figure 3 (this is the prototype example used to witness the effects of parametric resonance, see e.g.  or
 ), where: and, then, the Hill equation
which models the system is
Or after linearization around the lower equilibrium point
Note: Even though this paper developed the stabilization method for linear systems, we will show that it also works for the nonlinear system.
The operation point, i.e., the point produces an unstable response as can be seen in Figure 4.
The results of this paper show the existence of a constant gain and a periodic function such that the new system given by
Figure 3. Kapitza pendulum.
Figure 4. Arnold Tongues for the example. The operating point is located at .
Figure 5. Solution of the non-linear model.
has a stable solution. In this example, the resulting equation turns to be
The corresponding stability diagram of (42) is shown in Figure 7.
The response of (42) is depicted in Figure 8. The same stabilization scheme might be applied to the non-linear equation, this means
The corresponding response is shown in the Figure 9. Note that for small am-
Figure 6. Solution of the linear model.
Figure 7. Modified Arnold Tongues by the vibrational controller.
Figure 8. Stable solution of the linear model under the action of the vibrational control.
plitudes, the linear model resembles very good to the non-linear one. Also, note that in general, the stable solutions of a Linear Periodic System are not periodic, but almost- periodic  .
Figure 9. Stable solution of the non-linear model under the action of the vibrational control.
This paper presents an alternate and new method to compute the Arnold Tongues of a Hill equation, which is much faster than the traditional numerical integration method. Since the proposed algorithm calculates an approximation of the Monodromy matrix, it is possible to know, with some small errors due to the approximation, whether a given pair is a stable or an unstable operation point.
Moreover, a vibrational stabilization scheme for a given (unstable) Linear Periodic System is given. The approach presented is based on the capability of modifying the shape of the Arnold Tongues. Although the method proposed here was based on re- shaping the contiguous Tongue from which the pair belongs to, com- putationally speaking we may modify any Tongue as long as it is different from the Tongue to which the pair belongs originally. The proof follows the same lines developed in this paper.
H.J.K. thanks CONACyT and CINVESTAV-IPN for their support during his M.Sc. studies, while this document was written.
4Here, although the function q(mt) is not of period T (since m Î N) p(t) + q(mt) is of period T.
 Gel’fand, I.M. and Lidskii, V.B. (1958) On the Structure of the Regions of Stability of Linear Canonical Systems of Differential Equations with Periodic Coefficients. Uspekhi Matematicheskikh Nauk, 10, 3-40 (English Translation: AMS Translations, Series 2, 2, 143-181).
 Starzhinskii, V.M. (1955) A Survey of Works on the Conditions of Stability of the Trivial Solution of a System of Linear Differential Equations with Periodic Coefficients. American Mathematical Society Translations, 1, 189-237.
 Starzhinskii, V.M. (1963) On the Stability of Periodic Motions I and II. Buletinul Institutului Politehnic din Iasi, Serie Noua, 5, 51-100 (Buletinul Institutului Politehnic din Iasi, Serie Noua, 4, 9-68, 1958).
 Mailybaev, A.A. and Seyranian, A.P. (2009) Stabilization of Statically Unstable Systems by Parametric Excitation. Journal of Sound and Vibration, 323, 1016-1031.
 Jia, Y., Yan, J., Soga, K. and Seshia, A.A. (2014) Parametric Resonance for Vibration Energy Harvesting with Design Techniques to Passively Reduce the Initiation Threshold Amplitude. Smart Materials and Structures, 23, Article ID: 065011.
 Scapolan, M., Tehrani, M.G. and Bonisoli, E. (2016) Energy Harvesting Using Parametric Resonant System Due to Time-Varying Damping. Mechanical Systems and Signal Processing, 79, 149-165.
 Seyranian, A.P., Yabuno, H. and Tsumoto, K. (2005) Instability and Periodic Motion of a Physical Pendulum with Vibrating Suspension Point (Theoretical and Experimental Approach). Doklady Physics, 50, 467-472.