A Generic Ellipsoid Abstract Domain for Linear Time Invariant Systems
|
|
- Austen Chambers
- 5 years ago
- Views:
Transcription
1 A Generic Ellipsoid Abstract Domain for Linear Time Invariant Systems Pierre Roux 1 Romain Jobredeaux 2 Pierre-Loïc Garoche 1 Éric Féron 2 April 18, ONERA Toulouse, France 2 Georgia Tech Atlanta, Georgia, USA 1 / 26
2 Control-Command Critical Systems Critical systems rely on a control command computation core: flight commands of an aircraft; attitude control of a satellite; control of a car engine;... N 2 / 26
3 Stability of Linear Systems Those are reactive systems 3 / 26
4 Stability of Linear Systems Those are reactive systems periodically computing output as a combination of inputs and internal state variables. 3 / 26
5 Stability of Linear Systems Those are reactive systems periodically computing output as a combination of inputs and internal state variables. Goal Assuming a bound on inputs, we need to bound the outputs. 3 / 26
6 Stability of Linear Systems Those are reactive systems periodically computing output as a combination of inputs and internal state variables. Goal Assuming a bound on inputs, we need to bound the outputs. (for control theorists: open loop stability) 3 / 26
7 Stability of Linear Systems Those are reactive systems periodically computing output as a combination of inputs and internal state variables. Goal Assuming a bound on inputs, we need to bound the outputs. (for control theorists: open loop stability) We will focus on the (very common) case of linear systems. 3 / 26
8 Stability of linear systems, example // DiscreteStateSpace_A: A real_t A[16] = { , , , , , , , , , , , , , , , }; // DiscreteStateSpace_A: B real_t B[8] = { , , , , , , , }; static void MIMO_update(int_T tid) { static real_t xnew[4]; xnew[0] = (A[0])*St[0] + (A[1])*St[1] + (A[2])*St[2] + (A[3])*St[3]; xnew[0] += (B[0])*INPUT[0] + (B[1])*INPUT[1]; xnew[1] = (A[4])*St[0] + (A[5])*St[1] + (A[6])*St[2] + (A[7])*St[3]; xnew[1] += (B[2])*INPUT[0] + (B[3])*INPUT[1]; xnew[2] = (A[8])*St[0] + (A[9])*St[1] + (A[10])*St[2] + (A[11])*St[3]; xnew[2] += (B[4])*INPUT[0] + (B[5])*INPUT[1]; xnew[3] = (A[12])* St[0] + (A[13])*St[1] + (A[14])*St[2] + (A[15])*St[3]; xnew[3] += (B[6])*INPUT[0] + (B[7])*INPUT[1]; (void) memcpy(st, xnew, sizeof(real_t)*4); } 4 / 26
9 Stability of linear systems, example // DiscreteStateSpace_A: A real_t A[16] = { , , , , , , , , , , , , , , , }; // DiscreteStateSpace_A: B real_t B[8] = { , , , , , , , }; Result: if INPUT[0..1] remains in [ 1, 1] then St[0..3] remains in [-5, 5] static void MIMO_update(int_T tid) { static real_t xnew[4]; xnew[0] = (A[0])*St[0] + (A[1])*St[1] + (A[2])*St[2] + (A[3])*St[3]; xnew[0] += (B[0])*INPUT[0] + (B[1])*INPUT[1]; xnew[1] = (A[4])*St[0] + (A[5])*St[1] + (A[6])*St[2] + (A[7])*St[3]; xnew[1] += (B[2])*INPUT[0] + (B[3])*INPUT[1]; xnew[2] = (A[8])*St[0] + (A[9])*St[1] + (A[10])*St[2] + (A[11])*St[3]; xnew[2] += (B[4])*INPUT[0] + (B[5])*INPUT[1]; xnew[3] = (A[12])* St[0] + (A[13])*St[1] + (A[14])*St[2] + (A[15])*St[3]; xnew[3] += (B[6])*INPUT[0] + (B[7])*INPUT[1]; (void) memcpy(st, xnew, sizeof(real_t)*4); } 4 / 26
10 Quadratic invariants Linear invariants commonly used in static analysis are not well suited: at best costly; at worst no result. 5 / 26
11 Quadratic invariants Linear invariants commonly used in static analysis are not well suited: at best costly; at worst no result. y x 5 / 26
12 Quadratic invariants Linear invariants commonly used in static analysis are not well suited: at best costly; at worst no result. y x 5 / 26
13 Quadratic invariants Linear invariants commonly used in static analysis are not well suited: at best costly; at worst no result. Control theorists know for long that quadratic invariants are a good fit for linear systems. y y x x 5 / 26
14 Quadratic invariants Linear invariants commonly used in static analysis are not well suited: at best costly; at worst no result. Control theorists know for long that quadratic invariants are a good fit for linear systems. y y x x 5 / 26
15 1 Overall Method 2 Shape of the Ellipsoid 3 Floating Point Issues 4 Experimental Results 5 Related Work and Perspectives 6 / 26
16 1 Overall Method 2 Shape of the Ellipsoid 3 Floating Point Issues 4 Experimental Results 5 Related Work and Perspectives 7 / 26
17 Overall Method Lyapunov Stability Theorem For A R n n, B R n p, the sequence { x0 R n x k+1 = Ax k + Bu k is bounded for all u (R p ) N such that for all k N, u k 1 if and only if there exist P R n n positive definite such that A T P A P 0 where P 0 means that for all x R n : x 0 x T Px < 0. 8 / 26
18 Overall Method Lyapunov Stability, Invariant Invariant ellipsoid Moreover, there exist a λ { > 0 such that } x remains in the ellipsoid x R n x T P x λ. In Computer Science Language The property x T P x λ is a loop invariant. 9 / 26
19 { x Overall Method Lyapunov Stability, Illustration } x T Px 1 { Ax } x T Px 1 {Ax k + Bu u 1} Ax k x k 10 / 26
20 Overall Method Method Overall Method 1 First determine the shape of the ellipsoid by choosing a matrix P; 2 then find the smallest possible ratio λ such that x T P x λ is an invariant. 11 / 26
21 Overall Method Tools Definition (Semidefinite Programming) Minimize a linear objective function of variables y i under constraint k A 0 + y i A i 0 i=1 where the A i are known matrices and P 0 means x T P x 0 for all vector x. 12 / 26
22 Overall Method Tools Definition (Semidefinite Programming) Minimize a linear objective function of variables y i under constraint k A 0 + y i A i 0 i=1 where the A i are known matrices and P 0 means x T P x 0 for all vector x. Remark Efficient tools do exist. 12 / 26
23 1 Overall Method 2 Shape of the Ellipsoid 3 Floating Point Issues 4 Experimental Results 5 Related Work and Perspectives 13 / 26
24 Shape of the Ellipsoid What are we looking for? We look for a positive definite matrix P such that A T P A P / 26
25 Shape of the Ellipsoid What are we looking for? We look for a positive definite matrix P such that A T P A P 0. There are lot of very different solutions / 26
26 Shape of the Ellipsoid What are we looking for? We look for a positive definite matrix P such that A T P A P 0. There are lot of very different solutions......but we will then look for a λ > 0 such that x T P x λ is an invariant and the ellipsoid { x x T P x λ } should be small enough. Example y y x x is better than 14 / 26
27 Shape of the Ellipsoid What are we looking for? We look for a positive definite matrix P such that A T P A P 0. There are lot of very different solutions......but we will then look for a λ > 0 such that x T P x λ is an invariant and the ellipsoid { x x T P x λ } should be small enough. Example y y x x is better than Tested three heuristics. 14 / 26
28 1 Overall Method 2 Shape of the Ellipsoid 3 Floating Point Issues 4 Experimental Results 5 Related Work and Perspectives 15 / 26
29 Floating Point Issues Rounding Errors Actual computations are carried out with floating point numbers leading to rounding errors. Example int i = 0; float x = 0; while (i < ) { x += 0.1; ++i; } printf("%.0f\n", x); 16 / 26
30 Floating Point Issues Rounding Errors Actual computations are carried out with floating point numbers leading to rounding errors. Example int i = 0; float x = 0; while (i < ) { x += 0.1; ++i; } printf("%.0f\n", x); Gives ! 16 / 26
31 Floating Point Issues Rounding Errors Actual computations are carried out with floating point numbers leading to rounding errors. Example int i = 0; float x = 0; while (i < ) { x += 0.1; ++i; } printf("%.0f\n", x); Gives ! We have to distinguish two problems: rounding errors in the analyzed program; and rounding errors in the analyzer itself; 16 / 26
32 Floating Point Issues Rounding Errors Actual computations are carried out with floating point numbers leading to rounding errors. Example int i = 0; float x = 0; while (i < ) { x += 0.1; ++i; } printf("%.0f\n", x); Gives ! We have to distinguish two problems: rounding errors in the analyzed program; and rounding errors in the analyzer itself; both are addressed efficiently. 16 / 26
33 1 Overall Method 2 Shape of the Ellipsoid 3 Floating Point Issues 4 Experimental Results 5 Related Work and Perspectives 17 / 26
34 Experimental Results Our Prototype OCaml front-end to Scilab (lmisolver function). Only 2.5 kloc. Available under GPL at 18 / 26
35 Experimental Results Experimental Results Shape Bounds Valid [140.4; 189.9] Ex [22.2; 26.5] n=2, 1 input 0.23 [16.2; 17.6] [18.1; 25.2; 24.3; 33.7] Ex [6.3; 7.7; 2.2; 3.4] n=4, 1 input 0.40 [1.7; 2.0; 2.2; 2.5] Ex. 3 lead-lag 0.07 controller 0.17 [36.2; 36.1] n=2, 1 input 0.20 [38.8; 20.3] Ex. 4 LQG 0.09 [1.2; 0.9; 0.5] regulator 0.19 [0.9; 0.9; 0.9] n=3, 1 input 0.24 [0.7; 0.4; 0.3] Analysis times (in s) and bounds compared for the three heuristics. 19 / 26
36 Experimental Results Experimental Results, continued Shape Bounds Valid. Ex. 5 coupled 0.09 [9.8; 8.9; 11.0; 16.8] mass system 0.24 [5.7; 5.6; 6.4; 10.1] n=4, 2 inputs 0.48 [5.0; 4.9; 4.8; 4.7] Ex. 6 Butterworth 0.10 [7.5; 8.7; 6.1; 7.0; 6.5] low-pass filter 0.32 [3.6; 5.0; 4.7; 8.1; 8.9] n=5, 1 input 0.78 [2.3; 1.1; 1.9; 2.0; 2.9] Ex. 7 Dampened 0.07 [1.7; 2.1] oscillator 0.15 [2.0; 2.0] 0.20 n=2, no input 0.27 [1.5; 1.5] Ex. 8 Harmonic 0.08 [1.5; 1.5] oscillator 0.24 [1.5; 1.5] 0.20 n=2, no input 0.15 [1.5; 1.5] Analysis times (in s) and bounds compared for the three heuristics. 20 / 26
37 Experimental Results Experimental Results, continued (a) Ex. 1 (b) Ex. 2 (c) Ex. 3 (d) Ex. 4 Figure: Comparison of obtained ellipsoids by the three methods from lighter to darker, plus a random simulation trace ((b) and (d), being of dimension greater than 2, are cuts along planes containing the origin and two vectors of the canonical base, to show how the three different templates compare together). 21 / 26
38 Experimental Results Experimental Results, end (a) Ex. 5 (b) Ex. 6 (c) Ex. 7 (d) Ex. 8 Figure: Comparison of obtained ellipsoids by the three methods from lighter to darker, plus a random simulation trace (a) and (b), being of dimension greater than 2, are cuts along planes containing the origin and two vectors of the canonical base, to show how the three different templates compare together). 22 / 26
39 1 Overall Method 2 Shape of the Ellipsoid 3 Floating Point Issues 4 Experimental Results 5 Related Work and Perspectives 23 / 26
40 Related Work and Perspectives Related Work [Feret04] Adresses a more restricted class of systems. + Extremely precise thanks to formal expansion (unrolling). + Handles actual C code. + Proved industrially useful. [AdjéGaubertGoubault10] and [GawlitzaSeidl10] (policy iteration) Need to be given template ellipsoids. Does not address floating point issues. + Adresses a broader class of systems. + Can handle actual code. 24 / 26
41 Related Work and Perspectives Perspectives Adapt policy iteration. Include this new domain in our Lustre static analyzer. 25 / 26
42 Related Work and Perspectives Questions Thank you for your attention!? 26 / 26
43 Related Work and Perspectives Minimizing Condition Number We can look for the roundest possible ellipsoid by minimizing a new variable r in the extra constraint I P ri. Rationale: flat ellipsoids can yield a very bad bound on one of the variables. y x 27 / 31
44 Related Work and Perspectives Preserving the Shape Another heuristic is to minimize r (0, 1) in A T P A rp 0. Intuitively, the shape of ellipsoid that gets preserved the best through the update x k+1 = Ax k. y rp P x 28 / 31
45 Related Work and Perspectives All in One The two previous heuristics did not take B into account. 29 / 31
46 Related Work and Perspectives All in One The two previous heuristics did not take B into account. We could look for P maximizing r in P ri such that ( A T PA A T ) ( ) PBe i P 0 ei T B T PA 1 ei T B T τ PBe i 0 i 0 1 for all the vertices e i of the hypercube of dimension p. 29 / 31
47 Related Work and Perspectives All in One The two previous heuristics did not take B into account. We could look for P maximizing r in P ri such that ( A T PA A T ) ( ) PBe i P 0 ei T B T PA 1 ei T B T τ PBe i 0 i 0 1 for all the vertices e i of the hypercube of dimension p. + Previous equations amounts to the smallest ellipsoid such that x, u, u 1 x T P x 1 (Ax + Bu) T P (Ax + Bu) 1, i.e. theoretically the best solution. 29 / 31
48 Related Work and Perspectives All in One The two previous heuristics did not take B into account. We could look for P maximizing r in P ri such that ( A T PA A T ) ( ) PBe i P 0 ei T B T PA 1 ei T B T τ PBe i 0 i 0 1 for all the vertices e i of the hypercube of dimension p. + Previous equations amounts to the smallest ellipsoid such that x, u, u 1 x T P x 1 (Ax + Bu) T P (Ax + Bu) 1, i.e. theoretically the best solution. But in practice, it s hard to find a good solution. 29 / 31
49 Related Work and Perspectives Rounding Errors in the Program Knowing the precision of the floating point system used and the dimensions of matrices A and B of the analyzed system, we can compute two reals a and b such that if then (Ax + Bu) T P (Ax + Bu) λ fl(ax + Bu) T P fl(ax + Bu) a 2 λ + 2ab λ + b 2 with fl(e) the computation of e in any order and with any IEEE754 rounding mode (in practice a is near from 1 and b from 0). 30 / 31
50 Related Work and Perspectives Soundness of the result Checking the soundness of the result basically amounts to checking positive definiteness of a matrix. This is done by carefully bounding the rounding errors in a Cholesky decomposition. Hence an efficient soundness check (in O ( n 3) for an n n matrix). 31 / 31
Abstract Interpretation with Higher-Dimensional Ellipsoids and Conic Extrapolation
Abstract Interpretation with Higher-Dimensional Ellipsoids and Conic Extrapolation or potatoes & ice cream cones Mendes Oulamara, Arnaud Venet Computer Aided Verification, 2015 July 22, 2015 ÉCOLE NORMALE
More informationSparsity-Preserving Difference of Positive Semidefinite Matrix Representation of Indefinite Matrices
Sparsity-Preserving Difference of Positive Semidefinite Matrix Representation of Indefinite Matrices Jaehyun Park June 1 2016 Abstract We consider the problem of writing an arbitrary symmetric matrix as
More informationFormally Analyzing Adaptive Flight Control
Formally Analyzing Adaptive Flight Control Ashish Tiwari SRI International 333 Ravenswood Ave Menlo Park, CA 94025 Supported in part by NASA IRAC NRA grant number: NNX08AB95A Ashish Tiwari Symbolic Verification
More informationProving Unsatisfiability in Non-linear Arithmetic by Duality
Proving Unsatisfiability in Non-linear Arithmetic by Duality [work in progress] Daniel Larraz, Albert Oliveras, Enric Rodríguez-Carbonell and Albert Rubio Universitat Politècnica de Catalunya, Barcelona,
More informationLyapunov Stability Analysis: Open Loop
Copyright F.L. Lewis 008 All rights reserved Updated: hursday, August 8, 008 Lyapunov Stability Analysis: Open Loop We know that the stability of linear time-invariant (LI) dynamical systems can be determined
More informationEE221A Linear System Theory Final Exam
EE221A Linear System Theory Final Exam Professor C. Tomlin Department of Electrical Engineering and Computer Sciences, UC Berkeley Fall 2016 12/16/16, 8-11am Your answers must be supported by analysis,
More informationLecture 10 Linear Quadratic Stochastic Control with Partial State Observation
EE363 Winter 2008-09 Lecture 10 Linear Quadratic Stochastic Control with Partial State Observation partially observed linear-quadratic stochastic control problem estimation-control separation principle
More informationSUCCESSIVE POLE SHIFTING USING SAMPLED-DATA LQ REGULATORS. Sigeru Omatu
SUCCESSIVE POLE SHIFING USING SAMPLED-DAA LQ REGULAORS oru Fujinaka Sigeru Omatu Graduate School of Engineering, Osaka Prefecture University, 1-1 Gakuen-cho, Sakai, 599-8531 Japan Abstract: Design of sampled-data
More informationBalanced Truncation 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 2004: MODEL REDUCTION Balanced Truncation This lecture introduces balanced truncation for LTI
More informationTowards formal validation of aerospace systems
Towards formal validation of aerospace systems Eric Feron Dutton/Ducoffe Professor of Aerospace Engineering, Georgia Institute of Technology feron@gatech.edu Take-Home Message It changes every time I give
More informationTransformation of a PID Controller for Numerical Accuracy
Replace this file with prentcsmacro.sty for your meeting, or with entcsmacro.sty for your meeting. Both can be found at the ENTCS Macro Home Page. Transformation of a PID Controller for Numerical Accuracy
More informationCSCI 1951-G Optimization Methods in Finance Part 10: Conic Optimization
CSCI 1951-G Optimization Methods in Finance Part 10: Conic Optimization April 6, 2018 1 / 34 This material is covered in the textbook, Chapters 9 and 10. Some of the materials are taken from it. Some of
More informationAutomatic determination of numerical properties of software and systems
Automatic determination of numerical properties of software and systems Eric Goubault and Sylvie Putot Modelling and Analysis of Interacting Systems, CEA LIST MASCOT-NUM 2012 Meeting, March 21-23, 2012
More informationMath 411 Preliminaries
Math 411 Preliminaries Provide a list of preliminary vocabulary and concepts Preliminary Basic Netwon s method, Taylor series expansion (for single and multiple variables), Eigenvalue, Eigenvector, Vector
More informationEE363 homework 7 solutions
EE363 Prof. S. Boyd EE363 homework 7 solutions 1. Gain margin for a linear quadratic regulator. Let K be the optimal state feedback gain for the LQR problem with system ẋ = Ax + Bu, state cost matrix Q,
More informationScientific Computing with Case Studies SIAM Press, Lecture Notes for Unit VII Sparse Matrix
Scientific Computing with Case Studies SIAM Press, 2009 http://www.cs.umd.edu/users/oleary/sccswebpage Lecture Notes for Unit VII Sparse Matrix Computations Part 1: Direct Methods Dianne P. O Leary c 2008
More informationMODERN CONTROL DESIGN
CHAPTER 8 MODERN CONTROL DESIGN The classical design techniques of Chapters 6 and 7 are based on the root-locus and frequency response that utilize only the plant output for feedback with a dynamic controller
More informationPartially Observable Markov Decision Processes (POMDPs) Pieter Abbeel UC Berkeley EECS
Partially Observable Markov Decision Processes (POMDPs) Pieter Abbeel UC Berkeley EECS Many slides adapted from Jur van den Berg Outline POMDPs Separation Principle / Certainty Equivalence Locally Optimal
More informationMPC Infeasibility Handling
MPC Handling Thomas Wiese, TU Munich, KU Leuven supervised by H.J. Ferreau, Prof. M. Diehl (both KUL) and Dr. H. Gräb (TUM) October 9, 2008 1 / 42 MPC General MPC Strategies 2 / 42 Linear Discrete-Time
More information1. Find the solution of the following uncontrolled linear system. 2 α 1 1
Appendix B Revision Problems 1. Find the solution of the following uncontrolled linear system 0 1 1 ẋ = x, x(0) =. 2 3 1 Class test, August 1998 2. Given the linear system described by 2 α 1 1 ẋ = x +
More informationMCE/EEC 647/747: Robot Dynamics and Control. Lecture 8: Basic Lyapunov Stability Theory
MCE/EEC 647/747: Robot Dynamics and Control Lecture 8: Basic Lyapunov Stability Theory Reading: SHV Appendix Mechanical Engineering Hanz Richter, PhD MCE503 p.1/17 Stability in the sense of Lyapunov A
More informationLearning Methods for Linear Detectors
Intelligent Systems: Reasoning and Recognition James L. Crowley ENSIMAG 2 / MoSIG M1 Second Semester 2011/2012 Lesson 20 27 April 2012 Contents Learning Methods for Linear Detectors Learning Linear Detectors...2
More informationMiscellaneous. Regarding reading materials. Again, ask questions (if you have) and ask them earlier
Miscellaneous Regarding reading materials Reading materials will be provided as needed If no assigned reading, it means I think the material from class is sufficient Should be enough for you to do your
More informationarxiv: v1 [cs.lo] 24 Mar 2015
Property-based Polynomial Invariant Generation using Sums-of-Squares Optimization Assalé Adjé 1,a and Pierre-Loïc Garoche 1,a and Victor Magron 2,b arxiv:1503.07025v1 [cs.lo] 24 Mar 2015 1 Onera, the French
More informationRobust Control 5 Nominal Controller Design Continued
Robust Control 5 Nominal Controller Design Continued Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University 4/14/2003 Outline he LQR Problem A Generalization to LQR Min-Max
More informationAS computer hardware technology advances, both
1 Best-Harmonically-Fit Periodic Task Assignment Algorithm on Multiple Periodic Resources Chunhui Guo, Student Member, IEEE, Xiayu Hua, Student Member, IEEE, Hao Wu, Student Member, IEEE, Douglas Lautner,
More informationUnderstanding IC3. Aaron R. Bradley. ECEE, CU Boulder & Summit Middle School. Understanding IC3 1/55
Understanding IC3 Aaron R. Bradley ECEE, CU Boulder & Summit Middle School Understanding IC3 1/55 Further Reading This presentation is based on Bradley, A. R. Understanding IC3. In SAT, June 2012. http://theory.stanford.edu/~arbrad
More informationforms Christopher Engström November 14, 2014 MAA704: Matrix factorization and canonical forms Matrix properties Matrix factorization Canonical forms
Christopher Engström November 14, 2014 Hermitian LU QR echelon Contents of todays lecture Some interesting / useful / important of matrices Hermitian LU QR echelon Rewriting a as a product of several matrices.
More informationCourse 495: Advanced Statistical Machine Learning/Pattern Recognition
Course 495: Advanced Statistical Machine Learning/Pattern Recognition Lecturer: Stefanos Zafeiriou Goal (Lectures): To present discrete and continuous valued probabilistic linear dynamical systems (HMMs
More informationAnalysis of Discrete-Time Systems
TU Berlin Discrete-Time Control Systems 1 Analysis of Discrete-Time Systems Overview Stability Sensitivity and Robustness Controllability, Reachability, Observability, and Detectabiliy TU Berlin Discrete-Time
More informationStatic Output Feedback Stabilisation with H Performance for a Class of Plants
Static Output Feedback Stabilisation with H Performance for a Class of Plants E. Prempain and I. Postlethwaite Control and Instrumentation Research, Department of Engineering, University of Leicester,
More informationLecture 2: Discrete-time Linear Quadratic Optimal Control
ME 33, U Berkeley, Spring 04 Xu hen Lecture : Discrete-time Linear Quadratic Optimal ontrol Big picture Example onvergence of finite-time LQ solutions Big picture previously: dynamic programming and finite-horizon
More informationLecture 5: September 15
10-725/36-725: Convex Optimization Fall 2015 Lecture 5: September 15 Lecturer: Lecturer: Ryan Tibshirani Scribes: Scribes: Di Jin, Mengdi Wang, Bin Deng Note: LaTeX template courtesy of UC Berkeley EECS
More informationAutonomous Mobile Robot Design
Autonomous Mobile Robot Design Topic: Guidance and Control Introduction and PID Loops Dr. Kostas Alexis (CSE) Autonomous Robot Challenges How do I control where to go? Autonomous Mobile Robot Design Topic:
More informationPrincipal Component Analysis
Principal Component Analysis Laurenz Wiskott Institute for Theoretical Biology Humboldt-University Berlin Invalidenstraße 43 D-10115 Berlin, Germany 11 March 2004 1 Intuition Problem Statement Experimental
More informationThe ϵ-capacity of a gain matrix and tolerable disturbances: Discrete-time perturbed linear systems
IOSR Journal of Mathematics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 11, Issue 3 Ver. IV (May - Jun. 2015), PP 52-62 www.iosrjournals.org The ϵ-capacity of a gain matrix and tolerable disturbances:
More informationLECTURE NOTE #10 PROF. ALAN YUILLE
LECTURE NOTE #10 PROF. ALAN YUILLE 1. Principle Component Analysis (PCA) One way to deal with the curse of dimensionality is to project data down onto a space of low dimensions, see figure (1). Figure
More informationConvex relaxation. In example below, we have N = 6, and the cut we are considering
Convex relaxation The art and science of convex relaxation revolves around taking a non-convex problem that you want to solve, and replacing it with a convex problem which you can actually solve the solution
More informationAdvanced Optimization
Advanced Optimization Lecture 3: 1: Randomized Algorithms for for Continuous Discrete Problems Problems November 22, 2016 Master AIC Université Paris-Saclay, Orsay, France Anne Auger INRIA Saclay Ile-de-France
More informationOptimal control and estimation
Automatic Control 2 Optimal control and estimation Prof. Alberto Bemporad University of Trento Academic year 2010-2011 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 2010-2011
More informationContents lecture 5. Automatic Control III. Summary of lecture 4 (II/II) Summary of lecture 4 (I/II) u y F r. Lecture 5 H 2 and H loop shaping
Contents lecture 5 Automatic Control III Lecture 5 H 2 and H loop shaping Thomas Schön Division of Systems and Control Department of Information Technology Uppsala University. Email: thomas.schon@it.uu.se,
More informationThe Steepest Descent Algorithm for Unconstrained Optimization
The Steepest Descent Algorithm for Unconstrained Optimization Robert M. Freund February, 2014 c 2014 Massachusetts Institute of Technology. All rights reserved. 1 1 Steepest Descent Algorithm The problem
More informationAcknowledgements. Control System. Tracking. CS122A: Embedded System Design 4/24/2007. A Simple Introduction to Embedded Control Systems (PID Control)
Acknowledgements A Simple Introduction to Embedded Control Systems (PID Control) The material in this lecture is adapted from: F. Vahid and T. Givargis, Embedded System Design A Unified Hardware/Software
More informationX 2 3. Derive state transition matrix and its properties [10M] 4. (a) Derive a state space representation of the following system [5M] 1
QUESTION BANK 6 SIDDHARTH GROUP OF INSTITUTIONS :: PUTTUR Siddharth Nagar, Narayanavanam Road 5758 QUESTION BANK (DESCRIPTIVE) Subject with Code :SYSTEM THEORY(6EE75) Year &Sem: I-M.Tech& I-Sem UNIT-I
More informationAnalysis of Discrete-Time Systems
TU Berlin Discrete-Time Control Systems TU Berlin Discrete-Time Control Systems 2 Stability Definitions We define stability first with respect to changes in the initial conditions Analysis of Discrete-Time
More informationAlgebra Workshops 10 and 11
Algebra Workshops 1 and 11 Suggestion: For Workshop 1 please do questions 2,3 and 14. For the other questions, it s best to wait till the material is covered in lectures. Bilinear and Quadratic Forms on
More informationParallel Singular Value Decomposition. Jiaxing Tan
Parallel Singular Value Decomposition Jiaxing Tan Outline What is SVD? How to calculate SVD? How to parallelize SVD? Future Work What is SVD? Matrix Decomposition Eigen Decomposition A (non-zero) vector
More informationMath123 Lecture 1. Dr. Robert C. Busby. Lecturer: Office: Korman 266 Phone :
Lecturer: Math1 Lecture 1 Dr. Robert C. Busby Office: Korman 66 Phone : 15-895-1957 Email: rbusby@mcs.drexel.edu Course Web Site: http://www.mcs.drexel.edu/classes/calculus/math1_spring0/ (Links are case
More informationDecentralized Control of Stochastic Systems
Decentralized Control of Stochastic Systems Sanjay Lall Stanford University CDC-ECC Workshop, December 11, 2005 2 S. Lall, Stanford 2005.12.11.02 Decentralized Control G 1 G 2 G 3 G 4 G 5 y 1 u 1 y 2 u
More informationEvaluation of different wind estimation methods in flight tests with a fixed-wing UAV
Evaluation of different wind estimation methods in flight tests with a fixed-wing UAV Julian Sören Lorenz February 5, 2018 Contents 1 Glossary 2 2 Introduction 3 3 Tested algorithms 3 3.1 Unfiltered Method
More informationESC794: Special Topics: Model Predictive Control
ESC794: Special Topics: Model Predictive Control Discrete-Time Systems Hanz Richter, Professor Mechanical Engineering Department Cleveland State University Discrete-Time vs. Sampled-Data Systems A continuous-time
More informationA conjecture on sustained oscillations for a closed-loop heat equation
A conjecture on sustained oscillations for a closed-loop heat equation C.I. Byrnes, D.S. Gilliam Abstract The conjecture in this paper represents an initial step aimed toward understanding and shaping
More informationThe Multivariate Gaussian Distribution [DRAFT]
The Multivariate Gaussian Distribution DRAFT David S. Rosenberg Abstract This is a collection of a few key and standard results about multivariate Gaussian distributions. I have not included many proofs,
More informationU.C. Berkeley CS294: Spectral Methods and Expanders Handout 11 Luca Trevisan February 29, 2016
U.C. Berkeley CS294: Spectral Methods and Expanders Handout Luca Trevisan February 29, 206 Lecture : ARV In which we introduce semi-definite programming and a semi-definite programming relaxation of sparsest
More informationNeural Map. Structured Memory for Deep RL. Emilio Parisotto
Neural Map Structured Memory for Deep RL Emilio Parisotto eparisot@andrew.cmu.edu PhD Student Machine Learning Department Carnegie Mellon University Supervised Learning Most deep learning problems are
More informationORF 363/COS 323 Final Exam, Fall 2018
Name: Princeton University ORF 363/COS 323 Final Exam, Fall 2018 January 16, 2018 Instructor: A.A. Ahmadi AIs: Dibek, Duan, Gong, Khadir, Mirabelli, Pumir, Tang, Yu, Zhang 1. Please write out and sign
More informationMulti-Model Adaptive Regulation for a Family of Systems Containing Different Zero Structures
Preprints of the 19th World Congress The International Federation of Automatic Control Multi-Model Adaptive Regulation for a Family of Systems Containing Different Zero Structures Eric Peterson Harry G.
More informationComputing (the smallest) fixed points of monotone maps arising in static analysis by AI
Computing (the smallest fixed points of monotone maps arising in static analysis by AI Joint work, over the last 4 years: Stéphane Gaubert 1, Assalé Adjé, Eric Goubault 3, Sylvie Putot, Alexandru Costan,
More informationProgram verification. 18 October 2017
Program verification 18 October 2017 Example revisited // assume(n>2); void partition(int a[], int n) { int pivot = a[0]; int lo = 1, hi = n-1; while (lo
More informationEEE582 Homework Problems
EEE582 Homework Problems HW. Write a state-space realization of the linearized model for the cruise control system around speeds v = 4 (Section.3, http://tsakalis.faculty.asu.edu/notes/models.pdf). Use
More informationStochastic Models, Estimation and Control Peter S. Maybeck Volumes 1, 2 & 3 Tables of Contents
Navtech Part #s Volume 1 #1277 Volume 2 #1278 Volume 3 #1279 3 Volume Set #1280 Stochastic Models, Estimation and Control Peter S. Maybeck Volumes 1, 2 & 3 Tables of Contents Volume 1 Preface Contents
More informationGrothendieck s Inequality
Grothendieck s Inequality Leqi Zhu 1 Introduction Let A = (A ij ) R m n be an m n matrix. Then A defines a linear operator between normed spaces (R m, p ) and (R n, q ), for 1 p, q. The (p q)-norm of A
More informationRural/Urban Migration: The Dynamics of Eigenvectors
* Analysis of the Dynamic Structure of a System * Rural/Urban Migration: The Dynamics of Eigenvectors EGR 326 April 11, 2019 1. Develop the system model and create the Matlab/Simulink model 2. Plot and
More information14 Singular Value Decomposition
14 Singular Value Decomposition For any high-dimensional data analysis, one s first thought should often be: can I use an SVD? The singular value decomposition is an invaluable analysis tool for dealing
More informationNumber Theory. Zachary Friggstad. Programming Club Meeting
Number Theory Zachary Friggstad Programming Club Meeting Outline Factoring Sieve Multiplicative Functions Greatest Common Divisors Applications Chinese Remainder Theorem Throughout, problems to try are
More informationOn the Power of Robust Solutions in Two-Stage Stochastic and Adaptive Optimization Problems
MATHEMATICS OF OPERATIONS RESEARCH Vol. 35, No., May 010, pp. 84 305 issn 0364-765X eissn 156-5471 10 350 084 informs doi 10.187/moor.1090.0440 010 INFORMS On the Power of Robust Solutions in Two-Stage
More informationELEC4631 s Lecture 2: Dynamic Control Systems 7 March Overview of dynamic control systems
ELEC4631 s Lecture 2: Dynamic Control Systems 7 March 2011 Overview of dynamic control systems Goals of Controller design Autonomous dynamic systems Linear Multi-input multi-output (MIMO) systems Bat flight
More informationLU Factorization. LU factorization is the most common way of solving linear systems! Ax = b LUx = b
AM 205: lecture 7 Last time: LU factorization Today s lecture: Cholesky factorization, timing, QR factorization Reminder: assignment 1 due at 5 PM on Friday September 22 LU Factorization LU factorization
More informationTutorials in Optimization. Richard Socher
Tutorials in Optimization Richard Socher July 20, 2008 CONTENTS 1 Contents 1 Linear Algebra: Bilinear Form - A Simple Optimization Problem 2 1.1 Definitions........................................ 2 1.2
More informationConstrained optimization: direct methods (cont.)
Constrained optimization: direct methods (cont.) Jussi Hakanen Post-doctoral researcher jussi.hakanen@jyu.fi Direct methods Also known as methods of feasible directions Idea in a point x h, generate a
More informationLMI MODELLING 4. CONVEX LMI MODELLING. Didier HENRION. LAAS-CNRS Toulouse, FR Czech Tech Univ Prague, CZ. Universidad de Valladolid, SP March 2009
LMI MODELLING 4. CONVEX LMI MODELLING Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ Prague, CZ Universidad de Valladolid, SP March 2009 Minors A minor of a matrix F is the determinant of a submatrix
More informationSTABILITY ANALYSIS OF DAMPED SDOF SYSTEMS WITH TWO TIME DELAYS IN STATE FEEDBACK
Journal of Sound and Vibration (1998) 214(2), 213 225 Article No. sv971499 STABILITY ANALYSIS OF DAMPED SDOF SYSTEMS WITH TWO TIME DELAYS IN STATE FEEDBACK H. Y. HU ANDZ. H. WANG Institute of Vibration
More informationTOPOLOGICAL EQUIVALENCE OF LINEAR ORDINARY DIFFERENTIAL EQUATIONS
TOPOLOGICAL EQUIVALENCE OF LINEAR ORDINARY DIFFERENTIAL EQUATIONS ALEX HUMMELS Abstract. This paper proves a theorem that gives conditions for the topological equivalence of linear ordinary differential
More informationSemidefinite Programming Duality and Linear Time-invariant Systems
Semidefinite Programming Duality and Linear Time-invariant Systems Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University 2 July 2004 Workshop on Linear Matrix Inequalities in Control LAAS-CNRS,
More informationLQR, Kalman Filter, and LQG. Postgraduate Course, M.Sc. Electrical Engineering Department College of Engineering University of Salahaddin
LQR, Kalman Filter, and LQG Postgraduate Course, M.Sc. Electrical Engineering Department College of Engineering University of Salahaddin May 2015 Linear Quadratic Regulator (LQR) Consider a linear system
More informationAssignment #10: Diagonalization of Symmetric Matrices, Quadratic Forms, Optimization, Singular Value Decomposition. Name:
Assignment #10: Diagonalization of Symmetric Matrices, Quadratic Forms, Optimization, Singular Value Decomposition Due date: Friday, May 4, 2018 (1:35pm) Name: Section Number Assignment #10: Diagonalization
More information1 GSW Sets of Systems
1 Often, we have to solve a whole series of sets of simultaneous equations of the form y Ax, all of which have the same matrix A, but each of which has a different known vector y, and a different unknown
More informationLearning the Linear Dynamical System with ASOS ( Approximated Second-Order Statistics )
Learning the Linear Dynamical System with ASOS ( Approximated Second-Order Statistics ) James Martens University of Toronto June 24, 2010 Computer Science UNIVERSITY OF TORONTO James Martens (U of T) Learning
More informationConvex relaxation. In example below, we have N = 6, and the cut we are considering
Convex relaxation The art and science of convex relaxation revolves around taking a non-convex problem that you want to solve, and replacing it with a convex problem which you can actually solve the solution
More information8 A First Glimpse on Design with LMIs
8 A First Glimpse on Design with LMIs 8.1 Conceptual Design Problem Given a linear time invariant system design a linear time invariant controller or filter so as to guarantee some closed loop indices
More information1 An Overview and Brief History of Feedback Control 1. 2 Dynamic Models 23. Contents. Preface. xiii
Contents 1 An Overview and Brief History of Feedback Control 1 A Perspective on Feedback Control 1 Chapter Overview 2 1.1 A Simple Feedback System 3 1.2 A First Analysis of Feedback 6 1.3 Feedback System
More informationBasic Math for
Basic Math for 16-720 August 23, 2002 1 Linear Algebra 1.1 Vectors and Matrices First, a reminder of a few basic notations, definitions, and terminology: Unless indicated otherwise, vectors are always
More informationLecture Notes: Axiomatic Semantics and Hoare-style Verification
Lecture Notes: Axiomatic Semantics and Hoare-style Verification 17-355/17-665/17-819O: Program Analysis (Spring 2018) Claire Le Goues and Jonathan Aldrich clegoues@cs.cmu.edu, aldrich@cs.cmu.edu It has
More informationDO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED.
EE 537 Homewors Friedland Text Updated: Wednesday November 8 Some homewor assignments refer to Friedland s text For full credit show all wor. Some problems require hand calculations. In those cases do
More informationTemporal Modeling and Basic Speech Recognition
UNIVERSITY ILLINOIS @ URBANA-CHAMPAIGN OF CS 498PS Audio Computing Lab Temporal Modeling and Basic Speech Recognition Paris Smaragdis paris@illinois.edu paris.cs.illinois.edu Today s lecture Recognizing
More informationOn and Off-Policy Relational Reinforcement Learning
On and Off-Policy Relational Reinforcement Learning Christophe Rodrigues, Pierre Gérard, and Céline Rouveirol LIPN, UMR CNRS 73, Institut Galilée - Université Paris-Nord first.last@lipn.univ-paris13.fr
More informationCertified Roundoff Error Bounds using Semidefinite Programming
Certified Roundoff Error Bounds using Semidefinite Programming Victor Magron, CNRS VERIMAG joint work with G. Constantinides and A. Donaldson INRIA Mescal Team Seminar 19 November 2015 Victor Magron Certified
More informationAn Exact Stability Analysis Test for Single-Parameter. Polynomially-Dependent Linear Systems
An Exact Stability Analysis Test for Single-Parameter Polynomially-Dependent Linear Systems P. Tsiotras and P.-A. Bliman Abstract We provide a new condition for testing the stability of a single-parameter,
More informationAppendix I: Aircraft Flowchart Heuristic 1 pseudocode 1
1 Appendix I: Aircraft Flowchart Heuristic 1 pseudocode 1 Step 0: Initialization C res = C; Sa m1 = 0; Sy m1 = 0; Sg m1 = 0; a mn1 = A1 mn ; Sa m1 = Sa m1 + a mn1 ; C res = C res 1 + a mn1 ; y mn1 = Y
More informationANALYSIS OF LOAD PATTERNS IN RUBBER COMPONENTS FOR VEHICLES
ANALYSIS OF LOAD PATTERNS IN RUBBER COMPONENTS FOR VEHICLES Jerome Merel formerly at Hutchinson Corporate Research Center Israël Wander Apex Technologies Pierangelo Masarati, Marco Morandini Dipartimento
More informationHybrid Systems - Lecture n. 3 Lyapunov stability
OUTLINE Focus: stability of equilibrium point Hybrid Systems - Lecture n. 3 Lyapunov stability Maria Prandini DEI - Politecnico di Milano E-mail: prandini@elet.polimi.it continuous systems decribed by
More informationEconomic sensor/actuator selection and its application to flexible structure control
To appear in Proc. SPIE Int. Soc. Opt. Eng. 2004 Economic sensor/actuator selection and its application to flexible structure control Robert E. Skelton a and Faming Li a a Department of Mechanical and
More informationEE/ACM Applications of Convex Optimization in Signal Processing and Communications Lecture 4
EE/ACM 150 - Applications of Convex Optimization in Signal Processing and Communications Lecture 4 Andre Tkacenko Signal Processing Research Group Jet Propulsion Laboratory April 12, 2012 Andre Tkacenko
More informationCSE 241 Class 1. Jeremy Buhler. August 24,
CSE 41 Class 1 Jeremy Buhler August 4, 015 Before class, write URL on board: http://classes.engineering.wustl.edu/cse41/. Also: Jeremy Buhler, Office: Jolley 506, 314-935-6180 1 Welcome and Introduction
More informationSAMPLE SOLUTION TO EXAM in MAS501 Control Systems 2 Autumn 2015
FACULTY OF ENGINEERING AND SCIENCE SAMPLE SOLUTION TO EXAM in MAS501 Control Systems 2 Autumn 2015 Lecturer: Michael Ruderman Problem 1: Frequency-domain analysis and control design (15 pt) Given is a
More informationCHAPTER 11. A Revision. 1. The Computers and Numbers therein
CHAPTER A Revision. The Computers and Numbers therein Traditional computer science begins with a finite alphabet. By stringing elements of the alphabet one after another, one obtains strings. A set of
More informationNumerical Algorithms. IE 496 Lecture 20
Numerical Algorithms IE 496 Lecture 20 Reading for This Lecture Primary Miller and Boxer, Pages 124-128 Forsythe and Mohler, Sections 1 and 2 Numerical Algorithms Numerical Analysis So far, we have looked
More informationProgramming Languages and Compilers (CS 421)
Programming Languages and Compilers (CS 421) Sasa Misailovic 4110 SC, UIUC https://courses.engr.illinois.edu/cs421/fa2017/cs421a Based in part on slides by Mattox Beckman, as updated by Vikram Adve, Gul
More informationPhysicsAndMathsTutor.com. Advanced/Advanced Subsidiary. You must have: Mathematical Formulae and Statistical Tables (Blue)
Write your name here Surname Other names Pearson Edexcel International Advanced Level Centre Number Statistics S2 Advanced/Advanced Subsidiary Candidate Number Monday 27 June 2016 Morning Time: 1 hour
More informationFundamentals of Dynamical Systems / Discrete-Time Models. Dr. Dylan McNamara people.uncw.edu/ mcnamarad
Fundamentals of Dynamical Systems / Discrete-Time Models Dr. Dylan McNamara people.uncw.edu/ mcnamarad Dynamical systems theory Considers how systems autonomously change along time Ranges from Newtonian
More information