Control Using Logic & Switching: Part III Supervisory Control

Size: px
Start display at page:

Download "Control Using Logic & Switching: Part III Supervisory Control"

Transcription

1 Control Using Logic & Switching: Part III Spervisor Control Ttorial for the 40th CDC João P. Hespanha Universit of Sothern California Universit of California at Santa Barbara Otline Spervisor control overview Estimator-based linear spervisor control Estimator-based nonlinear spervisor control Examples 1

2 Spervisor control Motivation: in the control of complex and highl ncertain sstems, traditional methodologies based on a single controller do not provide satisfactor performance. spervisor controller 1 bank of candidate controllers controller n switching control exogenos distrbance/ noise w measred otpt Ke ideas: 1. Bild a bank of alternative controllers 2. Switch among them online based on measrements For simplicit we assme a stabilization problem, otherwise controllers shold have a reference inpt r Spervisor control Motivation: in the control of complex and highl ncertain sstems, traditional methodologies based on a single controller do not provide satisfactor performance. spervisor controller 1 bank of candidate controllers controller n switching control exogenos distrbance/ noise w measred otpt Spervisor: places in the feedback loop the controller that seems more promising based on the available measrements tpicall logic-based/hbrid sstem 2

3 Mlti-controller measred otpt controller 1 bank of candidate controllers controller n switching control Conceptal diagram: not efficient for man controllers & not possible for nstable controllers Mlti-controller measred otpt controller 1 bank of candidate controllers controller n switching control Given a famil of (n-dimensional) candidate controllers switching measred otpt control 3

4 (t) switching = 1 Mlti-controller switching times = 1 = 3 = 2 t Given a famil of (n-dimensional) candidate controllers switching measred otpt control Spervisor spervisor switching measred otpt control Tpicall an hbrid sstem: ϕ continos state δ discrete state continos vector field discrete transition fnction otpt fnction 4

5 Tpes of spervision Pre-roted spervision = 1 = 2 = 3 tr one controllers after another in a pre-defined seqence stop when the performance seems acceptable not effective when the nmber of controllers is large Performance-based spervision (direct) keep controller while observed performance is acceptable when performance of crrent controller becomes nacceptable, switch to controller that leads to best expected performance based on available data Estimator-based spervision (indirect) estimate model from observed data select controller based on crrent estimate Certaint Eqivalence Estimator-based spervision s setp w exogenos distrbance/ noise control measred otpt nmodeled Process is assmed to be in a famil dnamics M p small famil of sstems arond a nominal model N p parametric ncertaint for each in a famil M p, at least one candidate controller C q, q Q provides adeqate performance. controller selection fnction in M p, p P controller C q with q =χ(p) provides adeqate performance 5

6 measred otpt control Estimator-based spervisor decision logic switching Process is assmed to be in famil in M p,p P controller C q, q =χ(p) provides adeqate performance Mlti-estimator p estimate of the otpt that wold be correct if the was N p e p otpt estimation error that wold be small if the was N p Decision logic: e p small likel in M p shold se C q, q = χ(p) set = χ(p) Certaint eqivalence inspired measred otpt control Estimator-based spervisor Process is assmed to be in famil in M p, p P Mlti-estimator p estimate of the otpt that wold be correct if the was N p e p otpt estimation error that wold be small if the was N p Decision logic: mltiestimator mltiestimator decision logic switching A stabilit argment cannot be based on this becase tpicall in M controller p e p small C q, q =χ(p) bt not the converse provides adeqate performance e p small likel in M p shold se C q, q = χ(p) set = χ(p) Certaint eqivalence inspired 6

7 measred otpt control Estimator-based spervisor mltiestimator decision logic switching Process is assmed to be in famil in M p, p P controller C q, q =χ(p) provides adeqate performance Mlti-estimator p estimate of the otpt that wold be correct if the was N p e p otpt estimation error that wold be small if the was N p Decision logic: e p small set overall sstem overall state = χ(p) is detectable is small throgh e p Certaint eqivalence inspired, bt formall jstified b detectabilit Performance-based spervision measred otpt control performance monitor decision logic switching Candidate controllers: Performance monitor: π q measre of the expected performance of controller C q inferred from past data Decision logic: π is acceptable π is nacceptable keep crrent controller switch to controller C q corresponding to best π q 7

8 Abstract spervision Estimator and performance-based architectres share the same common architectre decision logic mlti-est. or perf. monitor switching w mlticontroller control measred otpt In this talk we will focs mostl on an estimator-based spervisor Abstract spervision Estimator and performance-based architectres share the same common architectre decision logic mltiestimator switching w mlticontroller control measred otpt switched sstem 8

9 The for basic properties (1-2) decision logic Matching propert: At least one of the e p is small Wh? in p * P: in M p* e p* is small switching essentiall a reqirement on the mlti-estimator Detectabilit propert: For each p P, the switched sstem is detectable throgh e p when = χ(p) index of controller that stabilizes es in M p essentiall a reqirement on the candidate controllers This propert jstifies sing the candidate controller that corresponds to a small estimation error. Wh? Certaint eqivalence stabilization theorem The for basic properties (3-4) switching decision logic Small error propert: There is a switching ρ : [0, ) P for which e ρ is small compared to an fixed e p and that is consistent with, i.e., = χ(ρ) ρ(t) can be viewed as crrent parameter estimate controller consistent with parameter estimate Non-destabilization propert: Detectabilit is preserved for the time-varing switched sstem (not jst for constant ) Tpicall reqires some form of slow switching Both are essentiall (conflicting) properties of the decision logic 9

10 Analsis otline (linear case, w = 0) switching decision logic 1st b the Matching propert: p * P sch that e p* is small 2nd b the Small error propert: ρ sch that = χ(ρ) and e ρ is small (when compared with e p* ) 3rd b the Detectabilit propert: there exist matrices K p sch that the matrices A q K p C p, q = χ(p) are asmptoticall stable injected 4th the switched sstem can be written as sstem asmptoticall stable b non-destabilization propert small b 2nd step x is small (and converges to zero if, e.g., e ρ L 2 ) Otline Spervisor control overview Estimator-based linear spervisor control Estimator-based nonlinear spervisor control Examples 10

11 Estimator-based linear spervisor control decision logic mltiestimator switching LINEAR w mlticontroller control measred otpt Class of admissible es d exogenos distrbance n measrement noise control measred otpt The transfer fnction from to is an nknown element of parametric ncertaint M p small famil of sstems arond a nominal transfer fnction ν p Tpicall mltiplicative nmod. dnamics additive nmod. dnamics or co-prime factorization of ν p (SISO) 11

12 A word on norms Given a and λ 0 Given a transfer fnction ν and λ 0 e λ t weighted L 2 norm of trncated to [0, t ) e λ t weighted H norm of ν (finite if all poles of ν have real part smaller than λ) The e λ t weighted L 2 indced norm of transfer fnction ν is nmericall eqal to the e λ t weighted H norm of ν ν ν,λ = γ de to initial conditions we will sometimes need a stabilit margin λ > 0... Candidate controllers Class of admissible es M p small famil of sstems arond a nominal transfer fnction ν p Assme given a famil of candidate controller transfer fnctions and a controller selection fnction χ : P Q sch that p P controller κ q, q = χ(p) stabilizes all es in M p χ maps parameter vales with the corresponding stabilizing controller No other constrain is posed on the candidate controllers: κ q can be designed sing an (nonadaptive) techniqe (e.g., pole placement, LQG/LQR, H-infinit, etc.) 12

13 measred otpt control Mlti-estimator mltiestimator with A E asmptoticall stable Matching propert: There exist positive constants c 0, c w, c ε, λ and some p* P sch that noise/distrb. initial cond. e p* is L 2 when noise/distrb. and nmodeled dnamics are L 2 recall nmodeled dnamics A simple mlti-estimator Class of admissible es (2 elements) Assming realizations are detectable, we cold make Lenberger observers asmptoticall stable Mlti-estimator: A E x E D E B E c p* (s I - A p* ) -1 b p* e p* = p* - 0 (exp. fast) with noise & nmodeled dnamics 13

14 Mlti-controller Given a famil of candidate controller transfer fnctions Compte (n C -dimensional) stabilizable and detectable realization Mlti-controller switching measred otpt control detectabilit propert? Switched sstem decision logic mltiestimator w mlticontroller switched sstem The switched sstem can be seen as the interconnection of the with the injected sstem essentiall the mlti-controller & mlti-estimator bt now qite 14

15 Constrcting the injected sstem 1st Take an (arbitrar) switching ρ : [0, ) P. 2nd Define the v ú e ρ = ρ 3rd Replace in the eqations of the mlti-estimator and mlti-controller b ρ v. v injected sstem Constrcting the injected sstem 1st Take an (arbitrar) switching ρ : [0, ) P 2nd Define the v ú e ρ = ρ 3rd Replace in the eqations of the mlti-estimator and mlti-controller b ρ v. mlticontroller mltiestimator ρ v v injected sstem 15

16 The injected sstem mlticontroller mltiestimator ρ v B inspection: for each p P, q Q eigenvales of A p q { sbset of eigenvales of A E } U { poles of the feedback interconnection of ν p with κ q } q =χ( p ) κ q stabilizes ν p A p q is asmptoticall stable If we choose ρ sch that =χ(ρ) then A ρ is alwas stable Switched sstem = injected sstem w Wh? becase v ú e ρ = ρ v injected sstem ρ ρ Using this diagram one can prove the detectabilit propert b inspection 16

17 v injected sstem ρ w Detectabilit propert ρ Sppose ρ = p P w = 0 = χ( p ) Q e p = 0 ( e ρ = v = 0 ) stabilit of injected sstem state & otpts, p of injected sstem 0 inpt & otpt, = p e p 0 state 0 whole state of switched sstem 0 detectabilit of The state of the switched sstem converges to zero along an soltion compatible with zero inpt w & otpt e p detectabilit Detectabilit propert w Detectabilit propert: Given an p P, setting ρ = p P and = χ( p ) Q : 1. The injected sstem is asmptoticall stable 2. The switched sstem is detectable throgh e p v injected sstem ρ ρ Also known as the Certaint Eqivalence Stabilization Theorem Stabilit of the injected sstem is not the onl mechanism to achieve detectabilit: e.g.: injected sstem i/o stable min. phase detectabilit of switched sstem (Certaint Eqivalence Otpt Stabilization Theorem) 17

18 Decision logic switching decision logic switched sstem 1. For bondedness one wants e ρ small for some ρ consistent with (i.e., = χ (ρ)) small error 2. To recover the static detectabilit of the time-varing switched sstem one wants slow switching. non-destabilization These are conflicting reqirements: 1. ρ shold follow smallest e p 2. = χ(ρ) shold not var estimation errors injected sstem Dwell-time switching start monitoring s p P measre of the size of e p over a window of length 1/λ forgetting factor wait τ D seconds Non-destabilizing propert: The minimm interval between consective discontinities of is τ D > 0. 18

19 Small error propert Assme P finite and p * P : (e.g., L 2 noise and no nmodeled dnamics) when we select ρ = p at time t we mst have Two possible cases: 1. Switching will stop in finite time T at some p P: < C* Small error propert Assme P finite and p * P : (e.g., L 2 noise and no nmodeled dnamics) Two possible cases: when we select ρ = p at time t we mst have 2. After some finite time T switching will occr onl among elements of a sbset P* of P, each appearing in ρ infinitel man times: < C* 19

20 Small error propert Assme P finite and p * P : (e.g., L 2 noise and no nmodeled dnamics) when we select ρ = p at time t we mst have Small error propert: (L 2 case) Assme that P is a finite set. If p * P for which then at least one error L 2 switched error will be L 2 Small error propert Small error propert: (L 2 case) Assme that P is a finite set. If p * P for which then at least one error L 2 switched error will be L 2 (for switching ρ defined b the logic) Small error propert: (general case) Assme that Q is a finite set with m element. For ever p P, t 0, switching ρ t : [0, t ) P sch that: 1. = χ(ρ t ) except at most on m time intervals of length τ D 2. althogh the bond ma not hold for e ρ, it will hold for another switching ρ t that is almost alwas consistent with The small error propert can still be generalized for the case when Q is not finite (i.e., infinitel man controllers) 20

21 Implementation isses start monitoring s How to efficientl compte a large nmber of monitoring s? It is alwas possible to write: appropriatel defined fnction wait τ D seconds From this and the definition of µ p : So, b linearit, all the µ p can be generated b: dimension is independent of the nmber of element in P Implementation isses start monitoring s When P is a continm (or ver large), it ma be isses with respect to the optimization for ρ. Things are eas, e.g., wait τ D seconds 1. P has a small nmber of elements 2. model is linearl parameterized on p (leads to k(p) qadratic) 3. there are closed form soltions (e.g., k(p) polnomial) 4. k(p) is convex on p sal reqirement in adaptive control reslts still hold if there exists a comptational dela τ C in performing the optimization, i.e. 21

22 Matching propert: p* P sch that Analsis (w = 0, ε = 0) Detectabilit propert: for frozen ρ = p P and = χ( p ) Q the injected sstem is asmptoticall stable Non-destabilizing propert: The minimm interval between consective discontinities of is τ D > 0. Small error propert (L 2 case, P finite): Matching propert: p* P sch that Analsis (w = 0, ε = 0) Detectabilit propert: for frozen ρ = p P and = χ( p ) Q the injected sstem is asmptoticall stable Non-destabilizing propert: The minimm interval between consective discontinities of is τ D > 0. Small error propert (L 2 case, P finite): Assmption (slow switching): τ D is large enogh and λ is small enogh so that the injected sstem is nif. exp. stable state transition matrix of injected sstem an switching with interval between consective discontinities no smaller than τ D 22

23 Analsis (w = 0, ε = 0) Matching propert: p* P sch that Detectabilit propert: for frozen ρ = p P and = χ( p ) Q the injected ρ sstem is asmptoticall stable v injected Non-destabilizing propert: The minimm sstem interval between consective discontinities of is τ D > 0. Small error propert (L 2 case, P finite): Assmption (slow switching): τ D is large enogh and λ is small enogh so that the injected sstem is nif. exp. stable state transition matrix of injected sstem an switching with interval between consective discontinities no smaller than τ D Analsis (w = 0, ε = 0) v injected sstem ρ ρ 1st b the Matching propert: p * P sch that 2nd b the Small error propert: 3rd b the Non-destabilization propert & assmption the injected sstem is nif. exp. stable (state transition matrix decas faster then e λ t ) 23

24 Analsis (w = 0, ε = 0) v injected sstem ρ ρ 2nd b the Small error propert: 3rd b the Non-destabilization propert & assmption the injected sstem is nif. exp. stable (state transition matrix decas faster then e λ t ) 4th b 2nd and 3rd (same for and ρ ) Analsis (w = 0, ε = 0) v injected sstem ρ ρ 4th b 2nd and 3rd (same for and ρ ) 5th b the detectabilit: state the is also L 2 and converges to zero 24

25 Analsis (w = 0, ε = 0) v injected sstem ρ ρ Theorem: Assming that P is finite and in the absence of noise and nmodeled dnamics (i.e., ε = 0, w(t) = 0, t 0) the states of the, the mlti-estimator, and the mlti-controller are all (e λt -weighted) L 2 and converge to zero as t. Matching propert: p* P sch that Analsis (general case) Detectabilit propert: for ρ = p P and = χ( p ) Q the injected sstem is asmptoticall stable Non-destabilizing propert: The minimm interval between consective discontinities of is τ D > 0. Small error propert (Q finite): p P, t 0, ρ t : [0, t ) P sch that: 1. = χ(ρ t ) except at most on m time intervals of length τ D 2. we will start b cheating and assming that = χ(ρ t ) all the time... Assmption (slow switching): τ D is large enogh and λ is small enogh so that state transition matrix of injected sstem an switching with interval between consective discontinities no smaller than τ D 25

26 Analsis (general case) v injected sstem ρ ρ Consider a fixed interval [0, t ) 1st b the Matching propert: p * P sch that 2nd b the Small error propert: ρ t sch that = χ( ρ t ) & 3rd se ρ t from Small error propert to constrct the injected sstem Analsis (general case) v injected sstem ρ t ρ t 2nd b the Small error propert: ρ t sch that = χ( ρ t ) & 3rd se ρ t from Small error propert to constrct the injected sstem since = χ( ρ t ) the injected sstem switches among stabilit matrices 4th b the Non-destabilization propert & assmption the injected sstem is nif. exp. stable (state transition matrix decas faster then e λ t ) 26

27 Analsis (general case) v injected sstem ρ t ρ t 2nd b the Small error propert: ρ t sch that = χ( ρ t ) & 4th b the Non-destabilization propert & assmption the injected sstem is nif. exp. stable (state transition matrix decas faster then e λ t ) finite λ,[0, t ) indced norm from v to : Analsis (general case) v injected sstem ρ t ρ t 2nd b the Small error propert: ρ t sch that = χ( ρ t ) & 4th b the Non-destabilization propert & assmption: 5th b small-gain argment (sing 2nd and 4th) jst as before: v bonded & injected sstem stable all s bonded 27

28 Analsis (general case) v injected sstem ρ t ρ t Where did we se = χ( ρ t )? 3rd se ρ t from Small error propert to constrct the injected sstem since = χ( ρ t ) the injected sstem switches among stabilit matrices 4th b the Non-destabilization propert & assmption the injected sstem is nif. exp. stable (state transition matrix decas faster then e λ t ) Analsis (general case) v injected sstem ρ t ρ t Where did we se = χ( ρ t )? 3rd se ρ t from Small error propert to constrct the injected sstem since = χ(ρ t ) except at most on m time intervals of length τ D, the injected sstem switches among stabilit matrices except at most on m time intervals of length τ D 4th b the Non-destabilization propert & assmption the injected sstem is still nif. exp. stable (state transition matrix decas faster then e λ t ) 28

29 Analsis (general case) v injected sstem ρ t ρ t Theorem: Assming that Q is a finite set with m element and the λ,[0, t ) norm of all s can be bonded b expressions of the form (finite indced norms) Moreover: w(t) is niforml bonded for t [0, ) all s niforml bonded w(t) 0 as t all s converge to zero as t So far Fast switching Assmption (slow switching): τ D is large enogh and λ is small enogh so that state transition matrix of injected sstem an switching with interval between consective discontinities no smaller than τ D Can be relaxed to Assmption (fast switching): λ is small enogh so that all matrices A p χ(p) λ I, p P are asmptoticall stable (an dwell-time τ D will do) 29

30 Fast switching Assmption (fast switching): λ is small enogh so that all matrices A p χ(p) λ I, p P are asmptoticall stable (an dwell-time τ D will do) Theorem: Assming that the is SISO and that the mlti-controller is realized as there exists a constant sch that when (no loss of generalit) the λ,[0, t ) norm of all s can be bonded b expressions of the form (finite indced norms) Moreover: w(t) is niforml bonded for t [0, ) all s niforml bonded w(t) 0 as t all s converge to zero as t proof: tilize internal strctre of injected sstem & inject more errors to boost rate of deca Other logics Dwell-time switching Scale-independent hsteresis switching start start hsteresis constant wait τ D seconds wait fixed amont of time wait ntil crrent monitoring becomes significantl larger than some other one n 30

31 Other logics Scale-independent hsteresis switching Hierarchical hsteresis switching start start wait ntil crrent monitoring becomes significantl larger than some other one n wait ntil crrent monitoring, or another one consistent with, becomes significantl larger than some other one n Other logics For both logics we have: Small error propert: For ever p P, t 0, switching ρ t : [0, t ) P sch that: all the time! Non-destabilizing propert: For ever p P nmber of switchings in the interval [τ, t ) nmber of elements in P (scale-indep.) or in Q (hierarchical) average dwell-time tpe growth either large h (hsteresis constant) or small λ (forgetting factor) leads to stabilit of injected sstem 31

32 Otline Spervisor control overview Estimator-based linear spervisor control Estimator-based nonlinear spervisor control Examples Estimator-based nonlinear spervisor control decision logic mltiestimator NON LINEAR switching w mlticontroller control measred otpt 32

33 Class of admissible es w exogenos distrbance/ noise control measred otpt nmodeled Process is assmed to be in a famil dnamics M p small famil of sstems arond a nominal model N p parametric ncertaint Tpicall metric on set of state-space model (?) Most reslts presented here: independent of metric d (e.g., detectabilit) or restricted to case ε p = 0 (e.g., matching) Candidate controllers Class of admissible es M p small famil of sstems arond a nominal model N p Assme given a famil of candidate controllers (withot loss of generalit all with same dimension) Mlti-controller: switching measred otpt control 33

34 Mlti-estimator measred otpt control mltiestimator How to design a mlti-estimator? we want: Matching propert: there exist some p* P sch that e p* is small Tpicall obtained b: in p * P: in M p* e p* is small when matches M p* the corresponding error mst be small Designing mlti-estimators - I Sppose nominal models N p, p P are of the form (state accessible) no exogenos inpt w state accessible Mlti-estimator: asmptoticall stable A When matches the nominal model N p* Matching propert: Assme M = { N p : p P } exponentiall p* P, c 0, λ * >0 : e p* (t) c 0 e -λ* t t 0 34

35 State-sharing Mlti-estimator: state of mlti-estimator is x E ú { z p : p P } can be large when P has man elements Sppose A p (, ) is separable, i.e., matrix vector B linearit, the p can be generated b: tre, e.g., is is linearl parameterized matrix with the size of M with ever colmn eqal to The dimension of the mlti-estimator is independent of the nmber of elements in P (cold even b infinit) Designing mlti-estimators - II Sppose nominal models N p, p P are of the form (otpt-injection awa from stable linear sstem) asmptoticall stable A p Mlti-estimator: nonlinear otpt injection (generalization of case I) When matches the nominal model N p* Matching propert: Assme M = { N p : p P } p* P, c 0, c w, λ * >0 : e p* (t) c 0 e -λ* t c w t 0 with c w = 0 in case w(t) = 0, t 0 State-sharing is possible when all A p are eqal an H p (, ) is separable: 35

36 Designing mlti-estimators - III Sppose nominal models N p, p P are of the form (otpt-inj. and coord. transf. awa from stable linear sstem) asmptoticall stable A p ζ p ú ξ p ξ p -1 (generalization of case I & II) cont. diff. coordinate transformation with continos inverse ξ p -1 (ma depend on nknown parameter p) The Matching propert is an inpt/otpt propert so the same mlti-estimator can be sed: Matching propert: Assme M = { N p : p P } p* P, c 0, c w, λ * >0 : e p* (t) c 0 e -λ* t c w t 0 with c w = 0 in case w(t) = 0, t 0 detectabilit propert? Switched sstem decision logic mltiestimator w mlticontroller switched sstem Also now the switched sstem can be seen as the interconnection of the with the injected sstem 36

37 Constrcting the injected sstem 1st Take a switching ρ : [0, ) P. 2nd Define the v ú e ρ = ρ 3rd Replace in the eqations of the mlti-estimator and mlti-controller b ρ v. mlticontroller mltiestimator ρ v Switched sstem = injected sstem w Q: How to get detectabilit on the switched sstem? A: Stabilit of the injected sstem v injected sstem ρ ρ 37

38 Stabilit & detectabilit of nonlinear sstems Stabilit: inpt small state x small Inpt-to-state stable (ISS) if β KL, γ K Integral inpt-to-state stable (iiss) if α K, β KL, γ K strictl weaker Notation: α:[0, ) [0, ) is class K continos, strictl increasing, α(0) = 0 is class K class K and nbonded β:[0, ) [0, ) [0, ) is class KL β(,t) K for fixed t & lim t β(s,t) = 0 (monotonicall) for fixed s Stabilit & detectabilit of nonlinear sstems Stabilit: inpt small state x small Inpt-to-state stable (ISS) if β KL, γ K Integral inpt-to-state stable (iiss) if α K, β KL, γ K strictl weaker One can show: 1. for ISS sstems: 0 soltion exist globall & x 0 2. for iiss sstems: 0 γ( ) < soltion exist globall & x 0 38

39 Stabilit & detectabilit of nonlinear sstems Detectabilit: inpt & otpt small state x small Detectabilit (or inpt/otpt-to-state stabilit IOSS) if β KL, γ, γ K Integral detectable (iioss) if α K, β KL, γ, γ K strictl weaker One can show: 1. for IOSS sstems:, 0 x 0 2. for iioss sstems: 0 γ ( ), 0 γ ( ) < x 0 Interconnecting stable & detectable sstems cascade (detectable) sstem 1 (stable) sstem 2 (detectable) Lemma 1 (cascade): i. sstem 1 ISS & sstem 2 detectable cascade detectable ii. sstem 1 integral ISS & sstem 2 detectable cascade integral detectable feedback (detectable) sstem 1 (detectable) Lemma 2 (feedback): i. sstem 1 detectable feedback detectable ii. sstem 1 integral detectable feedback integral detectable 39

40 Certaint Eqivalence Stabilization Theorem switched sstem (version 1) v injected sstem ρ ρ switched sstem (version 2) Σ 1 ρ Σ 2 injected sstem e ρ ρ Certaint Eqivalence Stabilization Theorem 1st detectable sstem Σ 2 detectable 2nd injected sstem ISS & 1st 3rd cascade Σ 1 Σ 2 detectable cascade Σ 1 Σ 2 detectable switched sstem detectable or 2nd injected sstem integral ISS & 1st cascade Σ 1 Σ 2 integral detectable 3rd cascade Σ 1 Σ 2 integral detectable switched sstem integral detectable switched sstem (version 2) Σ 1 ρ Σ 2 injected sstem e ρ ρ 40

41 Certaint Eqivalence Stabilization Theorem switched sstem v injected sstem ρ ρ Theorem: (Certaint Eqivalence Stabilization Theorem) Sppose the is detectable and take fixed ρ = p P and = q Q 1. injected sstem ISS switched sstem detectable. 2. injected sstem integral ISS switched sstem integral detectable Stabilit of the injected sstem is not the onl mechanism to achieve detectabilit: e.g., injected sstem i/o stable min. phase detectabilit of switched sstem (Nonlinear Certaint Eqivalence Otpt Stabilization Theorem) Achieving ISS for the injected sstem Theorem: (Certaint Eqivalence Stabilization Theorem) Sppose the is detectable and take fixed ρ = p P and = q Q 1. injected sstem ISS switched sstem detectable. 2. injected sstem integral ISS switched sstem integral detectable injected sstem mlticontroller mltiestimator ρ v ρ = p P = q Q We want to design candidate controllers that make the injected sstem (at least) integral ISS with respect to the distrbance inpt v Nonlinear robst control design problem, bt distrbance inpt v can be measred (v = e ρ = ρ ) the whole state of the injected sstem is measrable (x C, x E ) 41

42 Designing candidate controllers - I (feedback linearizable) Sppose the mlti-estimator is of the form with ( AD C p, B ) stabilizable To obtain the injected sstem, we make = p v : mst be inpt-to-state stabilized b the q = χ( p ) candidate controller (with respect to distrbance v) Candidate controller q = χ( p ): with AD C p B F p asmptoticall stable injected sstem ISS: Also applicable if the mlti-estimator is a coordinate transformation awa from this form Designing candidate controllers - II (inpt/otpt feedback linearizable) Sppose for each p P we can 1. partition the state of the mlti-estimator as 2. write its dnamics as with ( AD p C p, B p ) stabilizable ISS with respect to inpt ( x p,, ) Also applicable if the mlti-estimator is a coordinate transformation (possibl p-dependent) awa from this form This form is qite common becase tpicall one can linearize the sstem with respect to each individal p bt not with respect to all the p simltaneosl 42

43 Designing candidate controllers - II (inpt/otpt feedback linearizable) Sppose for each p P we can 1. partition the state of the mlti-estimator as 2. write its dnamics as with ( AD p C p, B p ) stabilizable ISS with respect to inpt ( x p,, ) Also applicable if the mlti-estimator is a coordinate transformation (possibl p-dependent) awa from this form Candidate controller q = χ( p ): with AD p C p B p F p asmptoticall stable injected sstem ISS: (cascade of ISS sstems) Detectabilit propert Detectabilit propert: For an of the previos mlti-estimators and candidate controller 1. The injected sstem is ISS (and also integral ISS) 2. The switched sstem is detectable throgh e p (and also integral detectable) Remarks: a. Other controllers also lead to detectabilit, e.g., one cold 1st se the feedback linearizing controller to find an ISS control Lapnov fnction 2nd se the ISS control Lapnov to constrct a more robst controller (e.g., sing an inverse optimal design) b. It is possible to achieve iiss for mch larger classes of sstems (e.g., sstems that cannot even be controlled b smooth time-invariant feedback) 43

44 Decision logic decision logic For nonlinear sstems dwell-time logics do not work becase of finite escape switching switched sstem estimation errors injected sstem Scale-independent hsteresis switching start monitoring s class K fnction from detectabilit propert hsteresis constant measre of the size of e p over a window of length 1/λ p P forgetting factor wait ntil crrent monitoring becomes significantl larger than some other one n All the µ p can be generated b a sstem with small dimension if γ p ( e p ) is separable. i.e., 44

45 Scale-independent hsteresis switching Theorem: Let P be finite with m elements. For ever p P and nmber of switchings in [τ, t ) Assme P is finite, the γ p are locall Lipschitz and maximm interval of existence of soltion niforml bonded on [0, T max ) niforml bonded on [0, T max ) Scale-independent hsteresis switching Theorem: Let P be finite with m elements. For ever p P and nmber of switchings in [τ, t ) Assme P is finite, the γ p are locall Lipschitz and maximm interval of existence of soltion Non-destabilizing propert: Switching will stop at some finite time T * [0, T max ) Small error propert: 45

46 Analsis (w = 0, no nmodeled dnamics) 1st b the Matching propert: p * P sch that e p* (t) c 0 e -λ* t t 0 2nd b the Non-destabilization propert: switching stops at a finite time T * [0, T max ) ρ(t) = p & (t) = χ(p) t [T *,T max ) 3rd b the Small error propert: 4th b the Detectabilit propert: the state x of the switched sstem is bonded on [T *,T max ) soltion exists globall T max = & x 0 as t Theorem: Assme that P is finite and all the γ p are locall Lipschitz. The state of the, mlti-estimator, mlti-controller, and all other s converge to zero as t. Otline Spervisor control overview Estimator-based linear spervisor control Estimator-based nonlinear spervisor control Examples 46

47 Example: One-link flexible maniplator x mass at the tip m t (x, t) deviation with respect to rigid bod torqe applied at the base T I H θ axis s inertia (.023) PDE (small bending): beam s elasticit transversal slice s inertia Bondar conditions: beam s length (113 cm) beam s mass densit (.68Kg total mass) Example: One-link flexible maniplator x mass at the tip m t (x, t) deviation with respect to rigid bod torqe applied at the base T I H θ axis s inertia (.023) Series expansion and trncation: eigenfnctions of the beam 47

48 Example: One-link flexible maniplator x mass at the tip m t (x, t) torqe applied at the base T Control measrements: I H θ axis s inertia (.023) Assmed not known a priori: m t [0,.1Kg] x sg [40cm, 60cm] base angle base anglar velocit tip position bending at position x sg (measred b a strain gage attached to the beam at position x sg ) Example: One-link flexible maniplator torqe Bode Diagrams Bode Diagrams 50 From: torqe 50 From: torqe 0 0 Phase (deg); Magnitde (db) To: base angle Phase (deg); Magnitde (db) To: base anglar velocit Freqenc (rad/sec) Freqenc (rad/sec) transfer fnctions as m t ranges over [0,.1Kg] and x sg ranges over [40cm, 60cm] Bode Diagrams Bode Diagrams 50 From: torqe 40 From: torqe 0 20 Phase (deg); Magnitde (db) To: tip position Phase (deg); Magnitde (db) To: bending Freqenc (rad/sec) Freqenc (rad/sec) 48

49 Example: One-link flexible maniplator torqe Class of admissible es: nknown parameter p ú (m t, x sg ) parameter set: grid of 18 points in [0,.1] [40, 60] M p famil arond a nominal transfer fnction corresponding to parameters p ú (m t, x sg ) For this problem it is not possible to write the coefficients of the nominal transfer fnctions as a fnction of the parameters becase these coefficients are the soltions to transcendental eqations that mst be compted nmericall. Famil of candidate controllers: 18 controllers designed sing LQR/LQE, one for each nominal model Example: One-link flexible maniplator m t x sg t tip position torqe set point t t 49

50 Example: One-link flexible maniplator (open-loop) Example: One-link flexible maniplator (closed-loop with fixed controller) 50

51 Example: One-link flexible maniplator (closed-loop with spervisor control) Example: Uncertain gain C 1 k < 4 The maximm gain margin achievable b a single linear time-invariant controller is 4 Dole, Francis, Tannenbam, Feedback Control Theor,

52 Example: Uncertain gain mlticontroller 1 k < 40 Anderson et. al., Mltiple Model Adaptive Control, Part 1: Finite Controller Coverings, Special George Zames Isse of IJRNC, Example: Uncertain gain otpt reference tre parameter vale monitoring s 52

53 Example: 2-dim SISO linear Class of admissible es: nominal transfer fnction nonlinear parameterized on p 1 1 nstable zero-pole cancellations An re-parameterization that makes the coefficients of the transfer fnction lie in a convex set will introdce an nstable zero-pole cancellation Bt the mlti-estimator is still separable and state-sharing can be sed Example: 2-dim SISO linear otpt reference tre parameter vale (withot noise) 53

54 Example: 2-dim SISO linear otpt reference tre parameter vale (with noise) Example: Distrbance Rejection (nknown freqenc) distrbance freqenc estimate otpt (rejection of one sinsoid) 54

55 Example: Distrbance Rejection (nknown freqenc) distrbance freqenc estimate otpt (rejection of two sinsoids) Example: Distrbance Rejection (nknown freqenc) distrbance freqenc estimate otpt (rejection of a sqare wave) 55

56 Example: Sstem in strict-feedback form Sppose nominal models N p, p P are of the form state accessible Mlti-estimator ( option I ): it is separable so we can do state-sharing When matches the nominal model N p* Matching propert: Assme M = { N p : p P } exponentiall p* P, c 0, λ * >0 : e p* (t) c 0 e -λ* t t 0 Example: Sstem in strict-feedback form Sppose nominal models N p, p P are of the form state accessible To facilitate the controller design, one can first back-step the sstem to simplif its stabilization: after the coordinate transformation the new state is no longer accessible now the control law stabilizes the sstem 56

57 Example: Sstem in strict-feedback form Sppose nominal models N p, p P are of the form Mlti-estimator ( option II ): it is separable so we can do state-sharing When matches the nominal model N p* Candidate controller q = χ(p): makes injected sstem ISS exponentiall Matching propert Detectabilit propert Example: Sstem in strict-feedback form p 2 p 1 α reference α 57

58 Example: Sstem in strict-feedback form Sppose nominal models N p, p P are of the form state accessible In the previos back-stepping procedre: the controller drives γ 0 nonlinearit is cancelled (even when p 1 < 0 and it introdces damping) One cold instead make pointwise min-norm design still leads to exponential decrease of α (withot canceling nonlinearit when p 1 < 0 ) Example: Sstem in strict-feedback form Sppose nominal models N p, p P are of the form state accessible A different recrsive procedre: In this case γ 0 exponentiall pointwise min-norm recrsive design α 0 exponentiall 58

59 Example: Sstem in strict-feedback form Sppose nominal models N p, p P are of the form Mlti-estimator ( option III ): When matches the nominal model N p* Candidate controller q = χ(p): exponentiall Matching propert Detectabilit propert Example: Sstem in strict-feedback form p 2 p 1 α reference α 59

60 Example: Sstem in strict-feedback form p 1 p 2 p 1 p 2 α α pointwise min-norm design feedback linearization design Example: Unstable-zero dnamics nknown parameter estimate otpt (stabilization) 60

61 Example: Unstable-zero dnamics nknown parameter estimate otpt (stabilization with noise) Example: Unstable-zero dnamics nknown parameter reference otpt (set-point with noise) 61

62 Example: Kinematic niccle robot x 2 θ 1 forward velocit 2 anglar velocit x 1 p 1, p 2 nknown parameters determined b the radis of the driving wheel and the distance between them This sstem cannot be stabilized b a continos time-invariant controller. The candidate controllers were themselves hbrid Example: Kinematic niccle robot 62

63 Example: Indction motor in crrent-fed mode ω is the onl measrable otpt λ 2 rotor flx 2 stator crrents ω rotor anglar velocit τ torqe generated Unknown parameters: τ L [τ min, τ max ] load torqe R [R min, R max ] rotor resistance Off-the-shelf field-oriented candidate controllers: 63

FRTN10 Exercise 12. Synthesis by Convex Optimization

FRTN10 Exercise 12. Synthesis by Convex Optimization FRTN Exercise 2. 2. We want to design a controller C for the stable SISO process P as shown in Figre 2. sing the Yola parametrization and convex optimization. To do this, the control loop mst first be

More information

Nonparametric Identification and Robust H Controller Synthesis for a Rotational/Translational Actuator

Nonparametric Identification and Robust H Controller Synthesis for a Rotational/Translational Actuator Proceedings of the 6 IEEE International Conference on Control Applications Mnich, Germany, October 4-6, 6 WeB16 Nonparametric Identification and Robst H Controller Synthesis for a Rotational/Translational

More information

Model Predictive Control Lecture VIa: Impulse Response Models

Model Predictive Control Lecture VIa: Impulse Response Models Moel Preictive Control Lectre VIa: Implse Response Moels Niet S. Kaisare Department of Chemical Engineering Inian Institte of Technolog Maras Ingreients of Moel Preictive Control Dnamic Moel Ftre preictions

More information

1. State-Space Linear Systems 2. Block Diagrams 3. Exercises

1. State-Space Linear Systems 2. Block Diagrams 3. Exercises LECTURE 1 State-Space Linear Sstems This lectre introdces state-space linear sstems, which are the main focs of this book. Contents 1. State-Space Linear Sstems 2. Block Diagrams 3. Exercises 1.1 State-Space

More information

Prototype Angle Domain Repetitive Control - Affine Parameterization Approach

Prototype Angle Domain Repetitive Control - Affine Parameterization Approach Prototpe Angle Domain Repetitive Control - Affine Parameterization Approach Perr Y. Li Professor Department of Mechanical Engineering Universit of Minnesota Minneapolis, Minnesota 55455 U.S.A. Email: lixxx99@mn.ed

More information

Chapter 4 Supervised learning:

Chapter 4 Supervised learning: Chapter 4 Spervised learning: Mltilayer Networks II Madaline Other Feedforward Networks Mltiple adalines of a sort as hidden nodes Weight change follows minimm distrbance principle Adaptive mlti-layer

More information

Classify by number of ports and examine the possible structures that result. Using only one-port elements, no more than two elements can be assembled.

Classify by number of ports and examine the possible structures that result. Using only one-port elements, no more than two elements can be assembled. Jnction elements in network models. Classify by nmber of ports and examine the possible strctres that reslt. Using only one-port elements, no more than two elements can be assembled. Combining two two-ports

More information

INPUT-OUTPUT APPROACH NUMERICAL EXAMPLES

INPUT-OUTPUT APPROACH NUMERICAL EXAMPLES INPUT-OUTPUT APPROACH NUMERICAL EXAMPLES EXERCISE s consider the linear dnamical sstem of order 2 with transfer fnction with Determine the gain 2 (H) of the inpt-otpt operator H associated with this sstem.

More information

3 2D Elastostatic Problems in Cartesian Coordinates

3 2D Elastostatic Problems in Cartesian Coordinates D lastostatic Problems in Cartesian Coordinates Two dimensional elastostatic problems are discssed in this Chapter, that is, static problems of either plane stress or plane strain. Cartesian coordinates

More information

Theoretical and Experimental Implementation of DC Motor Nonlinear Controllers

Theoretical and Experimental Implementation of DC Motor Nonlinear Controllers Theoretical and Experimental Implementation of DC Motor Nonlinear Controllers D.R. Espinoza-Trejo and D.U. Campos-Delgado Facltad de Ingeniería, CIEP, UASLP, espinoza trejo dr@aslp.mx Facltad de Ciencias,

More information

1 Differential Equations for Solid Mechanics

1 Differential Equations for Solid Mechanics 1 Differential Eqations for Solid Mechanics Simple problems involving homogeneos stress states have been considered so far, wherein the stress is the same throghot the component nder std. An eception to

More information

Chapter 3 MATHEMATICAL MODELING OF DYNAMIC SYSTEMS

Chapter 3 MATHEMATICAL MODELING OF DYNAMIC SYSTEMS Chapter 3 MATHEMATICAL MODELING OF DYNAMIC SYSTEMS 3. System Modeling Mathematical Modeling In designing control systems we mst be able to model engineered system dynamics. The model of a dynamic system

More information

Move Blocking Strategies in Receding Horizon Control

Move Blocking Strategies in Receding Horizon Control Move Blocking Strategies in Receding Horizon Control Raphael Cagienard, Pascal Grieder, Eric C. Kerrigan and Manfred Morari Abstract In order to deal with the comptational brden of optimal control, it

More information

Material. Lecture 8 Backlash and Quantization. Linear and Angular Backlash. Example: Parallel Kinematic Robot. Backlash.

Material. Lecture 8 Backlash and Quantization. Linear and Angular Backlash. Example: Parallel Kinematic Robot. Backlash. Lectre 8 Backlash and Qantization Material Toda s Goal: To know models and compensation methods for backlash Lectre slides Be able to analze the effect of qantization errors Note: We are sing analsis methods

More information

Controllers. Yale University, New Haven, CT Nagoya Works, Japan

Controllers.  Yale University, New Haven, CT Nagoya Works, Japan Servisor Control of Families of Noise Sressing Controllers Shogo Fjii Jo~ao P. Hesanha z A. Stehen Morse z shogo.fjii@ale.ed joao.hesanha@ale.ed morse@ssc.eng.ale.ed Det. of Electrical Engineering Yale

More information

Simple robustness measures for control of MISO and SIMO plants

Simple robustness measures for control of MISO and SIMO plants Preprints of the 8th IFAC World Congress Milano Ital) Agst 28 - September 2, 2 Simple robstness measres for control of MISO and SIMO plants W. P. Heath Sandira Gaadeen Control Sstems Centre, School of

More information

Geometric Image Manipulation. Lecture #4 Wednesday, January 24, 2018

Geometric Image Manipulation. Lecture #4 Wednesday, January 24, 2018 Geometric Image Maniplation Lectre 4 Wednesda, Janar 4, 08 Programming Assignment Image Maniplation: Contet To start with the obvios, an image is a D arra of piels Piel locations represent points on the

More information

The Linear Quadratic Regulator

The Linear Quadratic Regulator 10 The Linear Qadratic Reglator 10.1 Problem formlation This chapter concerns optimal control of dynamical systems. Most of this development concerns linear models with a particlarly simple notion of optimality.

More information

FEA Solution Procedure

FEA Solution Procedure EA Soltion Procedre (demonstrated with a -D bar element problem) EA Procedre for Static Analysis. Prepare the E model a. discretize (mesh) the strctre b. prescribe loads c. prescribe spports. Perform calclations

More information

A fundamental inverse problem in geosciences

A fundamental inverse problem in geosciences A fndamental inverse problem in geosciences Predict the vales of a spatial random field (SRF) sing a set of observed vales of the same and/or other SRFs. y i L i ( ) + v, i,..., n i ( P)? L i () : linear

More information

Control Systems

Control Systems 6.5 Control Systems Last Time: Introdction Motivation Corse Overview Project Math. Descriptions of Systems ~ Review Classification of Systems Linear Systems LTI Systems The notion of state and state variables

More information

Introdction In the three papers [NS97], [SG96], [SGN97], the combined setp or both eedback and alt detection lter design problem has been considered.

Introdction In the three papers [NS97], [SG96], [SGN97], the combined setp or both eedback and alt detection lter design problem has been considered. Robst Falt Detection in Open Loop s. losed Loop Henrik Niemann Jakob Stostrp z Version: Robst_FDI4.tex { Printed 5h 47m, Febrar 9, 998 Abstract The robstness aspects o alt detection and isolation (FDI)

More information

1 Introduction. r + _

1 Introduction. r + _ A method and an algorithm for obtaining the Stable Oscillator Regimes Parameters of the Nonlinear Sstems, with two time constants and Rela with Dela and Hsteresis NUŢU VASILE, MOLDOVEANU CRISTIAN-EMIL,

More information

Simplified Identification Scheme for Structures on a Flexible Base

Simplified Identification Scheme for Structures on a Flexible Base Simplified Identification Scheme for Strctres on a Flexible Base L.M. Star California State University, Long Beach G. Mylonais University of Patras, Greece J.P. Stewart University of California, Los Angeles

More information

BIOSTATISTICAL METHODS

BIOSTATISTICAL METHODS BIOSTATISTICAL METHOS FOR TRANSLATIONAL & CLINICAL RESEARCH ROC Crve: IAGNOSTIC MEICINE iagnostic tests have been presented as alwas having dichotomos otcomes. In some cases, the reslt of the test ma be

More information

Estimating models of inverse systems

Estimating models of inverse systems Estimating models of inverse systems Ylva Jng and Martin Enqvist Linköping University Post Print N.B.: When citing this work, cite the original article. Original Pblication: Ylva Jng and Martin Enqvist,

More information

Lecture Notes: Finite Element Analysis, J.E. Akin, Rice University

Lecture Notes: Finite Element Analysis, J.E. Akin, Rice University 9. TRUSS ANALYSIS... 1 9.1 PLANAR TRUSS... 1 9. SPACE TRUSS... 11 9.3 SUMMARY... 1 9.4 EXERCISES... 15 9. Trss analysis 9.1 Planar trss: The differential eqation for the eqilibrim of an elastic bar (above)

More information

Outline. Model Predictive Control: Current Status and Future Challenges. Separation of the control problem. Separation of the control problem

Outline. Model Predictive Control: Current Status and Future Challenges. Separation of the control problem. Separation of the control problem Otline Model Predictive Control: Crrent Stats and Ftre Challenges James B. Rawlings Department of Chemical and Biological Engineering University of Wisconsin Madison UCLA Control Symposim May, 6 Overview

More information

Partial Differential Equations with Applications

Partial Differential Equations with Applications Universit of Leeds MATH 33 Partial Differential Eqations with Applications Eamples to spplement Chapter on First Order PDEs Eample (Simple linear eqation, k + = 0, (, 0) = ϕ(), k a constant.) The characteristic

More information

Prandl established a universal velocity profile for flow parallel to the bed given by

Prandl established a universal velocity profile for flow parallel to the bed given by EM 0--00 (Part VI) (g) The nderlayers shold be at least three thicknesses of the W 50 stone, bt never less than 0.3 m (Ahrens 98b). The thickness can be calclated sing Eqation VI-5-9 with a coefficient

More information

Essentials of optimal control theory in ECON 4140

Essentials of optimal control theory in ECON 4140 Essentials of optimal control theory in ECON 4140 Things yo need to know (and a detail yo need not care abot). A few words abot dynamic optimization in general. Dynamic optimization can be thoght of as

More information

Gradient Projection Anti-windup Scheme on Constrained Planar LTI Systems. Justin Teo and Jonathan P. How

Gradient Projection Anti-windup Scheme on Constrained Planar LTI Systems. Justin Teo and Jonathan P. How 1 Gradient Projection Anti-windp Scheme on Constrained Planar LTI Systems Jstin Teo and Jonathan P. How Technical Report ACL1 1 Aerospace Controls Laboratory Department of Aeronatics and Astronatics Massachsetts

More information

Solving a System of Equations

Solving a System of Equations Solving a System of Eqations Objectives Understand how to solve a system of eqations with: - Gass Elimination Method - LU Decomposition Method - Gass-Seidel Method - Jacobi Method A system of linear algebraic

More information

Math 263 Assignment #3 Solutions. 1. A function z = f(x, y) is called harmonic if it satisfies Laplace s equation:

Math 263 Assignment #3 Solutions. 1. A function z = f(x, y) is called harmonic if it satisfies Laplace s equation: Math 263 Assignment #3 Soltions 1. A fnction z f(x, ) is called harmonic if it satisfies Laplace s eqation: 2 + 2 z 2 0 Determine whether or not the following are harmonic. (a) z x 2 + 2. We se the one-variable

More information

Robust Tracking and Regulation Control of Uncertain Piecewise Linear Hybrid Systems

Robust Tracking and Regulation Control of Uncertain Piecewise Linear Hybrid Systems ISIS Tech. Rept. - 2003-005 Robst Tracking and Reglation Control of Uncertain Piecewise Linear Hybrid Systems Hai Lin Panos J. Antsaklis Department of Electrical Engineering, University of Notre Dame,

More information

The Real Stabilizability Radius of the Multi-Link Inverted Pendulum

The Real Stabilizability Radius of the Multi-Link Inverted Pendulum Proceedings of the 26 American Control Conference Minneapolis, Minnesota, USA, Jne 14-16, 26 WeC123 The Real Stabilizability Radis of the Mlti-Link Inerted Pendlm Simon Lam and Edward J Daison Abstract

More information

Lecture Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2

Lecture Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2 BIJU PATNAIK UNIVERSITY OF TECHNOLOGY, ODISHA Lectre Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2 Prepared by, Dr. Sbhend Kmar Rath, BPUT, Odisha. Tring Machine- Miscellany UNIT 2 TURING MACHINE

More information

A Decomposition Method for Volume Flux. and Average Velocity of Thin Film Flow. of a Third Grade Fluid Down an Inclined Plane

A Decomposition Method for Volume Flux. and Average Velocity of Thin Film Flow. of a Third Grade Fluid Down an Inclined Plane Adv. Theor. Appl. Mech., Vol. 1, 8, no. 1, 9 A Decomposition Method for Volme Flx and Average Velocit of Thin Film Flow of a Third Grade Flid Down an Inclined Plane A. Sadighi, D.D. Ganji,. Sabzehmeidani

More information

AMS subject classifications. 47N70, 47A48, 47A20, 47A56, 47A62, 93B28, 93C55, 93D15

AMS subject classifications. 47N70, 47A48, 47A20, 47A56, 47A62, 93B28, 93C55, 93D15 OPTIMAL INPUT-OUTPUT STABILIZATION OF INFINITE-DIMENSIONAL DISCRETE TIME-INVARIANT LINEAR SYSTEMS BY OUTPUT INJECTION MARK R. OPMEER AND OLOF J. STAFFANS Abstract. We std the optimal inpt-otpt stabilization

More information

Technical Note. ODiSI-B Sensor Strain Gage Factor Uncertainty

Technical Note. ODiSI-B Sensor Strain Gage Factor Uncertainty Technical Note EN-FY160 Revision November 30, 016 ODiSI-B Sensor Strain Gage Factor Uncertainty Abstract Lna has pdated or strain sensor calibration tool to spport NIST-traceable measrements, to compte

More information

Look-Ahead Cyclic Pitch Control Using LIDAR

Look-Ahead Cyclic Pitch Control Using LIDAR Look-Ahead Cclic Pitch Control Using LIDAR David Schlipf 1 Simone Schler 2 1 Endowed Chair of Windenerg Universität Stttgart Allmandring 5B, D-7569 Stttgart schlipf@ifbni-stttgartde Patrick Gra 2 Frank

More information

A Model-Free Adaptive Control of Pulsed GTAW

A Model-Free Adaptive Control of Pulsed GTAW A Model-Free Adaptive Control of Plsed GTAW F.L. Lv 1, S.B. Chen 1, and S.W. Dai 1 Institte of Welding Technology, Shanghai Jiao Tong University, Shanghai 00030, P.R. China Department of Atomatic Control,

More information

Stability of Model Predictive Control using Markov Chain Monte Carlo Optimisation

Stability of Model Predictive Control using Markov Chain Monte Carlo Optimisation Stability of Model Predictive Control sing Markov Chain Monte Carlo Optimisation Elilini Siva, Pal Golart, Jan Maciejowski and Nikolas Kantas Abstract We apply stochastic Lyapnov theory to perform stability

More information

Sareban: Evaluation of Three Common Algorithms for Structure Active Control

Sareban: Evaluation of Three Common Algorithms for Structure Active Control Engineering, Technology & Applied Science Research Vol. 7, No. 3, 2017, 1638-1646 1638 Evalation of Three Common Algorithms for Strctre Active Control Mohammad Sareban Department of Civil Engineering Shahrood

More information

10.2 Solving Quadratic Equations by Completing the Square

10.2 Solving Quadratic Equations by Completing the Square . Solving Qadratic Eqations b Completing the Sqare Consider the eqation ( ) We can see clearl that the soltions are However, What if the eqation was given to s in standard form, that is 6 How wold we go

More information

THE foundation of many optimal controller synthesis. A System Level Approach to Controller Synthesis. arxiv: v2 [cs.

THE foundation of many optimal controller synthesis. A System Level Approach to Controller Synthesis. arxiv: v2 [cs. SUBMITTED TO IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. XX, NO. XX, OCTOBER 2016 1 A Sstem Level Approach to Controller Snthesis Yh-Shang Wang, Member, IEEE, Nikolai Matni, Member, IEEE, and John C.

More information

The Oscillatory Stable Regime of Nonlinear Systems, with two time constants

The Oscillatory Stable Regime of Nonlinear Systems, with two time constants 6th WSES International Conference on CIRCUITS SYSTEMS ELECTRONICSCONTROL & SIGNL PROCESSING Cairo Egpt Dec 9-3 7 5 The Oscillator Stable Regime of Nonlinear Sstems with two time constants NUŢU VSILE *

More information

Nonlinear parametric optimization using cylindrical algebraic decomposition

Nonlinear parametric optimization using cylindrical algebraic decomposition Proceedings of the 44th IEEE Conference on Decision and Control, and the Eropean Control Conference 2005 Seville, Spain, December 12-15, 2005 TC08.5 Nonlinear parametric optimization sing cylindrical algebraic

More information

Formal Methods for Deriving Element Equations

Formal Methods for Deriving Element Equations Formal Methods for Deriving Element Eqations And the importance of Shape Fnctions Formal Methods In previos lectres we obtained a bar element s stiffness eqations sing the Direct Method to obtain eact

More information

System identification of buildings equipped with closed-loop control devices

System identification of buildings equipped with closed-loop control devices System identification of bildings eqipped with closed-loop control devices Akira Mita a, Masako Kamibayashi b a Keio University, 3-14-1 Hiyoshi, Kohok-k, Yokohama 223-8522, Japan b East Japan Railway Company

More information

Reducing Conservatism in Flutterometer Predictions Using Volterra Modeling with Modal Parameter Estimation

Reducing Conservatism in Flutterometer Predictions Using Volterra Modeling with Modal Parameter Estimation JOURNAL OF AIRCRAFT Vol. 42, No. 4, Jly Agst 2005 Redcing Conservatism in Fltterometer Predictions Using Volterra Modeling with Modal Parameter Estimation Rick Lind and Joao Pedro Mortaga University of

More information

Chapter 2 Difficulties associated with corners

Chapter 2 Difficulties associated with corners Chapter Difficlties associated with corners This chapter is aimed at resolving the problems revealed in Chapter, which are cased b corners and/or discontinos bondar conditions. The first section introdces

More information

Nonlinear predictive control of dynamic systems represented by Wiener Hammerstein models

Nonlinear predictive control of dynamic systems represented by Wiener Hammerstein models Nonlinear Dn (26) 86:93 24 DOI.7/s7-6-2957- ORIGINAL PAPER Nonlinear predictive control of dnamic sstems represented b Wiener Hammerstein models Maciej Ławrńcz Received: 7 December 25 / Accepted: 2 Jl

More information

The Basic Feedback Loop. Design and Diagnosis. of the Basic Feedback Loop. Parallel form. The PID Controller. Tore Hägglund. The textbook version:

The Basic Feedback Loop. Design and Diagnosis. of the Basic Feedback Loop. Parallel form. The PID Controller. Tore Hägglund. The textbook version: The Basic Feedback Loop Design and Diagnosis l n of the Basic Feedback Loop sp Σ e x Controller Σ Process Σ Tore Hägglnd Department of Atomatic Control Lnd Institte of Technolog Lnd, Sweden The PID Controller

More information

Curves - Foundation of Free-form Surfaces

Curves - Foundation of Free-form Surfaces Crves - Fondation of Free-form Srfaces Why Not Simply Use a Point Matrix to Represent a Crve? Storage isse and limited resoltion Comptation and transformation Difficlties in calclating the intersections

More information

Vectors in Rn un. This definition of norm is an extension of the Pythagorean Theorem. Consider the vector u = (5, 8) in R 2

Vectors in Rn un. This definition of norm is an extension of the Pythagorean Theorem. Consider the vector u = (5, 8) in R 2 MATH 307 Vectors in Rn Dr. Neal, WKU Matrices of dimension 1 n can be thoght of as coordinates, or ectors, in n- dimensional space R n. We can perform special calclations on these ectors. In particlar,

More information

Calibration of the Characteristic Frequency of an Optical Tweezer using a Recursive Least-Squares Approach

Calibration of the Characteristic Frequency of an Optical Tweezer using a Recursive Least-Squares Approach Calibration of the Characteristic Freqenc of an Optical Tweezer sing a Recrsive Least-Sqares Approach Arna Ranaweera ranawera@engineering.csb.ed Bassam Bamieh bamieh@engineering.csb.ed Department of Mechanical

More information

Optimization via the Hamilton-Jacobi-Bellman Method: Theory and Applications

Optimization via the Hamilton-Jacobi-Bellman Method: Theory and Applications Optimization via the Hamilton-Jacobi-Bellman Method: Theory and Applications Navin Khaneja lectre notes taken by Christiane Koch Jne 24, 29 1 Variation yields a classical Hamiltonian system Sppose that

More information

Multivariable Ripple-Free Deadbeat Control

Multivariable Ripple-Free Deadbeat Control Scholarly Jornal of Mathematics and Compter Science Vol. (), pp. 9-9, October Available online at http:// www.scholarly-jornals.com/sjmcs ISS 76-8947 Scholarly-Jornals Fll Length Research aper Mltivariable

More information

Sources of Non Stationarity in the Semivariogram

Sources of Non Stationarity in the Semivariogram Sorces of Non Stationarity in the Semivariogram Migel A. Cba and Oy Leangthong Traditional ncertainty characterization techniqes sch as Simple Kriging or Seqential Gassian Simlation rely on stationary

More information

Elements of Coordinate System Transformations

Elements of Coordinate System Transformations B Elements of Coordinate System Transformations Coordinate system transformation is a powerfl tool for solving many geometrical and kinematic problems that pertain to the design of gear ctting tools and

More information

Control Performance Monitoring of State-Dependent Nonlinear Processes

Control Performance Monitoring of State-Dependent Nonlinear Processes Control Performance Monitoring of State-Dependent Nonlinear Processes Lis F. Recalde*, Hong Ye Wind Energy and Control Centre, Department of Electronic and Electrical Engineering, University of Strathclyde,

More information

An extremum seeking approach to parameterised loop-shaping control design

An extremum seeking approach to parameterised loop-shaping control design Preprints of the 19th World Congress The International Federation of Atomatic Control An extremm seeking approach to parameterised loop-shaping control design Chih Feng Lee Sei Zhen Khong Erik Frisk Mattias

More information

Bridging the Gap Between Multigrid, Hierarchical, and Receding-Horizon Control

Bridging the Gap Between Multigrid, Hierarchical, and Receding-Horizon Control Bridging the Gap Between Mltigrid, Hierarchical, and Receding-Horizon Control Victor M. Zavala Mathematics and Compter Science Division Argonne National Laboratory Argonne, IL 60439 USA (e-mail: vzavala@mcs.anl.gov).

More information

Complex Variables. For ECON 397 Macroeconometrics Steve Cunningham

Complex Variables. For ECON 397 Macroeconometrics Steve Cunningham Comple Variables For ECON 397 Macroeconometrics Steve Cnningham Open Disks or Neighborhoods Deinition. The set o all points which satis the ineqalit

More information

Setting The K Value And Polarization Mode Of The Delta Undulator

Setting The K Value And Polarization Mode Of The Delta Undulator LCLS-TN-4- Setting The Vale And Polarization Mode Of The Delta Undlator Zachary Wolf, Heinz-Dieter Nhn SLAC September 4, 04 Abstract This note provides the details for setting the longitdinal positions

More information

Dynamic Optimization of First-Order Systems via Static Parametric Programming: Application to Electrical Discharge Machining

Dynamic Optimization of First-Order Systems via Static Parametric Programming: Application to Electrical Discharge Machining Dynamic Optimization of First-Order Systems via Static Parametric Programming: Application to Electrical Discharge Machining P. Hgenin*, B. Srinivasan*, F. Altpeter**, R. Longchamp* * Laboratoire d Atomatiqe,

More information

Convergence analysis of ant colony learning

Convergence analysis of ant colony learning Delft University of Technology Delft Center for Systems and Control Technical report 11-012 Convergence analysis of ant colony learning J van Ast R Babška and B De Schtter If yo want to cite this report

More information

Subcritical bifurcation to innitely many rotating waves. Arnd Scheel. Freie Universitat Berlin. Arnimallee Berlin, Germany

Subcritical bifurcation to innitely many rotating waves. Arnd Scheel. Freie Universitat Berlin. Arnimallee Berlin, Germany Sbcritical bifrcation to innitely many rotating waves Arnd Scheel Institt fr Mathematik I Freie Universitat Berlin Arnimallee 2-6 14195 Berlin, Germany 1 Abstract We consider the eqation 00 + 1 r 0 k2

More information

HADAMARD-PERRON THEOREM

HADAMARD-PERRON THEOREM HADAMARD-PERRON THEOREM CARLANGELO LIVERANI. Invariant manifold of a fixed point He we will discss the simplest possible case in which the existence of invariant manifolds arises: the Hadamard-Perron theorem.

More information

Optimal Control of a Heterogeneous Two Server System with Consideration for Power and Performance

Optimal Control of a Heterogeneous Two Server System with Consideration for Power and Performance Optimal Control of a Heterogeneos Two Server System with Consideration for Power and Performance by Jiazheng Li A thesis presented to the University of Waterloo in flfilment of the thesis reqirement for

More information

DIGITAL LINEAR QUADRATIC SMITH PREDICTOR

DIGITAL LINEAR QUADRATIC SMITH PREDICTOR DIGITAL LINEAR QUADRATIC SMITH PREDICTOR Vladimír Bobál,, Marek Kbalčík, Petr Dostál, and Stanislav Talaš Tomas Bata Universit in Zlín Centre of Polmer Sstems, Universit Institte Department of Process

More information

Model Discrimination of Polynomial Systems via Stochastic Inputs

Model Discrimination of Polynomial Systems via Stochastic Inputs Model Discrimination of Polynomial Systems via Stochastic Inpts D. Georgiev and E. Klavins Abstract Systems biologists are often faced with competing models for a given experimental system. Unfortnately,

More information

IJSER. =η (3) = 1 INTRODUCTION DESCRIPTION OF THE DRIVE

IJSER. =η (3) = 1 INTRODUCTION DESCRIPTION OF THE DRIVE International Jornal of Scientific & Engineering Research, Volme 5, Isse 4, April-014 8 Low Cost Speed Sensor less PWM Inverter Fed Intion Motor Drive C.Saravanan 1, Dr.M.A.Panneerselvam Sr.Assistant Professor

More information

Lecture: Corporate Income Tax - Unlevered firms

Lecture: Corporate Income Tax - Unlevered firms Lectre: Corporate Income Tax - Unlevered firms Ltz Krschwitz & Andreas Löffler Disconted Cash Flow, Section 2.1, Otline 2.1 Unlevered firms Similar companies Notation 2.1.1 Valation eqation 2.1.2 Weak

More information

Hybrid modelling and model reduction for control & optimisation

Hybrid modelling and model reduction for control & optimisation Hybrid modelling and model redction for control & optimisation based on research done by RWTH-Aachen and TU Delft presented by Johan Grievink Models for control and optimiation market and environmental

More information

Chapter 3. Preferences and Utility

Chapter 3. Preferences and Utility Chapter 3 Preferences and Utilit Microeconomics stdies how individals make choices; different individals make different choices n important factor in making choices is individal s tastes or preferences

More information

Fuzzy Control of a Nonlinear Deterministic System for Different Operating Points

Fuzzy Control of a Nonlinear Deterministic System for Different Operating Points International Jornal of Electrical & Compter Sciences IJECS-IJE Vol: No: 9 Fzz Control of a Nonlinear Deterministic Sstem for Different Operating Points Gonca Ozmen Koca, Cafer Bal, Firat Universit, Technical

More information

Chapter 1: Differential Form of Basic Equations

Chapter 1: Differential Form of Basic Equations MEG 74 Energ and Variational Methods in Mechanics I Brendan J. O Toole, Ph.D. Associate Professor of Mechanical Engineering Howard R. Hghes College of Engineering Universit of Nevada Las Vegas TBE B- (7)

More information

FREQUENCY DOMAIN FLUTTER SOLUTION TECHNIQUE USING COMPLEX MU-ANALYSIS

FREQUENCY DOMAIN FLUTTER SOLUTION TECHNIQUE USING COMPLEX MU-ANALYSIS 7 TH INTERNATIONAL CONGRESS O THE AERONAUTICAL SCIENCES REQUENCY DOMAIN LUTTER SOLUTION TECHNIQUE USING COMPLEX MU-ANALYSIS Yingsong G, Zhichn Yang Northwestern Polytechnical University, Xi an, P. R. China,

More information

Linear and Nonlinear Model Predictive Control of Quadruple Tank Process

Linear and Nonlinear Model Predictive Control of Quadruple Tank Process Linear and Nonlinear Model Predictive Control of Qadrple Tank Process P.Srinivasarao Research scholar Dr.M.G.R.University Chennai, India P.Sbbaiah, PhD. Prof of Dhanalaxmi college of Engineering Thambaram

More information

Simulation Based Analysis of Two Different Control Strategies for PMSM

Simulation Based Analysis of Two Different Control Strategies for PMSM International Jornal of Engineering Trends and Technology (IJETT) - Volme4Isse4- April 23 Simlation Based Analysis of Two Different Control Strategies for PMSM Lindita Dhamo #, Aida Spahi #2 # Department

More information

Linear Strain Triangle and other types of 2D elements. By S. Ziaei Rad

Linear Strain Triangle and other types of 2D elements. By S. Ziaei Rad Linear Strain Triangle and other tpes o D elements B S. Ziaei Rad Linear Strain Triangle (LST or T6 This element is also called qadratic trianglar element. Qadratic Trianglar Element Linear Strain Triangle

More information

A NEW EVALUATION PLATFORM FOR NAVIGATION SYSTEMS

A NEW EVALUATION PLATFORM FOR NAVIGATION SYSTEMS A NEW EVALUATION PLATFORM FOR NAVIGATION SYSTEMS Thomas Hanefeld Sejerøe Niels Kjølstad polsen Ole Ravn Ørsted DTU, Atomation, The Technical Universit of Denmar, Bilding 6, DK-8 Kgs. Lngb, Denmar E-mail:9977@stdent.dt.d,

More information

Sensitivity Analysis in Bayesian Networks: From Single to Multiple Parameters

Sensitivity Analysis in Bayesian Networks: From Single to Multiple Parameters Sensitivity Analysis in Bayesian Networks: From Single to Mltiple Parameters Hei Chan and Adnan Darwiche Compter Science Department University of California, Los Angeles Los Angeles, CA 90095 {hei,darwiche}@cs.cla.ed

More information

Lab Manual for Engrd 202, Virtual Torsion Experiment. Aluminum module

Lab Manual for Engrd 202, Virtual Torsion Experiment. Aluminum module Lab Manal for Engrd 202, Virtal Torsion Experiment Alminm modle Introdction In this modle, o will perform data redction and analsis for circlar cross section alminm samples. B plotting the torqe vs. twist

More information

Robust H 2 Synthesis for Dual-stage Multi-sensing Track-following Servo Systems in HDDs

Robust H 2 Synthesis for Dual-stage Multi-sensing Track-following Servo Systems in HDDs Proceedings of the 6 American Control Conference Minneapolis Minnesota USA Jne 14-16 6 WeB18.1 Robst H Snthesis for al-stage Mlti-sensing Track-following Servo Sstems in Hs Rozo Nagamne Xinghi Hang and

More information

4.2 First-Order Logic

4.2 First-Order Logic 64 First-Order Logic and Type Theory The problem can be seen in the two qestionable rles In the existential introdction, the term a has not yet been introdced into the derivation and its se can therefore

More information

Lecture: Corporate Income Tax

Lecture: Corporate Income Tax Lectre: Corporate Income Tax Ltz Krschwitz & Andreas Löffler Disconted Cash Flow, Section 2.1, Otline 2.1 Unlevered firms Similar companies Notation 2.1.1 Valation eqation 2.1.2 Weak atoregressive cash

More information

Deakin Research Online

Deakin Research Online Deakin Research Online his is the pblished version: Herek A. Samir S Steven rinh Hie and Ha Qang P. 9 Performance of first and second-order sliding mode observers for nonlinear systems in AISA 9 : Proceedings

More information

MULTIOBJECTIVE OPTIMAL STRUCTURAL CONTROL OF THE NOTRE DAME BUILDING MODEL BENCHMARK *

MULTIOBJECTIVE OPTIMAL STRUCTURAL CONTROL OF THE NOTRE DAME BUILDING MODEL BENCHMARK * Thi d d i h F M k 4 0 4 Earthqake Engineering and Strctral Dynamics, in press. MULTIOBJECTIVE OPTIMAL STRUCTURAL CONTROL OF THE NOTRE DAME BUILDING MODEL BENCHMARK * E. A. JOHNSON, P. G. VOULGARIS, AND

More information

EE2 Mathematics : Functions of Multiple Variables

EE2 Mathematics : Functions of Multiple Variables EE2 Mathematics : Fnctions of Mltiple Variables http://www2.imperial.ac.k/ nsjones These notes are not identical word-for-word with m lectres which will be gien on the blackboard. Some of these notes ma

More information

Lecture 6 : Linear Fractional Transformations (LFTs) Dr.-Ing. Sudchai Boonto

Lecture 6 : Linear Fractional Transformations (LFTs) Dr.-Ing. Sudchai Boonto Lectre 6 : (LFTs) Dr-Ing Sdchai Boonto Department of Control System and Instrmentation Engineering King Mongkts Unniversity of Technology Thonbri Thailand Feedback Strctre d i d r e y z K G g n The standard

More information

STEP Support Programme. STEP III Hyperbolic Functions: Solutions

STEP Support Programme. STEP III Hyperbolic Functions: Solutions STEP Spport Programme STEP III Hyperbolic Fnctions: Soltions Start by sing the sbstittion t cosh x. This gives: sinh x cosh a cosh x cosh a sinh x t sinh x dt t dt t + ln t ln t + ln cosh a ln ln cosh

More information

AN ALTERNATIVE DECOUPLED SINGLE-INPUT FUZZY SLIDING MODE CONTROL WITH APPLICATIONS

AN ALTERNATIVE DECOUPLED SINGLE-INPUT FUZZY SLIDING MODE CONTROL WITH APPLICATIONS AN ALTERNATIVE DECOUPLED SINGLE-INPUT FUZZY SLIDING MODE CONTROL WITH APPLICATIONS Fang-Ming Y, Hng-Yan Chng* and Chen-Ning Hang Department of Electrical Engineering National Central University, Chngli,

More information

AN ALTERNATIVE DECOUPLED SINGLE-INPUT FUZZY SLIDING MODE CONTROL WITH APPLICATIONS

AN ALTERNATIVE DECOUPLED SINGLE-INPUT FUZZY SLIDING MODE CONTROL WITH APPLICATIONS AN ALTERNATIVE DECOUPLED SINGLE-INPUT FUZZY SLIDING MODE CONTROL WITH APPLICATIONS Fang-Ming Y, Hng-Yan Chng* and Chen-Ning Hang Department of Electrical Engineering National Central University, Chngli,

More information

III. Demonstration of a seismometer response with amplitude and phase responses at:

III. Demonstration of a seismometer response with amplitude and phase responses at: GG5330, Spring semester 006 Assignment #1, Seismometry and Grond Motions De 30 Janary 006. 1. Calibration Of A Seismometer Using Java: A really nifty se of Java is now available for demonstrating the seismic

More information

1 Undiscounted Problem (Deterministic)

1 Undiscounted Problem (Deterministic) Lectre 9: Linear Qadratic Control Problems 1 Undisconted Problem (Deterministic) Choose ( t ) 0 to Minimize (x trx t + tq t ) t=0 sbject to x t+1 = Ax t + B t, x 0 given. x t is an n-vector state, t a

More information

DESIGN AND SIMULATION OF SELF-TUNING PREDICTIVE CONTROL OF TIME-DELAY PROCESSES

DESIGN AND SIMULATION OF SELF-TUNING PREDICTIVE CONTROL OF TIME-DELAY PROCESSES DESIGN AND SIMULAION OF SELF-UNING PREDICIVE CONROL OF IME-DELAY PROCESSES Vladimír Bobál,, Marek Kbalčík and Petr Dostál, omas Bata University in Zlín Centre of Polymer Systems, University Institte Department

More information

Stability analysis of two predator-one stage-structured prey model incorporating a prey refuge

Stability analysis of two predator-one stage-structured prey model incorporating a prey refuge IOSR Jornal of Mathematics (IOSR-JM) e-issn: 78-578 p-issn: 39-765X. Volme Isse 3 Ver. II (Ma - Jn. 5) PP -5 www.iosrjornals.org Stabilit analsis of two predator-one stage-strctred pre model incorporating

More information