Adaptive Control for Singularly Perturbed Systems

Singularly perturbed systems are a class of mathematical systems that are not well approximated by their limits and can be used to model plants with multiple fast and slow states. Multiple-timescale systems are very common in engineering applications, but adaptive control can be sensitive to timescale effects. Recently a method called [K]control of Adaptive Multiple-timescale Systems (KAMS) has shown improved performance and increased robustness for singularly perturbed systems, but it has only been studied on systems using adaptive control for the slow states. This article extends KAMS to the general case when adaptive control is used to stabilize both the slow and fast states simultaneously. This causes complex interactions between the fast state reference model and the manifold to which the fast states converge. It is proven that under certain conditions the system still converges to the reference model despite these complex interactions. This method is demonstrated on a nonlinear, nonstandard, numerical example.


I. INTRODUCTION
Singularly perturbed differential equations can be used to model systems with elements that evolve at different rates.For example, singularly perturbed models have been published for aircraft [1], spacecraft [1], electric motors [2], nuclear reactors [2], factory logistics [3], and pandemics [4].Whereas adaptive controllers have sometimes been developed for the systems listed above, these adaptive methods have largely ignored the timescale separation (termed Full-Order Adaptive Control (FOAC)) or used sequential loop closure (e.g.[5], [6]).Singular perturbation theory is a more precise method of dealing with timescale behavior, but adaptive control research to date in the literature lacks a rigorous analytical method to check for stability in the presence of singularly perturbed plants.
Singular perturbation theory is a broad mathematical field that has been used in adaptive control design [7], [8], [9], but relatively little research addresses plants that are modeled with singularly perturbed differential equations.Researchers who have used adaptive control on singularly perturbed plants have primarily applied their methods to only a subset of the states and ignored the other dynamics [4], [10], [11].This method is called Reduced-Order Adaptive Control (ROAC) and it fails when the ignored dynamics are unstable [12].
Multiple-timescale control is a branch of control theory that specifically addresses singularly perturbed plants.However, adaptive control has yet to be considered by multipletimescale control researchers.Saha and Valasek designed controllers for uncertain singularly perturbed plants [13], [14], [15], [16], but their method derives the adaption laws using a full-order Lyapunov analysis.This makes their method difficult to generalize.
This article extends the [K]control of Adaptive Multipletimescale Systems (KAMS) methodology which was first introduced and developed in [12].KAMS provides a flexible framework that enables a wide class of modern adaptive methods to be applied to singularly perturbed systems.Compared to FOAC, ROAC, and sequential loop closure KAMS is more robust and rigorous.Compared to Saha and Valasek's method, KAMS is more general.
The singularly perturbed nature of the plant causes a subset of the states to evolve significantly faster than the other states.The general premise of KAMS is to use geometric singular perturbation theory to fully decouple the fast and slow states [17].Two different adaptive controllers can then be designed in isolation for these two independent subsystems.The independent control signals are fused using a wide class of methods from the field of multiple-timescale control.These multiple-timescale control fusion techniques have not been studied in the presence of adaptive control.KAMS addresses this gap in the literature.
Allowing adaptive control in both the fast and slow states is a challenging problem because of complex interactions between the slow timescale trajectory of the fast states and the fast state reference model.The present work builds upon the author's prior work [12] which discusses the much simpler case of adaptive control for only the slow states.Unlike [12], the present work makes no prior assumptions about the stability of the plant subsystems.The novel contribution of the present work is formal proof that under certain conditions the coupling present in the more accurate full-order model is insufficient to destabilize these adaptive controllers even though they are designed in isolation.
Section II details the KAMS control framework and associated singular perturbation analysis that is used to decouple the subsystems.In Section III, a set of conditions are derived that are sufficient to show that the states converge to their reference models.Finally, in Section IV an example of KAMS on a nonlinear nonstandard system is given.This example demonstrates how methods common in the literature -Sequential Control and Adaptive Nonlinear Dynamic Inversion (ANDI) -can be used on singularly perturbed systems within the framework of KAMS.

II. CONTROL SYNTHESIS
This section introduces KAMS, explains the assumptions, and describes the notation.For more details, the reader is referred to [18], [19], [20] for adaptive control, [17], [21] for multiple-timescale control, [22] for singular perturbation theory, and [23] for differential geometry in the context of control theory.

A. SYSTEM DESCRIPTION
This work addresses singularly perturbed systems that model multiple-timescale plants.A singularly perturbed system is a system that is a function of a small scalar but not well approximated by the limit as that scalar approaches zero.This scalar is called the singular perturbation parameter.The timescale of a system is a measure of how quickly a system's states evolve.The systems considered in this article have two timescales.The slow states (x ∈ D n x x ⊆ R n x ) evolve on the slow timescale (t s ) and the fast states (z ∈ D n z z ⊆ R n z ) evolve on the fast timescale (t f ).Conversion between fast time and slow time is a change of units.Let the timescale separation parameter be the ratio of the two timescales t s /t f .It can be shown that 0 < 1.The derivative with respect to the fast timescale is denoted d (•)/dt f ( •) and the derivative with respect to the slow timescale is denoted d (•)/dt s ( •).Using the above definitions it can be shown that ( As a general rule ź x and ź ≈ x.Whereas these relationships are not always true, they provide good intuition behind the meaning of the timescale separation parameter.Multipletimescale plants can be modeled using singular perturbation theory by making the timescale separation parameter a singular perturbation parameter. This work is generalized to the class of systems which are uncertain, nonlinear, multiple-input multiple-output (MIMO) plants of the form where u ∈ R n u is the system input.This system is singularly perturbed because because 0 < 1 and the functions f x and f z are defined such that O( (1).The order of a function (i.e. the output of the O operator) is a measure of the rate of change of that function as → 0. See [17,Appendix A.2] for a more formal definition.
Remark 1: Single-timescale systems can be written in the format described by (1).However, applying a multipletimescale control technique to a single timescale system comes with a performance penalty.The resulting closed-loop responses will be slower.This effect was identified and explored in [24] and [25].Oliveira et al. demonstrated this effect on during a neuromuscular electrical stimulation experiment [26].

B. SINGULAR PERTURBATION ANALYSIS
Geometric singular perturbation theory shows that the system can be approximated by two different asymptotic solutions.The first system is found by taking the limit as → 0 and is called the reduced slow subsystem.It is only a valid approximation when t 0. Note that the fast states are constrained to a subset of their domain z where z s is the root of (2b).In multiple-timescale control, z s is called the manifold.If (2b) can be solved for z s then the system is called standard.The second asymptotic solution for (1) is found by performing a change of timescales (recall that ( •) = ( •)) and again taking the limit as → 0 This is called the reduced fast subsystem and is only a valid approximation when t is very close to 0.

C. ADAPTIVE CONTROL
The control objective of this work is to determine the input as a function of the states that permits the full-order system to track a reference model asymptotically.The first step in this process is to design two different adaptive control algorithms that stabilize the reduced subsystems individually.Many adaptive control algorithms are available in the literature for this purpose.This article addresses a wide class of algorithms that fit the format described in this section.The input to the slow subsystem is u s ∈ R n u and the input to the fast subsystems is u f ∈ R n u .The variables x m ∈ D n x x and z m ∈ D Define r x ∈ R n rx to be the bounded input to the slow state reference model and a function of time.This function and its derivative are The reference models and adaptation laws must be selected in tandem with the control input so that the control objective is achieved.The differential equations describing the evolution of the reference models and parameter estimates are of the form Remark 2: Note that r x is implicitly included as a possible input to these functions because it is fully described by time.A wide array of adaptive methods fit this format (e.g.[18], [27]).
The role of the timescale separation parameter is important in these equations.If the control input is incorrectly designed then the timescale analysis in the previous section could be invalidated.The following two assumptions are made to prevent that.
Assumption 1: The manifold is an asymptotically stable equilibrium of the fast reference model in the reduced fast subsystem.
Assumption 2: The timescale of the reference models, the slow state reference model input, and the adaptation laws all match the timescale of the subsystem to which they are applied.Mathematically this means that O( These assumptions are intuitive.For example, if the reference model for the slow states evolved on the fast timescale then the slow states would not be able to keep up -or, more precisely, their evolution could not be decoupled from the fast states.

D. MULTIPLE-TIMESCALE FUSION
The inputs to the reduced-order subsystems have been defined and will form the building blocks of the full-order system input.Let the full-order input take the form u = g u (x, x m , θx , z, z m , θz , t s ) ( 7 ) The stability analysis in the next section depends upon the reduced-order models being stabilized by their inputs u s and u f .The control objective is to select u so that a singular perturbation analysis reduces u to u s and u f .Thus the reduced-order systems are simultaneously stabilized by a single input u.Various multiple-timescale control techniques accomplish this objective by fusing the control signals for the two reduced subsystems.Three candidate methods are summarized.See [17] for more information on each of these methods.

1) COMPOSITE CONTROL
Composite Control [21, p. 94-102] selects the control input to be u = u s + u f where u f = 0 when z = z s .The engineer first selects u s so that the reduced slow model is stable.Then the engineer selects u f so that u s + u f drives the fast states to z = z s .This requires prior knowledge of the system's open-loop manifold so the system must be of standard form.Composite Control (non-adapting) has seen many industrial applications including flexible robotics [28], electrical circuits [29], and chemical reactions [30].

2) SEQUENTIAL CONTROL
In Sequential Control [17] the fast states are used as the input to the slow subsystem.The manifold is selected such that the slow states converge to their reference model by setting z s = u s Then the input u can be selected to drive the fast states to the manifold.Thus Sequential Control uses u = u f .Valasek, Narang-Siddarth, and Saha applied Sequential Control (non-adapting) to two-core coupled nuclear reactor stabilization [17, p. 122-127] and a nonlinear spring springmass-damper [31], [32].

3) SIMULTANEOUS SLOW AND FAST TRACKING
Simultaneous Slow and Fast Tracking [2] uses the input u = u s = u f to stabilize both reduced-order systems simultaneously.As such this method is not suitable for underactuated systems.The advantage of this method is that the slow states and the fast states can both be commanded to any arbitrary trajectory within the state space (constrained only by timescales and smoothness).Unlike Composite Control and Sequential Control, Simultaneous Slow and Fast Tracking allows an arbitrary manifold.Narang-Siddarth and Valasek used Simultaneous Slow and Fast Tracking (non-adapting) to control an F/A-18 A Hornet through an aggressive vertical climb and roll maneuver [2].
A block diagram for the KAMS control framework described in the previous section is given in Fig. 1.As defined in the previous section the fast adaptive control is allowed to be a function of the slow states.This uncommon case is excluded from the block diagram for readability.

III. STABILITY ANALYSIS
This section develops tools for stability analysis of the fullorder system.Whereas the adaptive controllers have been designed so that the reduced-order systems are well-behaved, these properties might not extend to the coupled full-order system.The system of equations is rewritten as a single augmented system in terms of the error coordinates for notational simplicity.Examining the differential geometric nature of the augmented system leads to the desired important insights into the behavior of the full-order system.

A. AUGMENTED ERROR DYNAMICS
Adaptive control adds additional states (i.e. the reference model and adapting parameters) to the closed-loop system.These states evolve (see (6)) and effectively create a coupled augmented closed-loop system with control states and system states.The augmented closed-loop system is defined in this section.
The variables which describe the state of the system are x, x m , θx , z, z m , and θz .For notational simplicity, these states are concatenated together.Define Note that the differential equations describing the evolution of φ are found in ( 1) and ( 6) are dependent upon the system state variables, the input, and time.However, the input is also a function of the system state variables and time (see ( 7)) so the system's dynamics are entirely described by the system's state and time.Similarly, the manifold is a function of the slow states and the input (see (2b)) so it too can be described by a function of the system state and time.Let that function and its time derivative be The following assumption is made: Assumption 3: g z s is a diffeomorphism and the manifold evolves in the slow timescale so O( f z s ) = O (1).
Section III-C discusses the manifold in more detail.
If the control objective for the full-order system is successfully achieved then two things occur as t → ∞.First z → z m → z s .This is followed by x → x m .These goals imply a set of error variables.Define where r e φ ∈ R + .Whereas x m , θx , and θz don't necessarily converge to r x , θ x , and θ z respectively their relationship is nonetheless important.Note that e z = z − zm (11) A change of variables is now performed to describe the system in terms of the error variables.The new system state variables are Let the mapping h : Equations ( 1) and ( 6) can be rewritten in terms of the new state variables as This system can be written simply as or even simpler as where f e ξ , f e η , and f e φ are defined such that ( 17) and ( 18) are identitically equal to the vector field in ( 16).This last representation obscures the timescale behavior because O( f e φ ) = O(1).Equation ( 17) is recognizable as a singularly perturbed system of the type typically studied by multiple-timescale control researchers.In fact, traditional analysis tools are applicable.However, because the form of ( 16) is available additional insights are available.Let the subscript s be used to denote a variable or vector field on the slow subsystem manifold, e.g. e η,s represents e η when z = z m = z s .Similarly, let the subscript f represent a variable or vector field on the fast subsystem manifold.The augmented reduced slow subsystem in error coordinates can be found by setting = 0 and z = z m = z s such that This reduced slow subsystem can be written simply as éξ = f e ξ ,s (e ξ , e η,s , t s ) (20) where f e ξ ,s is defined such that ( 20) is identitically equal to the vector field in (19).Equation ( 20) looks very similar to ( 17) and ( 18) but represents a fundamentally different vector field which is only defined on a subset of the full-order domain (i.e.where z = z m = z s ).The augmented reduced fast subsystem in error coordinates is This reduced fast subsystem can be written simply as or even simpler as where f e ξ , f , f e η , f , and f e φ , f are defined such that ( 22) and ( 23) are identitically equal to the vector field in (21).Finally, in the reduced fast subsystem by Assumptions 2 and 3.

B. DIFFERENTIAL GEOMETRY
Differential geometry is a natural fit for the analysis of singularly perturbed systems because the differential equations which describe these systems form nonautonomous vector fields on a topological manifold.The term manifold has been used somewhat informally and will continue to be used to refer to z s .But it is worth noting that the reduced subsystems form differential submanifolds embedded within the full-order system manifold in the topological sense.Thus it is clear that g s is a diffeomorphic chart between the full-order manifold (M) and the reduced slow manifold (M s ).The chart between the full-order manifold and the reduced fast manifold (M f ) is the trivial automorphism.The stability analysis to follow will involve the time derivative of Lyapunov functions along a vector field that is a subset of the tangent bundle of one of these topological manifolds.The notation L(•) is used to represent the Lie derivative along the vector field given in the parentheses (the traditional subscript notation is not used to ensure the subscripts on the functions are readable).Let |(•)| p be the l p norm of a vector or the induced l p norm of a matrix and let (•) p be the L p norm over time where p ∈ [1, ∞].If the L p norm is applied to a vector then it means the L p norm of each component of the vector.Unless otherwise specified, all sets are subsets of the Euclidean Hilbert space with the dimension given in the superscript.An integer subscript (•) i on a variable (not to be confused with the subscript on the p-norms) represents the ith element of the vector.

C. MANIFOLD AND THE REFERENCE MODEL
The stability proofs in the next section are significantly complicated by the relationship between the manifold and the fast reference model.Traditional multiple-timescale control and adaptive control both use a feedback loop to ensure closed-loop stability.These feedback loops still exist in the KAMS control architecture.Fig. 2 is the block diagram of KAMS from Fig. 1 except that the traditional feedback loop has been highlighted.All paths which contribute to this loop are bolded but the primary loop is blue.However, KAMS has another unconventional feedback loop because the fast reference model uses the manifold as an input (Fig. 1), the manifold is a function of the slow states (see ( 9)), the slow states are coupled with the fast states, and the control objective is for the fast states to track the fast reference model which is itself a function of the manifold.This creates a feedback loop that is typically not seen in adaptive control.Fig. 3 highlights this feedback loop.Again, all paths which contribute to this loop are bolded but the primary unconventional loop is red.
The reference model adds a complication that is not encountered in traditional multiple-timescale control.If the fast reference model is not asymptotically stable then the steady state trajectory for the slow states may not be the manifold.This calls into question the validity of the slow subsystem and means that the multiple-timescale fusion stability proofs in prior work are not applicable.These effects are unavoidable because the full-order stability analysis works by extending the stability of the reduced subsystems to the full-order system.The slow subsystem assumes that the fast states have reached their manifold, so if the stability of the reduced slow subsystem is to have any bearing on the full-order system then the fast reference model must converge to that manifold.This is the purpose of Assumption 1. Reference models are not usually asymptotically stable when their input is time-varying since they are typically Type 1 linear systems.Thus they are only capable of tracking a step input with zero steady-state error.However, closer examination reveals that Assumption 1 is not as restrictive as it appears.Recall that the manifold is assumed to evolve on the slow timescale.Equation (24c) shows that in the fast timescale zs = 0. Thus the manifold is stationary in the reduced fast subsystem, and even a type 1 reference model can be asymptotically stable.Assumption 1 is usually satisfied.However, the full-order system does not benefit from this simplification.The steady-state value of zm influences the form and function of the full-order stability proofs in Section III-D.Three cases are studied: Case 1: There exist no prior assumptions about the stability of the fast reference model in relation to the full-order manifold.This is the most general case considered, but also has the most restrictive conditions.This case often requires the control objective to be downgraded to a regulation problem.Case 2: The fast reference model is always on the manifold.This case most commonly occurs when adaptive control is not necessary for the fast subsystem.The fast control drives the fast states directly to the manifold.This type of control can be modeled by setting zm = 0 and θz = 0. Note that a parallel simplification exists where there the slow control is non-adaptive, xm = 0, and θx = 0 but this still falls within Case 1 above.Case 3: The manifold is an asymptotically stable equilibrium of the fast state reference model in the context of the full-order system.This is possible but requires an unusual reference model.This case is a slightly stricter version of Assumption 1 which only requires asymptotic stability in a subset of the domain.Remark 3: In the present work, stating that the slow subsystem does not require adaptive control will be equivalent to saying xm = 0 and θx = 0. Similarly, stating that the fast subsystem does not require adaptive control will be equivalent to saying zm = 0 and θz = 0.This terminology may be slightly misleading because there exist model-free adaptive control algorithms (e.g.[33]) and there exist non-adaptive control methods which require reference models (e.g.Feedback Linearization [34]).However, the intricacies of these methods are not in the scope of this work.

D. FULL-ORDER SYSTEM STABILITY
In this section, the stability of KAMS is analyzed in the context of the full-order system.The goal is to develop conditions that, if met, extend the stability of the reduced subsystems to the full-order system.To that end, four related theorems are proved.Each theorem belongs to one of the three cases described in the previous section.All of the theorems in this work will make use of the vector v ∈ R 4 ≥0 which is defined as

1) FOUNDATION OF REDUCED-ORDER STABILITY
The proofs in this section are similar to the proof proposed by [29], [35].However, they have been significantly altered to account for adaptive control.The general process begins by generating a composite Lyapunov function using Lyapunov functions for the reduced-order subsystems.This composite Lyapunov function is then differentiated along the vector field describing the evolution of the full-order subsystem.
Using the stability of the reduced subsystems it is shown that the differences between reduced subsystems and the full-order system are insufficient to violate the negative definiteness.This implies that e x , e z ∈ L ∞ by Lyapunov's direct method [18,Theorem 3.4.1].The following four Lyapunov functions form the basis of this approach: These Lyapunov functions are positive definite functions of class C 1 (i.e. the function and its derivative are continuous) where V e x (0, 0, t s ) = V xm (0, t s ) = V e z (0, 0, t f ) = V zm (0, t f ) = 0. Let the adaptive control for the reduced subsystems be defined such that for some α 1 , α 3 , α 4 ∈ R + .The following assumption is now made: Assumption 4: The Lyapunov functions V e x , V e z , and as needed V zm are known and exist such that ( 27) is satisfied.
Note that the existence of V zm such that (27c) holds is sufficient to guarantee that Assumption 1 is satisfied.After the Lyapunov analysis, Barbalet's Lemma is used to prove convergence [18, Lemma 3.2.5]so the following assumption is made to ensure that the conditions of Barbalet's Lemma are satisfied: Assumption 5: The functions defined in the present work are sufficiently smooth and bounded so that the function is continuously differentiable as many times as necessary.Sufficiently bounded means that, as necessary, the domain of a function being in L ∞ is sufficient to imply that the function's range is also in L ∞ .
The definitions above are a formal way of saying and indeed imply that the adaptive control for the reduced subsystems is well designed.This conclusion only applies to the reduced subsystems.

2) CASE 1
There exist no prior assumptions about the stability of the fast reference model in relation to the full-order manifold.
Theorem 1: Let the matrix K = K T be defined as If ∃d ∈ (0, 1) and d * (1 − d ) such that K is positive definite, then e x , e z → 0 as t → ∞.
Proof: Define a composite Lyapunov function ) Conceptually this is the derivative in the subsystems plus some errors due to inaccuracies caused by the model reduction.Rearranging gives Some of these terms can be simplified because each Lyapunov function is not a function of all state variables (i.e. its partial derivative is zero).In doing so, must be carefully accounted for.
Some of the vector fields are the same in the reduced subsystem and the full-order subsystem, which allows further simplification Substituting the conditions from ( 27) and ( 28) gives: The triangle inequality shows that Using this and rearranging gives V ≤ −v T Kv (37) Thus, by Lyapunov's direct method e φ ∈ L ∞ .The goal is now to show that the conditions of Barbalat's Lemma are satisfied.This is done by showing that e x , e z , éx , éz ∈ L ∞ and e x , e z ∈ L 2 .By the arguments in Table 1 it can be concluded that these conditions are met.Note that the order of the lines in this table is significant.Thus, via Barbalat's Lemma, it is known that e x , e z → 0 as t → ∞.Corollary 1: Let the plant exist such that the reduced slow subsystem does not require adaptive control (i.e.xm = 0 and θx = 0).Assume that conditions (28b), (28c), and (28d) of Theorem 1 are true.Let the matrix K = K T be defined as If ∃d ∈ (0, 1) and d * (1 − d ) such that K is positive definite, then e x , e z → 0 as t → ∞.
Each of the following proofs assumes that f z − f z, f = 0 which occurs when does not appear on the right side of equation (1b).This is very common and making this assumption will aid in interpreting the results.

3) CASE 2
The fast reference model is always on the manifold.
Corollary 2: Let the plant exist such that the reduced fast subsystem does not require adaptive control (z m = 0 and θz = 0) and does not appear on the right side of (1b).Assume that condition (28b) of Theorem 1 is true.Then ∀ it is true that e x , e z → 0 as t → ∞.
Proof: The proof proceeds exactly as Theorem 1 except for γ = 0 and zm = 0.This reduces the matrix where V xm has been dropped because all of the cross terms of xm have been removed.This is simple enough for additional conclusions.By Sylvester's Criterion K is positive definite if and only if the leading principle minors (LPM) are positive [37].This gives rise to the following two inequalities which, if satisfied, imply that K is positive definite.
Inequality (40a) is satisfied by definition.Rearranging inequality (40b) gives when β = 0.When β = 0 then inequality (40b) is satisfied by definition.Recall that d is arbitrary so ∀ ∃d such that inequality (41) is satisfied.Continuing with Barbalet's Lemma as in Theorem 1 gives that e x , e z → 0 as t → ∞.

4) CASE 3
The manifold is an asymptotically stable equilibrium of the fast state reference model in the context of the full-order system.
Corollary 3: Assume that does not appear on the righthand side of (1b).Assume that condition (28b) of Theorem 1 is true.If ∃α 4 ∈ R + such that then ∀ it is true that e x , e z → 0 as t → ∞.
Proof: The proof largely follows Theorem 1 except for V xm is removed from the composite Lyapunov function.Before reaching (35), (34) can be rewritten using L( This is simple enough for additional conclusions.By Sylvester's Criterion K is positive definite if and only if the LPMs are positive.This gives rise to the following three inequalities which, if satisfied, imply that K is positive definite. Inequality (44a) is satisfied by definition.Rearranging the other two inequalities gives when β = 0.When β = 0 then inequalities (44b) and (44c) are satisfied by definition.Recall that d is arbitrary and ∀ ∃d such that the inequalities in (45) are satisfied.Continuing with Barbalet's Lemma as in Theorem 1 gives that e x , e z → 0 as t → ∞.Remark 4: Assumption 2 places bounds on the acceptable range of the adaptation gains.Note that is not required to implement the control.This is advantageous because can be difficult to determine.However, a rough approximation of allows the engineer to design the adaptive laws and reference models so that they evolve on the correct timescale.Beyond these conditions, is allowed to be uncertain.
Remark 5: The condition that K be positive definite limits the range of acceptable timescale separation parameters (e.g.see [29]).However, for Theorem 1 and Corollary 1 the analytical bounds would be complex.Remark 7: Systems which use adaptive control are likely to be nonstandard because adaptive control is specifically designed for systems with model uncertainties.Thus it is common for the open-loop manifold to be uncertain even if the system is standard in the traditional sense.Let the term uncertain nonstandard refer to this condition.Recent multiple-timescale control research has addressed nonstandard systems [17].Both Sequential Control and Simultaneous Slow and Fast Tracking are nonstandard methods because the manifold is specified.By comparison, Composite Control requires the open-loop manifold to be known apriori, so the manifold must be measured or analytically available.Thus Composite Control is well suited for systems that do not require adaptive control in the fast subsystem.

5) SUMMARY OF THEOREMS
This section describes each of the theorems that are proven in this article and provides criteria that can be used to determine which of the theorems apply to a given system.Theorem 1 is the most general but also has the most restrictive conditions on stability.It requires that the slow reference model be asymptotically stable to the reference model input.In practice, this can often limit the theorem to regulation.Three special cases of Theorem 1 were studied that are less restrictive.Corollary 1 is applicable when adaptive control is only used for the fast subsystem.Corollary 2 is applicable when adaptive control is only used for the slow subsystem.Corollary 3 allows adaptive control in both subsystems, but the manifold must be an asymptotically stable equilibrium of the fast reference model in the context of the full-order system.
Theorem 1 and Corollary 1 both allow the timescale separation parameter to appear on the right side of the fast states' differential equations and require checking the positive definiteness of a matrix.Corollaries 2 and 3 do not.KAMS typically requires differentiation of the manifold.In Theorem 1 and Corollary 1 the derivative of the manifold is used to ensure that condition (28d) is satisfied.The derivative of the manifold is not explicitly required for Corollary 3, but it is often required to ensure the manifold is an asymptotically stable equilibrium of the fast reference model.It is therefore significant that Corollary 2 does not require differentiating the manifold.

IV. VALIDATION
An example demonstrates and validates this method.Consider the following nonlinear, nonstandard, uncertain dynamical system where θ ∈ R + is an uncertain parameter.The control objective is for x to track the reference model

D. ALTERNATIVE APPROACH
Corollary 1 is also applicable because adaptive control is only required for the fast subsystem.To demonstrate this the problem is revised to a regulation problem and the fast reference model is redefined so that it is no longer asymptotically stable about the manifold zm = −a z zm (63) Note that the manifold is still a stable equilibrium.It is even asymptotically stable when the manifold is constant with respect to time, but it is not asymptotically stable in the context of the full-order system.Reference models such as this are useful if ŕx is not known apriori or the manifold is difficult to differentiate.From (57) it can be shown that where the domain has been restricted to x < 1 and z < 1.Note that δ 1 = 2k x + k 2 x , δ 2 = 0, δ 3 = k x , and δ 4 = k x .Substituting the values from the previous numerical example gives From Fig. 7 it can also be seen that when = 0.1 (as in the previous example) ∃d such that all of the LPMs are    positive.Thus by Corollary 1 e x , e z → 0 as t → ∞.The time evolution of the slow state for this alternative approach is shown in Fig. 8; the fast state is shown in Fig. 9; and the time evolution of the adapting gain is shown in Fig. 10.This example demonstrates how the manifold evolves on the slow timescale per Assumption 3.
Remark 8: As mentioned previously a common approach to these problems is to apply sequential loop closure to the subsystems.This example uses Sequential Control.To clarify Sequential Control is an extension of sequential loop closure.Applying sequential loop closure to this example would yield numerically indistinguishable results.Sequential Control was first published by Narang-Siddarth and Valasek in [17] where they showed that singular perturbation techniques could be used to rigorously show stability and obtain specific bounds on the time scale separation parameter.By extension, this same advantage is available to KAMS.The insights from Fig. 7 are a unique contribution of KAMS that is not available to traditional adaptive sequential loop closure implementations.Furthermore, unlike sequential loop closure, KAMS allows the use of Composite Control and Simultaneous Slow and Fast Tracking.
Remark 9: The MATLAB code used for both of the examples has been made open source and is available on Code Ocean [38].

V. CONCLUSION
This article extended the [K]control of Adaptive Multipletimescale Systems (KAMS) methodology to singularly perturbed systems with adaptive control in both the fast and slow subsystems; a wide class of adaptive control and multipletimescale control methods fit within this framework.Sufficient conditions for asymptotic stability were proven and coupling effects between the manifold and the fast reference model were identified.The stability of the full-order system was connected to the stability of the reduced-order systems through Theorem 1 and its corollaries.A nonlinear nonstandard system was used to demonstrate KAMS.
This article identified complex interactions between the fast reference model and the manifold which occur when adaptive control is used to stabilize the fast subsystem.This makes traditional multiple-timescale control proofs insufficient when adaptive control is used in the fast subsystem.The theorems proved in this article account for these complex interactions by carefully formatting the augmented error dynamics and by judiciously selecting sufficient conditions.The primary limitation of KAMS is the requirement to verify the conditions given in Theorem 1 and its corollaries.These conditions restrict the set of systems to which the theorems in this article can be applied.Lyapunov functions may not be known for some systems and adaptive control methods.Suitable Lyapunov functions are known for several popular adaptive control methods (e.g.Model Reference Adaptive Control and Adaptive Nonlinear Dynamic Inversion).Another limitation is that many applications will require differentiation of the manifold.This can be a complicated calculation, but it can sometimes be avoided by judicious selection of control objectives, careful system modeling, and the use of the correct corollary.See the alternative approach example above for a demonstration of this.Based upon the results presented in the article KAMS is judged to be a feasible control approach for uncertain nonstandard singularly perturbed systems regardless of which subsystem (fast, slow, or both) the uncertainty appears in.Further, KAMS is more capable than traditional sequential loop closure because it can be used to determine the minimum allowable timescale separation and it allows for the use of Composite Control and Simultaneous Slow and Fast Tracking.
There are several potential avenues for future research.First, future research could consider adapting laws that do not adapt in the same timescale as the subsystem to which they are applied.Second, future research could determine alternate formats for the upper bounds in Theorem 1. Finally, experimentally validating the performance of KAMS on physical systems would be insightful.
states.The parameters θx ∈ P n θx θ x ⊆ R n θx and θz ∈ P n θz θ z ⊆ R n θz are adaptive estimates of the true parameters θ x and θ z respectively.The true parameters are allowed to be time-varying.Let

Remark 6 :
Corollaries 1 and 2 study the case where only one subsystem requires adaptive control.If neither subsystem requires adaptive control then Theorem 1 reduces to [35, Theorem 1].

FIGURE 6 .
FIGURE 6. Evolution of the adapting gain.

FIGURE 7 .
FIGURE 7. Effects of varying and d on the applicablility of Corollary 1.

FIGURE 8 .
FIGURE 8. Evolution of the slow state for the alternative approach.

FIGURE 9 .
FIGURE 9. Evolution of the fast state for the alternative approach.

FIGURE 10 .
FIGURE 10.Evolution of the adapting gain for the alternative approach.

TABLE 1 . Proof that the conditions of Barbalat's Lemma are met.
From Corollary 1 it is known that if ∃d ∈ (0, 1) and d * (1 − d ) such that K is positive definite then e x , e z → 0 as t → ∞.By Sylvester's Criterion K is positive definite if and only if the LPMs are positive.The first LPM is positive by definition.The second and third LPMs depend upon and d.Fig.7plots this relationship.If both LPMs are positive for a given then the conditions of Corollary 1 are satisfied.In Fig.7note that as increases there is a point after which d such that both LPMs are positive simultaneously.At this point, the timescale separation is insufficient for Corollary 1 to guarantee convergence.