On Linear-Quadratic Control Theory of Implicit Difference Equations
|
|
- Patrick Poole
- 6 years ago
- Views:
Transcription
1 On Linear-Quadratic Control Theory of Implicit Difference Equations Daniel Bankmann Technische Universität Berlin 10. Elgersburg Workshop 2016 February 9, 2016 D. Bankmann (TU Berlin) Control of IDEs February 9, / 18
2 Linear IDEs Consider the IDE (implicit difference equation) Eσx j = Ax j + B, Ex 0 = Ex 0, (1) σ denotes the shift operator, i. e., σx j = x j+1, (x j ) (K n ) N0, ( ) (K m ) N0 are some sequences. x 0 is some consistent shift variable. The IDE is assumed to be regular, i. e., det(ze A) 0. Given x 0, minimize J σ (x, u) := j=0 [ xj ] [ Q ] [ ] S xj S R subject to (1) with Q = Q, R = R and lim j Ex j = 0. D. Bankmann (TU Berlin) Control of IDEs February 9, / 18
3 Motivation Optimal control problems have a variety of applications, e. g., Trajectory-following problems, e. g., in Robotics Reaching a goal with minimal energy consumption Discrete-time models are obtained by discretization methods from the continuous-time models. D. Bankmann (TU Berlin) Control of IDEs February 9, / 18
4 Discretized Stokes Problem Solve [ ] M 0 σ 0 0 [ ] v = p [ K + M D T D 0 ] [ ] v + p [ ] B u, Mv 0 0 = Mv 0. v describes the velocity p describes the pressure M is the mass matrix K is the stiffness matrix Minimize, j=0 i. e., Q = 0, S = 0, R = I m, such that lim j Mv j = 0. D. Bankmann (TU Berlin) Control of IDEs February 9, / 18
5 Discrete Time Algebraic Riccati Equations Solve A XA X + Q (A XB + S)(B XB + R) 1 (B XA + S ) = 0 for Hermitian X fulfilling B XB + R 0. Under some circumstances invertibility of B XB + R is not guaranteed, e. g., when R is indefinite (e. g., differential games) Can be rewritten as [ A XA X + Q A ] [ ] XB + S K [K ] B XA + S B = L. XB + R L D. Bankmann (TU Berlin) Control of IDEs February 9, / 18
6 Discrete Time Algebraic Riccati Equations Solve A XA X + Q (A XB + S)(B XB + R) 1 (B XA + S ) = 0 for Hermitian X fulfilling B XB + R 0. Under some circumstances invertibility of B XB + R is not guaranteed, e. g., when R is indefinite (e. g., differential games) In the IDE case [ A XA E X E + Q A XB + S ] [ ] K [K ] L, B XA + S B XB + R = V Σ L where V Σ denotes the system space, i. e., the space where all solutions (x, u) evolve. D. Bankmann (TU Berlin) Control of IDEs February 9, / 18
7 Discrete Time Algebraic Riccati Equations Solve A XA X + Q (A XB + S)(B XB + R) 1 (B XA + S ) = 0 for Hermitian X fulfilling B XB + R 0. Under some circumstances invertibility of B XB + R is not guaranteed, e. g., when R is indefinite (e. g., differential games) In the IDE case [ (V Σ ) A XA E XE + Q A ] [ ] XB + S B XA + S B V Σ K [K ] = L, XB + R L where V Σ is a basis matrix of V Σ. D. Bankmann (TU Berlin) Control of IDEs February 9, / 18
8 Discrete-Time KYP Lemma Solve (discrete-time) Kalman-Yakubovich-Popov (KYP) inequality [Kalman 1963; Popov 1961; Yakubovich 1962] [ A PA E PE + Q A ] PB + S B PA + S B PB + R V Σ 0, P = P Popov function [ (ze A) Φ(z) := 1 ] [ ] [ B Q S (ze A) 1 ] B S K m m (z), R I m where G(z) := G(z 1 ) I m D. Bankmann (TU Berlin) Control of IDEs February 9, / 18
9 Discrete-Time KYP Lemma Theorem (KYP Lemma for IDEs) Let (E, A, B, Q, S, R) Σ w m,n(k) and the system space V Σ be given with corresponding Popov function Φ(z) K m m (z). If there exists some P K n n that is a solution of the KYP inequality, then Φ(e iω ) 0 for all ω R with det(e iω E A) 0. If on the other hand (E, A, B) is R-controllable and Φ(e iω ) 0 for all ω R with det(e iω E A) 0, then there exists a solution P K n n of the KYP inequality. R-controllability R-controllability means controllability on the set of reachable states R. In the IDE case R K n. D. Bankmann (TU Berlin) Control of IDEs February 9, / 18
10 Lur e Equations For q := rk K(z) Φ(z) find X = X K n n, K K q n, and L K q m such that ] [ ] K [K ] L. [ A XA E XE + Q A XB + S B XA + S B XB + R = V Σ L If X is a solution of the KYP inequality, we can always find K K p n and L K p m solving the Lur e equation, p q. From now on E = I n, i. e., solve [ A XA X + Q A ] [ ] XB + S K [K ] B XA + S B = L. XB + R L D. Bankmann (TU Berlin) Control of IDEs February 9, / 18
11 Construction of Deflating Subspace Consider the pencil 0 I n 0 0 A B ze A = z A 0 0 I n Q S K 2n+m 2n+m [z]. B S R Find deflating subspace, i. e., Y K 2n+m n+m, Z K 2n+m n+q, Ẽ, Ã K n+q n+m with rk K(z) (zẽ Ã) = n + q such that (ze A)Y = Z(zẼ Ã). D. Bankmann (TU Berlin) Control of IDEs February 9, / 18
12 Construction of Deflating Subspace Consider the pencil 0 I n 0 0 A B ze A = z A 0 0 I n Q S K 2n+m 2n+m [z]. B S R Find deflating subspace, i. e., Y K 2n+m n+m, Z K 2n+m n+q, Ẽ, Ã K n+q n+m with rk K(z) (zẽ Ã) = n + q such that (ze A)Y = Z(zẼ Ã). Ansatz: Y µ X 0 Z µ I n 0 Y = Y x = I n 0, Z = Z x = A X K, Y u 0 I m Z u B X L and [ ] zin A B zẽ Ã =. K L D. Bankmann (TU Berlin) Control of IDEs February 9, / 18
13 Application to Infinite Horizon Optimal Control Suppose, we have a solution (X, K, L) of the Lur e equation, j 2 j 1 and optimal (x, u) with lim j x j = 0. Then x j 2 Xx j2 x j 1 Xx j1 j 2 1 = σ(xk Xx k) xk Xx k, σx k = x k+1 k=j 1 j 2 1 [ ] [ xj A = XA X A ] [ ] XB xj B XA B XB k=j 1 D. Bankmann (TU Berlin) Control of IDEs February 9, / 18
14 Application to Infinite Horizon Optimal Control Suppose, we have a solution (X, K, L) of the Lur e equation, j 2 j 1 and optimal (x, u) with lim j x j = 0. Then x0 Xx 0 = σ(xk Xx k) xk Xx k, σx k = x k+1 = k=0 [ xj k=0 ] [ A XA X A ] [ ] XB xj B XA B XB D. Bankmann (TU Berlin) Control of IDEs February 9, / 18
15 Application to Infinite Horizon Optimal Control Suppose, we have a solution (X, K, L) of the Lur e equation, j 2 j 1 and optimal (x, u) with lim j x j = 0. Then x0 Xx 0 = σ(xk Xx k) xk Xx k, σx k = x k+1 = = = k=0 [ xj k=0 [ xj k=0 [ xj k=0 ] [ A XA X A ] [ ] XB xj B XA B XB ] ([ A XA X + Q A XB + S B XA + S B XB + R ] [ ] K [K ] [ x L j L ] [ ]) [ ] Q S xj S R ] J σ (x, u) = Kx + Lu 2 l J σ (x, u) 2 D. Bankmann (TU Berlin) Control of IDEs February 9, / 18
16 Application to Infinite Horizon Optimal Control Suppose, optimal (x, u) fulfills Kx + Lu 2 l 2 = 0. Using the deflating subspace it holds that 0 σi n A B X 0 [ ] σa I n Q S I n 0 xj σb S u R 0 I j m [ ] [ ] σin A B xj = Z, Z K 2n+m n+q. K L D. Bankmann (TU Berlin) Control of IDEs February 9, / 18
17 Application to Infinite Horizon Optimal Control Suppose, optimal (x, u) fulfills Kx + Lu 2 l 2 = 0. Using the deflating subspace it holds that 0 σi n A B X 0 [ ] σa I n Q S I n 0 xj σb S u R 0 I j m [ ] [ ] σin A B xj = Z = 0, Z K 2n+m n+q. K L D. Bankmann (TU Berlin) Control of IDEs February 9, / 18
18 Application to Infinite Horizon Optimal Control Suppose, optimal (x, u) fulfills Kx + Lu 2 l 2 = 0. Using the deflating subspace it holds that 0 σi n A B X 0 [ ] σa I n Q S I n 0 xj σb S u R 0 I j m [ ] [ ] σin A B xj = Z = 0, Z K 2n+m n+q. K L Thus, µ j := Xx j solves the boundary value problem 0 I n 0 µ 0 A B µ A 0 0 σ x = I n Q S x, B 0 0 u 0 S R u x 0 = x 0, lim µ j = 0. j D. Bankmann (TU Berlin) Control of IDEs February 9, / 18
19 Palindromic Matrix Pencils Reformulation of the BVP [Schröder 2008] with yields m j := µ j µ j+1 0 I n 0 m 0 A B m A Q S σ x = I n Q S x, B S R u 0 S R u x 0 = x 0, m j = µ j with corresponding palindromic matrix pencil 0 I n 0 0 A B za A = z A Q S I n Q S K 2n+m 2n+m [z]. B S R 0 S R j=0 D. Bankmann (TU Berlin) Control of IDEs February 9, / 18
20 Deflating Subspace for Palindromic Matrix Pencils Find deflating subspace Y of 0 I n 0 0 A B za A = z A Q S I n Q S B S R 0 S R given a solution of [ A XA X + Q A ] [ ] XB + S K [K ] B XA + S B = L, XB + R L i. e., (za A)Y = Z(zẼ Ã) for Ẽ, Ã K n+q n+m. D. Bankmann (TU Berlin) Control of IDEs February 9, / 18
21 Deflating Subspace for Palindromic Matrix Pencils Find deflating subspace Y of 0 I n 0 0 A B za A = z A Q S I n Q S B S R 0 S R given a solution of [ A XA X + Q A ] [ ] XB + S K [K ] B XA + S B = L, XB + R L i. e., (za A)Y = Z(zẼ Ã) for Ẽ, Ã K n+q n+m. ] Y = does not work due to the variable transformation. [ X 0 I n 0 0 I m D. Bankmann (TU Berlin) Control of IDEs February 9, / 18
22 Deflating Subspace for Palindromic Matrix Pencils From µ j = Xx j = X [ I n 0 ] [ ] x j it follows that m j = µ j µ j+1 = X (x j+1 x j ) = X [ A I n B ] [ ] x j. D. Bankmann (TU Berlin) Control of IDEs February 9, / 18
23 Deflating Subspace for Palindromic Matrix Pencils From µ j = Xx j = X [ I n 0 ] [ ] x j it follows that m j = µ j µ j+1 = X (x j+1 x j ) = X [ A I n B ] [ ] x j. D. Bankmann (TU Berlin) Control of IDEs February 9, / 18
24 Deflating Subspace for Palindromic Matrix Pencils From µ j = Xx j = X [ I n 0 ] [ ] x j it follows that m j = µ j µ j+1 = X (x j+1 x j ) = X [ A I n B ] [ ] x j. Leads to the deflating subspace Y µ X (A I n ) XB Y = Y x = Y u, Z = I n 0 0 I m Z µ Z x Z u = I n 0 (I n A )X K, B X L and [ ] zẽ Ã = zin A B. (z 1)K (z 1)L D. Bankmann (TU Berlin) Control of IDEs February 9, / 18
25 Deflating Subspace for Palindromic Matrix Pencils Theorem Let (I n, A, B, Q, S, R) Σ w m,n(k) and consider the palindromic pencil za A. Further, let q = rk K(z) Φ(z) and assume that rk [ I n A B ] = n. Then the following are equivalent: There exists a solution of the (discrete time) Lur e equation. The KYP [ ] inequality [ has a solution and there exist Yµ Zµ Y =, Z = Y x Y u Z x Z u the matrix [ I n A ],Ẽ, Ã, such that (ze A)Y = Z(zẼ Ã) and B ] [ ] Y x has full row rank n; Y u the space Y = im Y is the largest space such that Y (A A)Y = 0. D. Bankmann (TU Berlin) Control of IDEs February 9, / 18
26 Deflating Subspace for IDEs In the IDE case find deflating subspace Y of 0 E 0 0 A B za A = z A Q S E Q S B S R 0 S R given a solution of [ A XA E XE + Q A XB + S ] [ ] K [K ] L, B XA + S B XB + R = V Σ L i. e., (za A)Y = Z(zẼ Ã) for Ẽ, Ã K n+q n+m. D. Bankmann (TU Berlin) Control of IDEs February 9, / 18
27 Deflating Subspace for IDEs In the IDE case find deflating subspace Y of 0 E 0 0 A B za A = z A Q S E Q S B S R 0 S R given a solution of [ A XA E XE + Q A XB + S ] [ ] K [K ] L, B XA + S B XB + R = V Σ L i. e., (za A)Y = Z(zẼ Ã) for Ẽ, Ã K n+q n+m. Part of the Solution: [ ] X (A E) + G1 XB + G Y = 2 V1 Σ V2 Σ, where im [ ] G 1 G 2 ker E, and [ V1 Σ V2 Σ ] [ ] x = u [ ] x, u [ ] x V u Σ. D. Bankmann (TU Berlin) Control of IDEs February 9, / 18
28 Concluding Remarks Conclusions We have seen the discrete time counterpart of the KYP Lemma for DAEs. We have seen the discrete time counterpart of Lur e equations and corresponding deflating subspaces for (explicit) difference equations. The results on Lur e equations for DAEs can also be adapted to Lur e equations on IDEs. D. Bankmann (TU Berlin) Control of IDEs February 9, / 18
29 Concluding Remarks Conclusions We have seen the discrete time counterpart of the KYP Lemma for DAEs. We have seen the discrete time counterpart of Lur e equations and corresponding deflating subspaces for (explicit) difference equations. The results on Lur e equations for DAEs can also be adapted to Lur e equations on IDEs. What was not shown The spectral structure of palindromic matrix pencils can be characterized w.r.t. the solvability of the KYP inequality. Certain discretizations display relations between continuous time results and their discrete time counterpart. D. Bankmann (TU Berlin) Control of IDEs February 9, / 18
30 Concluding Remarks Conclusions We have seen the discrete time counterpart of the KYP Lemma for DAEs. We have seen the discrete time counterpart of Lur e equations and corresponding deflating subspaces for (explicit) difference equations. The results on Lur e equations for DAEs can also be adapted to Lur e equations on IDEs. What was not shown The spectral structure of palindromic matrix pencils can be characterized w.r.t. the solvability of the KYP inequality. Certain discretizations display relations between continuous time results and their discrete time counterpart. Thank you for your attention! D. Bankmann (TU Berlin) Control of IDEs February 9, / 18
31 References Byers, R., D. S. Mackey, V. Mehrmann, and H. Xu (2008). Symplectic, BVD, and palindromic approaches to discrete-time control problems. Preprint. Kalman, R. E. (1963). Lyapunov functions for the problem of Lur e in automatic control. Proceedings of the National Academy of Sciences of the United States of America 49.2, p Mehrmann, V. (1991). The Autonomous Linear Quadratic Control Problem. Ed. by M. Thoma and A. Wyner. Vol Lecture Notes in Control and Information Sciences. Heidelberg: Springer. Popov, V. M. (1961). On absolute stability of non-linear automatic control systems. Automatika i Telemekhanika 12, pp Reis, T. (2011). Lur e equations and even matrix pencils. Linear Algebra and its Applications 434.1, pp Reis, T., O. Rendel, and M. Voigt (2015). The Kalman Yakubovich Popov inequality for differential-algebraic systems. Linear Algebra and its Applications 485, pp Schröder, C. (2008). Palindromic and even eigenvalue problems-analysis and numerical methods. Dissertation. Institut für Mathematik, Technische Universität Berlin. Voigt, M. (2015). On Linear-Quadratic Optimal Control and Robustness of Differential-Algebraic Systems. Doctoral Dissertation. Fakultät für Mathematik, Otto-von-Guericke-Universität Magdeburg. Berlin: Logos. Yakubovich, V. (1962). Solution of certain matrix inequalities in the stability theory of nonlinear control systems. 143, pp D. Bankmann (TU Berlin) Control of IDEs February 9, / 18
The Kalman-Yakubovich-Popov Lemma for Differential-Algebraic Equations with Applications
MAX PLANCK INSTITUTE Elgersburg Workshop Elgersburg February 11-14, 2013 The Kalman-Yakubovich-Popov Lemma for Differential-Algebraic Equations with Applications Timo Reis 1 Matthias Voigt 2 1 Department
More informationarxiv: v4 [math.oc] 27 Oct 2017
On linear-quadratic optimal control of implicit difference equations Daniel Bankmann a,1,, Matthias Voigt a, a Technische Universität Berlin, Institut für Mathematik, Sekretariat MA 4 5, Straße des 17.
More informationLINEAR-QUADRATIC OPTIMAL CONTROL OF DIFFERENTIAL-ALGEBRAIC SYSTEMS: THE INFINITE TIME HORIZON PROBLEM WITH ZERO TERMINAL STATE
LINEAR-QUADRATIC OPTIMAL CONTROL OF DIFFERENTIAL-ALGEBRAIC SYSTEMS: THE INFINITE TIME HORIZON PROBLEM WITH ZERO TERMINAL STATE TIMO REIS AND MATTHIAS VOIGT, Abstract. In this work we revisit the linear-quadratic
More informationPerturbation theory for eigenvalues of Hermitian pencils. Christian Mehl Institut für Mathematik TU Berlin, Germany. 9th Elgersburg Workshop
Perturbation theory for eigenvalues of Hermitian pencils Christian Mehl Institut für Mathematik TU Berlin, Germany 9th Elgersburg Workshop Elgersburg, March 3, 2014 joint work with Shreemayee Bora, Michael
More informationStability and Inertia Theorems for Generalized Lyapunov Equations
Published in Linear Algebra and its Applications, 355(1-3, 2002, pp. 297-314. Stability and Inertia Theorems for Generalized Lyapunov Equations Tatjana Stykel Abstract We study generalized Lyapunov equations
More informationContour integral solutions of Sylvester-type matrix equations
Contour integral solutions of Sylvester-type matrix equations Harald K. Wimmer Mathematisches Institut, Universität Würzburg, 97074 Würzburg, Germany Abstract The linear matrix equations AXB CXD = E, AX
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik On Combining Deflation and Iteration to Low-Rank Approximative Solution of Lur e Equations Federico Poloni and Timo Reis Nr. 2012-01 January 2012 ON COMBINING
More informationDefinite versus Indefinite Linear Algebra. Christian Mehl Institut für Mathematik TU Berlin Germany. 10th SIAM Conference on Applied Linear Algebra
Definite versus Indefinite Linear Algebra Christian Mehl Institut für Mathematik TU Berlin Germany 10th SIAM Conference on Applied Linear Algebra Monterey Bay Seaside, October 26-29, 2009 Indefinite Linear
More informationNonlinear palindromic eigenvalue problems and their numerical solution
Nonlinear palindromic eigenvalue problems and their numerical solution TU Berlin DFG Research Center Institut für Mathematik MATHEON IN MEMORIAM RALPH BYERS Polynomial eigenvalue problems k P(λ) x = (
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 informationIterative Solution of a Matrix Riccati Equation Arising in Stochastic Control
Iterative Solution of a Matrix Riccati Equation Arising in Stochastic Control Chun-Hua Guo Dedicated to Peter Lancaster on the occasion of his 70th birthday We consider iterative methods for finding the
More informationonly nite eigenvalues. This is an extension of earlier results from [2]. Then we concentrate on the Riccati equation appearing in H 2 and linear quadr
The discrete algebraic Riccati equation and linear matrix inequality nton. Stoorvogel y Department of Mathematics and Computing Science Eindhoven Univ. of Technology P.O. ox 53, 56 M Eindhoven The Netherlands
More informationOn Positive Real Lemma for Non-minimal Realization Systems
Proceedings of the 17th World Congress The International Federation of Automatic Control On Positive Real Lemma for Non-minimal Realization Systems Sadaaki Kunimatsu Kim Sang-Hoon Takao Fujii Mitsuaki
More informationStability of Parameter Adaptation Algorithms. Big picture
ME5895, UConn, Fall 215 Prof. Xu Chen Big picture For ˆθ (k + 1) = ˆθ (k) + [correction term] we haven t talked about whether ˆθ(k) will converge to the true value θ if k. We haven t even talked about
More informationModel order reduction of electrical circuits with nonlinear elements
Model order reduction of electrical circuits with nonlinear elements Andreas Steinbrecher and Tatjana Stykel 1 Introduction The efficient and robust numerical simulation of electrical circuits plays a
More informationStructured eigenvalue/eigenvector backward errors of matrix pencils arising in optimal control
Electronic Journal of Linear Algebra Volume 34 Volume 34 08) Article 39 08 Structured eigenvalue/eigenvector backward errors of matrix pencils arising in optimal control Christian Mehl Technische Universitaet
More informationSTRUCTURE PRESERVING DEFLATION OF INFINITE EIGENVALUES IN STRUCTURED PENCILS
Electronic Transactions on Numerical Analysis Volume 44 pp 1 24 215 Copyright c 215 ISSN 168 9613 ETNA STRUCTURE PRESERVING DEFLATION OF INFINITE EIGENVALUES IN STRUCTURED PENCILS VOLKER MEHRMANN AND HONGGUO
More informationAn Observation on the Positive Real Lemma
Journal of Mathematical Analysis and Applications 255, 48 49 (21) doi:1.16/jmaa.2.7241, available online at http://www.idealibrary.com on An Observation on the Positive Real Lemma Luciano Pandolfi Dipartimento
More informationThe Solvability Conditions for the Inverse Eigenvalue Problem of Hermitian and Generalized Skew-Hamiltonian Matrices and Its Approximation
The Solvability Conditions for the Inverse Eigenvalue Problem of Hermitian and Generalized Skew-Hamiltonian Matrices and Its Approximation Zheng-jian Bai Abstract In this paper, we first consider the inverse
More information6. Linear Quadratic Regulator Control
EE635 - Control System Theory 6. Linear Quadratic Regulator Control Jitkomut Songsiri algebraic Riccati Equation (ARE) infinite-time LQR (continuous) Hamiltonian matrix gain margin of LQR 6-1 Algebraic
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Numerical Solution Of Matrix Equations Arising In Linear-Quadratic Optimal Control Of Differential-Algebraic Equations Timo Reis, Olaf Rendel Nr. 216-36 December
More informationDiscrete Riccati equations and block Toeplitz matrices
Discrete Riccati equations and block Toeplitz matrices André Ran Vrije Universiteit Amsterdam Leonid Lerer Technion-Israel Institute of Technology Haifa André Ran and Leonid Lerer 1 Discrete algebraic
More informationISOLATED SEMIDEFINITE SOLUTIONS OF THE CONTINUOUS-TIME ALGEBRAIC RICCATI EQUATION
ISOLATED SEMIDEFINITE SOLUTIONS OF THE CONTINUOUS-TIME ALGEBRAIC RICCATI EQUATION Harald K. Wimmer 1 The set of all negative-semidefinite solutions of the CARE A X + XA + XBB X C C = 0 is homeomorphic
More informationStructured Krylov Subspace Methods for Eigenproblems with Spectral Symmetries
Structured Krylov Subspace Methods for Eigenproblems with Spectral Symmetries Fakultät für Mathematik TU Chemnitz, Germany Peter Benner benner@mathematik.tu-chemnitz.de joint work with Heike Faßbender
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 informationModel reduction of nonlinear circuit equations
Model reduction of nonlinear circuit equations Tatjana Stykel Technische Universität Berlin Joint work with T. Reis and A. Steinbrecher BIRS Workshop, Banff, Canada, October 25-29, 2010 T. Stykel. Model
More informationThe Newton-ADI Method for Large-Scale Algebraic Riccati Equations. Peter Benner.
The Newton-ADI Method for Large-Scale Algebraic Riccati Equations Mathematik in Industrie und Technik Fakultät für Mathematik Peter Benner benner@mathematik.tu-chemnitz.de Sonderforschungsbereich 393 S
More informationOptimal control of nonstructured nonlinear descriptor systems
Optimal control of nonstructured nonlinear descriptor systems TU Berlin DFG Research Center Institut für Mathematik MATHEON Workshop Elgersburg 19.02.07 joint work with Peter Kunkel Overview Applications
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 informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik The Kalman-Yakubovich-Popov inequality for differential-algebraic systems: existence of nonpositive solutions Timo Reis, Matthias Voigt Nr. 015-09 February
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 informationOSCILLATION AND SPECTRAL THEORY FOR LINEAR HAMILTONIAN SYSTEMS WITH NONLINEAR DEPENDENCE ON THE SPECTRAL PARAMETER
OSCILLATION AND SPECTRAL THEORY FOR LINEAR HAMILTONIAN SYSTEMS WITH NONLINEAR DEPENDENCE ON THE SPECTRAL PARAMETER Martin Bohner, Werner Kratz und Roman Simon Hilscher Preprint Series: 2011-09 Fakultät
More informationInput-to-state stability and interconnected Systems
10th Elgersburg School Day 1 Input-to-state stability and interconnected Systems Sergey Dashkovskiy Universität Würzburg Elgersburg, March 5, 2018 1/20 Introduction Consider Solution: ẋ := dx dt = ax,
More informationBlock companion matrices, discrete-time block diagonal stability and polynomial matrices
Block companion matrices, discrete-time block diagonal stability and polynomial matrices Harald K. Wimmer Mathematisches Institut Universität Würzburg D-97074 Würzburg Germany October 25, 2008 Abstract
More informationAdvanced Probability
Advanced Probability Perla Sousi October 10, 2011 Contents 1 Conditional expectation 1 1.1 Discrete case.................................. 3 1.2 Existence and uniqueness............................ 3 1
More informationReduced Basis Method for Parametric
1/24 Reduced Basis Method for Parametric H 2 -Optimal Control Problems NUMOC2017 Andreas Schmidt, Bernard Haasdonk Institute for Applied Analysis and Numerical Simulation, University of Stuttgart andreas.schmidt@mathematik.uni-stuttgart.de
More informationOn the Kalman-Yacubovich-Popov lemma and common Lyapunov solutions for matrices with regular inertia
On the Kalman-Yacubovich-Popov lemma and common Lyapunov solutions for matrices with regular inertia Oliver Mason, Robert Shorten and Selim Solmaz Abstract In this paper we extend the classical Lefschetz
More informationRank-one LMIs and Lyapunov's Inequality. Gjerrit Meinsma 4. Abstract. We describe a new proof of the well-known Lyapunov's matrix inequality about
Rank-one LMIs and Lyapunov's Inequality Didier Henrion 1;; Gjerrit Meinsma Abstract We describe a new proof of the well-known Lyapunov's matrix inequality about the location of the eigenvalues of a matrix
More informationModern Optimal Control
Modern Optimal Control Matthew M. Peet Arizona State University Lecture 21: Optimal Output Feedback Control connection is called the (lower) star-product of P and Optimal Output Feedback ansformation (LFT).
More informationKrylov Subspace Methods for the Evaluation of Matrix Functions. Applications and Algorithms
Krylov Subspace Methods for the Evaluation of Matrix Functions. Applications and Algorithms 2. First Results and Algorithms Michael Eiermann Institut für Numerische Mathematik und Optimierung Technische
More informationOn the choice of abstract projection vectors for second level preconditioners
On the choice of abstract projection vectors for second level preconditioners C. Vuik 1, J.M. Tang 1, and R. Nabben 2 1 Delft University of Technology 2 Technische Universität Berlin Institut für Mathematik
More informationMATH4406 (Control Theory) Unit 6: The Linear Quadratic Regulator (LQR) and Model Predictive Control (MPC) Prepared by Yoni Nazarathy, Artem
MATH4406 (Control Theory) Unit 6: The Linear Quadratic Regulator (LQR) and Model Predictive Control (MPC) Prepared by Yoni Nazarathy, Artem Pulemotov, September 12, 2012 Unit Outline Goal 1: Outline linear
More informationStabilität differential-algebraischer Systeme
Stabilität differential-algebraischer Systeme Caren Tischendorf, Universität zu Köln Elgersburger Arbeitstagung, 11.-14. Februar 2008 Tischendorf (Univ. zu Köln) Stabilität von DAEs Elgersburg, 11.-14.02.2008
More informationMathematical Optimisation, Chpt 2: Linear Equations and inequalities
Mathematical Optimisation, Chpt 2: Linear Equations and inequalities Peter J.C. Dickinson p.j.c.dickinson@utwente.nl http://dickinson.website version: 12/02/18 Monday 5th February 2018 Peter J.C. Dickinson
More informationEE C128 / ME C134 Feedback Control Systems
EE C128 / ME C134 Feedback Control Systems Lecture Additional Material Introduction to Model Predictive Control Maximilian Balandat Department of Electrical Engineering & Computer Science University of
More informationAn SVD-Like Matrix Decomposition and Its Applications
An SVD-Like Matrix Decomposition and Its Applications Hongguo Xu Abstract [ A matrix S C m m is symplectic if SJS = J, where J = I m I m. Symplectic matrices play an important role in the analysis and
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 informationSome recent results on controllability of coupled parabolic systems: Towards a Kalman condition
Some recent results on controllability of coupled parabolic systems: Towards a Kalman condition F. Ammar Khodja Clermont-Ferrand, June 2011 GOAL: 1 Show the important differences between scalar and non
More informationIntroduction to Nonlinear Control Lecture # 4 Passivity
p. 1/6 Introduction to Nonlinear Control Lecture # 4 Passivity È p. 2/6 Memoryless Functions ¹ y È Ý Ù È È È È u (b) µ power inflow = uy Resistor is passive if uy 0 p. 3/6 y y y u u u (a) (b) (c) Passive
More informationMATH 304 Linear Algebra Lecture 10: Linear independence. Wronskian.
MATH 304 Linear Algebra Lecture 10: Linear independence. Wronskian. Spanning set Let S be a subset of a vector space V. Definition. The span of the set S is the smallest subspace W V that contains S. If
More informationResearch Article Stabilization Analysis and Synthesis of Discrete-Time Descriptor Markov Jump Systems with Partially Unknown Transition Probabilities
Research Journal of Applied Sciences, Engineering and Technology 7(4): 728-734, 214 DOI:1.1926/rjaset.7.39 ISSN: 24-7459; e-issn: 24-7467 214 Maxwell Scientific Publication Corp. Submitted: February 25,
More informationThe norms can also be characterized in terms of Riccati inequalities.
9 Analysis of stability and H norms Consider the causal, linear, time-invariant system ẋ(t = Ax(t + Bu(t y(t = Cx(t Denote the transfer function G(s := C (si A 1 B. Theorem 85 The following statements
More informationAnalysis Preliminary Exam Workshop: Hilbert Spaces
Analysis Preliminary Exam Workshop: Hilbert Spaces 1. Hilbert spaces A Hilbert space H is a complete real or complex inner product space. Consider complex Hilbert spaces for definiteness. If (, ) : H H
More informationDifferential-Algebraic Equations Forum
Differential-Algebraic Equations Forum Editors-in-Chief Achim Ilchmann (TU Ilmenau, Ilmenau, Germany) Timo Reis (Universität Hamburg, Hamburg, Germany) Editorial Board Larry Biegler (Carnegie Mellon University,
More informationw T 1 w T 2. w T n 0 if i j 1 if i = j
Lyapunov Operator Let A F n n be given, and define a linear operator L A : C n n C n n as L A (X) := A X + XA Suppose A is diagonalizable (what follows can be generalized even if this is not possible -
More informationClosed-Loop Structure of Discrete Time H Controller
Closed-Loop Structure of Discrete Time H Controller Waree Kongprawechnon 1,Shun Ushida 2, Hidenori Kimura 2 Abstract This paper is concerned with the investigation of the closed-loop structure of a discrete
More informationControllability of switched differential-algebraic equations
Controllability of switched differential-algebraic equations Ferdinand Küsters Markus G-M Ruppert Stephan Trenn Fachbereich Mathematik Technische Universität Kaiserslautern Postfach 349 67653 Kaiserslautern
More informationPort-Hamiltonian Realizations of Linear Time Invariant Systems
Port-Hamiltonian Realizations of Linear Time Invariant Systems Christopher Beattie and Volker Mehrmann and Hongguo Xu December 16 215 Abstract The question when a general linear time invariant control
More informationRational Covariance Extension for Boundary Data and Positive Real Lemma with Positive Semidefinite Matrix Solution
Preprints of the 18th FAC World Congress Milano (taly) August 28 - September 2, 2011 Rational Covariance Extension for Boundary Data and Positive Real Lemma with Positive Semidefinite Matrix Solution Y
More information4F3 - Predictive Control
4F3 Predictive Control - Lecture 2 p 1/23 4F3 - Predictive Control Lecture 2 - Unconstrained Predictive Control Jan Maciejowski jmm@engcamacuk 4F3 Predictive Control - Lecture 2 p 2/23 References Predictive
More informationA Study of the Duality between Kalman Filters and LQR Problems
Purdue University Purdue e-pubs Department of Electrical and Computer Engineering Technical Reports Department of Electrical and Computer Engineering 11-3-2016 A Study of the Duality between Kalman Filters
More informationA Characterization of the Hurwitz Stability of Metzler Matrices
29 American Control Conference Hyatt Regency Riverfront, St Louis, MO, USA June -2, 29 WeC52 A Characterization of the Hurwitz Stability of Metzler Matrices Kumpati S Narendra and Robert Shorten 2 Abstract
More informationFinite dimensional indefinite inner product spaces and applications in Numerical Analysis
Finite dimensional indefinite inner product spaces and applications in Numerical Analysis Christian Mehl Technische Universität Berlin, Institut für Mathematik, MA 4-5, 10623 Berlin, Germany, Email: mehl@math.tu-berlin.de
More informationNumerical Treatment of Unstructured. Differential-Algebraic Equations. with Arbitrary Index
Numerical Treatment of Unstructured Differential-Algebraic Equations with Arbitrary Index Peter Kunkel (Leipzig) SDS2003, Bari-Monopoli, 22. 25.06.2003 Outline Numerical Treatment of Unstructured Differential-Algebraic
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 informationFinite Element Decompositions for Stable Time Integration of Flow Equations
MAX PLANCK INSTITUT August 13, 2015 Finite Element Decompositions for Stable Time Integration of Flow Equations Jan Heiland, Robert Altmann (TU Berlin) ICIAM 2015 Beijing DYNAMIK KOMPLEXER TECHNISCHER
More informationOn the -Sylvester Equation AX ± X B = C
Department of Applied Mathematics National Chiao Tung University, Taiwan A joint work with E. K.-W. Chu and W.-W. Lin Oct., 20, 2008 Outline 1 Introduction 2 -Sylvester Equation 3 Numerical Examples 4
More informationCyclic reduction and index reduction/shifting for a second-order probabilistic problem
Cyclic reduction and index reduction/shifting for a second-order probabilistic problem Federico Poloni 1 Joint work with Giang T. Nguyen 2 1 U Pisa, Computer Science Dept. 2 U Adelaide, School of Math.
More informationOn the linear quadratic optimization problems and associated Riccati equations for systems modeled by Ito linear differential equations
Innovativity in Modeling and Analytics Journal of Research vol. 1, 216, pp.13-33, http://imajor.info/ c imajor ISSN XXXX-XXXX On the linear quadratic optimization problems and associated Riccati equations
More informationRECURSIVE ESTIMATION AND KALMAN FILTERING
Chapter 3 RECURSIVE ESTIMATION AND KALMAN FILTERING 3. The Discrete Time Kalman Filter Consider the following estimation problem. Given the stochastic system with x k+ = Ax k + Gw k (3.) y k = Cx k + Hv
More informationControl, Stabilization and Numerics for Partial Differential Equations
Paris-Sud, Orsay, December 06 Control, Stabilization and Numerics for Partial Differential Equations Enrique Zuazua Universidad Autónoma 28049 Madrid, Spain enrique.zuazua@uam.es http://www.uam.es/enrique.zuazua
More informationNORMS ON SPACE OF MATRICES
NORMS ON SPACE OF MATRICES. Operator Norms on Space of linear maps Let A be an n n real matrix and x 0 be a vector in R n. We would like to use the Picard iteration method to solve for the following system
More informationLinear algebra properties of dissipative Hamiltonian descriptor systems
Linear algebra properties of dissipative Hamiltonian descriptor systems C. Mehl V. Mehrmann M. Wojtylak January 6, 28 Abstract A wide class of matrix pencils connected with dissipative Hamiltonian descriptor
More informationTheorem Let X be a symmetric solution of DR(X) = 0. Let X b be a symmetric approximation to Xand set V := X? X. b If R b = R + B T XB b and I + V B R
Initializing Newton's method for discrete-time algebraic Riccati equations using the buttery SZ algorithm Heike Fabender Zentrum fur Technomathematik Fachbereich - Mathematik und Informatik Universitat
More informationControl engineering sample exam paper - Model answers
Question Control engineering sample exam paper - Model answers a) By a direct computation we obtain x() =, x(2) =, x(3) =, x(4) = = x(). This trajectory is sketched in Figure (left). Note that A 2 = I
More informationEE363 homework 2 solutions
EE363 Prof. S. Boyd EE363 homework 2 solutions. Derivative of matrix inverse. Suppose that X : R R n n, and that X(t is invertible. Show that ( d d dt X(t = X(t dt X(t X(t. Hint: differentiate X(tX(t =
More informationContents 1 Introduction and Preliminaries 1 Embedding of Extended Matrix Pencils 3 Hamiltonian Triangular Forms 1 4 Skew-Hamiltonian/Hamiltonian Matri
Technische Universitat Chemnitz Sonderforschungsbereich 393 Numerische Simulation auf massiv parallelen Rechnern Peter Benner Ralph Byers Volker Mehrmann Hongguo Xu Numerical Computation of Deating Subspaces
More informationDOUBLING ALGORITHMS WITH PERMUTED LAGRANGIAN GRAPH BASES
DOUBLNG ALGORTHMS WTH PERMUTED LAGRANGAN GRAPH BASES VOLKER MEHRMANN AND FEDERCO POLON Abstract. We derive a new representation of Lagrangian subspaces in the form m Π T, X where Π is a symplectic matrix
More informationEML5311 Lyapunov Stability & Robust Control Design
EML5311 Lyapunov Stability & Robust Control Design 1 Lyapunov Stability criterion In Robust control design of nonlinear uncertain systems, stability theory plays an important role in engineering systems.
More informationOUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative methods ffl Krylov subspace methods ffl Preconditioning techniques: Iterative methods ILU
Preconditioning Techniques for Solving Large Sparse Linear Systems Arnold Reusken Institut für Geometrie und Praktische Mathematik RWTH-Aachen OUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative
More informationStability analysis of differential algebraic systems
Stability analysis of differential algebraic systems Volker Mehrmann TU Berlin, Institut für Mathematik with Vu Hoang Linh and Erik Van Vleck DFG Research Center MATHEON Mathematics for key technologies
More informationModern Optimal Control
Modern Optimal Control Matthew M. Peet Arizona State University Lecture 19: Stabilization via LMIs Optimization Optimization can be posed in functional form: min x F objective function : inequality constraints
More informationLecture 5: Control Over Lossy Networks
Lecture 5: Control Over Lossy Networks Yilin Mo July 2, 2015 1 Classical LQG Control The system: x k+1 = Ax k + Bu k + w k, y k = Cx k + v k x 0 N (0, Σ), w k N (0, Q), v k N (0, R). Information available
More informationSufficiency Conditions for Finite Escape Times in Systems of Quadratic Differential Equations
/. Inst. Maths Applies (1977) 19, 377-383 Sufficiency Conditions for Finite Escape Times in Systems of Quadratic Differential Equations W. M. GETZ AND D. H. JACOBSON National Research Institute for Mathematical
More informationx 3y 2z = 6 1.2) 2x 4y 3z = 8 3x + 6y + 8z = 5 x + 3y 2z + 5t = 4 1.5) 2x + 8y z + 9t = 9 3x + 5y 12z + 17t = 7
Linear Algebra and its Applications-Lab 1 1) Use Gaussian elimination to solve the following systems x 1 + x 2 2x 3 + 4x 4 = 5 1.1) 2x 1 + 2x 2 3x 3 + x 4 = 3 3x 1 + 3x 2 4x 3 2x 4 = 1 x + y + 2z = 4 1.4)
More informationBalanced realization and model order reduction for nonlinear systems based on singular value analysis
Balanced realization and model order reduction for nonlinear systems based on singular value analysis Kenji Fujimoto a, and Jacquelien M. A. Scherpen b a Department of Mechanical Science and Engineering
More informationRobust port-hamiltonian representations of passive systems
Robust port-hamiltonian representations of passive systems Christopher Beattie and Volker Mehrmann and Paul Van Dooren January 19, 2018 Abstract We discuss the problem of robust representations of stable
More informationFast Algorithms for SDPs derived from the Kalman-Yakubovich-Popov Lemma
Fast Algorithms for SDPs derived from the Kalman-Yakubovich-Popov Lemma Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University 8 September 2003 European Union RTN Summer School on Multi-Agent
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 informationEG4321/EG7040. Nonlinear Control. Dr. Matt Turner
EG4321/EG7040 Nonlinear Control Dr. Matt Turner EG4321/EG7040 [An introduction to] Nonlinear Control Dr. Matt Turner EG4321/EG7040 [An introduction to] Nonlinear [System Analysis] and Control Dr. Matt
More informationModel estimation through matrix equations in financial econometrics
Model estimation through matrix equations in financial econometrics Federico Poloni 1 Joint work with Giacomo Sbrana 2 1 Technische Universität Berlin (A. Von Humboldt postdoctoral fellow) 2 Rouen Business
More informationMathematics for Control Theory
Mathematics for Control Theory Outline of Dissipativity and Passivity Hanz Richter Mechanical Engineering Department Cleveland State University Reading materials Only as a reference: Charles A. Desoer
More informationLecture 4. Chapter 4: Lyapunov Stability. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University.
Lecture 4 Chapter 4: Lyapunov Stability Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 4 p. 1/86 Autonomous Systems Consider the autonomous system ẋ
More informationController Design for Robust Output Regulation of Regular Linear Systems
Controller Design for Robust Output Regulation of Regular Linear Systems L. Paunonen We present three dynamic error feedbac controllers for robust output regulation of regular linear systems. These controllers
More informationKrylov Subspace Methods for the Evaluation of Matrix Functions. Applications and Algorithms
Krylov Subspace Methods for the Evaluation of Matrix Functions. Applications and Algorithms 4. Monotonicity of the Lanczos Method Michael Eiermann Institut für Numerische Mathematik und Optimierung Technische
More informationSolution for Homework 5
Solution for Homework 5 ME243A/ECE23A Fall 27 Exercise 1 The computation of the reachable subspace in continuous time can be handled easily introducing the concepts of inner product, orthogonal complement
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis Volume 31, pp 306-330, 2008 Copyright 2008, ISSN 1068-9613 ETNA STRUCTURED POLYNOMIAL EIGENPROBLEMS RELATED TO TIME-DELAY SYSTEMS HEIKE FASSBENDER, D STEVEN
More informationSTRUCTURED CONDITION NUMBERS FOR INVARIANT SUBSPACES
STRUCTURED CONDITION NUMBERS FOR INVARIANT SUBSPACES RALPH BYERS AND DANIEL KRESSNER Abstract. Invariant subspaces of structured matrices are sometimes better conditioned with respect to structured perturbations
More informationControl Systems Design
ELEC4410 Control Systems Design Lecture 14: Controllability Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science Lecture 14: Controllability p.1/23 Outline
More informationLinear Quadratic Zero-Sum Two-Person Differential Games
Linear Quadratic Zero-Sum Two-Person Differential Games Pierre Bernhard To cite this version: Pierre Bernhard. Linear Quadratic Zero-Sum Two-Person Differential Games. Encyclopaedia of Systems and Control,
More information