A fast algorithm to compute the controllability, decentralized fixed-mode, and minimum-phase radius of LTI systems
|
|
- Moris Bruce
- 5 years ago
- Views:
Transcription
1 Proceedings of the 7th IEEE Conference on Decision Control Cancun, Mexico, Dec. 9-, 8 TuA6. A fast algorithm to compute the controllability, decentralized fixed-mode, minimum-phase radius of LTI systems Simon Lam Edward J. Davison Abstract In this paper, an efficient algorithm is presented for solving the -D optimization problem which is associated with computing the robustness measures of a number of linear time-invariant LTI) system properties such as the controllability radius, the decentralized fixed-mode radius, the minimum-phase radius. Unlike methods such as gradient search methods, the proposed algorithm works well independent of the initial trial point, obtains the minimum with a high level of confidence that it is indeed global. A numerical example is included. I. INTRODUCTION Consider the following LTI multivariable system ẋ = Ax + Bu y = Cx + Du ) where x R n, u R m, y R r are respectively the state, input, output vectors, A, B, C, D are constant matrices with the appropriate dimensions for n, m, r, maxr,m) n. It is well known that there exists a LTI controller that can assign the eigenvalues of the closed-loop system to any arbitrary spectrum if only if the system is controllable observable. However, when a controllable system is subject to parametric perturbations i.e. A A + A B B + B ), the system may be very close to becoming uncontrollable. Hence, a continuous controllability measure is more informative than the traditional yes/no controllability metric, which simply determines whether a system is controllable or not. The same can be said about other system properties such as observability, stability, minimum-phase, etc. In the current literature, various continuous measures have been proposed to measure the robustness of various system properties with respect to parametric perturbations. In particular, the controllability radius introduced in [, [ measures how close a controllable system is to being uncontrollable, the stability radius in [ measures how close a stable system is to an unstable one, the decentralized fixed-mode DFM) radius, [) measures how close a system with no DFMs is to having one. More recently, such robustness measures have been extended to characterize the robustness of a system s transmission zero properties. For instance, the minimum-phase radius in [6 measures how close a minimum-phase system is to becoming non-minimum phase. This work has been supported by NSERC under grant No. A96. S. Lam E. J. Davison is with the Systems Control Group, Department of Electrical Computer Engineering, University of Toronto, Toronto, ON, Canada MS G {simon,ted}@control.utoronto.ca Computing the stability radius mentioned above ) requires solving only a -D optimization problem, a fast algorithm to solve this is proposed in [7. On the other h, computing the other radii mentioned above i.e. except the stability radius) requires solving a -D optimization problem. As of now, the -D problem is solved using nonlinear search methods such as gradient search methods. A disadvantage with using such nonlinear techniques is that the minimum obtained may be a local minimum. Hence selecting the initial trial points are critical, a large number of initial points are generally tested in order to achieve a certain level of confidence that the obtained minimum is global. In this paper, we propose a more efficient algorithm for solving the -D optimization problem. The proposed algorithm has some similar features as the one in [7 in particular, the use of a so-called minimizing set ), has the advantage that it is not dependent on the initial test values, obtains the global minimum with a high degree of confidence. The paper is organized as follows. In Section II, the controllability radius is reviewed the -D optimization problem to be solved is described. Then in Section III, the few useful tools are developed, followed by a thorough description of the algorithm in Section IV. Finally, Section V gives a numerical example to demonstrate the effectiveness of the algorithm. II. NOTATION AND BACKGROUND In this paper, the field of real complex numbers are denoted by R C respectively. C + C U denotes the closed right half the closed upper half of the complex plane respectively, R + denotes the nonnegative real line [,+ ). The i-th singular value of a matrix M C p m is denoted by σ i M), where σ M) σ M). M denotes the spectral norm of a matrix M is equal to σ M). Also, M, M T, M H, M + denote respectively the complex conjugate, transpose, complex conjugate transpose, Moore-Penrose pseudoinverse of M. The real imaginary components of the matrix M are given by RM IM respectively. Finally, the set of eigenvalues of a square matrix A C n n is denoted by λa). A. Controllability radius The following definition is made. Definition.: Given a LTI system ), the controllability radius, rf c, is defined to be: r c FA,B) = inf{ [ A, B A F n n, B F n m, A + A,B + B ) is uncontrollable} ) /8/$. 8 IEEE 8
2 7th IEEE CDC, Cancun, Mexico, Dec. 9-, 8 TuA6. where F {C, R}. Similarly, the stabilizability radius, r c+ F, is defined as: r c+ F A,B) = inf{ [ A, B A F n n, B F n m, A + A,B + B ) is unstabilizable } ) In [ [, the complex real controllability radii are respectively shown to be given by: r c CA,B) = min s C σ n A sin B ) ) rra,b) c = min τ n A sin B ) ) s C where for W C p q, ) Re W γ Im W τ n W) := sup σ n γ, γ Im W Re W 6) For notational convenience in the paper, denote [ Re W γ ImW Pγ,W) := γ 7) ImW Re W Similarly, the complex real stabilizability radii are respectively shown to be: r c+ C A,B) = min σ n A sin B ) 8) s C + r c+ R A,B) = min τ n A sin B ) 9) s C + B. -D optimization problem It can be seen from ) 9) that a -D optimization problem in the complex plane is required to be solved when computing the complex real controllability/stabilizability radii. A similar -D problem is to be solved when computing the DFM radius, [), the minimum-phase radius 6). In this paper, a fast algorithm is proposed to solve the two stard general problems associated with such -D optimization problems; i.e. for given constant real matrices A,B,C,D) i {,...,n + minr,m)}: Problem : Find rcc,a,b,d) i A si B = min σ i ) s C Problem : Find rrc,a,b,d) i A si B = min τ i ) s C [ A si B For the sake of simplicity, denote Gs) :=. To avoid the trivial case where the radius is zero for all s C, we will assume that the two problems ) ) satisfy the following assumption: A B rank i ) To distinguish which field is consider, rc c i.e. F = C) is called the complex controllability radius rr c is called the real controllability radius. τ nw) is referred to as the n-th real perturbation value of W see [8). III. PRELIMINARY RESULTS Before presenting an efficient algorithm for solving ) ), some preliminary tools are needed; in particular Theorem. &. given below. Definition. The H matrix): Given real x, A R n n, B R n m, C R r n, D R r m, nonsingular matrices C n n, E C n+r) n+r), F C n+m) n+m), define: Hx,A,B,C,D,,E,F) := à B D ) C ) [ [ A B where à = A T + T, B = C T + T, [ [ C D C = B T + T, D = D T + T, [ [ T T = H, where T := satisfies T T [ [ T T x E = P H E ) T T l x FF H) P r I n for the permutation matrices P l = I n I r I m P r = I n I r I n I m. The following preliminary result [ is obtained. A s B Theorem.: Let M := E F, where s R, E C n+r) n+r), F C n+m) n+m), C n n are all nonsingular. Then, for given x R +, x σm) s λhx,a,b,c,d,,e,f)) ) Proof: M x σm) x λ M H xin+r M det M H = xi n+m [ ) xin+r M det P l P M H xi r = n+m A s B det A T s H C T C D B T D T [ ) T T + = T T à s B det = ) C D ) ) det D detã s B D C = s λ à B D )) C 9
3 7th IEEE CDC, Cancun, Mexico, Dec. 9-, 8 TuA6. Remark.: If D in ) is singular, then by ), ) can be modified such that s in ) is[ a generalized ) eigenvalue of the matrix pair à B C D, instead; [ à B i.e. s is in the set of λ C such that C D x = [ λ x for some non-zero eigenvector x C n+m+r. The following result is also obtained. Theorem.: Let γ,, θ [,π real x be given. Then, for all s R +, x σ P γ,g se iθ))) )) s λ H x,â, ˆB,Ĉ, ˆD, ˆ,T γ,n+rp nr,p nm T γ,n+m [ [ [ A B C where  =, ˆB =, Ĉ =, [ A [ B C D e ˆD =, ˆ iθ I = n D e iθ, P I nr = n I n I r I n I r T γ,n := γ, P mr = I n I n I m I m [ + γ I n iγ + γ I n + γ I n iγ + γ I n., Proof: It can easily be verified e.g. see [8) that T γ,n+r P γ,g se iθ)) [ ) G se Tγ,n+m iθ = Gse iθ ) A se iθ I n B = P nr A se iθ I n B C D P nm C D Therefore, P γ,g se iθ)) [ = Tγ,n+rP  sˆ nr Ĉ ˆB ˆD P nm T γ,n+m the proof follows immediately from Theorem.. IV. ALGORITHM An efficient iterative algorithm for solving ) ) for given constant real matrices A,B,C,D) i {,...,n + minm,r)} is now presented in this section. The following development focuses mainly on solving ), but it will be shown later that with a slight modification, the algorithm also applies to solving ). Firstly, let the global minimum of ) be denoted by: r = min τ igs)) = min s C s C sup γ, σ i Pγ,Gs))) is achieved at s C γ,, i.e. r = σ i Pγ,Gs ))). Now suppose that at the k-th iteration, for k =,,..., we are given r k, γ k, s k i.e. obtained from the previous iteration), which are approximations of r, γ, s respectively, where r k = σ i Pγ k,gs k ))). Furthermore, suppose we are also given the so-called minimizing set, S k, which is a closed) set in the complex plane that contains the global minimizer s. The basic idea of the algorithm is then to use the information of r k, γ k, s k, together with Theorem., to reduce the size of S k ; i.e. to reduce the region containing the global minimizer s. Denoting the reduced set by S k S k ), we then search within S k for a new point s k S k such that τ i Gs k )) < r k assign r k = τ i Gs k )). Also, assign γ k to be the value of γ that achieves r k ; i.e. r k = σ i Pγ k,gs k ))). Note that since r k < r k, then r k is a new better approximation of the global minimum r. The procedure is then repeated until either the size of S k is smaller than an user-specified tolerance, or until r k becomes very close to zero more on this later). A. The minimizing set, S k In the algorithm, the minimizing set is defined to be the union of particular sectors of interest in the closed upper half of the complex plane, can be described by S k [,π. The reason why the minimizing set is constrained to the closed upper half of the complex plane is due to the following e.g. see [6): A si σ i C A si τ i C ) B D ) B D A si B = σ i A si B = τ i ) ) In other words, the radii in ) ) are both symmetrical with respect to the ; hence one only needs to search within either the closed upper or lower half of the complex plane for the global minimum. It should be pointed out that S k is a subset of [,π, that the underlying minimizing set actually consists of the sectors in the complex plane described by S k ; i.e. the minimizing set is {s C s S k }, where s denotes the angle of s C. For convenience, however, we will sometimes refer to the minimizing set as S k. To reduce the size of S k based on given r k, γ k, s k, where r k = σ i Pγ k,gs k ))), we first construct the set R k : R k = {s C U σ i Pγ k,gs))) < r k } The significance of R k is that for all points not in R k, the radius can never be smaller than r k. This is because for all s / R k, τ i Gs)) = sup γ, σ i Pγ,Gs))) σ i Pγ k,gs))) r k Hence, one only needs to search within R k to find points in the complex plane that achieve a radius smaller than r k. Note that R k can also be a tighter) minimizing set. For ease of implementation, however, we chose the minimizing set to be the sectors containing R k i.e. as described by S k ) instead.
4 7th IEEE CDC, Cancun, Mexico, Dec. 9-, 8 TuA6. To obtain R k, we first fix a particular θ [,π, then find the set R θ k = { w R + σ i P γk,g we iθ))) < r k } 6) Therefore, R k = θ [,π R θ k 7) Note that Rk θ for a particular θ [,π can be obtained by applying Theorem.. In particular, by { Theorem., } the values along the ray, R θ := se iθ C s R +, that achieve a radius equal to r k are among the real nonnegative eigenvalues) of H r k,â, ˆB,Ĉ, ˆD, ˆ,T γ k,n+rp nr,p nm T γk,n+m. Hence, by solving for Λ )) = λ H r k,â, ˆB,Ĉ, ˆD, ˆ,T γ k,n+rp nr,p nm T γk,n+m then finding all real nonnegative s Λ such that σ i P γk,g se iθ))) = r k 8) we obtain the endpoints of the intervals along the ray R θ where σ i P γk,g se iθ))) < r k where σ i P γk,g se iθ))) > r k. To determine which interval is the former i.e. the one of interest), one can simply pick a trial point, w, within a particular interval evaluate σ i P γk,g we iθ))). If there are no real nonnegative eigenvalues, s Λ, such that σ i P γk,g se iθ))) = r k, then either Rk θ = R +, or Rk θ =. Again, a straightforward test using a trial point can confirm which is true. If Rk θ = is true for a particular θ [,π, then this implies that all the values along the ray R θ achieve a radius larger than r k, hence θ can be eliminated from the minimizing set, S k. Therefore the minimizing set for the next iteration is updated as follows. Let Θ k [,π be the set of intervals that describes the smallest union of sectors that contains R k ; i.e. Then update S k = S k Θ k. Θ k = { θ [,π R θ k } 9) B. Updating the current minimum, r k To obtain a radius smaller than r k i.e. a better approximation of the global minimum r ), one can perform a search e.g. a grid or rom search) within the set { R θ k θ S k }. From experience though, it is found that one can often obtain r k by letting s k = w m e iθm i.e. r k = τ i Gs k ))), where θ m is the midpoint of the largest interval of S k, w m is the midpoint of the largest interval of R θm k. C. Stopping criteria The algorithm has two main stopping criteria. Firstly, the algorithm stops when the size of the minimizing set, S k, the size of R k in 7) are both respectively smaller than user-specified tolerances, TOL Sk TOL Rk, where the size of S k is chosen to be the length of the largest interval in S k, the size of R k is chosen to be: sizer k ) = max θ { length of largest interval in R θ k } The second stopping criteria hles the special case when the radius is zero in spite of the fact that assumption ) is satisfied. For example, in [6, it is shown that the complex real minimum-phase radius are both when D = ; i.e. r n+minr,m) F A,B,C,) = where F {C, R}, this is achieved as real) s. Therefore, to prevent the algorithm from searching off into infinity, the second stopping criteria is added to quit when r k is smaller than an user-specified tolerance, TOL rk. D. Algorithm outline The algorithm can be summarized as follows: Algorithm.: Input: A,B,C,D) i, where ) is satisfied Input tolerances: TOL Sk, TOL Rk TOL rk Output: r, s, γ, where r = σ i Pγ,Gs ))) ) Initialization: Set S = [,π. Choose an arbitrary s C U, compute r γ, where r = σ i Pγ,Gs ))). ) Iteration k=,,...) a) Given r k, γ k, s k, S k. b) Reset Θ k =. c) For θ S k θ can be discretized steps of S k ) Compute Rk θ by finding all real nonnegative s λ H r k,â, ˆB,Ĉ, ˆD, ˆ,T γ k,n+rp nr, )) P nm T γk,n+m such that σ i P γk,g se iθ))) = r k. If Rk θ, then update Θ k Θ k θ. d) Update S k = S k Θ k. e) Quit if the length of the largest interval of S k < TOL Sk sizer k ) < TOL Rk. f) Find one point s k Rk θ such that r k < r k, where θ S k, r k = τ i Gs k )). ) Quit if r k < TOL rk. ) Update k k + goto step. E. Algorithm for solving problem ) To solve problem ), only two slight modifications of Algorithm. are needed, which actually results in a simpler algorithm. Firstly, instead of computing the eigenvalues ) of H r k,â, ˆB,Ĉ, ˆD, ˆ,T γ k,n+rp nr,p nm T γk,n+m in step -c) to obtain Rk θ, one is only required to compute the eigenvalues of H ) r k,a,b,c,d,e iθ I n,i n+r,i n+m the proof is trivial follows directly from Theorem.). Furthermore, when solving problem ), all the values within R k actually achieve a smaller radius than r k all the values outside R k do not). This is not true when solving problem ), where R k only achieves a bound. Hence, step -f) for solving problem ) can very easily be accomplished by selecting any point in R k.
5 7th IEEE CDC, Cancun, Mexico, Dec. 9-, 8 F. Algorithm for computing the real stabilizability radius To compute the real stabilizability radius, the real unstable DFM radius, the real minimum-phase radius, where the search is for a minimum in the closed right half of the complex plane i.e. by solving the following problem r RA,B,C,D) i A si B = min τ i ) s C + for given A,B,C,D) i {,...,n + minr,m)} only one very simply modification of Algorithm. is needed. In particular, one only needs to change the initialization step ) of Algorithm. from S = [,π to S = [, π the proof is trivial). G. Operations count To provide an idea of the computational requirements of using Algorithm., the number of operations required in terms of the number of singular value σ), real perturbation value τ), eigenvalue λ) problems are noted in Table I. It is to be noted that every time the radius is evaluated at a particular s C, a real perturbation value problem is solved. This occurs in step ) -f). In step ), the real perturbation value problem is evaluated only once, but step -f) may require multiple evaluations, depending on the number of trials needed before r k < r k is obtained. However, it is found from experiment that the number of trials is typically small i.e. close to when using the method outlined in Section IV-B). It is also to be noted that the main portion of the total operations count is from computing Rk θ for all θ S k in step -c). In terms of implementation, S k is discretized into a number of points. From experiment, it is found that it is generally sufficient to discretize S k into to points, or a minimum of. i.e. a maximum of 6 = π. ) points), whichever results in a smaller angular division. For a given θ S k, computing Rk θ requires: i) solving for the eigenvalues, Λ, of the n n matrix ) H r k,â, ˆB,Ĉ, ˆD, ˆ,T γ k,n+rp nr,p nm T γk,n+m ; ii) evaluating a number of n + r) n + m)-sized singular value problems to verify which real nonnegative eigenvalues, s Λ, there are a maximum of n such values) have the property that they satisfy 8). Furthermore, since there are at most n real nonnegative eigenvalues s Λ that may satisfy 8), Rk θ contains at most n + intervals. Hence, to determine which intervals in Rk θ achieve a radius less than r k, at most n + points are tested, resulting in n + additional n + r) n + m)-sized singular value problems. Remark.: It should be noted that Algorithm. solves only a limited number of real perturbation value problems, which by itself is a -D optimization problem involving multiple n+r) n+m)-sized singular value problems. Hence, this proves to be an advantage of Algorithm., as compared to say, a gradient search method, which requires solving a real perturbation value problem at each trial point, at each step involved with computing the steepest gradient. TABLE I SUMMARY OF OPERATIONS COUNT OF ALGORITHM. Step # of σ,τ,λ) Size of matrix τ n + r) n + m) -c for each λ n n θ S k ) at most n + ) σ n + r) n + m) -f at least τ n + r) n + m) Here, we denote σ, τ, λ as singular value, real perturbation value, eigenvalue problems. V. NUMERICAL EXAMPLE The following example can be found in [, where the real controllability radius of A =. B =. is found to be.9, is achieved at s =.97+j.98. Using Algorithm. with the arbitrary) starting point s = j, the same radius global minimizer are found to within significant figures in iterations to 6 significant figures in 7 iterations. Table II lists r k, s k, the minimizing set, S k for each iteration k =,...,9. Figure provides a grid plot of the real controllability radius with respect to a given s C. Figures to plot R k outlined by the dots) obtained at Iterations,,, 8 respectively, superimposed on a contour plot of Figure. The straight lines originating from the origin depict S k. It is interesting to note that on a computer with a Pentium IV.GHz processor, MB of RAM, MATLAB 7., Algorithm. took a total of about 6 sec. to complete, which is approximately the same amount of time a gradient search method requires to obtain a possibly local) minimum from a single starting point. Hence, if the gradient search method has to test initial points in order to achieve a certain level of confidence that the obtained minimum is global, then the gradient search method will take approximately times longer to run than Algorithm.. VI. CONCLUSIONS In this paper, an efficient algorithm is presented for solving the general -D optimization problems ) ), which are essential for computing the complex TABLE II ESTIMATES OF THE GLOBAL MINIMUM RADIUS r k ), THE MINIMIZER s k ), AND THE MINIMIZING SET S k ) AT EACH ITERATION k TuA6. Iter. k) r k s k S k.767 j [., j. [., [., j.86 [., j.97 [., j.98 [.7, j.9879 [., j.98 [.7, j.989 [.9, j.9897 [.,.
6 7th IEEE CDC, Cancun, Mexico, Dec. 9-, 8 TuA6. Grid plot of real controllability radius with respect to s Iteration radius Fig.. Grid plot of the controllability radius with respect to a given s C Iteration..... Fig.. Iteration Iteration Fig.. Iteration real controllability radius, the minimum-phase radius, the DFM radius, etc. The algorithm works by iteratively reducing the size of the so-called minimizing set which contains the global minimizer. Unlike methods such as the gradient search method, the choice of the initial point of the proposed Fig.. Iteration Iteration 8 Fig.. Iteration 8 algorithm is not crucial, the minimum is obtained with a high degree of confidence that it is indeed global. REFERENCES [ R. Eising, Between controllable uncontrollable, Systems Control Letters, vol., pp. 6 6, 98. [ G. Hu E. J. Davison, Real controllability/stabilizability radius of LTI systems, IEEE Trans. Automat. Contr., vol. 9, no., pp. 7,. [ L. Qiu, B. Bernhardsson, A. Rantzer, E. J. Davison, P. M. Young, J. C. Doyle, A formula for computation of the real stability radius, Automatica, vol., no. 6, pp , 99. [ A. F. Vaz E. J. Davison, A measure for the decentralized assignability of eigenvalues, Systems & Control Letters, vol., pp. 9 99, 988. [ S. Lam E. J. Davison, The real decentralized fixed mode radius of LTI systems, in 6 th IEEE Conference on Decision Control, New Orleans, LA, USA, December - 7, pp. 6. [6, The transmission zero at s radius the minimum phase radius of LTI systems, in 7 th International Federation of Automatic Control IFAC) world congress, Seoul, Korea, July 6-8, pp [7 J. Sreedhar, P. V. Dooren, A. L. Tits, A fast algorithm to compute the real structured stability radius, in International Series of Numerical Mathematics. Birkhäuser Verlag Basel, 996, vol., pp. 9. [8 B. Bernhardsson, A. Rantzer, L. Qiu, Real perturbation values real quadratic forms in a complex vector space, Linear Algebra its Applications, vol. 7, pp., 998.
An efficient algorithm to compute the real perturbation values of a matrix
An efficient algorithm to compute the real perturbation values of a matrix Simon Lam and Edward J. Davison 1 Abstract In this paper, an efficient algorithm is presented for solving the nonlinear 1-D optimization
More informationThe Transmission Zero at s Radius and the Minimum Phase Radius of LTI Systems
Proceedings of the 17th World ongress The International ederation of Automatic ontrol The Transmission Zero at s adius and the Minimum Phase adius of LTI Systems Simon Lam and Edward J. Davison Systems
More informationReal Robustness Radii and Performance Limitations of LTI Control Systems. Simon Sai-Ming Lam
Real Robustness Radii and Performance Limitations of LTI Control Systems by Simon Sai-Ming Lam A thesis submitted in conformity with the requirements for the degree of Doctor of Philosophy Graduate Department
More informationStability radii and spectral value sets for real matrix perturbations Introduction D. Hinrichsen and B. Kelb Institut fur Dynamische Systeme Universitat Bremen D-859 Bremen, F. R. G. dh@mathematik.uni-bremen.de
More informationRecent Advances in Positive Systems: The Servomechanism Problem
Recent Advances in Positive Systems: The Servomechanism Problem 47 th IEEE Conference on Decision and Control December 28. Bartek Roszak and Edward J. Davison Systems Control Group, University of Toronto
More informationSUCCESSIVE POLE SHIFTING USING SAMPLED-DATA LQ REGULATORS. Sigeru Omatu
SUCCESSIVE POLE SHIFING USING SAMPLED-DAA LQ REGULAORS oru Fujinaka Sigeru Omatu Graduate School of Engineering, Osaka Prefecture University, 1-1 Gakuen-cho, Sakai, 599-8531 Japan Abstract: Design of sampled-data
More informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product
More informationIMPULSIVE CONTROL OF DISCRETE-TIME NETWORKED SYSTEMS WITH COMMUNICATION DELAYS. Shumei Mu, Tianguang Chu, and Long Wang
IMPULSIVE CONTROL OF DISCRETE-TIME NETWORKED SYSTEMS WITH COMMUNICATION DELAYS Shumei Mu Tianguang Chu and Long Wang Intelligent Control Laboratory Center for Systems and Control Department of Mechanics
More informationPreliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012
Instructions Preliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012 The exam consists of four problems, each having multiple parts. You should attempt to solve all four problems. 1.
More information22.3. Repeated Eigenvalues and Symmetric Matrices. Introduction. Prerequisites. Learning Outcomes
Repeated Eigenvalues and Symmetric Matrices. Introduction In this Section we further develop the theory of eigenvalues and eigenvectors in two distinct directions. Firstly we look at matrices where one
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 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 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 informationZeros and zero dynamics
CHAPTER 4 Zeros and zero dynamics 41 Zero dynamics for SISO systems Consider a linear system defined by a strictly proper scalar transfer function that does not have any common zero and pole: g(s) =α p(s)
More informationarxiv: v1 [cs.sy] 29 Dec 2018
ON CHECKING NULL RANK CONDITIONS OF RATIONAL MATRICES ANDREAS VARGA Abstract. In this paper we discuss possible numerical approaches to reliably check the rank condition rankg(λ) = 0 for a given rational
More informationModel reduction via tangential interpolation
Model reduction via tangential interpolation K. Gallivan, A. Vandendorpe and P. Van Dooren May 14, 2002 1 Introduction Although most of the theory presented in this paper holds for both continuous-time
More informationẋ n = f n (x 1,...,x n,u 1,...,u m ) (5) y 1 = g 1 (x 1,...,x n,u 1,...,u m ) (6) y p = g p (x 1,...,x n,u 1,...,u m ) (7)
EEE582 Topical Outline A.A. Rodriguez Fall 2007 GWC 352, 965-3712 The following represents a detailed topical outline of the course. It attempts to highlight most of the key concepts to be covered and
More information1 Continuous-time Systems
Observability Completely controllable systems can be restructured by means of state feedback to have many desirable properties. But what if the state is not available for feedback? What if only the output
More informationDecentralized control with input saturation
Decentralized control with input saturation Ciprian Deliu Faculty of Mathematics and Computer Science Technical University Eindhoven Eindhoven, The Netherlands November 2006 Decentralized control with
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 informationIntroduction to Matrices
POLS 704 Introduction to Matrices Introduction to Matrices. The Cast of Characters A matrix is a rectangular array (i.e., a table) of numbers. For example, 2 3 X 4 5 6 (4 3) 7 8 9 0 0 0 Thismatrix,with4rowsand3columns,isoforder
More informationComputing generalized inverse systems using matrix pencil methods
Computing generalized inverse systems using matrix pencil methods A. Varga German Aerospace Center DLR - Oberpfaffenhofen Institute of Robotics and Mechatronics D-82234 Wessling, Germany. Andras.Varga@dlr.de
More informationOn Identification of Cascade Systems 1
On Identification of Cascade Systems 1 Bo Wahlberg Håkan Hjalmarsson Jonas Mårtensson Automatic Control and ACCESS, School of Electrical Engineering, KTH, SE-100 44 Stockholm, Sweden. (bo.wahlberg@ee.kth.se
More informationRepeated Eigenvalues and Symmetric Matrices
Repeated Eigenvalues and Symmetric Matrices. Introduction In this Section we further develop the theory of eigenvalues and eigenvectors in two distinct directions. Firstly we look at matrices where one
More informationGraph and Controller Design for Disturbance Attenuation in Consensus Networks
203 3th International Conference on Control, Automation and Systems (ICCAS 203) Oct. 20-23, 203 in Kimdaejung Convention Center, Gwangju, Korea Graph and Controller Design for Disturbance Attenuation in
More informationLinear algebra issues in Interior Point methods for bound-constrained least-squares problems
Linear algebra issues in Interior Point methods for bound-constrained least-squares problems Stefania Bellavia Dipartimento di Energetica S. Stecco Università degli Studi di Firenze Joint work with Jacek
More informationPerturbation theory for eigenvalues of Hermitian pencils. Christian Mehl Institut für Mathematik TU Berlin, Germany. 9th Elgersburg Workshop
Perturbation theory for eigenvalues of Hermitian pencils Christian Mehl Institut für Mathematik TU Berlin, Germany 9th Elgersburg Workshop Elgersburg, March 3, 2014 joint work with Shreemayee Bora, Michael
More informationMath Camp Lecture 4: Linear Algebra. Xiao Yu Wang. Aug 2010 MIT. Xiao Yu Wang (MIT) Math Camp /10 1 / 88
Math Camp 2010 Lecture 4: Linear Algebra Xiao Yu Wang MIT Aug 2010 Xiao Yu Wang (MIT) Math Camp 2010 08/10 1 / 88 Linear Algebra Game Plan Vector Spaces Linear Transformations and Matrices Determinant
More information[k,g,gfin] = hinfsyn(p,nmeas,ncon,gmin,gmax,tol)
8 H Controller Synthesis: hinfsyn The command syntax is k,g,gfin = hinfsyn(p,nmeas,ncon,gmin,gmax,tol) associated with the general feedback diagram below. e y P d u e G = F L (P, K) d K hinfsyn calculates
More informationA Cross-Associative Neural Network for SVD of Nonsquared Data Matrix in Signal Processing
IEEE TRANSACTIONS ON NEURAL NETWORKS, VOL. 12, NO. 5, SEPTEMBER 2001 1215 A Cross-Associative Neural Network for SVD of Nonsquared Data Matrix in Signal Processing Da-Zheng Feng, Zheng Bao, Xian-Da Zhang
More informationMulti-Model Adaptive Regulation for a Family of Systems Containing Different Zero Structures
Preprints of the 19th World Congress The International Federation of Automatic Control Multi-Model Adaptive Regulation for a Family of Systems Containing Different Zero Structures Eric Peterson Harry G.
More informationThe ϵ-capacity of a gain matrix and tolerable disturbances: Discrete-time perturbed linear systems
IOSR Journal of Mathematics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 11, Issue 3 Ver. IV (May - Jun. 2015), PP 52-62 www.iosrjournals.org The ϵ-capacity of a gain matrix and tolerable disturbances:
More informationComputational Methods CMSC/AMSC/MAPL 460. Eigenvalues and Eigenvectors. Ramani Duraiswami, Dept. of Computer Science
Computational Methods CMSC/AMSC/MAPL 460 Eigenvalues and Eigenvectors Ramani Duraiswami, Dept. of Computer Science Eigen Values of a Matrix Recap: A N N matrix A has an eigenvector x (non-zero) with corresponding
More informationThe Solvability Conditions for the Inverse Eigenvalue Problem of Hermitian and Generalized Skew-Hamiltonian Matrices and Its Approximation
The Solvability Conditions for the Inverse Eigenvalue Problem of Hermitian and Generalized Skew-Hamiltonian Matrices and Its Approximation Zheng-jian Bai Abstract In this paper, we first consider the inverse
More 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 informationThe Kalman-Yakubovich-Popov Lemma for Differential-Algebraic Equations with Applications
MAX PLANCK INSTITUTE Elgersburg Workshop Elgersburg February 11-14, 2013 The Kalman-Yakubovich-Popov Lemma for Differential-Algebraic Equations with Applications Timo Reis 1 Matthias Voigt 2 1 Department
More informationAPPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.
APPENDIX A Background Mathematics A. Linear Algebra A.. Vector algebra Let x denote the n-dimensional column vector with components 0 x x 2 B C @. A x n Definition 6 (scalar product). The scalar product
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 informationIterative Methods for Solving A x = b
Iterative Methods for Solving A x = b A good (free) online source for iterative methods for solving A x = b is given in the description of a set of iterative solvers called templates found at netlib: http
More informationELA THE MINIMUM-NORM LEAST-SQUARES SOLUTION OF A LINEAR SYSTEM AND SYMMETRIC RANK-ONE UPDATES
Volume 22, pp. 480-489, May 20 THE MINIMUM-NORM LEAST-SQUARES SOLUTION OF A LINEAR SYSTEM AND SYMMETRIC RANK-ONE UPDATES XUZHOU CHEN AND JUN JI Abstract. In this paper, we study the Moore-Penrose inverse
More informationModel reduction for linear systems by balancing
Model reduction for linear systems by balancing Bart Besselink Jan C. Willems Center for Systems and Control Johann Bernoulli Institute for Mathematics and Computer Science University of Groningen, Groningen,
More informationDiscrete Simulation of Power Law Noise
Discrete Simulation of Power Law Noise Neil Ashby 1,2 1 University of Colorado, Boulder, CO 80309-0390 USA 2 National Institute of Standards and Technology, Boulder, CO 80305 USA ashby@boulder.nist.gov
More informationand the nite horizon cost index with the nite terminal weighting matrix F > : N?1 X J(z r ; u; w) = [z(n)? z r (N)] T F [z(n)? z r (N)] + t= [kz? z r
Intervalwise Receding Horizon H 1 -Tracking Control for Discrete Linear Periodic Systems Ki Baek Kim, Jae-Won Lee, Young Il. Lee, and Wook Hyun Kwon School of Electrical Engineering Seoul National University,
More informationRECURSIVE SUBSPACE IDENTIFICATION IN THE LEAST SQUARES FRAMEWORK
RECURSIVE SUBSPACE IDENTIFICATION IN THE LEAST SQUARES FRAMEWORK TRNKA PAVEL AND HAVLENA VLADIMÍR Dept of Control Engineering, Czech Technical University, Technická 2, 166 27 Praha, Czech Republic mail:
More informationDECENTRALIZED CONTROL DESIGN USING LMI MODEL REDUCTION
Journal of ELECTRICAL ENGINEERING, VOL. 58, NO. 6, 2007, 307 312 DECENTRALIZED CONTROL DESIGN USING LMI MODEL REDUCTION Szabolcs Dorák Danica Rosinová Decentralized control design approach based on partial
More informationFast Computation of Moore-Penrose Inverse Matrices
Fast Computation of Moore-Penrose Inverse Matrices Pierre Courrieu To cite this version: Pierre Courrieu. Fast Computation of Moore-Penrose Inverse Matrices. Neural Information Processing - Letters and
More informationPOLE PLACEMENT. Sadegh Bolouki. Lecture slides for ECE 515. University of Illinois, Urbana-Champaign. Fall S. Bolouki (UIUC) 1 / 19
POLE PLACEMENT Sadegh Bolouki Lecture slides for ECE 515 University of Illinois, Urbana-Champaign Fall 2016 S. Bolouki (UIUC) 1 / 19 Outline 1 State Feedback 2 Observer 3 Observer Feedback 4 Reduced Order
More informationLecture Note 7: Iterative methods for solving linear systems. Xiaoqun Zhang Shanghai Jiao Tong University
Lecture Note 7: Iterative methods for solving linear systems Xiaoqun Zhang Shanghai Jiao Tong University Last updated: December 24, 2014 1.1 Review on linear algebra Norms of vectors and matrices vector
More informationRepresent this system in terms of a block diagram consisting only of. g From Newton s law: 2 : θ sin θ 9 θ ` T
Exercise (Block diagram decomposition). Consider a system P that maps each input to the solutions of 9 4 ` 3 9 Represent this system in terms of a block diagram consisting only of integrator systems, represented
More informationA Multi-Step Hybrid Method for Multi-Input Partial Quadratic Eigenvalue Assignment with Time Delay
A Multi-Step Hybrid Method for Multi-Input Partial Quadratic Eigenvalue Assignment with Time Delay Zheng-Jian Bai Mei-Xiang Chen Jin-Ku Yang April 14, 2012 Abstract A hybrid method was given by Ram, Mottershead,
More informationCOMP 558 lecture 18 Nov. 15, 2010
Least squares We have seen several least squares problems thus far, and we will see more in the upcoming lectures. For this reason it is good to have a more general picture of these problems and how to
More informationNetwork Reconstruction from Intrinsic Noise: Non-Minimum-Phase Systems
Preprints of the 19th World Congress he International Federation of Automatic Control Network Reconstruction from Intrinsic Noise: Non-Minimum-Phase Systems David Hayden, Ye Yuan Jorge Goncalves Department
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 informationA = 3 B = A 1 1 matrix is the same as a number or scalar, 3 = [3].
Appendix : A Very Brief Linear ALgebra Review Introduction Linear Algebra, also known as matrix theory, is an important element of all branches of mathematics Very often in this course we study the shapes
More informationMCE/EEC 647/747: Robot Dynamics and Control. Lecture 8: Basic Lyapunov Stability Theory
MCE/EEC 647/747: Robot Dynamics and Control Lecture 8: Basic Lyapunov Stability Theory Reading: SHV Appendix Mechanical Engineering Hanz Richter, PhD MCE503 p.1/17 Stability in the sense of Lyapunov A
More informationStructured singular values for Bernoulli matrices
Int. J. Adv. Appl. Math. and Mech. 44 2017 41 49 ISSN: 2347-2529 IJAAMM Journal homepage: www.ijaamm.com International Journal of Advances in Applied Mathematics and Mechanics Structured singular values
More informationNumerical Methods - Numerical Linear Algebra
Numerical Methods - Numerical Linear Algebra Y. K. Goh Universiti Tunku Abdul Rahman 2013 Y. K. Goh (UTAR) Numerical Methods - Numerical Linear Algebra I 2013 1 / 62 Outline 1 Motivation 2 Solving Linear
More informationA VERY BRIEF LINEAR ALGEBRA REVIEW for MAP 5485 Introduction to Mathematical Biophysics Fall 2010
A VERY BRIEF LINEAR ALGEBRA REVIEW for MAP 5485 Introduction to Mathematical Biophysics Fall 00 Introduction Linear Algebra, also known as matrix theory, is an important element of all branches of mathematics
More informationDistance bounds for prescribed multiple eigenvalues of matrix polynomials
Distance bounds for prescribed multiple eigenvalues of matrix polynomials Panayiotis J Psarrakos February 5, Abstract In this paper, motivated by a problem posed by Wilkinson, we study the coefficient
More informationUNCERTAINTY MODELING VIA FREQUENCY DOMAIN MODEL VALIDATION
AIAA 99-3959 UNCERTAINTY MODELING VIA FREQUENCY DOMAIN MODEL VALIDATION Martin R. Waszak, * NASA Langley Research Center, Hampton, Virginia Dominick Andrisani II, Purdue University, West Lafayette, Indiana
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 informationStructural Properties of LTI Singular Systems under Output Feedback
Structural Properties of LTI Singular Systems under Output Feedback Runyi Yu Dept. of Electrical and Electronic Eng. Eastern Mediterranean University Gazimagusa, Mersin 10, Turkey runyi.yu@emu.edu.tr Dianhui
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 23: GMRES and Other Krylov Subspace Methods; Preconditioning
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 23: GMRES and Other Krylov Subspace Methods; Preconditioning Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 18 Outline
More informationState Feedback and State Estimators Linear System Theory and Design, Chapter 8.
1 Linear System Theory and Design, http://zitompul.wordpress.com 2 0 1 4 2 Homework 7: State Estimators (a) For the same system as discussed in previous slides, design another closed-loop state estimator,
More informationRobust Stabilization of the Uncertain Linear Systems. Based on Descriptor Form Representation t. Toru ASAI* and Shinji HARA**
Robust Stabilization of the Uncertain Linear Systems Based on Descriptor Form Representation t Toru ASAI* and Shinji HARA** This paper proposes a necessary and sufficient condition for the quadratic stabilization
More informationLinear Quadratic Zero-Sum Two-Person Differential Games Pierre Bernhard June 15, 2013
Linear Quadratic Zero-Sum Two-Person Differential Games Pierre Bernhard June 15, 2013 Abstract As in optimal control theory, linear quadratic (LQ) differential games (DG) can be solved, even in high dimension,
More informationOn the Stabilization of Neutrally Stable Linear Discrete Time Systems
TWCCC Texas Wisconsin California Control Consortium Technical report number 2017 01 On the Stabilization of Neutrally Stable Linear Discrete Time Systems Travis J. Arnold and James B. Rawlings Department
More informationPOSITIVE SEMIDEFINITE INTERVALS FOR MATRIX PENCILS
POSITIVE SEMIDEFINITE INTERVALS FOR MATRIX PENCILS RICHARD J. CARON, HUIMING SONG, AND TIM TRAYNOR Abstract. Let A and E be real symmetric matrices. In this paper we are concerned with the determination
More informationA Novel Integral-Based Event Triggering Control for Linear Time-Invariant Systems
53rd IEEE Conference on Decision and Control December 15-17, 2014. Los Angeles, California, USA A Novel Integral-Based Event Triggering Control for Linear Time-Invariant Systems Seyed Hossein Mousavi 1,
More informationChapter 5. Linear Algebra. A linear (algebraic) equation in. unknowns, x 1, x 2,..., x n, is. an equation of the form
Chapter 5. Linear Algebra A linear (algebraic) equation in n unknowns, x 1, x 2,..., x n, is an equation of the form a 1 x 1 + a 2 x 2 + + a n x n = b where a 1, a 2,..., a n and b are real numbers. 1
More informationA Globally Stabilizing Receding Horizon Controller for Neutrally Stable Linear Systems with Input Constraints 1
A Globally Stabilizing Receding Horizon Controller for Neutrally Stable Linear Systems with Input Constraints 1 Ali Jadbabaie, Claudio De Persis, and Tae-Woong Yoon 2 Department of Electrical Engineering
More informationReal Eigenvalue Extraction and the Distance to Uncontrollability
Real Eigenvalue Extraction and the Distance to Uncontrollability Emre Mengi Computer Science Department Courant Institute of Mathematical Sciences New York University mengi@cs.nyu.edu May 22nd, 2006 Emre
More informationLOW ORDER H CONTROLLER DESIGN: AN LMI APPROACH
LOW ORDER H CONROLLER DESIGN: AN LMI APPROACH Guisheng Zhai, Shinichi Murao, Naoki Koyama, Masaharu Yoshida Faculty of Systems Engineering, Wakayama University, Wakayama 640-8510, Japan Email: zhai@sys.wakayama-u.ac.jp
More informationISOLATED SEMIDEFINITE SOLUTIONS OF THE CONTINUOUS-TIME ALGEBRAIC RICCATI EQUATION
ISOLATED SEMIDEFINITE SOLUTIONS OF THE CONTINUOUS-TIME ALGEBRAIC RICCATI EQUATION Harald K. Wimmer 1 The set of all negative-semidefinite solutions of the CARE A X + XA + XBB X C C = 0 is homeomorphic
More informationState estimation of uncertain multiple model with unknown inputs
State estimation of uncertain multiple model with unknown inputs Abdelkader Akhenak, Mohammed Chadli, Didier Maquin and José Ragot Centre de Recherche en Automatique de Nancy, CNRS UMR 79 Institut National
More informationFoundations of Matrix Analysis
1 Foundations of Matrix Analysis In this chapter we recall the basic elements of linear algebra which will be employed in the remainder of the text For most of the proofs as well as for the details, the
More informationResearch Article Partial Pole Placement in LMI Region
Control Science and Engineering Article ID 84128 5 pages http://dxdoiorg/11155/214/84128 Research Article Partial Pole Placement in LMI Region Liuli Ou 1 Shaobo Han 2 Yongji Wang 1 Shuai Dong 1 and Lei
More informationProjection of state space realizations
Chapter 1 Projection of state space realizations Antoine Vandendorpe and Paul Van Dooren Department of Mathematical Engineering Université catholique de Louvain B-1348 Louvain-la-Neuve Belgium 1.0.1 Description
More informationApplications of Controlled Invariance to the l 1 Optimal Control Problem
Applications of Controlled Invariance to the l 1 Optimal Control Problem Carlos E.T. Dórea and Jean-Claude Hennet LAAS-CNRS 7, Ave. du Colonel Roche, 31077 Toulouse Cédex 4, FRANCE Phone : (+33) 61 33
More informationEnforcing Passivity for Admittance Matrices Approximated by Rational Functions
IEEE TRANSACTIONS ON POWER SYSTEMS, VOL. 16, NO. 1, FEBRUARY 2001 97 Enforcing Passivity for Admittance Matrices Approximated by Rational Functions Bjørn Gustavsen, Member, IEEE and Adam Semlyen, Life
More informationEigenvalues, Eigenvectors. Eigenvalues and eigenvector will be fundamentally related to the nature of the solutions of state space systems.
Chapter 3 Linear Algebra In this Chapter we provide a review of some basic concepts from Linear Algebra which will be required in order to compute solutions of LTI systems in state space form, discuss
More informationLinear System Theory
Linear System Theory Wonhee Kim Chapter 6: Controllability & Observability Chapter 7: Minimal Realizations May 2, 217 1 / 31 Recap State space equation Linear Algebra Solutions of LTI and LTV system Stability
More informationBALANCING-RELATED MODEL REDUCTION FOR DATA-SPARSE SYSTEMS
BALANCING-RELATED Peter Benner Professur Mathematik in Industrie und Technik Fakultät für Mathematik Technische Universität Chemnitz Computational Methods with Applications Harrachov, 19 25 August 2007
More informationFirst, we review some important facts on the location of eigenvalues of matrices.
BLOCK NORMAL MATRICES AND GERSHGORIN-TYPE DISCS JAKUB KIERZKOWSKI AND ALICJA SMOKTUNOWICZ Abstract The block analogues of the theorems on inclusion regions for the eigenvalues of normal matrices are given
More informationSTAT 309: MATHEMATICAL COMPUTATIONS I FALL 2017 LECTURE 5
STAT 39: MATHEMATICAL COMPUTATIONS I FALL 17 LECTURE 5 1 existence of svd Theorem 1 (Existence of SVD) Every matrix has a singular value decomposition (condensed version) Proof Let A C m n and for simplicity
More informationVISUALIZING PSEUDOSPECTRA FOR POLYNOMIAL EIGENVALUE PROBLEMS. Adéla Klimentová *, Michael Šebek ** Czech Technical University in Prague
VSUALZNG PSEUDOSPECTRA FOR POLYNOMAL EGENVALUE PROBLEMS Adéla Klimentová *, Michael Šebek ** * Department of Control Engineering Czech Technical University in Prague ** nstitute of nformation Theory and
More informationSimultaneous State and Fault Estimation for Descriptor Systems using an Augmented PD Observer
Preprints of the 19th World Congress The International Federation of Automatic Control Simultaneous State and Fault Estimation for Descriptor Systems using an Augmented PD Observer Fengming Shi*, Ron J.
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 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 informationA NONSMOOTH, NONCONVEX OPTIMIZATION APPROACH TO ROBUST STABILIZATION BY STATIC OUTPUT FEEDBACK AND LOW-ORDER CONTROLLERS
A NONSMOOTH, NONCONVEX OPTIMIZATION APPROACH TO ROBUST STABILIZATION BY STATIC OUTPUT FEEDBACK AND LOW-ORDER CONTROLLERS James V Burke,1 Adrian S Lewis,2 Michael L Overton,3 University of Washington, Seattle,
More informationApplied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots Programming in
More informationScientific Computing: Dense Linear Systems
Scientific Computing: Dense Linear Systems Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course MATH-GA.2043 or CSCI-GA.2112, Spring 2012 February 9th, 2012 A. Donev (Courant Institute)
More information1 Cricket chirps: an example
Notes for 2016-09-26 1 Cricket chirps: an example Did you know that you can estimate the temperature by listening to the rate of chirps? The data set in Table 1 1. represents measurements of the number
More informationNonlinear Control Design for Linear Differential Inclusions via Convex Hull Quadratic Lyapunov Functions
Nonlinear Control Design for Linear Differential Inclusions via Convex Hull Quadratic Lyapunov Functions Tingshu Hu Abstract This paper presents a nonlinear control design method for robust stabilization
More informationChap 3. Linear Algebra
Chap 3. Linear Algebra Outlines 1. Introduction 2. Basis, Representation, and Orthonormalization 3. Linear Algebraic Equations 4. Similarity Transformation 5. Diagonal Form and Jordan Form 6. Functions
More informationLinear Quadratic Zero-Sum Two-Person Differential Games
Linear Quadratic Zero-Sum Two-Person Differential Games Pierre Bernhard To cite this version: Pierre Bernhard. Linear Quadratic Zero-Sum Two-Person Differential Games. Encyclopaedia of Systems and Control,
More informationyou expect to encounter difficulties when trying to solve A x = b? 4. A composite quadrature rule has error associated with it in the following form
Qualifying exam for numerical analysis (Spring 2017) Show your work for full credit. If you are unable to solve some part, attempt the subsequent parts. 1. Consider the following finite difference: f (0)
More information2nd Symposium on System, Structure and Control, Oaxaca, 2004
263 2nd Symposium on System, Structure and Control, Oaxaca, 2004 A PROJECTIVE ALGORITHM FOR STATIC OUTPUT FEEDBACK STABILIZATION Kaiyang Yang, Robert Orsi and John B. Moore Department of Systems Engineering,
More informationUnconstrained optimization
Chapter 4 Unconstrained optimization An unconstrained optimization problem takes the form min x Rnf(x) (4.1) for a target functional (also called objective function) f : R n R. In this chapter and throughout
More informationZero controllability in discrete-time structured systems
1 Zero controllability in discrete-time structured systems Jacob van der Woude arxiv:173.8394v1 [math.oc] 24 Mar 217 Abstract In this paper we consider complex dynamical networks modeled by means of state
More information