Abstract. + n) variables with a rank of ( n2 (n+1) )+n 2 ] equations in ( n2 (n+1)
|
|
- Philomena Walters
- 5 years ago
- Views:
Transcription
1 SOLUTION OF SECOND ORDER LINEARIZATION R.Devanathan School of Electrical and Electronic Engineering (Block S) Nanyang Technological University Singapore 6998 Abstract For a nonlinear system with a control input, a generalized form of the homological equation can be formulated to reduce the system to its normal form. In this paper, it has been shown that for a linearly controllable system with an appropriate choice of state feedback, the generalized homological equation can be solved to give an explicit solution, of a reduced order, to the problem of second order linearization. Introduction Linearization of nonlinear dynamic system was originally investigated by Poincare[, ] in terms of homological equations. Krener et. al.[] have considered a nonlinear system with a control input and showed that a generalized form of the homological equation can be formulated in this case.devanathan[]has shown that, for a linearly controllable system with an appropriate choice of state feedback, the system matrix can be made non-resonant.this concept is further exploited in this paper in the case of second order linearization. First, the normalizing transformation is explicitly solved for in an unique way. This is followed by the solution of the parameters of the input transformation using a system of equations which can be reduced to ( n(n ) ) equations in n variables whose rank is (n ).This is in contrast to the result of [] which derives a system of [( n (n+) )+n ] equations in ( n (n+) )+ n(n+) + n) variables with a rank of ( n (n+) )+ n(n+) + n ). Moreover, the solution derived in this paper is in an explicit form while that of [] is not. Background Krener et.al. [] have extended the Poincare s result on the normalizing transformation to include the control input in a straightforward way. Consider the equation ẋ = Ax + Gu + f (x)+f (x)+ + g (x)u + g (x)u + (.) the eigenvalues of A are non-resonant, u corresponds to a scalar input and f m (x) (g m (x)) corresponds to a vector-valued polynomial containing terms of the form x q x q x qn n,q i,i =,,,n are integers and n i= q i = m( n i= q i = m ),m. Using a change of variable as in x = y + h(y) (.)
2 consider a change of input as in v =(+β (x))u + α (x), +β (x) (.) β (x) andα (x) correspond to terms linear and of second order in x respectively. Following Devanathan [], it can be shown that the second order term f (x) and the g (x)u term can be removed provided the following generalized homological equations are satisfied. h (y) (Ay) Ah (y)+gα (y) =f (y), (.) y h (y) (Gu)+Gβ (y)u = g (y)u, u, (.) y State Feedback for Non-resonance In this paper, we consider a class of systems of the form ξ = f(ξ)+g(ξ)ζ (.) ξ is an n-tuple vector and f(ξ)and g(ξ) are vector fields.for simplicity, we assume a scalar input ζ. We require that the system (.) be linearly controllable [],i.e., the pair (F, G) is controllable F = f () and G = g(), at the assumed equilibrium point at ξ the origin. One can consider, without loss of generality, F and G to be in Brunovsky form [6]. The power series expansion of (.) about the origin can then be written as ẋ = Fx+ Gφ + O (x) () + γ (x, φ) () (.) superscript () corresponds to terms in x of degree greater than one, superscript () corresponds to terms in x of degree greater than or equal to one and x and φ are the transformed state and input variables respectively. To put (.) in the form of (.), matrix A is non-resonant,we introduce state feedback as in (.) then becomes φ = Kx + u (.) K =[k n,k n,,k,k ] t (.) ẋ = Ax + Gu+ O(x) () + γ(x, u) () (.) A = F GK (.6) Remark. We can choose the eigenvalues of (.6), without loss of generality, to be real and distinct ([]). Since, matrix A in (.6) is in phase-variable canonical form, using a change of coordinate involving the Vandermonde matrix [], (.) can be put in the form of (.) matrix A is diagonal and G =[,,, ] t. The solution of the generalized homological equations (.) and (.), in this case, results in the second order linearization.
3 Equations for Second Order Linearization Let h (y) =[h (y),h (y),,h n (y)] t (.) h i (y) = Q h iq Y Q,i=,,,n (.) Q =[q,q,,q n ] (.) n q =, q, =,,,n (.) = Y Q = y q y q yn qn (.) α (y) = Q α Q Y Q (.6) f i (y) = Q f iq Y Q,i=,,,n (.) f =[f,f,,f n ] t (.8) Further, given Q as in (.),let Q =[q,q,,q, (q ),q +,,q n ], q =, (,,,n) (.9) and for, k (,,,n). [q,q,,q k, (q k +),q k+,,q, (q ),q +,,q n ], > k Q k = [q,q,,q, (q ),q +,,q k, (q k +),q k+,,q n ], k > Q, k =, (.) q =,q k =,,k (,,,n) Let β (y) = Q b Q Y Q (.) g =[g,g,,,g n ] t (.) g i = g iq Y Q,i=,,,n Q (.) Then, combining Equations (.) and (.), one can get for q =,q k =, =,,,n; k (,,,n); i =,,,n q (<Q,Λ > λ i ) (α Q )+ [(q k +)(<Q k, Λ > λ i ) (α Q k )] b Q k = d i (.)
4 d i is a known quantity given by d i = q (<Q,Λ > λ i ) (f iq )+ [(q k +)(<Q k, Λ > λ i) (f iq k )] g iq k (.) The unknown parameters α Q,α Q k and b Q in the input transformation are to be solved from the equation (.). Using the simplified notations as follows: α Q k = α k ; α Q = α and b Q = b ; q =,q k =,, k (,,,n);i =,,,n (.6) Equation (.) can be put in the form, for =,,,n, d = C α Lb (.) d =[d,d,,d n ] t (.8) α =[α,α,,α n ] t (.9) L =[,,,, ] t (.) (.) and C is an n n matrix whose (i, k)-th element is given by C (i, k) = λ +λ k λ i, ( λ λ i ), k = k (.) i, k =,,,n; =,,,n The following two results are stated without proof. Theorem. For each =,,,n,then n matrix C is nonsingular. Representing C in terms of its columns as C =[c,c,,c NJ ] (.) and noting that c k = c k,α k = α k,k, (,,,n) Equation (.) can now be put in the form d is a n -tuple vector given as d =[DB] [ α b ] (.) d =[d t,dt,,dt n ]t (.)
5 d, (,,,n)isann-tuple vector given as in (.8). α is a n(n+) -tuple vector whose elements are α k, =,,,n,k =, +,,n.bisn-tuple vector given as b =[b,b,,b n ] t (.6) Disa(n x n(n+) ) matrix which can be explicitly expressed in terms of the elements c k,,k (,,,n). B is a n n matrix put in block diagonal form as B = dig[ L, L,, L] (.) Theorem. The rank of the n ( n(n+) + n) matrix [D B] is ( n(n+) + n ). Moreover, equation (.) further reduces to a system of n(n ) linear equations in b,=,,,n with rank (n ) whose solution corresponds to the solution of second order linearization. Example We illustrate the complete procedure through an example. Consider the system ẋ = x + u +.x x x.x + x x x u + O(x)() + O (x) () u (.) x =[x,x,x ] t, (.) and superscript ()(superscript ()) corresponds to terms in x of degree (). The state feedback u = v Kx applied to the above system K =[, 9, 6 6 ]t results in the eigenvalues,, leading to the system ż = z + v + 8z + 8 z + z 96z z 8z z + 69z z z + z + 6 z z z 69 z z + 8z z z + 68 z + 8z 6z z 86z z + 8z z z +6z +z z +z +z z +6z +z v + O (z) () + O (z) () v (.)
6 Formulating the equation(.), =,, δ = (.) δ =[α,α,α,α,α,α,b,b,b ] t (.) The matrix on the OHS of equation (.) has rank 8 only.the parameters α k,,k (,, ) can be uniquely expressed in terms of b, =,,.Hence, for reasons of consistency of the values of b, =,,, we have, b b b = 6 (.6) In the above, one of the b, (,, ) is arbitrary. Putting b =,b =8andb =, the normalizing and the input transformations are computed respectively as z = h (y)+y (.) h (y) = y + y y y y y + y y 9y 6y + 6y y y y 6y + 8y y y y + 8y y 6y y 89y + 8y y 89y (.8) w =(+8z z )v +z + z z z z 96 z + 6z z z (.9) The original system (.) then reduces to ẏ = dig(,, )y + []t w + O () (y)+o () (y, w) (.) the second degree terms in y and the product terms of the form (linear term in y)w are completely removed. 6
7 6 Conclusions The problem of transforming a dynamic system into its normal form is an old one. In the reduction to the normal form, the role played by the resonance of the system matrix is well-known. In the context of linearly controllable system, however, the control input can be used to provide state feedback so as to ensure the non-resonance of the system matrix. This, in turn, can lead to the explicit solution of second order linearization as shown in this paper.the future work is to extend the linearization technique (for a linearly controllable system) to an arbitrary order. Acknowledgment: The author would like to thank Nanyang Technological University, Singapore for the financial support towards presentation of this paper. References [] V.I.Arnold, Geometrical Methods in the Theory of Ordinary Differential Equations, Springer-Verlag, New York,pp.-88, 98. [] J. Guckenheimer and P.Holmes, Nonlinear Oscillations, Dynamical Systems and Bifurcation of Vector Fields, Springer-Verlag, New York,98. [] A.J.Krener, S.Karahan, and M.Hubbard, Approximate Normal Forms of Nonlinear Systems, Proc. th Conf. Decision and Control, pp.-9, 988. [] R. Devanathan, Linearization Condition Through State Feedback, IEEE Transactions on Automatic Control, Vol.6, No.8, Aug., pp.-6. [] A. J. Krener and W.Kang, Extended Normal Forms of Quadratic Systems, Proc. 9th Conf. Decision and Control, pp.9-9,99. [6] Brunovsky, P., A Classification of Linear Controllable Systems, Kybernetica cislo, Vol.,pp.-88, 9. [] B.C.Kuo,Automatic Control Systems Fourth Edition,Prentice-Hall,98,pp.-.
On the convergence of normal forms for analytic control systems
Chapter 1 On the convergence of normal forms for analytic control systems Wei Kang Department of Mathematics, Naval Postgraduate School Monterey, CA 93943 wkang@nps.navy.mil Arthur J. Krener Department
More informationChap. 3. Controlled Systems, Controllability
Chap. 3. Controlled Systems, Controllability 1. Controllability of Linear Systems 1.1. Kalman s Criterion Consider the linear system ẋ = Ax + Bu where x R n : state vector and u R m : input vector. A :
More informationLIE TRANSFORM CONSTRUCTION OF EXPONENTIAL NORMAL FORM OBSERVERS. Suba Thomas Harry G. Kwatny Bor-Chin Chang
LIE TRANSFORM CONSTRUCTION OF EXPONENTIAL NORMAL FORM OBSERVERS Suba Thomas Harry G Kwatny Bor-Chin Chang Department of Mechanical Engineering and Mechanics, Drexel University, Philadelphia, PA 1914, USA
More informationMath 1060 Linear Algebra Homework Exercises 1 1. Find the complete solutions (if any!) to each of the following systems of simultaneous equations:
Homework Exercises 1 1 Find the complete solutions (if any!) to each of the following systems of simultaneous equations: (i) x 4y + 3z = 2 3x 11y + 13z = 3 2x 9y + 2z = 7 x 2y + 6z = 2 (ii) x 4y + 3z =
More informationA matching pursuit technique for computing the simplest normal forms of vector fields
Journal of Symbolic Computation 35 (2003) 59 65 www.elsevier.com/locate/jsc A matching pursuit technique for computing the simplest normal forms of vector fields Pei Yu, Yuan Yuan Department of Applied
More informationCSL361 Problem set 4: Basic linear algebra
CSL361 Problem set 4: Basic linear algebra February 21, 2017 [Note:] If the numerical matrix computations turn out to be tedious, you may use the function rref in Matlab. 1 Row-reduced echelon matrices
More informationNATIONAL UNIVERSITY OF SINGAPORE PC5215 NUMERICAL RECIPES WITH APPLICATIONS. (Semester I: AY ) Time Allowed: 2 Hours
NATIONAL UNIVERSITY OF SINGAPORE PC55 NUMERICAL RECIPES WITH APPLICATIONS (Semester I: AY 0-) Time Allowed: Hours INSTRUCTIONS TO CANDIDATES. This examination paper contains FIVE questions and comprises
More informationIMU-Laser Scanner Localization: Observability Analysis
IMU-Laser Scanner Localization: Observability Analysis Faraz M. Mirzaei and Stergios I. Roumeliotis {faraz stergios}@cs.umn.edu Dept. of Computer Science & Engineering University of Minnesota Minneapolis,
More informationNon-regular feedback linearization of nonlinear systems via a normal form algorithm
Available online at www.sciencedirect.com Automatica 4 24 439 447 www.elsevier.com/locate/automatica Brief paper Non-regular feedback linearization of nonlinear systems via a normal form algorithm Daizhan
More informationControl design using Jordan controllable canonical form
Control design using Jordan controllable canonical form Krishna K Busawon School of Engineering, Ellison Building, University of Northumbria at Newcastle, Newcastle upon Tyne NE1 8ST, UK email: krishnabusawon@unnacuk
More informationStability of Hybrid Control Systems Based on Time-State Control Forms
Stability of Hybrid Control Systems Based on Time-State Control Forms Yoshikatsu HOSHI, Mitsuji SAMPEI, Shigeki NAKAURA Department of Mechanical and Control Engineering Tokyo Institute of Technology 2
More informationHIGHER ORDER CRITERION FOR THE NONEXISTENCE OF FORMAL FIRST INTEGRAL FOR NONLINEAR SYSTEMS. 1. Introduction
Electronic Journal of Differential Equations, Vol. 217 (217), No. 274, pp. 1 11. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu HIGHER ORDER CRITERION FOR THE NONEXISTENCE
More informationAnalysis of the Takens-Bogdanov bifurcation on m parameterized vector fields
Analysis of the Takens-Bogdanov bifurcation on m parameterized vector fields Francisco Armando Carrillo Navarro, Fernando Verduzco G., Joaquín Delgado F. Programa de Doctorado en Ciencias (Matemáticas),
More informationSolving x 2 Dy 2 = N in integers, where D > 0 is not a perfect square. Keith Matthews
Solving x 2 Dy 2 = N in integers, where D > 0 is not a perfect square. Keith Matthews Abstract We describe a neglected algorithm, based on simple continued fractions, due to Lagrange, for deciding the
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 informationFirst Integrals/Invariants & Symmetries for Autonomous Difference Equations
Proceedings of Institute of Mathematics of NAS of Ukraine 2004, Vol. 50, Part 3, 1253 1260 First Integrals/Invariants & Symmetries for Autonomous Difference Equations Leon ARRIOLA Department of Mathematical
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 informationAN EFFICIENT METHOD FOR COMPUTING THE SIMPLEST NORMAL FORMS OF VECTOR FIELDS
International Journal of Bifurcation and Chaos, Vol. 3, No. (2003) 9 46 c World Scientific Publishing Company AN EFFICIENT METHOD FOR COMPUTING THE SIMPLEST NORMAL FORMS OF VECTOR FIELDS P. YU and Y. YUAN
More informationSolutions to Dynamical Systems 2010 exam. Each question is worth 25 marks.
Solutions to Dynamical Systems exam Each question is worth marks [Unseen] Consider the following st order differential equation: dy dt Xy yy 4 a Find and classify all the fixed points of Hence draw the
More informationBUMPLESS SWITCHING CONTROLLERS. William A. Wolovich and Alan B. Arehart 1. December 27, Abstract
BUMPLESS SWITCHING CONTROLLERS William A. Wolovich and Alan B. Arehart 1 December 7, 1995 Abstract This paper outlines the design of bumpless switching controllers that can be used to stabilize MIMO plants
More informationLinear differential equations with constant coefficients Method of undetermined coefficients
Linear differential equations with constant coefficients Method of undetermined coefficients e u+vi = e u (cos vx + i sin vx), u, v R, i 2 = -1 Quasi-polynomial: Q α+βi,k (x) = e αx [cos βx ( f 0 + f 1
More informationIntroduction to Geometric Control
Introduction to Geometric Control Relative Degree Consider the square (no of inputs = no of outputs = p) affine control system ẋ = f(x) + g(x)u = f(x) + [ g (x),, g p (x) ] u () h (x) y = h(x) = (2) h
More informationENERGY DECAY ESTIMATES FOR LIENARD S EQUATION WITH QUADRATIC VISCOUS FEEDBACK
Electronic Journal of Differential Equations, Vol. 00(00, No. 70, pp. 1 1. ISSN: 107-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp ENERGY DECAY ESTIMATES
More informationEntrance Exam, Differential Equations April, (Solve exactly 6 out of the 8 problems) y + 2y + y cos(x 2 y) = 0, y(0) = 2, y (0) = 4.
Entrance Exam, Differential Equations April, 7 (Solve exactly 6 out of the 8 problems). Consider the following initial value problem: { y + y + y cos(x y) =, y() = y. Find all the values y such that the
More informationIEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 49, NO. 8, AUGUST Control Bifurcations
IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 49, NO. 8, AUGUST 2004 1231 Control Bifurcations Arthur J. Krener, Fellow, IEEE, Wei Kang, Senior Member, IEEE, Dong Eui Chang Abstract A parametrized nonlinear
More informationHopf Bifurcation Control for Affine Systems
opf Bifurcation Control for Affine Systems Ferno Verduzco Joaquin Alvarez Abstract In this paper we establish conditions to control the opf bifurcation of nonlinear systems with two uncontrollable modes
More informationNonlinear Control Systems
Nonlinear Control Systems António Pedro Aguiar pedro@isr.ist.utl.pt 7. Feedback Linearization IST-DEEC PhD Course http://users.isr.ist.utl.pt/%7epedro/ncs1/ 1 1 Feedback Linearization Given a nonlinear
More informationCENTER MANIFOLD AND NORMAL FORM THEORIES
3 rd Sperlonga Summer School on Mechanics and Engineering Sciences 3-7 September 013 SPERLONGA CENTER MANIFOLD AND NORMAL FORM THEORIES ANGELO LUONGO 1 THE CENTER MANIFOLD METHOD Existence of an invariant
More information235 Final exam review questions
5 Final exam review questions Paul Hacking December 4, 0 () Let A be an n n matrix and T : R n R n, T (x) = Ax the linear transformation with matrix A. What does it mean to say that a vector v R n is an
More informationBisimilar Finite Abstractions of Interconnected Systems
Bisimilar Finite Abstractions of Interconnected Systems Yuichi Tazaki and Jun-ichi Imura Tokyo Institute of Technology, Ōokayama 2-12-1, Meguro, Tokyo, Japan {tazaki,imura}@cyb.mei.titech.ac.jp http://www.cyb.mei.titech.ac.jp
More informationExercise Sheet 1.
Exercise Sheet 1 You can download my lecture and exercise sheets at the address http://sami.hust.edu.vn/giang-vien/?name=huynt 1) Let A, B be sets. What does the statement "A is not a subset of B " mean?
More informationLinear Algebra: Matrix Eigenvalue Problems
CHAPTER8 Linear Algebra: Matrix Eigenvalue Problems Chapter 8 p1 A matrix eigenvalue problem considers the vector equation (1) Ax = λx. 8.0 Linear Algebra: Matrix Eigenvalue Problems Here A is a given
More informationSymmetry Methods for Differential and Difference Equations. Peter Hydon
Lecture 2: How to find Lie symmetries Symmetry Methods for Differential and Difference Equations Peter Hydon University of Kent Outline 1 Reduction of order for ODEs and O Es 2 The infinitesimal generator
More informationCHAPTER 3. Matrix Eigenvalue Problems
A SERIES OF CLASS NOTES FOR 2005-2006 TO INTRODUCE LINEAR AND NONLINEAR PROBLEMS TO ENGINEERS, SCIENTISTS, AND APPLIED MATHEMATICIANS DE CLASS NOTES 3 A COLLECTION OF HANDOUTS ON SYSTEMS OF ORDINARY DIFFERENTIAL
More informationYURI LEVIN, MIKHAIL NEDIAK, AND ADI BEN-ISRAEL
Journal of Comput. & Applied Mathematics 139(2001), 197 213 DIRECT APPROACH TO CALCULUS OF VARIATIONS VIA NEWTON-RAPHSON METHOD YURI LEVIN, MIKHAIL NEDIAK, AND ADI BEN-ISRAEL Abstract. Consider m functions
More informationarxiv: v1 [math.ds] 17 Nov 2015
arxiv:1511.05836v1 [math.ds] 17 Nov 2015 A Note On Topological Conjugacy For Perpetual Points Awadhesh Prasad 1 Department of Physics and Astrophysics, University of Delhi, Delhi 110007, India. Abstract:
More information2.10 Saddles, Nodes, Foci and Centers
2.10 Saddles, Nodes, Foci and Centers In Section 1.5, a linear system (1 where x R 2 was said to have a saddle, node, focus or center at the origin if its phase portrait was linearly equivalent to one
More information1 Determinants. 1.1 Determinant
1 Determinants [SB], Chapter 9, p.188-196. [SB], Chapter 26, p.719-739. Bellow w ll study the central question: which additional conditions must satisfy a quadratic matrix A to be invertible, that is to
More information1. Introduction. The problem of transforming the nonlinear control singleinput. ξ = f(ξ)+g(ξ)u. x = φ(ξ), u = α(ξ)+β(ξ)v
SIAM J CONTROL OPTIM Vol 0, No 0, pp 000 000 c 2003 Society for Industrial and Applied Mathematics FEEDBACK CLASSIFICATION OF NONLINEAR SINGLE-INPUT CONTROL SYSTEMS WITH CONTROLLABLE LINEARIZATION: NORMAL
More informationJordan normal form notes (version date: 11/21/07)
Jordan normal form notes (version date: /2/7) If A has an eigenbasis {u,, u n }, ie a basis made up of eigenvectors, so that Au j = λ j u j, then A is diagonal with respect to that basis To see this, let
More informationObservability. Dynamic Systems. Lecture 2 Observability. Observability, continuous time: Observability, discrete time: = h (2) (x, u, u)
Observability Dynamic Systems Lecture 2 Observability Continuous time model: Discrete time model: ẋ(t) = f (x(t), u(t)), y(t) = h(x(t), u(t)) x(t + 1) = f (x(t), u(t)), y(t) = h(x(t)) Reglerteknik, ISY,
More informationEconomics 204 Fall 2013 Problem Set 5 Suggested Solutions
Economics 204 Fall 2013 Problem Set 5 Suggested Solutions 1. Let A and B be n n matrices such that A 2 = A and B 2 = B. Suppose that A and B have the same rank. Prove that A and B are similar. Solution.
More informationOn Bifurcations in Nonlinear Consensus Networks
American Control Conference Marriott Waterfront, Baltimore, MD, USA June -July, WeC.4 On Bifurcations in Nonlinear Consensus Networks Vaibhav Srivastava Jeff Moehlis Francesco Bullo Abstract The theory
More informationDegree Master of Science in Mathematical Modelling and Scientific Computing Mathematical Methods I Thursday, 12th January 2012, 9:30 a.m.- 11:30 a.m.
Degree Master of Science in Mathematical Modelling and Scientific Computing Mathematical Methods I Thursday, 12th January 2012, 9:30 a.m.- 11:30 a.m. Candidates should submit answers to a maximum of four
More informationHopf Bifurcation and Limit Cycle Analysis of the Rikitake System
ISSN 749-3889 (print), 749-3897 (online) International Journal of Nonlinear Science Vol.4(0) No.,pp.-5 Hopf Bifurcation and Limit Cycle Analysis of the Rikitake System Xuedi Wang, Tianyu Yang, Wei Xu Nonlinear
More informationLinear Algebra Practice Problems
Linear Algebra Practice Problems Math 24 Calculus III Summer 25, Session II. Determine whether the given set is a vector space. If not, give at least one axiom that is not satisfied. Unless otherwise stated,
More informationSTABILITY. Phase portraits and local stability
MAS271 Methods for differential equations Dr. R. Jain STABILITY Phase portraits and local stability We are interested in system of ordinary differential equations of the form ẋ = f(x, y), ẏ = g(x, y),
More informationNPTEL Online Course: Control Engineering
NPTEL Online Course: Control Engineering Ramkrishna Pasumarthy Assignment-11 : s 1. Consider a system described by state space model [ ] [ 0 1 1 x + u 5 1 2] y = [ 1 2 ] x What is the transfer function
More informationKernel Methods. Charles Elkan October 17, 2007
Kernel Methods Charles Elkan elkan@cs.ucsd.edu October 17, 2007 Remember the xor example of a classification problem that is not linearly separable. If we map every example into a new representation, then
More information(0, 0), (1, ), (2, ), (3, ), (4, ), (5, ), (6, ).
1 Interpolation: The method of constructing new data points within the range of a finite set of known data points That is if (x i, y i ), i = 1, N are known, with y i the dependent variable and x i [x
More informationAPPPHYS217 Tuesday 25 May 2010
APPPHYS7 Tuesday 5 May Our aim today is to take a brief tour of some topics in nonlinear dynamics. Some good references include: [Perko] Lawrence Perko Differential Equations and Dynamical Systems (Springer-Verlag
More informationMinimum-Phase Property of Nonlinear Systems in Terms of a Dissipation Inequality
Minimum-Phase Property of Nonlinear Systems in Terms of a Dissipation Inequality Christian Ebenbauer Institute for Systems Theory in Engineering, University of Stuttgart, 70550 Stuttgart, Germany ce@ist.uni-stuttgart.de
More informationAppendix A Solving Linear Matrix Inequality (LMI) Problems
Appendix A Solving Linear Matrix Inequality (LMI) Problems In this section, we present a brief introduction about linear matrix inequalities which have been used extensively to solve the FDI problems described
More informationCONTROL DESIGN FOR SET POINT TRACKING
Chapter 5 CONTROL DESIGN FOR SET POINT TRACKING In this chapter, we extend the pole placement, observer-based output feedback design to solve tracking problems. By tracking we mean that the output is commanded
More informationMA22S3 Summary Sheet: Ordinary Differential Equations
MA22S3 Summary Sheet: Ordinary Differential Equations December 14, 2017 Kreyszig s textbook is a suitable guide for this part of the module. Contents 1 Terminology 1 2 First order separable 2 2.1 Separable
More informationLOWELL. MICHIGAN, THURSDAY, MAY 23, Schools Close. method of distribution. t o t h e b o y s of '98 a n d '18. C o m e out a n d see t h e m get
UY Y 99 U XXX 5 Q == 5 K 8 9 $8 7 Y Y 8 q 6 Y x) 5 7 5 q U «U YU x q Y U ) U Z 9 Y 7 x 7 U x U x 9 5 & U Y U U 7 8 9 x 5 x U x x ; x x [ U K»5 98 8 x q 8 q K x Y Y x K Y 5 ~ 8» Y x ; ; ; ; ; ; 8; x 7;
More informationOn the positivity of linear weights in WENO approximations. Abstract
On the positivity of linear weights in WENO approximations Yuanyuan Liu, Chi-Wang Shu and Mengping Zhang 3 Abstract High order accurate weighted essentially non-oscillatory (WENO) schemes have been used
More informationMathematical Induction Assignments
1 Mathematical Induction Assignments Prove the Following using Principle of Mathematical induction 1) Prove that for any positive integer number n, n 3 + 2 n is divisible by 3 2) Prove that 1 3 + 2 3 +
More informationOutput Feedback and State Feedback. EL2620 Nonlinear Control. Nonlinear Observers. Nonlinear Controllers. ẋ = f(x,u), y = h(x)
Output Feedback and State Feedback EL2620 Nonlinear Control Lecture 10 Exact feedback linearization Input-output linearization Lyapunov-based control design methods ẋ = f(x,u) y = h(x) Output feedback:
More information1 Assignment 1: Nonlinear dynamics (due September
Assignment : Nonlinear dynamics (due September 4, 28). Consider the ordinary differential equation du/dt = cos(u). Sketch the equilibria and indicate by arrows the increase or decrease of the solutions.
More informationMATH 3030, Abstract Algebra Winter 2012 Toby Kenney Sample Midterm Examination Model Solutions
MATH 3030, Abstract Algebra Winter 2012 Toby Kenney Sample Midterm Examination Model Solutions Basic Questions 1. Give an example of a prime ideal which is not maximal. In the ring Z Z, the ideal {(0,
More informationOn Bifurcations. in Nonlinear Consensus Networks
On Bifurcations in Nonlinear Consensus Networks Vaibhav Srivastava Jeff Moehlis Francesco Bullo Abstract The theory of consensus dynamics is widely employed to study various linear behaviors in networked
More informationA Sufficient Condition for Local Controllability
A Sufficient Condition for Local Controllability of Nonlinear Systems Along Closed Orbits Kwanghee Nam and Aristotle Arapostathis Abstract We present a computable sufficient condition to determine local
More informationPOWER systems are increasingly operated closer to their
1438 IEEE TRANSACTIONS ON POWER SYSTEMS, VOL. 19, NO. 3, AUGUST 2004 Scaling of Normal Form Analysis Coefficients Under Coordinate Change Ian Dobson, Senior Member, IEEE, and Emilio Barocio, Member, IEEE
More informationRadial eigenvectors of the Laplacian of the nonbinary hypercube
Radial eigenvectors of the Laplacian of the nonbinary hypercube Murali K. Srinivasan Department of Mathematics Indian Institute of Technology, Bombay Powai, Mumbai 400076, INDIA mks@math.iitb.ac.in murali.k.srinivasan@gmail.com
More informationBASIC MATRIX ALGEBRA WITH ALGORITHMS AND APPLICATIONS ROBERT A. LIEBLER CHAPMAN & HALL/CRC
BASIC MATRIX ALGEBRA WITH ALGORITHMS AND APPLICATIONS ROBERT A. LIEBLER CHAPMAN & HALL/CRC A CRC Press Company Boca Raton London New York Washington, D.C. Contents Preface Examples Major results/proofs
More informationCS145: INTRODUCTION TO DATA MINING
CS145: INTRODUCTION TO DATA MINING 5: Vector Data: Support Vector Machine Instructor: Yizhou Sun yzsun@cs.ucla.edu October 18, 2017 Homework 1 Announcements Due end of the day of this Thursday (11:59pm)
More informationCONVERGENCE PROOF FOR RECURSIVE SOLUTION OF LINEAR-QUADRATIC NASH GAMES FOR QUASI-SINGULARLY PERTURBED SYSTEMS. S. Koskie, D. Skataric and B.
To appear in Dynamics of Continuous, Discrete and Impulsive Systems http:monotone.uwaterloo.ca/ journal CONVERGENCE PROOF FOR RECURSIVE SOLUTION OF LINEAR-QUADRATIC NASH GAMES FOR QUASI-SINGULARLY PERTURBED
More informationA matrix is a rectangular array of. objects arranged in rows and columns. The objects are called the entries. is called the size of the matrix, and
Section 5.5. Matrices and Vectors A matrix is a rectangular array of objects arranged in rows and columns. The objects are called the entries. A matrix with m rows and n columns is called an m n matrix.
More informationINVERSION IN INDIRECT OPTIMAL CONTROL
INVERSION IN INDIRECT OPTIMAL CONTROL François Chaplais, Nicolas Petit Centre Automatique et Systèmes, École Nationale Supérieure des Mines de Paris, 35, rue Saint-Honoré 7735 Fontainebleau Cedex, France,
More informationChapter 2 Hopf Bifurcation and Normal Form Computation
Chapter 2 Hopf Bifurcation and Normal Form Computation In this chapter, we discuss the computation of normal forms. First we present a general approach which combines center manifold theory with computation
More informationEIGENVALUE PROBLEMS. EIGENVALUE PROBLEMS p. 1/4
EIGENVALUE PROBLEMS EIGENVALUE PROBLEMS p. 1/4 EIGENVALUE PROBLEMS p. 2/4 Eigenvalues and eigenvectors Let A C n n. Suppose Ax = λx, x 0, then x is a (right) eigenvector of A, corresponding to the eigenvalue
More informationNonlinear Discrete-Time Observer Design with Linearizable Error Dynamics
622 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 48, NO. 4, APRIL 2003 Nonlinear Discrete-Time Observer Design with Linearizable Error Dynamics MingQing Xiao, Nikolaos Kazantzis, Costas Kravaris, Arthur
More informationHomework 6 Solutions. Solution. Note {e t, te t, t 2 e t, e 2t } is linearly independent. If β = {e t, te t, t 2 e t, e 2t }, then
Homework 6 Solutions 1 Let V be the real vector space spanned by the functions e t, te t, t 2 e t, e 2t Find a Jordan canonical basis and a Jordan canonical form of T on V dened by T (f) = f Solution Note
More informationADAPTIVE FEEDBACK LINEARIZING CONTROL OF CHUA S CIRCUIT
International Journal of Bifurcation and Chaos, Vol. 12, No. 7 (2002) 1599 1604 c World Scientific Publishing Company ADAPTIVE FEEDBACK LINEARIZING CONTROL OF CHUA S CIRCUIT KEVIN BARONE and SAHJENDRA
More informationZERO-HOPF BIFURCATION FOR A CLASS OF LORENZ-TYPE SYSTEMS
This is a preprint of: Zero-Hopf bifurcation for a class of Lorenz-type systems, Jaume Llibre, Ernesto Pérez-Chavela, Discrete Contin. Dyn. Syst. Ser. B, vol. 19(6), 1731 1736, 214. DOI: [doi:1.3934/dcdsb.214.19.1731]
More informationCDS 101/110a: Lecture 2.1 Dynamic Behavior
CDS 11/11a: Lecture.1 Dynamic Behavior Richard M. Murray 6 October 8 Goals: Learn to use phase portraits to visualize behavior of dynamical systems Understand different types of stability for an equilibrium
More informationA matrix is a rectangular array of. objects arranged in rows and columns. The objects are called the entries. is called the size of the matrix, and
Section 5.5. Matrices and Vectors A matrix is a rectangular array of objects arranged in rows and columns. The objects are called the entries. A matrix with m rows and n columns is called an m n matrix.
More informationAn Even Order Symmetric B Tensor is Positive Definite
An Even Order Symmetric B Tensor is Positive Definite Liqun Qi, Yisheng Song arxiv:1404.0452v4 [math.sp] 14 May 2014 October 17, 2018 Abstract It is easily checkable if a given tensor is a B tensor, or
More informationIntroduction to Modern Control MT 2016
CDT Autonomous and Intelligent Machines & Systems Introduction to Modern Control MT 2016 Alessandro Abate Lecture 2 First-order ordinary differential equations (ODE) Solution of a linear ODE Hints to nonlinear
More informationDisturbance Attenuation for a Class of Nonlinear Systems by Output Feedback
Disturbance Attenuation for a Class of Nonlinear Systems by Output Feedback Wei in Chunjiang Qian and Xianqing Huang Submitted to Systems & Control etters /5/ Abstract This paper studies the problem of
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationOn the Solution of the Van der Pol Equation
On the Solution of the Van der Pol Equation arxiv:1308.5965v1 [math-ph] 27 Aug 2013 Mayer Humi Department of Mathematical Sciences, Worcester Polytechnic Institute, 100 Institute Road, Worcester, MA 0l609
More informationSECTION 3.3. PROBLEM 22. The null space of a matrix A is: N(A) = {X : AX = 0}. Here are the calculations of AX for X = a,b,c,d, and e. =
SECTION 3.3. PROBLEM. The null space of a matrix A is: N(A) {X : AX }. Here are the calculations of AX for X a,b,c,d, and e. Aa [ ][ ] 3 3 [ ][ ] Ac 3 3 [ ] 3 3 [ ] 4+4 6+6 Ae [ ], Ab [ ][ ] 3 3 3 [ ]
More informationx x2 2 + x3 3 x4 3. Use the divided-difference method to find a polynomial of least degree that fits the values shown: (b)
Numerical Methods - PROBLEMS. The Taylor series, about the origin, for log( + x) is x x2 2 + x3 3 x4 4 + Find an upper bound on the magnitude of the truncation error on the interval x.5 when log( + x)
More informationCONTROL SYSTEMS ANALYSIS VIA BLIND SOURCE DECONVOLUTION. Kenji Sugimoto and Yoshito Kikkawa
CONTROL SYSTEMS ANALYSIS VIA LIND SOURCE DECONVOLUTION Kenji Sugimoto and Yoshito Kikkawa Nara Institute of Science and Technology Graduate School of Information Science 896-5 Takayama-cho, Ikoma-city,
More informationMath113: Linear Algebra. Beifang Chen
Math3: Linear Algebra Beifang Chen Spring 26 Contents Systems of Linear Equations 3 Systems of Linear Equations 3 Linear Systems 3 2 Geometric Interpretation 3 3 Matrices of Linear Systems 4 4 Elementary
More informationDiagonalization of Matrix
of Matrix King Saud University August 29, 2018 of Matrix Table of contents 1 2 of Matrix Definition If A M n (R) and λ R. We say that λ is an eigenvalue of the matrix A if there is X R n \ {0} such that
More informationAn homotopy method for exact tracking of nonlinear nonminimum phase systems: the example of the spherical inverted pendulum
9 American Control Conference Hyatt Regency Riverfront, St. Louis, MO, USA June -, 9 FrA.5 An homotopy method for exact tracking of nonlinear nonminimum phase systems: the example of the spherical inverted
More information8.385 MIT (Rosales) Hopf Bifurcations. 2 Contents Hopf bifurcation for second order scalar equations. 3. Reduction of general phase plane case to seco
8.385 MIT Hopf Bifurcations. Rodolfo R. Rosales Department of Mathematics Massachusetts Institute of Technology Cambridge, Massachusetts MA 239 September 25, 999 Abstract In two dimensions a Hopf bifurcation
More informationControl Systems Design, SC4026. SC4026 Fall 2010, dr. A. Abate, DCSC, TU Delft
Control Systems Design, SC426 SC426 Fall 2, dr A Abate, DCSC, TU Delft Lecture 5 Controllable Canonical and Observable Canonical Forms Stabilization by State Feedback State Estimation, Observer Design
More informationChapter 4 Euclid Space
Chapter 4 Euclid Space Inner Product Spaces Definition.. Let V be a real vector space over IR. A real inner product on V is a real valued function on V V, denoted by (, ), which satisfies () (x, y) = (y,
More informationCDS 101/110a: Lecture 2.1 Dynamic Behavior
CDS 11/11a: Lecture 2.1 Dynamic Behavior Richard M. Murray 6 October 28 Goals: Learn to use phase portraits to visualize behavior of dynamical systems Understand different types of stability for an equilibrium
More informationTwo Results About The Matrix Exponential
Two Results About The Matrix Exponential Hongguo Xu Abstract Two results about the matrix exponential are given. One is to characterize the matrices A which satisfy e A e AH = e AH e A, another is about
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 informationFeedback Linearization
Feedback Linearization Peter Al Hokayem and Eduardo Gallestey May 14, 2015 1 Introduction Consider a class o single-input-single-output (SISO) nonlinear systems o the orm ẋ = (x) + g(x)u (1) y = h(x) (2)
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ẋ = f(x, y), ẏ = g(x, y), (x, y) D, can only have periodic solutions if (f,g) changes sign in D or if (f,g)=0in D.
4 Periodic Solutions We have shown that in the case of an autonomous equation the periodic solutions correspond with closed orbits in phase-space. Autonomous two-dimensional systems with phase-space R
More informationIntro. Computer Control Systems: F8
Intro. Computer Control Systems: F8 Properties of state-space descriptions and feedback Dave Zachariah Dept. Information Technology, Div. Systems and Control 1 / 22 dave.zachariah@it.uu.se F7: Quiz! 2
More informationReview: control, feedback, etc. Today s topic: state-space models of systems; linearization
Plan of the Lecture Review: control, feedback, etc Today s topic: state-space models of systems; linearization Goal: a general framework that encompasses all examples of interest Once we have mastered
More information