BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS. Mauro Boccadoro Magnus Egerstedt Paolo Valigi Yorai Wardi
|
|
- Ethel Hopkins
- 5 years ago
- Views:
Transcription
1 BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS Mauro Boccaoro Magnus Egerstet Paolo Valigi Yorai Wari Dipartimento i Ingegneria Elettronica e ell Informazione Università i Perugia 06125, Perugia Italy {magnus,ywari}@ece.gatech.eu Electrical an Computer Engineering Georgia Institute of Technology Atlanta, GA Abstract: In this paper we report of a technique to esign optimal feeback control laws for hybri systems with autonomous (continuous) moes. Existing techniques esign the optimal switching surfaces base on a singular sample evolution of the system; hence proviing a solution epenent on the initial conitions. On the other han, the optimal switching times can be foun, proviing an an open loop control to the system, but those also are epenent on the initial conitions. The technique presente relies on a variational approach, giving the erivative of the switching times with respect to the initial conitions, thus proviing a tool to esign programs/algorithms generating switching surfaces which are optimal for any possible execution of the system. Keywors: Hybri Systems, Switching Surfaces, Optimal Control, Variational Methos 1. INTRODUCTION Consier a switche system with autonomous continuous ynamics, ẋ(t) f q(t) (x(t)), (1) q + (t) s(x(t), q(t)). (2) where (1) escribes the continuous ynamics of the state variable x X R n an (2) escribes the iscrete event ynamics of the system. Given an initial conition x 0 : x(t 0 ), the switching law (2) etermines the switching instants t i, i 1, 2,..., an thus the intervals where a certain moal function is active, as well as the initial conition for the o..e. which efines the evolution uner the next moe. The iscrete variable q is piecewise constant in time an belongs to a finite or countable set Q, hence, it can be expresse in terms of the inex i as q(i). In terms of such inex the ynamics of a switche system is: ẋ(t) f i (x(t)), t (t i 1, t i ] (3) i + s(x(t), i, t). (4) with the unerstaning that f i : f q(i), for a given map q(i), i.e., in this case (4) only expresses the occurrence of the i th switch, the specification of the next active moe being given by the map q(i). Since the continuous moes are autonomous, the evolution of the system is etermine by the
2 active moes, accoring to (4). When the function s oes not epen by the (continuous) state variable x, the switching instants are etermine as exogenous inputs, an the system is controlle in open loop (timing control); when s is epenent only on the state variables, the switching law is given in a feeback form, an it may be efine by switching surfaces in the state space. To formulate the problem we are intereste with, consier a simple execution of (3,4) with only one switch, starting at x(t 0 ) x 0 with moe 1, switching to moe 2 at time, an exogenous switch, an terminating either at a fixe final time t 2 or in corresponence of a terminal manifol efine by a function g(x), so that t 2 satisfies g(x(t 2 )) 0. For ease of reference, enote such two sets of possible executions by χ t an χ g, respectively. To fix notation, let the explicit representation of the evolution etermine by moe i be given by x(t) ϕ i (t, s, x(s)), hence, { ϕ1 (t, t x(t) 0, x 0 ) t [t 0, ] (5) ϕ 2 (t,, x( )) t (, t 2 ] Also, let x i : x(t i ), an R : f 1 (x 1 ) f 2 (x 1 ). In this paper the following conventions will be use: 1) vectors are column vectors; 2) the erivative of a scalar, e.g. L, w.r.t. a vector x is a row vector: L x : L x [ L x 1,..., L x n ]. (6) (hence L T x is a column vector). The Hessian matrix is enote by L xx. If f is a (column) vector, function of the vector x i.e., then f [f (1) (x),..., f (n) (x)] T f x : f x f (1) f (1)... x 1 x n f (n) (n) f... x 1 x n Accoring to this convention, for the scalars c, t an the vectors x, y, z, the usual chain rule applies to c(x(t)) an c(x(y)), i.e. t c xẋ, c y c x x y ( v stays for v t ); also: x T y z 1.1 Problem formulation y T x z + x T y z, (7) When the optimal control problem to minimize a cost function J t 0 L(x(t))t (8) is formulate, for some continuously ifferentiable function L, an such that L xx is symmetric, then it is known that when t 1, a (locally) optimal switching time, it satisfies the following conition, see e.g. (Egerstet et al., 2003): c(t 1 ) : pt (t 1 )R(x 1 ) 0 (9) where p T (t), for t [t 1, t 2 ] is given by: p T (t) t L x (x(s))φ 2 (s, t)s + p T (t 2 )Φ 2 (t 2, t) (10) with Φ i the transition matrix of the linearize time-varying system ż(t) fi(x(t)) x z(t), an p T (t 2 ) 0 for fixe final time an p T (t 2 ) L(x 2 )g x (x 2 )/L 2, for an evolution ening at a terminal manifol, where L 2 : g x (x 2 )f(x 2 ), the Lie erivative of g along f 2 evaluate at x 2. Assuming to start from a perturbe initial conition x 0 x 0 + δx 0 ; it is possible to use the information of optimality of t 1, as a switching time, to etermine t 1; in other wors: what is the epenence of the optimal switching time on the initial conitions? This problem is motivate by the etermination of optimal switching surfaces, which ten to solve optimal control problems for autonomous system via the synthesis of feeback laws, which may be pursue for specifications of stability or optimal control. Relevant application of such technique may arise in many areas such as behavior base robotics (Arkin, 1998), or manufacturing systems (Khmelnitsky an Caramanis, 1998) to cite a few. Computational methos exist an are base on the optimization of parametrize switching surfaces (Boccaoro et al., 2005). However, the choice of the optimal values for such parameters epen on the particular trajectory chosen to run an optimization program, an thus, funamentally, on the initial conitions (remin that the we are consiering a system with no continuous inputs). An interesting reference for this type of approach is (Giua et al., 2001), which aresse a timing optimization problem, an iscovere the special structure of the solution for linear quaratic problems. Inee, in that case it is possible to ientify homogeneous regions in the continuous state space, whose bounaries, when reache, etermine the optimal switches, thus proviing a feeback solution to a problem which is formulate in terms of an open loop strategy. Here we explicitly investigate the relation existing between optimal switching times an initial conitions, stuying how the conition of optimality (9) that switching times must satisfy, vary in epenence of the initial conitions.
3 This paper reports the work in progress towar this goal, which is still being pursue. 2. OPTIMAL SWITCHING TIMES V/S INITIAL CONDITIONS It is well known that, uner mil assumptions, executions of switche systems are continuous w.r.t. the initial conitions (Broucke an Arapostathis, 2002). If we assume that also the epenence of c on t 1 as well as t 1 on x 0 is such, we may characterize function t 1 by eriving (9) w.r.t. x 0 an setting this erivative to zero. In fact, if starting from x 0 x 0 + δx 0, it results t 1 t 1 + δt 1 ; then, by continuity, 0 c( t 1 ) c(t 1 ) + δx 0 + o(δx 0 ). Hence, set 0, to satisfy optimality conition for t 1. As we will see this yiels a formula for the variational epenence of t 1 on x 0. To go further, the superscript will be roppe (hence assuming that, x 1 etc. are relative to optimal executions) in orer to reuce the notational buren. By (7) we have that R T p() + p T ( ) R (11) To calculate p(t1), account for the following result, which is reaily verifie: x a t(x) f(s, x)s a t f x (s, x)s f(t, t)t x (12) Then, consiering first the simpler case of fixe final time, by (10, 7, 12) p( ) Φ T 2 (s, )L xx (x(s)) x(s) + Φ T 2 (s, ) L T x (x s) s L T x (x 1)Φ 2 (, ) (13) To get x(s) notice that x( ) ϕ 1 (, t 0, x 0 ), hence x(s) ϕ 2 (s,, ϕ 1 (, t 0, x 0 )) for s [, t 2 ], thus, x(s) x(s) + x(s) x 1 + x(s) x 1 x 1 x 1 x 0 (14) Now, x(s)/ f 2 (x(s)) 1, x(s)/ x 1 Φ 2 (s, ), x 1 / x 0 Φ 1 (, t 0 ), x 1 / f 1 (x 1 ), Φ 2 (, ) I, Φ 2 (s, ) Φ 2 (s, ) f 2(x 1 ) x (15) 1 For time invariant ynamics, [ϕ(s, t + h, x) ϕ(s, t, x)]/h [ϕ(s h, t, x) ϕ(s, t, x)]/h f(x(s)) + o(h). (to be transpose). It results: where I 1 I 3 p( ) I 1 I 2 + I 3 I 4 K (16) I 2 Φ T 2 (s, )L xx (x(s))φ 2 (s, )f 1 (x 1 ) I 4 Φ T 2 (s, )L xx (x(s))f 2 (x(s)) Φ T 2 (s, )L xx (x(s))φ 2 (s, )Φ 1 (, t 0 )s f T 2x (x 1)Φ T 2 (s, )L T x (x(s)) s K L T x (x 1 ) (17) To hanle these, integrate by parts I 2 (letting ), taking into account that L xx (x(s))f 2 (x(s))s L T x (x(s)) we have t2 Φ T 2 (s, )L xx (x(s))f 2 (x(s))s I 4 + Φ T 2 (s, )L T x (x(s)) t2 I 4 + Φ T 2 (t 2, )L T x (x 2) K (18) This leas to the cancellation of I 4 an K in (16). To complete, let s compute R(x 1 )/. Again, notice that x 1 x( ) x[ (x 0 ), x 0 ], hence, R(x 1 ) R [ x (x x1 1) + x ] 1 x [ 0 R x (x 1) f 1 (x 1 ) t ] 1 + Φ 1 (, t 0 ) (19) Multiplying this by p T ( ), (16) by R T from the left an summing up we finally obtain: ( ) [ R T (Qf 1 Φ T 2 (t 2, )L T x (x 2 )) + p T ( )R x f 1 ] + [ R T Q p T ( )R x ] Φ1 (, t 0 ) (20) where f 1 : f 1 (x 1 ), an Q : Φ T 2 (s, )L xx (x(s))φ 2 (s, )s (21) which is a kin of quaratic form co-costate. Notice that the term multiplying t1 above, is a scalar. So, if we know that t 1 is a local optimum for an evolution starting from x 0, then, assuming to start from x 0 x 0 +δx 0, we simply must switch at t 1 + δt 1 + o(δx 0 ). Accoring to (20), δt 1 [R T Q p T ( )R x ]Φ 1 (, t 0 ) δx 0 R T (Qf 1 Φ T 2 (t 2, )L T x (x 2)) + p T ( )R x f 1 (22)
4 3. TOWARD THE CONSTRUCTION OF THE OPTIMAL SWITCHING SURFACES To put in use Eq. (22) assume that one optimal switching time has been erive for a certain sample evolution of the system, e.g. one starting in ˆx 0. Then the optimal switching surfaces are efine by the optimal switching states yiele by the variation on the optimal switching times when initial conitions ifferent than ˆx 0 are consiere. However, it must be pai attention to the fact that the formula erive above works for a fixe final time: inee for the case of evolution ening at a terminal manifol the following result hols, Theorem 1. Consier a nominal an a perturbe execution of the set χ g, x( ) an y( ), respectively, the first starting at x 0 an the latter starting from a point y 0 which lies on the nominal trajectory; i.e., assume that it exists an interval δt 0 such that y 0 ϕ 1 (t 0 + δt 0, t 0, x 0 ). Then, the optimal switching time for all δt 0 < t 1 t 0 t 1(y 0 ) t 1(x 0 ) δt 0 (23) Proof Denote by a b a trajectory from point a to b, an let x(t 1 ) x 1 an x(t 2 ) x 2 where t 2 is the terminal time if the switching time from moe 1 to moe 2 is t 1. If (23) i not hol then assume t 1 (y 0) t 1 (x 0) δt 0 + ɛ (24) for some ɛ (assume with no loss of generality ɛ > 0). Let the nominal trajectory that switches at t 1 (x 0)+ɛ terminate at x ɛ 2. Denote A x(t 0) x(t 0 + δt 0 ), B x(t 0 + δt 0 ) x 1, C x 1 x 2, D x(t 1 ) x(t 1 + ɛ) an E x(t 1 + ɛ) xɛ 2. Accoring to (23) the optimal nominal trajectory is A B C paying for this the cost J(A) +J(B) + J(C). On the other han by (24) the perturbe trajectory B C incurs in a greater cost than B D E, which implies that J(C) > J(D)+J(E). This in turn means that if the nominal trajectory switches at t 1(x 0 ) + ɛ then it pays less than if the switch take place at t 1 Notice that Theorem 1 easily extens to negative δt 0, i.e., if y 0 is chosen such that the evolution starting from y 0 will reach x 0 we must a the time neee to reach x 0 from y 0 to the optimal (nominal) switching time. In case of fixe terminal time the optimal switching state may vary because the perturbe trajectory escribe in Theorem 1 above, switching at t 1 δt 0, reaches the point x(t 2 ) (of the nominal trajectory) at time instant t 2 δt 0, thence visits aitional states from t 2 δt 0 to t 2 (in other wors x( ) (t2 δt 0,t 2] is a set of states not visite by x( )). Such remnants of the perturbe trajectory a further costs, so that two ifferent trajectories, even if the starting point of one of them lies in the trajectory of the other, cannot really be properly compare, in terms of optimal switching states. This can be actually seen: take an i.c. y 0 ϕ 1 (t 0 + δt 0, t 0, x 0 ) very close to x 0, so that δx 0 f 1 (x 0 )δt 0 + o(δt 0 ). Multiplying (22) by such δx 0, we have that its numerator (plus higher orer terms) is: [R T Q + p T ( )R x ]Φ 1 (, t 0 )δx 0 [R T Qf 1 p T ( )R x f 1 ]δt 0 (25) where Φ 1 (, t 0 )f 1 (x 0 ) f 1 (x1) is ue to the fact that vector fiels obey their variational ynamics 2. Hence in this case δt 1 [R T Qf 1 p T ( )R x f 1 ] δt 0 R T (Qf 1 Φ T 2 (t 2, )L T x (x 2)) + p T ( )R x f 1 (26) In this case conition (23) is equivalent to δt 1 δt 0, so that to be verifie, enominator an numerator shoul have ha the same terms, oppose in sign. Here, the only term making the ifference, preventing (23) to hol (as expecte) is R T Φ T 2 (t 2, )L T x (x 2 ). Accoring to Theorem 1, an to the above verification, the proceure to buil optimal switching surfaces shoul be better pursue consiering evolution ening at terminal manifols, since variations in the switching times efine sounly optimal switching states as well. Theorem 1 also gives an hint about the set of initial conitions that shoul be consiere to set such proceure. Inee, it seems reasonable account only for that set of initial conitions which are transversal to the flow efine by the vector fiel of the initial ynamics (here f 1 ) which contains ˆx 0. Such set of initial conition is a surface itself an can be escribe by s(x) 0 where s is a R-value function such that s(ˆx 0 ) 0 an such that s x (x) is collinear with f 1 (x), so that s woul be a kin of potential of the vector fiel. This choice is justifie by Theorem 1, since the components of the variation δx 0 on some x 0 which are tangent to the flow yiel no ifference on the optimal switching state, hence giving no contribution to the construction of an optimal switching surface which is optimal for the executions etermine by any possible initial conition. 4. CONCLUSION AND FUTURE WORKS This paper presents the first steps to etermine optimal switching surfaces for hybri systems 2 Inee, the variational system ż(t) f(x(t)) z(t) has x the solution z(t) f(x(t)), which can be seen from the chain rule f f xf
5 with autonomous moes. Future work will be evote to the erivation of an analogue formula of (22) for executions ening at a terminal manifol. This in orer to pursue the program outline above, about the investigation of the impact of transverse variations in the initial conition on the switching states. REFERENCES Arkin, R.C. (1998). Behavior Base Robotics. The MIT Press. Cambrige, MA. Boccaoro, M., Y. Wari, M. Egerstet an E. Verriest (2005). Optimal control of switching surfaces in hybri ynamical systems. JD- EDS 15(4), Broucke, M. an A. Arapostathis (2002). Continuous selections of trajectories of hybri systems. Systems an Control Letters 47, Egerstet, M., Y. Wari an F. Delmotte (2003). Optimal control of switching times in switche ynamical systems. In: 42n IEEE Conference on Decision an Control (CDC 03). Maui, Hawaii, USA. Giua, A., C. Seatzu an C. Van Der Mee (2001). Optimal control of switche autonomous linear systems. In: 40 th IEEE Conf. on Decision an Control (CDC 2001). Orlano, FL, USA. pp Khmelnitsky, E. an M. Caramanis (1998). Onemachine n-part-type optimal setup scheuling: analytical characterization of switching surfaces. IEEE Trans. on Automatic Control 43(11),
OPTIMAL CONTROL OF SWITCHING SURFACES IN HYBRID DYNAMIC SYSTEMS. Mauro Boccadoro Magnus Egerstedt,1 Yorai Wardi,1
OPTIMAL CONTROL OF SWITCHING SURFACES IN HYBRID DYNAMIC SYSTEMS Mauro Boccadoro Magnus Egerstedt,1 Yorai Wardi,1 boccadoro@diei.unipg.it Dipartimento di Ingegneria Elettronica e dell Informazione Università
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 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 informationMinimum-time constrained velocity planning
7th Meiterranean Conference on Control & Automation Makeonia Palace, Thessaloniki, Greece June 4-6, 9 Minimum-time constraine velocity planning Gabriele Lini, Luca Consolini, Aurelio Piazzi Università
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 informationOptimal Control of Switching Surfaces
Optimal Control of Switching Surfaces Y. Wardi, M. Egerstedt, M. Boccadoro, and E. Verriest {ywardi,magnus,verriest}@ece.gatech.edu School of Electrical and Computer Engineering Georgia Institute of Technology
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 informationLie symmetry and Mei conservation law of continuum system
Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More informationThe Exact Form and General Integrating Factors
7 The Exact Form an General Integrating Factors In the previous chapters, we ve seen how separable an linear ifferential equations can be solve using methos for converting them to forms that can be easily
More informationEuler equations for multiple integrals
Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................
More informationA Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential
Avances in Applie Mathematics an Mechanics Av. Appl. Math. Mech. Vol. 1 No. 4 pp. 573-580 DOI: 10.4208/aamm.09-m0946 August 2009 A Note on Exact Solutions to Linear Differential Equations by the Matrix
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 informationDefine each term or concept.
Chapter Differentiation Course Number Section.1 The Derivative an the Tangent Line Problem Objective: In this lesson you learne how to fin the erivative of a function using the limit efinition an unerstan
More informationIPA Derivatives for Make-to-Stock Production-Inventory Systems With Backorders Under the (R,r) Policy
IPA Derivatives for Make-to-Stock Prouction-Inventory Systems With Backorers Uner the (Rr) Policy Yihong Fan a Benamin Melame b Yao Zhao c Yorai Wari Abstract This paper aresses Infinitesimal Perturbation
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 informationFrom Local to Global Control
Proceeings of the 47th IEEE Conference on Decision an Control Cancun, Mexico, Dec. 9-, 8 ThB. From Local to Global Control Stephen P. Banks, M. Tomás-Roríguez. Automatic Control Engineering Department,
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 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 information6 General properties of an autonomous system of two first order ODE
6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x
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 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 informationMarkov Chains in Continuous Time
Chapter 23 Markov Chains in Continuous Time Previously we looke at Markov chains, where the transitions betweenstatesoccurreatspecifietime- steps. That it, we mae time (a continuous variable) avance in
More informationθ x = f ( x,t) could be written as
9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)
More informationCHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold
CHAPTER 1 : DIFFERENTIABLE MANIFOLDS 1.1 The efinition of a ifferentiable manifol Let M be a topological space. This means that we have a family Ω of open sets efine on M. These satisfy (1), M Ω (2) the
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 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 informationThe total derivative. Chapter Lagrangian and Eulerian approaches
Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function
More informationMath Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors
Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+
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 informationOptimal 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 informationLinear and quadratic approximation
Linear an quaratic approximation November 11, 2013 Definition: Suppose f is a function that is ifferentiable on an interval I containing the point a. The linear approximation to f at a is the linear function
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 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 informationAdjoint Transient Sensitivity Analysis in Circuit Simulation
Ajoint Transient Sensitivity Analysis in Circuit Simulation Z. Ilievski 1, H. Xu 1, A. Verhoeven 1, E.J.W. ter Maten 1,2, W.H.A. Schilers 1,2 an R.M.M. Mattheij 1 1 Technische Universiteit Einhoven; e-mail:
More informationDiscrete Hamilton Jacobi Theory and Discrete Optimal Control
49th IEEE Conference on Decision an Control December 15-17, 2010 Hilton Atlanta Hotel, Atlanta, GA, USA Discrete Hamilton Jacobi Theory an Discrete Optimal Control Tomoi Ohsawa, Anthony M. Bloch, an Melvin
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 informationApproximate reduction of dynamic systems
Systems & Control Letters 57 2008 538 545 www.elsevier.com/locate/sysconle Approximate reuction of ynamic systems Paulo Tabuaa a,, Aaron D. Ames b, Agung Julius c, George J. Pappas c a Department of Electrical
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 informationA. Incorrect! The letter t does not appear in the expression of the given integral
AP Physics C - Problem Drill 1: The Funamental Theorem of Calculus Question No. 1 of 1 Instruction: (1) Rea the problem statement an answer choices carefully () Work the problems on paper as neee (3) Question
More information1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7.
Lectures Nine an Ten The WKB Approximation The WKB metho is a powerful tool to obtain solutions for many physical problems It is generally applicable to problems of wave propagation in which the frequency
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 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 informationThe effect of dissipation on solutions of the complex KdV equation
Mathematics an Computers in Simulation 69 (25) 589 599 The effect of issipation on solutions of the complex KV equation Jiahong Wu a,, Juan-Ming Yuan a,b a Department of Mathematics, Oklahoma State University,
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 informationRelation between the propagator matrix of geodesic deviation and the second-order derivatives of the characteristic function
Journal of Electromagnetic Waves an Applications 203 Vol. 27 No. 3 589 60 http://x.oi.org/0.080/0920507.203.808595 Relation between the propagator matrix of geoesic eviation an the secon-orer erivatives
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 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 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 informationQuantum Mechanics in Three Dimensions
Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.
More informationHow the potentials in different gauges yield the same retarded electric and magnetic fields
How the potentials in ifferent gauges yiel the same retare electric an magnetic fiels José A. Heras a Departamento e Física, E. S. F. M., Instituto Politécnico Nacional, México D. F. México an Department
More information2Algebraic ONLINE PAGE PROOFS. foundations
Algebraic founations. Kick off with CAS. Algebraic skills.3 Pascal s triangle an binomial expansions.4 The binomial theorem.5 Sets of real numbers.6 Surs.7 Review . Kick off with CAS Playing lotto Using
More informationSolutions to Practice Problems Tuesday, October 28, 2008
Solutions to Practice Problems Tuesay, October 28, 2008 1. The graph of the function f is shown below. Figure 1: The graph of f(x) What is x 1 + f(x)? What is x 1 f(x)? An oes x 1 f(x) exist? If so, what
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 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 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 informationImproved Rate-Based Pull and Push Strategies in Large Distributed Networks
Improve Rate-Base Pull an Push Strategies in Large Distribute Networks Wouter Minnebo an Benny Van Hout Department of Mathematics an Computer Science University of Antwerp - imins Mielheimlaan, B-00 Antwerp,
More informationReal-time economic optimization for a fermentation process using Model Predictive Control
Downloae from orbit.tu.k on: Nov 5, 218 Real-time economic optimization for a fermentation process using Moel Preictive Control Petersen, Lars Norbert; Jørgensen, John Bagterp Publishe in: Proceeings of
More informationThe Press-Schechter mass function
The Press-Schechter mass function To state the obvious: It is important to relate our theories to what we can observe. We have looke at linear perturbation theory, an we have consiere a simple moel for
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 informationIntroduction to variational calculus: Lecture notes 1
October 10, 2006 Introuction to variational calculus: Lecture notes 1 Ewin Langmann Mathematical Physics, KTH Physics, AlbaNova, SE-106 91 Stockholm, Sween Abstract I give an informal summary of variational
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 informationLeft-invariant extended Kalman filter and attitude estimation
Left-invariant extene Kalman filter an attitue estimation Silvere Bonnabel Abstract We consier a left-invariant ynamics on a Lie group. One way to efine riving an observation noises is to make them preserve
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 informationStability region estimation for systems with unmodeled dynamics
Stability region estimation for systems with unmoele ynamics Ufuk Topcu, Anrew Packar, Peter Seiler, an Gary Balas Abstract We propose a metho to compute invariant subsets of the robust region-of-attraction
More informationThermal conductivity of graded composites: Numerical simulations and an effective medium approximation
JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University
More informationarxiv: v1 [physics.class-ph] 20 Dec 2017
arxiv:1712.07328v1 [physics.class-ph] 20 Dec 2017 Demystifying the constancy of the Ermakov-Lewis invariant for a time epenent oscillator T. Pamanabhan IUCAA, Post Bag 4, Ganeshkhin, Pune - 411 007, Inia.
More informationEuler Equations: derivation, basic invariants and formulae
Euler Equations: erivation, basic invariants an formulae Mat 529, Lesson 1. 1 Derivation The incompressible Euler equations are couple with t u + u u + p = 0, (1) u = 0. (2) The unknown variable is the
More informationSYNCHRONOUS SEQUENTIAL CIRCUITS
CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents
More informationIntegration Review. May 11, 2013
Integration Review May 11, 2013 Goals: Review the funamental theorem of calculus. Review u-substitution. Review integration by parts. Do lots of integration eamples. 1 Funamental Theorem of Calculus In
More informationFall 2016: Calculus I Final
Answer the questions in the spaces provie on the question sheets. If you run out of room for an answer, continue on the back of the page. NO calculators or other electronic evices, books or notes are allowe
More information05 The Continuum Limit and the Wave Equation
Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,
More informationON THE OPTIMALITY SYSTEM FOR A 1 D EULER FLOW PROBLEM
ON THE OPTIMALITY SYSTEM FOR A D EULER FLOW PROBLEM Eugene M. Cliff Matthias Heinkenschloss y Ajit R. Shenoy z Interisciplinary Center for Applie Mathematics Virginia Tech Blacksburg, Virginia 46 Abstract
More informationChapter 3 Definitions and Theorems
Chapter 3 Definitions an Theorems (from 3.1) Definition of Tangent Line with slope of m If f is efine on an open interval containing c an the limit Δy lim Δx 0 Δx = lim f (c + Δx) f (c) = m Δx 0 Δx exists,
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 informationTRACKING CONTROL OF MULTIPLE MOBILE ROBOTS: A CASE STUDY OF INTER-ROBOT COLLISION-FREE PROBLEM
265 Asian Journal of Control, Vol. 4, No. 3, pp. 265-273, September 22 TRACKING CONTROL OF MULTIPLE MOBILE ROBOTS: A CASE STUDY OF INTER-ROBOT COLLISION-FREE PROBLEM Jurachart Jongusuk an Tsutomu Mita
More informationEnergy behaviour of the Boris method for charged-particle dynamics
Version of 25 April 218 Energy behaviour of the Boris metho for charge-particle ynamics Ernst Hairer 1, Christian Lubich 2 Abstract The Boris algorithm is a wiely use numerical integrator for the motion
More informationSurvey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013
Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing
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 informationComputing Derivatives
Chapter 2 Computing Derivatives 2.1 Elementary erivative rules Motivating Questions In this section, we strive to unerstan the ieas generate by the following important questions: What are alternate notations
More informationINTRODUCTION TO SYMPLECTIC MECHANICS: LECTURE IV. Maurice de Gosson
INTRODUCTION TO SYMPLECTIC MECHANICS: LECTURE IV Maurice e Gosson ii 4 Hamiltonian Mechanics Physically speaking, Hamiltonian mechanics is a paraphrase (an generalization!) of Newton s secon law, popularly
More informationDiscrete Mathematics
Discrete Mathematics 309 (009) 86 869 Contents lists available at ScienceDirect Discrete Mathematics journal homepage: wwwelseviercom/locate/isc Profile vectors in the lattice of subspaces Dániel Gerbner
More informationarxiv: v1 [math-ph] 5 May 2014
DIFFERENTIAL-ALGEBRAIC SOLUTIONS OF THE HEAT EQUATION VICTOR M. BUCHSTABER, ELENA YU. NETAY arxiv:1405.0926v1 [math-ph] 5 May 2014 Abstract. In this work we introuce the notion of ifferential-algebraic
More informationMake graph of g by adding c to the y-values. on the graph of f by c. multiplying the y-values. even-degree polynomial. graph goes up on both sides
Reference 1: Transformations of Graphs an En Behavior of Polynomial Graphs Transformations of graphs aitive constant constant on the outsie g(x) = + c Make graph of g by aing c to the y-values on the graph
More informationarxiv: v2 [math.dg] 16 Dec 2014
A ONOTONICITY FORULA AND TYPE-II SINGULARITIES FOR THE EAN CURVATURE FLOW arxiv:1312.4775v2 [math.dg] 16 Dec 2014 YONGBING ZHANG Abstract. In this paper, we introuce a monotonicity formula for the mean
More informationApproximate Reduction of Dynamical Systems
Proceeings of the 4th IEEE Conference on Decision & Control Manchester Gran Hyatt Hotel San Diego, CA, USA, December 3-, 6 FrIP.7 Approximate Reuction of Dynamical Systems Paulo Tabuaa, Aaron D. Ames,
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 informationLagrangian and Hamiltonian Mechanics
Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical
More informationNested Saturation with Guaranteed Real Poles 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
More informationTopic 7: Convergence of Random Variables
Topic 7: Convergence of Ranom Variables Course 003, 2016 Page 0 The Inference Problem So far, our starting point has been a given probability space (S, F, P). We now look at how to generate information
More informationConvergence rates of moment-sum-of-squares hierarchies for optimal control problems
Convergence rates of moment-sum-of-squares hierarchies for optimal control problems Milan Kora 1, Diier Henrion 2,3,4, Colin N. Jones 1 Draft of September 8, 2016 Abstract We stuy the convergence rate
More informationTEMPORAL AND TIME-FREQUENCY CORRELATION-BASED BLIND SOURCE SEPARATION METHODS. Yannick DEVILLE
TEMPORAL AND TIME-FREQUENCY CORRELATION-BASED BLIND SOURCE SEPARATION METHODS Yannick DEVILLE Université Paul Sabatier Laboratoire Acoustique, Métrologie, Instrumentation Bât. 3RB2, 8 Route e Narbonne,
More informationCode_Aster. Detection of the singularities and calculation of a map of size of elements
Titre : Détection es singularités et calcul une carte [...] Date : 0/0/0 Page : /6 Responsable : DLMAS Josselin Clé : R4.0.04 Révision : Detection of the singularities an calculation of a map of size of
More informationPredictive Control of a Laboratory Time Delay Process Experiment
Print ISSN:3 6; Online ISSN: 367-5357 DOI:0478/itc-03-0005 Preictive Control of a aboratory ime Delay Process Experiment S Enev Key Wors: Moel preictive control; time elay process; experimental results
More informationPerturbation Analysis and Optimization of Stochastic Flow Networks
IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. XX, NO. Y, MMM 2004 1 Perturbation Analysis an Optimization of Stochastic Flow Networks Gang Sun, Christos G. Cassanras, Yorai Wari, Christos G. Panayiotou,
More informationCentrum voor Wiskunde en Informatica
Centrum voor Wiskune en Informatica Moelling, Analysis an Simulation Moelling, Analysis an Simulation Conservation properties of smoothe particle hyroynamics applie to the shallow water equations J.E.
More informationCode_Aster. Detection of the singularities and computation of a card of size of elements
Titre : Détection es singularités et calcul une carte [...] Date : 0/0/0 Page : /6 Responsable : Josselin DLMAS Clé : R4.0.04 Révision : 9755 Detection of the singularities an computation of a car of size
More informationInterpolated Rigid-Body Motions and Robotics
Interpolate Rigi-Boy Motions an Robotics J.M. Selig Faculty of Business, Computing an Info. Management. Lonon South Bank University, Lonon SE AA, U.K. seligjm@lsbu.ac.uk Yaunquing Wu Dept. Mechanical Engineering.
More informationConservation laws a simple application to the telegraph equation
J Comput Electron 2008 7: 47 51 DOI 10.1007/s10825-008-0250-2 Conservation laws a simple application to the telegraph equation Uwe Norbrock Reinhol Kienzler Publishe online: 1 May 2008 Springer Scienceusiness
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 information