On one Application of Newton s Method to Stability Problem
|
|
- Willis Houston
- 5 years ago
- Views:
Transcription
1 Journal of Multidciplinary Engineering Science Technology (JMEST) ISSN: Vol. Issue 5, December On one Application of Newton s Method to Stability Problem Şerife Yılmaz Department of Mathematics, Faculty of Science, Anadolu University, Eskehir 26470, Turkey serifeyilmaz@anadolu.edu.tr Abstract In th paper we consider stability problem for switched linear systems. Th problem can be formulated as a convex minimization problem. By modifying the cost functions we apply the vector-valued Newton s method. Keywords Switched system; Hurwitz stability; common quadratic Lyapunov function; Newton s method I. INTRODUCTION Let AA be nn nn real matrix. If all eigenvalues of AA lie in the open left half plane then AA said to be Hurwitz stable. Hurwitz stability of AA equivalent to the following: There exts positive definite symmetric matrix PP such that AA TT PP + PPPP < 0 () the symbol T sts for the transpose, the symbol < for negative definiteness. Hurwitz stability of AA implies the asymptotic stability of the zero solution of the linear system xx = AAAA (2) xx = xx(tt) R nn. If the matrix AA switches between NN matrices AA, AA 2,, AA NN, i.e. AA {AA, AA 2,, AA NN } then the obtained system xx = AAAA (3) called a switched system. Sufficient condition for the asymptotic stability of the zero solution of (3) the extence of quadratic Lyapunov function of the form VV(xx) = xx TT PPPP PP > 0 AA ii TT PP + PPAA ii < 0 (ii =,2,, NN). (4) The matrix PP called a common solution to the Lyapunov inequalities (4). The stability problem of linear switched systems has been investigated in a lot of works (see [-] references therein). The papers [3-] study theoretical results for the extence of a common solution to (4). The papers [2-4] consider numerical algorithms for a common positive definite solution in the case of extence. In th paper we apply Newton s root finding method for the numerical generating of a common positive definite solution to (4). II. LYAPUNOV EQUATIONS In th section, we consider the Lyapunov inequality () which equivalent to the following equation. AA TT PP + PPPP = QQ (5) QQ > 0. We are looking for a positive definite solution PP of (5). In the iteration steps, the obtained PP guaranteed to be symmetric. The following theorem shows that in the case of Hurwitz stability of AA th implies the positive definiteness of PP. Theorem. Assume that AA Hurwitz stable. If there exts a symmetric solution PP to (5) then PP > 0. Proof: Define PP = ee AATTtt QQee AAAA dddd 0 ee AAAA sts for the matrix exponential. (6) Since AA Hurwitz stable, the matrix AA TT also Hurwitz stable. ee AAAA ee AATT tt define exponential functions with exponents Re(λλ ii ) tt < 0 λλ ii are the eigenvalues of AA. Th implies that the integral in ( 6 ) well defined. The matrix PP symmetric, positive definite satfies the following relation AA TT PP + PP AA = QQ (see [5]). Then AA TT PP PP + PP PP AA = 0. Multiplying by ee AATT tt ee AAAA give 0 = ee AATTtt AA TT PP PP + PP PP AA ee AAAA = dd dddd eeaatttt PP PP ee AAAA JMESTN
2 Journal of Multidciplinary Engineering Science Technology (JMEST) ISSN: Vol. Issue 5, December The integration from 0 to yields ee AATTtt PP PP ee AAAA 0 = 0. Using the fact that ee AAAA 0, ee AATTtt 0 as tt we obtain 0 PP PP = 0 PP = PP > 0. We are looking for a iterative procedure for a common PP satfying ( 4 ). Theorem allows to guarantee positive definiteness of PP obtained at each step of iteration. III. MODIFIED NEWTON S METHOD Consider a differentiable function FF: R nn R nn the following equation FF(xx) = 0. (7) Here xx = (xx, xx 2,, xx nn ) TT R nn, FF(xxff (xx), ff 2 (xx),, ff nn (xx) TT. Denote the Jacobian matrix by JJ(xx), i.e. JJ(xx ff ii(xx) (ii, jj =,2,, nn). xx jj The Newton method a method for an approximate solution of (7) starting from a suitable initial point xx 0 the iteration defined by xx kk = xx kk JJ(xx kk ) FF(xx kk ) (kk =,2, ) Define rr = nn(nn+) 2 xx xx 2 xx nn xx 2 xx nn+ xx 2nn PP = PP(xx xx nn xx 2nn xx rr The matrix inequalities (4) are equivalent to ff ii (xx) = λλ max (AA ii TT PP + PPAA ii ) < 0 (ii =,2,, NN) (8) xx = (xx, xx 2,, xx rr ) TT R rr, λλ max ( ) sts for the maximal eigenvalue. In the case of simple maximum eigenvalue of AA TT PP(xx) + PP(xx)AA, the gradient of ff(xx) = λλ max (AA TT PP(xx) + PP(xx)AA) should be easily calculated. Indeed since the function xx AA TT PP(xx) + PP(xx)AA linear then AA TT PP(xx) + PP(xx)AA = xx jj QQ jj. rr jj= Then ff(xx) = (uu TT QQ uu,, uu TT QQ rr uu), sts for the gradient, uu the unit eigenvector corresponding to the maximum eigenvalues of AA TT PP(xx) + PP(xx)AA (see [2]). Proposition. The function ff ii (xx) convex for each ii. Proof: The relation PP AA ii TT PP + PPAA ii linear. On the other h for symmetric CC, the function CC λλ max (CC) convex [6]. ff ii (xx) convex as a composition of linear convex functions. The system (4) has a common solution PP > 0 if only if there exts xx R rr such that ff ii (xx ) < 0 (ii =,2,, NN). (9) In order to apply Newton s method instead of the minimization of the functions ff ii (xx), we consider the system of equations ff ii (xx) = 0 (ii =,2,, NN). Without loss of generality we can set rr = NN. Indeed if NN > rr, we can combine some function by using the operation maximum. For example if NN = rr + then define gg (xx) = max{ff (xx), ff 2 (xx)}, gg ii (xx) = ff ii+ (xx) (ii = 2,3,, NN). Th operation preserves convexity. If NN < rr we use the operation of duplication. Thus from the now we assume that rr = NN. Define FF = (ff, ff 2,, ff rr ) TT consider the equation FF(xx) = 0 (0) xx R rr. If we apply the classical Newton s method to (0) we obtain the trivial sequence PP kk 0, since the functions ff ii (xx) are positive homogenous. To avoid th we impose the condition trace(pp) =. The following proposition shows that th does not violate the generality. Proposition 2. Assume that PP > 0 AA TT PP + PPPP < 0. Then AA TT PP + PP AA < 0 PP = trace(pp) PP. Proof: From PP = pp iiii > 0 it follows that for all xx R nn, xx 0, xx TT PPPP > 0. Taking xx ii = (0,,0,,0,,0) TT we obtain pp iiii > 0. trace(pp) > 0 AA TT PP + PP AA = trace(pp) [AATT PP + PPPP] < 0. The condition trace(pp) = reduces the number of variables from rr to rr. To solve (0) the following algorithm suggested. Algorithm. ) Consider the equation (0). Take initial matrix QQ = diag(,2,2,,2) consider AA TT PP + PPAA = QQ. JMESTN
3 The solution of th matrix equation let be PP 0. Dividing PP 0 by trace(pp 0 ) gives the initial iteration xx 0. 2) Replace the functions ff ii (xx) by ff ii(xx) + (ii =,2,, rr ). Apply Newton s iteration 2trace(PP 0 ) xx kk = xx kk JJ(xx kk ) FF (xx kk ) FF = (ff, ff 2,, ff rr ). 3) If ff ii(xx kk ) < 0 (ii =,2,, rr ) for some kk then stop. Otherwe continue. Example. Consider the Hurwitz stable matrices AA = 4 8 AA 2 = ff (xx) = λλ max (AA TT PP(xx) + PP(xx)AA ), ff 2 (xx) = λλ max (AA 2 TT PP(xx) + PP(xx)AA 2 ) Hence PP(xx xx xx 2 xx 2 xx QQ = AA TT PP + PPAA = QQ (ii =,2) PP 0 = trace(pp 0 ) PP = take the initial point xx 0 = (0.729, 0.354) TT. For th point calculations give the following maximum eigenvalues its corresponding unit eigenvectors: λλ max (AA TT PP(xx 0 ) + PP(xx 0 )AA ) = 0.75 the maximum eigenvector: uu = (,0) TT, λλ max (AA 2 TT PP(xx 0 ) + PP(xx 0 )AA 2 ) = the maximum eigenvector: uu 2 = ( , 0.705) TT. ff (xx 0 ) = 0.75, ff 2 (xx 0 ) = ff (xx) xx=xx 0 = ( 2, 2), ff 2 (xx) xx=xx 0 = (2.988,0.96). the Jacobian matrix of FF(xx) = ff (xx), ff 2 (xx) TT at xx 0 Journal of Multidciplinary Engineering Science Technology (JMEST) ISSN: Vol. Issue 5, December JJ(xx xx = = After 3 steps, we get ff (xx 3 ) < 0 ff 2 (xx 3 ) < 0 (see Table I). Hence for the matrix PP = PP(xx , AA TT ii PP + PPAA ii < 0 (ii =,2) are satfied. Table I kk xx kk ff (xx kk ) ff 22 (xx kk ) (0.30,0.9) TT (0.367, 0.35) TT (0.390, 0.247) TT Example 2. Consider the Hurwitz stable matrices AA = 2 2, AA = 2 AA 3 = ff (xx) = max(λλ max (AA TT PP(xx) + PP(xx)AA ), λλ max (AA 2 TT PP(xx) + PP(xx)AA 2 )), ff 2 (xx) = λλ max (AA 3 TT PP(xx) + PP(xx)AA 3 ). PP(xx xx xx 2 xx 2 xx QQ = AA TT PP + PPAA = QQ (ii =,2) PP 0 = trace(pp 0 ) =. Take the initial point xx 0 = (0.46,0.083) TT. For th point calculations give the following maximum eigenvalues its corresponding unit eigenvectors: λλ max (AA TT PP(xx 0 ) + PP(xx 0 )AA ) =, uu = (,0) TT, λλ max (AA 2 TT PP(xx 0 ) + PP(xx 0 )AA 2 ) = 0.680, uu 2 = ( 0.082,0.996) TT, λλ max (AA 3 TT PP(xx 0 ) + PP(xx 0 )AA 3 ) = 0.567, uu 3 = ( 0.987,0.60) TT. JMESTN
4 ff (xx 0 ) = max( 2.05,0.074) = 0.680, ff 2 (xx 0 ) = ff (xx) xx=xx 0 = (.9, 5.767), ff 2 (xx) xx=xx 0 = (0.786, 3.478). The Jacobian matrix of FF(xxff (xx), ff 2 (xx) TT at xx 0 JJ(xx xx = = After 6 steps, we get ff (xx 6 ) < 0 ff 2 (xx 6 ) < 0 (see Table II). Hence for the matrix 5 PP = PP(xx , AA ii TT PP + PPAA ii < 0 (ii =,2,3) are satfied. Table II kk xx kk ff (xx kk ) ff 22 (xx kk ) (0.80,0.336) TT (0.58,2.025) TT (0.553,0.33) TT (0.293,0.99) TT (0.46,0.38) TT Example 3. Consider the Hurwitz stable matrices AA = 0 2, AA 2 = 6 3, AA 3 = 2 4 2, AA 4 = 3 3 4, AA 5 = ff ii (xx) = λλ max (AA ii TT PP(xx) + PP(xx)AA ii ) Journal of Multidciplinary Engineering Science Technology (JMEST) ISSN: Vol. Issue 5, December xx xx 2 xx 3 PP(xxxx 2 xx 3 xx 4 xx 5 xx 5 xx xx 4 Hence 0 0 QQ = AA TT PP + PPAA = QQ (ii =,2) PP 0 = trace(pp 0 ) PP 0 = the initial point xx 0 = (0.069, 0.88,0.037,0.803, 0.74) TT. For th point calculations give the following: FF(xx 0 ) = (.363,.952,0.408,.843,.909) TT, JJ(xx 0 ) = xx = ( 0.090,0.252,0.079,0.67, 0.26) TT. After 36 steps, we get xx 36 = (0.42, 0.38, 0.044,0.28, 0.028) TT FF(xx 36 ) = ( 0.590, 0.360, 0.608, 0.447, 0.64) TT. Hence AA ii TT PP + PPAA ii < 0 (ii =,2,3,4,5) PP = PP(xx REFERENCES [] H. Lin P.J. Antsakl, Stability Stabilizability of Switched Linear Systems: A Survey of Recent Results, IEEE Transactions on Automatic Control, 54(2): , [2] D. Liberzon, Switching in System Control. Birkhäuser, Boston, [3] K.S. Narendra J. Balakrhnan, A common Lyapunov function for stable LTI systems with commuting A-matrices, IEEE Transactions on Automatic Control, 39, , 994. [4] Y. Mori, T. Mori Y. Kuroe, A solution to the common Lyapunov function problem for continuous-time systems, In Proceedings of the 36 th JMESTN
5 [5] Conference on Decion Control, vol. 3, pp , 997. [6] D. Liberzon, J.P. Hespanha A.S. Morse, Stability of switched linear systems: a Lie-algebraic condition, Systems & Control Letters, 37(3):7 22, 999. [7] R. Shorten K. Narendra, Necessary sufficient conditions for the extence of a common quadratic Lyapunov function for two stable second order linear time-invariant systems, in Proc. Amer. Control Conf., 40 44, 999. [8] R. Shorten, O. Mason, F.O. Cairbre, P. Curran, A unifying framework for the SISO circle criterion other quadratic stability criteria, Int. J. Control, 77(): 8, [9] R. Shorten, K. Narendra, O. Mason, A result on common quadratic Lyapunov functions, IEEE Transactions on Automatic Control, 48():0 3, [0] L. Gurvits, R. Shorten, O. Mason, On the stability of switched positive linear systems, IEEE Transactions on Automatic Control, 52(6):099 03, Journal of Multidciplinary Engineering Science Technology (JMEST) ISSN: Vol. Issue 5, December [] T.J. Laffey H. Šmigoc, Tensor conditions for the extence of a common solution to the Lyapunov equation, Linear Algebra Appl., vol. 420, pp , [2] R. Shorten, F. Wirth, O. Mason, K. Wulff C. King, Stability criteria for switched hybrid systems SIAM Rev. 49(4): , [3] S. Boyd Q. Yang, Structured simultaneous Lyapunov functions for system stability problems, International Journal of Control, 49, , 989. [4] D. Liberzon R. Tempo, Common Lyapunov Functions Gradient Algorithms, IEEE Transactions on Automatic Control, 49, , [5] T. Büyükköroğlu, Ö. Esen V. Dzhafarov, Common Lyapunov functions for some special classes of stable systems, IEEE Transactions on Automatic Control, 56, , 20. [6] H.K. Khalil, Nonlinear Systems. Prentice Hall, New Jersey, [7] R.A. Horn C.R. Johnson, Matrix Analys. Cambridge University Press, 992. JMESTN
A note on simultaneous block diagonally stable matrices
Electronic Journal of Linear Algebra Volume 26 Volume 26 (2013) Article 7 2013 A note on simultaneous block diagonally stable matrices I. Bras A. C. Carapito carapito@ubi.pt P. Rocha Follow this and additional
More informationQuadratic and Copositive Lyapunov Functions and the Stability of Positive Switched Linear Systems
Proceedings of the 2007 American Control Conference Marriott Marquis Hotel at Times Square New York City, USA, July 11-13, 2007 WeA20.1 Quadratic and Copositive Lyapunov Functions and the Stability of
More informationOn common quadratic Lyapunov functions for stable discrete-time LTI systems
On common quadratic Lyapunov functions for stable discrete-time LTI systems Oliver Mason, Hamilton Institute, NUI Maynooth, Maynooth, Co. Kildare, Ireland. (oliver.mason@may.ie) Robert Shorten 1, Hamilton
More informationOn the simultaneous diagonal stability of a pair of positive linear systems
On the simultaneous diagonal stability of a pair of positive linear systems Oliver Mason Hamilton Institute NUI Maynooth Ireland Robert Shorten Hamilton Institute NUI Maynooth Ireland Abstract In this
More informationChapter 22 : Electric potential
Chapter 22 : Electric potential What is electric potential? How does it relate to potential energy? How does it relate to electric field? Some simple applications What does it mean when it says 1.5 Volts
More informationWorksheets for GCSE Mathematics. Algebraic Expressions. Mr Black 's Maths Resources for Teachers GCSE 1-9. Algebra
Worksheets for GCSE Mathematics Algebraic Expressions Mr Black 's Maths Resources for Teachers GCSE 1-9 Algebra Algebraic Expressions Worksheets Contents Differentiated Independent Learning Worksheets
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 informationHybrid Systems Course Lyapunov stability
Hybrid Systems Course Lyapunov stability OUTLINE Focus: stability of an equilibrium point continuous systems decribed by ordinary differential equations (brief review) hybrid automata OUTLINE Focus: stability
More informationA Delay-dependent Condition for the Exponential Stability of Switched Linear Systems with Time-varying Delay
A Delay-dependent Condition for the Exponential Stability of Switched Linear Systems with Time-varying Delay Kreangkri Ratchagit Department of Mathematics Faculty of Science Maejo University Chiang Mai
More information7.3 The Jacobi and Gauss-Seidel Iterative Methods
7.3 The Jacobi and Gauss-Seidel Iterative Methods 1 The Jacobi Method Two assumptions made on Jacobi Method: 1.The system given by aa 11 xx 1 + aa 12 xx 2 + aa 1nn xx nn = bb 1 aa 21 xx 1 + aa 22 xx 2
More informationGramians based model reduction for hybrid switched systems
Gramians based model reduction for hybrid switched systems Y. Chahlaoui Younes.Chahlaoui@manchester.ac.uk Centre for Interdisciplinary Computational and Dynamical Analysis (CICADA) School of Mathematics
More informationMathematics Ext 2. HSC 2014 Solutions. Suite 403, 410 Elizabeth St, Surry Hills NSW 2010 keystoneeducation.com.
Mathematics Ext HSC 4 Solutions Suite 43, 4 Elizabeth St, Surry Hills NSW info@keystoneeducation.com.au keystoneeducation.com.au Mathematics Extension : HSC 4 Solutions Contents Multiple Choice... 3 Question...
More informationSTABILITY OF PLANAR NONLINEAR SWITCHED SYSTEMS
LABORATOIRE INORMATIQUE, SINAUX ET SYSTÈMES DE SOPHIA ANTIPOLIS UMR 6070 STABILITY O PLANAR NONLINEAR SWITCHED SYSTEMS Ugo Boscain, régoire Charlot Projet TOpModel Rapport de recherche ISRN I3S/RR 2004-07
More informationEXPONENTIAL STABILITY OF SWITCHED LINEAR SYSTEMS WITH TIME-VARYING DELAY
Electronic Journal of Differential Equations, Vol. 2007(2007), No. 159, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) EXPONENTIAL
More informationDaniel Liberzon. Abstract. 1. Introduction. TuC11.6. Proceedings of the European Control Conference 2009 Budapest, Hungary, August 23 26, 2009
Proceedings of the European Control Conference 2009 Budapest, Hungary, August 23 26, 2009 TuC11.6 On New On new Sufficient sufficient Conditions conditionsfor forstability stability of switched Switched
More informationWorksheets for GCSE Mathematics. Quadratics. mr-mathematics.com Maths Resources for Teachers. Algebra
Worksheets for GCSE Mathematics Quadratics mr-mathematics.com Maths Resources for Teachers Algebra Quadratics Worksheets Contents Differentiated Independent Learning Worksheets Solving x + bx + c by factorisation
More informationOn Linear Copositive Lyapunov Functions and the Stability of Switched Positive Linear Systems
1 On Linear Copositive Lyapunov Functions and the Stability of Switched Positive Linear Systems O. Mason and R. Shorten Abstract We consider the problem of common linear copositive function existence for
More informationAnalysis and design of switched normal systems
Nonlinear Analysis 65 (2006) 2248 2259 www.elsevier.com/locate/na Analysis and design of switched normal systems Guisheng Zhai a,, Xuping Xu b, Hai Lin c, Anthony N. Michel c a Department of Mechanical
More informationSupport Vector Machines. CSE 4309 Machine Learning Vassilis Athitsos Computer Science and Engineering Department University of Texas at Arlington
Support Vector Machines CSE 4309 Machine Learning Vassilis Athitsos Computer Science and Engineering Department University of Texas at Arlington 1 A Linearly Separable Problem Consider the binary classification
More informationIndependent Component Analysis and FastICA. Copyright Changwei Xiong June last update: July 7, 2016
Independent Component Analysis and FastICA Copyright Changwei Xiong 016 June 016 last update: July 7, 016 TABLE OF CONTENTS Table of Contents...1 1. Introduction.... Independence by Non-gaussianity....1.
More informationOn Dwell Time Minimization for Switched Delay Systems: Free-Weighting Matrices Method
On Dwell Time Minimization for Switched Delay Systems: Free-Weighting Matrices Method Ahmet Taha Koru Akın Delibaşı and Hitay Özbay Abstract In this paper we present a quasi-convex minimization method
More informationSECTION 5: POWER FLOW. ESE 470 Energy Distribution Systems
SECTION 5: POWER FLOW ESE 470 Energy Distribution Systems 2 Introduction Nodal Analysis 3 Consider the following circuit Three voltage sources VV sss, VV sss, VV sss Generic branch impedances Could be
More informationStability of Equilibrium Positions of Mechanical Systems with Switched Force Fields
SCIETIFIC PUBLICATIOS OF THE STATE UIVERSITY OF OVI PAZAR SER. A: APPL. MATH. IFORM. AD MECH. vol. 4, 2 2012, 35-39 Stability of Equilibrium Positions of Mechanical Systems with Switched Force Fields A.
More informationLogarithmic Functions
Name Student ID Number Group Name Group Members Logarithmic Functions 1. Solve the equations below. xx = xx = 5. Were you able solve both equations above? If so, was one of the equations easier to solve
More informationExam Programme VWO Mathematics A
Exam Programme VWO Mathematics A The exam. The exam programme recognizes the following domains: Domain A Domain B Domain C Domain D Domain E Mathematical skills Algebra and systematic counting Relationships
More informationHybrid Systems - Lecture n. 3 Lyapunov stability
OUTLINE Focus: stability of equilibrium point Hybrid Systems - Lecture n. 3 Lyapunov stability Maria Prandini DEI - Politecnico di Milano E-mail: prandini@elet.polimi.it continuous systems decribed by
More informationAn alternative proof of the Barker, Berman, Plemmons (BBP) result on diagonal stability and extensions - Corrected Version
1 An alternative proof of the Barker, Berman, Plemmons (BBP) result on diagonal stability and extensions - Corrected Version Robert Shorten 1, Oliver Mason 1 and Christopher King 2 Abstract The original
More information2.4 Error Analysis for Iterative Methods
2.4 Error Analysis for Iterative Methods 1 Definition 2.7. Order of Convergence Suppose {pp nn } nn=0 is a sequence that converges to pp with pp nn pp for all nn. If positive constants λλ and αα exist
More informationMath 171 Spring 2017 Final Exam. Problem Worth
Math 171 Spring 2017 Final Exam Problem 1 2 3 4 5 6 7 8 9 10 11 Worth 9 6 6 5 9 8 5 8 8 8 10 12 13 14 15 16 17 18 19 20 21 22 Total 8 5 5 6 6 8 6 6 6 6 6 150 Last Name: First Name: Student ID: Section:
More informationVariations. ECE 6540, Lecture 02 Multivariate Random Variables & Linear Algebra
Variations ECE 6540, Lecture 02 Multivariate Random Variables & Linear Algebra Last Time Probability Density Functions Normal Distribution Expectation / Expectation of a function Independence Uncorrelated
More information10.4 The Cross Product
Math 172 Chapter 10B notes Page 1 of 9 10.4 The Cross Product The cross product, or vector product, is defined in 3 dimensions only. Let aa = aa 1, aa 2, aa 3 bb = bb 1, bb 2, bb 3 then aa bb = aa 2 bb
More informationMathematical Methods 2019 v1.2
Examination This sample has been compiled by the QCAA to model one possible approach to allocating marks in an examination. It matches the examination mark allocations as specified in the syllabus (~ 60%
More informationLesson 1: Successive Differences in Polynomials
Lesson 1 Lesson 1: Successive Differences in Polynomials Classwork Opening Exercise John noticed patterns in the arrangement of numbers in the table below. 2.4 3.4 4.4 5.4 6.4 5.76 11.56 19.36 29.16 40.96
More informationP.3 Division of Polynomials
00 section P3 P.3 Division of Polynomials In this section we will discuss dividing polynomials. The result of division of polynomials is not always a polynomial. For example, xx + 1 divided by xx becomes
More informationPolynomials and Polynomial Functions
1 Polynomials and Polynomial Functions One of the simplest types of algebraic expressions are polynomials. They are formed only by addition and multiplication of variables and constants. Since both addition
More information10.1 Three Dimensional Space
Math 172 Chapter 10A notes Page 1 of 12 10.1 Three Dimensional Space 2D space 0 xx.. xx-, 0 yy yy-, PP(xx, yy) [Fig. 1] Point PP represented by (xx, yy), an ordered pair of real nos. Set of all ordered
More informationMath 3 Unit 4: Rational Functions
Math Unit : Rational Functions Unit Title Standards. Equivalent Rational Expressions A.APR.6. Multiplying and Dividing Rational Expressions A.APR.7. Adding and Subtracting Rational Expressions A.APR.7.
More informationOn the stability of switched positive linear systems
On the stability of switched positive linear systems L. Gurvits, R. Shorten and O. Mason Abstract It was recently conjectured that the Hurwitz stability of the convex hull of a set of Metzler matrices
More informationStability and Stabilizability of Switched Linear Systems: A Short Survey of Recent Results
Proceedings of the 2005 IEEE International Symposium on Intelligent Control Limassol, Cyprus, June 27-29, 2005 MoA01-5 Stability and Stabilizability of Switched Linear Systems: A Short Survey of Recent
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 informationLecture No. 1 Introduction to Method of Weighted Residuals. Solve the differential equation L (u) = p(x) in V where L is a differential operator
Lecture No. 1 Introduction to Method of Weighted Residuals Solve the differential equation L (u) = p(x) in V where L is a differential operator with boundary conditions S(u) = g(x) on Γ where S is a differential
More informationLarge Scale Data Analysis Using Deep Learning
Large Scale Data Analysis Using Deep Learning Linear Algebra U Kang Seoul National University U Kang 1 In This Lecture Overview of linear algebra (but, not a comprehensive survey) Focused on the subset
More informationL 2 -induced Gains of Switched Systems and Classes of Switching Signals
L 2 -induced Gains of Switched Systems and Classes of Switching Signals Kenji Hirata and João P. Hespanha Abstract This paper addresses the L 2-induced gain analysis for switched linear systems. We exploit
More informationApproximate Second Order Algorithms. Seo Taek Kong, Nithin Tangellamudi, Zhikai Guo
Approximate Second Order Algorithms Seo Taek Kong, Nithin Tangellamudi, Zhikai Guo Why Second Order Algorithms? Invariant under affine transformations e.g. stretching a function preserves the convergence
More informationP.2 Multiplication of Polynomials
1 P.2 Multiplication of Polynomials aa + bb aa + bb As shown in the previous section, addition and subtraction of polynomials results in another polynomial. This means that the set of polynomials is closed
More informationIntegrated Attitude Determination and Control System via Magnetic Measurements and Actuation
Proceedings of the th WSEAS International Conference on Automatic Control, Modelling and Simulation Integrated Attitude Determination and Control System via Magnetic Measurements and Actuation MOHAMMAD
More informationSWITCHED SYSTEMS THAT ARE PERIODICALLY STABLE MAY BE UNSTABLE
SWITCHED SYSTEMS THAT ARE PERIODICALLY STABLE MAY BE UNSTABLE VINCENT D. BLONDEL, JACQUES THEYS, AND ALEXANDER A. VLADIMIROV 1. Introduction In this contribution, we prove the existence of switched linear
More informationA Posteriori Error Estimates For Discontinuous Galerkin Methods Using Non-polynomial Basis Functions
Lin Lin A Posteriori DG using Non-Polynomial Basis 1 A Posteriori Error Estimates For Discontinuous Galerkin Methods Using Non-polynomial Basis Functions Lin Lin Department of Mathematics, UC Berkeley;
More informationOn optimal quadratic Lyapunov functions for polynomial systems
On optimal quadratic Lyapunov functions for polynomial systems G. Chesi 1,A.Tesi 2, A. Vicino 1 1 Dipartimento di Ingegneria dell Informazione, Università disiena Via Roma 56, 53100 Siena, Italy 2 Dipartimento
More informationEstimation of cumulative distribution function with spline functions
INTERNATIONAL JOURNAL OF ECONOMICS AND STATISTICS Volume 5, 017 Estimation of cumulative distribution function with functions Akhlitdin Nizamitdinov, Aladdin Shamilov Abstract The estimation of the cumulative
More information3.2 A2 - Just Like Derivatives but Backwards
3. A - Just Like Derivatives but Backwards The Definite Integral In the previous lesson, you saw that as the number of rectangles got larger and larger, the values of Ln, Mn, and Rn all grew closer and
More informationWork, Energy, and Power. Chapter 6 of Essential University Physics, Richard Wolfson, 3 rd Edition
Work, Energy, and Power Chapter 6 of Essential University Physics, Richard Wolfson, 3 rd Edition 1 With the knowledge we got so far, we can handle the situation on the left but not the one on the right.
More informationRow and Column Representatives in Qualitative Analysis of Arbitrary Switching Positive Systems
ROMANIAN JOURNAL OF INFORMATION SCIENCE AND TECHNOLOGY Volume 19, Numbers 1-2, 2016, 127 136 Row and Column Representatives in Qualitative Analysis of Arbitrary Switching Positive Systems Octavian PASTRAVANU,
More informationSTABILITY OF 2D SWITCHED LINEAR SYSTEMS WITH CONIC SWITCHING. 1. Introduction
STABILITY OF D SWITCHED LINEAR SYSTEMS WITH CONIC SWITCHING H LENS, M E BROUCKE, AND A ARAPOSTATHIS 1 Introduction The obective of this paper is to explore necessary and sufficient conditions for stability
More informationM.5 Modeling the Effect of Functional Responses
M.5 Modeling the Effect of Functional Responses The functional response is referred to the predation rate as a function of the number of prey per predator. It is recognized that as the number of prey increases,
More informationF.1 Greatest Common Factor and Factoring by Grouping
section F1 214 is the reverse process of multiplication. polynomials in algebra has similar role as factoring numbers in arithmetic. Any number can be expressed as a product of prime numbers. For example,
More informationAnalysis of different Lyapunov function constructions for interconnected hybrid systems
Analysis of different Lyapunov function constructions for interconnected hybrid systems Guosong Yang 1 Daniel Liberzon 1 Andrii Mironchenko 2 1 Coordinated Science Laboratory University of Illinois at
More informationOn The Cauchy Problem For Some Parabolic Fractional Partial Differential Equations With Time Delays
Journal of Mathematics and System Science 6 (216) 194-199 doi: 1.17265/2159-5291/216.5.3 D DAVID PUBLISHING On The Cauchy Problem For Some Parabolic Fractional Partial Differential Equations With Time
More informationTime Domain Analysis of Linear Systems Ch2. University of Central Oklahoma Dr. Mohamed Bingabr
Time Domain Analysis of Linear Systems Ch2 University of Central Oklahoma Dr. Mohamed Bingabr Outline Zero-input Response Impulse Response h(t) Convolution Zero-State Response System Stability System Response
More informationGLOBAL ANALYSIS OF PIECEWISE LINEAR SYSTEMS USING IMPACT MAPS AND QUADRATIC SURFACE LYAPUNOV FUNCTIONS
GLOBAL ANALYSIS OF PIECEWISE LINEAR SYSTEMS USING IMPACT MAPS AND QUADRATIC SURFACE LYAPUNOV FUNCTIONS Jorge M. Gonçalves, Alexandre Megretski y, Munther A. Dahleh y California Institute of Technology
More informationSwitched systems: stability
Switched systems: stability OUTLINE Switched Systems Stability of Switched Systems OUTLINE Switched Systems Stability of Switched Systems a family of systems SWITCHED SYSTEMS SWITCHED SYSTEMS a family
More informationDefinition: A sequence is a function from a subset of the integers (usually either the set
Math 3336 Section 2.4 Sequences and Summations Sequences Geometric Progression Arithmetic Progression Recurrence Relation Fibonacci Sequence Summations Definition: A sequence is a function from a subset
More informationEE 341 A Few Review Problems Fall 2018
EE 341 A Few Review Problems Fall 2018 1. A linear, time-invariant system with input, xx(tt), and output, yy(tt), has the unit-step response, cc(tt), shown below. a. Find, and accurately sketch, the system
More informationDisturbance Attenuation Properties for Discrete-Time Uncertain Switched Linear Systems
Disturbance Attenuation Properties for Discrete-Time Uncertain Switched Linear Systems Hai Lin Department of Electrical Engineering University of Notre Dame Notre Dame, IN 46556, USA Panos J. Antsaklis
More informationConvergence Rate of Nonlinear Switched Systems
Convergence Rate of Nonlinear Switched Systems Philippe JOUAN and Saïd NACIRI arxiv:1511.01737v1 [math.oc] 5 Nov 2015 January 23, 2018 Abstract This paper is concerned with the convergence rate of the
More informationOn the stability of switched positive linear systems
1 On the stability of switched positive linear systems L. Gurvits, R. Shorten and O. Mason Abstract It was recently conjectured that the Hurwitz stability of the convex hull of a set of Metzler matrices
More informationSolving Fuzzy Nonlinear Equations by a General Iterative Method
2062062062062060 Journal of Uncertain Systems Vol.4, No.3, pp.206-25, 200 Online at: www.jus.org.uk Solving Fuzzy Nonlinear Equations by a General Iterative Method Anjeli Garg, S.R. Singh * Department
More informationClassical RSA algorithm
Classical RSA algorithm We need to discuss some mathematics (number theory) first Modulo-NN arithmetic (modular arithmetic, clock arithmetic) 9 (mod 7) 4 3 5 (mod 7) congruent (I will also use = instead
More informationIntroduction to the General operation of Matrices
Introduction to the General operation of Matrices CLAUDE ZIAD BAYEH 1, 1 Faculty of Engineering II, Lebanese University EGRDI transaction on mathematics (004) LEBANON Email: claude_bayeh_cbegrdi@hotmail.com
More informationOn Several Composite Quadratic Lyapunov Functions for Switched Systems
On Several Composite Quadratic Lyapunov Functions for Switched Systems Tingshu Hu, Liqiang Ma and Zongli Lin Abstract Three types of composite quadratic Lyapunov funtions are used for deriving conditions
More informationdue date: third day of class estimated time: 10 hours (for planning purposes only; work until you finish)
Honors PreCalculus Summer Work 016 due date: third day of class estimated time: 10 hours (for planning purposes only; work until you finish) Dear Honors PreCalculus Students, This assignment is designed
More informationA note on the Adomian decomposition method for generalized Burgers-Fisher equation
INTERNATIONAL JOURNAL OF MATHEMATICS AND COMPUTERS IN SIMULATION Volume 11, 017 A note on the Adomian decomposition method for generalized Burgers-Fisher equation M. Meštrović, E. Ocvirk and D. Kunštek
More information2. Linear algebra. matrices and vectors. linear equations. range and nullspace of matrices. function of vectors, gradient and Hessian
FE661 - Statistical Methods for Financial Engineering 2. Linear algebra Jitkomut Songsiri matrices and vectors linear equations range and nullspace of matrices function of vectors, gradient and Hessian
More informationMath 3 Unit 3: Polynomial Functions
Math 3 Unit 3: Polynomial Functions Unit Title Standards 3.1 End Behavior of Polynomial Functions F.IF.7c 3.2 Graphing Polynomial Functions F.IF.7c, A.APR3 3.3 Writing Equations of Polynomial Functions
More informationLimit Cycles in High-Resolution Quantized Feedback Systems
Limit Cycles in High-Resolution Quantized Feedback Systems Li Hong Idris Lim School of Engineering University of Glasgow Glasgow, United Kingdom LiHonIdris.Lim@glasgow.ac.uk Ai Poh Loh Department of Electrical
More informationPreliminary results on the stability of switched positive linear systems
1 Preliminary results on the stability of switched positive linear systems L. Gurvits, R. Shorten and O. Mason Abstract It was recently conjectured that the Hurwitz stability of a polytope of Metzler matrices
More informationSection 2: Equations and Inequalities
Topic 1: Equations: True or False?... 29 Topic 2: Identifying Properties When Solving Equations... 31 Topic 3: Solving Equations... 34 Topic 4: Solving Equations Using the Zero Product Property... 36 Topic
More informationInterval Matrix Systems and Invariance of Non-Symmetrical Contractive Sets
ROMANIAN JOURNAL OF INFORMATION SCIENCE AND TECHNOLOGY Volume 19, Numbers 1-2, 2016, 137 147 Interval Matrix Systems and Invariance of Non-Symmetrical Contractive Sets Octavian PASTRAVANU, Mihaela-Hanako
More informationLesson 14: Solving Inequalities
Hart Interactive Algebra 1 Lesson 14 Classwork 1. Consider the inequality xx 2 + 4xx 5. a. Think about some possible values to assign to xx that make this inequality a true statement. Find at least two
More informationQuadratic Equations and Functions
50 Quadratic Equations and Functions In this chapter, we discuss various ways of solving quadratic equations, aaxx 2 + bbbb + cc 0, including equations quadratic in form, such as xx 2 + xx 1 20 0, and
More informationStability of Switched Linear Hyperbolic Systems by Lyapunov Techniques
2196 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 59, NO. 8, AUGUST 2014 Stability of Switched Linear Hyperbolic Systems by Lyapunov Techniques Christophe Prieur, Antoine Girard, Emmanuel Witrant Abstract
More informationADAPTIVE NEURAL NETWORK CONTROLLER DESIGN FOR BLENDED-WING UAV WITH COMPLEX DAMAGE
ADAPTIVE NEURAL NETWORK CONTROLLER DESIGN FOR BLENDED-WING UAV WITH COMPLEX DAMAGE Kijoon Kim*, Jongmin Ahn**, Seungkeun Kim*, Jinyoung Suk* *Chungnam National University, **Agency for Defense and Development
More informationMath 95--Review Prealgebra--page 1
Math 95--Review Prealgebra--page 1 Name Date In order to do well in algebra, there are many different ideas from prealgebra that you MUST know. Some of the main ideas follow. If these ideas are not just
More informationControl of Mobile Robots
Control of Mobile Robots Regulation and trajectory tracking Prof. Luca Bascetta (luca.bascetta@polimi.it) Politecnico di Milano Dipartimento di Elettronica, Informazione e Bioingegneria Organization and
More information20-Apr-2009, 15:21. The conference website can be found at
This is a preprint of an article accepted for presentation at POSTA 09 Third Multidisciplinary International Symposium on Positive Systems: Theory and Application c 2009 Springer. The conference website
More informationMath 3 Unit 3: Polynomial Functions
Math 3 Unit 3: Polynomial Functions Unit Title Standards 3.1 End Behavior of Polynomial Functions F.IF.7c 3.2 Graphing Polynomial Functions F.IF.7c, A.APR3 3.3 Writing Equations of Polynomial Functions
More informationCS Lecture 8 & 9. Lagrange Multipliers & Varitional Bounds
CS 6347 Lecture 8 & 9 Lagrange Multipliers & Varitional Bounds General Optimization subject to: min ff 0() R nn ff ii 0, h ii = 0, ii = 1,, mm ii = 1,, pp 2 General Optimization subject to: min ff 0()
More informationHeat, Work, and the First Law of Thermodynamics. Chapter 18 of Essential University Physics, Richard Wolfson, 3 rd Edition
Heat, Work, and the First Law of Thermodynamics Chapter 18 of Essential University Physics, Richard Wolfson, 3 rd Edition 1 Different ways to increase the internal energy of system: 2 Joule s apparatus
More informationLesson 7: Algebraic Expressions The Commutative and Associative Properties
Algebraic Expressions The Commutative and Associative Properties Classwork Exercise 1 Suzy draws the following picture to represent the sum 3 + 4: Ben looks at this picture from the opposite side of the
More informationLyapunov Stability of Linear Predictor Feedback for Distributed Input Delays
IEEE TRANSACTIONS ON AUTOMATIC CONTROL VOL. 56 NO. 3 MARCH 2011 655 Lyapunov Stability of Linear Predictor Feedback for Distributed Input Delays Nikolaos Bekiaris-Liberis Miroslav Krstic In this case system
More informationLyapunov-based control of quantum systems
Lyapunov-based control of quantum systems Symeon Grivopoulos Bassam Bamieh Department of Mechanical and Environmental Engineering University of California, Santa Barbara, CA 936-57 symeon,bamieh@engineering.ucsb.edu
More informationTransition to College Math and Statistics
Transition to College Math and Statistics Summer Work 016 due date: third day of class estimated time: 10 hours (for planning purposes only; work until you finish) Dear College Algebra Students, This assignment
More informationA Lie-algebraic condition for stability of switched nonlinear systems
43rd IEEE Conference on Decision and Control December 14-17, 2004 Atlantis, Paradise Island, Bahamas FrB01.6 A Lie-algebraic condition for stability of switched nonlinear systems Michael Margaliot and
More informationInternational Journal of PharmTech Research CODEN (USA): IJPRIF, ISSN: Vol.8, No.7, pp , 2015
International Journal of PharmTech Research CODEN (USA): IJPRIF, ISSN: 0974-4304 Vol.8, No.7, pp 99-, 05 Lotka-Volterra Two-Species Mutualistic Biology Models and Their Ecological Monitoring Sundarapandian
More informationUnit II. Page 1 of 12
Unit II (1) Basic Terminology: (i) Exhaustive Events: A set of events is said to be exhaustive, if it includes all the possible events. For example, in tossing a coin there are two exhaustive cases either
More informationFast Linear Iterations for Distributed Averaging 1
Fast Linear Iterations for Distributed Averaging 1 Lin Xiao Stephen Boyd Information Systems Laboratory, Stanford University Stanford, CA 943-91 lxiao@stanford.edu, boyd@stanford.edu Abstract We consider
More informationSecondary 3H Unit = 1 = 7. Lesson 3.3 Worksheet. Simplify: Lesson 3.6 Worksheet
Secondary H Unit Lesson Worksheet Simplify: mm + 2 mm 2 4 mm+6 mm + 2 mm 2 mm 20 mm+4 5 2 9+20 2 0+25 4 +2 2 + 2 8 2 6 5. 2 yy 2 + yy 6. +2 + 5 2 2 2 0 Lesson 6 Worksheet List all asymptotes, holes and
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 informationGlobal stabilization of feedforward systems with exponentially unstable Jacobian linearization
Global stabilization of feedforward systems with exponentially unstable Jacobian linearization F Grognard, R Sepulchre, G Bastin Center for Systems Engineering and Applied Mechanics Université catholique
More informationMath, Stats, and Mathstats Review ECONOMETRICS (ECON 360) BEN VAN KAMMEN, PHD
Math, Stats, and Mathstats Review ECONOMETRICS (ECON 360) BEN VAN KAMMEN, PHD Outline These preliminaries serve to signal to students what tools they need to know to succeed in ECON 360 and refresh their
More information