On Necessary and Sufficient Conditions for Incremental Stability of Hybrid Systems Using the Graphical Distance Between Solutions

Size: px
Start display at page:

Download "On Necessary and Sufficient Conditions for Incremental Stability of Hybrid Systems Using the Graphical Distance Between Solutions"

Transcription

1 On Necessary and Sufficien Condiions for Incremenal Sabiliy of Hybrid Sysems Using he Graphical Disance Beween Soluions Yuchun Li and Ricardo G. Sanfelice Absrac This paper inroduces new incremenal sabiliy noions for a class of hybrid dynamical sysems given in erms of differenial equaions and difference equaions wih sae consrains. Incremenal sabiliy is defined as he propery ha he disance beween every pair of soluions o he sysem has sable behavior (incremenal sabiliy) and approaches zero asympoically (incremenal araciviy) in erms of graphical convergence. Basic properies of he class of graphically incremenally sable sysems are considered as well as hose implied by he new noions are revealed. Moreover, several sufficien and necessary condiions for a hybrid sysem wih such a propery are esablished. Examples are presened hroughou he paper o illusrae he noions and resuls. I. INTRODUCTION In recen years, incremenal sabiliy-like properies have been used in he sudy of synchronizaion [], [2], [3], observer design [4], [5], conrol design [6], as well as he sudy of convergen sysems [7]. Unforunaely, he incremenal sabiliy noions based on Euclidean and Riemannian disance canno be applied direcly o sysems wih variables ha can change coninuously and, a imes, jump discreely. These sysems, known as hybrid sysems, are capable of modeling a wide range of complex dynamical sysems, including roboic, auomoive, and power sysems as well as naural processes. Alhough se sabiliy heory in erms of Lyapunov funcions is available (see [8], [9]), incremenal sabiliy noions for such sysems would enable he sudy of properies similar o he curren noion for coninuousime sysems. While he iniial effor in [] defines an incremenal sabiliy noion ha prioriizes ordinary ime, he noion herein migh be difficul o guaranee for hybrid sysems ha conain he so called peaking phenomenon as recognized in, e.g., [], [2]. In his paper, we inroduce conceps for he sudy of incremenal graphical sabiliy of hybrid sysems ha may have he peaking phenomenon. In paricular, for a class of hybrid sysems in he framework of [8], we inroduce new noions of incremenal sabiliy and incremenal araciviy based on he graphical disance beween soluions. Moreover, sufficien condiions based on hose found in conracion heory and finie-ime sabiliy heory wihin a neighborhood of he jump se are esablished. Furhermore, i is shown ha finie-ime sabiliy wihin a neighborhood of he jump se is also a necessary condiion for a hybrid sysem o be graphically incremenally sable or aracive. To he bes Y. Li and R. G. Sanfelice are wih he Deparmen of Compuer Engineering, Universiy of California, Sana Cruz, CA yuchunli,ricardo@ucsc.edu This research has been parially suppored by he Naional Science Foundaion under CAREER Gran no. ECS- 536 and by he AFOSR under Gran no. FA of our knowledge, he noion of incremenal sabiliy and is properies for hybrid sysems have no been horoughly sudied, only discussed briefly in [3] for a class of ransiion sysems in he conex of bisimulaions, and in [] for a paricular class of hybrid sysems prioriizing ordinary ime. The remainder of his paper is organized as follows. Secion II briefly reviews he hybrid sysem framework used and is basic properies. Secion III inroduces he noion of graphical incremenal sabiliy and araciviy, and illusraes hem in wo examples. Secion IV presens he main resuls, which consis of several sufficien and necessary condiions. Examples are discussed hroughou he paper o illusrae he resuls. Due o he space limiaions, complee proofs will be published elsewhere. A. Noaion II. PRELIMINARIES ON HYBRID SYSTEMS AND GRAPH NOTIONS Given a se S R n, he closure of S is he inersecion of all closed ses conaining S, denoed by S; S is said o be discree if here exiss δ > such ha for each x S, (x+δb) S={x}; cons is he closure of he convex hull of he se S. R :=[, ) and N := {,,...}. Given vecors ν R n, w R m, ν defines he Euclidean vecor norm ν = ν ν and [ν w ] is equivalen o (ν,w). Given a funcion f : R m R n, is domain of definiion is denoed by dom f, i.e., dom f := {x R m : f(x) is defined}. The range of f is denoed by rge f, i.e., rge f := {f(x) : x domf}. The righ limi of he funcion f is defined as f + (x) := lim ν + f(x+ν) if i exiss. Given a poin y R n and a closed se A R n, y A := inf x A x y. A funcion α : R R is a class-k funcion, also wrien α K, if α is zero a zero, coninuous, sricly increasing, and unbounded; α is posiive definie, also wrien α PD, if α(s)> for all s> and α()=. Given a real number x R, floor(x) is he closes ineger o x from below. A funcion V : R n R is called a Lyapunov funcion wih respec o a se A if V is coninuously differeniable and such ha c ( x A ) V(x) c 2 ( x A ) for all x R n and some funcions c,c 2 K. B. Preliminaries In his paper, a hybrid sysem H has daa (C,f,D,g) and is defined by ż = f(z) z C, z + () = g(z) z D,

2 where z R n is he sae, he map f defines he flow map capuring he coninuous dynamics and C defines he flow se on whichf is effecive. The map g defines he jump map and models he discree behavior, while D defines he jump se, which is he se of poins from where jumps are allowed. A soluion φ o H is paramerized by (,j) R N, where denoes ordinary ime and j denoes jump ime. The domain dom φ R N is a hybrid ime domain when i saisfies [8, Definiion 2.3]. A soluion o H is called maximal if i canno be exended, i.e., i is no a runcaed version of anoher soluion. I is called complee if is domain is unbounded. A soluion is Zeno if i is complee and is domain is bounded in he direcion. A soluion is precompac if i is complee and bounded. The se S H conains all maximal soluions o H, and he se S H (ξ) conains all maximal soluions o H from ξ. A hybrid sysem H is said o saisfy he hybrid basic condiions if i saisfies [8, Definiion 6.5]. The propery of pre-forward compleeness of a hybrid sysem H is characerized in [8, Definiion 6.2]. The graph of a hybrid arc is defined in [8, Definiion 5.2], and he disance of graphs of wo hybrid arcs is measured by ε-closeness noion in [8, Definiion 4.]. In order o characerize he propery of hybrid arcs graphically converging o each oher, we inroduce he following noion. Definiion 2.: Given ε >, wo hybrid arcs φ and φ 2 are evenually ε-close if here exiss T > such ha (a) for each (,j) domφ such ha + j > T, here exiss (s,j) domφ 2 saisfying s < ε and φ (,j) φ 2 (s,j) < ε, (2) (b) for each (,j) domφ 2 such ha + j > T, here exiss (s,j) domφ saisfying s < ε and φ 2 (,j) φ (s,j) < ε. (3) We refer he reader o [8] and [9] for more deails on hese noions and he hybrid sysems framework. III. INCREMENTAL GRAPHICAL STABILITY In his paper, for hybrid sysems H as in (), we are ineresed in characerizing he incremenal sabiliy propery, namely, he noion ha he graphical disance beween every pair of maximal soluions o he sysem has sable behavior and approaches zero asympoically. To highligh he inricacies of his propery in he hybrid seing, a canonical example of hybrid sysems, he so-called bouncing ball sysem, is considered. Every maximal soluion o he bouncing ball sysem is Zeno and converges o he origin; see [8, Example. and 2.2] for more deails. Consider wo soluions o his sysem, given by φ and φ 2, from iniial condiions φ (,) = (5,) (ball iniialized a a posiive heigh wih zero velociy) and φ 2 (,)=(,3) (ball iniialized a he ground wih a posiive velociy). Figure (a) shows he posiion (firs) componen (φ i for i {,2}) of hese wo soluions, and Figure (b) shows he velociy (second) componen (φ 2 i for i {,2}) of hem. The Zeno behavior of he soluions makes i A soluion is Zeno if i is complee and is domain is bounded in he direcion j φ φ 2 (a) The firs componen (heigh φ i ) of soluions from φ (,) = (5,) and φ 2 (,) = (,3). 6 8 d( φ 2, φ 2 2 ) (b) The projecion of he second componen (velociy φ 2 i ) of soluions from φ (,) = (3,3) and φ 2 (,) = (3,3.) on he direcion. Fig.. Soluions o he bouncing ball sysem. The Euclidean disance, which is d( φ 2 (,j ), φ 2 2 (,j 2)) := φ 2 (,j ) φ 2 2 (,j 2) for all (,j ) dom φ and (,j 2 ) dom φ 2, has repeiive large peaks. exremely difficul o analyze incremenal propery of hese wo soluions. In fac, since graphical incremenal sabiliy requires comparing he graphs beween wo soluions, if one soluion (φ ) reaches he Zeno ime sooner han he second soluion (φ 2 ), hen he graphical disance beween hese wo soluions canno be evaluaed afer one of hem has approached he Zeno ime. This siuaion is shown in Figure (a), where φ 2 approaches Zeno a abou = 6sec while φ is sill describing he moion of he ball bouncing. Such exreme difference in he domains of he soluions makes i difficul (if no impossible) o compare φ and φ 2. Furhermore, a similar siuaion is encounered if, insead, he poinwise disance is used o measure he disance beween soluions. As shown in Figure (b), he poinwise disance of velociy (second) componens of wo soluions ( φ 2 i for i {,2}) has repeiive large peaks, even hough hey are iniialized very close o each oher. To avoid such possible dramaic differences on he hybrid ime domains and disances beween soluions o a hybrid sysem, hroughou he paper, we consider a class of hybrid sysems ha saisfies he following assumpion. Assumpion 3.: The hybrid sysem H = (C,f,D,g) is such ha ) each maximal soluion o H has a hybrid ime domain ha is unbounded in he direcion; 2) he flow map f is coninuously differeniable; 3) here exiss γ > such ha for each maximal soluion o H, he flow ime beween wo consecuive jumps is lower bounded by γ. Lemma 2.7 in [4] provides a sufficien condiion for condiion 3). As we show laer on, hese assumpions are no resricive, in fac, some of hem are necessary for esablishing sufficien condiions for a hybrid sysem o be graphically incremenally sable. The noion of incremenal sabiliy used in his work measures he graph disance beween soluions o hybrid sysems. I is defined as follows. Definiion 3.2: Consider a hybrid sysem H wih sae z R n. The hybrid sysem H is said o be ) incremenally graphically sable (δs) if for every ε > here exiss δ > such ha, for any wo maximal φ 2 φ 2 2

3 soluions φ,φ 2 o H, φ (,) φ 2 (,) δ implies φ and φ 2 are ε-close; 2) incremenally graphically locally aracive (δla) if here exissµ > such ha for everyε >, for any wo maximal soluionsφ,φ 2 oh, φ (,) φ 2 (,) µ implies ha φ and φ 2 are evenually ε-close; 3) incremenally graphically globally aracive (δga) if for every ε >, and for any wo maximal soluions φ,φ 2 o H, φ and φ 2 are evenually ε-close; 4) incremenally graphically locally asympoically sable (δlas) if i is boh δs and δla. Remark 3.3: Noe ha he δla noion is differen from he δs noion. The former requires ha every pair of maximal soluions o H iniialized close converge o each oher graphically when complee, while he laer requires ha every pair of wo maximal soluions o H iniialized close say close graphically. The following examples illusrae some of he incremenal sabiliy noions (δs and δlas) in Definiion 3.2. Example 3.4: (A imer) Consider he hybrid sysem H ż = z [,], (4) z + = z =. (5) Noe ha every maximal soluion o H is complee. Consider wo maximal soluions φ, φ 2 o he sysem. To show δs, for a given ε >, le < δ < ε and assume φ (,) φ 2 (,) < δ. Wihou loss of generaliy, we furher suppose φ (,) > φ 2 (,). Then, soluion φ jumps before φ 2. For each j N \{}, le j = max (,j ) domφ domφ 2 and j = min (,j) domφ domφ 2. Then, we have ha for each [, ], here exiss (s,) domφ 2 such ha s = and φ (,) φ 2 (,) = φ (,)+ φ 2 (,) δ<ε. (6) For each [, ], φ (,) φ 2 (,) = φ 2 (,) (7) φ (,) φ 2 (,) δ < ε, where we used he fac ha φ (,)+ =. Moreover, for each [, ], φ (,) φ 2 (,) = δ < ε. (8) Noe ha δ < ε holds for all [, ]. In fac, for each [ i, i ], where i N\{,}, φ (,i ) φ 2 (,i ) = (φ ( i,i ) φ 2( i,i )) δ < ε. Moreover, for each [ i, i ], where i N\{,}, φ ( i,i ) φ 2 (,i ) = φ ( i,i ) φ 2 ( i,i ) ( i ) < ε, and φ (,i) φ 2 ( i,i) = i δ < ε. Therefore, he sysem is δs. Moreover, since he disance beween φ and φ 2 does no converge o zero, φ,φ 2 are no evenually ε-close and hus he sysem is neiher δla nor δga. As shown in Figure 2(a), he domains of wo soluions o he imer sysem are differen from each oher. In such case, he Euclidean disance may no be a good candidae of a disance funcion for he sudy of incremenal properies, φ φ (a) The projecions of wo soluions from φ (,) = and φ 2 (,) =.2 on he direcion. d(φ, φ2) dg(φ, φ2) (b) Comparison beween Euclidean disance (op) and graphical disance (boom) beween φ and φ 2. Fig. 2. Two soluions φ and φ 2 o he imer sysem in Example 3.4. Unlike he Euclidean disance, which is d(φ (,j ),φ 2 (,j 2 )) := φ (,j ) φ 2 (,j 2 ) for all (,j ) domφ and (,j 2 ) domφ 2, which assumes he value.8 for.2 seconds afer every.8 seconds, he graphical disance saisfies d g(φ (,j),φ 2 ) := min{ φ (,j) φ 2 (s,j) : s.2,(,j) domφ,(s,j) domφ 2 }, which is upper bounded by.2. as shown in he op subfigure of Figure 2(b). Noe ha no maer how close he wo maximal soluions are iniialized, he peak always exiss for he Euclidean disance beween hem. However, he graphical disance beween he soluions, as shown in he boom subfigure of Figure 2(b), is bounded by.2 for all hybrid ime in he respecive domain of he soluions. We will presen a sysemaic way o check he propery laer on in his paper. 2 Example 3.5: Consider a hybrid sysem H wih daa f(z) : = z z C := [i,i+] i {2k:k N} (9) g(z) : = z z D := {2k : k N\{}}. Noe ha every maximal soluion o H is complee. Given ε >, consider wo maximal soluions φ and φ 2 such ha φ (,) φ 2 (,) < δ, where δ = min{,ε}. Then, i is guaraneed ha J := sup (,j) domφ j = sup (,j) domφ2 j <. For eachj N\{}, le j = max (,j ) domφ domφ 2 and j = min (,j) domφ domφ 2. Wihou loss of generaliy, assume φ 2 (,) > φ (,) 2, hen φ jumps firs. Then, we have ha for each [, ], here exiss (s,) domφ 2 such ha s = and φ (,) φ 2 (,) = φ (,)e φ 2 (,)e δ<ε. () For each [, ], φ (,) φ 2 (,) = e φ (,) e φ 2 (,) e φ (,) e φ 2 (,) δ < ε, () where we used he propery e φ (,)=floor(φ (,))= floor(φ 2 (,)) = e φ2 (,). Noe ha = ln(φ (,)) ln(floor(φ (,))) and = ln(φ 2 (,)) ln(floor(φ 2 (,))), herefore, = ln(φ 2 (,)) ln(φ (,)). Furhermore, by mean value heorem, here exiss φ [φ (,),φ 2 (,)] such ha = φ φ (,) φ 2 (,) φ (,) φ 2 (,) δ < ε. Similarly, for each [, ], φ (,) φ 2 (,) φ 2(,) φ (,) δ. (2) If φ (,),φ 2 (,), we have ha φ (,) 2 I may be possible o consruc an alernaive disance funcion ha is decreasing along soluions using he ideas in [5].

4 φ φ (a) The projecions of wo soluions from φ (,) = 4.5 and φ 2 (,) = 4.2 on he direcion. d(φ, φ2) dg(φ, φ2) (b) Comparison beween graphical disance and Euclidean disance beween φ and φ 2. Fig. 3. Two soluions φ and φ 2 o he sysem H. Unlike he Euclidean disance, which is d(φ (,j ),φ 2 (,j 2 )) := φ (,j ) φ 2 (,j 2 ) for all (,j ) domφ and (,j 2 ) domφ 2, which does no decrease along soluions, he graphical disance saisfies d g(φ (,j),φ 2 ) := min{ φ (,j) φ 2 (s,j) : s.2,(,j) domφ,(s,j) domφ 2 }, which decreases along soluions. φ 2 (,) e φ (,) φ 2 (,) δ for all (,) domφ = domφ 2. In fac, lim,(,) domφ φ (,) φ 2 (,) =. Therefore, he sysem H is δlas. As shown in Figure 3(a), he domains of wo soluions o he sysem H are differen from each oher. The Euclidean disance beween φ and φ 2 is shown in he op subfigure of Figure 3(b), while he graphical disance beween hem, as shown in he boom subfigure of Figure 3(b), converges o zero asympoically. IV. SUFFICIENT AND NECESSARY CONDITIONS FOR GRAPHICAL INCREMENTAL STABILITY NOTIONS In his secion, we presen several sufficien and necessary condiions of incremenal graphical sabiliy properies for hybrid sysems ha saisfy Assumpion 3.. The following resul highlighs a necessary propery of he hybrid ime domains of wo hybrid arcs ha are graphically close. In paricular, i holds for soluions o a hybrid sysem H ha is δs or δga. Lemma 4.: Given ε > and wo hybrid arcs φ and φ 2, he following hold: ) If φ and φ 2 are ε-close, hen sup j = sup j. (3) (,j) domφ (,j) domφ 2 2) If φ and φ 2 are complee and evenually ε-close, hen (3) holds. Example 4.2: Consider he Example 3.5, and wo soluions φ and φ 2 o H from φ (,)=4.5 and φ 2 (,)=, respecively. Soluionφ jumps wice while soluionφ 2 never jumps. Therefore, φ and φ 2 is no evenually ε-close which prevens he sysem from being δga even hough i is shown o be δlas in Example 3.5. The following resul esablishes ha uniqueness is a necessary condiion for δs. In urn, in paricular, i jusifies he choice of single-valued flow and jump maps in he definiion of H in (). Proposiion 4.3: (uniqueness of soluions under δs) Consider a hybrid sysem H wih sae z R n. Suppose H saisfies Assumpion 3. and H is δs. Then, every maximal soluion o H is unique. Based on Proposiion 4.3, assuming uniqueness of soluions o H is no a all resricive when sudying incremenal graphical sabiliy. Hence, in he following resuls we impose he following uniqueness of soluions assumpion. Assumpion 4.4: The hybrid sysem H = (C,f,D,g) is such ha each maximal soluion φ o H is unique. A sufficien condiion for guaraneeing uniqueness of maximal soluions requires f o be locally Lipschiz and no flow from C D, see [8, Proposiion 2.]. When he jump se D is a discree se, he following sufficien condiions for a hybrid sysem H o be δlas are esablished. In paricular, he graphical disance beween any wo maximal soluions o a hybrid sysem H sricly decreases during flows. Theorem 4.5: (δlas hrough flow wih D being a discree se) Consider a hybrid sysem H = (C,f,D,g) wih sae z R n. Suppose H saisfies Assumpion 3., Assumpion 4.4, and he hybrid basic condiions. Moreover, suppose D is a discree se. If here exis β > and δ > such ha H saisfies ) f(z)+ f (z) 2βI for all z conc; 2) for each δ [,δ ], each maximal soluion φ o H from φ(, ) saisfying φ(,) C, φ(,) D = δ (4) is such ha here exiss s [,δ] saisfying φ(s,) D =, φ(,) D δ [,s], (5) and he maximal soluion φ o H from g(φ(s,)) saisfies hen, H is δlas. φ(,) g(φ(s,)) δ [,s]; (6) Skech of he Proof: Given ε >, and using δ as in he iem 2) of assumpion and γ as in Assumpion 3., consider wo maximal soluions φ, φ 2 o H such ha φ (,) φ 2 (,) < δ, where δ is chosen such ha < δ min{ε,δ,γ} and for each z D, (z+δb) D = {z}. Ifφ (,),φ 2 (,) C and no jump occurs o eiherφ or φ 2, hen, by he generalized mean value heorem for vecorvalued funcions and iem ), for all (, ), we have d d φ (,) φ 2 (,) 2 2β φ (,) φ 2 (,) 2. (7) Then, by he comparison lemma, for all [, ), φ (,) φ 2 (,) exp( β) φ (,) φ 2 (,) δ. (8) Now consider he case when eiher φ or φ 2 jump. Assume φ jumps firs and J =. Furhermore, for each j N \ {}, le j = max (,j ) domφ domφ 2 and j = min (,j) domφ domφ 2, and =. Case I: When φ (,),φ 2 (,) C, similarly as in (7), for all [, ], we have ha φ (,) φ 2 (,) δ. (9) When =, since φ jumps firs, φ (,) D and φ (,) = g(φ (,)). Noe ha under iem 3) of Assumpion 3., g(d) D =. Therefore, by using (7) and iem 2), we obain

5 a) for each j N\{} and each [ j, j+ ]: φ (,j) φ 2 (,j) exp( β( j+ j )) φ (,) φ 2 (,) ε, where j := j k= ( k k ), b) for each j N\{} and each [ j, j ]: φ 2 (,j ) φ ( j,j ) exp( β j ) φ (,) φ 2 (,) ε, c) for each j N\{} and each [ j, j ]: (2) (2) φ (,j) φ 2 ( j,j) (22) exp( β j ) φ (,) φ 2 (,) ε. Therefore, φ and φ 2 are ε-close. Case II If φ (,),φ 2 (,) D, since D is a discree se and φ (,) φ 2 (,) δ, φ (,) = φ 2 (,). By he uniqueness of soluions, we have ha φ (,j) = φ 2 (,j) for all (,j) domφ. Thus, φ and φ 2 are ε-close. Case III If φ (,) C,φ 2 (,) D, he argumens follows similarly as in Case I. Therefore, by combining argumens in hree cases, i is proved ha φ and φ 2 are ε-close which implies H is δs. Noe ha he case when J < follows similarly wih [ J, ) {J} domφ domφ 2 and above. On he oher hand, he δla propery follows from (2), (2) and (22). Remark 4.6: The firs condiion in Theorem 4.5 guaranees sric decrease of he graphical disance beween wo maximal soluions on he inersecions of heir hybrid ime domains. The second condiion in Theorem 4.5 implies ha on he mismached pars of heir hybrid ime domains, he graphical disance beween hem does no grow. The condiion proposed in iem 2) of Theorem 4.5 can be guaraneed by he following sufficien condiion. Proposiion 4.7: Consider a hybrid sysem H = (C,f,D,g) wih sae z R n. Suppose H saisfies Assumpion 3., Assumpion 4.4, and he hybrid basic condiions. Moreover, suppose D is a discree se. Then, iem 2) of Theorem 4.5 holds if here exiss δ > such ha, for any z D, he following hold: here exis c,c 2 >, c 2 (,c ], and α (,) such ha ) for V (z) = z z 2, we have V (z),f(z) + c V α(z) < and z z 2α c ( α) for all z C ((z +δ B)\D), 2) for V 2 (z) = z g(z ) 2, we have V 2 (z),f(z) c 2 V2 α(z) < for all z C (g(z )+δ B). The following example illusraes he sufficien condiion in Theorem 4.5. Example 4.8: Consider he sysem in Example 3.5 and he se M = [,M], where M >. For he hybrid sysem H M = (C M,f,D M,g), he condiions in Theorem 4.5 can be verified as follows. Each maximal soluion o H M is complee and is domain is unbounded in he direcion. Moreover, he flow map is coninuously differeniable on con(c M). Furhermore, for any maximal soluion φ o he sysem H M from φ(,) (C D) M, denoe d := max{x : x C,x φ(,)}. If φ(,), hen φ never jumps and he jump ime beween wo consecuive jumps is bounded below by. If d 2, he flow ime beween wo consecuive jumps of φ is bounded below by ln d d ln floor(m+2) floor(m+). For all z con(c M), f(z)+ f(z) = 2, so iem ) in Theorem 4.5 is saisfied wih β =. Moreover, given z D M, he Lyapunov funcion V (z) = z z 2 saisfies V (z),f(z) = 2(z z )( z) 2z (z z ) = 2z V 2 (z) for z (M C) ((z + ln d d B) \ (M D)), where we used he propery ha z z for all z (M C) ((z + ln d d B)\(M D)). Furhermore, he Lyapunov funcion V 2 (z) = z g(z ) 2 saisfies V 2 (z),f(z) = 2(z g(z ))( z) 2z (g(z ) z) = 2g(z )V 2 2 (z) for z (g(z )+ln d d B) (C M), where we used he propery ha z g(z ) for all z (g(z ) + ln d d B) (C M) and g(z ) = z < z. Then, by Proposiion 4.7 and Theorem 4.5, we have ha H M is δlas. I follows ha he finie-ime convergence propery in iem 2) of Theorem 4.5 is a necessary condiion for a hybrid sysem H o be δs or δla. Indeed, wihou he finie-ime convergence propery nearby D and g(d), he graphs of he soluions would no be close. Theorem 4.9: (necessary condiion for δs and δla) Consider a hybrid sysem H wih sae z R n. Suppose H saisfies Assumpion 3., Assumpion 4.4, and he hybrid basic condiions. Furhermore, suppose H is δs or δla. Then, here exiss δ > such ha each maximal soluion φ o H from φ(,) saisfying φ(,) D δ andφ(,) C converges o D wihin finie ime, i.e., here exiss s > such ha φ(s,) D =. Under furher assumpions, he above resuls for a discree jump se can be exended o he case of a generic se D. For his purpose, we inroduce he following forward invarian noion. Definiion 4.: (forward invariance from away of D) A se A R n is said o be forward invarian for H from away of D if for each soluion φ o H from φ(,) A \ D, φ(,) A for all (,) domφ. Now, we are ready o presen sufficien condiions for a generic jump se. Theorem 4.: (δlas hrough flow for generic D) Consider a hybrid sysem H = (C,f,D,g) wih sae z R n. SupposeH saisfies Assumpion 3., Assumpion 4.4, and he hybrid basic condiions. If here exis β > and δ > such ha H saisfies ) f(z)+ f (z) 2βI for all z conc; 2) for each δ [,δ ], each maximal soluion φ o H saisfying φ(,) C, φ(,) D = δ (23) is such ha here exiss s [,δ] saisfying φ(s,) D = (24) and he se φ(s,)+δb is forward invarian from away

6 of D, and each maximal soluion φ o H from φ(,) g(φ(s,))+δb (25) saisfies φ(,) g(φ(s,))+δb (26) for all [,s]; 3) he jump map g is locally Lipschiz on D wih Lipschiz consan L [,], 3 i.e., g(z ) g(z 2 ) L z z 2 for all z,z 2 D such ha z z 2 δ ; hen, H is δlas. The following resul esablishes a sufficien condiion for a hybrid sysem H o be δlas hrough jumps. In paricular, he graphical disance beween any wo maximal soluions o a hybrid sysem H sricly decreases during jumps. Theorem 4.2: (δlas hrough jumps for generic D) Consider a hybrid sysem H = (C,f,D,g) wih sae z R n. SupposeH saisfies Assumpion 3., Assumpion 4.4, and he hybrid basic condiions. If here exis δ,l,l 2 > such ha ) f(z)+ f(z) for all z conc; 2) for each δ [,δ ], each maximal soluion φ o H from φ(, ) saisfying φ(,) C, φ(,) D = δ saisfies φ(s,) D = for some s [,δ]; 3) for each z D and each δ [,δ ], he se z +δb is forward invarian for H from away of D; 4) he jump map g is locally Lipschiz on D wih Lipschiz consan L, i.e., g(z ) g(z 2 ) L z z 2 for all z,z 2 D such ha z z 2 δ ; 5) f is bounded on conc wih bound L 2, i.e., f(z) L 2 for all z conc; 6) L +L 2 ; hen, H is δs. Furhermore, if he domain of each maximal soluion o H is unbounded in he j direcion, and L and L 2 can be chosen such ha L +L 2 <, hen, H is δlas. The following example illusraes he condiions in Theorem 4.2. Example 4.3: Consider he imer sysem in Example 3.4. Each maximal soluionφ o i has a domain ha is unbounded in he and j direcion. Moreover, he flow ime beween wo consecuive jumps of φ is lower bounded by. The condiion in iem ) of Theorem 4.2 can be verified as f(z)+ f(z) = for all z conc. The condiion in iem 2) can be verified according o Proposiion 4.7. Consider δ (,) and he Lyapunov funcion V(z) = z 2 D. For each z (D+δ B) C\D, we have V(z)=(z ) 2 and V(z),f(z) = 2( z)= 2V 2(z), where we used he propery ha z for all z (D+δ B)\D. Iem 3) of Theorem 4.2 follows from he fac D={} is a singleon and V(z),f(z) = 2( z)< for allz (D+δ B) C\D. Iem 4) of Theorem 4.2 is saisfied wih c =, and iem 5) of Theorem 4.2 is saisfied wih c 2 =. Therefore, he imer sysem in Example 3.4 is δs. 3 Such g is also known as a weak conracion map. V. CONCLUSION In his paper, we show ha graphical incremenal sabiliy is a key noion for he sudy of incremenal sabiliy for hybrid sysems. Oher noions based on poinwise Euclidean disance fall shor when applied o sysems ha exhibi he peaking phenomenon, which is a ypical behavior in racking and observer design for hybrid sysems. Several sufficien and necessary condiions for a hybrid sysem o be graphically incremenally sable and graphically incremenally aracive were provided and illusraed in examples. REFERENCES [] C. Cai and G. Chen. Synchronizaion of complex dynamical neworks by he incremenal ISS approach. Physica A: Saisical Mechanics and is Applicaions, 37(2): , 26. [2] Q. Pham, N. Tabareau, and J. J. Sloine. A conracion heory approach o sochasic incremenal sabiliy. IEEE Transacions on Auomaic Conrol, 54(4):86 82, April 29. [3] A Hamadeh, G. B San, and J. Goncalves. Consrucive synchronizaion of neworked feedback sysems. In Proc. 49h IEEE Conference on Decision and Conrol, pages , Dec 2. [4] R. G. Sanfelice and L. Praly. Convergence of nonlinear observers on R n wih a Riemannian meric (Par I). IEEE Transacions on Auomaic Conrol, 57(7):79 722, July 22. [5] A. P. Dani, S. J. Chung, and S. Huchinson. Observer design for sochasic nonlinear sysems using conracion analysis. In Proc. 5s Annual Conference on Decision and Conrol, pages , Dec 22. [6] M. Zamani and P. Tabuada. Backsepping design for incremenal sabiliy. IEEE Transacions on Auomaic Conrol, 56(9): , Sep 2. [7] A. Pavlov, N. van de Wouw, and H. Nijmeijer. Frequency response funcions for nonlinear convergen sysems. IEEE Transacions on Auomaic Conrol, 52(6):59 65, June 27. [8] R. Goebel, R. G. Sanfelice, and A. R. Teel. Hybrid Dynamical Sysems: Modeling, Sabiliy, and Robusness. Princeon Universiy Press, New Jersey, 22. [9] R. Goebel, R. G. Sanfelice, and A. R. Teel. Hybrid dynamical sysems. IEEE Conrol Sysems Magazine, 29(2):28 93, April 29. [] Y. Li, S. Phillips, and R. G. Sanfelice. Resuls on incremenal sabiliy for a class of hybrid sysems. In Proc. of IEEE 53rd Annual Conference on Decision and Conrol, pages , Dec 24. [] R. G. Sanfelice, J. J. B. Biemond, N. van de Wouw, and W. P. Maurice. H. Heemels. An embedding approach for he design of saefeedback racking conrollers for references wih jumps. In. J. Robus Nonlinear Conrol, 24:585 68, 24. [2] S. Galeani, L. Menini, and A. Poini. Robus rajecory racking for a class of hybrid sysems: An inernal model principle approach. IEEE Transacions on Auomaic Conrol, 57(2): , Feb 22. [3] P. Prabhakar, J. Liu, and R. M. Murray. Pre-orders for reasoning abou sabiliy properies wih respec o inpu of hybrid sysems. In Proc. of he Inernaional Conference on Embedded Sofware (EMSOFT), pages, Sep 23. [4] R. G. Sanfelice, R. Goebel, and A. R. Teel. Invariance principles for hybrid sysems wih connecions o deecabiliy and asympoic sabiliy. IEEE Transacions on Auomaic Conrol, 52(2): , 27. [5] J. J. B. Biemond, N. van de Wouw, W. P. M. H. Heemels, and H. Nijmeijer. Tracking conrol for hybrid sysems wih sae-riggered jumps. IEEE Transacions on Auomaic Conrol, 58(4):876 89, April 23.

The Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales

The Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales Advances in Dynamical Sysems and Applicaions. ISSN 0973-5321 Volume 1 Number 1 (2006, pp. 103 112 c Research India Publicaions hp://www.ripublicaion.com/adsa.hm The Asympoic Behavior of Nonoscillaory Soluions

More information

2. Nonlinear Conservation Law Equations

2. Nonlinear Conservation Law Equations . Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear

More information

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions Muli-Period Sochasic Models: Opimali of (s, S) Polic for -Convex Objecive Funcions Consider a seing similar o he N-sage newsvendor problem excep ha now here is a fixed re-ordering cos (> 0) for each (re-)order.

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,

More information

Optimality Conditions for Unconstrained Problems

Optimality Conditions for Unconstrained Problems 62 CHAPTER 6 Opimaliy Condiions for Unconsrained Problems 1 Unconsrained Opimizaion 11 Exisence Consider he problem of minimizing he funcion f : R n R where f is coninuous on all of R n : P min f(x) x

More information

STABILITY OF NONLINEAR NEUTRAL DELAY DIFFERENTIAL EQUATIONS WITH VARIABLE DELAYS

STABILITY OF NONLINEAR NEUTRAL DELAY DIFFERENTIAL EQUATIONS WITH VARIABLE DELAYS Elecronic Journal of Differenial Equaions, Vol. 217 217, No. 118, pp. 1 14. ISSN: 172-6691. URL: hp://ejde.mah.xsae.edu or hp://ejde.mah.un.edu STABILITY OF NONLINEAR NEUTRAL DELAY DIFFERENTIAL EQUATIONS

More information

15. Vector Valued Functions

15. Vector Valued Functions 1. Vecor Valued Funcions Up o his poin, we have presened vecors wih consan componens, for example, 1, and,,4. However, we can allow he componens of a vecor o be funcions of a common variable. For example,

More information

Hybrid Control and Switched Systems. Lecture #3 What can go wrong? Trajectories of hybrid systems

Hybrid Control and Switched Systems. Lecture #3 What can go wrong? Trajectories of hybrid systems Hybrid Conrol and Swiched Sysems Lecure #3 Wha can go wrong? Trajecories of hybrid sysems João P. Hespanha Universiy of California a Sana Barbara Summary 1. Trajecories of hybrid sysems: Soluion o a hybrid

More information

Lecture 20: Riccati Equations and Least Squares Feedback Control

Lecture 20: Riccati Equations and Least Squares Feedback Control 34-5 LINEAR SYSTEMS Lecure : Riccai Equaions and Leas Squares Feedback Conrol 5.6.4 Sae Feedback via Riccai Equaions A recursive approach in generaing he marix-valued funcion W ( ) equaion for i for he

More information

(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4)

(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4) Phase Plane Analysis of Linear Sysems Adaped from Applied Nonlinear Conrol by Sloine and Li The general form of a linear second-order sysem is a c b d From and b bc d a Differeniaion of and hen subsiuion

More information

A Necessary and Sufficient Condition for the Solutions of a Functional Differential Equation to Be Oscillatory or Tend to Zero

A Necessary and Sufficient Condition for the Solutions of a Functional Differential Equation to Be Oscillatory or Tend to Zero JOURNAL OF MAEMAICAL ANALYSIS AND APPLICAIONS 24, 7887 1997 ARICLE NO. AY965143 A Necessary and Sufficien Condiion for he Soluions of a Funcional Differenial Equaion o Be Oscillaory or end o Zero Piambar

More information

4 Sequences of measurable functions

4 Sequences of measurable functions 4 Sequences of measurable funcions 1. Le (Ω, A, µ) be a measure space (complee, afer a possible applicaion of he compleion heorem). In his chaper we invesigae relaions beween various (nonequivalen) convergences

More information

Boundedness and Exponential Asymptotic Stability in Dynamical Systems with Applications to Nonlinear Differential Equations with Unbounded Terms

Boundedness and Exponential Asymptotic Stability in Dynamical Systems with Applications to Nonlinear Differential Equations with Unbounded Terms Advances in Dynamical Sysems and Applicaions. ISSN 0973-531 Volume Number 1 007, pp. 107 11 Research India Publicaions hp://www.ripublicaion.com/adsa.hm Boundedness and Exponenial Asympoic Sabiliy in Dynamical

More information

Asymptotic instability of nonlinear differential equations

Asymptotic instability of nonlinear differential equations Elecronic Journal of Differenial Equaions, Vol. 1997(1997), No. 16, pp. 1 7. ISSN: 172-6691. URL: hp://ejde.mah.sw.edu or hp://ejde.mah.un.edu fp (login: fp) 147.26.13.11 or 129.12.3.113 Asympoic insabiliy

More information

Global Synchronization of Directed Networks with Fast Switching Topologies

Global Synchronization of Directed Networks with Fast Switching Topologies Commun. Theor. Phys. (Beijing, China) 52 (2009) pp. 1019 1924 c Chinese Physical Sociey and IOP Publishing Ld Vol. 52, No. 6, December 15, 2009 Global Synchronizaion of Direced Neworks wih Fas Swiching

More information

STATE-SPACE MODELLING. A mass balance across the tank gives:

STATE-SPACE MODELLING. A mass balance across the tank gives: B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing

More information

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still. Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in

More information

Ann. Funct. Anal. 2 (2011), no. 2, A nnals of F unctional A nalysis ISSN: (electronic) URL:

Ann. Funct. Anal. 2 (2011), no. 2, A nnals of F unctional A nalysis ISSN: (electronic) URL: Ann. Func. Anal. 2 2011, no. 2, 34 41 A nnals of F uncional A nalysis ISSN: 2008-8752 elecronic URL: www.emis.de/journals/afa/ CLASSIFICAION OF POSIIVE SOLUIONS OF NONLINEAR SYSEMS OF VOLERRA INEGRAL EQUAIONS

More information

10. State Space Methods

10. State Space Methods . Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he

More information

Stability and Bifurcation in a Neural Network Model with Two Delays

Stability and Bifurcation in a Neural Network Model with Two Delays Inernaional Mahemaical Forum, Vol. 6, 11, no. 35, 175-1731 Sabiliy and Bifurcaion in a Neural Nework Model wih Two Delays GuangPing Hu and XiaoLing Li School of Mahemaics and Physics, Nanjing Universiy

More information

4. Advanced Stability Theory

4. Advanced Stability Theory Applied Nonlinear Conrol Nguyen an ien - 4 4 Advanced Sabiliy heory he objecive of his chaper is o presen sabiliy analysis for non-auonomous sysems 41 Conceps of Sabiliy for Non-Auonomous Sysems Equilibrium

More information

Vehicle Arrival Models : Headway

Vehicle Arrival Models : Headway Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where

More information

di Bernardo, M. (1995). A purely adaptive controller to synchronize and control chaotic systems.

di Bernardo, M. (1995). A purely adaptive controller to synchronize and control chaotic systems. di ernardo, M. (995). A purely adapive conroller o synchronize and conrol chaoic sysems. hps://doi.org/.6/375-96(96)8-x Early version, also known as pre-prin Link o published version (if available):.6/375-96(96)8-x

More information

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x WEEK-3 Reciaion PHYS 131 Ch. 3: FOC 1, 3, 4, 6, 14. Problems 9, 37, 41 & 71 and Ch. 4: FOC 1, 3, 5, 8. Problems 3, 5 & 16. Feb 8, 018 Ch. 3: FOC 1, 3, 4, 6, 14. 1. (a) The horizonal componen of he projecile

More information

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation: M ah 5 7 Fall 9 L ecure O c. 4, 9 ) Hamilon- J acobi Equaion: Weak S oluion We coninue he sudy of he Hamilon-Jacobi equaion: We have shown ha u + H D u) = R n, ) ; u = g R n { = }. ). In general we canno

More information

Sliding Mode Extremum Seeking Control for Linear Quadratic Dynamic Game

Sliding Mode Extremum Seeking Control for Linear Quadratic Dynamic Game Sliding Mode Exremum Seeking Conrol for Linear Quadraic Dynamic Game Yaodong Pan and Ümi Özgüner ITS Research Group, AIST Tsukuba Eas Namiki --, Tsukuba-shi,Ibaraki-ken 5-856, Japan e-mail: pan.yaodong@ais.go.jp

More information

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details! MAT 257, Handou 6: Ocober 7-2, 20. I. Assignmen. Finish reading Chaper 2 of Spiva, rereading earlier secions as necessary. handou and fill in some missing deails! II. Higher derivaives. Also, read his

More information

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t M ah 5 2 7 Fall 2 0 0 9 L ecure 1 0 O c. 7, 2 0 0 9 Hamilon- J acobi Equaion: Explici Formulas In his lecure we ry o apply he mehod of characerisics o he Hamilon-Jacobi equaion: u + H D u, x = 0 in R n

More information

arxiv: v1 [math.fa] 9 Dec 2018

arxiv: v1 [math.fa] 9 Dec 2018 AN INVERSE FUNCTION THEOREM CONVERSE arxiv:1812.03561v1 [mah.fa] 9 Dec 2018 JIMMIE LAWSON Absrac. We esablish he following converse of he well-known inverse funcion heorem. Le g : U V and f : V U be inverse

More information

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

Course Notes for EE227C (Spring 2018): Convex Optimization and Approximation

Course Notes for EE227C (Spring 2018): Convex Optimization and Approximation Course Noes for EE7C Spring 018: Convex Opimizaion and Approximaion Insrucor: Moriz Hard Email: hard+ee7c@berkeley.edu Graduae Insrucor: Max Simchowiz Email: msimchow+ee7c@berkeley.edu Ocober 15, 018 3

More information

Convergence of the Neumann series in higher norms

Convergence of the Neumann series in higher norms Convergence of he Neumann series in higher norms Charles L. Epsein Deparmen of Mahemaics, Universiy of Pennsylvania Version 1.0 Augus 1, 003 Absrac Naural condiions on an operaor A are given so ha he Neumann

More information

Oscillation of an Euler Cauchy Dynamic Equation S. Huff, G. Olumolode, N. Pennington, and A. Peterson

Oscillation of an Euler Cauchy Dynamic Equation S. Huff, G. Olumolode, N. Pennington, and A. Peterson PROCEEDINGS OF THE FOURTH INTERNATIONAL CONFERENCE ON DYNAMICAL SYSTEMS AND DIFFERENTIAL EQUATIONS May 4 7, 00, Wilmingon, NC, USA pp 0 Oscillaion of an Euler Cauchy Dynamic Equaion S Huff, G Olumolode,

More information

POSITIVE SOLUTIONS OF NEUTRAL DELAY DIFFERENTIAL EQUATION

POSITIVE SOLUTIONS OF NEUTRAL DELAY DIFFERENTIAL EQUATION Novi Sad J. Mah. Vol. 32, No. 2, 2002, 95-108 95 POSITIVE SOLUTIONS OF NEUTRAL DELAY DIFFERENTIAL EQUATION Hajnalka Péics 1, János Karsai 2 Absrac. We consider he scalar nonauonomous neural delay differenial

More information

The motions of the celt on a horizontal plane with viscous friction

The motions of the celt on a horizontal plane with viscous friction The h Join Inernaional Conference on Mulibody Sysem Dynamics June 8, 18, Lisboa, Porugal The moions of he cel on a horizonal plane wih viscous fricion Maria A. Munisyna 1 1 Moscow Insiue of Physics and

More information

EE 315 Notes. Gürdal Arslan CLASS 1. (Sections ) What is a signal?

EE 315 Notes. Gürdal Arslan CLASS 1. (Sections ) What is a signal? EE 35 Noes Gürdal Arslan CLASS (Secions.-.2) Wha is a signal? In his class, a signal is some funcion of ime and i represens how some physical quaniy changes over some window of ime. Examples: velociy of

More information

Variational Iteration Method for Solving System of Fractional Order Ordinary Differential Equations

Variational Iteration Method for Solving System of Fractional Order Ordinary Differential Equations IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 1, Issue 6 Ver. II (Nov - Dec. 214), PP 48-54 Variaional Ieraion Mehod for Solving Sysem of Fracional Order Ordinary Differenial

More information

5. Stochastic processes (1)

5. Stochastic processes (1) Lec05.pp S-38.45 - Inroducion o Teleraffic Theory Spring 2005 Conens Basic conceps Poisson process 2 Sochasic processes () Consider some quaniy in a eleraffic (or any) sysem I ypically evolves in ime randomly

More information

Mean-square Stability Control for Networked Systems with Stochastic Time Delay

Mean-square Stability Control for Networked Systems with Stochastic Time Delay JOURNAL OF SIMULAION VOL. 5 NO. May 7 Mean-square Sabiliy Conrol for Newored Sysems wih Sochasic ime Delay YAO Hejun YUAN Fushun School of Mahemaics and Saisics Anyang Normal Universiy Anyang Henan. 455

More information

MODULE 3 FUNCTION OF A RANDOM VARIABLE AND ITS DISTRIBUTION LECTURES PROBABILITY DISTRIBUTION OF A FUNCTION OF A RANDOM VARIABLE

MODULE 3 FUNCTION OF A RANDOM VARIABLE AND ITS DISTRIBUTION LECTURES PROBABILITY DISTRIBUTION OF A FUNCTION OF A RANDOM VARIABLE Topics MODULE 3 FUNCTION OF A RANDOM VARIABLE AND ITS DISTRIBUTION LECTURES 2-6 3. FUNCTION OF A RANDOM VARIABLE 3.2 PROBABILITY DISTRIBUTION OF A FUNCTION OF A RANDOM VARIABLE 3.3 EXPECTATION AND MOMENTS

More information

International Journal of Scientific & Engineering Research, Volume 4, Issue 10, October ISSN

International Journal of Scientific & Engineering Research, Volume 4, Issue 10, October ISSN Inernaional Journal of Scienific & Engineering Research, Volume 4, Issue 10, Ocober-2013 900 FUZZY MEAN RESIDUAL LIFE ORDERING OF FUZZY RANDOM VARIABLES J. EARNEST LAZARUS PIRIYAKUMAR 1, A. YAMUNA 2 1.

More information

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux Gues Lecures for Dr. MacFarlane s EE3350 Par Deux Michael Plane Mon., 08-30-2010 Wrie name in corner. Poin ou his is a review, so I will go faser. Remind hem o go lisen o online lecure abou geing an A

More information

On Oscillation of a Generalized Logistic Equation with Several Delays

On Oscillation of a Generalized Logistic Equation with Several Delays Journal of Mahemaical Analysis and Applicaions 253, 389 45 (21) doi:1.16/jmaa.2.714, available online a hp://www.idealibrary.com on On Oscillaion of a Generalized Logisic Equaion wih Several Delays Leonid

More information

SOLUTIONS TO ECE 3084

SOLUTIONS TO ECE 3084 SOLUTIONS TO ECE 384 PROBLEM 2.. For each sysem below, specify wheher or no i is: (i) memoryless; (ii) causal; (iii) inverible; (iv) linear; (v) ime invarian; Explain your reasoning. If he propery is no

More information

Heat kernel and Harnack inequality on Riemannian manifolds

Heat kernel and Harnack inequality on Riemannian manifolds Hea kernel and Harnack inequaliy on Riemannian manifolds Alexander Grigor yan UHK 11/02/2014 onens 1 Laplace operaor and hea kernel 1 2 Uniform Faber-Krahn inequaliy 3 3 Gaussian upper bounds 4 4 ean-value

More information

On-line Adaptive Optimal Timing Control of Switched Systems

On-line Adaptive Optimal Timing Control of Switched Systems On-line Adapive Opimal Timing Conrol of Swiched Sysems X.C. Ding, Y. Wardi and M. Egersed Absrac In his paper we consider he problem of opimizing over he swiching imes for a muli-modal dynamic sysem when

More information

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes Some common engineering funcions 2.7 Inroducion This secion provides a caalogue of some common funcions ofen used in Science and Engineering. These include polynomials, raional funcions, he modulus funcion

More information

6.2 Transforms of Derivatives and Integrals.

6.2 Transforms of Derivatives and Integrals. SEC. 6.2 Transforms of Derivaives and Inegrals. ODEs 2 3 33 39 23. Change of scale. If l( f ()) F(s) and c is any 33 45 APPLICATION OF s-shifting posiive consan, show ha l( f (c)) F(s>c)>c (Hin: In Probs.

More information

18 Biological models with discrete time

18 Biological models with discrete time 8 Biological models wih discree ime The mos imporan applicaions, however, may be pedagogical. The elegan body of mahemaical heory peraining o linear sysems (Fourier analysis, orhogonal funcions, and so

More information

11!Hí MATHEMATICS : ERDŐS AND ULAM PROC. N. A. S. of decomposiion, properly speaking) conradics he possibiliy of defining a counably addiive real-valu

11!Hí MATHEMATICS : ERDŐS AND ULAM PROC. N. A. S. of decomposiion, properly speaking) conradics he possibiliy of defining a counably addiive real-valu ON EQUATIONS WITH SETS AS UNKNOWNS BY PAUL ERDŐS AND S. ULAM DEPARTMENT OF MATHEMATICS, UNIVERSITY OF COLORADO, BOULDER Communicaed May 27, 1968 We shall presen here a number of resuls in se heory concerning

More information

An Introduction to Malliavin calculus and its applications

An Introduction to Malliavin calculus and its applications An Inroducion o Malliavin calculus and is applicaions Lecure 5: Smoohness of he densiy and Hörmander s heorem David Nualar Deparmen of Mahemaics Kansas Universiy Universiy of Wyoming Summer School 214

More information

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems. Mah 2250-004 Week 4 April 6-20 secions 7.-7.3 firs order sysems of linear differenial equaions; 7.4 mass-spring sysems. Mon Apr 6 7.-7.2 Sysems of differenial equaions (7.), and he vecor Calculus we need

More information

1 Solutions to selected problems

1 Solutions to selected problems 1 Soluions o seleced problems 1. Le A B R n. Show ha in A in B bu in general bd A bd B. Soluion. Le x in A. Then here is ɛ > 0 such ha B ɛ (x) A B. This shows x in B. If A = [0, 1] and B = [0, 2], hen

More information

A Finite-time Convergent Observer with Robustness to Piecewise-constant Measurement Noise

A Finite-time Convergent Observer with Robustness to Piecewise-constant Measurement Noise A Finie-ime Convergen Observer wih Robusness o Piecewise-consan Measuremen Noise Yuchun Li a, Ricardo G. Sanfelice a a Deparmen of Compuer Engineering, Universiy of California, 56 High Sree, Sana Cruz,

More information

6. Stochastic calculus with jump processes

6. Stochastic calculus with jump processes A) Trading sraegies (1/3) Marke wih d asses S = (S 1,, S d ) A rading sraegy can be modelled wih a vecor φ describing he quaniies invesed in each asse a each insan : φ = (φ 1,, φ d ) The value a of a porfolio

More information

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should Cambridge Universiy Press 978--36-60033-7 Cambridge Inernaional AS and A Level Mahemaics: Mechanics Coursebook Excerp More Informaion Chaper The moion of projeciles In his chaper he model of free moion

More information

Solutions from Chapter 9.1 and 9.2

Solutions from Chapter 9.1 and 9.2 Soluions from Chaper 9 and 92 Secion 9 Problem # This basically boils down o an exercise in he chain rule from calculus We are looking for soluions of he form: u( x) = f( k x c) where k x R 3 and k is

More information

EECE251. Circuit Analysis I. Set 4: Capacitors, Inductors, and First-Order Linear Circuits

EECE251. Circuit Analysis I. Set 4: Capacitors, Inductors, and First-Order Linear Circuits EEE25 ircui Analysis I Se 4: apaciors, Inducors, and Firs-Order inear ircuis Shahriar Mirabbasi Deparmen of Elecrical and ompuer Engineering Universiy of Briish olumbia shahriar@ece.ubc.ca Overview Passive

More information

Robust estimation based on the first- and third-moment restrictions of the power transformation model

Robust estimation based on the first- and third-moment restrictions of the power transformation model h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,

More information

A DELAY-DEPENDENT STABILITY CRITERIA FOR T-S FUZZY SYSTEM WITH TIME-DELAYS

A DELAY-DEPENDENT STABILITY CRITERIA FOR T-S FUZZY SYSTEM WITH TIME-DELAYS A DELAY-DEPENDENT STABILITY CRITERIA FOR T-S FUZZY SYSTEM WITH TIME-DELAYS Xinping Guan ;1 Fenglei Li Cailian Chen Insiue of Elecrical Engineering, Yanshan Universiy, Qinhuangdao, 066004, China. Deparmen

More information

Signal and System (Chapter 3. Continuous-Time Systems)

Signal and System (Chapter 3. Continuous-Time Systems) Signal and Sysem (Chaper 3. Coninuous-Time Sysems) Prof. Kwang-Chun Ho kwangho@hansung.ac.kr Tel: 0-760-453 Fax:0-760-4435 1 Dep. Elecronics and Informaion Eng. 1 Nodes, Branches, Loops A nework wih b

More information

Class Meeting # 10: Introduction to the Wave Equation

Class Meeting # 10: Introduction to the Wave Equation MATH 8.5 COURSE NOTES - CLASS MEETING # 0 8.5 Inroducion o PDEs, Fall 0 Professor: Jared Speck Class Meeing # 0: Inroducion o he Wave Equaion. Wha is he wave equaion? The sandard wave equaion for a funcion

More information

MATH 5720: Gradient Methods Hung Phan, UMass Lowell October 4, 2018

MATH 5720: Gradient Methods Hung Phan, UMass Lowell October 4, 2018 MATH 5720: Gradien Mehods Hung Phan, UMass Lowell Ocober 4, 208 Descen Direcion Mehods Consider he problem min { f(x) x R n}. The general descen direcions mehod is x k+ = x k + k d k where x k is he curren

More information

BOUNDEDNESS OF MAXIMAL FUNCTIONS ON NON-DOUBLING MANIFOLDS WITH ENDS

BOUNDEDNESS OF MAXIMAL FUNCTIONS ON NON-DOUBLING MANIFOLDS WITH ENDS BOUNDEDNESS OF MAXIMAL FUNCTIONS ON NON-DOUBLING MANIFOLDS WITH ENDS XUAN THINH DUONG, JI LI, AND ADAM SIKORA Absrac Le M be a manifold wih ends consruced in [2] and be he Laplace-Belrami operaor on M

More information

Logic in computer science

Logic in computer science Logic in compuer science Logic plays an imporan role in compuer science Logic is ofen called he calculus of compuer science Logic plays a similar role in compuer science o ha played by calculus in he physical

More information

EXERCISES FOR SECTION 1.5

EXERCISES FOR SECTION 1.5 1.5 Exisence and Uniqueness of Soluions 43 20. 1 v c 21. 1 v c 1 2 4 6 8 10 1 2 2 4 6 8 10 Graph of approximae soluion obained using Euler s mehod wih = 0.1. Graph of approximae soluion obained using Euler

More information

Empirical Process Theory

Empirical Process Theory Empirical Process heory 4.384 ime Series Analysis, Fall 27 Reciaion by Paul Schrimpf Supplemenary o lecures given by Anna Mikusheva Ocober 7, 28 Reciaion 7 Empirical Process heory Le x be a real-valued

More information

Predator - Prey Model Trajectories and the nonlinear conservation law

Predator - Prey Model Trajectories and the nonlinear conservation law Predaor - Prey Model Trajecories and he nonlinear conservaion law James K. Peerson Deparmen of Biological Sciences and Deparmen of Mahemaical Sciences Clemson Universiy Ocober 28, 213 Ouline Drawing Trajecories

More information

Correspondence should be addressed to Nguyen Buong,

Correspondence should be addressed to Nguyen Buong, Hindawi Publishing Corporaion Fixed Poin Theory and Applicaions Volume 011, Aricle ID 76859, 10 pages doi:101155/011/76859 Research Aricle An Implici Ieraion Mehod for Variaional Inequaliies over he Se

More information

Linear Response Theory: The connection between QFT and experiments

Linear Response Theory: The connection between QFT and experiments Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and

More information

Electrical and current self-induction

Electrical and current self-induction Elecrical and curren self-inducion F. F. Mende hp://fmnauka.narod.ru/works.hml mende_fedor@mail.ru Absrac The aricle considers he self-inducance of reacive elemens. Elecrical self-inducion To he laws of

More information

THE MYSTERY OF STOCHASTIC MECHANICS. Edward Nelson Department of Mathematics Princeton University

THE MYSTERY OF STOCHASTIC MECHANICS. Edward Nelson Department of Mathematics Princeton University THE MYSTERY OF STOCHASTIC MECHANICS Edward Nelson Deparmen of Mahemaics Princeon Universiy 1 Classical Hamilon-Jacobi heory N paricles of various masses on a Euclidean space. Incorporae he masses in he

More information

t 2 B F x,t n dsdt t u x,t dxdt

t 2 B F x,t n dsdt t u x,t dxdt Evoluion Equaions For 0, fixed, le U U0, where U denoes a bounded open se in R n.suppose ha U is filled wih a maerial in which a conaminan is being ranspored by various means including diffusion and convecion.

More information

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals

More information

IMPLICIT AND INVERSE FUNCTION THEOREMS PAUL SCHRIMPF 1 OCTOBER 25, 2013

IMPLICIT AND INVERSE FUNCTION THEOREMS PAUL SCHRIMPF 1 OCTOBER 25, 2013 IMPLICI AND INVERSE FUNCION HEOREMS PAUL SCHRIMPF 1 OCOBER 25, 213 UNIVERSIY OF BRIISH COLUMBIA ECONOMICS 526 We have exensively sudied how o solve sysems of linear equaions. We know how o check wheher

More information

Chapter 3 Boundary Value Problem

Chapter 3 Boundary Value Problem Chaper 3 Boundary Value Problem A boundary value problem (BVP) is a problem, ypically an ODE or a PDE, which has values assigned on he physical boundary of he domain in which he problem is specified. Le

More information

KINEMATICS IN ONE DIMENSION

KINEMATICS IN ONE DIMENSION KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec

More information

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1.

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1. Roboics I April 11, 017 Exercise 1 he kinemaics of a 3R spaial robo is specified by he Denavi-Harenberg parameers in ab 1 i α i d i a i θ i 1 π/ L 1 0 1 0 0 L 3 0 0 L 3 3 able 1: able of DH parameers of

More information

) were both constant and we brought them from under the integral.

) were both constant and we brought them from under the integral. YIELD-PER-RECRUIT (coninued The yield-per-recrui model applies o a cohor, bu we saw in he Age Disribuions lecure ha he properies of a cohor do no apply in general o a collecion of cohors, which is wha

More information

!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)

!!#$%&#'()!#&'(*%)+,&',-)./0)1-*23) "#"$%&#'()"#&'(*%)+,&',-)./)1-*) #$%&'()*+,&',-.%,/)*+,-&1*#$)()5*6$+$%*,7&*-'-&1*(,-&*6&,7.$%$+*&%'(*8$&',-,%'-&1*(,-&*6&,79*(&,%: ;..,*&1$&$.$%&'()*1$$.,'&',-9*(&,%)?%*,('&5

More information

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution Physics 7b: Saisical Mechanics Fokker-Planck Equaion The Langevin equaion approach o he evoluion of he velociy disribuion for he Brownian paricle migh leave you uncomforable. A more formal reamen of his

More information

An Introduction to Backward Stochastic Differential Equations (BSDEs) PIMS Summer School 2016 in Mathematical Finance.

An Introduction to Backward Stochastic Differential Equations (BSDEs) PIMS Summer School 2016 in Mathematical Finance. 1 An Inroducion o Backward Sochasic Differenial Equaions (BSDEs) PIMS Summer School 2016 in Mahemaical Finance June 25, 2016 Chrisoph Frei cfrei@ualbera.ca This inroducion is based on Touzi [14], Bouchard

More information

Chapter 7 Response of First-order RL and RC Circuits

Chapter 7 Response of First-order RL and RC Circuits Chaper 7 Response of Firs-order RL and RC Circuis 7.- The Naural Response of RL and RC Circuis 7.3 The Sep Response of RL and RC Circuis 7.4 A General Soluion for Sep and Naural Responses 7.5 Sequenial

More information

Ordinary Differential Equations

Ordinary Differential Equations Ordinary Differenial Equaions 5. Examples of linear differenial equaions and heir applicaions We consider some examples of sysems of linear differenial equaions wih consan coefficiens y = a y +... + a

More information

Product of Fuzzy Metric Spaces and Fixed Point Theorems

Product of Fuzzy Metric Spaces and Fixed Point Theorems In. J. Conemp. Mah. Sciences, Vol. 3, 2008, no. 15, 703-712 Produc of Fuzzy Meric Spaces and Fixed Poin Theorems Mohd. Rafi Segi Rahma School of Applied Mahemaics The Universiy of Noingham Malaysia Campus

More information

Bifurcation Analysis of a Stage-Structured Prey-Predator System with Discrete and Continuous Delays

Bifurcation Analysis of a Stage-Structured Prey-Predator System with Discrete and Continuous Delays Applied Mahemaics 4 59-64 hp://dx.doi.org/.46/am..4744 Published Online July (hp://www.scirp.org/ournal/am) Bifurcaion Analysis of a Sage-Srucured Prey-Predaor Sysem wih Discree and Coninuous Delays Shunyi

More information

STABILITY OF RESET SWITCHING SYSTEMS

STABILITY OF RESET SWITCHING SYSTEMS STABILITY OF RESET SWITCHING SYSTEMS J.P. Paxman,G.Vinnicombe Cambridge Universiy Engineering Deparmen, UK, jpp7@cam.ac.uk Keywords: Sabiliy, swiching sysems, bumpless ransfer, conroller realizaion. Absrac

More information

arxiv: v1 [math.pr] 19 Feb 2011

arxiv: v1 [math.pr] 19 Feb 2011 A NOTE ON FELLER SEMIGROUPS AND RESOLVENTS VADIM KOSTRYKIN, JÜRGEN POTTHOFF, AND ROBERT SCHRADER ABSTRACT. Various equivalen condiions for a semigroup or a resolven generaed by a Markov process o be of

More information

Families with no matchings of size s

Families with no matchings of size s Families wih no machings of size s Peer Franl Andrey Kupavsii Absrac Le 2, s 2 be posiive inegers. Le be an n-elemen se, n s. Subses of 2 are called families. If F ( ), hen i is called - uniform. Wha is

More information

Cash Flow Valuation Mode Lin Discrete Time

Cash Flow Valuation Mode Lin Discrete Time IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728,p-ISSN: 2319-765X, 6, Issue 6 (May. - Jun. 2013), PP 35-41 Cash Flow Valuaion Mode Lin Discree Time Olayiwola. M. A. and Oni, N. O. Deparmen of Mahemaics

More information

Random Walk with Anti-Correlated Steps

Random Walk with Anti-Correlated Steps Random Walk wih Ani-Correlaed Seps John Noga Dirk Wagner 2 Absrac We conjecure he expeced value of random walks wih ani-correlaed seps o be exacly. We suppor his conjecure wih 2 plausibiliy argumens and

More information

14 Autoregressive Moving Average Models

14 Autoregressive Moving Average Models 14 Auoregressive Moving Average Models In his chaper an imporan parameric family of saionary ime series is inroduced, he family of he auoregressive moving average, or ARMA, processes. For a large class

More information

Matlab and Python programming: how to get started

Matlab and Python programming: how to get started Malab and Pyhon programming: how o ge sared Equipping readers he skills o wrie programs o explore complex sysems and discover ineresing paerns from big daa is one of he main goals of his book. In his chaper,

More information

A New Perturbative Approach in Nonlinear Singularity Analysis

A New Perturbative Approach in Nonlinear Singularity Analysis Journal of Mahemaics and Saisics 7 (: 49-54, ISSN 549-644 Science Publicaions A New Perurbaive Approach in Nonlinear Singulariy Analysis Ta-Leung Yee Deparmen of Mahemaics and Informaion Technology, The

More information

Existence of multiple positive periodic solutions for functional differential equations

Existence of multiple positive periodic solutions for functional differential equations J. Mah. Anal. Appl. 325 (27) 1378 1389 www.elsevier.com/locae/jmaa Exisence of muliple posiive periodic soluions for funcional differenial equaions Zhijun Zeng a,b,,libi a, Meng Fan a a School of Mahemaics

More information

An Invariance for (2+1)-Extension of Burgers Equation and Formulae to Obtain Solutions of KP Equation

An Invariance for (2+1)-Extension of Burgers Equation and Formulae to Obtain Solutions of KP Equation Commun Theor Phys Beijing, China 43 2005 pp 591 596 c Inernaional Academic Publishers Vol 43, No 4, April 15, 2005 An Invariance for 2+1-Eension of Burgers Equaion Formulae o Obain Soluions of KP Equaion

More information

Chapter 8 The Complete Response of RL and RC Circuits

Chapter 8 The Complete Response of RL and RC Circuits Chaper 8 The Complee Response of RL and RC Circuis Seoul Naional Universiy Deparmen of Elecrical and Compuer Engineering Wha is Firs Order Circuis? Circuis ha conain only one inducor or only one capacior

More information

arxiv: v1 [math.gm] 4 Nov 2018

arxiv: v1 [math.gm] 4 Nov 2018 Unpredicable Soluions of Linear Differenial Equaions Mara Akhme 1,, Mehme Onur Fen 2, Madina Tleubergenova 3,4, Akylbek Zhamanshin 3,4 1 Deparmen of Mahemaics, Middle Eas Technical Universiy, 06800, Ankara,

More information

On the probabilistic stability of the monomial functional equation

On the probabilistic stability of the monomial functional equation Available online a www.jnsa.com J. Nonlinear Sci. Appl. 6 (013), 51 59 Research Aricle On he probabilisic sabiliy of he monomial funcional equaion Claudia Zaharia Wes Universiy of Timişoara, Deparmen of

More information

Problemas das Aulas Práticas

Problemas das Aulas Práticas Mesrado Inegrado em Engenharia Elecroécnica e de Compuadores Conrolo em Espaço de Esados Problemas das Aulas Práicas J. Miranda Lemos Fevereiro de 3 Translaed o English by José Gaspar, 6 J. M. Lemos, IST

More information