Nested Saturation with Guaranteed Real Poles 1
|
|
- Kelley Hood
- 5 years ago
- Views:
Transcription
1 Neste Saturation with Guarantee Real Poles Eric N Johnson 2 an Suresh K Kannan 3 School of Aerospace Engineering Georgia Institute of Technology, Atlanta, GA 3332 Abstract The global stabilization of asymptotically nullcontrollable linear systems with boune controls has been stuie extensively An early contribution was by Teel [6] who propose a set of neste saturators to globally asymptotically stabilize the special case of n- integrators with one input Using this law however, the close loop system pole locations epen on the choice of coorinate transformation use to arrive at the control law In this paper we suggest an approach that allows the esigner to pick transformations that facilitate the placement of the close loop poles on the negative real axis Introuction The problem aresse involves the global stabilization of a chain of integrators ẋ = x 2,, ẋ n = u () The system given by () is a subset of a class of systems that are sai to be asymptotically null-controllable with boune controls [3, ] This property was shown in [2] to be equivalent to the system being stabilizable an having all open-loop poles in the close left-half plane It was shown in [4] that it is not possible to globally stabilize integrator chains of orer n > 2 using a boune linear feeback law However, it was shown by Teel in [6] that a nonlinear law consisting of neste saturators can guarantee global asymptotic stability for integrator chains of any orer n This control law may be expresse as u = σ n (h n (x) + σ n (h n (x) + + σ (h (x)))) where h i are linear combinations of the state (feeback) an the saturation functions σ i satisfy certain properties The existence of such a globally stabilizing control law was establishe in [6] by choosing one set of h i s This work was supporte in part by DARPA contract # C-34 2 Lockhee Martin Assistant Professor of Avionics Integration ericjohnson@aegatecheu 3 Grauate Research Assistant suresh kannan@aegatecheu /3/$7 23 IEEE 497 such that global asymptotic stability coul be proven The choice of h i is a esign egree of freeom an may be exercise to prescribe pole locations an the linear ynamics when ifferent elements of the control law are saturate We observe that the h i chosen by Teel with conventional saturation functions (see Definition 2) results in all the poles of the close loop system resiing at when none of the saturation elements in the control law are saturate If the k th saturator is the outermost element to be saturate, then the resulting close loop system has poles at with multiplicity n k an poles at with multiplicity k, at least until the element comes out of saturation A iscussion on the prescription of performance by pole placement (both real an complex) is provie in [5], however no explicit transformation is provie Another aspect is the behavior of these poles as ifferent elements of the control law saturate Ieally, these poles shoul not change when saturation occurs Both these properties (pole placement an movement when saturate) are useful if the neste saturation control law is to be employe in practice We believe that the simple an elegant neste saturation law can benefit greatly from these properties Hence, the effort here is to evelop a transformation, ie, a way to select h i such that close loop poles for the unsaturate system may be prescribe as { a, a 2, a n }, where a i R \ an a i > for stability Aitionally, it will be shown that when the outermost saturate element is σ k, the poles of resulting linear system resie at { a, a 2, a n k,, 2, k } 2 Main Result Definition (Linear saturation) Define constants (L, M) R+ such that < L M Now, efine a function σ : R R σ is sai to be a linear saturation if it is continuous, nonecreasing an satisfies a sσ(s) > s b σ(s) = s when s L c σ(s) M s R Denver, Colorao June 4-6, 23
2 Definition 2 (Conventional saturation) σ : R R is sai to be a conventional saturation if it has a limit M R+ such that a sσ(s) > s b σ(s) = s when s M c σ(s) = M when s > M Remark σ is sai to be saturate when its argument is not in its linear region For linear saturation this occurs when s > L For conventional saturation this occurs when s > M Remark 2 Conventional saturation is a special case of linear saturation with L = M an a constant saturation value M Lemma Consier a chain of n-integrators, given by (), which may be represente as ẋ = A x x + B x u, with x R n, u R an A x =, B x = (2) then there exists a linear transformation y = T yx x which transforms () into ẏ = A y y + B y u where, a n a n a n a n A y =, B y = a 2 an the elements a i R \ with i = n Proof: Given a set of coefficients a n a n a (3) A = {a, a 2,, a n } (4) Let A l A represent a subset containing the first l elements of A Define a function Fk m(a l) which acts over the set A l Fk m is use to generate the prouct of combinations of elements taken m at a time from A l The number of ( such) combinations is given by the binomial l coefficient Hence, F m k m(a l) may be treate as a generating function that outputs the k th combination of the prouct of m elements taken from the set A l without repetition an isregaring orer Note that Fk = 498 In orer to generate the transformation T yx, efine the function C(l, m), with l [,, n], m [,, l] an m l, over the set of coefficients A given by (4) C(l, m) = C l m k= F m k (A l ) (5) C(l, ) = (6) ( l C m l is the binomial coefficient The new coori- m) nate system is characterize by y n i = a i+ i C(i, j)x n j, i [,, n ] an the transformation T yx is explicitly given by T yx(n i)(n j) = a i+ C(i, j) i j (7) T yx(n i)(n j) = i < j (8) for i, j [,, n ] Aitionally, T yx is an upper iagonal matrix with non-zero iagonal entries Hence, T xy = Tyx exists Finally, observing that ẋ = T yx A x T yx y + T yx B x u = A y y + B y u it is enough to verify that A y T yx = T yx A x an that T yx B x = B y This may be carrie out using Equations 2, 3 an 8 Theorem For the system given by () Given any set of positive constants {(L i, M i )}, where L i M i for i =,, n an M i < 2 L i+ for i =,, n, an for any set of functions {σ i } that are linear saturations for {(L i, M i )}, there exists a linear coorinate transformation y = T yx x such that the boune control u = σ n (y n + σ n (y n + + σ (y ))) (9) results in a globally asymptotically stable system Proof: In short, use the transformation given by Lemma in the proof of Theorem 2 in [6] It is however restate here for completeness Use the coorinate transformation y = T yx x given by Lemma an choose the set of coefficients a i > Substituting the neste saturation law given by Eq (9) into Eq () an expaning yiels the close loop system ẏ = a n [y 2 + +y n σ n (y n + σ n ( σ (y )))] ẏ 2 = a n [y 3 + +y n σ n (y n + σ n ( σ (y )))] () ẏ n = a 2 [ y n σ n (y n + σ n ( σ (y )))] ẏ n = a σ n (y n + σ n ( σ (y ))) Denver, Colorao June 4-6, 23
3 The trajectory of y n is examine first Choosing a Lyapunov function V n = y 2 n, with y n R Its erivative V n may be written as V n = 2a y n [σ n (y n + σ n (y n + + σ (y )))] Noting that a i > Definition, conitions (a), (b), imply that y n an σ n ( ) are the same sign only if y n + σ n ( ) is the same sign as y n Conition (c) of Definition applie to σ n an having chosen M n < 2 L n, it can be seen that V n < for all y n / Q n = {y n : y n 2 L n} If starting outsie Q n, the trajectory of y n eventually enters Q n in finite time Since the RHS of Eq () is globally Lipschitz, the erivatives are boune resulting in the remaining states y y n remaining boune for any given finite time Once y n has entere Q n, conition (b) of Definition implies σ n operates in its linear region because the argument to σ n is boune as y n + σ n ( ) 2 L n + M n L n The equation for the evolution of y n is now given by ẏ n = a 2 y n a 2 y n a 2 σ n (y n + + σ (y )) = a 2 σ n (y n + + σ (y )) which is similar to the expression for ẏ n Using similar arguments as that use for the evolution of y n, it can be shown that y n enters a set Q n in finite time an remains in Q n thereafter with all remaining states being boune Continuing in the same fashion, it can be shown that every state y i for i [,, n], enters a set Q i = {y i : y i 2 L i} in finite time an all saturation functions σ i are operating in their linear regions Hence after a certain finite amount of time the governing equations, Eq (), becomes ẏ = a n y ẏ 2 = a n (y + y 2 ) ẏ n = a (y + y y n ) which is exponentially stable Corollary (Pole location) If the saturators use are Conventional saturation, an none of the σ i are saturate, the poles of the linearize close loop system resie at { a, a 2,, a n } During perios when the outermost saturate element is the k th saturator, σ k, the poles of the resulting close loop linear system resie at { a, a 2,, a n k,, 2,, k } 499 Proof: Using the neste saturation law, the closeloop n-integrator system may be expresse as ẋ n + σ n (y n + σ n (y n + + σ (y ))) = When the k th saturator is saturate, an σ k+ σ n are not saturate, the close loop system is given by ẋ n + y n + y n + y k+ ± M k = This represents a force linear system where the forcing function is the constant M k Examining the homogeneous part = ẋ n + y n + y n + + y k+ Using Eq (7) to expan y i = ẋ n + a C(, j)x n j + a 2 C(, j)x n j + n (k+) + a n k C (n (k + ), j) x n j Noting that x = x, ẋ = x 2,, x (n ) = x n, x (n) = ẋ n, an substituting p = n k for clarity the characteristic equation may be written as Υ(λ) = λ n + a C(, )λ n + a 2 C(, )λ n + a 2 C(, )λ n 2 + a p C(p, )λ n + + a p C(p, p )λ k Factoring out λ k Υ(λ) =λ k [λ p + a C(, )λ p + a 2 C(, )λ p + a 2 C(, )λ p 2 + a p C(p, )λ p + + a p C(p, p )] an may be written in its final form as Υ(λ) = λ k (λ + a )(λ + a 2 )(λ + a p ) which has k zeros an p = n k non-zero stable poles at known locations Corollary 2 During perios when σ k is the outermost saturate element in the control law of Theorem an the coorinate transformation use is given Denver, Colorao June 4-6, 23
4 by Lemma, then, in steay-state, the magnitue of the k th erivative, ẋ k, is given by ẋ k = M k a n k C(n (k + ), n (k + )) () for k [,, n ] for k = n ẋ k = M k (2) Proof: If σ k is saturate, the close loop system may be written as Using Eq (7) ẋ n + y n + y n + y k+ ± M k = = ẋ n + a C(, j)x n j + a 2 C(, j)x n j + (3) n (k+) + a n k C (n (k + ), j) x n j ± M k When the outermost saturate element is σ k, the ynamics eventually reach a saturate-equilibrium region where higher-orer erivatives reach zero So, x (n) x (k+), ie, ẋ n, x n x k+2 go to zero The only term left from Eq (3) is a n k C(n (k + ), n (k + ))x k+ ± M k = (4) Noting that ẋ k = x k+, rearranging Eq (4) an taking the absolute value of both sies results in Eq () Finally, when k = n, the outermost saturator σ n is saturate an Eq (3) reuces to If x (n) (t) L n+ ɛ for all t t an for some ɛ > an given linear saturation functions σ i with parameters (L i, M i ) satisfying, L i M i i =,, n + M i < 2 L i+ i =,, n M n ɛ then, the feeback u = x (n) σ n (y n + σ n (y n + + σ (y ))) with y = T yx e given by Lemma, where, e i = x i x (i ) for i = n, results in a globally asymptotically stable system Aitionally if conventional saturation elements are use, the error ynamics are governe by Corollary an quasi-steay rates governe by Corollary 2 Proof: The ynamics of Eq (6) may be expresse in terms of the error e ė = e 2,, ė n = x (n) + σ n+ (u) With the given control law, if the magnitue of the n th erivative of the comman x is always such that x (n) (t) L n+ ɛ for all t t an M n ɛ, then the magnitue of the argument of σ n+ is x (n) σ n ( ) L n+ an σ n+ is always in its linear region, resulting in the close loop error ynamics becoming ė = e 2,, ė n = σ n (y n + σ n (y n + + σ (y ))) (8) The conitions of this corollary an form of Eq (8) satisfy the requirements of Theorem This implies that the ynamics of e are asymptotically stable an hence x tracks x asymptotically The form of Eq (8) also allows Corollary an Corollary 2 to be applie irectly ẋ n ± M n = (5) Rearranging Eq (5) an taking magnitues of both sies results in Eq (2) Corollary 3 (Restricte Tracking) Consier a nonlinear system with magnitue saturation at the input u given by ẋ = x 2,, ẋ n = σ n+ (u) (6) an a compatible reference signal given by [ x (t), ẋ (t), x (n) (t) ] (7) 5 3 Examples Global Stabilization: Consier the problem of stabilizing the 3 r orer system ẋ = x 2, ẋ 2 = x 3, ẋ 3 = u using boune control u [, ] (conventional saturation) with poles at {, 3, 2} Then, {a, a 2, a 3 } = {, 3, 2} The transformation require to achieve these poles may be expresse as y y 2 = a 3(a a 2 ) a 3 (a + a 2 ) a 3 a 2 (a ) a 2 x y 3 x 2 a x 3 Denver, Colorao June 4-6, 23
5 σ x States Saturation Functions time /s Figure : Initial conition response of a 3 r orer system Using the neste saturation law given by Theorem an choosing the saturation element parameters as follows M 3 =, L 3 = M 3 M 2 = 2 L 3 ɛ, L 2 = M 2 M = 2 L 2 ɛ, L = M where ɛ is a small positive number, that is use to satisfy the inequality M i < 2 L i+ Aitionally the saturation element parameters are chosen L i = M i (Conventional saturation) Then, the close loop system is given by ẋ 3 + σ 3 (y 3 + σ 2 (y 2 + σ (y ))) = An initial conition response with x = [, 5, ] is shown in Figure The figure also shows the outputs of the saturation elements - 5s, σ 3 is saturate 5-3s, σ 2 is saturate 3-46s, σ is saturate 46-5s, control law is unsaturate The only region where the system practically reaches a saturate-equilibrium is when σ is saturate, between an 4 secons The equilibrium value for ẋ is given by Corollary 2 ẋ = x 2 = M a 2 a = 833 an matches the slope of x in Figure x x 2 x 3 σ σ 2 σ 3 5 Restricte Tracking: Consier a chain of 4 integrators where, σ 5 (u) represents a magnitue saturate actuator ẋ = x 2, ẋ 2 = x 3, ẋ 3 = x 4, ẋ 4 = σ 5 (u) σ 5 is a conventional saturation function with parameters (L 5, M 5 ) A compatible comman may be represente as [x, ẋ, ẍ, x, x ] Defining the error as, e = x x, the error erivatives may be written as ė = x 2 ẋ ë = x 3 ẍ e = x 4 x e = σ 5 (u) x The control is given by Corollary 3 u = x σ 4 (y 4 + σ 3 (y 3 + σ 2 (y 2 + σ (y )))) with x L 5 ɛ an M 4 ɛ, for some ɛ >, the saturation function parameters (L i, M i ) chosen to satisfy the conitions given by Corollary 3 an y i given by Lemma The coorinate transformation use is y = T yx e where T yx is given by a a 2 a 3 a 4 (a a 3 + a 2 a 3 + a a 2 )a 4 (a + a 2 + a 3 )a 4 a 4 a a 2 a 3 (a + a 2 )a 3 a 3 a a 2 a 2 a Here, the poles were taken to be at { a, a 2, a 3, a 4 } = { 5,, 2, 3} The saturation function parameters were chosen as L 5 = ɛ = 2 L 5 L 4 = ɛ L 3 = 2 L 4 L 2 = 2 L 3 L = 2 L 2 M 5 = L 5 ɛ M 4 = ɛ ɛ M 3 = L 3 ɛ M 2 = L 2 ɛ M = L ɛ where ɛ is a small positive number chosen to satisfy M i < L i If Corollary 2 is evaluate for various saturation elements being saturate ė 4 = M 4 when, σ 4 is saturate ė 3 = M 3 a when, σ 3 is saturate ė 2 = M 2 a a 2 when, σ 2 is saturate ė = M a a 2 a 3 when, σ is saturate (9) Denver, Colorao June 4-6, 23
6 The response of this system to a sinusoial comman compatible with x = 5 sin(5t) an zero initial conitions is illustrate in Figure 2 From Eq (9) notice that as the banwith ie, a i is increase, the error rates in saturate-equilibrium ecrease Hence for higher banwith, the overall settling time can be higher, which is perhaps counter-intuitive This aspect is further illustrate in Figure 3 where it is observe that the control law with faster poles (all at -5) takes longer to be regulate back to than the system with slower poles (all at -5) The initial conition use was x = [,,, 2] T with zero comman x, x x vs x Saturation Functions x x 5 σ σ 2 4 Conclusion The extensions presente thus far provies a transformation that allows placement of poles on the non-zero real axis Assuming the poles chosen are stable, global asymptotic stability is guarantee These chosen poles are guarantee to remain constant, apart from going to zero when the respective saturation element saturates Finally, when in saturate-equilibrium ue to σ k saturating, the quasi-steay rate of state change is given by a Corollary σ time /s Figure 2: Response to a sinusoial comman for a 4 th orer system σ 3 σ 4 References [] Tingshu Hu an Zongli Lin Control Systems with Actuator Saturation : Analysis an Design Birkhaüser, 2 [2] Euaro D Sontag An algebraic approach to boune controllability of linear systems 39:8 88, 984 [3] Hector J Sussmann, Euaro D Sontag, an Yui Yang A general result on the stabilization of linear systems using boune controls IEEE Transactions on Automatic Control, 39(2), ec 994 [4] Hector J Sussmann an Yui Yang On the stabilizability of multiple integrators by means of boune feeback controls In Proceeings of the 3th IEEE Conference on Decision an Control, Brighton, UK, 99 [5] Anrew R Teel Feeback stabilization : Nonlinear solutions to inherently nonlinear problems Technical Report UCB/ERL M92/65, Electronics Research Laboratory, University of California, Berkeley, CA 9472, June 992 [6] Anrew R Teel Global stabilization an restricte tracking for multiple integrators with boune controls Systems & Control Letters, 8:65 7, 992 x x reponse for fast an slow poles Poles at 5 Poles at time /s Figure 3: Comparison of the initial conition response, for a 4 th orer system The soli curve settles faster an has all poles at -5 whilst the ashe-curve settles slower an has poles at Denver, Colorao June 4-6, 23
Optimal Variable-Structure Control Tracking of Spacecraft Maneuvers
Optimal Variable-Structure Control racking of Spacecraft Maneuvers John L. Crassiis 1 Srinivas R. Vaali F. Lanis Markley 3 Introuction In recent years, much effort has been evote to the close-loop esign
More informationExperimental Robustness Study of a Second-Order Sliding Mode Controller
Experimental Robustness Stuy of a Secon-Orer Sliing Moe Controller Anré Blom, Bram e Jager Einhoven University of Technology Department of Mechanical Engineering P.O. Box 513, 5600 MB Einhoven, The Netherlans
More informationEE 370L Controls Laboratory. Laboratory Exercise #7 Root Locus. Department of Electrical and Computer Engineering University of Nevada, at Las Vegas
EE 370L Controls Laboratory Laboratory Exercise #7 Root Locus Department of Electrical an Computer Engineering University of Nevaa, at Las Vegas 1. Learning Objectives To emonstrate the concept of error
More informationExponential Tracking Control of Nonlinear Systems with Actuator Nonlinearity
Preprints of the 9th Worl Congress The International Feeration of Automatic Control Cape Town, South Africa. August -9, Exponential Tracking Control of Nonlinear Systems with Actuator Nonlinearity Zhengqiang
More informationLaplacian Cooperative Attitude Control of Multiple Rigid Bodies
Laplacian Cooperative Attitue Control of Multiple Rigi Boies Dimos V. Dimarogonas, Panagiotis Tsiotras an Kostas J. Kyriakopoulos Abstract Motivate by the fact that linear controllers can stabilize the
More informationAdaptive Control with a Nested Saturation Reference Model
Adaptive Control with a Nested Saturation Reference Model Suresh K Kannan and Eric N Johnson School of Aerospace Engineering Georgia Institute of Technology, Atlanta, GA 3332 This paper introduces a neural
More informationTotal Energy Shaping of a Class of Underactuated Port-Hamiltonian Systems using a New Set of Closed-Loop Potential Shape Variables*
51st IEEE Conference on Decision an Control December 1-13 212. Maui Hawaii USA Total Energy Shaping of a Class of Uneractuate Port-Hamiltonian Systems using a New Set of Close-Loop Potential Shape Variables*
More informationNonlinear Adaptive Ship Course Tracking Control Based on Backstepping and Nussbaum Gain
Nonlinear Aaptive Ship Course Tracking Control Base on Backstepping an Nussbaum Gain Jialu Du, Chen Guo Abstract A nonlinear aaptive controller combining aaptive Backstepping algorithm with Nussbaum gain
More information19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control
19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior
More informationImplicit Differentiation
Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,
More informationELEC3114 Control Systems 1
ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.
More informationVIRTUAL STRUCTURE BASED SPACECRAFT FORMATION CONTROL WITH FORMATION FEEDBACK
AIAA Guiance, Navigation, an Control Conference an Exhibit 5-8 August, Monterey, California AIAA -9 VIRTUAL STRUCTURE BASED SPACECRAT ORMATION CONTROL WITH ORMATION EEDBACK Wei Ren Ranal W. Bear Department
More informationThe derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)
Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)
More informationLinear First-Order Equations
5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)
More informationOn Using Unstable Electrohydraulic Valves for Control
Kailash Krishnaswamy Perry Y. Li Department of Mechanical Engineering, University of Minnesota, 111 Church St. SE, Minneapolis, MN 55455 e-mail: kk,pli @me.umn.eu On Using Unstable Electrohyraulic Valves
More informationLecture 5. Symmetric Shearer s Lemma
Stanfor University Spring 208 Math 233: Non-constructive methos in combinatorics Instructor: Jan Vonrák Lecture ate: January 23, 208 Original scribe: Erik Bates Lecture 5 Symmetric Shearer s Lemma Here
More informationLecture Note 7: Switching Stabilization via Control-Lyapunov Function
ECE7850: Hybrid Systems:Theory and Applications Lecture Note 7: Switching Stabilization via Control-Lyapunov Function Wei Zhang Assistant Professor Department of Electrical and Computer Engineering Ohio
More informationTRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS
TRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS Francesco Bullo Richar M. Murray Control an Dynamical Systems California Institute of Technology Pasaena, CA 91125 Fax : + 1-818-796-8914 email
More informationAdaptive Gain-Scheduled H Control of Linear Parameter-Varying Systems with Time-Delayed Elements
Aaptive Gain-Scheule H Control of Linear Parameter-Varying Systems with ime-delaye Elements Yoshihiko Miyasato he Institute of Statistical Mathematics 4-6-7 Minami-Azabu, Minato-ku, okyo 6-8569, Japan
More informationLeast-Squares Regression on Sparse Spaces
Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction
More informationPhysics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1
Physics 5153 Classical Mechanics The Virial Theorem an The Poisson Bracket 1 Introuction In this lecture we will consier two applications of the Hamiltonian. The first, the Virial Theorem, applies to systems
More informationCalculus of Variations
16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t
More informationLectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs
Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent
More informationGlobal stabilization of feedforward systems with exponentially unstable Jacobian linearization
Global stabilization of feedforward systems with exponentially unstable Jacobian linearization F Grognard, R Sepulchre, G Bastin Center for Systems Engineering and Applied Mechanics Université catholique
More informationConsider for simplicity a 3rd-order IIR filter with a transfer function. where
Basic IIR Digital Filter The causal IIR igital filters we are concerne with in this course are characterie by a real rational transfer function of or, equivalently by a constant coefficient ifference equation
More informationNumerical Integrator. Graphics
1 Introuction CS229 Dynamics Hanout The question of the week is how owe write a ynamic simulator for particles, rigi boies, or an articulate character such as a human figure?" In their SIGGRPH course notes,
More informationBasic IIR Digital Filter Structures
Basic IIR Digital Filter Structures The causal IIR igital filters we are concerne with in this course are characterie by a real rational transfer function of or, equivalently by a constant coefficient
More information1 Lecture 20: Implicit differentiation
Lecture 20: Implicit ifferentiation. Outline The technique of implicit ifferentiation Tangent lines to a circle Derivatives of inverse functions by implicit ifferentiation Examples.2 Implicit ifferentiation
More informationSlide10 Haykin Chapter 14: Neurodynamics (3rd Ed. Chapter 13)
Slie10 Haykin Chapter 14: Neuroynamics (3r E. Chapter 13) CPSC 636-600 Instructor: Yoonsuck Choe Spring 2012 Neural Networks with Temporal Behavior Inclusion of feeback gives temporal characteristics to
More informationAn inductance lookup table application for analysis of reluctance stepper motor model
ARCHIVES OF ELECTRICAL ENGINEERING VOL. 60(), pp. 5- (0) DOI 0.478/ v07-0-000-y An inuctance lookup table application for analysis of reluctance stepper motor moel JAKUB BERNAT, JAKUB KOŁOTA, SŁAWOMIR
More informationProblem Description The problem we consider is stabilization of a single-input multiple-state system with simultaneous magnitude and rate saturations,
SEMI-GLOBAL RESULTS ON STABILIZATION OF LINEAR SYSTEMS WITH INPUT RATE AND MAGNITUDE SATURATIONS Trygve Lauvdal and Thor I. Fossen y Norwegian University of Science and Technology, N-7 Trondheim, NORWAY.
More informationSliding mode approach to congestion control in connection-oriented communication networks
JOURNAL OF APPLIED COMPUTER SCIENCE Vol. xx. No xx (200x), pp. xx-xx Sliing moe approach to congestion control in connection-oriente communication networks Anrzej Bartoszewicz, Justyna Żuk Technical University
More informationMath 1271 Solutions for Fall 2005 Final Exam
Math 7 Solutions for Fall 5 Final Eam ) Since the equation + y = e y cannot be rearrange algebraically in orer to write y as an eplicit function of, we must instea ifferentiate this relation implicitly
More informationMathematical Review Problems
Fall 6 Louis Scuiero Mathematical Review Problems I. Polynomial Equations an Graphs (Barrante--Chap. ). First egree equation an graph y f() x mx b where m is the slope of the line an b is the line's intercept
More informationAnalytic Scaling Formulas for Crossed Laser Acceleration in Vacuum
October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945
More information'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21
Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting
More informationA PAC-Bayesian Approach to Spectrally-Normalized Margin Bounds for Neural Networks
A PAC-Bayesian Approach to Spectrally-Normalize Margin Bouns for Neural Networks Behnam Neyshabur, Srinah Bhojanapalli, Davi McAllester, Nathan Srebro Toyota Technological Institute at Chicago {bneyshabur,
More information6.003 Homework #7 Solutions
6.003 Homework #7 Solutions Problems. Secon-orer systems The impulse response of a secon-orer CT system has the form h(t) = e σt cos(ω t + φ)u(t) where the parameters σ, ω, an φ are relate to the parameters
More informationCalculus of Variations
Calculus of Variations Lagrangian formalism is the main tool of theoretical classical mechanics. Calculus of Variations is a part of Mathematics which Lagrangian formalism is base on. In this section,
More informationNeural Network Controller for Robotic Manipulator
MMAE54 Robotics- Class Project Paper Neural Network Controller for Robotic Manipulator Kai Qian Department of Biomeical Engineering, Illinois Institute of echnology, Chicago, IL 666 USA. Introuction Artificial
More informationSYSTEMS OF DIFFERENTIAL EQUATIONS, EULER S FORMULA. where L is some constant, usually called the Lipschitz constant. An example is
SYSTEMS OF DIFFERENTIAL EQUATIONS, EULER S FORMULA. Uniqueness for solutions of ifferential equations. We consier the system of ifferential equations given by x = v( x), () t with a given initial conition
More informationChapter 2 Lagrangian Modeling
Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie
More informationFORMATION INPUT-TO-STATE STABILITY. Herbert G. Tanner and George J. Pappas
Copyright 2002 IFAC 5th Triennial Worl Congress, Barcelona, Spain FORMATION INPUT-TO-STATE STABILITY Herbert G. Tanner an George J. Pappas Department of Electrical Engineering University of Pennsylvania
More informationA new approach to explicit MPC using self-optimizing control
28 American Control Conference Westin Seattle Hotel, Seattle, Washington, USA June -3, 28 WeA3.2 A new approach to explicit MPC using self-optimizing control Henrik Manum, Sriharakumar Narasimhan an Sigur
More informationFinal Exam Study Guide and Practice Problems Solutions
Final Exam Stuy Guie an Practice Problems Solutions Note: These problems are just some of the types of problems that might appear on the exam. However, to fully prepare for the exam, in aition to making
More informationAdaptive Optimal Path Following for High Wind Flights
Milano (Italy) August - September, 11 Aaptive Optimal Path Following for High Win Flights Ashwini Ratnoo P.B. Sujit Mangal Kothari Postoctoral Fellow, Department of Aerospace Engineering, Technion-Israel
More informationA Review of Feedforward Control Approaches in Nanopositioning for High-Speed SPM
Garrett M. Clayton Department of Mechanical Engineering, Villanova University, Villanova, PA 1985 e-mail: garrett.clayton@villanova.eu Szuchi Tien Department of Mechanical Engineering, National Cheng Kung
More informationA New Backstepping Sliding Mode Guidance Law Considering Control Loop Dynamics
pp. 9-6 A New Backstepping liing Moe Guiance Law Consiering Control Loop Dynamics V. Behnamgol *, A. Vali an A. Mohammai 3, an 3. Department of Control Engineering, Malek Ashtar University of Technology
More informationRobust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k
A Proof of Lemma 2 B Proof of Lemma 3 Proof: Since the support of LL istributions is R, two such istributions are equivalent absolutely continuous with respect to each other an the ivergence is well-efine
More informationLecture 6: Calculus. In Song Kim. September 7, 2011
Lecture 6: Calculus In Song Kim September 7, 20 Introuction to Differential Calculus In our previous lecture we came up with several ways to analyze functions. We saw previously that the slope of a linear
More informationPDE Notes, Lecture #11
PDE Notes, Lecture # from Professor Jalal Shatah s Lectures Febuary 9th, 2009 Sobolev Spaces Recall that for u L loc we can efine the weak erivative Du by Du, φ := udφ φ C0 If v L loc such that Du, φ =
More informationAgmon Kolmogorov Inequalities on l 2 (Z d )
Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,
More informationThe dynamics of the simple pendulum
.,, 9 G. Voyatzis, ept. of Physics, University of hessaloniki he ynamics of the simple penulum Analytic methos of Mechanics + Computations with Mathematica Outline. he mathematical escription of the moel.
More informationTime Headway Requirements for String Stability of Homogeneous Linear Unidirectionally Connected Systems
Joint 48th IEEE Conference on Decision an Control an 8th Chinese Control Conference Shanghai, PR China, December 6-8, 009 WeBIn53 Time Heaway Requirements for String Stability of Homogeneous Linear Uniirectionally
More informationSome Remarks on the Boundedness and Convergence Properties of Smooth Sliding Mode Controllers
International Journal of Automation an Computing 6(2, May 2009, 154-158 DOI: 10.1007/s11633-009-0154-z Some Remarks on the Bouneness an Convergence Properties of Smooth Sliing Moe Controllers Wallace Moreira
More informationLecture 2 Lagrangian formulation of classical mechanics Mechanics
Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,
More informationEquilibrium in Queues Under Unknown Service Times and Service Value
University of Pennsylvania ScholarlyCommons Finance Papers Wharton Faculty Research 1-2014 Equilibrium in Queues Uner Unknown Service Times an Service Value Laurens Debo Senthil K. Veeraraghavan University
More informationRobust Adaptive Control for a Class of Systems with Deadzone Nonlinearity
Intelligent Control an Automation, 5, 6, -9 Publishe Online February 5 in SciRes. http://www.scirp.org/journal/ica http://x.oi.org/.436/ica.5.6 Robust Aaptive Control for a Class of Systems with Deazone
More informationThe Principle of Least Action
Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of
More informationTable of Common Derivatives By David Abraham
Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec
More informationMath 1B, lecture 8: Integration by parts
Math B, lecture 8: Integration by parts Nathan Pflueger 23 September 2 Introuction Integration by parts, similarly to integration by substitution, reverses a well-known technique of ifferentiation an explores
More informationOn classical orthogonal polynomials and differential operators
INSTITUTE OF PHYSICS PUBLISHING JOURNAL OF PHYSICS A: MATHEMATICAL AND GENERAL J. Phys. A: Math. Gen. 38 2005) 6379 6383 oi:10.1088/0305-4470/38/28/010 On classical orthogonal polynomials an ifferential
More informationImplicit Lyapunov control of closed quantum systems
Joint 48th IEEE Conference on Decision an Control an 28th Chinese Control Conference Shanghai, P.R. China, December 16-18, 29 ThAIn1.4 Implicit Lyapunov control of close quantum systems Shouwei Zhao, Hai
More informationDynamics of the synchronous machine
ELEC0047 - Power system ynamics, control an stability Dynamics of the synchronous machine Thierry Van Cutsem t.vancutsem@ulg.ac.be www.montefiore.ulg.ac.be/~vct These slies follow those presente in course
More informationPermanent vs. Determinant
Permanent vs. Determinant Frank Ban Introuction A major problem in theoretical computer science is the Permanent vs. Determinant problem. It asks: given an n by n matrix of ineterminates A = (a i,j ) an
More informationUnit #6 - Families of Functions, Taylor Polynomials, l Hopital s Rule
Unit # - Families of Functions, Taylor Polynomials, l Hopital s Rule Some problems an solutions selecte or aapte from Hughes-Hallett Calculus. Critical Points. Consier the function f) = 54 +. b) a) Fin
More informationComputing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions
Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5
More information23 Implicit differentiation
23 Implicit ifferentiation 23.1 Statement The equation y = x 2 + 3x + 1 expresses a relationship between the quantities x an y. If a value of x is given, then a corresponing value of y is etermine. For
More informationRobust Tracking Control of Robot Manipulator Using Dissipativity Theory
Moern Applie Science July 008 Robust racking Control of Robot Manipulator Using Dissipativity heory Hongrui Wang Key Lab of Inustrial Computer Control Engineering of Hebei Province Yanshan University Qinhuangao
More informationBEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS. Mauro Boccadoro Magnus Egerstedt Paolo Valigi Yorai Wardi
BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS Mauro Boccaoro Magnus Egerstet Paolo Valigi Yorai Wari {boccaoro,valigi}@iei.unipg.it Dipartimento i Ingegneria Elettronica
More informationNOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,
NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which
More informationIEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. XX, NO. XX, MONTH YEAR 1
IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. XX, NO. XX, MONTH YEAR 1 Non-linear complementary filters on the special orthogonal group Robert Mahony, Member, IEEE, Tarek Hamel, Member, IEEE, an Jean-Michel
More informationMath 115 Section 018 Course Note
Course Note 1 General Functions Definition 1.1. A function is a rule that takes certain numbers as inputs an assigns to each a efinite output number. The set of all input numbers is calle the omain of
More informationAn Optimal Algorithm for Bandit and Zero-Order Convex Optimization with Two-Point Feedback
Journal of Machine Learning Research 8 07) - Submitte /6; Publishe 5/7 An Optimal Algorithm for Banit an Zero-Orer Convex Optimization with wo-point Feeback Oha Shamir Department of Computer Science an
More informationAssignment 1. g i (x 1,..., x n ) dx i = 0. i=1
Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition
More informationThe canonical controllers and regular interconnection
Systems & Control Letters ( www.elsevier.com/locate/sysconle The canonical controllers an regular interconnection A.A. Julius a,, J.C. Willems b, M.N. Belur c, H.L. Trentelman a Department of Applie Mathematics,
More information16.30/31, Fall 2010 Recitation # 1
6./, Fall Recitation # September, In this recitation we consiere the following problem. Given a plant with open-loop transfer function.569s +.5 G p (s) = s +.7s +.97, esign a feeback control system such
More informationSchrödinger s equation.
Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of
More informationCriteria for Global Stability of Coupled Systems with Application to Robust Output Feedback Design for Active Surge Control
Criteria for Global Stability of Couple Systems with Application to Robust Output Feeback Design for Active Surge Control Shiriaev, Anton; Johansson, Rolf; Robertsson, Aners; Freiovich, Leoni 9 Link to
More informationinitial configuration initial configuration end
Design of trajectory stabilizing feeback for riftless at systems M. FLIESS y J. L EVINE z P. MARTIN x P. ROUCHON { ECC95 Abstract A esign metho for robust stabilization of riftless at systems aroun trajectories
More informationLecture 3 Basic Feedback
Lecture 3 Basic Feeback Simple control esign an analsis Linear moel with feeback control Use simple moel esign control valiate Simple P loop with an integrator Velocit estimation Time scale Cascae control
More informationOptimal Control of Spatially Distributed Systems
Optimal Control of Spatially Distribute Systems Naer Motee an Ali Jababaie Abstract In this paper, we stuy the structural properties of optimal control of spatially istribute systems. Such systems consist
More informationIntroduction to the Vlasov-Poisson system
Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its
More informationA Sketch of Menshikov s Theorem
A Sketch of Menshikov s Theorem Thomas Bao March 14, 2010 Abstract Let Λ be an infinite, locally finite oriente multi-graph with C Λ finite an strongly connecte, an let p
More informationSystems & Control Letters
Systems & ontrol Letters ( ) ontents lists available at ScienceDirect Systems & ontrol Letters journal homepage: www.elsevier.com/locate/sysconle A converse to the eterministic separation principle Jochen
More informationChaos, Solitons and Fractals Nonlinear Science, and Nonequilibrium and Complex Phenomena
Chaos, Solitons an Fractals (7 64 73 Contents lists available at ScienceDirect Chaos, Solitons an Fractals onlinear Science, an onequilibrium an Complex Phenomena journal homepage: www.elsevier.com/locate/chaos
More informationState-Space Model for a Multi-Machine System
State-Space Moel for a Multi-Machine System These notes parallel section.4 in the text. We are ealing with classically moele machines (IEEE Type.), constant impeance loas, an a network reuce to its internal
More informationAverage value of position for the anharmonic oscillator: Classical versus quantum results
verage value of position for the anharmonic oscillator: Classical versus quantum results R. W. Robinett Department of Physics, The Pennsylvania State University, University Park, Pennsylvania 682 Receive
More informationInterconnected Systems of Fliess Operators
Interconnecte Systems of Fliess Operators W. Steven Gray Yaqin Li Department of Electrical an Computer Engineering Ol Dominion University Norfolk, Virginia 23529 USA Abstract Given two analytic nonlinear
More informationDiagonalization of Matrices Dr. E. Jacobs
Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is
More informationDetermine Power Transfer Limits of An SMIB System through Linear System Analysis with Nonlinear Simulation Validation
Determine Power Transfer Limits of An SMIB System through Linear System Analysis with Nonlinear Simulation Valiation Yin Li, Stuent Member, IEEE, Lingling Fan, Senior Member, IEEE Abstract This paper extens
More informationIndirect Adaptive Fuzzy and Impulsive Control of Nonlinear Systems
International Journal of Automation an Computing 7(4), November 200, 484-49 DOI: 0007/s633-00-053-7 Inirect Aaptive Fuzzy an Impulsive Control of Nonlinear Systems Hai-Bo Jiang School of Mathematics, Yancheng
More informationFree rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012
Free rotation of a rigi boy 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012 1 Introuction In this section, we escribe the motion of a rigi boy that is free to rotate
More informationSwitching Time Optimization in Discretized Hybrid Dynamical Systems
Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set
More informationTime-Optimal Motion Control of Piezoelectric Actuator: STM Application
Time-Optimal Motion Control of Piezoelectric Actuator: STM Application Yongai Xu, Peter H. Mecl Abstract This paper exaes the problem of time-optimal motion control in the context of Scanning Tunneling
More informationSeparation Principle for a Class of Nonlinear Feedback Systems Augmented with Observers
Proceeings of the 17th Worl Congress The International Feeration of Automatic Control Separation Principle for a Class of Nonlinear Feeback Systems Augmente with Observers A. Shiriaev, R. Johansson A.
More information3.7 Implicit Differentiation -- A Brief Introduction -- Student Notes
Fin these erivatives of these functions: y.7 Implicit Differentiation -- A Brief Introuction -- Stuent Notes tan y sin tan = sin y e = e = Write the inverses of these functions: y tan y sin How woul we
More informationVectors in two dimensions
Vectors in two imensions Until now, we have been working in one imension only The main reason for this is to become familiar with the main physical ieas like Newton s secon law, without the aitional complication
More informationA global Implicit Function Theorem without initial point and its applications to control of non-affine systems of high dimensions
J. Math. Anal. Appl. 313 (2006) 251 261 www.elsevier.com/locate/jmaa A global Implicit Function Theorem without initial point an its applications to control of non-affine systems of high imensions Weinian
More informationOnline Appendix for Trade Policy under Monopolistic Competition with Firm Selection
Online Appenix for Trae Policy uner Monopolistic Competition with Firm Selection Kyle Bagwell Stanfor University an NBER Seung Hoon Lee Georgia Institute of Technology September 6, 2018 In this Online
More information