JOURNAL OF GUIDANCE,CONTROL, AND DYNAMICS Vol. 37, No. 2, March–April 2014
Virtual Formation Method for Precise Autonomous Absolute Orbit Control
S. De Florio∗ University of Glasgow, Glasgow, Scotland G12 8QQ, United Kingdom S. D’Amico† DLR, German Aerospace Center, Wessling 82234, Germany and G. Radice‡ University of Glasgow, Glasgow, Scotland G12 8QQ, United Kingdom DOI: 10.2514/1.61575 This paper analyzes the problem of precise autonomous orbit control of a spacecraft in a low Earth orbit. Autonomous onboard orbit control means that the spacecraft maintains its ground track close to a reference trajectory, without operator intervention. The problem is formulated as a specific case of two spacecraft in formation in which one, the reference, is virtual and affected only by the Earth’s gravitational field. A new parametrization, the relative Earth-fixed elements, is introduced to describe the relative motion of the two spacecraft’s subsatellite points on the Earth surface. These relative elements enable the translation of absolute-into-relative orbit-control requirements and the straightforward use of modern control-theory techniques for orbit control. As a demonstration, linear and quadratic optimum regulators are designed and compared, by means of numerical simulations, with an analytical autonomous orbit-control algorithm that has been validated in flight. The differences of the numerical and analytical control methods are identified and discussed.
Nomenclature δe = relative eccentricity vector amplitude δ A = spacecraft reference area, m2 ej = relative eccentricity vector component A −1 δF = relative Earth-fixed elements vector = dynamic model matrix, s δ a = semimajor axis, m h = normalized altitude difference B δi = relative inclination vector amplitude, rad = control matrix δ C = output matrix ij = relative inclination vector component, rad δκ = relative orbital elements vector CD = drag coefficient δ e = eccentricity Lj = phase difference vector component, rad δr = normalized relative position vector ej = eccentricity vector component δ G rj = component of normalized relative position vector = gain matrix δ G = disturbance gain matrix u = relative argument of latitude, rad d δv = normalized relative velocity vector gj = gain δ i = inclination, rad vj = component of normalized relative velocity vector ϵ = normalized relative orbital elements vector, m J2 = geopotential second-order zonal coefficient η m = spacecraft mass, kg = normalized altitude ∕ κ = mean orbital elements vector n = reference spacecraft mean motion, rad s λ R = Earth’s equatorial radius, m = longitude, rad E μ = Earth’s gravitational coefficient, m3∕s2 sj = characteristic root
Downloaded by STANFORD UNIVERSITY on August 2, 2015 | http://arc.aiaa.org DOI: 10.2514/1.61575 ξ T = coordinates transformation matrix = thrusters execution error, % ρ ∕ 3 t = time, s = atmospheric density, kg m ϕ = relative perigee, rad u = argument of latitude, rad φ u = constant value of argument of latitude, rad = latitude, rad ∕ Ω = right ascension of the ascending node, rad v = spacecraft velocity, m s Ω_ ∕ Δh = altitude difference, m = secular rotation of the line of nodes, rad s Δ ω = argument of periapsis, rad vj = impulsive maneuver velocity increment component, ω ’ ∕ ∕ E = Earth s rotation rate, rad s m s θ δa = relative semimajor axis = relative ascending node, rad R− = set of negative real numbers R = set of positive real numbers Z = set of integer numbers Received 27 December 2012; revision received 3 June 2013; accepted for publication 3 June 2013; published online 5 February 2014. Copyright © 2013 by the American Institute of Aeronautics and Astronautics, Inc. All Subscripts rights reserved. Copies of this paper may be made for personal or internal use, on condition that the copier pay the $10.00 per-copy fee to the Copyright c = closed-loop Clearance Center, Inc., 222 Rosewood Drive, Danvers, MA 01923; include d = disturbance the code 1533-3884/14 and $10.00 in correspondence with the CCC. *Ph.D. Candidate, Space Advanced Research Team, School of E = Earth Engineering; [email protected]. g = gravity †Research Lead, GSOC/Space Flight Technology, Münchner Straβe 20. N = cross-track ‡Reader (Aerospace Sciences), Space Advanced Research Team, School of R = radial Engineering. R = reference orbit
425 426 DE FLORIO, D’AMICO, AND RADICE
r = reduced II. Virtual Formation Model: Parametrization T = along-track A. Relative Orbital Elements The absolute orbit-control system design can be formulated as the Superscripts control problem of two spacecraft in formation in which one is virtual T = vector transpose and not affected by nongravitational orbit perturbations. The most appropriate parametrization to represent the relative motion of the real spacecraft with respect to the reference is the set of relative orbital elements (ROEs) [26], shown in Eq. (1), where the subscript R refers to the reference orbit: 0 1 0 1 δ − ∕ I. Introduction a a aR aR B δ C B − C B ex C B ex exR C N AUTONOMOUS orbit-control system can fulfil very strict B δ C B − C δκ B ey C B ey eyR C A control requirements on different orbit parameters in real time B δ C B − C (1) B ix C B i iR C and with a significant reduction in flight dynamics ground operations @ δ A @ Ω − Ω A [1]. In the last two decades, several studies on the autonomous orbit iy R sin i δ − control of satellites in a low Earth orbit (LEO) have been conducted u u uR – [2 8]. The TacSat-2 [9], Demeter [10], and autonomous orbit- These parameters are obtained as a nonlinear combination of the keeping (AOK) PRISMA [11] in-flight experiments represent mean orbital elements κ a; ex;ey;i;Ω;u [27–29]. The relative milestones in demonstrating that this technology has now reached a orbit representation of Eq. (1) is based on the relative eccentricity and sufficient level of maturity to be applied routinely to LEO missions. inclination vectors [30] defined in Cartesian and polar notations as The next step is to define a general and rigorous formalization of the autonomous orbit-control problem and explore new control methods. δex cos ϕ δix cos θ This paper analyzes the problem of autonomous absolute orbit δe δe δi δi (2) δe sin ϕ δi sin θ control as a specific case of two spacecraft in formation in which one, y y ’ the reference, is virtual and is affected only by the Earth s The phases of the relative e∕i vectors are termed relative perigee ϕ gravitational field. The control actions are implemented by means of and relative ascending node θ because they characterize the relative in-plane and out-of-plane velocity increments. A new param- geometry and determine the angular locations of the perigee and etrization, the relative Earth-fixed elements (REFEs), analogous to ascending node of the relative orbit. The normalized position δr the relative orbital elements (ROEs) [12–14], is introduced to T T δrRδrT δrN ∕aR and velocity δv δvRδvT δvN ∕ naR vectors describe the relative motion of the real and reference subsatellite of the spacecraft relative to the reference orbit in the RTN orbital points on the Earth surface. The REFEs allow the general formali- frame (R pointing along the orbit’s radius, N pointing along the zation of the absolute orbit-control requirements, which are usually angular momentum vector, and T N × R pointing in the direction expressed through specific Earth-fixed operational parameters such of motion for a circular orbit) can be described through the relative as altitude deviation, phase difference, etc. [15]. A direct mapping orbital elements [26,30] as between REFE and ROE enables the direct translation of absolute δr T into relative orbit-control requirements. By means of this new p δκ formalization, the deviation between the actual and the reference δv T 0 v 10 1 nominal orbit [16] can be defined in an Earth-fixed coordinate system δ analogous to the orbital frame [17]. This approach allows moreover B 1 − cos u − sin u 000CB a C B CB C the straightforward use of modern control-theory techniques for orbit B CB δ C B − 3∕2 u 2 sin u −2 cos u 01∕ tan i 1 CB ex C control. Three types of control algorithms are compared. The first is B CB C B CB δ C the one used in the AOK [3,11] on the PRISMA [18,19] mission. The B 0 0 0 sin u − cos u 0 CB ey C B CB C AOK controller adopts a guidance law for the orbital longitude of the B 0 sin u − cos u 000CB δ C B CB ix C ascending node and implements an analytical feedback control B CB C B − 3∕2 2 cos u 2 sin u 000CB δi C algorithm using along- and anti-along-track velocity increments. The @ A@ y A second and third controllers are the linear and the quadratic optimum 000cosu sin u 0 δu regulators from modern control theory [20,21]. This research finds its
Downloaded by STANFORD UNIVERSITY on August 2, 2015 | http://arc.aiaa.org DOI: 10.2514/1.61575 natural place and exploitation in the steadily increasing orbit-control (3) accuracy requirements for current and future Earth observation missions such as TerraSAR-X [16,22] or Sentinel-1 [23], whose nominal frozen repeat cycle orbits have very strict control require- ments within an orbital tube defined in an Earth-fixed frame. Following a thorough theoretical explanation of the concepts just mentioned, the final part of this paper is dedicated to the presentation of numerical simulations in which the different types of control paradigms are compared. As a case study, the spacecraft MANGO of the PRISMA formation is used because this research is based on the experience and flight data acquired during that mission. The control algorithm implementation used for this research is, in fact, grounded on the expertise gained in the realization of the PRISMA mission navigation and control flight software [24]. The simulation results are evaluated from a performance and operational point of view to formulate a first conclusion about the advantages and disadvantages of the control techniques considered here. The DLR, German Aerosapce Center/German Space Operations Center formation- flying simulation platform [25] is used to perform the analyses in a test environment very representative of reality. This test platform includes an orbit propagator and control flight software, and it simulates actuators and navigation errors. Fig. 1 Geometry of the relative Earth-fixed elements. DE FLORIO, D’AMICO, AND RADICE 427
Fig. 2 Earth-fixed and orbital reference frames.
Equation (3) represents the first-order solution of the Clohessy– intersection points R1 and R tR of the reference orbit projection with Wiltshire equations [31] expressed in terms of relative orbital the circle of latitude φ at times t and tR (Fig. 1). The minus sign in elements. the second member of Eq. (4) is due to the convention that δλδt is positive when δt is negative and vice versa. The time difference δt can B. Relative Earth-Fixed Elements be written as δt −δφn sin iR, where vφ n sin iR is the φ The constraints on the relative position of the subsatellite points component of the velocity of the reference subsatellite point moving ’ (SSPs) of the real and reference spacecraft and on the difference of on the Earth s surface, and n is the reference spacecraft mean motion. their altitudes determine the virtual formation’s geometry to be Hence, Eq. (4) can be written as maintained in the RTN orbital frame. Referring to Fig. 1, the λφη _ jωE − ΩRj cos φ reference frame has the origin in the subsatellite point considered δLλ δλ − cot iR δφ (7) at the altitude of the satellite, the λ axis tangent to the local circle n sin iR of latitude and pointing eastward, the φ axis tangent to the local meridian and pointing northward, and the η axis pointing along the Because the absolute orbit-control requirements are formulated in orbit radius. The relative position of the real and reference SSPs is terms of REFE [Eqs. (4–6)] but the control is realized in terms of ROE defined in the λRφRηR frame of the reference spacecraft by the [Eq. (1)], a direct mapping between the two systems is required. T phase difference vector δL δLλ; δLφ and by δh Δh∕aR, the Referring to Figs. 2b and 2c for the transformation from RTN to T altitude difference normalized to aR: λφη and using Eqs. (5–7), the REFE vector δF δLλδLφδh can be mapped into RTN coordinates with the transformation
δLλ δλδ δλδ δλ − δφ cot iR − cos φjω − Ω_ Rjδt (4) t 0 t E δF TEFδr (8)
0 1 δLφ δφ (5) _ _ jωE−ΩRj jωE−ΩRj 1 B 0 cos φ cos φ cos iR − 1 C T B n n sin iR C EF @ 0 sin iR cos iR A δh δη (6) 10 0 (9) where δλ; δφ; δη λ − λR; φ − φR; η − ηR , iR is the reference π ∈ orbit inclination, φ is the real spacecraft’s latitude at time t, ωE Equation (7) is singular for iR 0 k (with k Z) because, for −5 −1 7.292115 × 10 rad s is the Earth’s rotation rate, Ω_ R is the secular those values of iR, the orbital and equatorial planes are parallel. It is δλ rotation of the reference orbit line of nodes, δt t − tR, and tR is the interesting to remark that the term δt 0 of Eq. (4) can be written RTN δφ time at which the reference spacecraft passes at latitude φ. Figure 1 directly in the coordinate system by imposing Downloaded by STANFORD UNIVERSITY on August 2, 2015 | http://arc.aiaa.org DOI: 10.2514/1.61575 δ δ δ depicts the case in which t>tR. The phase difference vector’s rT sin iR rN cos iR 0, and substituting the solution rT −δ ∕ π component δLλ δLλ φ; δλ; δφ; δt is the distance, normalized to rN cos iR sin iR, singular for iR 0 k , in the equation δλ δ − δ aR, of the real and reference ground tracks measured along the λ axis rT cos iR rNpsin iR. Finally, from Eqs. (3) and (8) and 2 at latitude φ. δLλ is also a function of δt because the real and reference substituting cos φ 1 − sin u sin i , vector δF and its time spacecraft pass at latitude φ at different times, and the coordinates derivative d δF ∕dt, evaluated fixing u u , can be written in terms systems λRφRηR and λφη move with the SSPs while the Earth is of relative orbital elements using the transformation matrix δλ T T T T rotating. The quantity δt 0 is the normalized distance, measured EF p, with p defined in Eq. (3): λ along the axis, between the intersection points S and R1 of the real and reference orbit projections with the circle of latitude φ at time t δF u; δκ T u δκ (10) (Figs. 1 and 2). Referring to Fig. 2a, as in the coordinate system λ φ η ≈ δφ π − −δφ −δλ −δφ R R R is R1 cot iR ; , and S ; ,it δλ − ≈ δλ − δφ δλ d δF d δκ results in δt 0 R1 S cot iR. The quantity δt is the T u (11) normalized distance measured along the λ axis between the dt dt 0 1 h i 3 su τ 1−cu ciR cu B − uτ 2τsu −2τcu τciR − 1 1 τ C B 2 siR cu siR C B C T u B 3 C (12) B − usiR 2susiR −2cusiR suciR 1 − cu ciR siR C @ 2 A 1 −cu −su 000 428 DE FLORIO, D’AMICO, AND RADICE p _ 2 ~ ~ with τ jωE − ΩRj∕n 1 − sin u sin i , su sin u, cu is the Jacobian evaluated at κR of the vectors sum Ag κ Ad κ , cos u, siR sin iR, and ciR cos iR. Because u t is periodic, the and δκ is the relative orbital elements vector. vectorial function δF u; δκ t is obtained from the function Making the proper modifications to matrix A~ κR for the δF u t ; δκ t by considering only the subset u; δκ t of the normalization in aR and the introduction of δiy [Eq. (1)], and adding function’s domain u t ; δκ t . This is why the term dT∕dt δκ the control term B κ Δv: dT∕du du ∕dt δκ 0 does not appear in Eq. (11) if the variation of ϵ A κ ϵ x B κ Δv i is considered negligible. This procedure is justified by the control _ R d (20) design approach explained at the end of Sec. III.A. Equations (10) A κ A κ A κ ϵ δκ and (11) are valid under the assumptions of near-circular orbits, and where R g R d R [Eqs. (A2) and (A4)], aR is the relative orbital elements vector [Eq. (1)] normalized to the separations between the real and reference spacecraft are small when x compared to the reference orbit radius; see Eq. (3). Table 1 shows the semimajor axis, d results from the direct (dyadic) vectors product A~ κ T Δv Δ Δ Δ T form assumed by the relative Earth-fixed elements when evaluated at d R 1;aR;aR; 1; 1; 1 , and vR; vT ; vN is the the ascending node (u 0). The phase difference vector component vector of impulsive velocity increments in the RTN orbital frame. B κ δLλ at the ascending node is commonly used as operational parameter Matrix [Eq. (A5)] represents the Gauss variational equations of for the maintenance of phased orbits. If the inclination remains equal motion adapted for near-circular nonequatorial orbits [15]. The to its nominal value, the control of the phase difference at the equator elements of matrix B κ are computed with good approximation [29] will be the most effective way to monitor the displacement between using the mean orbital elements. The Gauss equations provide the real and reference ground tracks. It can also be noticed that the relationships between the impulsive velocity increments in the RTN maintenance of the altitude deviation δh requires the control of the orbital frame and the increments of the orbital elements. semimajor axis as well as the eccentricity vector. Equation (20) can be written in the general form ϵ_ AI ϵ B Δv (21) x 0A x 0 III. Virtual Formation Model: Controller Design _ d 0 d The state–space representation approach is used here for the Equation (21) is the representation of a tracking system [20] in which development of the autonomous orbit control with techniques from the atmospheric drag perturbation vector xd is represented as an modern control theory. In general, the linear model used has the form additional state variable, the disturbance input, which is assumed to satisfy the model x_ d A0xd, and I is the 6 × 6 identity matrix. If the ϵ Aϵ BΔv x _ d (13) feedback controller is designed to compute the control inputs Δvj always in the same place of the orbit (u u ), xd can be assumed as y Cϵ (14) A ≡ constant [Eq. (A3)] (i.e., 0 0). This design approach stems from the consideration that the implementation of an orbit-control strategy Δv −Gy − G x d d (15) implies the specification not only of the magnitude of the corrective maneuvers but also the in-orbit location that maximizes their ϵ A B where is the state vector; and are the dynamic model and efficiency. The choice of the less-expensive maneuver’s in-orbit x control input matrix, respectively; d is the modelled state location depends on the operational parameter that has to be Δv G perturbation; is the impulsive velocity increment vector; and controlled. An out-of-plane maneuver Δv to change the orbit’s G N and d are the output and disturbance gain matrices, respectively. inclination, for example, according to the Gauss equations [Eq. (A5)] y The output vector is composed of relative Earth-fixed elements is most effective if placed at the node (u 0) while at the orbit’s introduced by Eq. (10). highest latitude (u π∕2)ifΩ has to be changed. On the other hand, while the semimajor axis can be changed with an along-track A. Linear Dynamic Model maneuver ΔvT with the same effectiveness anywhere along the orbit, The state–space model of the orbital motion of the real and the ΔvT most effective [11] to control δLλ (Table 1) has to be reference spacecraft is given by computed at the equator (ascending or descending node) for reasons of symmetry. κ_ A~ κ A~ κ (16) g d B. Reduced Model ~ κ_R Ag κR (17) In case no eccentricity or inclination are to be actively controlled, the model can be reduced to Eq. (22) by considering only the states δ δ δ Downloaded by STANFORD UNIVERSITY on August 2, 2015 | http://arc.aiaa.org DOI: 10.2514/1.61575 where κ a; ex;ey;i;Ω;u is the mean orbital elements vector. a, iy, and u. The elements aij agij adij are given by ~ ~ Vector functions Ag and Ad [Eqs. (A1) and (A3)] describe the Eqs. (A2) and A4). The use of this model allows the design of a linear κ δ δ behavior of the mean orbital elements under the influence of the J2 controller for the relative Earth-fixed elements Lλ and Lφ. The gravitational perturbation and atmospheric drag [29,32]. The mean passive control of δh can be achieved by a proper in-orbit placement orbital elements κR of the reference spacecraft are affected only by of the along-track maneuvers as explained in the next section: the Earth’s gravitational field as they define the nominal trajectory 0 1 0 1 δ [15,16]. The linearization around κR of the difference between a11 00 a 2 B C B C Eqs. (16) and (17) yields A 3 a nJ2 B C ϵ B δ C r 2 2 @ a51 00A r aR@ iy A 4 RE 1 − e d δκ a 00 δu A~ κR δκ A~ κR (18) 61 dt d 0 1 0 p 1 02 0 μa B C B C where 1 B C A B C Br @ 0 0 sin u A xdr − CDρ@ 0 A n m ∂ A~ κ A~ κ −20− sin u∕ tan i 0 A~ κ g d R ∂κ (19) κR (22)
Table 1 Relative Earth-fixed elements (specific case)
iR u δLλ δLφ δh _ iR 0 δiy∕ sin iR jωE − ΩRj∕n δu − 2δey δu − 2δey sin iR δa − δex DE FLORIO, D’AMICO, AND RADICE 429
– ∕ ω − Ω_ ∕ In this case, Eqs. (10 12) have the following form: where the terms c12 1 sin iR, c13 j E Rj n and c21 ∕ ω − Ω_ ∕ – ar21 sin iR ar31 j E Rj n are obtained from Eqs. (23 25) δF u; ϵr aR Tr u ϵr (23) imposing u 0. Δ − δ δ ∕ The control input is given by vT g1 Lλ g2d Lλ dt aR. The objective here is to find the gains g1 and g2 that place the poles of δF A A − B GC d T A ϵ the closed-loop dynamic matrix c r r at the locations aR r u r r (24) A dt desired. The characteristic polynomial of c is
I − A 2 − 0 h i 1 js cj s s b1c21g2 ar11 s b1 ar21 c12 ar31 c13 g1 τ 1−cu ciR − 3 uτ 1 cu τ 2 B 2 cu siR C s s a^ g s a^ g (29) B C 1 2 2 1 Tr u B − 3 − C (25) @ 2 usiR 1 cu ciR siR A One of the three poles of Ac is placed in the origin regardless of the 100 value of the gains. This is due to the fact that the part of the system depending on δiy is not controllable by ΔvT. Indeed the controllability test matrix of Eq. (30) has rank 2 < 3. However, δLλ C. Linear Control can be controlled by means of variations of δu, which compensate the The linear control law for the system Eq. (21) has the general form variations of δiy (Table 1): [20] of Eq. (15). G and Gd are the gain matrices, and y Cϵ is the 0 1 ’ C b a b a2 b system s output. The terms of matrix will be computed from B 1 r11 1 r11 1 C Eq. (12) because the goal of the proposed absolute orbit controller Q 2 B C ctr BAB A B @ 0 ar b1 ar ar b1 A (30) is the maintenance of one or more relative Earth-fixed elements r r r r 21 11 21
within their control windows. For the closed-loop system to be 0 ar31 b1 ar11 ar31 b1 asymptotically stable, the characteristic roots [20] of the closed-loop p A A − BGC 2 dynamics matrix c must have negative real parts. The closed-loop poles s −a^ a^ p− 4a^ ∕2 of Eq.p (29) are G 1 1 2 This can be accomplished by a suitable choice of the gain matrix if real (or complex-conjugate) if ja^1j > 2 a^2 (or ja^1j < 2 a^2) and the system is controllable. Once the gains and thus the poles of matrix a^ b c g > 0.Ifa^ > 0, the poles are placed on the left of A G 2 1 21 1 1 c have been set, matrix d is obtained substituting Eq. (15) in the imaginary axis of the complex plane, and the closed-loop system Eq. (20), imposing y Cϵ and the steady-state condition ϵ_ 0: is stable. These stability conditions impose constraints on the value of the gains as resumed by Table 2. If the poles chosen are G CA−1B −1CA−1I − δ ∕ d c c (26) complex-conjugate, the contribution of the term g2 d Lλ dt aR to the control action ΔvT will be negligible because aRd δLλ ∕dt has −1 The linear control system is designed by means of pole placement. an order of magnitude of 30∕86; 400 m · s [11], and the value of g2 The choice of which relative Earth-fixed element has to be controlled is limited in the range indicated in the second row of Table 2. This is determined by the mission requirements, whereas the best place means that to also control d δLλ ∕dt, the poles have to be on the and direction of the orbital maneuvers is also dictated by the Gauss negative real axis because, in this way, a suitably large value of g2 can equations [Eq. (A5)]. First, an in-plane orbit-control system will be be obtained. The gains are related to the poles by the following considered with the single control input (ΔvT ) and two outputs (δLλ equation: and d δLλ ∕dt), computed at the orbit’s ascending node. This is the basic control required for the maintenance of a repeat-track orbit [15]. s s s1 s2 ar g 1 2 ;g 11 (31) Second, an in-plane/out-of-plane controller with two control inputs 1 2 b1 ar21 c12 ar31 c13 b1c21 (ΔvT and ΔvN) and three outputs (δLλ, d δLλ ∕dt, and δiy)is designed. In this case, the in-orbit phasing δLφ can be restrained in a − with s1, s2 ∈ R . The gain values chosen as a first guess are g1 control window as well. The design of these regulators will be based Δ ∕ δ Δ ∕ δ ∕ sgn c21 vTδ aR Lλ and g2 sgn c21 vT δ aRd Lλ on the reduced model [Eq. (22)]. L Δ ΔMAX δ dδL ∕ ∈ dt MAX, where vTδL , vTdδL , aR LλMAX , and aRd Lλ dt MAX R are limits imposed by design, and sgn c21 is the sign of c21. The 1. In-Plane Control with One Input and Two Outputs dynamics of the closed-loop system can be verified and adjusted by In this case, the design of the feedback system is finalized to computing the poles with Eq. (31), adjusting their placement, and Downloaded by STANFORD UNIVERSITY on August 2, 2015 | http://arc.aiaa.org DOI: 10.2514/1.61575 control δLλ and d δLλ ∕dt, computed at the ascending node iterating the process. The control input is (Table 1), by means of along-track velocity increments. This means that the orbit controller is designed to work only once per orbit at Δ −GCϵ − δ δ δ vT r g2c21 a g1c12 iy g1c13 u aR (32) most. The system components are 0 1 0 1 0 1 The following subsystem of Eq. (27) is considered for the a 00 δ r11 a b1 G B C B C B C determination of the disturbance gain matrix d [Eq. (15)]: A B C ϵ B δ C B B C r @ ar21 00A r aR@ iy A r @ 0 A d δa ar 00 δu 0 aR a δa aR b Δv x y aRδa (33) 31 0 1 dt r11 1 T d1 xd1 B C ∕ G G x B C for which Eq. (26) yields g0 1 b1, and then g1 g2 d g0 00 dr @ 0 A (27) 0 Table 2 Gain constraints ∕ with ar , ar , ar , b1 2 n, and xd1 given by Eq. (22), and g1, Pole type c21 i > 0 and g1 > 0 c21 i < 0 and g1 < 0 ∈ 11 21 31 y Cϵ p p g2 R. The output r is a 2 a^ a 2 a^ Real g > r11 2 g < r11 2 0 1 2 b c 2 b c δa 1 21 p p 1 21 δLλ ar ar 2 a^2 ar 2 a^2 ar 0 c12 c13 @ δ A 11 11 11 11 aR d δLλ aR iy (28) Complex < g2 < < g2 < c 00 b1c21 b1c21 b1c21 b1c21 dt 21 δu 430 DE FLORIO, D’AMICO, AND RADICE