Linear Riccati Dynamics, Constant Feedback, and Controllability in Linear Quadratic Control Problems
|
|
- Stanley Booker
- 5 years ago
- Views:
Transcription
1 Linear Riccati Dynamics, Constant Feedback, and Controllability in Linear Quadratic Control Problems July 2001 Revised: December 2005 Ronald J. Balvers Douglas W. Mitchell Department of Economics Department of Economics P.O. Box 6025 P.O. Box 6025 West Virginia University West Virginia University Morgantown, WV Morgantown, WV Phone: (304) Phone: (304) ABSTRACT Conditions are derived for linear-quadratic control (LQC) problems to exhibit linear evolution of the Riccati matrix and constancy of the control feedback matrix. One of these conditions involves a matrix upon whose rank a necessary condition and a sufficient condition for controllability are based. Linearity of Riccati evolution allows for rapid iterative calculation, and constancy of the control feedback matrix allows for time-invariant comparative static analysis of policy reactions. JEL classification: C61 Keywords: Controllability, Riccati Equation, Linear Quadratic Control.
2 Linear Riccati Dynamics, Constant Feedback, and Controllability in Linear Quadratic Control Problems 1. Introduction This paper presents some conditions for linear evolution of the Riccati matrix, for constancy of the control-feedback matrix, and for controllability or lack thereof, in the standard linear quadratic control (LQC) problem. Computational aspects and simplification of the LQC problem have figured prominently in the recent literature on dynamic linear economies [see Amman (1997), Anderson et al. (1997), Anderson and Moore (1985), Binder and Pesaran (2000), Blanchard and Kahn (1980), Ehlgen (1999), Klein (2000), Ljungqvist and Sargent (2000), and Sims (2000)]. Linearity of Riccati evolution allows for rapid iterative calculation, and constancy of the control feedback matrix allows for time-invariant comparative static analysis of policy reactions. 2. The control problem A typical statement of the finite horizon linear quadratic control problem is: (1a), (1b) subject to, where K and A are n x n, C is n x k, y is n x 1, and u is k x 1. The cost matrix K is positive t definite, and C has full column rank; the transition matrix A need not have full rank. If the original problem statement has control costs, one can augment the state vector with the costed t
3 controls [see Chow (1975)], putting all costs on the state vector and thus giving the problem formulation in equations (1). It is well known [Chow (1975)] that the optimal controls are given by: (2), t T, (3),, t T, where the symmetric n x n matrix is positive definite. The nonlinear dynamic matrix equation (3) is the Riccati equation. The system is said to be controllable if and only if the state vector can be driven from any value to any value within n periods. The well-known necessary and sufficient condition for controllability [e.g., Anderson and Moore (1979), (1990), and Söderström (1994)] is that. 3. A Preliminary Lemma Write with C being q x k where q n - k, and with C being k x k. Since C 1 2 was assumed to be of full column rank, there is at least one k x k sub-matrix of C that is invertible. Proper prior arrangement of the vector (and concomitant arrangement of C, A, and K ) is thus sufficient to guarantee that C 2 is invertible. Define the n x q matrix M as, so that. 2
4 We can now prove the following identity: LEMMA (IDENTITY RE-EXPRESSING A FUNCTION OF C IN TERMS OF M). For C (which is n x k with full column rank) and M (which is n x q (q = n - k) with full column rank) such that, and D any (symmetric) positive definite n x n matrix, the following identity holds: (4). Proof. Consider a Cholesky decomposition such that, with full rank (such a decomposition exists for any positive definite matrix). Then we can premultiply equation (4) by, postmultiply by, and rearrange to obtain (5). Note that both sides of equation (5) express symmetric idempotent matrices. Define the n x n matrix. It is easy to confirm that. Hence, defining the n x n matrix, we have: (6). 3
5 Next show that Y has full rank: (i) Q has full rank n and C has full column rank k, so Y 1 = QC has rank k; (ii) Q has full rank n and M has full column rank q, so Y 2 = has rank q; (iii), so all the k columns of Y are 1 independent of all the q columns of Y2. Given (i), (ii), and (iii), Y has full rank. Thus, from equation (6), Z = 0 implying that by the definition of Z equation (5) must hold; but equation (5) is equivalent to equation (4) in the statement of the Lemma. 4. Conditions for Linear Riccati Evolution, Constant Feedback, and Controllability Apply the Lemma to the Riccati equation (3), setting which is in fact positive definite. This yields the following alternative statement of the Riccati equation: (7),. This allows us to obtain Theorem 1: THEOREM 1 (LINEAR DYNAMICS OF ). Linear evolution of, and constancy of the control feedback matrix for all, are both implied by each of (a), and (b). More strongly, (b) implies constancy of. Proof. (a) Post-multiply equation (7) by C. With this gives for all t T. Using this equation and its transpose in equations (2) and (3) gives 4
6 , and (8),. (b) Use a standard matrix inversion identity [e.g. Söderström (1994), pp ] on equation (7) to obtain (9),, and post-multiply this by M to obtain when. Using this in equation (7) gives which is constant for t T. Equation (2) with constant gives constant for t T. Linear evolution of in equation (8) gives obvious computational advantages. Constancy of the control feedback matrix implies that comparative static analysis of policy reactions can be conducted in straightforward fashion, with the results not dependent on time to horizon. This Theorem also implies a result for stabilizability [Ljungqvist and Sargent (2000), p.61]: 5
7 COROLLARY (SUFFICIENT CONDITION FOR STABILIZABILITY). is sufficient for stabilizability, and moreover, under this condition optimal stabilization drives the state variables to their target values of 0 in two periods. The Corollary follows trivially from the immediate constancy of. For period T the cost matrix is ; for period T-1 it equals its steady state value; and with more than two periods to the horizon, the cost is unchanged; thus is driven to its desired value 0 in two periods. For the latter to occur, two things are necessary: A and C, and hence M, must be such that it is feasible to do so; and K must be such that it is desirable to do so. The corollary shows that the condition meets both of these criteria. THEOREM 2 (RELATION OF TO CONTROLLABILITY). (a) Full row rank of implies controllability; (b) implies, but is not implied by, lack of controllability. Proof. (a) Consider a canonical form of the decision problem in equations (1) with such that and. Any possible form of the problem in equations (1) can be put in this canonical form by setting, where and. It is easy 6
8 to check that, full row rank (n) of which is equivalent to controllability (so controllability of the canonical problem is equivalent to controllability of the original problem). Now *, where we partition A as with being q x k. Thus we need to show that, if, the system is controllable. (Of course this requires q k -- equivalently, n 2k.) Let =. Clearly full row rank of * implies full row rank of W. So the canonical problem and hence the original problem are controllable. (b) We show that implies, where the latter defines noncontrollability. Consider again the canonical form referred to above. Since = which now equals zero, this directly implies lack of controllability for the canonical problem (sub-vector cannot be controlled directly since, nor indirectly by controlling since ). As shown above, lack of controllability of the canonical problem implies lack of controllability of the original problem. Part (b) further asserts that lack of controllability need not imply that. We show this by example: Let and take to be idempotent and such that. 7
9 Then the system is noncontrollable since 1, since C and AC each have k columns. Note that does not imply lack of stabilizability. We can show this by counterexample: let, in which case since by construction of M. If the state equation is stable in the absence of time-varying control, and hence trivially is stabilizable. 1 A numerical example is:. 8
10 References Amman, H. M., 1997, Numerical methods for linear-quadratic models, in H.M. Amman, D.A. Kendrick, and J. Rust, eds., Handbook of Computational Economics (North-Holland, Amsterdam). Anderson, B. D. A. and J. B. Moore, 1979.Optimal filtering (Prentice Hall, Englewood Cliffs). Anderson, B. D. A. and J. B. Moore, Optimal control: linear quadratic methods (Prentice Hall, Englewood Cliffs). Anderson, E., L. P. Hansen, E. R. McGrattan and T. Sargent, Mechanics of forming and estimating dynamical linear economies, in H.M. Amman, D.A. Kendrick, and J. Rust, eds., Handbook of Computational Economics (North-Holland, Amsterdam). Anderson, G. and G. Moore, A linear algebraic procedure for solving linear perfect foresight models, Economics Letters 17, Binder, M. and M. H. Pesaran, Solution of finite-horizon multivariate linear rational expectations models and sparse linear systems, Journal of Economic Dynamics and Control 24, Blanchard, O. J. and Kahn, C. M., The solution of linear difference models under rational expectations, Econometrica 48, Chow, G. C., Analysis and control of dynamic economic systems (Wiley, New York). Ehlgen, J., A nonrecursive solution method for the linear-quadratic optimal control problem with a singular transition matrix, Computational Economics 13, Klein, P., Using the generalized Schur form to solve a multivariate linear rational expectations model, Journal of Economic Dynamics and Control 24, Ljungqvist, L. and T. J. Sargent, Recursive macroeconomic theory (MIT Press, Cambridge). Sims, C. S., 2000, Solving linear rational expectations models, Working paper, Yale University. Söderström, T., Discrete-time stochastic systems: estimation and control (Prentice-Hall, New York). 9
Linear Riccati Dynamics, Constant Feedback, and Controllability in Linear Quadratic Control Problems
Linear Riccati Dynamics, Constant Feedback, and Controllability in Linear Quadratic Control Problems July 2001 Ronald J. Balvers Douglas W. Mitchell Department of Economics Department of Economics P.O.
More informationReducing the Dimensionality of Linear Quadratic Control Problems*
Reducing the Dimensionality of Linear Quadratic Control Problems* February 2001 Ronald J. Balvers Douglas W. Mitchell Department of Economics Department of Economics P.O. Box 6025 P.O. Box 6025 West Virginia
More informationVolume 30, Issue 3. A note on Kalman filter approach to solution of rational expectations models
Volume 30, Issue 3 A note on Kalman filter approach to solution of rational expectations models Marco Maria Sorge BGSE, University of Bonn Abstract In this note, a class of nonlinear dynamic models under
More informationA Simple Algorithm for Solving Ramsey Optimal Policy with Exogenous Forcing Variables
MPRA Munich Personal RePEc Archive A Simple Algorithm for Solving Ramsey Optimal Policy with Exogenous Forcing Variables Jean-Bernard Chatelain and Kirsten Ralf Paris School of Economics, Université Paris
More informationA Simple Algorithm for Solving Ramsey Optimal Policy with Exogenous Forcing Variables
A Simple Algorithm for Solving Ramsey Optimal Policy with Exogenous Forcing Variables Jean-Bernard Chatelain and Kirsten Ralf arxiv:1708.07996v1 [q-fin.ec] 26 Aug 2017 June 20, 2018 Abstract This algorithm
More informationSolving Linear Rational Expectation Models
Solving Linear Rational Expectation Models Dr. Tai-kuang Ho Associate Professor. Department of Quantitative Finance, National Tsing Hua University, No. 101, Section 2, Kuang-Fu Road, Hsinchu, Taiwan 30013,
More informationOptimization under Commitment and Discretion, the Recursive Saddlepoint Method, and Targeting Rules and Instrument Rules: Lecture Notes
CommDiscTRIR.tex Preliminary; comments welcome Optimization under Commitment and Discretion, the Recursive Saddlepoint Method, and Targeting Rules and Instrument Rules: Lecture Notes Lars E.O. Svensson
More informationMatrix Mathematics. Theory, Facts, and Formulas with Application to Linear Systems Theory. Dennis S. Bernstein
Matrix Mathematics Theory, Facts, and Formulas with Application to Linear Systems Theory Dennis S. Bernstein PRINCETON UNIVERSITY PRESS PRINCETON AND OXFORD Contents Special Symbols xv Conventions, Notation,
More informationEstimating and Identifying Vector Autoregressions Under Diagonality and Block Exogeneity Restrictions
Estimating and Identifying Vector Autoregressions Under Diagonality and Block Exogeneity Restrictions William D. Lastrapes Department of Economics Terry College of Business University of Georgia Athens,
More informationMeasurement Errors and the Kalman Filter: A Unified Exposition
Luiss Lab of European Economics LLEE Working Document no. 45 Measurement Errors and the Kalman Filter: A Unified Exposition Salvatore Nisticò February 2007 Outputs from LLEE research in progress, as well
More informationOPTIMAL CONTROL AND ESTIMATION
OPTIMAL CONTROL AND ESTIMATION Robert F. Stengel Department of Mechanical and Aerospace Engineering Princeton University, Princeton, New Jersey DOVER PUBLICATIONS, INC. New York CONTENTS 1. INTRODUCTION
More informationDecentralised economies I
Decentralised economies I Martin Ellison 1 Motivation In the first two lectures on dynamic programming and the stochastic growth model, we solved the maximisation problem of the representative agent. More
More informationNew Introduction to Multiple Time Series Analysis
Helmut Lütkepohl New Introduction to Multiple Time Series Analysis With 49 Figures and 36 Tables Springer Contents 1 Introduction 1 1.1 Objectives of Analyzing Multiple Time Series 1 1.2 Some Basics 2
More informationDepartment of Economics. Working Papers
10ISSN 1183-1057 SIMON FRASER UNIVERSITY Department of Economics Working Papers 08-07 "A Cauchy-Schwarz inequality for expectation of matrices" Pascal Lavergne November 2008 Economics A Cauchy-Schwarz
More informationMATRIX ALGEBRA AND SYSTEMS OF EQUATIONS. + + x 1 x 2. x n 8 (4) 3 4 2
MATRIX ALGEBRA AND SYSTEMS OF EQUATIONS SYSTEMS OF EQUATIONS AND MATRICES Representation of a linear system The general system of m equations in n unknowns can be written a x + a 2 x 2 + + a n x n b a
More informationMatrix Polynomial Conditions for the Existence of Rational Expectations.
Matrix Polynomial Conditions for the Existence of Rational Expectations. John Hunter Department of Economics and Finance Section, School of Social Science, Brunel University, Uxbridge, Middlesex, UB8 PH,
More informationLessons in Estimation Theory for Signal Processing, Communications, and Control
Lessons in Estimation Theory for Signal Processing, Communications, and Control Jerry M. Mendel Department of Electrical Engineering University of Southern California Los Angeles, California PRENTICE HALL
More informationAn Alternative Method to Obtain the Blanchard and Kahn Solutions of Rational Expectations Models
An Alternative Method to Obtain the Blanchard and Kahn Solutions of Rational Expectations Models Santiago Acosta Ormaechea Alejo Macaya August 8, 2007 Abstract In this paper we show that the method of
More informationRecursive Solutions of the Matrix Equations X + A T X 1 A = Q and X A T X 1 A = Q
Applied Mathematical Sciences, Vol. 2, 2008, no. 38, 1855-1872 Recursive Solutions of the Matrix Equations X + A T X 1 A = Q and X A T X 1 A = Q Nicholas Assimakis Department of Electronics, Technological
More informationIn: Proc. BENELEARN-98, 8th Belgian-Dutch Conference on Machine Learning, pp 9-46, 998 Linear Quadratic Regulation using Reinforcement Learning Stephan ten Hagen? and Ben Krose Department of Mathematics,
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 informationRiccati difference equations to non linear extended Kalman filter constraints
International Journal of Scientific & Engineering Research Volume 3, Issue 12, December-2012 1 Riccati difference equations to non linear extended Kalman filter constraints Abstract Elizabeth.S 1 & Jothilakshmi.R
More informationOn Solving Large Algebraic. Riccati Matrix Equations
International Mathematical Forum, 5, 2010, no. 33, 1637-1644 On Solving Large Algebraic Riccati Matrix Equations Amer Kaabi Department of Basic Science Khoramshahr Marine Science and Technology University
More informationProperties of Matrices and Operations on Matrices
Properties of Matrices and Operations on Matrices A common data structure for statistical analysis is a rectangular array or matris. Rows represent individual observational units, or just observations,
More informationUsing the generalized Schur form to solve a multivariate linear rational expectations model
Journal of Economic Dynamics & Control 24 (2000) 1405}1423 Using the generalized Schur form to solve a multivariate linear rational expectations model Paul Klein* Institute for International Economic Studies,
More informationMATRICES. knowledge on matrices Knowledge on matrix operations. Matrix as a tool of solving linear equations with two or three unknowns.
MATRICES After studying this chapter you will acquire the skills in knowledge on matrices Knowledge on matrix operations. Matrix as a tool of solving linear equations with two or three unknowns. List of
More informationThe matrix will only be consistent if the last entry of row three is 0, meaning 2b 3 + b 2 b 1 = 0.
) Find all solutions of the linear system. Express the answer in vector form. x + 2x + x + x 5 = 2 2x 2 + 2x + 2x + x 5 = 8 x + 2x + x + 9x 5 = 2 2 Solution: Reduce the augmented matrix [ 2 2 2 8 ] to
More informationSecond-order approximation of dynamic models without the use of tensors
Second-order approximation of dynamic models without the use of tensors Paul Klein a, a University of Western Ontario First draft: May 17, 2005 This version: January 24, 2006 Abstract Several approaches
More informationMacroeconomics Theory II
Macroeconomics Theory II Economies as ynamic Systems Francesco Franco Nova SBE February 11, 2014 Francesco Franco Macroeconomics Theory II 1/18 First-order The variable z t follows a first-order di erence
More informationChapter 6. Maximum Likelihood Analysis of Dynamic Stochastic General Equilibrium (DSGE) Models
Chapter 6. Maximum Likelihood Analysis of Dynamic Stochastic General Equilibrium (DSGE) Models Fall 22 Contents Introduction 2. An illustrative example........................... 2.2 Discussion...................................
More informationSolving a Non-Linear Model: The Importance of Model Specification for Deriving a Suitable Solution
Solving a Non-Linear Model: The Importance of Model Specification for Deriving a Suitable Solution Ric D. Herbert a and Peter J. Stemp b a School of Design, Communication and Information Technology, The
More informationIntroduction to Quantitative Techniques for MSc Programmes SCHOOL OF ECONOMICS, MATHEMATICS AND STATISTICS MALET STREET LONDON WC1E 7HX
Introduction to Quantitative Techniques for MSc Programmes SCHOOL OF ECONOMICS, MATHEMATICS AND STATISTICS MALET STREET LONDON WC1E 7HX September 2007 MSc Sep Intro QT 1 Who are these course for? The September
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 informationThe Kalman filter is arguably one of the most notable algorithms
LECTURE E NOTES «Kalman Filtering with Newton s Method JEFFREY HUMPHERYS and JEREMY WEST The Kalman filter is arguably one of the most notable algorithms of the 0th century [1]. In this article, we derive
More informationChapter 1. Matrix Algebra
ST4233, Linear Models, Semester 1 2008-2009 Chapter 1. Matrix Algebra 1 Matrix and vector notation Definition 1.1 A matrix is a rectangular or square array of numbers of variables. We use uppercase boldface
More informationAuxiliary signal design for failure detection in uncertain systems
Auxiliary signal design for failure detection in uncertain systems R. Nikoukhah, S. L. Campbell and F. Delebecque Abstract An auxiliary signal is an input signal that enhances the identifiability of a
More informationExample Linear Algebra Competency Test
Example Linear Algebra Competency Test The 4 questions below are a combination of True or False, multiple choice, fill in the blank, and computations involving matrices and vectors. In the latter case,
More informationHANDBOOK OF COMPUTATIONAL ECONOMICS
EIBElEIEIEIEIEIEIEIBBeEieBEieEIEHEIElBE HANDBOOK OF COMPUTATIONAL ECONOMICS VOLUME 1 Edited by HANS M. AMMAN University of Amsterdam DAVID A. KENDRICK University of Texas and JOHN RUST University of Wisconsin
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 information1 Solving Linear Rational Expectations Models
September 29, 2001 SCHUR.TEX Economics 731 International Monetary University of Pennsylvania Theory and Policy Martín Uribe Fall 2001 1 Solving Linear Rational Expectations Models Consider a vector x t
More informationE-Stability vis-a-vis Determinacy Results for a Broad Class. of Linear Rational Expectations Models. Bennett T. McCallum
E-Stability vis-a-vis Determinacy Results for a Broad Class of Linear Rational Expectations Models Bennett T. McCallum Carnegie Mellon University, Tepper School 256, Pittsburgh, PA 15213 USA and National
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 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 informationProperties of Zero-Free Spectral Matrices Brian D. O. Anderson, Life Fellow, IEEE, and Manfred Deistler, Fellow, IEEE
IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL 54, NO 10, OCTOBER 2009 2365 Properties of Zero-Free Spectral Matrices Brian D O Anderson, Life Fellow, IEEE, and Manfred Deistler, Fellow, IEEE Abstract In
More informationTHE LARGE number of components in a power system
IEEE TRANSACTIONS ON POWER SYSTEMS, VOL. 19, NO. 2, MAY 2004 857 Order Reduction of the Dynamic Model of a Linear Weakly Periodic System Part I: General Methodology Abner Ramirez, Member, IEEE, Adam Semlyen,
More information1300 Linear Algebra and Vector Geometry
1300 Linear Algebra and Vector Geometry R. Craigen Office: MH 523 Email: craigenr@umanitoba.ca May-June 2017 Matrix Inversion Algorithm One payoff from this theorem: It gives us a way to invert matrices.
More informationSolution of System of Linear Equations & Eigen Values and Eigen Vectors
Solution of System of Linear Equations & Eigen Values and Eigen Vectors P. Sam Johnson March 30, 2015 P. Sam Johnson (NITK) Solution of System of Linear Equations & Eigen Values and Eigen Vectors March
More informationHands-on Matrix Algebra Using R
Preface vii 1. R Preliminaries 1 1.1 Matrix Defined, Deeper Understanding Using Software.. 1 1.2 Introduction, Why R?.................... 2 1.3 Obtaining R.......................... 4 1.4 Reference Manuals
More informationEL2520 Control Theory and Practice
EL2520 Control Theory and Practice Lecture 8: Linear quadratic control Mikael Johansson School of Electrical Engineering KTH, Stockholm, Sweden Linear quadratic control Allows to compute the controller
More informationLecture 7: Linear-Quadratic Dynamic Programming Real Business Cycle Models
Lecture 7: Linear-Quadratic Dynamic Programming Real Business Cycle Models Shinichi Nishiyama Graduate School of Economics Kyoto University January 10, 2019 Abstract In this lecture, we solve and simulate
More informationResearch Article Minor Prime Factorization for n-d Polynomial Matrices over Arbitrary Coefficient Field
Complexity, Article ID 6235649, 9 pages https://doi.org/10.1155/2018/6235649 Research Article Minor Prime Factorization for n-d Polynomial Matrices over Arbitrary Coefficient Field Jinwang Liu, Dongmei
More informationLMI Based Model Order Reduction Considering the Minimum Phase Characteristic of the System
LMI Based Model Order Reduction Considering the Minimum Phase Characteristic of the System Gholamreza Khademi, Haniyeh Mohammadi, and Maryam Dehghani School of Electrical and Computer Engineering Shiraz
More informationReview of Controllability Results of Dynamical System
IOSR Journal of Mathematics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 13, Issue 4 Ver. II (Jul. Aug. 2017), PP 01-05 www.iosrjournals.org Review of Controllability Results of Dynamical System
More informationAn Alternative Proof of the Greville Formula
JOURNAL OF OPTIMIZATION THEORY AND APPLICATIONS: Vol. 94, No. 1, pp. 23-28, JULY 1997 An Alternative Proof of the Greville Formula F. E. UDWADIA1 AND R. E. KALABA2 Abstract. A simple proof of the Greville
More informationGeneric degree structure of the minimal polynomial nullspace basis: a block Toeplitz matrix approach
Generic degree structure of the minimal polynomial nullspace basis: a block Toeplitz matrix approach Bhaskar Ramasubramanian, Swanand R Khare and Madhu N Belur Abstract This paper formulates the problem
More informationCANONICAL LOSSLESS STATE-SPACE SYSTEMS: STAIRCASE FORMS AND THE SCHUR ALGORITHM
CANONICAL LOSSLESS STATE-SPACE SYSTEMS: STAIRCASE FORMS AND THE SCHUR ALGORITHM Ralf L.M. Peeters Bernard Hanzon Martine Olivi Dept. Mathematics, Universiteit Maastricht, P.O. Box 616, 6200 MD Maastricht,
More informationPARAMETERIZATION OF STATE FEEDBACK GAINS FOR POLE PLACEMENT
PARAMETERIZATION OF STATE FEEDBACK GAINS FOR POLE PLACEMENT Hans Norlander Systems and Control, Department of Information Technology Uppsala University P O Box 337 SE 75105 UPPSALA, Sweden HansNorlander@ituuse
More informationMarkov Chains, Stochastic Processes, and Matrix Decompositions
Markov Chains, Stochastic Processes, and Matrix Decompositions 5 May 2014 Outline 1 Markov Chains Outline 1 Markov Chains 2 Introduction Perron-Frobenius Matrix Decompositions and Markov Chains Spectral
More informationPOLI270 - Linear Algebra
POLI7 - Linear Algebra Septemer 8th Basics a x + a x +... + a n x n b () is the linear form where a, b are parameters and x n are variables. For a given equation such as x +x you only need a variable and
More informationThe dual simplex method with bounds
The dual simplex method with bounds Linear programming basis. Let a linear programming problem be given by min s.t. c T x Ax = b x R n, (P) where we assume A R m n to be full row rank (we will see in the
More informationStatistical and Adaptive Signal Processing
r Statistical and Adaptive Signal Processing Spectral Estimation, Signal Modeling, Adaptive Filtering and Array Processing Dimitris G. Manolakis Massachusetts Institute of Technology Lincoln Laboratory
More informationDiscrete-Time H Gaussian Filter
Proceedings of the 17th World Congress The International Federation of Automatic Control Discrete-Time H Gaussian Filter Ali Tahmasebi and Xiang Chen Department of Electrical and Computer Engineering,
More informationAdaptive State Feedback Nash Strategies for Linear Quadratic Discrete-Time Games
Adaptive State Feedbac Nash Strategies for Linear Quadratic Discrete-Time Games Dan Shen and Jose B. Cruz, Jr. Intelligent Automation Inc., Rocville, MD 2858 USA (email: dshen@i-a-i.com). The Ohio State
More informationPerformance assessment of MIMO systems under partial information
Performance assessment of MIMO systems under partial information H Xia P Majecki A Ordys M Grimble Abstract Minimum variance (MV) can characterize the most fundamental performance limitation of a system,
More informationA Simple Derivation of Right Interactor for Tall Transfer Function Matrices and its Application to Inner-Outer Factorization Continuous-Time Case
A Simple Derivation of Right Interactor for Tall Transfer Function Matrices and its Application to Inner-Outer Factorization Continuous-Time Case ATARU KASE Osaka Institute of Technology Department of
More informationReview Questions REVIEW QUESTIONS 71
REVIEW QUESTIONS 71 MATLAB, is [42]. For a comprehensive treatment of error analysis and perturbation theory for linear systems and many other problems in linear algebra, see [126, 241]. An overview of
More informationDefinition 2.3. We define addition and multiplication of matrices as follows.
14 Chapter 2 Matrices In this chapter, we review matrix algebra from Linear Algebra I, consider row and column operations on matrices, and define the rank of a matrix. Along the way prove that the row
More informationTo appear in IEEE Trans. on Automatic Control Revised 12/31/97. Output Feedback
o appear in IEEE rans. on Automatic Control Revised 12/31/97 he Design of Strictly Positive Real Systems Using Constant Output Feedback C.-H. Huang P.A. Ioannou y J. Maroulas z M.G. Safonov x Abstract
More informationApplications for solving DSGE models. October 25th, 2011
MATLAB Workshop Applications for solving DSGE models Freddy Rojas Cama Marola Castillo Quinto Preliminary October 25th, 2011 A model to solve The model The model A model is set up in order to draw conclusions
More informationUnbiased minimum variance estimation for systems with unknown exogenous inputs
Unbiased minimum variance estimation for systems with unknown exogenous inputs Mohamed Darouach, Michel Zasadzinski To cite this version: Mohamed Darouach, Michel Zasadzinski. Unbiased minimum variance
More informationLinear Algebra. The analysis of many models in the social sciences reduces to the study of systems of equations.
POLI 7 - Mathematical and Statistical Foundations Prof S Saiegh Fall Lecture Notes - Class 4 October 4, Linear Algebra The analysis of many models in the social sciences reduces to the study of systems
More informationComputational Methods for Feedback Control in Damped Gyroscopic Second-order Systems 1
Computational Methods for Feedback Control in Damped Gyroscopic Second-order Systems 1 B. N. Datta, IEEE Fellow 2 D. R. Sarkissian 3 Abstract Two new computationally viable algorithms are proposed for
More informationChapter 1 Theory and Econometrics
Chapter 1 Theory and Econometrics Complete market economies are all alike... Robert E. Lucas, Jr. (1989) 1.1. Introduction Economic theory identifies patterns that unite apparently diverse subjects. Consider
More informationOptimal Polynomial Control for Discrete-Time Systems
1 Optimal Polynomial Control for Discrete-Time Systems Prof Guy Beale Electrical and Computer Engineering Department George Mason University Fairfax, Virginia Correspondence concerning this paper should
More informationLinear-Quadratic Optimal Control: Full-State Feedback
Chapter 4 Linear-Quadratic Optimal Control: Full-State Feedback 1 Linear quadratic optimization is a basic method for designing controllers for linear (and often nonlinear) dynamical systems and is actually
More information1 With state-contingent debt
STOCKHOLM DOCTORAL PROGRAM IN ECONOMICS Helshögskolan i Stockholm Stockholms universitet Paul Klein Email: paul.klein@iies.su.se URL: http://paulklein.se/makro2.html Macroeconomics II Spring 2010 Lecture
More informationMATH 2331 Linear Algebra. Section 2.1 Matrix Operations. Definition: A : m n, B : n p. Example: Compute AB, if possible.
MATH 2331 Linear Algebra Section 2.1 Matrix Operations Definition: A : m n, B : n p ( 1 2 p ) ( 1 2 p ) AB = A b b b = Ab Ab Ab Example: Compute AB, if possible. 1 Row-column rule: i-j-th entry of AB:
More informationElementary Matrices. MATH 322, Linear Algebra I. J. Robert Buchanan. Spring Department of Mathematics
Elementary Matrices MATH 322, Linear Algebra I J. Robert Buchanan Department of Mathematics Spring 2015 Outline Today s discussion will focus on: elementary matrices and their properties, using elementary
More informationOn the adjacency matrix of a block graph
On the adjacency matrix of a block graph R. B. Bapat Stat-Math Unit Indian Statistical Institute, Delhi 7-SJSS Marg, New Delhi 110 016, India. email: rbb@isid.ac.in Souvik Roy Economics and Planning Unit
More informationA note on the equality of the BLUPs for new observations under two linear models
ACTA ET COMMENTATIONES UNIVERSITATIS TARTUENSIS DE MATHEMATICA Volume 14, 2010 A note on the equality of the BLUPs for new observations under two linear models Stephen J Haslett and Simo Puntanen Abstract
More informationLinear Algebra. Min Yan
Linear Algebra Min Yan January 2, 2018 2 Contents 1 Vector Space 7 1.1 Definition................................. 7 1.1.1 Axioms of Vector Space..................... 7 1.1.2 Consequence of Axiom......................
More informationInterlacing Inequalities for Totally Nonnegative Matrices
Interlacing Inequalities for Totally Nonnegative Matrices Chi-Kwong Li and Roy Mathias October 26, 2004 Dedicated to Professor T. Ando on the occasion of his 70th birthday. Abstract Suppose λ 1 λ n 0 are
More informationGeneralized Principal Pivot Transform
Generalized Principal Pivot Transform M. Rajesh Kannan and R. B. Bapat Indian Statistical Institute New Delhi, 110016, India Abstract The generalized principal pivot transform is a generalization of the
More informationSolving Homogeneous Systems with Sub-matrices
Pure Mathematical Sciences, Vol 7, 218, no 1, 11-18 HIKARI Ltd, wwwm-hikaricom https://doiorg/112988/pms218843 Solving Homogeneous Systems with Sub-matrices Massoud Malek Mathematics, California State
More informationLinear Algebra II. 2 Matrices. Notes 2 21st October Matrix algebra
MTH6140 Linear Algebra II Notes 2 21st October 2010 2 Matrices You have certainly seen matrices before; indeed, we met some in the first chapter of the notes Here we revise matrix algebra, consider row
More informationSome New Results on Linear Quadratic Regulator Design for Lossless Systems
Some New Results on Linear Quadratic Regulator Design for Lossless Systems Luigi Fortuna, Giovanni Muscato Maria Gabriella Xibilia Dipartimento Elettrico Elettronico e Sistemistico Universitá degli Studi
More informationNext topics: Solving systems of linear equations
Next topics: Solving systems of linear equations 1 Gaussian elimination (today) 2 Gaussian elimination with partial pivoting (Week 9) 3 The method of LU-decomposition (Week 10) 4 Iterative techniques:
More informationStat 159/259: Linear Algebra Notes
Stat 159/259: Linear Algebra Notes Jarrod Millman November 16, 2015 Abstract These notes assume you ve taken a semester of undergraduate linear algebra. In particular, I assume you are familiar with the
More informationConsistent Bivariate Distribution
A Characterization of the Normal Conditional Distributions MATSUNO 79 Therefore, the function ( ) = G( : a/(1 b2)) = N(0, a/(1 b2)) is a solu- tion for the integral equation (10). The constant times of
More informationContents. 1 State-Space Linear Systems 5. 2 Linearization Causality, Time Invariance, and Linearity 31
Contents Preamble xiii Linear Systems I Basic Concepts 1 I System Representation 3 1 State-Space Linear Systems 5 1.1 State-Space Linear Systems 5 1.2 Block Diagrams 7 1.3 Exercises 11 2 Linearization
More informationSOLVING LINEAR RATIONAL EXPECTATIONS MODELS. Three ways to solve a linear model
SOLVING LINEAR RATIONAL EXPECTATIONS MODELS KRISTOFFER P. NIMARK Three ways to solve a linear model Solving a model using full information rational expectations as the equilibrium concept involves integrating
More informationOn the Principle of Optimality for Nonstationary Deterministic Dynamic Programming
On the Principle of Optimality for Nonstationary Deterministic Dynamic Programming Takashi Kamihigashi January 15, 2007 Abstract This note studies a general nonstationary infinite-horizon optimization
More informationDepartment of Econometrics and Business Statistics
ISSN 440-77X Australia Department of Econometrics and Business Statistics http://wwwbusecomonasheduau/depts/ebs/pubs/wpapers/ The Asymptotic Distribution of the LIML Estimator in a artially Identified
More informationDenis ARZELIER arzelier
COURSE ON LMI OPTIMIZATION WITH APPLICATIONS IN CONTROL PART II.2 LMIs IN SYSTEMS CONTROL STATE-SPACE METHODS PERFORMANCE ANALYSIS and SYNTHESIS Denis ARZELIER www.laas.fr/ arzelier arzelier@laas.fr 15
More informationTitle without the persistently exciting c. works must be obtained from the IEE
Title Exact convergence analysis of adapt without the persistently exciting c Author(s) Sakai, H; Yang, JM; Oka, T Citation IEEE TRANSACTIONS ON SIGNAL 55(5): 2077-2083 PROCESS Issue Date 2007-05 URL http://hdl.handle.net/2433/50544
More informationStochastic convexity in dynamic programming
Economic Theory 22, 447 455 (2003) Stochastic convexity in dynamic programming Alp E. Atakan Department of Economics, Columbia University New York, NY 10027, USA (e-mail: aea15@columbia.edu) Received:
More informationA method for computing quadratic Brunovsky forms
Electronic Journal of Linear Algebra Volume 13 Volume 13 (25) Article 3 25 A method for computing quadratic Brunovsky forms Wen-Long Jin wjin@uciedu Follow this and additional works at: http://repositoryuwyoedu/ela
More informationGENERALIZED POWER SYMMETRIC STOCHASTIC MATRICES
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 27, Number 7, Pages 987 994 S 0002-9939(99)0526-6 Article electronically published on March 7, 999 GENERALIZED POWER SYMMETRIC STOCHASTIC MATRICES
More informationOptimal Control. McGill COMP 765 Oct 3 rd, 2017
Optimal Control McGill COMP 765 Oct 3 rd, 2017 Classical Control Quiz Question 1: Can a PID controller be used to balance an inverted pendulum: A) That starts upright? B) That must be swung-up (perhaps
More informationSTABILITY OF SECOND-ORDER ASYMMETRIC LINEAR MECHANICAL SYSTEMS WITH APPLICATION TO ROBOT GRASPING
STABILITY OF SECOND-ORDER ASYMMETRIC LINEAR MECHANICAL SYSTEMS WITH APPLICATION TO ROBOT GRASPING Amir Shapiro Department of Mechanical Engineering Ben Gurion University of the Negev P.O.B. 653 Beer Sheva
More information