On external semi-global stochastic stabilization of linear systems with input saturation

Size: px
Start display at page:

Download "On external semi-global stochastic stabilization of linear systems with input saturation"

Transcription

1 1 On exernal semi-global sochasic sabilizaion of linear sysems wih inpu sauraion Anon A. Soorvogel 1,AliSaberi 2 and Siep Weiland 3 Absrac This paper inroduces he noion of exernal semi-global sochasic sabilizaion for linear plans wih sauraing acuaors, driven by a sochasic exernal disurbance, and having random Gaussiandisribued iniial condiions. The aim of his sabilizaion is o conrol such plans by a possibly nonlinear saic sae feedback law ha achieves global asympoic sabiliy in he absence of disurbances, while guaraneeing a bounded variance of he sae vecor for all ime in he presence of disurbances and Gaussian disribued iniial condiions. We repor complee resuls for boh coninuous- and discree-ime open-loop criically sable plans. I. INTRODUCTION Inernal and exernal sabilizaion of linear plans wih acuaors subjec o sauraion has been he subjec of inense renewed ineres among he conrol research communiy for he pas wo decades. The number of recen books and special issues of conrol journals devoed o his subjec maer evidence his inense research focus, see for insance [1], [11], [4], [15] and he references herein. I is now considered a classical fac ha boh coninuous- and discree-ime linear plans wih sauraing acuaors can be globally inernally sabilized if and only if all of he open-loop poles are locaed in he closed lef half plane in he coninuous-ime case), and are in inside 1 Deparmen of Mahemaics, and Compuing Science, Eindhoven Univ. of Technology, P.O. Box 513, 5600 MB Eindhoven, The Neherlands. Deparmen of Elecrical Engineering, Mahemaics and Compuer Science, Delf Univ. of Technology, P.O. Box 5031, 2600 GA Delf, The Neherlands. A.A.Soorvogel@ue.nl. 2 School of Elecrical Engineering and Compuer Science, Washingon Sae Universiy, Pullman, WA , U.S.A. saberi@eecs.wsu.edu. 3 Deparmen of Elecrical Engineering, Eindhoven Universiy of Technology, P.O. Box 513, 5600 MB Eindhoven, The Neherlands, s.weiland@ue.nl.

2 2 and/or on he uni circle in he discree-ime case). These condiions on open-loop plans can be equivalenly saed as a requiremen ha he plan be asympoically null conrollable wih bounded conrol ANCBC). I is also a classical fac ha, in general, global inernal sabilizaion of ANCBC plans requires nonlinear feedback laws. See for insance. [10], [2], [16]. By weakening he noion of global inernal sabilizaion o a semiglobal framework, Saberi e al. showed ha ANCBC plans wih sauraing acuaors can be inernally sabilized in he semiglobal sense using only linear feedback laws. Such a relaxaion from a global o semiglobal framework is, from an engineering sandpoin, boh sensible and aracive. A body of work exiss proposing design mehodologies low-gain, low-high gain, and variaions) for boh coninuousand discree-ime ANCBC plans wih sauraing aciaors. See for insance [3], [17], [6], [5] Wih regard o exernal sabilizaion, he picure is complicaed. Unlike linear plans, inernal sabiliy does no necessarily guaranee exernal sabiliy when sauraion is presen. Hence, sabilizaion mus be done simulaneously in boh he inernal and exernal sense. There is a body of work ha examines he sandard noion of L p sabiliy for ANCBC plans wih sauraing acuaors when he exernal inpu disurbance) is inpu-addiive and a he same ime requiring eiher global or semi global inernal sabiliy in he absence of disurbnce, see [7]. For he more general case in which he disurbance is no necessarily inpu-addiive, here are surprising resuls in he lieraure which poin o he complexiy of he noion of exernal sabiliy for ANCBC plans wih sauraing acuaors. Noable among hese is he resul from [13] of Soorvogel e.al., ha for a double inegraor wih sauraing acuaors, any linear inernally sabilizing conrol law can achieve L p sabiliy wihou finie gain for p [1, 2], bu no linear inernally sabilizing conrol law can achieve L p sabiliy for any p 2, ]. Moreover, considering he recenly developed noion of inpu-o-sae sabiliy ISS) as a framework for simuaneous exernal and inernal sabiliy, such a double-inegraor sysem canno achieve ISS sabiliy wih a linear feedback law, see [8]. All hese resuls poin oward he delicacy of he exernal sabiliy concep for linear plans wih sauraing acuaors. Moreover, his lieraure also poins ou ha, for susained disurbances, he induced L norm does no yield a suiable problem formulaion. I conains mahemaical funcions ha are no reasonable models for disurbances, and i is impossible o ge a good sable response if we use hese funcions as models for disurbance. In fac one would need o consider a class of sensible susained disurbances from an engineering poin of view. All hese consideraions lead us o believe ha a suiable noion of exernal

3 3 sabiliy for linear plans wih sauraing acuaors, and indeed for general nonlinear sysems in he presence of disurbances and nonzero iniial condiions, is ye o be developed. In his paper we look a he simulaneous exernal and inernal sabilizaion of ANCBC plans wih sauraing acuaors when he exernal inpu is a sochasic disurbance. Specifically, we consider a linear ime-invarian sysem subjec o inpu sauraion, sochasic exernal disurbances and random Gaussian disribued iniial condiions, independen of he exernal disurbances. The aim will be o conrol his sysem by a possibly nonlinear saic sae feedback law ha achieves global asympoic sabiliy in he absence of disurbances, while guaraneeing a bounded variance of he sae vecor for all ime. This problem seems a naural exension of he resuls in [9], [12]. In [14], a resul was claimed. In his paper we give a complee proof for he case ha he sysem is neurally sable, i.e., he eigenvalues of he sysem marix of he coninuous-ime discree-ime) sysem are in he closed lef half plane open uni disc) and he eigenvalues on he imaginary axis uni circle) have equal algebraic and geomeric mulipliciy. The paper is organized as follows. For boh discree- and coninuous-ime sysems, a formal problem formulaion and he main resuls are saed in Secion II. A formal proof of he discreeime case for neurally sable sysems is given in Secion III. A discussion on he implicaions of he main resuls, some exensions and conclusions are colleced in Secion V. II. PROBLEM FORMULATION AND MAIN RESULTS A. The discree-ime case In discree ime, we consider sysems of he form xk +1)=Axk)+Bσuk)) + Ewk) 1) where he sae x, he conrol u and he disurbance w are vecor-valued signals of dimension n, m and l, respecively. Here, k Z +, w is a whie noise sochasic process wih variance Q and mean 0, he iniial condiion x 0 of 1) is a Gaussian random vecor independen of wk) for all k 0. Moreover σ is he sandard sauraion funcion given by: 1 u< 1 σu) = u 1 u 1 1 u>1

4 4 An admissible feedback is a nonlinear feedback of he form uk) =fxk)) 2) where f : R n R m is a coninuous map wih f0) = 0. We herefore consider nonlinear saic sae feedbacks. We will be ineresed in he following problem. Problem II.1 Given he sysem 1), he semiglobal exernal sochasic sabilizaion problem is o find an admissible feedback 2) such ha he following properies hold: 1) in he absence of he exernal inpu w, he equilibrium poin x =0of he conrolled sysem 1)-2) is globally asympoically sable. 2) he variance Varxk)) of he sae of he conrolled sysem 1)-2) is bounded over k 0. The fac ha conrollers exis ha achieve global asympoic sabiliy in he absence of disurbances as described in condiion 1) is well-known. The main objecive of his paper is o look a he addiional requiremen on he variance. The following is he main resul of his paper for discree ime sysems: Theorem II.2 Consider he sysem 1) and suppose ha A, B) is sabilizable while A is neurally sable, i.e. he eigenvalues of A are in he closed uni disc and he eigenvalues on he uni circle have equal geomeric and algebraic mullipliciy. In he above case, here exiss a linear feedback which solves he semiglobal exernal sochasic sabilizaion problem as defined in Problem II.1. We claim ha he condiion above is only sufficien and ha we can weaken he condiion of A neurally sable o he requiremen ha all eigenvalues of A are in he closed uni disc. However, his does, in general, require a nonlinear feedback. Conjecure II.1 Consider he sysem 1). Then for any given variance Q of he whie noise sochasic process w here exiss a feedback 2) which solves he semiglobal exernal sochasic sabilizaion problem as defined in Problem II.1 if and only if A, B) is sabilizable while he eigenvalues of A are in he closed uni disc.

5 5 B. The coninuous-ime case In coninuous ime we consider he sochasic differenial equaion dx) = Ax)d + Bσu))d + Edw) 3) where he sae x, he conrol u and he disurbance w are vecor-valued signals of dimension n, m and l, respecively. Here w is a Wiener process a process of l independen Brownian moions) wih mean 0 and rae Q, hais, Var[W )] = Q and he iniial condiion x 0 of 3) is a Gaussian random vecor which is independen of w. Is soluion x is rigorously defined hrough Wiener inegrals and is a Gauss-Markov process. Like in he discree ime case, σ denoes he sandard sauraion funcion and admissible feedbacks are possibly nonlinear saic sae feedbacks of he form 2) where f is a Lipschiz-coninuous mapping wih f0) = 0. Problem II.3 Given he sysem 3), he semiglobal exernal sochasic sabilizaion problem is o find an admissible feedback 2) such ha he following properies hold: 1) in he absence of he exernal inpu w, he equilibrium poin x =0of he conrolled sysem 3)-2) is globally asympoically sable. 2) he variance Varx)) of he sae of he conrolled sysem 3)-2) is bounded over 0. Like in he discree ime, he fac ha conrollers exis ha achieve global asympoic sabiliy in he absence of disurbances as described in condiion 1) is well-known. The main objecive of his paper is o look a he addiional requiremen on he variance. The following is he main resul of his paper for coninuous ime sysems: Theorem II.4 Consider he sysem 3) and suppose ha A, B) is sabilizable while A is neurally sable, i.e. he eigenvalues of A are in he closed lef half plane and he eigenvalues on he imaginary axis have equal geomeric and algebraic mullipliciy. Then for any given rae Q of he Wiener process w, here exiss a linear feedback which solves he semiglobal exernal sochasic sabilizaion problem as defined in Problem II.3.

6 6 We claim ha he condiion above is only sufficien and ha we can weaken he condiion of A neurally sable o he requiremen ha all eigenvalues of A are in he closed lef half plane. However, his does, in general, require a nonlinear feedback. Conjecure II.2 Consider he sysem 3). Then for any given variance Q of he whie noise sochasic process w) here exiss a feedback 2) which solves he semiglobal exernal sochasic sabilizaion problem as defined in Problem II.3 if and only if A, B) is sabilizable while he eigenvalues of A are in he closed lef half plane. III. PROOFS FOR THE DISCRETE-TIME CASE We firs presen a lile lemma ha we need laer: Lemma III.1 For any a R m and b R m and n>1 we have: σa + b) 2 ) 1 1 2n σa) 2 8n b 2 where σ is he sandard sauraion funcion. Proof: I is easily verified ha: σa + b) σa) b where he inequaliy is componenwise and b R m is defined as he absolue value of each componen of b. This implies: Using hen yields he required resul. σa + b) 2 σa) 2 2 b T σa)+ b 2 b T σa) 4n b n σa) 2 Using a simple basis ransformaion we can assume wihou loss of generaliy ha: A = A 11 0, B = B 1, E = E 1 0 A 22 B 2 E 2 wih A 11 asympoically sable while A T 22 A 22 = I. This is guaraneed by he fac ha A is neurally sable. Nex we consider he subsysem: x 2 +1)=A 22 x 2 )+B 2 u)+e 2 w) 4)

7 7 If we find an admissible feedback of he form u) =fx 2 )) 5) which solves he semiglobal sochasic sabilizaion problem for he sysem 4) hen i is easily verified ha his conroller also solves he semiglobal sochasic sabilizaion problem for he original sysem 1). In oher words we can resric aenion o he sysem 4). We will use a conroller of he form: uk) = B2 T A 22x 2 k) and we will esablish ha for any given variance Q of he whie noise here exiss small enough such ha he resuling sysem has a bounded variance for he sae. The choice of is independen of he variance of he iniial condiion x 0. Proof: [Proof of Theorem II.2] To keep noaion simple, we assume he original sysem saisfies A T A = I and we consider a feedback of he form: uk) = B T Axk) 6) The general case of he heorem can hen easily be esablished using he reducion o he sysem 4) as described earlier. If he sysem is no affeced by noise, i is well-known ha he feedback 6) achieves asympoic sabiliy which is also easily verified by noing ha V x) =x T x is a suiable Lyapunov funcion. To sudy he effec of he noise, we firs noe ha due o he Markov propery of he sae, we have: E [xk + r) T xk + r) xk + s)] = E [E [xk + r) T xk + r) xk + )] xk + s)] for r s. Nex we consider: E [xk +1) T xk +1) xk)]

8 8 Using he feedback 6), he dynamics 1) and he fac ha A T A = I we obain: E [xk +1) T xk +1) xk)] = xk) T xk) 2xk) T A T BσB T Axk)) + σb T Axk)) T B T BσB T Axk)) + race EQE T xk) T xk) 2 σbt Axk)) T σb T Axk)) + σb T Axk)) T B T BσB T Axk)) + race EQE T xk) T xk) 1 σbt Axk)) 2 + race EQE T 7) provided is small enough such ha B T B<I. Nex, we consider: [ σbt E A xk +1)) ] 2 xk) for some ineger 1. Wege: [ σbt E A xk +1)) ] 2 xk) [ σbt = E A +1 xk) B T A BσB T Axk)) + B T A Ewk)) ] 2 xk) 1 1 2n ) σbt A +1 xk) 2 8n E [ B T A BσB T Axk)) +B T A Ewk)) 2 xk) ] 1 1 2n ) σbt A +1 xk)) 2 8n 2 B T A B 2 σb T Axk)) B T A 2 race EQE T where he firs inequaliy follows from Lemma III.1. If we choose such ha B T B<I[as used earlier in deriving 7)] and 8n B T A 2 1, 8n B T A B 2 1 for all, hen i is easy o see ha for all < and any posiive ineger we have: [ σbt E A xk +1)) ] 2 xk) 1 1 2n ) σbt A +1 xk)) 2 3/2 σb T Axk)) 2 3/2 race EQE T The following lemma presens a crucial inequaliy:

9 9 Lemma III.2 We have for small enough ha E [xk + s) T xk + s) xk)] xk) T xk) 1 s =1 1 ) σbt A xk)) 2 + s 1+ s ) race EQE T 8) 2n n Proof: We will esablish his inequaliy recursively. Noe ha we already esablished his inequaliy for s =1given 7). Nex, consider any s>1 and assume he above inequaliy is rue for s and any k). Then we obain: E [xk + s +1) T xk + s +1) xk)] = E [E [xk + s +1) T xk + s +1) xk +1)] xk)] [ E xk +1) T xk +1) 1 s 1 ) σbt A xk +1)) 2 2n =1 +s 1+ s ) ] race EQE T xk) n Using he previously obained inequaliy we have: [ s E 1 ) σbt A xk +1)) 2 xk)] 2n =1 1 1 ) [ s 1 ) σbt A +1 xk)) ] 2 s 3/2 σb T Axk)) 2n 2n =1 s 3/2 race EQE T s 1 +1 ) σbt A +1 xk)) 2 s 3/2 σb T Axk)) s 3/2 race EQE T 2n =1

10 10 We hen obain: E [xk + s +1) T xk + s +1) xk)] = E [E [xk + s +1) T xk + s +1) xk +1)] xk)] xk) T xk) 1 σbt Axk)) 2 + race EQE T 1 s 1 +1 ) σbt A +1 xk)) 2 2n =1 xk) T xk) 1 + s σb T Axk)) 2 + s race EQE T + s s+1 =1 provided < is such ha 1+ s ) race EQE T n 1 ) σbt A xk)) 2 +s +1) 1+ s +1 ) race EQE T 2n n 1 n, s3/2 1 2n. Nex we noe ha since A, B) is conrollable while A is inverible, he marix B T A B T A 2. B T A n is injecive which implies ha here exiss α he smalles singular value of his marix divided by n) such ha for any x here exiss a posiive ineger n such ha and hence B T A x α x 1 { } σbt A x) 2 min α 2 x 2, 1 Bu hen 8) yields: } E [xk + n) T xk + n) xk)] xk) T xk) {α 1 min 2 xk) 2, 1 +2nrace EQE T 2 We define M = 1 2 α 2

11 11 and yk) =max{ xk) T xk),m} The following lemma presens anoher crucial resul: Lemma III.3 We have: E [yk + n) yk)] yk) 9) provided is small enough such ha provided is small enough such ha he inequaliy of Lemma III.2 is valid for s =1,...,n and: 1 2 > 2n race EQET and α 2 < 1 Moreover, have: E[yk + pn)] E[yk)], p =1, 2,... 10) Proof: This follows from he fac ha if yk) =xk) T xk) and hence xk) T xk) M we E [xk + n) T xk + n) xk)] } yk) {α 1 min 2 xk) 2, 1 +2nrace EQE T 2 yk) 1 +2nrace EQET 2 yk) while for he case yk) =M and hence xk) T xk) M we have: E [xk + n) T xk + n) xk)] } xk) T xk) {α 1 min 2 xk) 2, 1 +2nrace EQE T α2 )xk) T xk)+2n race EQE T α2 )yk)+2n race EQE T yk) 1 2 α2 M +2n race EQE T yk) bu hen using ha yk) M): E [yk + n) xk)] yk)

12 12 and finally, E [yk + n) yk)] = E [E [yk + n) xk)] yk)] E [yk) yk)] = yk) Similarlywherewehaveobecarefulsinceyk) is no Markov), E [yk +2n) yk)] = E [E [yk +2n) xk + n),xk)] yk)] = E [E [E [yk +2n) xk + n),xk)] xk)] yk)] = E [E [E [yk +2n) xk + n)] xk)] yk)] E [E [yk + n) xk)] yk)] E [yk) yk)] = yk) Using a simple recursion, we obain 10). We have, using 7): E[xk +1) T xk +1) xk)] xk) T xk) + race EQE T which yields: E[xk) T xk)] Var[x0)] + k race EQE T for k =1,...,n 1 and hence: E[yk)] max { M,Var[x0)] + k race EQE T } for k =1,...,n 1. This implies ha he expecaion of y0),...,yn 1) is bounded. Using 10) we find ha he expecaion of yk) is bounded in k. Usinghayk) xk) T xk) his yields ha he variance of xk) is bounded in k. IV. PROOFS FOR THE CONTINUOUS-TIME CASE Using a simple basis ransformaion we can assume wihou loss of generaliy ha: A = A 11 0, B = B 1, E = E 1 0 A 22 B 2 E 2

13 13 wih A 11 asympoically sable while A 22 + A T 22 =0. This is guaraneed by he fac ha A is neurally sable. Nex we consider he subsysem: dx 2 ) =A 22 x 2 )d + B 2 σu))d + E 2 dw) 11) If we find an admissible feedback of he form u) =fx 2 )) 12) which solves he semiglobal exernal sochasic sabilizaion problem for he sysem 11) hen i is easily verified ha his conroller also solves he semiglobal sochasic sabilizaion problem for he original sysem 3). In oher words we can resric aenion o he sysem 11). We will use a conroller of he form: u) = B2x T 2 ) and we will esablish ha for any given rae Q of he Wiener process here exiss small enough such ha he resuling sysem has a bounded variance for he sae. The choice of is independen of he variance of he iniial condiion x 0. Proof: [Proof of Theorem II.4] To keep noaion simple, we assume he original sysem saisfies A + A T =0and we consider a feedback of he form: u) = B T x) 13) The general case of he heorem can hen easily be esablished using he reducion o he sysem 11) as described earlier. If he sysem is no affeced by noise, i is well-known ha he feedback 13) achieves asympoic sabiliy which is also easily verified by noing ha V x) =x T x is a suiable Lyapunov funcion. To sudy he effec of he noise, we firs noe ha: τ τ xτ) =e Aτ ) x)+ e Aτ s) Bσus))ds + e Aτ s) Edws) = e Aτ ) x)+v τ) We noe ha σu) is bounded by 1 and he second inegral has bounded variance for all τ such ha τ [, +1]. This implies ha here exiss a consan N such ha: Var[v τ)] <N 14)

14 14 for all τ such ha τ [, +1]. Nex, we noe ha σb T xτ)) 2 = σ B T e Aτ ) x)+b T v τ) ) σbt e Aτ ) x)) B 2 v τ) 2 15) where we applied Lemma III.1 wih n =1. Nex we noe ha: dx T x =2x T Bσu) + race E T E)d +2x T Edw using ha A + A T =0and Iô s lemma. This yields ha x +1) T x +1) x) T x) = +1 2xτ) T Bσuτ))dτ + race E T E xτ) T Edwτ) and hence: E[x +1) T x +1) x)] On he oher hand, we have: = x) T x)+e +1 2xτ) T Bσuτ))dτ x) + race E T E Combining he las wo equaions we obain: 2xτ) T Bσuτ)) = 2xτ) T Bσ B T xτ)) 2 σbt xτ)) 2 E[x +1) T x +1) x)] x) T x) + race E T E 2 E +1 σb T xτ)) 2 dτ x) Using 14), 15) and he fac ha x) and v τ) are independen for τ [, +1] we obain a crucial inequaliy: +1 E[x +1) T x +1) x)] x) T x)+v 1 σb T e Aτ ) x)) 2 dτ 16)

15 15 where V =16 B 2 N + race E T E. Lemma IV.1 Assume ha A, B) is conrollable. Then here exis suiably defined consans α and R such ha for all < we have ha: 1 +1 σbt e Aτ ) x)) 2 dτ min{α x) 2, R} Proof: If he sauraion does no ge acivaed and hence we ge: +1 σbt e Aτ ) x)) 2 dτ = 2 x) T Qx) where: Q = 1 e Aτ BB T e Aτ dτ Conrollabiliy of A, B) implies ha Q is inverible and hence we obain: 0 x)qx) α x) 2 where α is he smalles singular value of Q. Hence we obain: 1 +1 σbt e Aτ ) x)) 2 dτ α x) 2 17) On he oher hand, in case he sauraion does ge acivaed hen i is easily seen ha here mus exis an N independen of ) such ha: We obain: +1 We consider he inegral: x) > N σbt e Aτ ) x)) 2 dτ +1 σ B T e +1 σ B T e Aτ ) Nx) x) Aτ ) Nx) x) ) 2 dτ ) 2 dτ 18) In case sauraion does no arise in his inegral we immediaely find ha he erm is larger or equal o 2 α Nx) x) 2 αn 2

16 16 using he same echniques as for he case ha he original inegral did no conain a erm which sauraes. Bu hen we find 1 +1 σbt e Aτ ) x)) 2 dτ 1 αn 2 19) Finally we look a he case where he sauraion in he inegral 18) does ge acivaed. In ha case we noe ha here exiss τ such ha: ) σ B T Aτ ) Nx) 2 e 1 x) Bu he derivaive wih respec o τ of his erm is bounded by some consan γ independen of our choice for τ [, +1] and independen of our choice for x). Thisimpliesha This yields σ B T e Aτ ) Nx) σbt e Aτ ) x)) 2 dτ 1 x) +1 ) 2 dτ 1 γ σ B T e Aτ ) Nx) x) ) 2 dτ 1 γ Combining 17) wih 19) and 20), he resul of he lemma follows immediaely for R =min{n 2, 1 γ }. 20) The above lemma yields, wih he help of 16), ha: E[x +1) T x +1) x)] x) T x) min{α x) 2, R } + V We define: and M = R 2 α y) =max{ x) T x),m} The following lemma presens anoher crucial resul:

17 17 Lemma IV.2 We have: E [y +1) y)] y) 21) provided is small enough such ha: < R V and α < 1. Moreover, have: Proof: E[y + p)] E[y)], p =1, 2,... 22) This follows from he fac ha if y) =x) T x) and hence x) T x) M we E [x +1) T x +1) x)] { } y) min α x) 2, R + V y) R + V y) while for he case y) =M and hence x) T x) M we have: bu hen using ha y) M): and finally, E [x +1) T x +1) x)] { } x) T x) min α x) 2, R + V 1 α)x) T x)+v 1 α)y)+v y) αm + V y) E [y +1) x)] y) E [y +1) y)] = E [E [y +1) x)] y)] E [y) y)] = y)

18 18 The inequaliy 22) hen follows using he same argumens as in he proof of Lemma III.3. I is easy o esablish ha: E[y)] is bounded on he inerval [0, 1]. Combined wih 22) we find ha he variance of y) is bounded. Using ha y) x) T x) his yields ha he variance of x) is bounded in. V. DISCUSSION AND CONCLUSIONS This paper is he coninuaion of our earlier effors oward undersanding and suiably formulaing he noion of simulaneous exernal and inernal sabilizaion for linear plans wih sauraing acuaors. Here we have shown ha i is possible o conrol an open-loop criically sable linear plan in coninuous or discree ime, which is subjec o inpu sauraion, sochasic exernal disurbances and Gaussian-disribued iniial condiions, via a linear sae feedback law ha achieves global asympoic sabiliy of he closed-loop sysem in he absence of disurbances, while guaraneeing bounded variance of he sae vecor for all ime in he presence of disurbances. We have ermed his conrol problem he exernal semi-global sabilizaion problem. Furhermore, we have conjecured ha such a problem is also solvable for open-loop criically unsable linear plans, albei via a nonlinear sae feedback law. REFERENCES [1] D.S. BERNSTEIN AND A.N. MICHEL, Gues Eds., Special Issue on sauraing acuaors, In. J. Robus & Nonlinear Conrol, 55), 1995, pp [2] A.T. FULLER, In-he-large sabiliy of relay and sauraing conrol sysems wih linear conroller, In. J. Conr., 104), 1969, pp [3] P. HOU, A. SABERI, Z. LIN, AND P. SANNUTI, Simulaneously exernal and inernal sabilizaion for coninuous and discree-ime criically unsable sysems wih sauraing acuaors, Auomaica, 3412), 1998, pp [4] V. KAPILA AND G. GRIGORIADIS, Eds., Acuaor sauraion conrol, Marcel Dekker, [5] Z. LIN AND A. SABERI, Semi-global exponenial sabilizaion of linear discree-ime sysems subjec o inpu sauraion via linear feedbacks, Sys. & Conr. Leers, 24, 1995, pp [6] Z. LIN, A. SABERI, AND A.A. STOORVOGEL, Semi-global sabilizaion of linear discree-ime sysems subjec o inpu sauraion via linear feedback - an ARE-based approach, IEEE Trans. Au. Conr., 418), 1996, pp [7] A. SABERI, P. HOU, AND A.A. STOORVOGEL, On simulaneous global exernal and global inernal sabilizaion of criically unsable linear sysems wih sauraing acuaors, IEEE Trans. Au. Conr., 456), 2000, pp [8] G. SHI AND A. SABERI, On he inpu-o-sae sabiliy ISS) of a double inegraor wih sauraed linear conrol laws, in Proceedings of he American Conrol Conference, Anchorage, Alaska, 2002, pp

19 19 [9] G. SHI, A.SABERI, AND A.A. STOORVOGEL, OnL p l p) performance wih global inernal sabiliy for linear sysems wih acuaors subjec o ampliude and rae sauraion, in American Conrol Conference, Chicago, IL, 2000, pp [10] E.D. SONTAG AND H.J. SUSSMANN, Nonlinear oupu feedback design for linear sysems wih sauraing conrols, in Proc. 29h CDC, Honolulu, 1990, pp [11] A. SABERI AND A.A. STOORVOGEL, Gues Eds., Special Issue on conrol problems wih consrains, In. J. Robus & Nonlinear Conrol, 910), 1999, pp [12] A.A. STOORVOGEL,A.SABERI, AND G. SHI, On achieving L p l p) performance wih global inernal sabiliy for linear plans wih sauraing acuaors, in Proc. 38h CDC, Phoenix, AZ, 1999, pp [13] A.A. STOORVOGEL, G. SHI, AND A. SABERI, Exernal sabiliy of a double inegraor wih sauraed linear conrol laws, Dynamics of Coninuous Discree and Impulsive Sysems, Series B: Applicaions & Algorihms, 114-5), 2004, pp [14] A.A. STOORVOGEL,S.WEILAND, AND A. SABERI, On sabilizaion of linear sysems wih sochasic disurbances and inpu sauraion, in Proc. 43rd CDC, The Bahamas, 2004, pp [15] S. TARBOURIECH AND G. GARCIA, Eds., Conrol of uncerain sysems wih bounded inpus, vol. 227 of Lecure noes in conrol and informaion sciences, Springer Verlag, [16] A.R. TEEL, Global sabilizaion and resriced racking for muliple inegraors wih bounded conrols, Sys. & Conr. Leers, 183), 1992, pp [17] Y. YANG, Global sabilizaion of linear sysems wih bounded feedback, PhD hesis, Rugers Universiy, New Brunswick, 1993.

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

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

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

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

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

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

arxiv: v1 [math.ca] 15 Nov 2016

arxiv: v1 [math.ca] 15 Nov 2016 arxiv:6.599v [mah.ca] 5 Nov 26 Counerexamples on Jumarie s hree basic fracional calculus formulae for non-differeniable coninuous funcions Cheng-shi Liu Deparmen of Mahemaics Norheas Peroleum Universiy

More information

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3 and d = c b - b c c d = c b - b c c This process is coninued unil he nh row has been compleed. The complee array of coefficiens is riangular. Noe ha in developing he array an enire row may be divided or

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

Recursive Least-Squares Fixed-Interval Smoother Using Covariance Information based on Innovation Approach in Linear Continuous Stochastic Systems

Recursive Least-Squares Fixed-Interval Smoother Using Covariance Information based on Innovation Approach in Linear Continuous Stochastic Systems 8 Froniers in Signal Processing, Vol. 1, No. 1, July 217 hps://dx.doi.org/1.2266/fsp.217.112 Recursive Leas-Squares Fixed-Inerval Smooher Using Covariance Informaion based on Innovaion Approach in Linear

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

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

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

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

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

On a Fractional Stochastic Landau-Ginzburg Equation

On a Fractional Stochastic Landau-Ginzburg Equation Applied Mahemaical Sciences, Vol. 4, 1, no. 7, 317-35 On a Fracional Sochasic Landau-Ginzburg Equaion Nguyen Tien Dung Deparmen of Mahemaics, FPT Universiy 15B Pham Hung Sree, Hanoi, Vienam dungn@fp.edu.vn

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

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

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

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

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

LAPLACE TRANSFORM AND TRANSFER FUNCTION

LAPLACE TRANSFORM AND TRANSFER FUNCTION CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dep. of Chemical and Biological Engineering 5-1 Road Map of he Lecure V Laplace Transform and Transfer funcions

More information

A Note on the Equivalence of Fractional Relaxation Equations to Differential Equations with Varying Coefficients

A Note on the Equivalence of Fractional Relaxation Equations to Differential Equations with Varying Coefficients mahemaics Aricle A Noe on he Equivalence of Fracional Relaxaion Equaions o Differenial Equaions wih Varying Coefficiens Francesco Mainardi Deparmen of Physics and Asronomy, Universiy of Bologna, and he

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

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

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

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

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

DISCRETE GRONWALL LEMMA AND APPLICATIONS

DISCRETE GRONWALL LEMMA AND APPLICATIONS DISCRETE GRONWALL LEMMA AND APPLICATIONS JOHN M. HOLTE MAA NORTH CENTRAL SECTION MEETING AT UND 24 OCTOBER 29 Gronwall s lemma saes an inequaliy ha is useful in he heory of differenial equaions. Here is

More information

Chapter 6. Systems of First Order Linear Differential Equations

Chapter 6. Systems of First Order Linear Differential Equations Chaper 6 Sysems of Firs Order Linear Differenial Equaions We will only discuss firs order sysems However higher order sysems may be made ino firs order sysems by a rick shown below We will have a sligh

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

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

dy dx = xey (a) y(0) = 2 (b) y(1) = 2.5 SOLUTION: See next page

dy dx = xey (a) y(0) = 2 (b) y(1) = 2.5 SOLUTION: See next page Assignmen 1 MATH 2270 SOLUTION Please wrie ou complee soluions for each of he following 6 problems (one more will sill be added). You may, of course, consul wih your classmaes, he exbook or oher resources,

More information

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid

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

The L p -Version of the Generalized Bohl Perron Principle for Vector Equations with Infinite Delay

The L p -Version of the Generalized Bohl Perron Principle for Vector Equations with Infinite Delay Advances in Dynamical Sysems and Applicaions ISSN 973-5321, Volume 6, Number 2, pp. 177 184 (211) hp://campus.ms.edu/adsa The L p -Version of he Generalized Bohl Perron Principle for Vecor Equaions wih

More information

Solutions to Assignment 1

Solutions to Assignment 1 MA 2326 Differenial Equaions Insrucor: Peronela Radu Friday, February 8, 203 Soluions o Assignmen. Find he general soluions of he following ODEs: (a) 2 x = an x Soluion: I is a separable equaion as we

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

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality Marix Versions of Some Refinemens of he Arihmeic-Geomeric Mean Inequaliy Bao Qi Feng and Andrew Tonge Absrac. We esablish marix versions of refinemens due o Alzer ], Carwrigh and Field 4], and Mercer 5]

More information

Two Coupled Oscillators / Normal Modes

Two Coupled Oscillators / Normal Modes Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own

More information

The expectation value of the field operator.

The expectation value of the field operator. The expecaion value of he field operaor. Dan Solomon Universiy of Illinois Chicago, IL dsolom@uic.edu June, 04 Absrac. Much of he mahemaical developmen of quanum field heory has been in suppor of deermining

More information

BBP-type formulas, in general bases, for arctangents of real numbers

BBP-type formulas, in general bases, for arctangents of real numbers Noes on Number Theory and Discree Mahemaics Vol. 19, 13, No. 3, 33 54 BBP-ype formulas, in general bases, for arcangens of real numbers Kunle Adegoke 1 and Olawanle Layeni 2 1 Deparmen of Physics, Obafemi

More information

Approximation Algorithms for Unique Games via Orthogonal Separators

Approximation Algorithms for Unique Games via Orthogonal Separators Approximaion Algorihms for Unique Games via Orhogonal Separaors Lecure noes by Konsanin Makarychev. Lecure noes are based on he papers [CMM06a, CMM06b, LM4]. Unique Games In hese lecure noes, we define

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

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004 ODEs II, Lecure : Homogeneous Linear Sysems - I Mike Raugh March 8, 4 Inroducion. In he firs lecure we discussed a sysem of linear ODEs for modeling he excreion of lead from he human body, saw how o ransform

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

Properties Of Solutions To A Generalized Liénard Equation With Forcing Term

Properties Of Solutions To A Generalized Liénard Equation With Forcing Term Applied Mahemaics E-Noes, 8(28), 4-44 c ISSN 67-25 Available free a mirror sies of hp://www.mah.nhu.edu.w/ amen/ Properies Of Soluions To A Generalized Liénard Equaion Wih Forcing Term Allan Kroopnick

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

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

Continuous Time. Time-Domain System Analysis. Impulse Response. Impulse Response. Impulse Response. Impulse Response. ( t) + b 0.

Continuous Time. Time-Domain System Analysis. Impulse Response. Impulse Response. Impulse Response. Impulse Response. ( t) + b 0. Time-Domain Sysem Analysis Coninuous Time. J. Robers - All Righs Reserved. Edied by Dr. Rober Akl 1. J. Robers - All Righs Reserved. Edied by Dr. Rober Akl 2 Le a sysem be described by a 2 y ( ) + a 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

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

Undetermined coefficients for local fractional differential equations

Undetermined coefficients for local fractional differential equations Available online a www.isr-publicaions.com/jmcs J. Mah. Compuer Sci. 16 (2016), 140 146 Research Aricle Undeermined coefficiens for local fracional differenial equaions Roshdi Khalil a,, Mohammed Al Horani

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

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

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

Existence of positive solution for a third-order three-point BVP with sign-changing Green s function

Existence of positive solution for a third-order three-point BVP with sign-changing Green s function Elecronic Journal of Qualiaive Theory of Differenial Equaions 13, No. 3, 1-11; hp://www.mah.u-szeged.hu/ejqde/ Exisence of posiive soluion for a hird-order hree-poin BVP wih sign-changing Green s funcion

More information

Robust Stability and Stabilization for Singular Systems With State Delay and Parameter Uncertainty

Robust Stability and Stabilization for Singular Systems With State Delay and Parameter Uncertainty 11 IEEE RANSACIONS ON AUOMAIC CONROL, VOL. 47, NO. 7, JULY Now, we choose 8= 1; ; 1 9= 19; 5 5; 7 he corresponding Riccai equaion (13) admis he maximal soluion 5= 1; ; 7 : By heorem 2, all opimal conrols

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

SOME MORE APPLICATIONS OF THE HAHN-BANACH THEOREM

SOME MORE APPLICATIONS OF THE HAHN-BANACH THEOREM SOME MORE APPLICATIONS OF THE HAHN-BANACH THEOREM FRANCISCO JAVIER GARCÍA-PACHECO, DANIELE PUGLISI, AND GUSTI VAN ZYL Absrac We give a new proof of he fac ha equivalen norms on subspaces can be exended

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

Linear Quadratic Regulator (LQR) - State Feedback Design

Linear Quadratic Regulator (LQR) - State Feedback Design Linear Quadrai Regulaor (LQR) - Sae Feedbak Design A sysem is expressed in sae variable form as x = Ax + Bu n m wih x( ) R, u( ) R and he iniial ondiion x() = x A he sabilizaion problem using sae variable

More information

LECTURE 1: GENERALIZED RAY KNIGHT THEOREM FOR FINITE MARKOV CHAINS

LECTURE 1: GENERALIZED RAY KNIGHT THEOREM FOR FINITE MARKOV CHAINS LECTURE : GENERALIZED RAY KNIGHT THEOREM FOR FINITE MARKOV CHAINS We will work wih a coninuous ime reversible Markov chain X on a finie conneced sae space, wih generaor Lf(x = y q x,yf(y. (Recall ha q

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

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is UNIT IMPULSE RESPONSE, UNIT STEP RESPONSE, STABILITY. Uni impulse funcion (Dirac dela funcion, dela funcion) rigorously defined is no sricly a funcion, bu disribuion (or measure), precise reamen requires

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

Reading from Young & Freedman: For this topic, read sections 25.4 & 25.5, the introduction to chapter 26 and sections 26.1 to 26.2 & 26.4.

Reading from Young & Freedman: For this topic, read sections 25.4 & 25.5, the introduction to chapter 26 and sections 26.1 to 26.2 & 26.4. PHY1 Elecriciy Topic 7 (Lecures 1 & 11) Elecric Circuis n his opic, we will cover: 1) Elecromoive Force (EMF) ) Series and parallel resisor combinaions 3) Kirchhoff s rules for circuis 4) Time dependence

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

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

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

Stabilization of Exponentially Unstable Linear Systems with Saturating Actuators

Stabilization of Exponentially Unstable Linear Systems with Saturating Actuators IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 46, NO. 6, JUNE 2 973 Sabilizaion of Exponenially Unsable Linear Sysems wih Sauraing Acuaors Tingshu Hu, Zongli Lin, and Li Qiu Absrac We sudy he problem of

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

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

Chapter 1 Fundamental Concepts

Chapter 1 Fundamental Concepts Chaper 1 Fundamenal Conceps 1 Signals A signal is a paern of variaion of a physical quaniy, ofen as a funcion of ime (bu also space, disance, posiion, ec). These quaniies are usually he independen variables

More information

CONTRIBUTION TO IMPULSIVE EQUATIONS

CONTRIBUTION TO IMPULSIVE EQUATIONS European Scienific Journal Sepember 214 /SPECIAL/ ediion Vol.3 ISSN: 1857 7881 (Prin) e - ISSN 1857-7431 CONTRIBUTION TO IMPULSIVE EQUATIONS Berrabah Faima Zohra, MA Universiy of sidi bel abbes/ Algeria

More information

Chapter #1 EEE8013 EEE3001. Linear Controller Design and State Space Analysis

Chapter #1 EEE8013 EEE3001. Linear Controller Design and State Space Analysis Chaper EEE83 EEE3 Chaper # EEE83 EEE3 Linear Conroller Design and Sae Space Analysis Ordinary Differenial Equaions.... Inroducion.... Firs Order ODEs... 3. Second Order ODEs... 7 3. General Maerial...

More information

Sliding Mode Controller for Unstable Systems

Sliding Mode Controller for Unstable Systems S. SIVARAMAKRISHNAN e al., Sliding Mode Conroller for Unsable Sysems, Chem. Biochem. Eng. Q. 22 (1) 41 47 (28) 41 Sliding Mode Conroller for Unsable Sysems S. Sivaramakrishnan, A. K. Tangirala, and M.

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

Four Generations of Higher Order Sliding Mode Controllers. L. Fridman Universidad Nacional Autonoma de Mexico Aussois, June, 10th, 2015

Four Generations of Higher Order Sliding Mode Controllers. L. Fridman Universidad Nacional Autonoma de Mexico Aussois, June, 10th, 2015 Four Generaions of Higher Order Sliding Mode Conrollers L. Fridman Universidad Nacional Auonoma de Mexico Aussois, June, 1h, 215 Ouline 1 Generaion 1:Two Main Conceps of Sliding Mode Conrol 2 Generaion

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

Differential Equations

Differential Equations Mah 21 (Fall 29) Differenial Equaions Soluion #3 1. Find he paricular soluion of he following differenial equaion by variaion of parameer (a) y + y = csc (b) 2 y + y y = ln, > Soluion: (a) The corresponding

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

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

AN OPTIMAL CONTROL PROBLEM FOR SPLINES ASSOCIATED TO LINEAR DIFFERENTIAL OPERATORS

AN OPTIMAL CONTROL PROBLEM FOR SPLINES ASSOCIATED TO LINEAR DIFFERENTIAL OPERATORS CONTROLO 6 7h Poruguese Conference on Auomaic Conrol Insiuo Superior Técnico, Lisboa, Porugal Sepember -3, 6 AN OPTIMAL CONTROL PROBLEM FOR SPLINES ASSOCIATED TO LINEAR DIFFERENTIAL OPERATORS Rui C. Rodrigues,

More information

EXISTENCE OF NON-OSCILLATORY SOLUTIONS TO FIRST-ORDER NEUTRAL DIFFERENTIAL EQUATIONS

EXISTENCE OF NON-OSCILLATORY SOLUTIONS TO FIRST-ORDER NEUTRAL DIFFERENTIAL EQUATIONS Elecronic Journal of Differenial Equaions, Vol. 206 (206, No. 39, pp.. ISSN: 072-669. URL: hp://ejde.mah.xsae.edu or hp://ejde.mah.un.edu fp ejde.mah.xsae.edu EXISTENCE OF NON-OSCILLATORY SOLUTIONS TO

More information

Supplement for Stochastic Convex Optimization: Faster Local Growth Implies Faster Global Convergence

Supplement for Stochastic Convex Optimization: Faster Local Growth Implies Faster Global Convergence Supplemen for Sochasic Convex Opimizaion: Faser Local Growh Implies Faser Global Convergence Yi Xu Qihang Lin ianbao Yang Proof of heorem heorem Suppose Assumpion holds and F (w) obeys he LGC (6) Given

More information

TO our knowledge, most exciting results on the existence

TO our knowledge, most exciting results on the existence IAENG Inernaional Journal of Applied Mahemaics, 42:, IJAM_42 2 Exisence and Uniqueness of a Periodic Soluion for hird-order Delay Differenial Equaion wih wo Deviaing Argumens A. M. A. Abou-El-Ela, A. I.

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

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients Secion 3.5 Nonhomogeneous Equaions; Mehod of Undeermined Coefficiens Key Terms/Ideas: Linear Differenial operaor Nonlinear operaor Second order homogeneous DE Second order nonhomogeneous DE Soluion o homogeneous

More information

Fractional Method of Characteristics for Fractional Partial Differential Equations

Fractional Method of Characteristics for Fractional Partial Differential Equations Fracional Mehod of Characerisics for Fracional Parial Differenial Equaions Guo-cheng Wu* Modern Teile Insiue, Donghua Universiy, 188 Yan-an ilu Road, Shanghai 51, PR China Absrac The mehod of characerisics

More information

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon 3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of

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

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

A Primal-Dual Type Algorithm with the O(1/t) Convergence Rate for Large Scale Constrained Convex Programs

A Primal-Dual Type Algorithm with the O(1/t) Convergence Rate for Large Scale Constrained Convex Programs PROC. IEEE CONFERENCE ON DECISION AND CONTROL, 06 A Primal-Dual Type Algorihm wih he O(/) Convergence Rae for Large Scale Consrained Convex Programs Hao Yu and Michael J. Neely Absrac This paper considers

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

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

Georey E. Hinton. University oftoronto. Technical Report CRG-TR February 22, Abstract

Georey E. Hinton. University oftoronto.   Technical Report CRG-TR February 22, Abstract Parameer Esimaion for Linear Dynamical Sysems Zoubin Ghahramani Georey E. Hinon Deparmen of Compuer Science Universiy oftorono 6 King's College Road Torono, Canada M5S A4 Email: zoubin@cs.orono.edu Technical

More information

Non-uniform Small-gain Theorems for Systems with Unstable Invariant Sets

Non-uniform Small-gain Theorems for Systems with Unstable Invariant Sets Proceedings of he 47h IEEE Conference on Decision and Conrol Cancun, Mexico, Dec. 9-11, 2008 Non-uniform Small-gain Theorems for Sysems wih Unsable Invarian Ses Ivan Y. Tyukin, Erik Seur, Henk Nijmeijer,

More information

Most Probable Phase Portraits of Stochastic Differential Equations and Its Numerical Simulation

Most Probable Phase Portraits of Stochastic Differential Equations and Its Numerical Simulation Mos Probable Phase Porrais of Sochasic Differenial Equaions and Is Numerical Simulaion Bing Yang, Zhu Zeng and Ling Wang 3 School of Mahemaics and Saisics, Huazhong Universiy of Science and Technology,

More information