Trajectory Smoothing as a Linear Optimal Control Problem
|
|
- Lee Ellis
- 5 years ago
- Views:
Transcription
1 Trajectory Smoothing as a Linear Optimal Control Problem Biswadip Dey & P. S. Krishnaprasad Allerton Conference Monticello, IL, USA. October 4, 212
2 Background and Motivation 1 Background and Motivation 2 Problem Generative Model Inverse Problem Relationship between Linear and Non-linear Generative Models 3 Optimal Control Based Approach for Trajectory Reconstruction Path Independence Lemma Existence of Optimal Initial Condition Optimal Reconstruction as a Linear Smoother Co-State Based Approach 4 Cross-validation Approach to Inverse Problem 5 Numerical Results
3 Background and Motivation Background and Motivation To explore underlying strategies and motion (pursuit, collective motion etc.) governing control laws, by extracting parameters of motion (namely curvature, speed, lateral acceleration etc.) from sampled observations of trajectories. To extract control inputs from sampled data.
4 Problem 1 Background and Motivation 2 Problem Generative Model Inverse Problem Relationship between Linear and Non-linear Generative Models 3 Optimal Control Based Approach for Trajectory Reconstruction Path Independence Lemma Existence of Optimal Initial Condition Optimal Reconstruction as a Linear Smoother Co-State Based Approach 4 Cross-validation Approach to Inverse Problem 5 Numerical Results
5 Problem Generative Model Generative Models for a Curve inr 3 (Non-linear and Linear) Natural Frenet Frame ṙ = νt T = ν(k 1 M 1 +k 2 M 2 ) Ṁ 1 = νk 1 T Ṁ 2 = νk 2 T Trajectory T r M2 M1 (1) The natural curvatures are the steering inputs and the speed is a time function dictated by propulsive/lift/drag mechanisms. Linear Generative Model ṙ = v v = a ȧ = u (2) Jerk, i.e. the third-derivative of position, is viewed as the control. LTI representation ẋ = Ax+Bu r = Cx with, x = [ r T v T a T ] T ; A = I I ;B = I C = [ I ] ; (3) Controllable and Observable
6 Problem Inverse Problem Regularized Inverse Problem T M2 : Data Points r M1 Trajectory Given a time series of observed positions, generate a smooth trajectory to fit the data points. The inverse problem is ill-posed. Highly sensitive to noise. Non-unique. A regularization parameter is introduced to control the amount of smoothing. Ordinary cross validation is a standard approach to choose an optimal value for the regularization parameter.
7 Problem Inverse Problem Extracting Curvature (Inverse Problem) Non-linear Optimization Minimize Observations r i Regularized Inversion States, Controls r, [T,M 1,M 2],k 1,k 2,ν ( N ) T r(t i ) r i 2 +λ ( k k ν 2 )dt i= subject to Dynamics in (1), Initial Condition, and Input (4) Linear-Quadratic Control Minimize Observations r i ( N i= Regularized Inversion States, Controls r,v,a,u ) T r(t i ) r i 2 +λ u T udt subject to Dynamics in (3), Initial Condition, and Input (5)
8 Problem Relationship between Linear and Non-linear Generative Models Relationship between Two Approaches for Modelling a Curve Natural-Frenet Frame Linear Model (Triple Integrator) v = νt a = νt +ν 2 k 1 M 1 +ν 2 k 2 M 2 u = ( ν ν 3 (k1 2 +k2 2 ))T +(3ν νk 1 +ν 2 k1 )M 1 +(3ν νk 2 +ν 2 k2 )M 2 Linear Model (Triple Integrator) Natural-Frenet Frame ν = v T = v v T = 1 ν (a (a T)T) κ = T ν τ = v (a u) v a 2 k 1,k 2,M 1,M 2 can be computed by assuming suitable intial conditions. ( t ) k 1 (t) = κcos θ + τ(σ)dσ ( k 2 (t) = κsin θ + t ) τ(σ)dσ t M 1 (t) = M 1 () ν(σ)k 1 (σ)t(σ)dσ t M 2 (t) = M 2 () ν(σ)k 2 (σ)t(σ)dσ
9 Optimal Control Based Approach for Trajectory Reconstruction 1 Background and Motivation 2 Problem Generative Model Inverse Problem Relationship between Linear and Non-linear Generative Models 3 Optimal Control Based Approach for Trajectory Reconstruction Path Independence Lemma Existence of Optimal Initial Condition Optimal Reconstruction as a Linear Smoother Co-State Based Approach 4 Cross-validation Approach to Inverse Problem 5 Numerical Results
10 Optimal Control Based Approach for Trajectory Reconstruction Path Independence Lemma Application of Path Independence Lemma Optimal Control Problem: Minimize x(t ),u J(x(t ),u) = subject to x(t ) R n, u U, N i= Dynamics in (3) T r(t i ) r i 2 +λ u T udt (6) Path Independence: Along trajectories of (3) = x T (t i )K(t + i )x(t i) x T (t i+1 )K(t i+1 )x(t i+1) t [ ] T [ i+1 x K +A + T ][ K +KA KB x u B T K u t + i = x T (t i )η(t + i ) xt (t i+1 )η(t i+1 )+ t t + i i+1 [ x u ] dt ] T [ η +A T η B T η ] dt for all i {,1,,N 1}
11 Optimal Control Based Approach for Trajectory Reconstruction Path Independence Lemma Application of Path Independence Lemma Assumptions on the the dynamics and boundary values of K and η: K = A T K KA+KBB T K, K(t + N ) =, K(t + i ) K(t i ) = 1 λ CT C. (7) With the assumptions (7) and (8), we obtain η = ( A T KBB T) η, η(t + N ) =, η(t + i ) η(t i ) = 2 λ CT r i. ] N J(x(t ),u) = λ [x T (t )K(t )x(t )+x T (t )η(t ) + ri T r i 1 T 4 λ B T η(t) 2 dt i= T ( +λ u(t) +B T K(t)x(t)+ 1 ) 2 η(t) 2 dt. (9) Optimal control input: u opt (t) = B T ( K(t)x(t)+ 1 2 η(t) ) (1) Optimal initial condition: ] [K(t ) x opt (t )+ 1 2 η(t ) =. (11) (8)
12 Optimal Control Based Approach for Trajectory Reconstruction Existence of Optimal Initial Condition Existence of Solution for (11) - Sketch of Proof Proposition 1 The solution of the Riccati equation (7) assumes the form K(t i ) = 1 N Φ Σ (t i,t k )C T CΦ T Σ λ (t i,t k ) k=i for any i {,1,,N} where Σ(t) = (A 1 2 BBT K(t)) T and Φ Σ is the transition matrix of Σ. Holds true for i = N. Apply mathematical induction. Proposition 2 ( Σ T,C) forms an observable pair for the problem of our interest (3). Apply Silverman-Meadows rank condition.
13 Optimal Control Based Approach for Trajectory Reconstruction Existence of Optimal Initial Condition Existence of Solution for (11) - Sketch of Proof Theorem 1 The equation [K(t ) ] x opt (t )+ 1 2 η(t ) =. is uniquely solvable for almost any time index set {t i } N i=. Observe K(t ) can be represented as K(t ) = 1 λ CT C, with C CΦ Σ T(t 1,t ) C =.. CΦ Σ T(t N,t ) Consider the system ξ = Σ T ξ; γ = Cξ. The outputs, corresponding to two different initial conditions, do not match identically over any interval. ξ a ξ b Cξ a Cξ b (almost surely) Otherwise, consider an arbitrary close perturbation of the original time index set {t i } N i=, to obtain full rank for C.
14 Optimal Control Based Approach for Trajectory Reconstruction Optimal Reconstruction as a Linear Smoother Linearity in the Reconstructed Trajectory Closed loop dynamics: ẋ(t) = Σ T x(t) 1 2 BBT η(t) with Σ = [ A BB T K(t) ] T. x opt(t ) and η( ) are linear in observed data {r i } N i=. where r(t k ) = 1 N [ CF λ λ (k,i)c T] r i (12) i= F λ (k,i) = Φ T Σ(t [ ] 1Φ Σ(t,t k ) K(t ),t i ) + min{i,k} j=1 ( tj t j 1 Φ T Σ(σ,t k )BB T Φ Σ(σ,t i )dσ ) Can be be viewed as a global alternative to Savitzky-Golay smoothing filters. Can be used as a building block to obtain a fixed lag smoothing algorithm.
15 Optimal Control Based Approach for Trajectory Reconstruction Co-State Based Approach An Alternative Co-State Based Approach Co-state variables: p(t) K(t)x(t)+ 1 2 η(t) An optimal trajectory between two observation times can be viewed as the base integral curve of the following Hamiltonian dynamics [ ] [ ][ ] d x(t) A BB T x(t) = dt p(t) A T p(t) Jump condition for the co-state variables: p(t + i ) p(t i ) = 1 λ CT (r i r(t i )) Terminal condition for the co-state variables: p(t + N ) = p(t ) =
16 Optimal Control Based Approach for Trajectory Reconstruction Co-State Based Approach An Alternative Co-State Based Approach Forward-propagation of x(t i ) and p(t + i ): [ ] [ x(ti+1 ) e A i e A iw i p(t + i+1 ) = ] 1 λ CT Ce A i [e AT i + 1λ CT Ce A iw i [ ] + 1 r i+1 λ CT where W i is defined as ] [ x(ti ) p(t + i ) ] W i = i e Aσ BB T e ATσ dσ ( i = t i+1 t i ) Optimal initial condition is obtained by solving ( N 1 ) [ ] I N N 1 [ I] Λ i 1 x(t i= λ CT C ) = [ I] Λ j Γr i (13) i= j=i where, [ Λ i = e A i 1 λ CT Ce A i e A iw i ] [e AT i + 1 λ CT Ce A iw i ] [ ;Γ = 1 λ CT ]
17 Cross-validation Approach to Inverse Problem 1 Background and Motivation 2 Problem Generative Model Inverse Problem Relationship between Linear and Non-linear Generative Models 3 Optimal Control Based Approach for Trajectory Reconstruction Path Independence Lemma Existence of Optimal Initial Condition Optimal Reconstruction as a Linear Smoother Co-State Based Approach 4 Cross-validation Approach to Inverse Problem 5 Numerical Results
18 Cross-validation Approach to Inverse Problem Cross-validation Approach to Determination of Penalty Parameter We use leaving-out-one version of the Ordinary Cross Validation (OCV) technique. Let, {x [λ,k] opt,u[λ,k] } be a minimizer of: N T r(t i ) r i 2 +λ u T udt i= i k Let the reconstructed trajectory be r [λ,k] ( ). Then the OCV cost is defined as: V (λ) = 1 N +1 N r [λ,k] (t k ) r k 2 k= Hence, OCV estimate for λ is defined as: λ = argmin V (λ) λ R +
19 Numerical Results 1 Background and Motivation 2 Problem Generative Model Inverse Problem Relationship between Linear and Non-linear Generative Models 3 Optimal Control Based Approach for Trajectory Reconstruction Path Independence Lemma Existence of Optimal Initial Condition Optimal Reconstruction as a Linear Smoother Co-State Based Approach 4 Cross-validation Approach to Inverse Problem 5 Numerical Results
20 Numerical Results Numerical Result - Spherical Curve Avg. Fit Error/Radius:
21 Numerical Results Numerical Result - Circular Helix Avg. Fit Error/Radius:
22 References References E. Justh and P. S. Krishnaprasad, Optimal Natural Frames, Comm. Inf. Syst., 11(1):17-34, 211. T. Flash and N. Hogan, The coordination of arm movements: An experimentally confirmed mathematical model, The Journal of Neuroscience, 5(7): , R. L. Bishop, There is more than one way to frame a curve, The American Mathematical Monthly, 82(3): , L. M. Silverman and H. E. Meadows, Controllability and observability in time-variable linear systems, SIAM J. on Control and Optimization, 5(1):64-73, P.V. Reddy, Steering laws for pursuit, M. S. Thesis, University of Maryland, College Park, 27. B. Dey and P. S. Krishnaprasad, Trajectory smoothing as a linear optimal control problem (extended version),
23 ???
Homework Solution # 3
ECSE 644 Optimal Control Feb, 4 Due: Feb 17, 4 (Tuesday) Homework Solution # 3 1 (5%) Consider the discrete nonlinear control system in Homework # For the optimal control and trajectory that you have found
More informationExtremal Collective Behavior
49th IEEE Conference on Decision and Control December 5-7, Hilton Atlanta Hotel, Atlanta, GA, USA Extremal Collective Behavior E W Justh and P S Krishnaprasad Abstract Curves and natural frames can be
More information1 Controllability and Observability
1 Controllability and Observability 1.1 Linear Time-Invariant (LTI) Systems State-space: Dimensions: Notation Transfer function: ẋ = Ax+Bu, x() = x, y = Cx+Du. x R n, u R m, y R p. Note that H(s) is always
More information5th-order differentiation
ARBITRARY-ORDER REAL-TIME EXACT ROBUST DIFFERENTIATION A. Levant Applied Mathematics Dept., Tel-Aviv University, Israel E-mail: levant@post.tau.ac.il Homepage: http://www.tau.ac.il/~levant/ 5th-order differentiation
More informationMotion Camouflage for Coverage
1 American Control Conference Marriott Waterfront, Baltimore, MD, USA June 3-July, 1 FrC5.4 Motion Camouflage for Coverage M. Mischiati and P. S. Krishnaprasad Abstract Pursuit strategies can lead to cohesive
More informationGramians based model reduction for hybrid switched systems
Gramians based model reduction for hybrid switched systems Y. Chahlaoui Younes.Chahlaoui@manchester.ac.uk Centre for Interdisciplinary Computational and Dynamical Analysis (CICADA) School of Mathematics
More informationModel reduction for linear systems by balancing
Model reduction for linear systems by balancing Bart Besselink Jan C. Willems Center for Systems and Control Johann Bernoulli Institute for Mathematics and Computer Science University of Groningen, Groningen,
More informationTopic # Feedback Control Systems
Topic #17 16.31 Feedback Control Systems Deterministic LQR Optimal control and the Riccati equation Weight Selection Fall 2007 16.31 17 1 Linear Quadratic Regulator (LQR) Have seen the solutions to the
More informationDeterministic Dynamic Programming
Deterministic Dynamic Programming 1 Value Function Consider the following optimal control problem in Mayer s form: V (t 0, x 0 ) = inf u U J(t 1, x(t 1 )) (1) subject to ẋ(t) = f(t, x(t), u(t)), x(t 0
More informationThree-Dimensional Motion Coordination in a Spatiotemporal Flowfield
IEEE TRANSACTIONS ON AUTOMATIC CONTROL 1 Three-Dimensional Motion Coordination in a Spatiotemporal Flowfield Sonia Hernandez and Dere A. Paley, Member, IEEE Abstract Decentralized algorithms to stabilize
More informationNonlinear Control Systems
Nonlinear Control Systems António Pedro Aguiar pedro@isr.ist.utl.pt 7. Feedback Linearization IST-DEEC PhD Course http://users.isr.ist.utl.pt/%7epedro/ncs1/ 1 1 Feedback Linearization Given a nonlinear
More informationDuality and dynamics in Hamilton-Jacobi theory for fully convex problems of control
Duality and dynamics in Hamilton-Jacobi theory for fully convex problems of control RTyrrell Rockafellar and Peter R Wolenski Abstract This paper describes some recent results in Hamilton- Jacobi theory
More informationChap. 3. Controlled Systems, Controllability
Chap. 3. Controlled Systems, Controllability 1. Controllability of Linear Systems 1.1. Kalman s Criterion Consider the linear system ẋ = Ax + Bu where x R n : state vector and u R m : input vector. A :
More informationNew ideas in the non-equilibrium statistical physics and the micro approach to transportation flows
New ideas in the non-equilibrium statistical physics and the micro approach to transportation flows Plenary talk on the conference Stochastic and Analytic Methods in Mathematical Physics, Yerevan, Armenia,
More informationRobotics. Dynamics. Marc Toussaint U Stuttgart
Robotics Dynamics 1D point mass, damping & oscillation, PID, dynamics of mechanical systems, Euler-Lagrange equation, Newton-Euler recursion, general robot dynamics, joint space control, reference trajectory
More informationAn introduction to Mathematical Theory of Control
An introduction to Mathematical Theory of Control Vasile Staicu University of Aveiro UNICA, May 2018 Vasile Staicu (University of Aveiro) An introduction to Mathematical Theory of Control UNICA, May 2018
More informationThe Frenet Serret formulas
The Frenet Serret formulas Attila Máté Brooklyn College of the City University of New York January 19, 2017 Contents Contents 1 1 The Frenet Serret frame of a space curve 1 2 The Frenet Serret formulas
More informationA new statistical framework to infer gene regulatory networks with hidden transcription factors
A new statistical framework to infer gene regulatory networks with hidden transcription factors Ivan Vujačić, Javier González, Ernst Wit University of Groningen Johann Bernoulli Institute PEDS II, June
More informationEvent-triggered stabilization of linear systems under channel blackouts
Event-triggered stabilization of linear systems under channel blackouts Pavankumar Tallapragada, Massimo Franceschetti & Jorge Cortés Allerton Conference, 30 Sept. 2015 Acknowledgements: National Science
More informationStationary trajectories, singular Hamiltonian systems and ill-posed Interconnection
Stationary trajectories, singular Hamiltonian systems and ill-posed Interconnection S.C. Jugade, Debasattam Pal, Rachel K. Kalaimani and Madhu N. Belur Department of Electrical Engineering Indian Institute
More informationII. Unit Speed Curves
The Geometry of Curves, Part I Rob Donnelly From Murray State University s Calculus III, Fall 2001 note: This material supplements Sections 13.3 and 13.4 of the text Calculus with Early Transcendentals,
More informationProblem 1 Cost of an Infinite Horizon LQR
THE UNIVERSITY OF TEXAS AT SAN ANTONIO EE 5243 INTRODUCTION TO CYBER-PHYSICAL SYSTEMS H O M E W O R K # 5 Ahmad F. Taha October 12, 215 Homework Instructions: 1. Type your solutions in the LATEX homework
More informationInfeasibility Detection and an Inexact Active-Set Method for Large-Scale Nonlinear Optimization
Infeasibility Detection and an Inexact Active-Set Method for Large-Scale Nonlinear Optimization Frank E. Curtis, Lehigh University involving joint work with James V. Burke, University of Washington Daniel
More informationarxiv: v1 [cs.sy] 29 Sep 2017
Beacon-referenced Mutual Pursuit in Three Dimensions Kevin S. Galloway and Biswadip Dey arxiv:709.037v cs.sy 9 Sep 07 Abstract Motivated by station-keeping applications in various unmanned settings, this
More informationMin-Max Output Integral Sliding Mode Control for Multiplant Linear Uncertain Systems
Proceedings of the 27 American Control Conference Marriott Marquis Hotel at Times Square New York City, USA, July -3, 27 FrC.4 Min-Max Output Integral Sliding Mode Control for Multiplant Linear Uncertain
More informationProject Optimal Control: Optimal Control with State Variable Inequality Constraints (SVIC)
Project Optimal Control: Optimal Control with State Variable Inequality Constraints (SVIC) Mahdi Ghazaei, Meike Stemmann Automatic Control LTH Lund University Problem Formulation Find the maximum of the
More information6.241 Dynamic Systems and Control
6.241 Dynamic Systems and Control Lecture 24: H2 Synthesis Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology May 4, 2011 E. Frazzoli (MIT) Lecture 24: H 2 Synthesis May
More informationPresenter: Noriyoshi Fukaya
Y. Martel, F. Merle, and T.-P. Tsai, Stability and Asymptotic Stability in the Energy Space of the Sum of N Solitons for Subcritical gkdv Equations, Comm. Math. Phys. 31 (00), 347-373. Presenter: Noriyoshi
More informationState Space Representation of Gaussian Processes
State Space Representation of Gaussian Processes Simo Särkkä Department of Biomedical Engineering and Computational Science (BECS) Aalto University, Espoo, Finland June 12th, 2013 Simo Särkkä (Aalto University)
More informationNonlinear Analysis 71 (2009) Contents lists available at ScienceDirect. Nonlinear Analysis. journal homepage:
Nonlinear Analysis 71 2009 2744 2752 Contents lists available at ScienceDirect Nonlinear Analysis journal homepage: www.elsevier.com/locate/na A nonlinear inequality and applications N.S. Hoang A.G. Ramm
More information13 Path Planning Cubic Path P 2 P 1. θ 2
13 Path Planning Path planning includes three tasks: 1 Defining a geometric curve for the end-effector between two points. 2 Defining a rotational motion between two orientations. 3 Defining a time function
More informationObjective. 1 Specification/modeling of the controlled system. 2 Specification of a performance criterion
Optimal Control Problem Formulation Optimal Control Lectures 17-18: Problem Formulation Benoît Chachuat Department of Chemical Engineering Spring 2009 Objective Determine the control
More informationTutorial on Control and State Constrained Optimal Control Pro. Control Problems and Applications Part 3 : Pure State Constraints
Tutorial on Control and State Constrained Optimal Control Problems and Applications Part 3 : Pure State Constraints University of Münster Institute of Computational and Applied Mathematics SADCO Summer
More informationSubject: Optimal Control Assignment-1 (Related to Lecture notes 1-10)
Subject: Optimal Control Assignment- (Related to Lecture notes -). Design a oil mug, shown in fig., to hold as much oil possible. The height and radius of the mug should not be more than 6cm. The mug must
More informationLinear-Quadratic-Gaussian (LQG) Controllers and Kalman Filters
Linear-Quadratic-Gaussian (LQG) Controllers and Kalman Filters Emo Todorov Applied Mathematics and Computer Science & Engineering University of Washington Winter 204 Emo Todorov (UW) AMATH/CSE 579, Winter
More informationTheorema Egregium, Intrinsic Curvature
Theorema Egregium, Intrinsic Curvature Cartan s second structural equation dω 12 = K θ 1 θ 2, (1) is, as we have seen, a powerful tool for computing curvature. But it also has far reaching theoretical
More informationOHSx XM521 Multivariable Differential Calculus: Homework Solutions 14.1
OHSx XM5 Multivariable Differential Calculus: Homework Solutions 4. (8) Describe the graph of the equation. r = i + tj + (t )k. Solution: Let y(t) = t, so that z(t) = t = y. In the yz-plane, this is just
More informationOPTIMAL CONTROL OF A CLASS OF PARABOLIC
OPTIMAL CONTROL OF A CLASS OF PARABOLIC TIME-VARYING PDES James Ng a, Ilyasse Aksikas a,b, Stevan Dubljevic a a Department of Chemical & Materials Engineering, University of Alberta, Canada b Department
More information5.2 The Levi-Civita Connection on Surfaces. 1 Parallel transport of vector fields on a surface M
5.2 The Levi-Civita Connection on Surfaces In this section, we define the parallel transport of vector fields on a surface M, and then we introduce the concept of the Levi-Civita connection, which is also
More informationOPTIMAL CONTROL. Sadegh Bolouki. Lecture slides for ECE 515. University of Illinois, Urbana-Champaign. Fall S. Bolouki (UIUC) 1 / 28
OPTIMAL CONTROL Sadegh Bolouki Lecture slides for ECE 515 University of Illinois, Urbana-Champaign Fall 2016 S. Bolouki (UIUC) 1 / 28 (Example from Optimal Control Theory, Kirk) Objective: To get from
More informationComplex Laplacians and Applications in Multi-Agent Systems
1 Complex Laplacians and Applications in Multi-Agent Systems Jiu-Gang Dong, and Li Qiu, Fellow, IEEE arxiv:1406.186v [math.oc] 14 Apr 015 Abstract Complex-valued Laplacians have been shown to be powerful
More informationInput to state Stability
Input to state Stability Mini course, Universität Stuttgart, November 2004 Lars Grüne, Mathematisches Institut, Universität Bayreuth Part IV: Applications ISS Consider with solutions ϕ(t, x, w) ẋ(t) =
More informationCHAPTER 4 DIFFERENTIAL VECTOR CALCULUS
CHAPTER 4 DIFFERENTIAL VECTOR CALCULUS 4.1 Vector Functions 4.2 Calculus of Vector Functions 4.3 Tangents REVIEW: Vectors Scalar a quantity only with its magnitude Example: temperature, speed, mass, volume
More informationCALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems. CDS 110b
CALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems CDS 110b R. M. Murray Kalman Filters 14 January 2007 Reading: This set of lectures provides a brief introduction to Kalman filtering, following
More information16.31 Fall 2005 Lecture Presentation Mon 31-Oct-05 ver 1.1
16.31 Fall 2005 Lecture Presentation Mon 31-Oct-05 ver 1.1 Charles P. Coleman October 31, 2005 1 / 40 : Controllability Tests Observability Tests LEARNING OUTCOMES: Perform controllability tests Perform
More informationMULTI-AGENT TRACKING OF A HIGH-DIMENSIONAL ACTIVE LEADER WITH SWITCHING TOPOLOGY
Jrl Syst Sci & Complexity (2009) 22: 722 731 MULTI-AGENT TRACKING OF A HIGH-DIMENSIONAL ACTIVE LEADER WITH SWITCHING TOPOLOGY Yiguang HONG Xiaoli WANG Received: 11 May 2009 / Revised: 16 June 2009 c 2009
More informationLinear-quadratic control problem with a linear term on semiinfinite interval: theory and applications
Linear-quadratic control problem with a linear term on semiinfinite interval: theory and applications L. Faybusovich T. Mouktonglang Department of Mathematics, University of Notre Dame, Notre Dame, IN
More informationFeedback Optimal Control of Low-thrust Orbit Transfer in Central Gravity Field
Vol. 4, No. 4, 23 Feedback Optimal Control of Low-thrust Orbit Transfer in Central Gravity Field Ashraf H. Owis Department of Astronomy, Space and Meteorology, Faculty of Science, Cairo University Department
More informationGame Theory Extra Lecture 1 (BoB)
Game Theory 2014 Extra Lecture 1 (BoB) Differential games Tools from optimal control Dynamic programming Hamilton-Jacobi-Bellman-Isaacs equation Zerosum linear quadratic games and H control Baser/Olsder,
More informationProbabilistic Graphical Models
Probabilistic Graphical Models Brown University CSCI 2950-P, Spring 2013 Prof. Erik Sudderth Lecture 12: Gaussian Belief Propagation, State Space Models and Kalman Filters Guest Kalman Filter Lecture by
More informationObservability. It was the property in Lyapunov stability which allowed us to resolve that
Observability We have seen observability twice already It was the property which permitted us to retrieve the initial state from the initial data {u(0),y(0),u(1),y(1),...,u(n 1),y(n 1)} It was the property
More information9 Controller Discretization
9 Controller Discretization In most applications, a control system is implemented in a digital fashion on a computer. This implies that the measurements that are supplied to the control system must be
More informationRegularization of the Inverse Laplace Transform with Applications in Nuclear Magnetic Resonance Relaxometry Candidacy Exam
Regularization of the Inverse Laplace Transform with Applications in Nuclear Magnetic Resonance Relaxometry Candidacy Exam Applied Mathematics, Applied Statistics, & Scientific Computation University of
More informationAlberto Bressan. Department of Mathematics, Penn State University
Non-cooperative Differential Games A Homotopy Approach Alberto Bressan Department of Mathematics, Penn State University 1 Differential Games d dt x(t) = G(x(t), u 1(t), u 2 (t)), x(0) = y, u i (t) U i
More informationAn introduction to Birkhoff normal form
An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an
More informationAn Extension of. Friedmann-Robertson-Walker Theory beyond Big Bang
c 2009 International Press Adv. Theor. Math. Phys. 12 (2009) 1 28 An Extension of arxiv:0906.5504v1 [gr-qc] 30 Jun 2009 Friedmann-Robertson-Walker Theory beyond Big Bang Joachim Schröter University of
More informationTrajectory Tracking Control of Bimodal Piecewise Affine Systems
25 American Control Conference June 8-1, 25. Portland, OR, USA ThB17.4 Trajectory Tracking Control of Bimodal Piecewise Affine Systems Kazunori Sakurama, Toshiharu Sugie and Kazushi Nakano Abstract This
More informationAn Optimal Control Problem for Rigid Body Motions in Minkowski Space
Applied Mathematical Sciences, Vol. 5, 011, no. 5, 559-569 An Optimal Control Problem for Rigid Body Motions in Minkowski Space Nemat Abazari Department of Mathematics, Ardabil Branch Islamic Azad University,
More informationInput-output finite-time stabilization for a class of hybrid systems
Input-output finite-time stabilization for a class of hybrid systems Francesco Amato 1 Gianmaria De Tommasi 2 1 Università degli Studi Magna Græcia di Catanzaro, Catanzaro, Italy, 2 Università degli Studi
More information= m(0) + 4e 2 ( 3e 2 ) 2e 2, 1 (2k + k 2 ) dt. m(0) = u + R 1 B T P x 2 R dt. u + R 1 B T P y 2 R dt +
ECE 553, Spring 8 Posted: May nd, 8 Problem Set #7 Solution Solutions: 1. The optimal controller is still the one given in the solution to the Problem 6 in Homework #5: u (x, t) = p(t)x k(t), t. The minimum
More informationHIGHER ORDER SLIDING MODES AND ARBITRARY-ORDER EXACT ROBUST DIFFERENTIATION
HIGHER ORDER SLIDING MODES AND ARBITRARY-ORDER EXACT ROBUST DIFFERENTIATION A. Levant Institute for Industrial Mathematics, 4/24 Yehuda Ha-Nachtom St., Beer-Sheva 843, Israel Fax: +972-7-232 and E-mail:
More informationLinear Quadratic Zero-Sum Two-Person Differential Games Pierre Bernhard June 15, 2013
Linear Quadratic Zero-Sum Two-Person Differential Games Pierre Bernhard June 15, 2013 Abstract As in optimal control theory, linear quadratic (LQ) differential games (DG) can be solved, even in high dimension,
More informationTangent and Normal Vectors
Tangent and Normal Vectors MATH 311, Calculus III J. Robert Buchanan Department of Mathematics Fall 2011 Navigation When an observer is traveling along with a moving point, for example the passengers in
More informationStochastic Differential Equations.
Chapter 3 Stochastic Differential Equations. 3.1 Existence and Uniqueness. One of the ways of constructing a Diffusion process is to solve the stochastic differential equation dx(t) = σ(t, x(t)) dβ(t)
More informationInterpolated Rigid-Body Motions and Robotics
Interpolated Rigid-Body Motions and Robotics J.M. Selig London South Bank University and Yuanqing Wu Shanghai Jiaotong University. IROS Beijing 2006 p.1/22 Introduction Interpolation of rigid motions important
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science : Dynamic Systems Spring 2011
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science 6.4: Dynamic Systems Spring Homework Solutions Exercise 3. a) We are given the single input LTI system: [
More informatione 2 = e 1 = e 3 = v 1 (v 2 v 3 ) = det(v 1, v 2, v 3 ).
3. Frames In 3D space, a sequence of 3 linearly independent vectors v 1, v 2, v 3 is called a frame, since it gives a coordinate system (a frame of reference). Any vector v can be written as a linear combination
More informationLinear System Theory
Linear System Theory Wonhee Kim Chapter 6: Controllability & Observability Chapter 7: Minimal Realizations May 2, 217 1 / 31 Recap State space equation Linear Algebra Solutions of LTI and LTV system Stability
More informationTAM 212 Worksheet 3: Track design. Solutions
Name: Group members: TAM Worksheet 3: Track design Solutions In this worksheet we will design one curve on the track of the Chicago O Hare Airport Transit System ATS):. On the diagram below, sketch what
More informationarxiv: v1 [math.pr] 18 Oct 2016
An Alternative Approach to Mean Field Game with Major and Minor Players, and Applications to Herders Impacts arxiv:161.544v1 [math.pr 18 Oct 216 Rene Carmona August 7, 218 Abstract Peiqi Wang The goal
More informationTangent and Normal Vector - (11.5)
Tangent and Normal Vector - (.5). Principal Unit Normal Vector Let C be the curve traced out by the vector-valued function rt vector T t r r t t is the unit tangent vector to the curve C. Now define N
More informationThere is a function, the arc length function s(t) defined by s(t) = It follows that r(t) = p ( s(t) )
MATH 20550 Acceleration, Curvature and Related Topics Fall 2016 The goal of these notes is to show how to compute curvature and torsion from a more or less arbitrary parametrization of a curve. We will
More informationSpezielle Funktionen 4 Skalare Kugelfunktionen
Spezielle Funktionen 4 Skalare Kugelfunktionen M. Gutting 28. Mai 2015 4.1 Scalar Spherical Harmonics We start by introducing some basic spherical notation. S 2 = {ξ R 3 ξ = 1} is the unit sphere in R
More informationEXAM 2 ANSWERS AND SOLUTIONS, MATH 233 WEDNESDAY, OCTOBER 18, 2000
EXAM 2 ANSWERS AND SOLUTIONS, MATH 233 WEDNESDAY, OCTOBER 18, 2000 This examination has 30 multiple choice questions. Problems are worth one point apiece, for a total of 30 points for the whole examination.
More informationLecture 1 From Continuous-Time to Discrete-Time
Lecture From Continuous-Time to Discrete-Time Outline. Continuous and Discrete-Time Signals and Systems................. What is a signal?................................2 What is a system?.............................
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More informationOn feedback stabilizability of time-delay systems in Banach spaces
On feedback stabilizability of time-delay systems in Banach spaces S. Hadd and Q.-C. Zhong q.zhong@liv.ac.uk Dept. of Electrical Eng. & Electronics The University of Liverpool United Kingdom Outline Background
More informationAC&ST AUTOMATIC CONTROL AND SYSTEM THEORY SYSTEMS AND MODELS. Claudio Melchiorri
C. Melchiorri (DEI) Automatic Control & System Theory 1 AUTOMATIC CONTROL AND SYSTEM THEORY SYSTEMS AND MODELS Claudio Melchiorri Dipartimento di Ingegneria dell Energia Elettrica e dell Informazione (DEI)
More informationOn the stochastic nonlinear Schrödinger equation
On the stochastic nonlinear Schrödinger equation Annie Millet collaboration with Z. Brzezniak SAMM, Paris 1 and PMA Workshop Women in Applied Mathematics, Heraklion - May 3 211 Outline 1 The NL Shrödinger
More informationAdaptive Piecewise Polynomial Estimation via Trend Filtering
Adaptive Piecewise Polynomial Estimation via Trend Filtering Liubo Li, ShanShan Tu The Ohio State University li.2201@osu.edu, tu.162@osu.edu October 1, 2015 Liubo Li, ShanShan Tu (OSU) Trend Filtering
More informationOn Solvability of an Indefinite Riccati Equation
Working Paper Series National Centre of Competence in Research Financial Valuation and Risk Management Working Paper No. 883 On Solvability of an Indefinite Riccati Equation Kai Du First version: December
More informationarxiv: v3 [math.ds] 22 Feb 2012
Stability of interconnected impulsive systems with and without time-delays using Lyapunov methods arxiv:1011.2865v3 [math.ds] 22 Feb 2012 Sergey Dashkovskiy a, Michael Kosmykov b, Andrii Mironchenko b,
More informationBrownian Motion. 1 Definition Brownian Motion Wiener measure... 3
Brownian Motion Contents 1 Definition 2 1.1 Brownian Motion................................. 2 1.2 Wiener measure.................................. 3 2 Construction 4 2.1 Gaussian process.................................
More informationPythagorean-hodograph curves
1 / 24 Pythagorean-hodograph curves V. Vitrih Raziskovalni matematični seminar 20. 2. 2012 2 / 24 1 2 3 4 5 3 / 24 Let r : [a, b] R 2 be a planar polynomial parametric curve ( ) x(t) r(t) =, y(t) where
More informationContinuous Curvature Path Generation Based on Bézier Curves for Autonomous Vehicles
Continuous Curvature Path Generation Based on Bézier Curves for Autonomous Vehicles Ji-wung Choi, Renwick E. Curry, Gabriel Hugh Elkaim Abstract In this paper we present two path planning algorithms based
More informationLimit solutions for control systems
Limit solutions for control systems M. Soledad Aronna Escola de Matemática Aplicada, FGV-Rio Oktobermat XV, 19 e 20 de Outubro de 2017, PUC-Rio 1 Introduction & motivation 2 Limit solutions 3 Commutative
More informationMA 351 Fall 2007 Exam #1 Review Solutions 1
MA 35 Fall 27 Exam # Review Solutions THERE MAY BE TYPOS in these solutions. Please let me know if you find any.. Consider the two surfaces ρ 3 csc θ in spherical coordinates and r 3 in cylindrical coordinates.
More informationNonlinear Observers. Jaime A. Moreno. Eléctrica y Computación Instituto de Ingeniería Universidad Nacional Autónoma de México
Nonlinear Observers Jaime A. Moreno JMorenoP@ii.unam.mx Eléctrica y Computación Instituto de Ingeniería Universidad Nacional Autónoma de México XVI Congreso Latinoamericano de Control Automático October
More informationQuadratic estimates and perturbations of Dirac type operators on doubling measure metric spaces
Quadratic estimates and perturbations of Dirac type operators on doubling measure metric spaces Lashi Bandara maths.anu.edu.au/~bandara Mathematical Sciences Institute Australian National University February
More information5.1 2D example 59 Figure 5.1: Parabolic velocity field in a straight two-dimensional pipe. Figure 5.2: Concentration on the input boundary of the pipe. The vertical axis corresponds to r 2 -coordinate,
More informationOn Irregular Linear Quadratic Control: Deterministic Case
1 On Irregular Linear Quadratic Control: Deterministic Case Huanshui Zhang and Juanjuan Xu arxiv:1711.9213v2 math.oc 1 Mar 218 Abstract This paper is concerned with fundamental optimal linear quadratic
More informationDEPARTMENT OF ELECTRICAL AND ELECTRONIC ENGINEERING EXAMINATIONS 2010
[E2.5] IMPERIAL COLLEGE LONDON DEPARTMENT OF ELECTRICAL AND ELECTRONIC ENGINEERING EXAMINATIONS 2010 EEE/ISE PART II MEng. BEng and ACGI SIGNALS AND LINEAR SYSTEMS Time allowed: 2:00 hours There are FOUR
More informationOptimal Linear Feedback Control for Incompressible Fluid Flow
Optimal Linear Feedback Control for Incompressible Fluid Flow Miroslav K. Stoyanov Thesis submitted to the Faculty of the Virginia Polytechnic Institute and State University in partial fulfillment of the
More informationLecture Note 13:Continuous Time Switched Optimal Control: Embedding Principle and Numerical Algorithms
ECE785: Hybrid Systems:Theory and Applications Lecture Note 13:Continuous Time Switched Optimal Control: Embedding Principle and Numerical Algorithms Wei Zhang Assistant Professor Department of Electrical
More informationKalman Filter and Parameter Identification. Florian Herzog
Kalman Filter and Parameter Identification Florian Herzog 2013 Continuous-time Kalman Filter In this chapter, we shall use stochastic processes with independent increments w 1 (.) and w 2 (.) at the input
More informationAn efficient approach to stochastic optimal control. Bert Kappen SNN Radboud University Nijmegen the Netherlands
An efficient approach to stochastic optimal control Bert Kappen SNN Radboud University Nijmegen the Netherlands Bert Kappen Examples of control tasks Motor control Bert Kappen Pascal workshop, 27-29 May
More informationINERTIALIZATION: HIGH-PERFORMANCE ANIMATION TRANSITIONS IN GEARS OF WAR DAVID BOLLO PRINCIPAL SOFTWARE ENGINEER THE COALITION (MICROSOFT STUDIOS)
INERTIALIZATION: HIGH-PERFORMANCE ANIMATION TRANSITIONS IN GEARS OF WAR DAVID BOLLO PRINCIPAL SOFTWARE ENGINEER THE COALITION (MICROSOFT STUDIOS) ANIMATION TRANSITIONS ANIMATION TRANSITIONS Transition
More informationTime-Invariant Linear Quadratic Regulators!
Time-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 17 Asymptotic approach from time-varying to constant gains Elimination of cross weighting
More informationReachability and Controllability
Capitolo. INTRODUCTION 4. Reachability and Controllability Reachability. The reachability problem is to find the set of all the final states x(t ) reachable starting from a given initial state x(t ) :
More informationCubic spline collocation for a class of weakly singular Fredholm integral equations and corresponding eigenvalue problem
Cubic spline collocation for a class of weakly singular Fredholm integral equations and corresponding eigenvalue problem ARVET PEDAS Institute of Mathematics University of Tartu J. Liivi 2, 549 Tartu ESTOIA
More information