It has been known for some time that proportional output feedback will stabilize MIMO, minimum-phase, linear time-invariant systems if the feedback gain is sufficiently large. High-gain adaptive controllers achieve stability by automatically driving up the feedback gain monotonically. More recently, it was demonstrated that sample-and-hold implementations of the high-gain adaptive controller also require adaptation of the sampling rate. In this paper, we use recent advances in the mathematical field of dynamic equations on time scales to unify and generalize the discrete and continuous versions of the high-gain adaptive controller. We prove the stability of high-gain adaptive controllers on a wide class of time scales.
1. Introduction
The concept of high-gain adaptive feedback arose from a desire to stabilize certain classes of linear continuous systems without the need to explicitly identify the unknown system parameters. This type of adaptive controller does not identify system parameters at all, but rather adapts the feedback gain itself in order to regulate the system. A number of papers examine the details of various kinds of high-gain adaptive controllers [1–4], among others. More recently, several papers have discussed one particularly practical angle on the high-gain adaptive controller, namely, how to cope with input/output sampling. In particular, Owens [5] showed that it is not generally possible to stabilize a linear system with adaptive high-gain feedback under uniform sampling. Thus, Owens et al., develop a mechanism to adapt the sampling rate as well as the gain, a notion subsequently improved upon by Ilchmann and Townley [6, 7], Ilchmann and Ryan [8], and Logemann [9].
In this paper, we employ results from the burgeoning new field of mathematics called dynamic equations on time scales to accomplish three principal objectives. First, we use time scales to unify the continuous and discrete versions of the high-gain controller, which have previously been treated separately. Next we give an upper bound on the system graininess to guarantee stabilizability for a much wider class of time scales than previously known, including mixed continuous/discrete time scales. Third, the paper represents the first application of several fairly recent advances in stability theory and Lyapunov theory for systems on time scales, and two new lemmas are presented in that vein. We also give a simulation of a high-gain controller on a mixed time scale.
2. Background
We first state two assumptions that are required in the subsequent text.
The system model and feedback law are given by the linear, time-invariant, minimum phase system
ẋ(t)=Ax(t)+Bu(t),y(t)=Cx(t),x(0)=x0,u(t)=-k(t)y(t)
for t≥0. System parameters A∈ℝn×n, B,CT∈ℝn×m, x0∈ℝn, and n are unknown. The feedback gain k:ℝ→ℝ+ is piecewise continuous, and nondecreasing as t→∞. By minimum phase, we mean that the polynomial
det[A-λIBC0]
with λ∈ℂ is Hurwitz (zeros in open left hand plane).
Furthermore,
(CB)T+(CB)>0,
that is, (CB)T+(CB) is positive definite. (In [3], it is pointed out that a nonsingular input/output transformation T always exists such that B̃=BT-1 and C̃=TC give (C̃B̃)T+(C̃B̃)>0.)
Under these conditions, it has been known for some time (e.g., [1]) that there are a wide class of gain adaptation laws k(t)=f(y(t)), f:ℝm→ℝ, that can asymptotically stabilize system (2.1) in the sense that
y∈L2[t0,∞),limt→∞k(t)<∞.
Subsequently, various authors [1, 5, 6, 10] assumed that the output is obtained via sample-and-hold, that is, yi:=y(ti) with ki:=k(ti) and i∈ℕ0. Thus it becomes necessary also to adapt the sample period hi:=ti+1-ti, so the closed-loop control objectives are
yi∈ℓ2[i0,∞),limi→∞ki<∞,limi→∞hi>0.
Though several variations on these results exist, these remain the basic control results for continuous and discrete high-gain adaptive controllers. The continuous and discrete cases have previously been treated quite differently, but we now construct a common framework for both using time-scale theory.
3. A Time-Scale Model
The system of (A1) can be replaced by
xΔ(t)=Â(t)x(t)+B̂(t)u(t),y(t)=Cx(t),x(0)=x0,u(t)=k(t)y(t),t∈𝕋,
where 𝕋 is any time scale unbounded above with 0∈𝕋. With a series expansion similar to [1], we see that
Â(t):=expc(μ(t)A)A,B̂(t)=expc(μ(t)A)B,
where expc is the matrix power series function
expc(X):=I+12X+16X2+⋯+1n!Xn-1+⋯
and μ is the time scale graininess. Implementing control law, (2.2) then gives
xΔ(t)=𝒜(t)x(t),𝒜(t):=expc(μ(t)A)(A-k(t)BC).
Note that Â,B̂, and 𝒜 may all be time-varying, but we will henceforth drop the explicit reference to t for these variables. For future reference, we also note that if supt∈𝕋μ(t)<∞, then ∥Â∥ and ∥B̂∥ are also bounded (c.f. Appendix).
The design objectives are to find graininess μ and feedback gain k as functions of the output y,
μ(t)=g(y(t)),g:ℝm→[0,∞),k(t)=f(y(t)),f:ℝm→ℝ+
with t∈𝕋 and k∈Crd(𝕋) nondecreasing, such that
y∈L2[t0,∞)𝕋:={w:∫t0∞|w(t)|2Δt<∞},limt→∞k(t)<∞.
It is important to keep in mind the generality of the expressions above. A great deal of mathematical machinery supports the existence of delta derivatives on arbitrary times scales, as well as the existence and characteristics of solutions to (3.4). See, for example, [11, 12].
4. Stability Preliminaries
We begin this section with a definition and theorem from the work of Pötzsche et al. [13].
Definition 4.1.
The set of exponential stability for the time-varying scalar equation zΔ(t)=λ(t)z(t), z0:=z(t0) with λ:𝕋→ℂ and z,λ∈Crd(𝕋), is given by
𝒮(𝕋):=𝒮ℂ(𝕋)∪𝒮ℝ(𝕋),
where
𝒮ℂ(𝕋):={η(t)∈ℂ:α=-lim supt→∞1t-t0∫t0tlog|1+μ(t)η(t)|μ(t)Δt>0},𝒮ℝ(𝕋):={η(t)∈ℝ:∀t∈𝕋,∃τ>twithτ∈𝕋such that1+μ(t)η(t)=0},
with η∈Crd(𝕋) and arbitrary t0∈𝕋.
Theorem 4.2 (see [13]).
Solutions of the scalar equation zΔ(t)=λ(t)z(t) are exponentially stable on an arbitrary 𝕋 if and only if λ(t)∈𝒮(𝕋).
We note here that Pötzche, Siegmund, and Wirth did not explicitly consider scenarios where η is time-varying, but their stability analysis remains unchanged for η(t).
The set 𝒮ℝ(𝕋) contains nonregressive eigenvalues λ(t), and a loose interpretation of 𝒮ℂ(𝕋) suggests that it is necessary for a regressive eigenvalue to reside in the area of the complex plane where ∥1+μ(t)λ(t)∥<1 “most” of the time. The contour ∥1+μ(t)λ(t)∥=1 is termed the Hilger Circle. Since the solution of zΔ(t)=λ(t)z(t) is x(t)=z0eλ(t,t0), Theorem 4.2 states that if λ(t)∈𝒮(𝕋) then some K(t0)>1 exists such that
∥z(t)∥=∥z0∥∥eλ(t,t0)∥≤∥z0∥Ke-α(t-t0),
where eλ is a generalized time scale exponential. The Hilger Circle will be important in the upcoming Lyapunov analysis, as will the following lemmas.
Lemma 4.3.
Let z∈Crd(𝕋) be a function which is known to satisfy the inequality zΔ(t)≤λ(t)z(t), z(0)=z0∈ℝ with λ:𝕋→ℝ and λ∈Crd(𝕋). If z(t)>0 for all t≥t0∈𝕋, then λ(t)∈ℛ+.
Proof.
Defining v(t)=zΔ(t)-λ(t)z(t) gives rise to the initial value problem
zΔ(t)-λ(t)z(t)=v(t),z(0)=z0,t∈𝕋,
where v(t)≤0.
First, suppose λ(t) is negatively regressive for t≥t0, that is, λ(t)<-1/μ(t) with μ(t)>0. Then (4.5) yields zσ=μλz+z+μv=z(1+μλ)+μv<0.
On the other hand, suppose λ(T) is nonregressive for some T>t0,T∈𝕋. If λ(t)<-1/μ(t) over T>t≥t0, then invoke the preceding argument. If λ(t)>-1/μ(t) over T>t≥t0, then solve (4.5) to get
z(t)=eλ(t,t0)z0+∫t0teλ(t,σ(τ))v(τ)Δτ,T>t≥t0.
Since eλ(t,t0)>0 for t<T, we see that ∫t0teλ(t,σ(τ))v(τ)Δτ<0. However, for t≥T, (4.6) becomes
z(t)=0+∫t0ρ(T)eλ(t,σ(τ))v(τ)Δτ,t≥T.
Thus, z(t)<0 for all t≥t0 for both negatively regressive and nonregressive λ(t), a contradiction of the Lemma's presupposition. This leaves only λ(t)∈ℛ+.
At this point we pause briefly to discuss Lyapunov theory on time scales. DaCunha produced two pivotal works [14, 15] on solutions P of the generalized time scale Lyapunov equation,
A(t)TP(t)+P(t)A(t)+μ(t)AT(t)P(t)A(t)=-Q(t),t∈𝕋,
where A(t),P(t),Q(t)∈ℝn×n, A and Q are known and Q(t)>0. Though it will not be necessary to solve (4.8) in this work, we will see that the form of (4.8) leads to an upper bound on the graininess that is generally applicable to MIMO systems, an advancement beyond previous works which gave an explicit bound only for SISO systems. DaCunha proves that a positive definite solution P(t) to the time scale Lyapunov equation with positive definite Q(t) exists if the eigenvalues of A(t) are in the Hilger circle for all t∈𝕋. Furthermore, P(t) is unique. As with the well-known result from continuous system theory (c.f. [16]), the solution is constructive, with
P(t)=∫𝕊tΦA(t)T(s,0)Q(t)ΦA(t)(s,0)Δs,
where ΦA(t,t0) denotes the transition matrix for the linear system xΔ(t)=A(t)x(t), x(t0)=I. The correct interpretation of this integral is crucial: for each t∈𝕋, the time scale over which the integration is performed is 𝕊t:=μ(t)ℕ0, which has constant graininess for each fixed t. (Since this paper was written, we have explored the subject on algebraic and dynamic Lyapunov equations in more detail; we refer the interested reader to [17].)
Before the next lemma, we define
𝒦:={k:𝕋→ℝ+∣k∈Crd(𝕋),kΔ(t)≥0∀t∈𝕋,limt→∞k(t)=∞}.
The next lemma follows directly.
Lemma 4.4.
Given assumptions (A1) and (A2) and k∈𝒦, there exists a nonzero graininess μ̅(t) and a time t* such that, for all μ(t)≤μ̅(t) and t>t*, the matrix -kCB̂ satisfies a time scale Lyapunov equation with P=I, Q(t)≥ε2 for small ε2>0, and B̂ from (3.2).
Proof.
We construct μ̅(t) as
μ̅(t):=1k(t)(λmin{CB+(CB)T}λmax{(CB)TCB}-ε1),t∈𝕋,
with ε1>0 sufficiently small so that μ̅(t)>0 on 𝕋. This holds if
(CB)T+CB-μ̅(t)k(t)(CB)TCB=ε1I.
Multiplying (4.12) by -k(t) gives
(-k(t)CB)T+(-k(t)CB)+μ̅(t)(-k(t)CB)T(-k(t)CB)=-k(t)ε1I.
(We now drop the explicit time-dependence for readability.) Set
Σ(μ̅):=expc(μ̅A)-I,
so that B̂(μ̅)=[I+Σ(μ̅)]B, yielding
(-kCB̂)T+(-kCB̂)+μ̅(-kCB̂)T(-kCB̂)=-k(ε1I-Z),
where each term of Z is a product of contants times Σ(μ̅). Since ∥Σ(μ̅)∥→0 as μ̅→0 (c.f. Appendix), there exists a time t*∈𝕋 when ∥Z∥<ε1. Because the preceding arguments admit any graininess μ(t)≤μ̅(t), it follows that
(-kCB̂)T+(-kCB̂)+μ(-kCB̂)T(-kCB̂)=-kQ,t>t*,μ(t)≤μ̅(t),
where Q(t)≥ε2 and ε2:=ε1-∥Z∥>0.
We comment briefly on the intuitive implication of Lemma 4.4. Equation (4.16) shows that there exist positive definite Q(t) and P=I to satisfy an equation of the form (4.8). According to DaCunha, this implies that the eigenvalues of -kCB̂ lie strictly within the Hilger circle for t>t*.
5. System Stability
We now come to the three central theorems of the paper. If BC is not known to be full rank, or cannot be full rank because of the input/output dimensions, then it must be assumed (or determined a priori) that the eigenvalues of A-k(t)BC attain negative real parts at some point in time. This phenomenon is investigated in-depth by other authors [1, 4]. We then make use of the observation by Owens et al. [3] that there must exist some k*>0 such that if k(t)≡k* the system of (A1) has a positive-real realization. This, together with the Kalman-Yakubovich Lemma [16], implies existence of P,Q>0 so that
(A-k*BC)TP+P(A-k*BC)=-Q,PB=CTCB.
Theorem 5.1 (Exponential Stability).
In addition to (A1) and (A2), suppose
t∈𝕋 where 𝕋 is a time scale which is unbounded above but with μ(t)≤μ̅(t),
k∈𝒦 (implying from (4.11) that μ(t)→0, but not necessarily monotonically),
BC is not necessarily full rank, but there exists a time t*∈𝕋 such that the eigenvalues of (A-k(t)BC) are strictly in the left-hand complex plane for t≥t*.
Then the system (3.1), is exponentially stable in the sense that there exists time t0∈𝕋 and constants K(t0)≥1, α≥0, such that
∥x(t)∥≤∥x(t0)∥Ke-(1/2)α(t-t0),t≥t0.
Proof.
Set k*:=k(t*). Then assumption (iii) is the prerequisite for (5.1). Again, we suppress the time-dependence of x, k, μ, Â, and B̂. Similarly to Lemma 4.4, terms containing Σ(μ) may be added to the first equality in (5.1) to obtain
(Â-k*B̂C)TP+P(Â-k*B̂C)≤-ε3,t≥t**,
for some small ε3>0. Note t=t** is the point at which terms involving Σ(μ) become small enough for (5.3) to hold. Defining Z1:=CTCΣB-PΣB, the second equality in (5.1) gives
PB̂=CTCB̂+Z1.
Consider the Lyapunov function V=x(t)TPx(t) with P from (5.1). Then, using (5.1), (5.3) and Lemma 4.4,
VΔ=xTΔPx+xσTPxΔ=xT[𝒜TP+P𝒜+μ𝒜TP𝒜]x=xT[(Â-k*B̂C)TP+P(Â-k*B̂C)-(k-k*)(CTB̂TP+PB̂C)+μ(Â-kB̂C)TP(Â-kB̂C)]x≤xT[-ε3-(k-k*)(CTB̂TCTC+CTCB̂C+2Z1C)+μÂTPÂ-2μkÂTPB̂C+μk2CTB̂TCTCB̂C+μk2CTB̂TZ1C]x=xT[-ε3+μÂTPÂ]x+yT[-k(CB̂)T-kCB̂+μk2(CB̂)TCB̂+μk2B̂TZ1+k*((CB̂)T+CB̂)]y+xT[-(k-k*)2Z1-2μkÂTPB̂]y≤-(ε3-μ∥ÂTPÂ∥)xTx-(kε2-μk∥B̂∥∥kZ1∥-2k*∥CB̂∥)yTy+(2∥kZ1∥+2μk∥ÂTPB̂∥)xTy.
At this point, we observe the following.
∥ÂTPÂ∥ is bounded because  is bounded. Set γ1:=supt∈𝕋∥ÂTPÂ∥.
μk is bounded by assumption (i) and (4.11).
∥kZ1∥ is bounded because k is proportional to 1/μ and (1/μ)∥Σ(μ)∥→const as μ→0. Set γ2:=supt∈𝕋2∥kZ1∥.
Set γ3:=supt∈𝕋μk∥B̂∥γ2.
Set γ4:=supt∈𝕋2k*∥CB̂∥.
Set γ5:=supt∈𝕋2μk∥ÂTPB̂∥.
Recalling the standard inequality xTy≤βxTx+(1/β)yTy for any β>0, we then have
VΔ≤-(ε3-μγ1-βγ2-βγ5)xTx-(kε2-γ3-γ4-1βγ2-1βγ5)yTy.
By assumption (ii), there exists a time t0≥max{t*,t**} such that, for sufficiently small β,
VΔ≤-(ε3-μγ1-βγ2-βγ5)γmin(P)V:=η(t)V,t≥t0,
with η(t)<0. Then, by [11, Theorem 6.1], Theorem 4.2, and Lemma 4.3 it follows that there exists K(t0)≥1 so that
∥x(t)∥≤K∥x(t0)∥e-(1/2)α(t-t0),t≥t0,
where α is defined in (4.2).
We point out that, when BC is full rank, assumption (iii) above is no longer necessary as there always exists a k* such that the eigenvalues of (A-k(t)BC) are strictly real-negative for k(t)≥k*. One more lemma is required before the next theorem.
Lemma 5.2 ..
If 𝕋 is a time scale with bounded graininess (i.e., μ∞:=supt∈𝕋μ(t)<∞), then
c1∫t0∞eαtdt≤∫t0∞eαtΔt≤c2∫t0∞eαtdt,
where c1, c2, α∈ℝ and c1, c2>0.
Proof.
Consider the case when α>0. The process of time scale integration is akin to the approximation of a continuous integral via a left-endpoint sum of (variable width) rectangles. If the function to be summed is increasing (as in this case), the sum of rectangular areas will be less than the continuous integral, meaning c2=1, c1<1. One estimate of the lower bound, then, follows by simply increasing c1 until c1eαt meets one of the rectangle right endpoints which are given by eαρ(t). Thus c1eαt≤eαρ(t), or equivalently, c1eασ(t)≤eαt. This in turn yields c1≤e-αμ(t). Therefore, the most conservative bound is given by c1≤e-αμ∞. The case for α<0 can be argued similarly, leading to the lemma's following conclusion:
c1={e-αμ∞,α>0,1,α≤0,c2={1,α≥0,e-αμ∞,α<0.
We are now in a position to state the main theorem of the paper.
Theorem 5.3.
In addition to (A1) and (A2), assume the prototypical update law, kΔ(t)=∥y(t)∥2 with k0:=k(0)>0. Then limt→∞k(t)<∞ and y∈L2[t0,∞)𝕋.
Proof.
For the sake of contradiction, assume k(t)→∞ as t→∞. Then k∈𝒦. Theorem 5.1 yields x∈L∞[t0,∞)𝕋 and therefore y∈L∞[t0,∞)𝕋. The solution for k(t) is (by [11, Theorem 2.77]),
k(t)=k(t0)+∫t0t∥y(τ)∥2Δτ≤k(t0)+∥C∥2∫t0t∥x(τ)∥2Δτ,t≥t0.
In conjunction with Lemma 5.2, this allows
limt→∞k(t)≤k(t0)+∥C∥2∫t0∞∥x(t)∥2Δt≤k(t0)+∥C∥2∥x(t0)∥2K2∫t0∞e-α(t-t0)Δt=k0+∥C∥2∥x(t0)∥2K2∫0∞e-αtΔt<∞.
This contradicts the assumption, so it must be that k(t)<∞ for t∈𝕋. It also immediately follows that ∫t0∞∥y∥2Δt<∞.
It seems possible that Theorem 5.3 may be improved to show that the output is convergent, that is, y(t)→0 as t→∞. This is left as an open problem.
6. Discussion
We remark here that there is a great amount of freedom in the choice of the update law for k(t) (c.f. [5]). We use the simplest choice for convenience (as do most of authors); the essential arguments remain unchanged for other choices. There is also freedom in the choice of update for μ(t). The expression (4.12) can be simplified to
(CB)T+CB-μ(t)k(t)(CB)TCB>0,t∈𝕋,
for any μ(t)≤μ̅(t). In the SISO case, this further reduces to the expression derived by Owens [5], that μkCB<2. It requires the graininess (which may interpreted as the system sampling step size for sample-and-hold systems) to share at least an inverse relationship with the gain, but is otherwise quite unrestrictive. Ilchmann and Townley [7] note that μ(t)=1/(k(t)logk(t)) meets μkCB<2 after sufficient time without knowledge of CB. While previous works have always constructed a monotonically decreasing step size, the time scale-based arguments in this paper reveal even greater freedom: μ(t) may actually increase, jump between continuous, and discrete intervals, or even exhibit bounded randomness. Two examples of the usefulness of this freedom come next.
For the first example, we point out that the notation in the previous sections somewhat belies the fact that the system's time domain (its time scale) may be fully or partially discrete, and thus there is no guarantee in Theorem 5.3 that the output has stabilized between samples. As pointed out in [7] and elsewhere, a sampled system with period μ is detectable if and only if
λk-λl2πjμ∉ℤfor anyλk≠λl,j=-1,
where λk,λl∈{0,spec(A)}. We next comment on how to circumvent the intrasample stabilization problem. Both of the following methods essentially permit the graininess μ to “wiggle” a bit so that an output sample must eventually occur away from a zero-crossing. Recalling that μ(t)≤μ̅(t) from Theorem 5.1(i), let μ(t)=μ̅(t)v(t), where 0<v≤1is one of the sequences below:
(1) Let v be consist of an infinitely repeated subsequence with n!+1 elements that are random numbers between 0 and 1. Let these elements, labeled v1,v2,…, be irrational multiples of each other. Assume μ(t) converges to μ̅v1>0 such that y(t)=0 but the true continuous output is nonzero. This implies that there exist integers k,l such that (λk-λl)μ̅v1/2πj∈ℤ. As sequence v advances, there may be at worst n! combinations of j,k such that (λj-λk)μ̅vr/2πj∈ℤ for r=1⋯n!. However, at the next instant in time, (λk-λl)μ̅vn!+1/2πj must be irrational and therefore not in ℤ. The controller will detect a nonzero output and continue to adapt k(t) and μ(t). In practice, of course, it is not possible to obtain a sequence of truly irrational numbers, but most modern computer controllers have enough accuracy to represent the ratio of two very large integers, so that this technique would only fail for impractically high magnitudes of ∥λk-λl∥.
(2) Let v be a sequence of random numbers in the specified range. Even in a computer with only 8-bit resolution for v, the probability of (λ1-λ2)μ̅v/2πi∈ℤ drops drastically after a few sample periods.
We remark here that if (6.2) holds, then (𝒜,C) is detectable because (A,C) is detectable. Thus, the stability of y implies the stability of x. We do not dwell on this here, but see, for example, [7] for a similar argument.
For the second example, we consider a problem posed by distributed control networks (c.f. [18, 19]). Here, one communication network supports many control loops as well as a certain volume of unrelated high-priority traffic. The (unpredictable) high-priority traffic may block the control traffic at times, forcing longer-than-anticipated sample periods. In normal operation, the controller may sample fast enough to behave, for all practical purposes, like a continuous control with μ(t)=0. At the blocking instant, tb, μ(tb) rises to some unpredictable level. The scheduling question is when should the blocked controller emit a communications packet of high enough priority to override the block? One answer is straightforward: just before μ(tb) exceeds μ̅(tb), the known maximum sample period that will guarantee plant stability. (In [19] we suggest that even longer delays are possible under certain conditions.) Intuitively, lower gains permit longer blocking delays.
The example of the previous paragraph is closely modeled using the variable ℙ[a,b] time-scale, which is continuous for an interval a, then has a gap for an interval b, then repeats. Figure 1 shows the regulation of a system implementing an adaptive gain controller in a blocking situation:
In graph (A), output y(t) (thick line) and gain k(t) are plotted. In graph (B), μ(t) (thick line) is kept at 90% of μ(t) during blocking. A block occurs after every 1 second of continuous control.
A=[01-11],B=[11],C=[1,0].
The gain begins at k0=1/2 and the sampling period at μ(0)=0. The bounding function for the graininess is μ̅(t)=1.9/k(t) (so that μ̅kCB<2).
7. Conclusions
In summary, the paper illustrates a new unified continuous/discrete approach to the high-gain adaptive controller. Using developments in the new field of time-scale theory, the unified results reveal that this type of feedback control works well on a much wider variety of time scales than explored in previous literature, including those that switch between continuous (or nearly continuous) and discrete domains or those without monotonically decreasing graininess. A simulation of an adaptive controller on a mixed continuous/discrete time scale is also given. It is our hope that time-scale theory may find wider application in the broad fields of signals and systems as it seems that many of the tools needed in those fields are beginning to appear in their generalized forms.
Appendix
We comment on the properties of the “expc” function referenced in the main body of the paper.
Lemma A.1.
The power series (3.3) has the following properties:
expc(X)X=Xexpc(X).
expc(X)=(eX-I)X-1 when X-1 exists.
For real, scalar arguments x, expc(x)=eix/2sinc(x/2), where sinc denotes the sine cardinal function. (This is the motivation for the expc notation.)
∥expc(X)∥≤exp(∥X∥).
∥expc(X)-I∥≤exp(∥X∥)-1.
Proof.
Parts 1–3 follow immediately from the definition. To verify 4, note
∥∑n=1∞Xn-1n!∥≤∑n=1∞∥X∥n-1n!≤∑n=1∞∥X∥n-1(n-1)!=exp(∥X∥).
Part 5 follows from a similar argument. Note that, by property 5, the decomposition expc(μX)=I+Σ(μ) gives ∥Σ(μ)∥→0 as μ→0, with uniform convergence.
Acknowledgments
This work was supported by NSF Grants no. EHS-0410685 and CMMI no. 726996 as well as a Baylor University Research Council grant. The authors thank their colleague, Robert J. Marks II, for his very helpful suggestions throughout this project.
IlchmannA.OwensD. H.Prätzel-WoltersD.High-gain robust adaptive controllers for multivariable systems19878539740410.1016/0167-6911(87)90077-6MR890075ZBL0632.93046KhalilH. K.SaberiA.Adaptive stabilization of a class of nonlinear systems using high-gain feedback198732111031103510.1109/TAC.1987.1104481MR909974ZBL0625.93040OwensD. H.Pratzel-WoltersD.IlchmannA.Positive-real structure and high-gain adaptive stabilization1987416718110.1093/imamci/4.2.167ZBL0631.93058WillemsJ. C.ByrnesC. I.Global adaptive stabilization in the absence of information on the sign of the high frequency gain198462Berlin, GermanySpringer4957Lecture Notes in Control and Information SciencesMR87668210.1007/BFb0004944ZBL0549.93043OwensD. H.Adaptive stabilization using a variable sampling rate19966311071192-s2.0-002334440410.1080/00207179608921834ZBL0845.93074IlchmannA.TownleyS.Adaptive high-gain λ-tracking with variable sampling rate1999364285293MR175202010.1016/S0167-6911(98)00101-7IlchmannA.TownleyS.Adaptive sampling control of high-gain stabilizable systems199944101961196610.1109/9.793786MR1716167ZBL0956.93056IlchmannA.RyanE. P.On gain adaptation in adaptive control200348589589910.1109/TAC.2003.811276MR1980604LogemannH.MårtenssonB.Adaptive stabilization of infinite-dimensional systems199237121869188310.1109/9.182474MR1200604ZBL0797.93041ÖzdemirN.TownleyS.Integral control by variable sampling based on steady-state data20033911351402-s2.0-0029771174ZBL1016.93044BohnerM.PetersonA.2003Boston, Mass, USABirkhäuserxii+348MR1962542BohnerM.PetersonA.2001Boston, Mass, USABirkhäuserx+358MR1843232PötzscheC.SiegmundS.WirthF.A spectral characterization of exponential stability for linear time-invariant systems on time scales2003951223124110.3934/dcds.2003.9.1223MR1974425ZBL1054.34086DaCunhaJ. J.Stability for time varying linear dynamic systems on time scales2005176238141010.1016/j.cam.2004.07.026MR2116401ZBL1064.39005DaCunhaJ. J.2004Waco, Tex, USABaylor UniversitySlotineJ. E.LiW.1991Upper Saddle River, NJ, USAPrentice HallDavisJ. M.GravagneI. A.MarksR. J.IIRamosA. A.Algebraic and dynamic Lyapunov equations on time scalessubmittedCervinA.HenrikssonD.LincolnB.EkerJ.ÅrzénK.-E.How does control timing affect performance? Analysis and simulation of timing using Jitterbug and TrueTime200323316302-s2.0-003269021510.1109/MCS.2003.1200240GravagneI. A.DavisJ. M.DaCunhaJ. J.MarksR. J.IIBandwidth reduction for controller area networks using adaptive samplingProceedings of the IEEE International Conference on Robotics and AutomationApril 2004New Orleans, La, USA525052552-s2.0-0031122310