#### Abstract

A systematic research on the structure-preserving controller is investigated in this paper, including its applications to the second-order, first-order, time-periodic, or degenerated astrodynamics, respectively. The general form of the controller is deduced for the typical Hamiltonian system in full feedback and position-only feedback modes, which is successful in changing the hyperbolic equilibrium to an elliptic one. With the poles assigned at any different positions on imaginary axis, the controlled Hamiltonian system is Lyapunov stable. The Floquet multiplier is employed to measure the stability of time-dependent Hamiltonian system, because the equilibrium of periodic system may be unstable even though the equilibrium is always elliptic. One type of periodic orbits is achieved by the resonant conditions of control gains, and another type is making judicious choice in the foundational motions with different frequencies. The control gains are selected from the viewpoint of both the local and global optimizations on fuel cost. This controller is applied to some astrodynamics to achieve some interesting conclusions, including stable lissajous orbits in solar sail’s three-body problem and degenerated two-body problem, quasiperiodic formation flying on a -perturbed mean circular orbit, and controlled frozen orbits for a spacecraft with a high area-to-mass ratio.

#### 1. Introduction

It is well known that most of the astrodynamical problems could be classified as hyperbolic Hamiltonian systems, for example, circular restricted 3-body problem (CR3BP). There are the hyperbolic equilibria which have stable, unstable, and centre manifolds for hyperbolic Hamiltonian system, just like the collinear libration points , and in CR3BP [1], of which the unstable manifolds may cause instability in the Lyapunov sense.

Some astronautical missions, requiring increasing the coverage rate to the object (ground station or target spacecraft) or avoiding communication signals being lost in the sun, are preferring bounded motions near hyperbolic equilibria, such as lissajous (or halo) orbit generated in CR3BP [1, 2]. There presented two different approaches to generate bounded relative trajectories, of which one is obtaining the bounded orbits by making judicious choice in initial conditions. Certain zero-steady controllers could be applied to stabilize the nominal trajectories, owing to the hyperbolicity inheriting from the topological type of the equilibrium.

However, for the other case without a natural trajectory to track, some station keeping maneuvers are necessary to implement modification of the topology of the equilibrium. This paper proposes a structure-preserving controller to generate bounded trajectories around the equilibrium, which has many potential applications in astrodynamics. Firstly, a simple structure-preserving controller was constructed by Scheeres et al. [3], in which involved only the instantaneous stable and unstable manifolds, and the manifold’s gains are constrained to be equal. By choosing the proper manifolds’ gains, M. Xu and S. Xu [4] have demonstrated that the poles of the typical second-order Hamiltonian system can be allocated to any position along the imaginary axis and then obtained a new quasiperiodic orbit type referred to as a stable lissajous orbit. Therefore, the controller was powerful enough to work out some famous and difficult problems in sail’s restricted three-body problem that there exist the bounded trajectories near the off-axis solar sail equilibrium [4], and in maintaining the eccentricity in the neighborhood of the unstable fixed point. Besides, the controller implemented the stabilization on a saddle-node bifurcation point for the degenerated case, which makes the modified elliptic equilibrium unique for the controlled system to change the instability of the equilibrium [5]. Further work was implemented on the stabilization of the equilibrium for time-periodic system, which has time varying topological types and no fixed-dimensional unstable/stable/center manifolds [6]. The 2-dimensions temporally independent systems possess fixed-dimensions of unstable/stable/center manifolds and a permanent pair of hyperbolic eigenvalues. Tabuada and Pappas [7] focused on local accessibility preserving abstractions, and provided conditions under which local accessibility properties of the abstracted Hamiltonian system are equivalent to the local accessibility properties of the original Hamiltonian control system.

A systematic investigation on the structure-preserving controller is presented in this paper, including some applications in second-order, first-order, time-periodic, and degenerated Hamiltonian astrodynamics, respectively. The general form of the controller is designed for the typical second-order Hamiltonian system, which can work as full feedback and position-only modes. The controller is successful in changing the hyperbolic equilibrium (saddle) to an elliptic one (center) with the poles on the imaginary axis (marginal stability), and then the controlled Hamiltonian system can achieve the Lyapunov stability by means of the Morse lemma. In contrast to the time-independent Hamiltonian system, the equilibrium of periodic system may be unstable even though the equilibrium is always elliptic during its period. Therefore, it is necessary to check whether the Floquet multipliers lie on the unit circle in the complex plane, or whether the moduli of all the multipliers are equal to 1. Two different approaches to generate the periodic orbits are discussed, one of which is achieved by the resonant conditions of control gains for the resonant periodic orbits, and another is making judicious choice in the foundational motions with different frequencies. Furthermore, some selection techniques on the gains are considered from the point view of both the local and global optimizations.

This controller is applied to some astrodynamics to achieve some interesting conclusions, including stable lissajous orbits in solar sail’s three-body problem and degenerated two-body problem, quasiperiodic formation flying on a -perturbed mean circular orbit, controlled frozen orbits for a spacecraft with high area-to-mass ratio.

#### 2. Hamiltonian System and Hyperbolic Equilibrium

##### 2.1. Typical Hamiltonian System in Astrodynamics

For the typical Hamiltonian system in astrodynamics, the Hamiltonian function has the following form as [4]: where is the general form of pseudopotential function relying only on the coordinate and time; is the characteristic parameter introduced by the Coriolis acceleration, such as the mean angular motion of the synodic frame used in CR3BP, and Local Vertical-Local Horizontal (LVLH) frame used in formation flying.

In general, the variables can be expressed by the physical variables from the Legendre transformation: where is the identical matrix, and is the symplectic matrix. The astrodynamics in physical variables can be expressed as For a hyperbolic Hamiltonian system, there exist hyperbolic equilibria that have stable, unstable, and center manifolds [1], which are similar to the collinear libration points , , and in CR3BP. Hence, the equilibrium of the typical 2-dimensional Hamiltonian system in astrodynamics can be obtained from the following equations: Denote as the second derivative matrix of the pseudopotential function about the position vector , and the elements of the matrix have the general expression as . And then a hyperbolic (saddle) equilibrium results from .

Generally, unstable manifolds of the hyperbolic equilibrium may cause instability. The spacecraft needs to be located in the unstable region around the hyperbolic equilibrium in order to meet the requirements of some missions. Thus, it is necessary to establish a stabilizing strategy for the hyperbolic equilibrium.

##### 2.2. Typical Hyperbolic Equilibrium in Astrodynamics

For the typical second-order Hamiltonian system, different and can be used to define different astrodynamical problems, illustrated in the following section.

and define the circular restricted three-body problem (CR3BP) [8], where is the mass ratio between the primaries, and , are, respectively, the distances between the spacecraft and the two primaries, that is, and . There are five libration points, three collinear ones of which are hyperbolic, labeled as , and shown in Figure 1.

and define the planar solar sail three body problem [4], where is the solar pressure acceleration and is the sail surface normal vector. There are infinite equilibria existing in this system, which can be parameterized by the sail lightness number : where is the sail lightness number which is the ratio of the radiation pressure acceleration to gravitational acceleration. This parameterization for Sun-Earth/Moon system generates level surfaces as demonstrated in Figure 2.

and define the solar sail’s two-body problem [1], where is the constant determined by solar pressure; the polar coordinates are parallel to the sun line and the minimal radius from the axis, respectively; is the constant angular momentum directed along the sun line; . indicates that this Hamiltonian system is degenerated for the absence of the Coriolis acceleration and a unique double equilibrium due to .

Considering the formation flying on a -perturbed mean circular orbit, the linearized relative dynamics can be deduced as [6, 9] where is the relative position in the chief’s LVLH rotating frame, is the gravitational potential function of the deputy including perturbation, is the coordinate transformation matrix from the LVLH to inertial frames, and is the gradient vector described in orthogonal or spherical coordinates. Due to the osculating orbital elements, the relative dynamics on a -perturbed mean circular orbit is time-periodic Hamiltonian system. This time-periodic system has the fixed equilibrium , that is, the location of chief; however, the equilibrium of the system has time-varying topological types and no fixed dimensional unstable/stable/center manifolds, which are quite different from the two-dimensional time-independent system with a permanent pair of hyperbolic eigenvalues and fixed dimensions of unstable/stable/center manifolds, as shown in Figure 3.

**(a)**

**(b)**

Some averaging techniques can be used to reduce the typical astrodynamics into a first-order Hamiltonian system. For a spacecraft with high area-to-mass ratio orbiting the Earth, its dynamics is strongly perturbed by the term of the gravitational field due to the Earth’s oblateness and by the effect of solar radiation pressure (SRP). The secular rate of the orbital elements due to SRP and term can be written as [10] where and are the constant parameters measuring the solar radiation pressure and the oblateness. The reduced system exists three equilibria, that is, and , where is hyperbolic and the others are elliptic, as shown in Figure 4.

#### 3. Structure-Preserving Stabilization for Hamiltonian System

##### 3.1. Controller Design Preserving Hamiltonian Structure

For the typical second-order Hamiltonian system, the variation equation around the equilibrium is: For this unstable system, there exist hyperbolic eigenvalues associated with the stable and unstable manifolds , and elliptic eigenvalues associated with the centre manifolds and .

The Hamiltonian structure-preserving controller then is constructed as where is used to change the Coriolis acceleration, and are, respectively, the gains of unstable, stable, and centre manifolds. The symmetry of matrix and the skew symmetry of matrix guarantee the linear feedback controller preserves the Hamiltonian structure [3].

##### 3.2. Stability Analysis

For the 2-dimensional and time-independent Hamiltonian system, some theorem and proposition about the stability of the controller can be deduced as follows.

Theorem 1. *The poles can be assigned at any different positions in imaginary axis, and , and required are not unique.*

*Proof. *Define the matrix as
where and are composing the stable and unstable manifolds of the hyperbolic eigenvalues as
and and are composing the center manifolds of the elliptic eigenvalues as
Thus , and will be transformed by the controller:
where , and can be assigned arbitrarily to satisfy the linearly stable conditions, so long as matrix is not regular. Obviously, cannot be regular due to the fact that is different from for the different eigenvalues and .

Denote and as the solutions to the characteristic equation of the second-order dynamics, and then and can be expressed as
Consequently, , and can be solved from and .

Furthermore, and are restricted by
but if one sets
then . Since this system has second-order elementary factors, the Jordan form of will have long-term dispersions because the matrix cannot be diagonalized.

Because of the nonuniqueness of , and solved from (9), , and are nonunique. Suppose that sets of and generate the same values of and in (9), and the two sets have the following relationships as
It is obtained from (7) that
then the real solution is:
we can fix to get two different sets generated from and . Besides, the fact that there exist different sets generating the same values of and indicate , and are nonunique.

Hence, the controlled Hamiltonian can be expressed as Expand around the equilibrium region to get where represents the Hamiltonian value at the equilibrium and the first-order polynomial disappears because of the equilibrium.

Proposition 2. *The Coriolis acceleration can be modified by , which cannot stabilize the system independently for does not change , and .*

The controller developed here can transform the hyperbolic equilibrium (saddle) to an elliptic one (center) according to the theorems and propositions. The elliptic equilibrium has a linear symplectic transformation, which transforms to the following form: According to Morse lemma, there exists an analytic diffeomorphism around the equilibrium, which transforms to the following form: thus the energy surfaces are locally diffeomorphic to a family of spheres which shrink down to the equilibria: . The Lyapunov stability follows because the trajectories are tangent to the energy surfaces [11]. Therefore, we have obtained that the nonlinear full dynamics is also Lyapunov stable.

Theorem 3. *The gains can be chosen large enough to guarantee the modified elliptic equilibrium is the only equilibrium for the controlled system.*

*Proof. *Since it is easy to demonstrate that is negative definite, the geometric structure of the Hill’s region for the controlled system will be changed by adding to the Hill’s region which depends on .

Another equilibrium of the controlled system (if exist) will stay in Island or Mainland [1] because it must stay in the region where and for the negative definition of .

Here it is just illustrated that the equilibrium cannot stay in Mainland, and the negative inclusion for Island can be obtained by the similar way. We can verify whether is true by increasing and , since is the bottom point of Mainland; moreover, we need to test which one descends faster between ( and . The comparison is very tedious but simply, thus it will be ignored here.

For the controlled system, there exists only the modified elliptic equilibrium, which means the equilibrium has the global minimum pseudopotential. So the controller’s stabilizing region is the whole phase space.

Furthermore, the equilibrium of time-periodic system may be unstable even if the equilibrium is always elliptic during its period, compared with the independent Hamiltonian system. Consequently, it is necessary to verify if the modules of all the multipliers are 1, or if the Floquet multipliers lie on the unit circle in the complex plane. Considering the relative dynamics on a -perturbed mean circular orbit, the Floquet stability of the periodic Hamiltonian system is obtained by the larger controller gain rather than the critical gain , as shown in Figure 5.

**(a)**

**(b)**

##### 3.3. Resonant Condition

A resonance happens when the eigenvalues of a controlled system have relation , where and and are reduced with each other. Thus the trajectories will be periodic for any initial condition of spacecraft. We will give a particular solution for the resonance in the following section.

Suppose , and then: hence, A particular solution can be constructed as If , a particular solution can be constructed as where is an arbitrary parameter.

Compared with the osculating frequencies of the controlled elliptic equilibrium for the global investigation during its period, the characteristic frequencies defined by the Floquet multipliers have more applications for the time-periodic Hamiltonian system. Floquet multipliers are the eigenvalues of the linear Poincaré mapping with the unit eigenvectors , and . Characteristic frequencies and are approximately equal to zeros inheriting from two of osculating frequencies around , while depends on the controller gain. The linear combinations of eigenvectors, , and , are utilized to generate six different types of near-periodic relative orbits with their periods of and , respectively. The six different types of bounded relative orbits generated by the linear combinations can be then regarded as the foundational motions around the equilibrium.

Particularly, the linear combinations and can be used as the initial values to generate the periodic relative orbits with their periods of when satisfies the resonant condition . To be stricter, all the relative periodic orbits, which produced by the judicious choice combined of the eigenvectors, are near periodic only for the nonlinear term that ignored by the linearized differential equations describing the relative dynamics.

It is worth mentioning that the six eigenvectors of the Floquet multipliers span the entire space of the relative position and velocity . Therefore, any initial relative positions and velocity vectors can be decomposed as where is constant coefficients, and is chosen among , and . Hence, all the initial relative position and velocity can generate a bounded trajectory around the controlled elliptic equilibrium, which depends on the topological type of the equilibrium. Moreover, all the general trajectories are quasiperiodic for the different frequencies between the foundational motions and all the bounded trajectories are involved in the foundational motions caused by the eigenvectors.

##### 3.4. Cost and Optimization

Because , and are nonunique, it is necessary to investigate the rule of choosing , and to stabilize the system. Denote as the collection of , and which will allocate the controlled system with the expected poles and .

The controller’s output can be then determined as Define to specify the sensitivity of the controller as .

The Frobenius norm which is consistent to the Euclidean norm of the vector, can be used to measure the controller’s sensitivity as

Define and then deduce the matrix from and as

which has the following relationship as: Thus, put it into (33) and (30) to obtain: where can be calculated directly from the constant matrix according to the invariant manifolds of the equilibrium.

For full feedback, can be treated as the optimization index to select more suitable values for the manifolds gains; however, for the position feedback, the optimization index can be set as . With the constraint equation (15), the selection for controller’s gains can be considered as a nonlinear constrained optimization, as: The programming functions “fmincon” and “confuneq” in MATLAB can solve this optimization problem.

As confirmed by zero real part of the eigenvalues in Figure 6, the equilibrium is elliptic during the period for time-periodic Hamiltonian system. Therefore, there exist three osculating imaginary eigenvalues , and for the controlled equilibrium. The quasiperiodic relative trajectories around the equilibrium comprise three different oscillating motions with the osculating frequencies of , and , which can be verified by the second-order term in the new Hamiltonian function yielded by a linear symplectic transformation:

**(a)**

**(b)**

According to the local optimization, less cost is obtained by small controller gain illustrated by (34); however, the smaller gain cannot maintain the relative position in a domain closer to the equilibrium, which may consume more fuels.

According to the global optimization, an averaging quadratic cost function is defined to measure the fuel consumption during for the continuous thrust engine [12]: where is the position component of the periodic orbit which is developed by the eigenvector . The component in the th row and th column of the measuring matrix is , which inherits the near-periodicity from the characteristic frequencies and depends only on the controlled elliptic equilibrium. Therefore, inherited from the longest period of the characteristic frequency , that is, , can serve as the foundational measuring matrix that can be obtained offline (or on the ground) from the linearized dynamics for the fast and approximative estimation. The consumption function during has the form as: where is symmetric and positive-definite and . Generally, there exists a unitary matrix from Singular Value Decomposition (SVD) to transform the function as , where represents the eigenvalues of and is the th column element of the intermediate variable .

The identified gain can be derived to minimize if we know the initial relative position and velocity. The specific gain satisfying the minimal trace of , that is, , is however preferable to minimize for the maximum likelihood, for an unknown case on the initial conditions.

#### 4. Application to Nonlinear Astrodynamics

##### 4.1. Stable Lissajous Orbits of Solar Sail: Application to Second-Order Hamiltonian System

A solar sail is a new type of spacecraft with no fuel, and its orbital control is realized by the solar radiation pressure by means of orientating its attitude relative to the Sun. Therefore, we can neglect the control consumption and implement the complicated control law, which is incomparable to the low-trust ion propulsion or the impulse propulsion.

We noticed the fact that there exists quasiperiodic or periodic orbits near the sail’s equilibria, and several periodic orbits have been proposed in the solar sail restricted three-body problem. A survey [13] aimed to summarize the combined wealth of literature concerned with the dynamics, stability and control of highly non-Keplerian orbits for various low thrust propulsion devices, and to demonstrate some of these potential applications, like Capturing near earth objects [14]. McInnes [15] is the first person who applied the classical libration point orbit theory to generate the halo orbit around the on-axis sail’s equilibrium. Baoyin and McInnes [16] achieved two different halo orbits types around the on-axis sail’s equilibrium. McInnes [15] has tried to generate lissajous orbits near the equilibrium as well, but failed to keep a bounded trajectory for a long time. Waters and McInnes [17] applied the classical theory to generate some halo orbits that is only available for the off-axis sail equilibria in some specified regions. What is more, all the interesting orbits presented by researchers are unstable and have one-dim unstable manifolds, hence the station-keeping strategies becomes quite necessary.

The sail associated with the Sun and Earth-Moon system, (where the Earth and Moon are regarded as a whole celestial body, denoted as Earth/Moon), is regarded as the solar sail restricted three-body system. In this section we suppose that the Sun-Earth/Moon system revolves around a circular orbit.

The structure-preserving controller is proposed to generate stable lissajous orbits and Lyapunov orbits for nonresonance shown in Figure 7, and resonant periodic orbits shown in Figure 8 (1 : 2, 1 : 3, 1 : 4, 2 : 3, resp.).

The feedback gains are initially chosen with and in view of the local optimization, and then the gains will be optimized so as to obtain the same poles. The gains refined by optimization are , and , with the optimization index improved from 600 to 483 (80.5%) for the position feedback. While for full feedback, the gains refined by optimization are , and , with the optimization index improved from to (77.71%).

##### 4.2. Bounded Orbits by Degenerated Solar Sail: Application to Degenerated Hamiltonian System

McInnes [18] has researched the displaced non-Keplerian orbits above the Sun. Similar steps have been taken by Bookless and McInnes [19, 20] to the displaced non-Keplerian orbits above the planet. Nevertheless, just as the stability analysis for motions around the equilibrium, Bookless’ analysis for the dynamics and control is linear and local. There presented the nonlinear dynamical analysis for a displaced orbit above a planet, and there also investigated the motions around the equilibria for the nonresonance case with the help of the Birkhoff normal form and dynamical system techniques. M. Xu and S. Xu [1] has obtained two important contributions to derive necessary and sufficiency conditions for the motion stability around the equilibria, and to illustrate numerically that the non-transition critical KAM torus is filled with the Lyapunov (1,1)-homoclinic orbits. defines the degradation case for the two body problem of solar sails, which is the saddle-node bifurcation point.

For the degenerated system, there exists only one double equilibrium, and the frequencies for the liberalized motion are, respectively, and . What is more, it is demonstrated from the variation equation around the equilibrium that the equilibrium has no stable or unstable manifolds but centre manifolds. Owing to Theorem 1, it guarantees that the stabilization for the degenerated can be achieved by the equilibrium centre manifolds.

Similar to the nondegradation case, the suitable gains can be chosen to guarantee the modified elliptic equilibrium is the only equilibrium for the controlled system. So the controller’s stabilizing region is all the phase space. The quasiperiodic orbits generated by the controller are shown in Figure 9, and their positions in 3D physical space are shown in Figure 10.

**(a)**

**(b)**

**(c)**

**(d)**

##### 4.3. Quasiperiodic Formation Flying: Application to Time-Periodic Hamiltonian System

Roughly, we can classify the relative motion control approaches into two branches: continuous control that using low-thrust electric propulsion and impulsive control that relying on chemical thrusters. While the former approaches used mostly full-state feedback of the Cartesian relative position and velocity to develop high-accuracy tracking control laws, thus capable to maintain small steady-state errors under a myriad of orbital perturbations [21–25]; however, the maintenance of the later assumes that the tracking reference orbital elements are known and attempts to generate control commands that will match the osculating orbital elements to the expected values [12, 26, 27].

The initial relative position of and the relative velocity generated from the classical equations for the formation flying on a -perturbed mean circular orbit, are applied to generate the quasiperiodic relative trajectories with the preferred gain 0.31 illustrated in Figure 11. More rapid changes in configuration are obtained by the controller developed in this paper, rather than by the change in invariant relative orbits caused by the secular term of perigee argument, from the projection of the configuration on the or plane. Therefore, the mission that requires rapid changes in formation configuration, like on-orbit inspection and repair could potentially use the addressed controller. Furthermore, the quasiperiodic relative trajectories which comprise six foundational motions with different frequencies have no reference relative configuration to track for the continuous thrust, make the trajectory prediction difficult by other nonallied spacecrafts.

**(a)**

**(b)**

**(c)**

**(d)**

Figure 12 illustrates the topological type of the equilibrium of the controlled system, which is elliptic during the period, as confirmed by zero real part of the eigenvalues, for the mean circular orbit whose mean inclination is 80°, mean semi-major axis is 7178.137 km, and initial argument of latitude is 0°.

**(a)**

**(b)**

should share its value with to minimize the controller outputs from view of local optimization, according to the contour curves of illustrated in Figure 13. The specific gain owning the minimal trace of , that is, , is preferable to minimize for the maximum likelihood if the initial conditions are unknown. The optimizing result obtained from local optimization is illustrated in Figure 14 that the optimal gain of 0.31 has a smaller trace than the critical gain of 0.27.

##### 4.4. Controlled Frozen Orbits: Application to First-Order Hamiltonian System

The high area-to-mass spacecraft’s orbital evolution behaves strangely under the influence of the solar radiation pressure and the perturbation due to the oblateness of the Earth. Since 1960, there have been observed long-term oscillations of the eccentricity in the satellites’ orbital behaviours, like Vanguard and the ECHO balloon. The two variables, the osculating orbital eccentricity and the solar angle between the orbit pericentre and the Sun-line [10, 28], can be used to describe this dynamical system with low inclinations through a Hamiltonian function.

All spacecrafts around the hyperbolic eccentricity cannot be maintained on such an orbit; on the contrary, they will move on the unstable manifold towards higher eccentricities, or on the unstable manifold towards lower eccentricities based on the quite small change of the orbital elements. Controlled frozen trajectories can be obtained in the phase space around the hyperbolic eccentricity by means of applying the Hamiltonian structure-preserving controller, such as the one illustrated by red dotted line in Figure 15.

**(a)**

**(b)**

The mission illustrated in Figure 15(a) comprises three phases: (i) transferring (blue line) from a circular orbit to the hyperbolic equilibrium, (ii) stabilizing near the hyperbolic equilibrium (red line), and (iii) transferring (dotted blue line) beyond the critical eccentricity. While the mission illustrated in Figure 15(b) comprises three phases as well: (i) transferring from a circular orbit to the hyperbolic equilibrium (blue line), (ii) stabilizing (red line) near the hyperbolic equilibrium, and (iii) transferring (the dotted blue line) back to the circular orbit.

#### 5. Conclusion

Bounded motions near hyperbolic equilibria have been broadly applied to various astronautical missions. The Hamiltonian structure-preserving controller can change the hyperbolic equilibrium (saddle) to an elliptic one (center), and then some quasiperiodic bounded trajectories will emerge naturally from the KAM tori generated by the controller.

In this paper, a systematic investigation is presented on the general form for full feedback and position-only feedback modes, Lyapunov and Floquet stability analysis, and control gain selection of the structure-preserving controller, in the context of second-order, first-order, time-periodic, and degenerated Hamiltonian system, respectively. One type of periodic trajectories is achieved by the resonant conditions of control gains for the resonant periodic orbits, and another type is making judicious choice in the foundational motions with different frequencies. Therefore, the controller has potential applications in stable lissajous orbits for solar sail’s three-body problem and degenerated two-body problem, quasiperiodic formation flying on a -perturbed mean circular orbit, or the controlled frozen orbits for a spacecraft with high area-to-mass ratio.

Furthermore, some further researches can be implemented, including how to use unfixed dimensional center manifolds to feedback the controller, the further investigation for the global Floquet stable controller gain for quasiperiodic Hamiltonian system, and the relationships between the different frequencies and the critical gain for Floquet stability or the optimal gain.

#### Acknowledgments

The research is supported by the National Natural Science Foundation of China (11172020), the National High Technology Research and Development Program of China (863 Program: 2012AA120601), Talent Foundation supported by the Fundamental Research Funds for the Central Universities, Aerospace Science and Technology Innovation Foundation of China Aerospace Science Corporation, and Innovation Fund of China Academy of Space Technology.