Stabilization and Performance Analysis for a Class of Switched Systems

Size: px
Start display at page:

Download "Stabilization and Performance Analysis for a Class of Switched Systems"

Transcription

1 43rd IEEE Conference on Decson and Control December 14-17, 2004 Atlants, Paradse Island, Bahamas ThB02.6 Stablzaton and Performance Analyss for a Class of Swtched Systems Le Fang, Ha Ln, and Panos J. Antsakls Abstract Ths paper nvestgates stablty and control desgn problems wth performance analyss for dscrete-tme swtched lnear systems. The swtched Lyapunov functon method s combned wth Fnsler s Lemma to generate varous tests n the enlarged space contanng both the state and ts tme dfference, allowng extra degree of freedom for stablty analyss and control desgn. Two performance measures beng consdered are the decay rate and the nputoutput performance. A new LMI based stablty test for the exstence of swtched Lyapunov functons s frst developed. If a swtched Lyapunov functon exsts, asymptotc stablty of the swtched system also mples ts exponental stablty. An LMI optmzaton problem s then formulated to fnd a bound on the decay rate of the system. To attan the bound, state feedback control gans are desgned. Usng the same framework and the well-known S-procedure, a generalzed suffcent LMI condton s obtaned whch guarantees a γ- performance of the closed-loop swtched systems subject to nput dsturbances. I. INTRODUCTION In recent years, consderable research efforts have been devoted to the study of swtched systems. The motvaton for studyng swtched systems comes from the fact that many practcal systems are nherently multmodal 8, and the fact that some of ntellgent control methods are based on the dea of swtchng between dfferent controllers 9, 11, 18, 20. The exstence of systems that cannot be asymptotcally stablzed by a sngle statc contnuous feedback controller 4 also motvates the study. A survey of basc problems n stablty and desgn of swtched systems s gven n 15. One of the basc problems s to fnd condtons whch guarantee that swtched systems are asymptotcally stable under arbtrary swtchng sequences. Stablty analyss of swtched systems s usually carred out n the Lyapunov framework 3, 9. For swtched lnear systems, stablty under arbtrary swtchng s equvalent to the exstence of a common Lyapunov functon 15. Although progress has been recently made 5, 16, 19, 25, fndng a common Lyapunov functon s stll an open problem. Qute often, a lnear matrx nequalty (LMI) problem formulaton s used to obtan suffcent stablty condtons by constructng a set of quadratc Lyapunov-lke functons 12, 23. More recently, the swtched Lyapunov functon (SLF) method and less conservatve LMI based condtons were developed n 6 for stablty analyss and control desgn for swtched lnear systems. The partal support of the Natonal Scence Foundaton (NSF CCR ) s gratefully acknowledged. L. Fang, H. Ln and P. Antsakls are wth Department of Electrcal Engneerng, Unversty of Notre Dame, Notre Dame, IN E-mal: {lfang, hln1, antsakls.1}@nd.edu /04/$ IEEE 3265 For swtched lnear systems, the exstence of SLF s a weaker condton than the solvablty of Le algebra, whch mples the exstence of a common quadratc Lyapunov functon 1. The Le algebra approach, however, can be generalzed to study swtched nonlnear systems 17 whle the extenson of SLF approach to the nonlnear context s not straghtforward. In the SLF method, we also assume that the SLF strctly decrease along the solutons of the systems for all tme nstances. Ths restrctve assumpton can be relaxed for certan classes of swtched lnear systems. One can deduce asymptotc stablty usng multple Lyapunov functons whose Le dervatve are only negatve semdefnte 10. The SLF method has been appled to solve dfferent problems of swtched lnear systems 7, 14, 28. In partcular, the nput-output performance problem was studed n 7. There are some related works n the lterature on analyzng the nput-output propertes of swtched systems. For example, the L 2 analyss 29 and the l dsturbance attenuaton problem 13 have been studed. Motvated by the work n 6, we combne the SLF method wth Fnsler s Lemma 21, 26 to study stablty and control desgn problems wth performance analyss for dscrete-tme arbtrarly swtchng lnear systems. Two performance measures consdered n ths paper are the decay rate and the nput-output performance. Frst, a new LMIbased necessary and suffcent condton, whch generalzes the one n 6, s obtaned to check the exstence of a SLF. Our new method s conceptually smple. Here, dfference equatons are consdered as constrants and these dynamcal constrants are ncorporated nto the stablty analyss condton through the use of matrx Lagrange multplers. The key dea s to ncrease the dmenson of the LMIs and to ntroduce new matrx varables, allowng extra degree of freedom for stablty analyss and control desgn. We show that f a SLF exsts, the swtched system s not only asymptotcally stable but also exponentally stable. Of partcular nterest s the formulaton of an LMI optmzaton problem to fnd a sharp estmate of the decay rate to the orgn for swtched lnear systems. To attan a bound on the decay rate, swtched state feedback control desgn s nvestgated. By swtched control desgn, we mean the desgn of state feedback control gans for each subsystem such that the closed-loop swtched system s asymptotcally (or exponentally) stable. The closed-loop swtched system s the one correspondng to the closed-loop subsystems under an arbtrary swtchng rule. Fnally, we study the robustness of the feedback control law when the system s subject to nput dsturbances from a γ-performance pont of vew, whch provdes an upper bound on the

2 worst case energy ampltude gan for swtched systems over all possble nputs and swtchng sgnals. Usng the same mathematcal tool and the well-known S-procedure 2, a new LMI-based suffcent condton s obtaned to ensure the asymptotc stablty of the swtched system whle satsfyng a γ-performance condton. Better performance level than the one n 7 s guaranteed due to extra matrx varables ntroduced n the new LMI condton. Most of the ntroductory materal on SLF can be found n 6. The paper s organzed as follows. Secton II gves the problem formulaton. Stablty analyss of swtched systems under arbtrary swtchng s addressed n Secton III usng the SLF method combned wth Fnsler s Lemma. Exponental stablty of swtched systems s studed n Secton IV. A bound on the decay rate to the orgn s found by solvng an LMI optmzaton problem. Swtched state feedback control gans are then desgned to attan the bound. In Secton V, the nput-output performance and synthess problems are nvestgated. Secton VI concludes the paper. Notaton: The notaton s standard. Z s the set of nteger numbers. R m n denotes the set of m n real matrces. S n denotes the set of n n real symmetrc matrces and S n +, the set of n n real symmetrc postve-defnte matrces. λ(m) stands for the egenvalue of a matrx M. M T s the transpose of the matrx M. M>0(M<0) means that M s postve defnte (negatve defnte). II. PROBLEM FORMULATION Consder the class of swtched lnear systems =A σ(t) +B σ(t) u(t), t Z + (1) where R n s the state, u(t) R m s the control nput, and Z + = {0, 1, 2,...}. The swtchng sgnal σ(t) :Z + I = {1,...,N} s a pecewse constant functon of tme t, whch s unknown a pror. (For notatonal smplcty, we may not explctly menton the tme-dependence of the swtchng sgnal below.) Here, (A,B ) ( I) are constant matrces of approprate dmensons denotng the subsystems, and N 2 s the number of subsystems. We are nterested n the followng mportant questons. Q1) Suffcent condtons for asymptotc stablty: Is t possble to specfy condtons such that the autonomous swtched system s asymptotcally stable? Q2) Decay rate: If the answer to Q1) s yes, can a bound on the decay rate to the orgn be found? Q3) Control desgn: If the answer to Q2) s yes, does a control u(t) exst such that the bound s actually attaned? Q4) Robustness of the control law: Is the control robust when subject to nput dsturbances? In the sequel we show that answers to all the above questons are ndeed yes. III. STABILITY ANALYSIS In ths secton, we nvestgate the stablty of the orgn of an autonomous swtched system. The asymptotc stablty of 3266 the system s verfed by means of a set of LMIs formulated n terms of subsystem A-matrces. If feasble, these LMIs provde a set of Lyapunov matrces that can be combned to form a swtched quadratc Lyapunov functon. The autonomous swtched system s gven by =A σ(t). (2) Defne the ndcaton functon ξ(t) =ξ 1 (t),...,ξ N (t) T (3) wth { 1, σ(t) = ξ (t) = 0, otherwse Then, the swtched system (2) can also be wrtten as N = ξ (t)a. (4) =1 To check asymptotc stablty of system (2), a SLF wth a structure smlar to that of the system descrpton was used 6: ( N ) V (t, ) = x T (t) ξ (t)p, (5) =1 where P S n +,( =1,...,N). If such a postve-defnte Lyapunov functon exsts and V (t, ) = V (t +1,) V (t, ) (6) s negatve defnte along the solutons of (2), then the orgn of the swtched system (2) s asymptotcally stable. In the followng, a new LMI-based necessary and suffcent condton s obtaned by combnng SLF method wth Fnsler s Lemma 21, 26. Ths condton s more general than those condtons of Theorem 2 n 6 and no matrx nverson s nvolved n the constructon of SLF. We frst ntroduce Fnsler s Lemma, whch has been prevously used n the control lterature manly wth the purpose of elmnatng desgn varables n matrx nequaltes. In ths context, Fnsler s Lemma s usually referred to as Elmnaton Lemma. Lemma 1 (Fnsler s Lemma): Let x R n, P S n, and H R m n such that rank(h) =r<n. The followng statements are equvalent: 1) x T Px < 0, Hx =0,x 0. 2) X R n m : P + XH + H T X T < 0. In Lemma 1, tem 1) has a constraned quadratc form n R n whle tem 2) provdes an unconstraned quadratc form, where the constrant s taken nto account by ntroducng multpler X. Recall that the requrement V (t, ) < 0, 0 can be stated as P,P j S n + such that T T P 0 0 P j A I =0, < 0, 0,

3 assumng that σ(t) = and σ(t +1)=j. It now becomes clear that Lemma 1 can be appled to the swtched system under study. In the new procedure, the dynamc dfference equatons that characterze the system are seen as constrants, whch are naturally ncorporated nto the stablty condton usng Fnsler s Lemma. In contrast wth standard space methods, where stablty s carred n the space of the state vector, the stablty test s generated n the enlarged space contanng both the state and ts tme dfference. Theorem 1: There exsts a Lyapunov functon of the form (5) whose dfference s negatve defnte, provng asymptotc stablty of (2) f and only f there exst P S n +, and matrces F,G R n n ( =1,...,N), satsfyng A F T + F A T P A G F G T AT F T P j G G T < 0, (7) (, j) I I. The Lyapunov functon s then gven by (5). Proof: Apply Lemma 1 wth P 0 x,p, 0 P j A H T T F,X I G T, P S n +,F,G R n n ( I). The equvalence between asymptotc stablty of (2) and the followng feasblty test s then establshed: P S n +,F,G R n n ( I): A T F T + F A P A T G F G T A F T P j G G T < 0. The result follows by transposng A n (8). The key dea behnd Theorem 1 s to ncrease the dmenson of the LMIs and to ntroduce new matrx varables F and G, here dentfed as Lagrange multplers, allowng some degree of freedom to verfy (5) and (6). Wth specal choces of F and G, we have the followng corollary. Corollary 1: The followng statements are equvalent. ) There exsts a Lyapunov functon of the form (5) whose dfference s negatve defnte, provng asymptotc stablty of (2). ) There exst P S n +,( =1,...,N), satsfyng P A T P j > 0, (, j) I I. (9) P j A P j The Lyapunov functon s then gven by (5). ) There exst P S n + and G R n n ( =1,...,N), satsfyng P A G G T AT P j G G T < 0, (, j) I I. (10) The Lyapunov functon s gven by (5). (8) 3267 TABLE I NUMERICAL COMPLEXITY ASSOCIATED WITH THREE STABILITY TESTS Stablty Tests K (scalar varables) L (LMI rows) Nn(n+1) Cor. 1-) 2 Nn+2N 2 n Nn(n+1) Cor. 1-) + Nn 2 2 Nn+2N 2 n Nn(n+1) Theorem 1 +2Nn 2 2 Nn+2N 2 n Proof: The equvalence of )-) follows from Theorem 1 by makng F = G =0for tem ) and F =0for tem ), respectvely n condton (8). It s not dffcult to see that Corollary 1 s essentally equvalent to Theorem 2 n 6 whle the proof here s more straghtforward. Moreover, wth no matrx nverson nvolved n the Lyapunov functon, Theorem 1 allows us to formulate an LMI optmzaton problem to fnd a sharp estmate of exponental convergence rate of (2) as llustrated n Secton IV, provded that (7) s feasble. Remark 1: The numercal complexty assocated wth the LMI condtons can be computed n terms of the number K of scalar varables and number L of LMI rows. As dscussed n 2, the number of floatng pont operaton or the tme requred to test the feasblty of the set of LMIs s proportonal to K 3 L. Table I shows K and L as a functon of n (states) and N (subsystems) for three tests presented here. (For a practcal purpose, only those nstances wth N 10, n 10 s tractable.) In the case of restrctve swtchng sgnals, we can modfy these condtons or nvoke the S-procedure to mprove the conservatsm of these condtons. Take for nstance, a system whch does not allow arbtrary transtons between subsystems wll have the set of all ordered pars (, j) of subsystem ndces much smaller than I I. IV. ATTAINABLE BOUNDS ON DECAY RATE A. A Bound on Decay Rate From Theorem 1, the asymptotc stablty of a swtched lnear system under arbtrary swtchng can be checked wth the feasblty test (7). Indeed, we can say more: f the swtched system (2) s asymptotcally stable, t s also exponentally stable about the orgn,.e., κ > 0 and 0 ξ<1, such that κ ξ t x(0). (11) for all ntal condtons x(0) and for all t 0. In fact, the Lyapunov functon (5) s postve defnte, decrescent, and radally unbounded snce V (t, 0) = 0, t 0, and η 2 V (t, ) ρ 2 (12) for all R n wth η = mn I λ mn (P ) and ρ = max I λ max (P ) postve scalars. Furthermore, V ν 2 (13)

4 wth ν =mn (,j) I I λ mn (P A T P ja ) <ρ. Wth the above observatons, the exponental stablty of (2) about the orgn mmedately follows from the well-known fact 24: Theorem 2: The sequence s exponentally stable about the orgn f there exsts a Lyapunov functon V (t, ) such that t 0 η 2 V (t, ) ρ 2, V (t +1,) V (t, ) ν 2, (14) η, ρ, ν > 0. Then κ ξ t x(0), where κ 2 = ρ/η and ξ 2 =1 ν/ρ. We obtan the followng corollary. Corollary 2: If the LMIs gven n Theorem 1 or Corollary 1 are satsfed, the autonomous swtched system (2) s exponentally stable about the orgn. In the case of verfyng exponental stablty of (2), t may be desrable not only fnd feasble solutons to (10) or (7) but to search for solutons that gve an estmate of the decay rate ξ n Theorem 2. To ths end, we want to maxmze the rato ν/ρ. Snce V = x T (t) ( A T P ja P ), constrants (12) and (13) can be rewrtten as: ηi < P <ρi, A F T + F A T P + νi A G F G T AT F T P j G G T (, j) I I, < 0, Note that nequaltes (15) represent strct LMIs but the constrants (12) and (13) are non-strct. Recall that mnmzaton under non-strct LMI constrants gves the same result as mnmzaton under strct LMI constrants when both strct and non-strct LMI constrants are feasble 2. Ths s the case for (15). To obtan a well-posed optmzaton problem, we should normalze ρ to 1 snce the rato ν/ρ can be made arbtrarly small by choosng suffcently large ρ wthout volatng constrants (15). Wth 1 ρ, ν ν/ρ, η η/ρ, G G /ρ, P P /ρ, the followng optmzaton problem s proposed: maxmze ν subject to ηi < P <I, AF T + F A T P + νi A G F G T A T F T P j G G T ν, η > 0, (, j) I I, < 0, (15) (16) Remark 2: Constrants ηi < P < I ( I) lmt the condton number of P to 1/η. The advantage of the condton number lmt s that t wll prevent the LMI soluton algorthm from convergng to P that could lead to roundoff problems. The followng theorem answers the Q2): 3268 Theorem 3: Swtched system (2) s exponentally stable wth a decay rate ξ = (1 ν) 1 2 f the problem (16) s feasble. The assured bound on the decay rate s gven by (1 ν opt ) 1 2 wth νopt the optmal value of optmzaton problem (16). Remark 3: It s possble to produce a better estmate of the decay rate of the swtched system. However, a fner parttonng s needed 27. A SLF havng the same swtchng sgnals as the swtched system may not be suffcent. B. Swtched State Feedback An mportant aspect of the new condtons (7) gven n Theorem 1 s that they are LMIs n P, F and G where there s no cross product between the matrx A and the Lyapunov matrx P ( I). Ths fact has an mpact on the synthess problem consdered below. Let us consder the swtched systems =A σ +B σ u(t) (17) where u(t) s the control and the swtchng sgnal s avalable n real-tme. The stablzng state feedback control problem s to fnd u(t) =K σ (18) such that the correspondng closed-loop swtched system =(A σ + B σ K σ ) (19) s stable. The followng theorem gves a suffcent condton to buld a swtched state feedback controller, whch ensures the exponental stablty of the closed-loop swtched system. Moreover, ths controller s optmal n the sense that the bound on the decay rate of the system s attaned. Q3) s thus solved. Theorem 4: If there s a soluton to maxmze ν subject to ηi < P <I, P + νi A G + B R G T AT + R T BT P j G G T < 0, ν, η > 0, (, j) I I, then the state feedback control gven by (18) wth K = R G 1, I (20) exponentally stablzes the system (17). The decay rate of the system s gven by ξ =(1 ν) 1 2. Proof: Condtons of Theorem 4 lead to ηi < P <I, P + νi (A + B K )G G T (A + B K ) T P j G G T < 0, ν, η > 0, (, j) I I, whch are equvalent to condton (16) wrtten for the closed-loop system (19) wth F = 0. The result then

5 follows from Theorem 3. Note that satsfyng (20) mples P j G G T < 0 and matrces G are non-sngular. Hence the followng feedback gan K = R G 1 s always avalable whenever (20) are feasble. Remark 4: The condton gven n Theorem 4 can also be adapted to the swtched statc output feedback control as shown n 6. V. INPUT-OUTPUT PERFORMANCE The next step s naturally to consder the robustness of the swtched state feedback control. Ths s the last queston, Q4) proposed n Secton II. We study the robustness from a γ-performance pont of vew, that s the desgned swtched feedback control ensures that the worst case energy ampltude gan of the closed loop system s less than or equal to some specfed postve level γ. A. γ-performance Consder an autonomous dscrete-tme swtched system gven by { = Aσ +Bσ w w(t) z(t) = Cσ z+dzw σ w(t) (21) where R n s the state, x(0) = 0, w(t) R q s the dsturbance and z(t) R p s an output vector. The swtchng rule σ s defned as prevously. Smlarly, the matrces (A σ,bσ w,cσ,d z σ zw ) are allowed to take values, at an arbtrary tme, n the fnte set {(A 1,B1 w,c1 z,d1 zw ),...,(A N,BN w,cn z,dn zw )}. Gven γ>0, the γ-performance for the swtched system (21) s defned as below. Defnton 1 (7): The autonomous system (21) s sad to have a γ-performance f t s asymptotcally stable and z T (t)z(t) <γ 2 t=0 w(t) L 2,.e., t=0 w T (t)w(t), (22) w T (t)w(t) <. t=0 Now, defne the same Lyapunov functon V (t, ) > 0, x 0consdered n Secton III, and the modfed Lyapunov stablty condtons V (t, ) < 0, γ 2 w T (t)w(t) z T (t)z(t), (,,w(t),z(t)) satsfyng (21), (,,w(t),z(t)) 0, (23) for a gven γ. The S-procedure 2 s nvoked to generate the equvalent condton V (t, ) <γ 2 w T (t)w(t) z T (t)z(t), (,,w(t),z(t)) satsfyng (21), (,,w(t),z(t)) 0. (24) If (24) s feasble for some 0 <γ< then t s possble to conclude that the system (21) s nternally asymptotcally 3269 stable and has a γ-performance snce 0 < V(t +1,) < γ 2 w T (k)w(k) t t z T (k)z(k), whch s vald for all t>0. In partcular, take t, 0 <V( ) <γ 2 w T (k)w(k) z T (k)z(k), (25) whch mples γ 2 wt (k)w(k) > zt (k)z(k). We are now ready to state the followng theorem whch gves a suffcent condton to check f the autonomous system (21) has a γ-performance. Theorem 5: The system (21) has a γ-performance f there exst P S n +, F 1, G 1 R n n, F 2,G 2 R n p, H 1 R p n, J 1 R q n, H 2 R p p, J 2 R q p, = {1,...,N}, such that where P = U = x P + U + U T < 0, (26) P P j I γ 2 I, (, j) (I I), (27) F 1A + F 2C z F 1 F 2 F 1B w + F 2D zw G 1A + G 2C z G 1 G 2 G 1B w + G 2D zw H 1A + H 2C z H 1 H 2 H 1B w + H 2D zw J 1A + J 2C z J 1 J 2 J 1B w + J 2D zw Proof: Assgn H T z(t) w(t) A T,P C T I 0 0 I B T D T,X and apply Lemma 1 on (24). P P j I γ 2 I F 1 G 1 H 1 J 1 F 2 G 2 H 2 J 2,, B. Swtched Feedback Control wth γ-performance In ths secton, we consder swtched lnear systems gven by { = Aσ +B σ u(t)+bσ w w(t) z(t) = Cσ z+dz σ u(t)+dzw σ w(t) (28) A γ-performance state feedback problem s to fnd a control (18) such that the correspondng closed-loop swtched system { = (Aσ + B σ K σ )+Bσ ww(t) z(t) = (Cσ z + DσK z σ )+Dσ zw (29) w(t) has a γ-performance. The followng theorem gves a suffcent condton to buld such a controller..

6 Theorem 6: There exsts a swtched state feedback control (18) such that the closed-loop swtched system (29) swtched systems, Proc. IEEE 41st CDC, Las Vegas, NE, 2002, pp. 7 J. Daafouz and J. Bernussou, Robust dynamc output feedback for has a γ-performance f there exst P S n +, G 1 R n n , F 2,G 2 R m n, H 2 R m m, J 2 R m q 8 W. P. Dayawansa and C. F. Martn, A converse Lyapunov theorem for, R a class of dynamcal systems whch undergo swtchng, IEEE Trans. R m n, = {1,...,N}, such that AC, vol. 44, no. 4, pp , R. A. Decarlo, M. S. Brancky, S. Pettersson, and B. Lennartson, P + U + U T < 0, (30) Perspectves and results on the stablty and stablzablty of hybrd systems, Proc. of the IEEE, vol. 88, no. 7, pp , where P s gven by (27) and 10 J. Hespanha, Extendng LaSalle s Invarance Prncple to Swtched B wf 2 A G 1 + B R + B wg 2 B wh 2 B wj Lnear Systems, Proc. of the 40th CDC, Orlando, FL, 2001, pp G U = J. P. Hespanha and A. S. Morse, Swtchng between stablzng F 2 G 2 H 2 J. 2 D zw F 2 C G 1 + D R + D zw G 2 D zw H 2 D zw J 2 The γ-performance state feedback control law s gven by (18) wth K = R G 1 1. Proof: Use a transposed verson of Theorem 5 wth F 1 =0,H 1 =0,J 1 =0, (31) and R = K G 1. In order to obtan a convex condton (30), we made the choce n (31) but t s not unque. There are other choces that can also lead to a convex condton. Remark 5: By settng F 2 =0, G 2 =0, H 2 = I, and J 2 = 0, we can recover the correspondng condton n 7. Wthout restrctons on F 2, G 2, H 2, and J 2, a better γ-performance level wth respect to the one n 7 can be obtaned. However, ths came at expense of a more ntensve computaton; see also 22. Remark 6: The best upper bound on the L 2 -nduced gan can be acheved by mnmzng γ subject to the constrants defned by (30), whch s a classcal egenvalue problem 2. VI. CONCLUSION In ths paper the SLF method has been combned wth Fnsler s Lemma to study swtched lnear systems. New and less conservatve LMI condtons are developed for stablty and control desgn problems wth performance analyss. Output feedback control problem can also be treated n a smlar way by usng the technque developed n ths paper. ACKNOWLEDGMENT The authors would lke to thank Prof. Leond Faybusovch at Unversty of Notre Dame for helpful dscussons. REFERENCES 1 A. Agrachev and Danel Lberzon, Le Algebra Condtons for Exponental Stablty of Swtched Systems, Proc. of the 38th CDC, Phoenx, AZ, 1999, pp S. Boyd, L. EL Ghaou, E. Feron, and V. Balakrshnan, Lnear Matrx Inequaltes n System and Control Theory, SIAM, M. S. Brancky, Multple Lyapunov functons and other analyss tools for swtched and hybrd systems, IEEE Trans. Automat. Contr., vol. 43, no. 4, pp , R. W. Brockett, Asymptotc stablty and feedback stablzaton, n Dfferental Geometrc Control Theory, R. W. Brockett, R. S. Mllman, and H. J. Sussmann, Eds. Boston, MA: Brkhuser, 1983, pp D. Cheng, L. Guo, and J. Huang, On quadratc Lyapunov functons, IEEE Trans. AC, vol. 48, no. 5, pp , J. Daafouz, P. Rednger, and C. Iung, Stablty analyss and control synthess for swtched systems: A swtched Lyapunov functon approach, IEEE Trans. AC, vol. 47, no. 11, pp , controllers, Automatca, vol. 38, no. 11, pp , M. Johansson and A. Rantzer, Computaton of pecewse Lyapunov functons for hybrd systems, IEEE Trans. AC, vol. 43, no. 4, pp , H. Ln and P. J. Antsakls, Dsturbance attenuaton propertes for dscrete-tme uncertan lnear swtched systems, n Proc. 42nd IEEE CDC, Mau, HI, 2003, pp Z. J, L. Wang, and G. Xe, Stablzng dscrete-tme swtched systems va observer-based statc output feedback, IEEE Int. Conf. SMC, 2003 Washngton, D.C., pp D. Lberzon and A. S. Morse, Basc problems n stablty and desgn of swtched systems, IEEE Contr. Syst. Mag., vol. 19, no. 5, pp , D. Lberzon and R. Tempo, Gradent algorthm for fndng common Lyapunov functons, Proc. 42nd IEEE CDC, Mau, HI, 2003, pp M. Margalot and D. Lberzon, Le-algebrac stablty condtons for nonlnear swtched systems and dfferental nclusons, submtted to Systems and Control Letters, Jan A. S. Morse, Supervsory control of famles of lnear set-pont controllers - Part 1: Exact Matchng, IEEE Trans. AC, vol. 41, no. 10, pp , K. S. Narendra and J. Balakrshnan, A common Lyapunov functon for stable LTI systems wth commutng A-matrces, IEEE Trans. AC, vol. 39, no. 12, pp , K. S. Narendra and J. Balakrshnan, Adaptve control usng multple models, IEEE Trans. Automat. Contr., vol. 42, pp , M. C. de Olvera and R. E. Skelton, Stablty tests for constraned lnear systems, n S. O. Reza Moheman, edtor, Perspectve n Robust Control, Lecture Notes n Control and Informaton Scence, pp , Sprnger-Verlag, P. J. de Olvera, R. C. L. F. Olvera, V. J. S. Lete, V. F. Montagner, and P. L. D. Peres, LMI based robust stablty condtons for lnear uncertan systems: A numercal comparson, Proc. 41st IEEE CDC, Las Vegas, 2002, pp S. Pettersson and B. Lnnartson, Hybrd system stablty and robustness verfcaton usng lnear matrx nequaltes, Int. J. Control, vol. 75, nos. 16/17, pp , W. J. Rugh, Lnear Systems Theory, Englewood Clffs, New Jersey, Prentce-Hall, R. N. Shorten and K. S. Narendra, On the stablty and exstence of common Lyapunov functons for stable lnear swtchng systems, Proc. 37th IEEE CDC, Tampa, FL, 1998, pp R. E. Skelton, T. Iwasak, and K. M. Grgorads, A Unfed Algebrac Approach to Lnear Control Desgn, London, UK: Taylor & Francs, P. Spångéus, A negatve result on pecewse quadratc Lyapunov functons for decay rate analyss, Proc. 38th IEEE CDC, Phoenx, AZ, 1999, pp D. Xe, L. Wang, F. Hao, and G. Xe, Robust stablty analyss and control synthess for dscrete-tme uncertan swtched systems, Proc. 42nd IEEE CDC, Mau, HI, 2003, pp G. Zha, B. Hu, K. Yasuda, and A. N. Mchel, Qualtatve analyss of dscrete-tme swtched systems, n Proc. ACC, Anchorage, AK, 2002, pp

Static Output Feedback Control for Discrete-time Switched Linear Systems under Arbitrary Switching

Static Output Feedback Control for Discrete-time Switched Linear Systems under Arbitrary Switching 29 Amercan Control Conference Hyatt Regency Rverfront, St. Lous, MO, USA June -2, 29 ThA2.6 Statc Output Feedback Control for Dscrete-tme Swtched Lnear Systems under Arbtrary Swtchng Da-We Dng and Guang-Hong

More information

Stability analysis for class of switched nonlinear systems

Stability analysis for class of switched nonlinear systems Stablty analyss for class of swtched nonlnear systems The MIT Faculty has made ths artcle openly avalable. Please share how ths access benefts you. Your story matters. Ctaton As Publshed Publsher Shaker,

More information

APPENDIX A Some Linear Algebra

APPENDIX A Some Linear Algebra APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,

More information

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017 U.C. Berkeley CS94: Beyond Worst-Case Analyss Handout 4s Luca Trevsan September 5, 07 Summary of Lecture 4 In whch we ntroduce semdefnte programmng and apply t to Max Cut. Semdefnte Programmng Recall that

More information

princeton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg

princeton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg prnceton unv. F 17 cos 521: Advanced Algorthm Desgn Lecture 7: LP Dualty Lecturer: Matt Wenberg Scrbe: LP Dualty s an extremely useful tool for analyzng structural propertes of lnear programs. Whle there

More information

MMA and GCMMA two methods for nonlinear optimization

MMA and GCMMA two methods for nonlinear optimization MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could

More information

Control of Uncertain Bilinear Systems using Linear Controllers: Stability Region Estimation and Controller Design

Control of Uncertain Bilinear Systems using Linear Controllers: Stability Region Estimation and Controller Design Control of Uncertan Blnear Systems usng Lnear Controllers: Stablty Regon Estmaton Controller Desgn Shoudong Huang Department of Engneerng Australan Natonal Unversty Canberra, ACT 2, Australa shoudong.huang@anu.edu.au

More information

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:

More information

THE GUARANTEED COST CONTROL FOR UNCERTAIN LARGE SCALE INTERCONNECTED SYSTEMS

THE GUARANTEED COST CONTROL FOR UNCERTAIN LARGE SCALE INTERCONNECTED SYSTEMS Copyrght 22 IFAC 5th rennal World Congress, Barcelona, Span HE GUARANEED COS CONROL FOR UNCERAIN LARGE SCALE INERCONNECED SYSEMS Hroak Mukadan Yasuyuk akato Yoshyuk anaka Koch Mzukam Faculty of Informaton

More information

Introduction. - The Second Lyapunov Method. - The First Lyapunov Method

Introduction. - The Second Lyapunov Method. - The First Lyapunov Method Stablty Analyss A. Khak Sedgh Control Systems Group Faculty of Electrcal and Computer Engneerng K. N. Toos Unversty of Technology February 2009 1 Introducton Stablty s the most promnent characterstc of

More information

The Order Relation and Trace Inequalities for. Hermitian Operators

The Order Relation and Trace Inequalities for. Hermitian Operators Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence

More information

Solutions to exam in SF1811 Optimization, Jan 14, 2015

Solutions to exam in SF1811 Optimization, Jan 14, 2015 Solutons to exam n SF8 Optmzaton, Jan 4, 25 3 3 O------O -4 \ / \ / The network: \/ where all lnks go from left to rght. /\ / \ / \ 6 O------O -5 2 4.(a) Let x = ( x 3, x 4, x 23, x 24 ) T, where the varable

More information

Lecture 10 Support Vector Machines II

Lecture 10 Support Vector Machines II Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed

More information

More metrics on cartesian products

More metrics on cartesian products More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of

More information

Solutions HW #2. minimize. Ax = b. Give the dual problem, and make the implicit equality constraints explicit. Solution.

Solutions HW #2. minimize. Ax = b. Give the dual problem, and make the implicit equality constraints explicit. Solution. Solutons HW #2 Dual of general LP. Fnd the dual functon of the LP mnmze subject to c T x Gx h Ax = b. Gve the dual problem, and make the mplct equalty constrants explct. Soluton. 1. The Lagrangan s L(x,

More information

COS 521: Advanced Algorithms Game Theory and Linear Programming

COS 521: Advanced Algorithms Game Theory and Linear Programming COS 521: Advanced Algorthms Game Theory and Lnear Programmng Moses Charkar February 27, 2013 In these notes, we ntroduce some basc concepts n game theory and lnear programmng (LP). We show a connecton

More information

Feature Selection: Part 1

Feature Selection: Part 1 CSE 546: Machne Learnng Lecture 5 Feature Selecton: Part 1 Instructor: Sham Kakade 1 Regresson n the hgh dmensonal settng How do we learn when the number of features d s greater than the sample sze n?

More information

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009 College of Computer & Informaton Scence Fall 2009 Northeastern Unversty 20 October 2009 CS7880: Algorthmc Power Tools Scrbe: Jan Wen and Laura Poplawsk Lecture Outlne: Prmal-dual schema Network Desgn:

More information

Stability and Stabilization for Discrete Systems with Time-varying Delays Based on the Average Dwell-time Method

Stability and Stabilization for Discrete Systems with Time-varying Delays Based on the Average Dwell-time Method Proceedngs of the 29 IEEE Internatonal Conference on Systems, an, and Cybernetcs San Antono, TX, USA - October 29 Stablty and Stablzaton for Dscrete Systems wth Tme-varyng Delays Based on the Average Dwell-tme

More information

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 )

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 ) Kangweon-Kyungk Math. Jour. 4 1996), No. 1, pp. 7 16 AN ITERATIVE ROW-ACTION METHOD FOR MULTICOMMODITY TRANSPORTATION PROBLEMS Yong Joon Ryang Abstract. The optmzaton problems wth quadratc constrants often

More information

Assortment Optimization under MNL

Assortment Optimization under MNL Assortment Optmzaton under MNL Haotan Song Aprl 30, 2017 1 Introducton The assortment optmzaton problem ams to fnd the revenue-maxmzng assortment of products to offer when the prces of products are fxed.

More information

Inexact Newton Methods for Inverse Eigenvalue Problems

Inexact Newton Methods for Inverse Eigenvalue Problems Inexact Newton Methods for Inverse Egenvalue Problems Zheng-jan Ba Abstract In ths paper, we survey some of the latest development n usng nexact Newton-lke methods for solvng nverse egenvalue problems.

More information

Difference Equations

Difference Equations Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information

The lower and upper bounds on Perron root of nonnegative irreducible matrices

The lower and upper bounds on Perron root of nonnegative irreducible matrices Journal of Computatonal Appled Mathematcs 217 (2008) 259 267 wwwelsevercom/locate/cam The lower upper bounds on Perron root of nonnegatve rreducble matrces Guang-Xn Huang a,, Feng Yn b,keguo a a College

More information

Errors for Linear Systems

Errors for Linear Systems Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch

More information

Computing Correlated Equilibria in Multi-Player Games

Computing Correlated Equilibria in Multi-Player Games Computng Correlated Equlbra n Mult-Player Games Chrstos H. Papadmtrou Presented by Zhanxang Huang December 7th, 2005 1 The Author Dr. Chrstos H. Papadmtrou CS professor at UC Berkley (taught at Harvard,

More information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and

More information

TAIL BOUNDS FOR SUMS OF GEOMETRIC AND EXPONENTIAL VARIABLES

TAIL BOUNDS FOR SUMS OF GEOMETRIC AND EXPONENTIAL VARIABLES TAIL BOUNDS FOR SUMS OF GEOMETRIC AND EXPONENTIAL VARIABLES SVANTE JANSON Abstract. We gve explct bounds for the tal probabltes for sums of ndependent geometrc or exponental varables, possbly wth dfferent

More information

6.854J / J Advanced Algorithms Fall 2008

6.854J / J Advanced Algorithms Fall 2008 MIT OpenCourseWare http://ocw.mt.edu 6.854J / 18.415J Advanced Algorthms Fall 2008 For nformaton about ctng these materals or our Terms of Use, vst: http://ocw.mt.edu/terms. 18.415/6.854 Advanced Algorthms

More information

Maximizing the number of nonnegative subsets

Maximizing the number of nonnegative subsets Maxmzng the number of nonnegatve subsets Noga Alon Hao Huang December 1, 213 Abstract Gven a set of n real numbers, f the sum of elements of every subset of sze larger than k s negatve, what s the maxmum

More information

Improved delay-dependent stability criteria for discrete-time stochastic neural networks with time-varying delays

Improved delay-dependent stability criteria for discrete-time stochastic neural networks with time-varying delays Avalable onlne at www.scencedrect.com Proceda Engneerng 5 ( 4456 446 Improved delay-dependent stablty crtera for dscrete-tme stochastc neural networs wth tme-varyng delays Meng-zhuo Luo a Shou-mng Zhong

More information

Perfect Competition and the Nash Bargaining Solution

Perfect Competition and the Nash Bargaining Solution Perfect Competton and the Nash Barganng Soluton Renhard John Department of Economcs Unversty of Bonn Adenauerallee 24-42 53113 Bonn, Germany emal: rohn@un-bonn.de May 2005 Abstract For a lnear exchange

More information

Some modelling aspects for the Matlab implementation of MMA

Some modelling aspects for the Matlab implementation of MMA Some modellng aspects for the Matlab mplementaton of MMA Krster Svanberg krlle@math.kth.se Optmzaton and Systems Theory Department of Mathematcs KTH, SE 10044 Stockholm September 2004 1. Consdered optmzaton

More information

Lecture 21: Numerical methods for pricing American type derivatives

Lecture 21: Numerical methods for pricing American type derivatives Lecture 21: Numercal methods for prcng Amercan type dervatves Xaoguang Wang STAT 598W Aprl 10th, 2014 (STAT 598W) Lecture 21 1 / 26 Outlne 1 Fnte Dfference Method Explct Method Penalty Method (STAT 598W)

More information

Robust observed-state feedback design. for discrete-time systems rational in the uncertainties

Robust observed-state feedback design. for discrete-time systems rational in the uncertainties Robust observed-state feedback desgn for dscrete-tme systems ratonal n the uncertantes Dmtr Peaucelle Yosho Ebhara & Yohe Hosoe Semnar at Kolloquum Technsche Kybernetk, May 10, 016 Unversty of Stuttgart

More information

Stability Analysis and Anti-windup Design of Switched Systems with Actuator Saturation

Stability Analysis and Anti-windup Design of Switched Systems with Actuator Saturation Internatonal Journal of Automaton and Computng 14(5), October 2017, 615-625 DOI: 101007/s11633-015-0920-z Stablty Analyss and Ant-wndup Desgn of Swtched Systems wth Actuator Saturaton Xn-Quan Zhang 1 Xao-Yn

More information

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve

More information

A CHARACTERIZATION OF ADDITIVE DERIVATIONS ON VON NEUMANN ALGEBRAS

A CHARACTERIZATION OF ADDITIVE DERIVATIONS ON VON NEUMANN ALGEBRAS Journal of Mathematcal Scences: Advances and Applcatons Volume 25, 2014, Pages 1-12 A CHARACTERIZATION OF ADDITIVE DERIVATIONS ON VON NEUMANN ALGEBRAS JIA JI, WEN ZHANG and XIAOFEI QI Department of Mathematcs

More information

MATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS

MATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS MATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS These are nformal notes whch cover some of the materal whch s not n the course book. The man purpose s to gve a number of nontrval examples

More information

A new construction of 3-separable matrices via an improved decoding of Macula s construction

A new construction of 3-separable matrices via an improved decoding of Macula s construction Dscrete Optmzaton 5 008 700 704 Contents lsts avalable at ScenceDrect Dscrete Optmzaton journal homepage: wwwelsevercom/locate/dsopt A new constructon of 3-separable matrces va an mproved decodng of Macula

More information

ECE559VV Project Report

ECE559VV Project Report ECE559VV Project Report (Supplementary Notes Loc Xuan Bu I. MAX SUM-RATE SCHEDULING: THE UPLINK CASE We have seen (n the presentaton that, for downlnk (broadcast channels, the strategy maxmzng the sum-rate

More information

The equation of motion of a dynamical system is given by a set of differential equations. That is (1)

The equation of motion of a dynamical system is given by a set of differential equations. That is (1) Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence

More information

Problem Set 9 Solutions

Problem Set 9 Solutions Desgn and Analyss of Algorthms May 4, 2015 Massachusetts Insttute of Technology 6.046J/18.410J Profs. Erk Demane, Srn Devadas, and Nancy Lynch Problem Set 9 Solutons Problem Set 9 Solutons Ths problem

More information

CHAPTER 4. Vector Spaces

CHAPTER 4. Vector Spaces man 2007/2/16 page 234 CHAPTER 4 Vector Spaces To crtcze mathematcs for ts abstracton s to mss the pont entrel. Abstracton s what makes mathematcs work. Ian Stewart The man am of ths tet s to stud lnear

More information

Decentralized robust control design using LMI

Decentralized robust control design using LMI Acta Montanstca Slovaca Ročník 3 (28) číslo -4 Decentralzed robust control desgn usng LMI Anna Flasová and Dušan Krokavec ávrh robustného decentralzovaného radena pomocou LMI he paper deals wth applcaton

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

Blowing Up Affine Hybrid Systems

Blowing Up Affine Hybrid Systems 43rd IEEE Conference on Decson and Control December 14-17, 004 Atlants, Paradse Island, Bahamas TuB01.3 Blowng Up Affne Hybrd Systems Aaron D. Ames* and Shankar Sastry* Department of EECS Unversty of Calforna

More information

FINITE-TIME BOUNDEDNESS AND STABILIZATION OF SWITCHED LINEAR SYSTEMS

FINITE-TIME BOUNDEDNESS AND STABILIZATION OF SWITCHED LINEAR SYSTEMS K Y BERNETIKA VOLUM E 46 21, NUMBER 5, P AGES 87 889 FINITE-TIME BOUNDEDNESS AND STABILIZATION OF SWITCHED LINEAR SYSTEMS Habo Du, Xangze Ln and Shhua L In ths paper, fnte-tme boundedness and stablzaton

More information

EEE 241: Linear Systems

EEE 241: Linear Systems EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they

More information

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,

More information

Module 9. Lecture 6. Duality in Assignment Problems

Module 9. Lecture 6. Duality in Assignment Problems Module 9 1 Lecture 6 Dualty n Assgnment Problems In ths lecture we attempt to answer few other mportant questons posed n earler lecture for (AP) and see how some of them can be explaned through the concept

More information

Outline. Communication. Bellman Ford Algorithm. Bellman Ford Example. Bellman Ford Shortest Path [1]

Outline. Communication. Bellman Ford Algorithm. Bellman Ford Example. Bellman Ford Shortest Path [1] DYNAMIC SHORTEST PATH SEARCH AND SYNCHRONIZED TASK SWITCHING Jay Wagenpfel, Adran Trachte 2 Outlne Shortest Communcaton Path Searchng Bellmann Ford algorthm Algorthm for dynamc case Modfcatons to our algorthm

More information

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction ECONOMICS 5* -- NOTE (Summary) ECON 5* -- NOTE The Multple Classcal Lnear Regresson Model (CLRM): Specfcaton and Assumptons. Introducton CLRM stands for the Classcal Lnear Regresson Model. The CLRM s also

More information

The Minimum Universal Cost Flow in an Infeasible Flow Network

The Minimum Universal Cost Flow in an Infeasible Flow Network Journal of Scences, Islamc Republc of Iran 17(2): 175-180 (2006) Unversty of Tehran, ISSN 1016-1104 http://jscencesutacr The Mnmum Unversal Cost Flow n an Infeasble Flow Network H Saleh Fathabad * M Bagheran

More information

A Local Variational Problem of Second Order for a Class of Optimal Control Problems with Nonsmooth Objective Function

A Local Variational Problem of Second Order for a Class of Optimal Control Problems with Nonsmooth Objective Function A Local Varatonal Problem of Second Order for a Class of Optmal Control Problems wth Nonsmooth Objectve Functon Alexander P. Afanasev Insttute for Informaton Transmsson Problems, Russan Academy of Scences,

More information

Economics 101. Lecture 4 - Equilibrium and Efficiency

Economics 101. Lecture 4 - Equilibrium and Efficiency Economcs 0 Lecture 4 - Equlbrum and Effcency Intro As dscussed n the prevous lecture, we wll now move from an envronment where we looed at consumers mang decsons n solaton to analyzng economes full of

More information

NP-Completeness : Proofs

NP-Completeness : Proofs NP-Completeness : Proofs Proof Methods A method to show a decson problem Π NP-complete s as follows. (1) Show Π NP. (2) Choose an NP-complete problem Π. (3) Show Π Π. A method to show an optmzaton problem

More information

Lecture 17: Lee-Sidford Barrier

Lecture 17: Lee-Sidford Barrier CSE 599: Interplay between Convex Optmzaton and Geometry Wnter 2018 Lecturer: Yn Tat Lee Lecture 17: Lee-Sdford Barrer Dsclamer: Please tell me any mstake you notced. In ths lecture, we talk about the

More information

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0 MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector

More information

Approximate D-optimal designs of experiments on the convex hull of a finite set of information matrices

Approximate D-optimal designs of experiments on the convex hull of a finite set of information matrices Approxmate D-optmal desgns of experments on the convex hull of a fnte set of nformaton matrces Radoslav Harman, Mára Trnovská Department of Appled Mathematcs and Statstcs Faculty of Mathematcs, Physcs

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION

More information

General viscosity iterative method for a sequence of quasi-nonexpansive mappings

General viscosity iterative method for a sequence of quasi-nonexpansive mappings Avalable onlne at www.tjnsa.com J. Nonlnear Sc. Appl. 9 (2016), 5672 5682 Research Artcle General vscosty teratve method for a sequence of quas-nonexpansve mappngs Cuje Zhang, Ynan Wang College of Scence,

More information

Welfare Properties of General Equilibrium. What can be said about optimality properties of resource allocation implied by general equilibrium?

Welfare Properties of General Equilibrium. What can be said about optimality properties of resource allocation implied by general equilibrium? APPLIED WELFARE ECONOMICS AND POLICY ANALYSIS Welfare Propertes of General Equlbrum What can be sad about optmalty propertes of resource allocaton mpled by general equlbrum? Any crteron used to compare

More information

Generalized Linear Methods

Generalized Linear Methods Generalzed Lnear Methods 1 Introducton In the Ensemble Methods the general dea s that usng a combnaton of several weak learner one could make a better learner. More formally, assume that we have a set

More information

Deriving the X-Z Identity from Auxiliary Space Method

Deriving the X-Z Identity from Auxiliary Space Method Dervng the X-Z Identty from Auxlary Space Method Long Chen Department of Mathematcs, Unversty of Calforna at Irvne, Irvne, CA 92697 chenlong@math.uc.edu 1 Iteratve Methods In ths paper we dscuss teratve

More information

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1 P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the

More information

LOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin

LOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin Proceedngs of the 007 Wnter Smulaton Conference S G Henderson, B Bller, M-H Hseh, J Shortle, J D Tew, and R R Barton, eds LOW BIAS INTEGRATED PATH ESTIMATORS James M Calvn Department of Computer Scence

More information

A New Refinement of Jacobi Method for Solution of Linear System Equations AX=b

A New Refinement of Jacobi Method for Solution of Linear System Equations AX=b Int J Contemp Math Scences, Vol 3, 28, no 17, 819-827 A New Refnement of Jacob Method for Soluton of Lnear System Equatons AX=b F Naem Dafchah Department of Mathematcs, Faculty of Scences Unversty of Gulan,

More information

STAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 16

STAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 16 STAT 39: MATHEMATICAL COMPUTATIONS I FALL 218 LECTURE 16 1 why teratve methods f we have a lnear system Ax = b where A s very, very large but s ether sparse or structured (eg, banded, Toepltz, banded plus

More information

Perron Vectors of an Irreducible Nonnegative Interval Matrix

Perron Vectors of an Irreducible Nonnegative Interval Matrix Perron Vectors of an Irreducble Nonnegatve Interval Matrx Jr Rohn August 4 2005 Abstract As s well known an rreducble nonnegatve matrx possesses a unquely determned Perron vector. As the man result of

More information

FAULT TOLERANT CONTROL OF SWITCHED NONLINEAR SYSTEMS WITH TIME DELAY UNDER ASYNCHRONOUS SWITCHING

FAULT TOLERANT CONTROL OF SWITCHED NONLINEAR SYSTEMS WITH TIME DELAY UNDER ASYNCHRONOUS SWITCHING Int. J. Appl. Math. Comput. Sc. 2010 Vol. 20 No. 3 497 506 DOI: 10.2478/v10006-010-0036-0 FAULT TOLERANT CONTROL OF SWITCHED NONLINEAR SYSTEMS WITH TIME DELAY UNDER ASYNCHRONOUS SWITCHING ZHENGRONG XIANG

More information

FREQUENCY DISTRIBUTIONS Page 1 of The idea of a frequency distribution for sets of observations will be introduced,

FREQUENCY DISTRIBUTIONS Page 1 of The idea of a frequency distribution for sets of observations will be introduced, FREQUENCY DISTRIBUTIONS Page 1 of 6 I. Introducton 1. The dea of a frequency dstrbuton for sets of observatons wll be ntroduced, together wth some of the mechancs for constructng dstrbutons of data. Then

More information

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Structure and Drive Paul A. Jensen Copyright July 20, 2003 Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.

More information

Kernel Methods and SVMs Extension

Kernel Methods and SVMs Extension Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general

More information

Application of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems

Application of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems Mathematca Aeterna, Vol. 1, 011, no. 06, 405 415 Applcaton of B-Splne to Numercal Soluton of a System of Sngularly Perturbed Problems Yogesh Gupta Department of Mathematcs Unted College of Engneerng &

More information

U.C. Berkeley CS294: Spectral Methods and Expanders Handout 8 Luca Trevisan February 17, 2016

U.C. Berkeley CS294: Spectral Methods and Expanders Handout 8 Luca Trevisan February 17, 2016 U.C. Berkeley CS94: Spectral Methods and Expanders Handout 8 Luca Trevsan February 7, 06 Lecture 8: Spectral Algorthms Wrap-up In whch we talk about even more generalzatons of Cheeger s nequaltes, and

More information

R R 2 fl 2 + F 2 z j z 2 R n, g Fgure : R s polytopc and the boundary of R and R j s characterzed as beng a subset of fl j + F j z j z 2 R n, g (here

R R 2 fl 2 + F 2 z j z 2 R n, g Fgure : R s polytopc and the boundary of R and R j s characterzed as beng a subset of fl j + F j z j z 2 R n, g (here Quadratc Stablzaton and Control of Pecewse-Lnear Systems Arash Hassb Stephen Boyd Informaton Systems Laboratory Stanford Unversty Stanford, CA 935-55 USA arash@sl.stanford.edu boyd@sl.stanford.edu Abstract

More information

Lecture 10 Support Vector Machines. Oct

Lecture 10 Support Vector Machines. Oct Lecture 10 Support Vector Machnes Oct - 20-2008 Lnear Separators Whch of the lnear separators s optmal? Concept of Margn Recall that n Perceptron, we learned that the convergence rate of the Perceptron

More information

Convexity preserving interpolation by splines of arbitrary degree

Convexity preserving interpolation by splines of arbitrary degree Computer Scence Journal of Moldova, vol.18, no.1(52), 2010 Convexty preservng nterpolaton by splnes of arbtrary degree Igor Verlan Abstract In the present paper an algorthm of C 2 nterpolaton of dscrete

More information

Output Feedback Stabilization of Networked Control Systems With Packet Dropouts

Output Feedback Stabilization of Networked Control Systems With Packet Dropouts IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 52, NO. 9, SEPTEMBER 2007 705 Output Feedback Stablzaton of Networked Control Systems Wth Packet Dropouts Wen-An Zhang and L Yu Abstract In ths paper, we dscuss

More information

Case A. P k = Ni ( 2L i k 1 ) + (# big cells) 10d 2 P k.

Case A. P k = Ni ( 2L i k 1 ) + (# big cells) 10d 2 P k. THE CELLULAR METHOD In ths lecture, we ntroduce the cellular method as an approach to ncdence geometry theorems lke the Szemeréd-Trotter theorem. The method was ntroduced n the paper Combnatoral complexty

More information

FUZZY GOAL PROGRAMMING VS ORDINARY FUZZY PROGRAMMING APPROACH FOR MULTI OBJECTIVE PROGRAMMING PROBLEM

FUZZY GOAL PROGRAMMING VS ORDINARY FUZZY PROGRAMMING APPROACH FOR MULTI OBJECTIVE PROGRAMMING PROBLEM Internatonal Conference on Ceramcs, Bkaner, Inda Internatonal Journal of Modern Physcs: Conference Seres Vol. 22 (2013) 757 761 World Scentfc Publshng Company DOI: 10.1142/S2010194513010982 FUZZY GOAL

More information

Neuro-Adaptive Design - I:

Neuro-Adaptive Design - I: Lecture 36 Neuro-Adaptve Desgn - I: A Robustfyng ool for Dynamc Inverson Desgn Dr. Radhakant Padh Asst. Professor Dept. of Aerospace Engneerng Indan Insttute of Scence - Bangalore Motvaton Perfect system

More information

Additional Codes using Finite Difference Method. 1 HJB Equation for Consumption-Saving Problem Without Uncertainty

Additional Codes using Finite Difference Method. 1 HJB Equation for Consumption-Saving Problem Without Uncertainty Addtonal Codes usng Fnte Dfference Method Benamn Moll 1 HJB Equaton for Consumpton-Savng Problem Wthout Uncertanty Before consderng the case wth stochastc ncome n http://www.prnceton.edu/~moll/ HACTproect/HACT_Numercal_Appendx.pdf,

More information

Lecture 4. Instructor: Haipeng Luo

Lecture 4. Instructor: Haipeng Luo Lecture 4 Instructor: Hapeng Luo In the followng lectures, we focus on the expert problem and study more adaptve algorthms. Although Hedge s proven to be worst-case optmal, one may wonder how well t would

More information

Zeros and Zero Dynamics for Linear, Time-delay System

Zeros and Zero Dynamics for Linear, Time-delay System UNIVERSITA POLITECNICA DELLE MARCHE - FACOLTA DI INGEGNERIA Dpartmento d Ingegnerua Informatca, Gestonale e dell Automazone LabMACS Laboratory of Modelng, Analyss and Control of Dynamcal System Zeros and

More information

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family IOSR Journal of Mathematcs IOSR-JM) ISSN: 2278-5728. Volume 3, Issue 3 Sep-Oct. 202), PP 44-48 www.osrjournals.org Usng T.O.M to Estmate Parameter of dstrbutons that have not Sngle Exponental Famly Jubran

More information

Quantum and Classical Information Theory with Disentropy

Quantum and Classical Information Theory with Disentropy Quantum and Classcal Informaton Theory wth Dsentropy R V Ramos rubensramos@ufcbr Lab of Quantum Informaton Technology, Department of Telenformatc Engneerng Federal Unversty of Ceara - DETI/UFC, CP 6007

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 12 10/21/2013. Martingale Concentration Inequalities and Applications

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 12 10/21/2013. Martingale Concentration Inequalities and Applications MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.65/15.070J Fall 013 Lecture 1 10/1/013 Martngale Concentraton Inequaltes and Applcatons Content. 1. Exponental concentraton for martngales wth bounded ncrements.

More information

Lecture Notes on Linear Regression

Lecture Notes on Linear Regression Lecture Notes on Lnear Regresson Feng L fl@sdueducn Shandong Unversty, Chna Lnear Regresson Problem In regresson problem, we am at predct a contnuous target value gven an nput feature vector We assume

More information

Stability of Switched Linear Systems on Cones: A Generating Function Approach

Stability of Switched Linear Systems on Cones: A Generating Function Approach 49th IEEE Conference on Decson and Control December 15-17, 2010 Hlton Atlanta Hotel, Atlanta, GA, USA Stablty of Swtched Lnear Systems on Cones: A Generatng Functon Approach Jngla Shen and Jangha Hu Abstract

More information

A RELAXED SUFFICIENT CONDITION FOR ROBUST STABILITY OF AFFINE SYSTEMS

A RELAXED SUFFICIENT CONDITION FOR ROBUST STABILITY OF AFFINE SYSTEMS Доклади на Българската академия на науките Comptes rendus de l Académe bulgare des Scences Tome 60 No 9 2007 SCIENCES ET INGENIERIE Théore des systèmes A RELAXED SUFFICIENT CONDITION FOR ROBUST STABILITY

More information

Real-Time Systems. Multiprocessor scheduling. Multiprocessor scheduling. Multiprocessor scheduling

Real-Time Systems. Multiprocessor scheduling. Multiprocessor scheduling. Multiprocessor scheduling Real-Tme Systems Multprocessor schedulng Specfcaton Implementaton Verfcaton Multprocessor schedulng -- -- Global schedulng How are tasks assgned to processors? Statc assgnment The processor(s) used for

More information

BOUNDEDNESS OF THE RIESZ TRANSFORM WITH MATRIX A 2 WEIGHTS

BOUNDEDNESS OF THE RIESZ TRANSFORM WITH MATRIX A 2 WEIGHTS BOUNDEDNESS OF THE IESZ TANSFOM WITH MATIX A WEIGHTS Introducton Let L = L ( n, be the functon space wth norm (ˆ f L = f(x C dx d < For a d d matrx valued functon W : wth W (x postve sem-defnte for all

More information

arxiv: v1 [math.ho] 18 May 2008

arxiv: v1 [math.ho] 18 May 2008 Recurrence Formulas for Fbonacc Sums Adlson J. V. Brandão, João L. Martns 2 arxv:0805.2707v [math.ho] 8 May 2008 Abstract. In ths artcle we present a new recurrence formula for a fnte sum nvolvng the Fbonacc

More information

Calculation of time complexity (3%)

Calculation of time complexity (3%) Problem 1. (30%) Calculaton of tme complexty (3%) Gven n ctes, usng exhaust search to see every result takes O(n!). Calculaton of tme needed to solve the problem (2%) 40 ctes:40! dfferent tours 40 add

More information