L 1 Adaptive Controller for Nonlinear Systems in the Presence of

Size: px
Start display at page:

Download "L 1 Adaptive Controller for Nonlinear Systems in the Presence of"

Transcription

1 28 American Conrol Conference Wesin Seale Hoel, Seale, Washingon, USA June -3, 28 FrA4.3 L Adapive Conroller for Nonlinear Sysems in he Presence of Unmodelled Dynamics: Par II Chengyu Cao and Naira Hovakimyan Absrac This paper presens a novel adapive conrol mehodology for a class of uncerain nonlinear sysems in he presence of unmodelled dynamics. The adapive conroller ensures uniformly bounded ransien response for sysem s boh inpu and oupu signals simulaneously. The performance bounds can be sysemaically improved by increasing he adapaion gain. I. INTRODUCTION Reference [] considers a class of uncerain sysems in he presence of ime and sae dependen unknown nonlineariies and develops an adapive conrol mehodology ha ensures uniformly bounded ransien response for sysem s inpu/oupu signals simulaneously. This paper eends he resuls of [] o a class of uncerain nonlinear sysems in he presence of unmodeled dynamics. The L norm bounds for he error signals beween he closed-loop adapive sysem and he closed-loop reference sysem can be sysemaically reduced by increasing he adapaion gain. The paper is organied as follows. Secion II gives he problem formulaion. In Secion III, he novel L adapive conrol archiecure is presened. Sabiliy and uniform performance bounds are presened in Secion IV. In Secion V, simulaion resuls are presened, while Secion VI concludes he paper. II. PROBLEM FORMULATION Consider he following sysem dynamics: ẋ = A m + bωu + f,, = g o,, ẋ = g,,, y = c, =, where R n is he sysem sae vecor measurable, u R is he conrol signal, y R is he regulaed oupu, b,c R n are known consan vecors, A m is a known n n Hurwi mari, ω R is he unknown conrol effeciveness, and are he oupu and he sae vecor of he unmodelled dynamics, while f, g o, g are unknown nonlinear funcions. Assumpion : [Known sign for conrol effeciveness] There eis ω u > ω l > such ha ω l ω ω u. Assumpion 2: [Sabiliy of inernal dynamics] The - dynamics are bounded inpu bounded oupu BIBO sable, i.e. here eis L > and L 2 > such ha L L L + L 2. 2 Research is suppored by AFOSR under Conrac No. FA The auhors are wih Aerospace & Ocean Engineering, Virginia Polyechnic Insiue & Sae Universiy, Blacksburg, VA , {chengyu, nhovakim@v.edu} Furher, le X [ ]. Assumpion 3: [Semiglobal Lipschi condiion] For any δ >, here eis posiive K δ and B such ha fx, fx 2, K δ X X 2, 3 f, B 4 for all X i δ, i =,2, uniformly in. Assumpion 4: [Semiglobal uniform boundedness of parial derivaives] For any δ >, here eis d f δ >, and d f δ > such ha for any δ, he parial derivaives of fx, are piece-wise coninuous and bounded fx, d f δ, fx, d f δ. 5 X The conrol objecive is o design a full-sae feedback adapive conroller o ensure ha y racks a given bounded reference signal r boh in ransien and seady sae, while all oher error signals remain bounded. III. L ADAPTIVE CONTROLLER The design of L adapive conroller involves a sricly proper ransfer funcion Ds and a gain k R +, which leads o a sricly proper sable Cs = ωkds + ωkds wih DC gain C =. The simples choice of Ds is 6 Ds = s, 7 which yields a firs order sricly proper Cs in he following form: Cs = ωk s + ωk. 8 Le Hs = si A m b, 9 and r be he signal wih is Laplace ransformaion si A m. Since A m is Hurwi and is finie, r L is finie. Furher, for every δ > le where L δ δ δ K δ, δ ma{δ + γ, L δ + γ + L 2 }, and γ is an arbirary posiive consan /8/$ AACC. 499 Auhoried licensed use limied o: UNIVERSITY OF CONNECTICUT. Downloaded on February 2, 29 a 9:44 from IEEE Xplore. Resricions apply.

2 For he proofs of sabiliy and performance bounds, he choice of Ds and k furher needs o ensure ha here eiss ρ r such ha: Gs L < ρ r k g CsHs L r L where r L /ρ r L ρr + B, 2 Gs = Hs Cs, and k g = c A m b. 3 We consider he following sae predicor or passive idenifier for generaion of he adapive laws: ˆ = A mˆ + b ˆωu + ˆθ L + ˆσ ŷ = c ˆ, ˆ =. 4 The adapive esimaes ˆω, ˆθ, ˆσ are defined as: ˆθ = ΓProjˆθ, L Pb, ˆθ = ˆθ ˆσ = ΓProjˆσ, Pb, ˆσ = ˆσ 5 ˆω = ΓProjˆω, Pbu, ˆω = ˆω, where = ˆ, Γ R + is he adapaion gain, P is he soluion of he algebraic equaion A mp +PA m = Q, Q >, and he projecion operaor ensures ha he adapive esimaes ˆω, ˆθ, ˆσ remain inside he compac ses [ω l, ω u ], [ θ b, θ b ], [ σ b, σ b ], respecively, wih θ b, σ b defined as follows: θ b = L ρ, σ b = B + ǫ, 6 where ǫ is arbirary posiive consan, and ρ = ρ r + β, 7 wih arbirary β > γ. Remark : In he following analysis, we demonsrae ha ρ r and ρ characerie he domain of aracion of he closed loop reference sysem ye o be defined and he sysem in respecively. Since γ and β can be se arbirarily small, ρ can approimae ρ r arbirarily closely. The conrol signal is generaed hrough gain feedback of he following sysem: χs = Ds rs, us = kχs, 8 where k R + was inroduced in 6, while rs is he Laplace ransformaion of he signal r = ˆωu + ˆθ L + ˆσ k g r. 9 The complee L adapive conroller consiss of 4, 5 and 8, subjec o he L -gain upper bound in 2. As compared o he corresponding L adapive conrol archiecure for sysems wihou unmodeled dynamics in [], he only difference is ha here we use L in 4 and 5 insead of. Consequenly, he subsequen analysis is also principally differen from he one in []. IV. ANALYSIS OF L ADAPTIVE CONTROLLER A. Reference Sysem We now consider he following closed-loop reference sysem wih is conrol signal and sysem response being defined as follows: ẋ ref = A m ref + b f ref,, +ωu ref, ref =, 2 u ref s = Cs/ωk g rs r ref s, 2 y ref = c ref, 22 where r ref s is he Laplace ransformaion of he signal r ref = f ref,,, and k g is inroduced in 3. The ne Lemma esablishes he sabiliy of he closed-loop sysem in Lemma : For he closed-loop reference sysem in 2-22, subjec o he L -gain upper bound in 2, if < ρ r and hen L L ref L + γ + L 2, 23 where ρ r is inroduced in 2 and ref L < ρ r, 24 u ref L < ρ ur, 25 ρ ur = Cs/ω L ρ r L ρr + B + k g r L. Proof. I follows from 2-22 ha ref s = Gs r ref s+hscsk g rs+si A m. 26 Eample 5.2 in [2] page 99 implies ha refτ L Gs L r refτ L + k g CsHs L r τ L + r L. 27 If 24 is no rue, since ref = < ρ r and ref is coninuous, here eiss τ [, ] such ha I follows from 23 and 28 ha and hence refτ L ρ r, 28 ref τ = ρ r. 29 τ L L ρ r + γ + L 2, X τ L ρ r ma{ρ r + γ,l ρ r + γ + L 2 }. 3 Furher, i follows from Assumpion 3 ha r refτ L K ρr ρ r + B, and he redefiniion in leads o he following upper bound r refτ L L ρr ρ r + B. 3 Since r τ L r L, i follows from 27 ha refτ L Gs L L ρr ρ r + r L + k g CsHs L r L + Gs L B Auhoried licensed use limied o: UNIVERSITY OF CONNECTICUT. Downloaded on February 2, 29 a 9:44 from IEEE Xplore. Resricions apply.

3 The condiion in 2 can be solved for ρ r o obain he following upper bound leading o Gs L L ρr ρ r + k g CsHs L r L + r L + Gs L B < ρ r, 33 refτ L < ρ r, 34 which conradics 29 and proves he upper bound in 24. This furher implies ha he upper bound in 3 holds for any τ, i.e. r ref L < ρ r L ρr + B. 35 Eample 5.2 in [2] page 99 leads o he following bound u ref L < Cs/ω L ρ r L ρr + B + k g r L, 36 which proves 25. B. Equivalen Linear Time-Varying Sysem In his secion we demonsrae ha he sysem wih unmodelled dynamics in can be ransformed ino a linear sysem wih unknown ime-varying parameers. To sreamline he subsequen analysis, we need o inroduce several noaions. Le γ be he desired performance bound for L, and β be an arbirary posiive consan verifying he following equaliy Cs L Gs L K ρr γ + β = γ, 37 where ρ r was defined in 3, and γ was inroduced in. I follows from 2 ha Gs L L ρr <. Since K ρr < L ρr, we have Gs L K ρr >, which implies ha he condiion in 37 can be always saisfied by selecing γ and β sufficienly small. Furher, le ρ u = ρ ur + γ 2, 38 Cs γ 2 = K ρr γ Cs ω + γ L ω c o H o s c o 39 L I follows from Lemma 4 in [3] ha here eiss c o R n such ha c o Hs = N ns N d s, 4 where he order of N d s is one more han he order of N n s, and boh N n s and N d s are sable polynomials. Lemma 2: For he sysem in, if L ρ, u L ρ u, 4 here eis differeniable θτ and στ wih bounded θτ and στ over τ [, ] such ha θτ < θ b, 42 στ < σ b, 43 fτ,τ,τ = θτ τ L + στ. 44 The proof is similar o he proof of Lemma 2 in [] and is herefore omied. If 4 holds, Lemma 2 implies ha he sysem in can be rewrien over τ [, ] as: ẋτ = A m τ + bθτ τ L + ωuτ + στ, yτ = c τ, =, 45 where θτ, στ are unknown ime-varying signals subjec o for all τ [, ], and heir derivaives are also uniformly bounded for all τ [, ]: θτ d θ ρ,ρ u <, στ d σ ρ,ρ u <. 46 Le θτ ˆθτ θτ, ωτ ˆωτ ω, στ ˆστ στ. 47 I follows from 4 and 45 ha for all τ [,] τ = A m τ + b ωτuτ + θτ τ L + στ 48 and =. C. Tracking error signal Lemma 3: For he sysem in and he L adapive conroller in 4, 5 and 8, if hen where L ρ, u ρ u, 49 L θ m ρ,ρ u λ min PΓ, 5 θ m ρ,ρ u 4θ 2 b + 4σ 2 b + ω u ω l 2 +4 λ map λ min Q θ bd θ ρ,ρ u + σ b d σ ρ,ρ u. 5 The proof is similar o Lemma 3 in []. D. Transien and Seady-Sae Performance The block diagram of he closed-loop sysem wih L adapive conroller and he reference sysem is illusraed in Figure. We noe ha he reference sysem is no implemenable since i uses he unknown signal and he unknown funcion f. This closed-loop sysem is only used for analysis purposes and does no affec he implemenaion of L adapive conroller. In he following Theorem, we prove he sabiliy and he ransien performance of he closed-loop sysem wih L adapive conroller. Theorem : Consider he closed-loop sysem wih L adapive conroller defined via 4, 5, 8, subjec o 2, and he reference sysem in If ρ r, 52 and he adapive gain verifies he lower bound Γ > θ mρ,ρ u λ min Pγ 2, 53 4 Auhoried licensed use limied o: UNIVERSITY OF CONNECTICUT. Downloaded on February 2, 29 a 9:44 from IEEE Xplore. Resricions apply.

4 Reference Sysem Fig.. sysem. we have: f ref,, u Cs + + ωu ref Sysem = g,, & g,, o = f,, + ω u L Adapive Conroller k g si A m si A m b b ref Closed-loop sysem wih L adapive conroller and reference ref L ρ r, 54 u ref L ρ ur, 55 r r L γ, 56 ref L < γ, 57 u u ref L < γ 2, 58 where γ 2 is defined in 39. Proof. The proof will be done by conradicion. Assume ha are no rue. Then, since ref = γ, u u ref =, and τ, ref τ, uτ, u ref τ are coninuous, here eiss such ha while ref = γ, or u u ref = γ 2, 59 ref L γ, u u ref L γ 2. 6 On he oher hand, i follows from Assumpion 2 ha L L ref L + γ + L 2. 6 Consequenly, Lemma implies ha ref L ρ r, u ref L ρ ur. 62 Taking ino consideraion he definiions from 7 and 38, i follows from 6 ha L ρ r + γ ρ, u L ρ ur + γ 2 = ρ u. 63 Choosing he adapive gain according o 53, Lemma 3 implies ha L γ. 64 Le rτ = ωτuτ + θτ τ L + στ, r τ = θτ τ L + στ. I follows from 8 ha χs = Dsωus+r s k g rs+ rs, where rs and r s are he Laplace ransformaions of signals rτ and r τ. Consequenly χs = Ds + kωds r s k g rs + rs, us = kds + kωds r s k g rs + rs. 65 Using he definiion of Cs from 6, we can wrie ωus = Csr s k g rs + rs, 66 and he sysem in consequenly akes he form: s = Hs Csr s + Csk g rs Cs rs + si A m. 67 Le eτ = τ ref τ. Then from 26 we have es = Hs Csr 2 s Cs rs, e =, 68 where r 2 s is he Laplace ransformaion of he signal r 2 τ = θτ τ L refτ L. 69 Since 44 implies ha fτ,τ,τ = θτ τ L + στ, f ref τ,τ,τ = θτ refτ L + στ, i follows from 69 ha r 2 τ = fτ,τ,τ f ref τ,τ,τ. 7 Eample 5.2 in [2] page 99 gives he following upper bound: e L Gs L r 2 L + r 3 L, 7 where r 3 τ is he signal wih is Laplace ransformaion being r 3 s = CsHs rs. From he error dynamics in 48 we have s = Hs rs, which leads o r 3 s = Cs s, and hence r 3 L Cs L L. Subsiuing 62 ino 6, we obain and hence L L ρ r + γ + L 2, X L ma{ρ r + γ,l ρ r + γ + L 2 }. 72 Similarly, we have X ref L X L ma{ρ r +γ,l ρ r +γ +L 2 }, 73 where X ref = [ ref ]. I follows from 72, 73 ha Xτ ρ r, X ref τ ρ r, over τ [, ], and hence Assumpion 3 implies ha r 2 τ K ρr eτ over τ [, ], which implies ha r 2 L K ρr e L. 74 From 7 we have e L Gs L K ρr e L + Cs L L. The upper bound in 64 and he L - gain upper bound in 2 lead o he following upper bound 42 Auhoried licensed use limied o: UNIVERSITY OF CONNECTICUT. Downloaded on February 2, 29 a 9:44 from IEEE Xplore. Resricions apply.

5 e L γ Cs L / Gs L K ρr, which along wih 37 leads o e L γ β < γ. 75 We noice ha from 2 and 66 ha one can derive us u ref s = Cs ω r 2s r 4 s, where r 4 s = Cs ω rs. Therefore, i follows from Eample 5.2 in [2] page 99 ha u u ref L Cs/ω L r 2 L + r 4 L, and hence u u ref L Cs/ω L K ρr e L + r 4 L. 76 Cs We have r 4 s = ω c o Hsc o Hs rs = Cs ω c o Hsc o s, where c o is inroduced in 4. Using he polynomials from 4, we can wrie ha Cs Cs N ω = d s c o Hs ω N ns. Since Cs is sable and sricly proper, he complee sysem Cs c o Hs is proper and sable, which implies ha is L gain eiss and is finie. Hence, we have r 4 L Cs ω c o Hsc o L. L The upper bound in 63 leads o r 4 L Cs ω c o Hs c o γ. 77 L I follows from 75, 76 and 77 and he definiion of γ 2 in 39 ha u u ref L K ρr Cs/ω L γ β + Cs ω c o Hs c o γ < γ L We noe ha he upper bounds in 75 and 78 conradic he equaliy in 59, which proves Since are uniform bounds for all, he upper bound in 56 follows from 64 direcly, while he upper bounds in follow from 62 correspondingly. I follows from 53 ha we can achieve arbirarily small γ by increasing he adapive gain. V. SIMULATIONS Consider he dynamics of a single-link robo arm roaing on a verical plane: I q + Fq, q,, = u, 79 where q and q are measured angular posiion and velociy, respecively, u is he inpu orque, I is he unknown momen of ineria, represens he oupu of unmodelled dynamics, Fq, q,, is an unknown nonlineariy ha lumps he forces and orques due o graviy, fricion, disurbance, oher eernal sources and unmodelled dynamics. The conrol objecive is o design u o achieve racking of bounded reference inpu r by q, where r L. Le = [ 2 ] = [q q]. The sysem in 79 can be presened in he canonical form: ẋ = A m + bωu + f,, y = c, [ ] where b = [ ], c = [ ], A m =, ω =.4 /I is he unknown conrol effeciveness, and he unknown funcion f is given by: f, = [.4] F 2,,,. Le he unknown conrol effeciveness, he unknown Fig. 2. Gs L L ρ wih respec o ωk. nonlineariy be given by ω = /I = and F, 2,, = , while is he oupu of he unmodelled dynamics, given by: s = s s 2 + 3s + 2 σs, σ = sin.2 + 2, so ha he compac ses can be conservaively chosen according o he following upper and lower bounds ω l =.5, ω u = 2, θ b = 2,σ b = Time a solid, ˆ dashed, and rdoed ime b Time-hisory of u Fig. 3. Performance of L adapive conroller for σ = sin In he implemenaion of[ he L adapive conroller, ] we se Q = 2I and hence P = Firs, we need o verify he condiion in 2. Leing s Ds = /s, we have Gs = s+ωkhs, Hs = s [ s 2 +.4s+ s 2 +.4s+ ]. We choose conservaive L ρ = 2. In Fig. 2, we plo Gs L L ρ as a funcion of ωk and compare i o. We noice ha for ωk > 3, we have 43 Auhoried licensed use limied o: UNIVERSITY OF CONNECTICUT. Downloaded on February 2, 29 a 9:44 from IEEE Xplore. Resricions apply.

6 Time a solid, ˆ dashed, and rdoed unmodelled dynamics. The conrol signal and he sysem response approimae he same signals of a closed-loop reference sysem, which can be designed o achieve desired specificaions. REFERENCES [] C. Cao and N. Hovakimyan. L adapive conroller for a class of sysems wih unknown nonlineariies: Par I. Submied o American Conrol Conference, 28. [2] H. K. Khalil. Nonlinear Sysems. Prenice Hall, Englewood Cliffs, NJ, 22. [3] C. Cao and N. Hovakimyan. Guaraneed ransien performance wih L adapive conroller for sysems wih unknown ime-varying parameers: Par I. American Conrol Conference, pages , ime b Time-hisory of u Fig Performance of L adapive conroller for σ = sin Time a solid, ˆ dashed, and rdoed ime b Time-hisory of u Fig. 5. Performance of L adapive conroller for σ = sin sin5. Gs L L ρ <. Since ω >.5, we se k = 6. We se he adapive gain Γ c =. The simulaion resuls of he L adapive conroller are shown in Figures 3a-3b for he reference inpu r = cos.5. Ne, we consider he same sysem in he presence of σ = sin The simulaion resuls, wihou any reuning of he conroller, are shown in 4a-4b. Finally, we change σ = sin sin5. The simulaion resuls are shown in 5a-5b. We noe ha he L adapive conroller guaranees smooh and uniform ransien performance in he presence of differen unknown nonlineariies and does no require any reuning. We also noice ha and ˆ are almos he same in Figs. 3a, 4a and 5a. VI. CONCLUSION A novel L adapive conrol archiecure is presened ha has guaraneed ransien response in addiion o sable racking for uncerain nonlinear sysems in he presence of 44 Auhoried licensed use limied o: UNIVERSITY OF CONNECTICUT. Downloaded on February 2, 29 a 9:44 from IEEE Xplore. Resricions apply.

L 1 Adaptive Controller for a Class of Systems with Unknown

L 1 Adaptive Controller for a Class of Systems with Unknown 28 American Control Conference Westin Seattle Hotel, Seattle, Washington, USA June -3, 28 FrA4.2 L Adaptive Controller for a Class of Systems with Unknown Nonlinearities: Part I Chengyu Cao and Naira Hovakimyan

More information

Module 4: Time Response of discrete time systems Lecture Note 2

Module 4: Time Response of discrete time systems Lecture Note 2 Module 4: Time Response of discree ime sysems Lecure Noe 2 1 Prooype second order sysem The sudy of a second order sysem is imporan because many higher order sysem can be approimaed by a second order model

More information

4. Advanced Stability Theory

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

More information

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

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

More information

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

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

More information

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

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

More information

Lecture 20: Riccati Equations and Least Squares Feedback Control

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

More information

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

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

More information

10. State Space Methods

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

More information

Stability and Bifurcation in a Neural Network Model with Two Delays

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

More information

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

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

More information

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

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

More information

15. Vector Valued Functions

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

More information

arxiv: v1 [math.gm] 4 Nov 2018

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

More information

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

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

More information

Sliding Mode Extremum Seeking Control for Linear Quadratic Dynamic Game

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

More information

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

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

More information

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

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

More information

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

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

More information

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

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

More information

LAPLACE TRANSFORM AND TRANSFER FUNCTION

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

More information

Haar Wavelet Operational Matrix Method for Solving Fractional Partial Differential Equations

Haar Wavelet Operational Matrix Method for Solving Fractional Partial Differential Equations Copyrigh 22 Tech Science Press CMES, vol.88, no.3, pp.229-243, 22 Haar Wavele Operaional Mari Mehod for Solving Fracional Parial Differenial Equaions Mingu Yi and Yiming Chen Absrac: In his paper, Haar

More information

Settling Time Design and Parameter Tuning Methods for Finite-Time P-PI Control

Settling Time Design and Parameter Tuning Methods for Finite-Time P-PI Control Journal of Conrol Science an Engineering (6) - oi:.765/8-/6.. D DAVID PUBLISHING Seling ime Design an Parameer uning Mehos for Finie-ime P-PI Conrol Keigo Hiruma, Hisaazu Naamura an Yasuyui Saoh. Deparmen

More information

Chapter 4. Truncation Errors

Chapter 4. Truncation Errors Chaper 4. Truncaion Errors and he Taylor Series Truncaion Errors and he Taylor Series Non-elemenary funcions such as rigonomeric, eponenial, and ohers are epressed in an approimae fashion using Taylor

More information

Week 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem)

Week 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem) Week 1 Lecure Problems, 5 Wha if somehing oscillaes wih no obvious spring? Wha is ω? (problem se problem) Sar wih Try and ge o SHM form E. Full beer can in lake, oscillaing F = m & = ge rearrange: F =

More information

L 1 Adaptive Output Feedback Controller to Systems of Unknown

L 1 Adaptive Output Feedback Controller to Systems of Unknown Proceedings of the 27 American Control Conference Marriott Marquis Hotel at Times Square New York City, USA, July 11-13, 27 WeB1.1 L 1 Adaptive Output Feedback Controller to Systems of Unknown Dimension

More information

1 Solutions to selected problems

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

More information

Designing Information Devices and Systems I Spring 2019 Lecture Notes Note 17

Designing Information Devices and Systems I Spring 2019 Lecture Notes Note 17 EES 16A Designing Informaion Devices and Sysems I Spring 019 Lecure Noes Noe 17 17.1 apaciive ouchscreen In he las noe, we saw ha a capacior consiss of wo pieces on conducive maerial separaed by a nonconducive

More information

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

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

More information

Fractional Method of Characteristics for Fractional Partial Differential Equations

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

More information

Chapter 2. First Order Scalar Equations

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

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

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

More information

Chapter 8 The Complete Response of RL and RC Circuits

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

More information

An Introduction to Malliavin calculus and its applications

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

More information

Anti-Disturbance Control for Multiple Disturbances

Anti-Disturbance Control for Multiple Disturbances Workshop a 3 ACC Ani-Disurbance Conrol for Muliple Disurbances Lei Guo (lguo@buaa.edu.cn) Naional Key Laboraory on Science and Technology on Aircraf Conrol, Beihang Universiy, Beijing, 9, P.R. China. Presened

More information

Linear Control System EE 711. Design. Lecture 8 Dr. Mostafa Abdel-geliel

Linear Control System EE 711. Design. Lecture 8 Dr. Mostafa Abdel-geliel Linear Conrol Sysem EE 7 MIMO Sae Space Analysis and Design Lecure 8 Dr. Mosafa Abdel-geliel Course Conens Review Sae Space SS modeling and analysis Sae feed back design Oupu feedback design Observer design

More information

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

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

More information

Online Appendix to Solution Methods for Models with Rare Disasters

Online Appendix to Solution Methods for Models with Rare Disasters Online Appendix o Soluion Mehods for Models wih Rare Disasers Jesús Fernández-Villaverde and Oren Levinal In his Online Appendix, we presen he Euler condiions of he model, we develop he pricing Calvo block,

More information

, the. L and the L. x x. max. i n. It is easy to show that these two norms satisfy the following relation: x x n x = (17.3) max

, the. L and the L. x x. max. i n. It is easy to show that these two norms satisfy the following relation: x x n x = (17.3) max ecure 8 7. Sabiliy Analyi For an n dimenional vecor R n, he and he vecor norm are defined a: = T = i n i (7.) I i eay o how ha hee wo norm aify he following relaion: n (7.) If a vecor i ime-dependen, hen

More information

Math 334 Fall 2011 Homework 11 Solutions

Math 334 Fall 2011 Homework 11 Solutions Dec. 2, 2 Mah 334 Fall 2 Homework Soluions Basic Problem. Transform he following iniial value problem ino an iniial value problem for a sysem: u + p()u + q() u g(), u() u, u () v. () Soluion. Le v u. Then

More information

SOLUTIONS TO ECE 3084

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

More information

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

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

More information

L p -L q -Time decay estimate for solution of the Cauchy problem for hyperbolic partial differential equations of linear thermoelasticity

L p -L q -Time decay estimate for solution of the Cauchy problem for hyperbolic partial differential equations of linear thermoelasticity ANNALES POLONICI MATHEMATICI LIV.2 99) L p -L q -Time decay esimae for soluion of he Cauchy problem for hyperbolic parial differenial equaions of linear hermoelasiciy by Jerzy Gawinecki Warszawa) Absrac.

More information

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

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

More information

Differential Equations

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

More information

Laplace Transforms. Examples. Is this equation differential? y 2 2y + 1 = 0, y 2 2y + 1 = 0, (y ) 2 2y + 1 = cos x,

Laplace Transforms. Examples. Is this equation differential? y 2 2y + 1 = 0, y 2 2y + 1 = 0, (y ) 2 2y + 1 = cos x, Laplace Transforms Definiion. An ordinary differenial equaion is an equaion ha conains one or several derivaives of an unknown funcion which we call y and which we wan o deermine from he equaion. The equaion

More information

Problemas das Aulas Práticas

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

More information

Let us start with a two dimensional case. We consider a vector ( x,

Let us start with a two dimensional case. We consider a vector ( x, Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our

More information

Computation of the Effect of Space Harmonics on Starting Process of Induction Motors Using TSFEM

Computation of the Effect of Space Harmonics on Starting Process of Induction Motors Using TSFEM Journal of elecrical sysems Special Issue N 01 : November 2009 pp: 48-52 Compuaion of he Effec of Space Harmonics on Saring Process of Inducion Moors Using TSFEM Youcef Ouazir USTHB Laboraoire des sysèmes

More information

Theory of! Partial Differential Equations!

Theory of! Partial Differential Equations! hp://www.nd.edu/~gryggva/cfd-course/! Ouline! Theory o! Parial Dierenial Equaions! Gréar Tryggvason! Spring 011! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and

More information

ESTIMATION OF AMPLITUDE OF OSCILLATIONS IN SLIDING MODE SYSTEMS CAUSED BY TIME DELAY

ESTIMATION OF AMPLITUDE OF OSCILLATIONS IN SLIDING MODE SYSTEMS CAUSED BY TIME DELAY Asian Journal of Conrol, Vol 6, No, pp 57-53, December 57 -Brief Paper- ESTIMATION OF AMPLITUDE OF OSCILLATIONS IN SLIDING MODE SYSTEMS CAUSED BY TIME DELAY P Acosa, A Polyakov, L Fridman, and V Srygin

More information

Some Basic Information about M-S-D Systems

Some Basic Information about M-S-D Systems Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,

More information

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

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

More information

Homotopy Perturbation Method for Solving Some Initial Boundary Value Problems with Non Local Conditions

Homotopy Perturbation Method for Solving Some Initial Boundary Value Problems with Non Local Conditions Proceedings of he World Congress on Engineering and Compuer Science 23 Vol I WCECS 23, 23-25 Ocober, 23, San Francisco, USA Homoopy Perurbaion Mehod for Solving Some Iniial Boundary Value Problems wih

More information

TO our knowledge, most exciting results on the existence

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

More information

On Oscillation of a Generalized Logistic Equation with Several Delays

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

More information

arxiv: v1 [math.pr] 19 Feb 2011

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

More information

Modal identification of structures from roving input data by means of maximum likelihood estimation of the state space model

Modal identification of structures from roving input data by means of maximum likelihood estimation of the state space model Modal idenificaion of srucures from roving inpu daa by means of maximum likelihood esimaion of he sae space model J. Cara, J. Juan, E. Alarcón Absrac The usual way o perform a forced vibraion es is o fix

More information

Asymptotic instability of nonlinear differential equations

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

More information

Challenge Problems. DIS 203 and 210. March 6, (e 2) k. k(k + 2). k=1. f(x) = k(k + 2) = 1 x k

Challenge Problems. DIS 203 and 210. March 6, (e 2) k. k(k + 2). k=1. f(x) = k(k + 2) = 1 x k Challenge Problems DIS 03 and 0 March 6, 05 Choose one of he following problems, and work on i in your group. Your goal is o convince me ha your answer is correc. Even if your answer isn compleely correc,

More information

Chapter 7 Response of First-order RL and RC Circuits

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

More information

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

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

More information

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

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

More information

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes Half-Range Series 2.5 Inroducion In his Secion we address he following problem: Can we find a Fourier series expansion of a funcion defined over a finie inerval? Of course we recognise ha such a funcion

More information

ENGI 9420 Engineering Analysis Assignment 2 Solutions

ENGI 9420 Engineering Analysis Assignment 2 Solutions ENGI 940 Engineering Analysis Assignmen Soluions 0 Fall [Second order ODEs, Laplace ransforms; Secions.0-.09]. Use Laplace ransforms o solve he iniial value problem [0] dy y, y( 0) 4 d + [This was Quesion

More information

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux

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

More information

Sliding Mode Controller for Unstable Systems

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

More information

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

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

More information

Optimal Path Planning for Flexible Redundant Robot Manipulators

Optimal Path Planning for Flexible Redundant Robot Manipulators 25 WSEAS In. Conf. on DYNAMICAL SYSEMS and CONROL, Venice, Ialy, November 2-4, 25 (pp363-368) Opimal Pah Planning for Flexible Redundan Robo Manipulaors H. HOMAEI, M. KESHMIRI Deparmen of Mechanical Engineering

More information

Two Coupled Oscillators / Normal Modes

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

More information

Expert Advice for Amateurs

Expert Advice for Amateurs Exper Advice for Amaeurs Ernes K. Lai Online Appendix - Exisence of Equilibria The analysis in his secion is performed under more general payoff funcions. Wihou aking an explici form, he payoffs of he

More information

Theory of! Partial Differential Equations-I!

Theory of! Partial Differential Equations-I! hp://users.wpi.edu/~grear/me61.hml! Ouline! Theory o! Parial Dierenial Equaions-I! Gréar Tryggvason! Spring 010! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and

More information

The bang-bang funnel controller: time delays and case study

The bang-bang funnel controller: time delays and case study 213 European Conrol Conference ECC) July 17-19, 213, Zürich, Swizerland. The bang-bang funnel conroller: ime delays and case sudy Daniel Liberzon and Sephan Trenn Absrac We invesigae he recenly inroduced

More information

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

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

More information

Dynamic Effects of Feedback Control!

Dynamic Effects of Feedback Control! Dynamic Effecs of Feedback Conrol! Rober Sengel! Roboics and Inelligen Sysems MAE 345, Princeon Universiy, 2017 Inner, Middle, and Ouer Feedback Conrol Loops Sep Response of Linear, Time- Invarian (LTI)

More information

Multi-Layer Switching Control

Multi-Layer Switching Control 5 American Conrol Conference June 8-1, 5. Porland, OR, USA FrC7.3 Muli-Layer Swiching Conrol Idin Karuei, Nader Meskin, Amir G. Aghdam Deparmen of Elecrical and Compuer Engineering, Concordia Universiy

More information

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

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

More information

L 1 Adaptive Controller for Multi Input Multi Output Systems in the Presence of Unmatched Disturbances

L 1 Adaptive Controller for Multi Input Multi Output Systems in the Presence of Unmatched Disturbances 28 American Control Conference Westin Seattle Hotel, Seattle, Washington, USA June -3, 28 FrA4.4 L Adaptive Controller for Multi Input Multi Output Systems in the Presence of Unmatched Disturbances Chengyu

More information

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

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

More information

Optimal Control of Dc Motor Using Performance Index of Energy

Optimal Control of Dc Motor Using Performance Index of Energy American Journal of Engineering esearch AJE 06 American Journal of Engineering esearch AJE e-issn: 30-0847 p-issn : 30-0936 Volume-5, Issue-, pp-57-6 www.ajer.org esearch Paper Open Access Opimal Conrol

More information

R.#W.#Erickson# Department#of#Electrical,#Computer,#and#Energy#Engineering# University#of#Colorado,#Boulder#

R.#W.#Erickson# Department#of#Electrical,#Computer,#and#Energy#Engineering# University#of#Colorado,#Boulder# .#W.#Erickson# Deparmen#of#Elecrical,#Compuer,#and#Energy#Engineering# Universiy#of#Colorado,#Boulder# Chaper 2 Principles of Seady-Sae Converer Analysis 2.1. Inroducion 2.2. Inducor vol-second balance,

More information

Laplace transfom: t-translation rule , Haynes Miller and Jeremy Orloff

Laplace transfom: t-translation rule , Haynes Miller and Jeremy Orloff Laplace ransfom: -ranslaion rule 8.03, Haynes Miller and Jeremy Orloff Inroducory example Consider he sysem ẋ + 3x = f(, where f is he inpu and x he response. We know is uni impulse response is 0 for

More information

A NEW TECHNOLOGY FOR SOLVING DIFFUSION AND HEAT EQUATIONS

A NEW TECHNOLOGY FOR SOLVING DIFFUSION AND HEAT EQUATIONS THERMAL SCIENCE: Year 7, Vol., No. A, pp. 33-4 33 A NEW TECHNOLOGY FOR SOLVING DIFFUSION AND HEAT EQUATIONS by Xiao-Jun YANG a and Feng GAO a,b * a School of Mechanics and Civil Engineering, China Universiy

More information

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

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

More information

Convergence of the Neumann series in higher norms

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

More information

Numerical Dispersion

Numerical Dispersion eview of Linear Numerical Sabiliy Numerical Dispersion n he previous lecure, we considered he linear numerical sabiliy of boh advecion and diffusion erms when approimaed wih several spaial and emporal

More information

Economics 8105 Macroeconomic Theory Recitation 6

Economics 8105 Macroeconomic Theory Recitation 6 Economics 8105 Macroeconomic Theory Reciaion 6 Conor Ryan Ocober 11h, 2016 Ouline: Opimal Taxaion wih Governmen Invesmen 1 Governmen Expendiure in Producion In hese noes we will examine a model in which

More information

Existence of multiple positive periodic solutions for functional differential equations

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

More information

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

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

More information

arxiv: v1 [math.fa] 3 Jan 2019

arxiv: v1 [math.fa] 3 Jan 2019 DAMPED AND DIVERGENCE EXACT SOLUTIONS FOR THE DUFFING EQUATION USING LEAF FUNCTIONS AND HYPERBOLIC LEAF FUNCTIONS A PREPRINT arxiv:9.66v [mah.fa] Jan 9 Kazunori Shinohara Deparmen of Mechanical Sysems

More information

CONTRIBUTION TO IMPULSIVE EQUATIONS

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

More information

Dynamic Analysis of Loads Moving Over Structures

Dynamic Analysis of Loads Moving Over Structures h Inernaional ongress of roaian ociey of echanics epember, 18-, 3 Bizovac, roaia ynamic nalysis of Loads oving Over rucures Ivica Kožar, Ivana Šimac Keywords: moving load, direc acceleraion mehod 1. Inroducion

More information

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

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

More information

The Contradiction within Equations of Motion with Constant Acceleration

The Contradiction within Equations of Motion with Constant Acceleration The Conradicion wihin Equaions of Moion wih Consan Acceleraion Louai Hassan Elzein Basheir (Daed: July 7, 0 This paper is prepared o demonsrae he violaion of rules of mahemaics in he algebraic derivaion

More information

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

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

More information

Econ107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8)

Econ107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8) I. Definiions and Problems A. Perfec Mulicollineariy Econ7 Applied Economerics Topic 7: Mulicollineariy (Sudenmund, Chaper 8) Definiion: Perfec mulicollineariy exiss in a following K-variable regression

More information

POSITIVE AND MONOTONE SYSTEMS IN A PARTIALLY ORDERED SPACE

POSITIVE AND MONOTONE SYSTEMS IN A PARTIALLY ORDERED SPACE Urainian Mahemaical Journal, Vol. 55, No. 2, 2003 POSITIVE AND MONOTONE SYSTEMS IN A PARTIALLY ORDERED SPACE A. G. Mazo UDC 517.983.27 We invesigae properies of posiive and monoone differenial sysems wih

More information

Multi-scale 2D acoustic full waveform inversion with high frequency impulsive source

Multi-scale 2D acoustic full waveform inversion with high frequency impulsive source Muli-scale D acousic full waveform inversion wih high frequency impulsive source Vladimir N Zubov*, Universiy of Calgary, Calgary AB vzubov@ucalgaryca and Michael P Lamoureux, Universiy of Calgary, Calgary

More information

2. Nonlinear Conservation Law Equations

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

More information

A Uniform Asymptotic Estimate for Discounted Aggregate Claims with Subexponential Tails

A Uniform Asymptotic Estimate for Discounted Aggregate Claims with Subexponential Tails A Uniform Asympoic Esimae for Discouned Aggregae Claims wih Subeponenial Tails Xuemiao Hao and Qihe Tang Deparmen of Saisics and Acuarial Science The Universiy of Iowa 241 Schaeffer Hall, Iowa Ciy, IA

More information