Semi-definite Programming Techniques for Quadratic Inverse Eigenvalue Problem
|
|
- Louisa Patrick
- 5 years ago
- Views:
Transcription
1 Semi-definite Programming Techniques for Quadratic Inverse Eigenvalue Problem Moody T. Chu North Carolina State University March 31, National Chiao Tung University
2 Outline Basic Ideas General Quadratic Model Constraints and Fundamental Questions Some Applications Brief History Model Updating Spectrum Modification Challenges Structure Preservation Symmetry Positive Semi-definiteness Other Difficulties Semi-definite Programming Quick Overview Solvability Test Model Updating Conclusion
3 Outline Basic Ideas General Quadratic Model Constraints and Fundamental Questions Some Applications Brief History Model Updating Spectrum Modification Challenges Structure Preservation Symmetry Positive Semi-definiteness Other Difficulties Semi-definite Programming Quick Overview Solvability Test Model Updating Conclusion
4 Outline Basic Ideas General Quadratic Model Constraints and Fundamental Questions Some Applications Brief History Model Updating Spectrum Modification Challenges Structure Preservation Symmetry Positive Semi-definiteness Other Difficulties Semi-definite Programming Quick Overview Solvability Test Model Updating Conclusion
5 Outline Basic Ideas General Quadratic Model Constraints and Fundamental Questions Some Applications Brief History Model Updating Spectrum Modification Challenges Structure Preservation Symmetry Positive Semi-definiteness Other Difficulties Semi-definite Programming Quick Overview Solvability Test Model Updating Conclusion
6 Outline Basic Ideas General Quadratic Model Constraints and Fundamental Questions Some Applications Brief History Model Updating Spectrum Modification Challenges Structure Preservation Symmetry Positive Semi-definiteness Other Difficulties Semi-definite Programming Quick Overview Solvability Test Model Updating Conclusion
7 Mass-Spring Vibration Equation of motion: m m m m ẍ+ c 1 + c 2 0 c c 2 0 c 2 + c 3 c c 3 c ẋ+ k 1 + k 2 + k 5 k 2 k 5 0 k 2 k 2 + k 3 k 3 0 k 5 k 3 k 3 + k 4 + k 5 k k 4 k x = f(t).
8 Second-Order Differential System Dynamical system with n-degree-of-freedom: Mẍ + (C + G)ẋ + (K + N)x = F. Some interpretations: M := Mass matrix M = M 0. C := Damping matrix C = C. K := Stiffness matrix K = K 0. G := Gyroscopic matrix G = G. N := Circulatory matrix N = N. F := External force.
9 Integrated Circuit Ohm s Law: V = RI = L di dt = 1 C I dt. Kirchhoff s Law: The algebraic sum of the current at any junction must be zero. The algebraic sum of voltage drops around any closed loop must be zero. Very large scale integration.
10 Quadratic Eigenvalue Problem Assume the homogeneous solution x(t): x = e λt u. Look for nontrivial solution to the QEP: Q(λ)u := (λ 2 M + λc + K )u = 0. If M is nonsingular, then there are 2n eigenpairs. Many applications. Many numerical techniques.
11 Direct Problem versus Inverse Problem Direct Problem: Given physical parameters, e.g., mass, length, elasticity, etc., analyze and derive the spectral information and, hence, induce the dynamical behavior. Inverse Problem: Validate, determine, or estimate the parameters of the system according to its observed or expected behavior.
12 A Simplified Interpretation Direct problem Express the behavior in terms of the parameters. Inverse problem Express the parameters in term of the behavior.
13 Structured Problem Be it from finite element discretization or from real physical configuration, the coefficient matrices M, C and K are often structured. The structure of symmetry and positive definiteness in the mass-spring system. The structure of asymmetry and connectivity in the RCL model. M = C = K = L 2 L L 2 L L L , 0 R 2 R 2 0 R 1 + R R 4 0 R 2 R 2 + R 3 0 R R C C3 C C C3 C , How to solve the structured problem?
14 Partial Spectral Information In a large and complicated physical system, it is often impossible to obtain the entire spectral information. Quantities related to high frequency terms generally are susceptible to measurement errors due to the finite bandwidth of measuring devices. Only partial spectral information are of interest and observable. How much eigeninformation is enough?
15 Solvability: Fundamental Questions Determine a necessary or a sufficient condition under which an IEP has a solution. Computability: Develop a scheme through which, knowing a priori that the given spectral data are feasible, the solution to an IEP can be constructed numerically. Sensitivity: Quantify how a solution to an IEP is subject to changes of the spectral data. Applicability: Distinguish whether the given data are exact or approximate, complete or incomplete, and whether only an estimation of the parameters of the system is sufficient. Decide between physical realizability and physical uncertainty which constraint of the problem should be enforced.
16 Vibrating String m 2 m 3 F m 1 x 1 x 2 x 3 m 4 x 4 F h Discretization: Equally spaced particles (with spacing h and mass m i ) on a string. Subject to a constant horizontal tension F.
17 Equation of motion for 4 particles: d 2 x 1 m 1 dt 2 = F x 1 h + F x 2 x 1 h d 2 x 2 m 2 dt 2 = F x 2 x 1 + F x 3 x 2 h h d 2 x 3 m 3 dt 2 = F x 3 x 2 + F x 4 x 3 h h d 2 x 4 m 4 dt 2 = F x 4 x 3 h In matrix form: d 2 x dt 2 = DAx x = [x 2 1, x 2, x 3, x 4 ] T A = D = diag(d 1, d 2, d 3, d 4 ) with d i = F m i h. F x 4 h
18 An MIEP Eigenvalues of DA: Necessarily real and nonnegative. Are the squares of the so called natural frequencies of the system. Inverse problem: Want to place the weight m i so that the system has a prescribed set of natural frequencies. A is symmetric and tridiagonal. D is diagonal. This is a multiplicative inverse eigenvalue problem. What is the range of possible natural frequencies? Can the system have arbitrary natural frequencies?
19 Control Theory ẋ(t) = Ax(t) + Bu(t). x(t) R n denotes the state of a certain physical system. u(t) R m denotes a control to the system. Want to select the input u(t) so that the dynamics of the resulting x(t) is driven into a certain desired state.
20 Feedback Control and Pole Assignment Problem State feedback. u(t) = }{{} F x(t) = ẋ(t) = (A + BF)x(t). gain matrix Output feedback. y(t) = }{{} C x(t) known u(t) = }{{} K y(t) output matrix = ẋ(t) = (A + BKC)x(t). Choose F or K so that the matrices A + BF or A + BKC has a prescribed set of eigenvalues.
21 Stochastic Inverse Spectrum Problem Construct a Markov chain with prescribed spectrum. Row-stochastic structure. No strings of symmetry. Eigenvalues can appear in complex conjugate pairs. A hard problem (Karpelevič 51, Minc 88). Need to know how the eigenvalues of n n stochastic matrices are distributed in C. The set Θ n consisting eigenvalues of n n stochastic matrices is completely characterized. The Karpelevič theorem characterizes only one complex value a time. When will two or more points in Θ n belong to the spectrum of the same stochastic matrix?
22 Basic Ideas Model Updating Structure Preservation Semi-definite Programming Conclusion Karpelevič s Theorem A number λ is an eigenvalue for a stochastic matrix if and only if it belongs to a region Θ n such as the one shown below for n =
23 Applications in Other Disciplines Applied Mechanics and structure design Construct a model of a mass-spring system with prescribed natural frequencies/modes. Applied physics Quantum mechanics. Geophysics. Numerical analysis Computing B-stable RK methods have real poles. Quadratures. Mathematical analysis Inverse Sturm-Liouville problems.
24 Brief History Inverse Sturm-Liouville problem: Ambartsumyan 29 Krein 33 Borg 46, Levinson 49 Gel fand&levitan 51 Kac 66 (Can one hear the shape of a drum?) Hochstadt 73, Barcilon 74, McLaughlin 76, Hald 78 Zayed 82, Issacson et al 83, McLaughlin 86, Andersson 88 Lowe et al 95, Rundell 97 Matrix theory: Downing&Householder 56, Mirsky 58 Hochstadt 67 de Oliveira 70, Hald 72, Golub 73, Friedland 77, de Boor&Golub 78 Biegler-König 81, Shapiro 83, Barcilon 86, Sun 86, Boley&Golub 87 Landau 94, Chu 98 (A list of 39 different types of IEPs.)
25 Applied Mechanics: Barcilon 74 Gottlieb 83, Gladwell 86 Ram 91, Gladwell 96, Nylen&Uhlig 97 Computation: Morel 76, Boley&Golud 78 Nocedal et al 83, Friedland et al 88, Laurie 88 Chu 90, Zhou&Dai 91, Trench 97, Xu 98 Chinese Scholars: Lots of contributions. Particularly strong in solvability and sensitivity analysis. Mostly published in Chinese, not accessible by the West. Quadratic Inverse Eigenvalue Problems: Real challenge. Research barely begins.
26 Establishing an Effective Model Precise mathematical models of physical systems are rarely available. Disturbance to the measurement. Imperfect characterization of the natural phenomenon. Discretization error.. Necessary to update a primitive model to attain consistency with empirical results.
27 The Updating Problem Given A structured quadratic pencil (M 0, C 0, K 0 ), A few of its associated eigenpairs {(λ i, u i )} k i=1 with k < 2n, Newly measured eigenpairs {(σ i, y i )} k i=1. Update (M 0, C 0, K 0 ) to (M, C, K ) so that The same structure is maintained, The subset {(λ i, u i )} k i=1 is replaced by {(σ i, y i )} k i=1 as k eigenpairs of (M, C, K ). Must there be some other relationship between the pencils (M, C, K ) and (M 0, C 0, K 0 )?
28 Structure Preserving The structure imposed on the model depends inherently on the underlying physical system. Physical feasibility requires that structure be preserved. Some examples: M is diagonal and positive; C and K are symmetric and banded; K is positive semi-definite. The inner-connectivity of an electronic circuit. The structure varies from problem to problem. The analysis of a solution varies from problem to problem? Would prefer to see a unified theory.
29 Spurious Eigenstructure Vibration parameters not related to the newly measured parameters should remain invariant. These parameters have already been proved to be acceptable in the original model. Do not wish to introduce new vibrations via updating. Simply do not know of any information about these parameters. No spill-over is desirable. Update (M 0, C 0, K 0 ) to (M, C, K ) so that The same structure is maintained, The subset {(λ i, u i )} k i=1 is replaced by {(σ i, y i )} k i=1 as k eigenpairs of (M, C, K ). The remaining 2n k eigenpairs, unmeasurable and unknown, are the same as those of the original (M 0, C 0, K 0 ). Can the no spill-over phenomenon be achieved?
30 Minimal or Robust Changes The solution to the updating problem is not unique. Wish to optimize the adjustment or robustness. minimize (M,C,K ) subject to M M 0 2 F + C C 0 2 F + K K 0 2 F structural and spectral constraints. Need to know the representation of a solution before we can search for the optimal solution. Would be difficult to parameterize the feasible set. No longer needed with SDP techniques.
31 Setup Describe the partial eigeninformation via (Λ, X) R k k R n k. Λ = diag{λ [2] 1,..., λ[2] l, λ 2l+1,..., λ k }, X = [x 1R, x 1I,..., x lr, x li, x 2l+1,..., x k ].» λ [2] αj β j = j R 2 2, β β j α j 0, for j = 1,..., l, j All given eigenvalues are simple. Look for real and structured matrices M, C and K so that the equation MXΛ 2 + CXΛ + KX = 0, (1) is satisfied. What structure and how many eigenpairs are needed?
32 Maintaining Symmetry (k n ) Assume the QR decomposition X = Q [ R 0 ]. Write Define S := RΛR 1. (Λ, X) (S, R). Q MQ = (This is a fixed matrix.) (Same eigeninformation.) [ ] [ ] [ ] M11 M 12, Q C11 C CQ = 12, Q K11 K KQ = 12. M 21 M 22 C 21 C 22 K 21 K 22 Look for quadratic pencils so that M 11 RΛ 2 + C 11 RΛ + K 11 R = 0, (a square pencil,) M 21 RΛ 2 + C 21 RΛ + K 21 R = 0, (a rectangular pencil.)
33 M = arbitrary SPD. Parametric Solution (Kuo-Lin-Xu 2005) M 11 = arbitrary SPD. (k(k + 1)/2) 2 M 21 = arbitrary = M 12. k(n k) M 22 = M 21 M 1 11 M 12 + arbitrary SPD. ((n k)(n k + 1)/2) 2 C 11 = (M 11 S + S M 11 + R ΩR 1 ). C 21 = arbitrary = C 12. k(n k) C 22 = arbitrary symmetric. ((n k)(n k + 1)/2) 2 K 11 = S M 11 S + R ΩΛR 1. K 21 = K 12 = (M 21S 2 + C 21 S). K 22 = arbitrary symmetric. ((n k)(n k + 1)/2) 2 Ω = Additional free parameters. k {[ ] [ ] } ξ1 η Ω := diag 1 ξl η,..., l, ξ η 1 ξ 1 η l ξ 2l+1,..., ξ k. l
34 Counting Freedom Symmetric M, C and K = 3 2n(n + 1) unknowns. Given k eigenpairs = nk equations. The system will be under-determined and, hence, solvable, so long as k < 3 2 (n + 1) eigenpairs are specified. If k n, Free choice of M, C 21, C 22, K 22 and Ω = 3 n(n + 1) nk free 2 parameters in the general solution.»» M11 M 12 C11 C, 12 M 21 M 22 C 21 C 22» K11 K, 12 K 21 K 22 What are solutions when k > n? Cannot guarantee positive semi-definiteness of C or K. Cannot deal with gyroscopic or other structured systems..
35 Maximal Allowable Eigenpairs (Chu-Datta-Lin-Xu 2008) Provides an upper bound of the allowable eigenpairs, { 3l + 1, if n = 2l, k max = 3l + 2, if n = 2l + 1. Parametric representation of all solutions up to k = k max. Shows that spill-over for damped system is generally unpreventable. No guarantee of positive semi-definiteness or other structures.
36 Maintain Positive Semi-definiteness Can do QIEP for M 0, C 0, and K 0, with complete information of eigenvalues and eigenvectors (Lancaster-Prells 2005). with all eigenvalues and eigenvectors corresponding to real eigenvalues (Lancaster 2007). Can do QIEP for M 0, and K 0, but no C, with partially prescribed eigenpairs (Chu-Kuo-Lin 2004). Cannot maintain SPD structure with fine grains of connectivity.
37 Eigenstructure Completion (Dong-Lin-Chu 2008) Spectral decomposition (Gohbert-Lancaster-Rodman 1980, Chu-Xu 2008) M = (XJH 1 X ) 1, C = MXJ 2 H 1 X M, K = MXJ 3 H 1 X M + CM 1 C. (J, X) = a standard pair eigeninformation. A particular nonsingular H. XH 1 X = 0. H J = (H J). H = H. Block diagonal with at most 2 2 blocks If all eigenvalues of J have negative real part, and if XJ 2 H 1 X 0, then generically M 0, K 0, and C 0.
38 Other Difficulties Most theory and algorithms are structure dependent. Unable to address finer grain of inner-connectivity. Unable to satisfy the requisite nonnegativity. Unable to characterize mixed-type structure. Unable to tackle local" model updating. Would like to have a unified theory and/or robust algorithms. Literature moves in diverse directions.
39 What is SDP? Primal problem: Minimize X C X Subject to A i X = b i, i = 1,..., m X 0, C, A 1,..., A m = given symmetric matrices. = Frobenius inner product. Dual problem: Maximize y,s Subject to b y m i=1 y ia i + S = C, S 0,
40 Powerful Mathematical Programming Tool SDP is versatile in modeling problems arising in broad discipline areas. Combinatorial optimization, Min-max eigenvalue problems,... Controller synthesis, CAD design, Network queueing, QIEPs. SDP has a considerable similarity to the linear programming. Vector of variables in LP Symmetric matrices. Nonnegativity constraints Positive semi-definite constraints. Duality theory, algorithms, convergence and polynomial time-complexity for LP can be extended mechanically" to SDP.
41 Research Activities Profusion of research results: The book, Handbook of Semi-definite Programming, has 877 references. Online bibliography has 1333 references, available at http: //liinwww.ira.uka.de/bbliography/math/psd.html Many well developed MATLAB-based software packages. SeDuMe (Sturm), downloadable at SDPT3 (Toh-Tütüncü-Todd), downloadable at Useful interface. YALMIP (Löfberb), downloadable at
42 Solvability Test % Define symmetric variables sm = sdpvar(n,n); sc = sdpvar(n,n); sk = sdpvar(n,n); % Specify equality and definiteness constraints F = set(sm*x*lambda^2+sc*x*lambda+sk*x == zeros(n,k)); F = F + set(sm >= 0) + set(sc >= 0) + set(sk >= 0); % Select SDPT3 as the solver ops = sdpsettings( solver, sdpt3 );
43 % Invoke SDP solver solvesdp(f,[],ops); % Retrieve numerical solution M = double(sm) C = double(sc) K = double(sk) % Check feasibility checkset(f);
44 Make-or-break Calculation For solvability test, no objective function is needed. The feasibility set is convex, containing at least the trivial solution. An appropriate normalization is needed to assess the final construction. SPD algorithm is simple but efficient. Can be exploited to answer many theoretical questions. How many arbitrarily prescribed eigenstructure (all eigenvalues must have negative real part) can be completed for a QIEP with positive semi-definite coefficient matrices?
45 Random Eigenstructure Completion (n = 10) k c k p k n k t
46 Handling Gyroscopic Constraints : % Define skew-symmetric variables sg = sdpvar{n,n, skew, real ); sn = sdpvar(n,n, skew, real ); % Modify equality constraints F = set(sm*x*lambda^2+(sc+sg)*x*lambda+(sk+sn)*x ==... zeros(n,k)); :
47 Versatility and Generality Any of the skew-symmetric G or N adds n(n 1) 2 extra variables. Would increase solvability by allowing a larger number k of prescribed eigenpairs. Other types of coefficients, such as Toeplitz, Hankel, palindromic, sparsity pattern, or fixed entries can easily be incorporated into the SDP formulations for QIEPs. Being linear in coefficient matrices, the SDP approach can be used for matrix polynomials of general degrees.
48 Connectivity and Nonnegativity Constraints % Define physical parameters sm = sdpvar(4,1); sc = sdpvar(3,1); sk = sdpvar(5,1); % Define connectivity constraints sm = diag(sm); sc = [sc(1)+sc(2) 0 -sc(2) 0; ; -sc(2) 0 sc(2)+sc(3) -sc(3); 0 0 -sc(3) sc(3)]; sk = [sk(1)+sk(2)+sk(5) -sk(2) -sk(5) 0; -sk(2) sk(2)+sk(3) -sk(3) 0; -sk(5) -sk(3) sk(3)+sk(4)+sk(5) -sk(4); 0 0 -sk(4) sk(4)]
49 % Define equality and nonnegativity constraints F = set(sm*x*lambda^2+sc*x*lambda+sk*x == zeros(n,k)); F = F + set(sm > 0) + set(sc > 0) + set(sk > 0); % Normalize the first mass (or others) solvesdp(f,(sm(1,1)-1)^2,ops);
50 Convex Programming The SDP formulation is a convex programming so long as the artificial objective function is convex. The SDP computation tells Either a completely successful reconstruction" with positive physical parameters, or An utterly disastrous failure" giving rise to either the trivial solution or a degenerated system.
51 Model Updating % Define variables/constraints, structured if needed sdpvar um uc uk : F =... % Include rotated Lorentz cones F = F + set(rcone(reshape((m_0-sm),n^2,1),um,1/2)); F = F + set(rcone(reshape((c_0-sc),n^2,1),uc,1/2)); F = F + set(rcone(reshape((k_0-sk),n^2,1),uk,1/2)); % Define objective function and call for SOCP solvers sdpsolve(f,um + uc + uk,ops);
52 Harwell-Boeing BCSST*01 Data Generalized eigenvalue problem of size (λm 0 K 0 )x = 0, Has 24 positive eigenvalues and the remaining 24 at positive infinity. Improperly scaling: M 0 F and K 0 F Eigenvalue computation would be particularly challenging. (λm 0 K 0 )x 2 of the 24 real-valued eigenvalues and eigenvectors computed by MATLAB ranges from to
53 Sparsity Pattern of BCSST*01 Data 0 M 0 0 K nz = nz = 400
54 Updating Procedures Update (λ, x) to (µ, x). µ stands for a one-sided" random perturbation of eigenvalue λ. with σ = 10 p, p = 1,, 10. µ = λ(1 + σ randn(1) ) Two steps involved in the updating procedure: (µ, x) feasibility (M, K ) optimization (M 0, K 0 ). Measure feasibility by (µm K )x 2. Search for minimal update.
55 Performance of Updating BCSST*01 Data Residules over 20 samples of μ Objective values over 20 samples of μ 4 2 log 10 ( (μm K)x 2) log 10 (J (M,K)) p p
56 Conclusion The inverse eigenvalue problem is a diverse area full of research interests and activities. The inverse problem, expressing the parameters in term of the behavior, is just as important as the direct problem in applications. (The AMS has recently added a new number 65F18 for inverse eigenvalue problem to its Subject Classifications.) Model updating can be embedded as a quadratic inverse eigenvalue problem. The general solution to the inverse eigenvalue problem for adjoint quadratic pencils with up to 3 2 (n + 1) given eigenpairs is completely known. Structured problems are more complicated (and needed!) SDP techniques offers a unified and effectual avenue of attack on the QIEPs in general and the MUPs in particular.
Inverse Eigenvalue Problem with Non-simple Eigenvalues for Damped Vibration Systems
Journal of Informatics Mathematical Sciences Volume 1 (2009), Numbers 2 & 3, pp. 91 97 RGN Publications (Invited paper) Inverse Eigenvalue Problem with Non-simple Eigenvalues for Damped Vibration Systems
More informationInverse Eigenvalue Problems: Theory, Algorithms, and Applications
Inverse Eigenvalue Problems: Theory, Algorithms, and Applications Moody T. Chu North Carolina State University Gene H. Golub Stanford University OXPORD UNIVERSITY PRESS List of Acronyms List of Figures
More informationInverse Eigenvalue Problems: Theory and Applications
Inverse Eigenvalue Problems: Theory and Applications A Series of Lectures to be Presented at IRMA, CNR, Bari, Italy Moody T. Chu (Joint with Gene Golub) Department of Mathematics North Carolina State University
More informationComputational Methods for Feedback Control in Damped Gyroscopic Second-order Systems 1
Computational Methods for Feedback Control in Damped Gyroscopic Second-order Systems 1 B. N. Datta, IEEE Fellow 2 D. R. Sarkissian 3 Abstract Two new computationally viable algorithms are proposed for
More informationOn an Inverse Problem for a Quadratic Eigenvalue Problem
International Journal of Difference Equations ISSN 0973-6069, Volume 12, Number 1, pp. 13 26 (2017) http://campus.mst.edu/ijde On an Inverse Problem for a Quadratic Eigenvalue Problem Ebru Ergun and Adil
More informationPartial Eigenvalue Assignment in Linear Systems: Existence, Uniqueness and Numerical Solution
Partial Eigenvalue Assignment in Linear Systems: Existence, Uniqueness and Numerical Solution Biswa N. Datta, IEEE Fellow Department of Mathematics Northern Illinois University DeKalb, IL, 60115 USA e-mail:
More informationQuadratic Inverse Eigenvalue Problems, Active Vibration Control and Model Updating
Quadratic Inverse Eigenvalue Problems, Active Vibration Control and Model Updating Biswa N. Datta,1,2 Northern Illinois University, Mathematical Sciences, Watson Hall 320 DeKalb, IL 60115-2888, USA Vadim
More informationUpdating quadratic models with no spillover effect on unmeasured spectral data
INSTITUTE OF PHYSICS PUBLISHING Inverse Problems 3 (007) 43 56 INVERSE PROBLEMS doi:0.088/066-56/3//03 Updating quadratic models with no spillover effect on unmeasured spectral data Moody T Chu,4, Wen-Wei
More informationOutline Introduction: Problem Description Diculties Algebraic Structure: Algebraic Varieties Rank Decient Toeplitz Matrices Constructing Lower Rank St
Structured Lower Rank Approximation by Moody T. Chu (NCSU) joint with Robert E. Funderlic (NCSU) and Robert J. Plemmons (Wake Forest) March 5, 1998 Outline Introduction: Problem Description Diculties Algebraic
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 informationDefinite versus Indefinite Linear Algebra. Christian Mehl Institut für Mathematik TU Berlin Germany. 10th SIAM Conference on Applied Linear Algebra
Definite versus Indefinite Linear Algebra Christian Mehl Institut für Mathematik TU Berlin Germany 10th SIAM Conference on Applied Linear Algebra Monterey Bay Seaside, October 26-29, 2009 Indefinite Linear
More informationThe Eigenvalue Shift Technique and Its Eigenstructure Analysis of a Matrix
The Eigenvalue Shift Technique and Its Eigenstructure Analysis of a Matrix Chun-Yueh Chiang Center for General Education, National Formosa University, Huwei 632, Taiwan. Matthew M. Lin 2, Department of
More informationResearch overview. Seminar September 4, Lehigh University Department of Industrial & Systems Engineering. Research overview.
Research overview Lehigh University Department of Industrial & Systems Engineering COR@L Seminar September 4, 2008 1 Duality without regularity condition Duality in non-exact arithmetic 2 interior point
More informationLecture Note 5: Semidefinite Programming for Stability Analysis
ECE7850: Hybrid Systems:Theory and Applications Lecture Note 5: Semidefinite Programming for Stability Analysis Wei Zhang Assistant Professor Department of Electrical and Computer Engineering Ohio State
More informationABSTRACT QUADRATIC INVERSE EIGENVALUE PROBLEMS: THEORY, METHODS, AND APPLICATIONS
ABSTRACT QUADRATIC INVERSE EIGENVALUE PROBLEMS: THEORY, METHODS, AND APPLICATIONS Vadim Olegovich Sokolov, Ph.D. Department of Mathematical Sciences Northern Illinois University, 2008 Biswa Nath Datta,
More informationAppendix A Solving Linear Matrix Inequality (LMI) Problems
Appendix A Solving Linear Matrix Inequality (LMI) Problems In this section, we present a brief introduction about linear matrix inequalities which have been used extensively to solve the FDI problems described
More informationCSC Linear Programming and Combinatorial Optimization Lecture 10: Semidefinite Programming
CSC2411 - Linear Programming and Combinatorial Optimization Lecture 10: Semidefinite Programming Notes taken by Mike Jamieson March 28, 2005 Summary: In this lecture, we introduce semidefinite programming
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 informationLinear Matrix Inequalities in Control
Linear Matrix Inequalities in Control Delft Center for Systems and Control (DCSC) Delft University of Technology The Netherlands Department of Electrical Engineering Eindhoven University of Technology
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 informationEE364 Review Session 4
EE364 Review Session 4 EE364 Review Outline: convex optimization examples solving quasiconvex problems by bisection exercise 4.47 1 Convex optimization problems we have seen linear programming (LP) quadratic
More information15. Conic optimization
L. Vandenberghe EE236C (Spring 216) 15. Conic optimization conic linear program examples modeling duality 15-1 Generalized (conic) inequalities Conic inequality: a constraint x K where K is a convex cone
More information6-1 The Positivstellensatz P. Parrilo and S. Lall, ECC
6-1 The Positivstellensatz P. Parrilo and S. Lall, ECC 2003 2003.09.02.10 6. The Positivstellensatz Basic semialgebraic sets Semialgebraic sets Tarski-Seidenberg and quantifier elimination Feasibility
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 informationA Control Methodology for Constrained Linear Systems Based on Positive Invariance of Polyhedra
A Control Methodology for Constrained Linear Systems Based on Positive Invariance of Polyhedra Jean-Claude HENNET LAAS-CNRS Toulouse, France Co-workers: Marina VASSILAKI University of Patras, GREECE Jean-Paul
More informationMultiobjective Optimization Applied to Robust H 2 /H State-feedback Control Synthesis
Multiobjective Optimization Applied to Robust H 2 /H State-feedback Control Synthesis Eduardo N. Gonçalves, Reinaldo M. Palhares, and Ricardo H. C. Takahashi Abstract This paper presents an algorithm for
More informationReal Symmetric Matrices and Semidefinite Programming
Real Symmetric Matrices and Semidefinite Programming Tatsiana Maskalevich Abstract Symmetric real matrices attain an important property stating that all their eigenvalues are real. This gives rise to many
More informationSecond-order cone programming
Outline Second-order cone programming, PhD Lehigh University Department of Industrial and Systems Engineering February 10, 2009 Outline 1 Basic properties Spectral decomposition The cone of squares The
More information1 Review of simple harmonic oscillator
MATHEMATICS 7302 (Analytical Dynamics YEAR 2017 2018, TERM 2 HANDOUT #8: COUPLED OSCILLATIONS AND NORMAL MODES 1 Review of simple harmonic oscillator In MATH 1301/1302 you studied the simple harmonic oscillator:
More informationLinear Algebra. P R E R E Q U I S I T E S A S S E S S M E N T Ahmad F. Taha August 24, 2015
THE UNIVERSITY OF TEXAS AT SAN ANTONIO EE 5243 INTRODUCTION TO CYBER-PHYSICAL SYSTEMS P R E R E Q U I S I T E S A S S E S S M E N T Ahmad F. Taha August 24, 2015 The objective of this exercise is to assess
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 informationON INVERSE QUADRATIC EIGENVALUE PROBLEMS WITH PARTIALLY PRESCRIBED EIGENSTRUCTURE
ON INVERSE QUADRATIC EIGENVALUE PROBLEMS WITH PARTIALLY PRESCRIBED EIGENSTRUCTURE MOODY T. CHU, YUEN-CHENG KUO, AND WEN-WEI LIN Abstract. The inverse eigenvalue problem of constructing real and symmetric
More informationarzelier
COURSE ON LMI OPTIMIZATION WITH APPLICATIONS IN CONTROL PART II.1 LMIs IN SYSTEMS CONTROL STATE-SPACE METHODS STABILITY ANALYSIS Didier HENRION www.laas.fr/ henrion henrion@laas.fr Denis ARZELIER www.laas.fr/
More information1 Number Systems and Errors 1
Contents 1 Number Systems and Errors 1 1.1 Introduction................................ 1 1.2 Number Representation and Base of Numbers............. 1 1.2.1 Normalized Floating-point Representation...........
More informationLinear Algebra Methods for Data Mining
Linear Algebra Methods for Data Mining Saara Hyvönen, Saara.Hyvonen@cs.helsinki.fi Spring 2007 1. Basic Linear Algebra Linear Algebra Methods for Data Mining, Spring 2007, University of Helsinki Example
More informationPeter J. Dukes. 22 August, 2012
22 August, 22 Graph decomposition Let G and H be graphs on m n vertices. A decompostion of G into copies of H is a collection {H i } of subgraphs of G such that each H i = H, and every edge of G belongs
More informationSemidefinite and Second Order Cone Programming Seminar Fall 2001 Lecture 5
Semidefinite and Second Order Cone Programming Seminar Fall 2001 Lecture 5 Instructor: Farid Alizadeh Scribe: Anton Riabov 10/08/2001 1 Overview We continue studying the maximum eigenvalue SDP, and generalize
More informationConvex Optimization of Graph Laplacian Eigenvalues
Convex Optimization of Graph Laplacian Eigenvalues Stephen Boyd Abstract. We consider the problem of choosing the edge weights of an undirected graph so as to maximize or minimize some function of the
More informationChapter 3 Parameterized Inverse Eigenvalue Problems
Chapter 3 Parameterized Inverse Eigenvalue Problems Overview. General Results. Additive Inverse Eigenvalue Problems. Multiplicative Inverse Eigenvalue Problems. 51 52 Parameterized Inverse Eigenvalue Problems
More informationInverse eigenvalue problems involving multiple spectra
Inverse eigenvalue problems involving multiple spectra G.M.L. Gladwell Department of Civil Engineering University of Waterloo Waterloo, Ontario, Canada N2L 3G1 ggladwell@uwaterloo.ca URL: http://www.civil.uwaterloo.ca/ggladwell
More informationInverse Eigenvalue Problems and Their Associated Approximation Problems for Matrices with J-(Skew) Centrosymmetry
Inverse Eigenvalue Problems and Their Associated Approximation Problems for Matrices with -(Skew) Centrosymmetry Zhong-Yun Liu 1 You-Cai Duan 1 Yun-Feng Lai 1 Yu-Lin Zhang 1 School of Math., Changsha University
More informationFast Algorithms for SDPs derived from the Kalman-Yakubovich-Popov Lemma
Fast Algorithms for SDPs derived from the Kalman-Yakubovich-Popov Lemma Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University 8 September 2003 European Union RTN Summer School on Multi-Agent
More informationCONTROL DESIGN FOR SET POINT TRACKING
Chapter 5 CONTROL DESIGN FOR SET POINT TRACKING In this chapter, we extend the pole placement, observer-based output feedback design to solve tracking problems. By tracking we mean that the output is commanded
More informationFast ADMM for Sum of Squares Programs Using Partial Orthogonality
Fast ADMM for Sum of Squares Programs Using Partial Orthogonality Antonis Papachristodoulou Department of Engineering Science University of Oxford www.eng.ox.ac.uk/control/sysos antonis@eng.ox.ac.uk with
More informationProblem Set 5 Solutions 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Problem Set 5 Solutions The problem set deals with Hankel
More informationNonlinear Eigenvalue Problems: An Introduction
Nonlinear Eigenvalue Problems: An Introduction Cedric Effenberger Seminar for Applied Mathematics ETH Zurich Pro*Doc Workshop Disentis, August 18 21, 2010 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction
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 informationControl Systems Design, SC4026. SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft
Control Systems Design, SC4026 SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft Lecture 4 Controllability (a.k.a. Reachability) vs Observability Algebraic Tests (Kalman rank condition & Hautus test) A few
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 informationm i=1 c ix i i=1 F ix i F 0, X O.
What is SDP? for a beginner of SDP Copyright C 2005 SDPA Project 1 Introduction This note is a short course for SemiDefinite Programming s SDP beginners. SDP has various applications in, for example, control
More informationIntroduction. Chapter One
Chapter One Introduction The aim of this book is to describe and explain the beautiful mathematical relationships between matrices, moments, orthogonal polynomials, quadrature rules and the Lanczos and
More informationConvex computation of the region of attraction for polynomial control systems
Convex computation of the region of attraction for polynomial control systems Didier Henrion LAAS-CNRS Toulouse & CTU Prague Milan Korda EPFL Lausanne Region of Attraction (ROA) ẋ = f (x,u), x(t) X, u(t)
More informationRegion of attraction approximations for polynomial dynamical systems
Region of attraction approximations for polynomial dynamical systems Milan Korda EPFL Lausanne Didier Henrion LAAS-CNRS Toulouse & CTU Prague Colin N. Jones EPFL Lausanne Region of Attraction (ROA) ẋ(t)
More informationDescriptor system techniques in solving H 2 -optimal fault detection problems
Descriptor system techniques in solving H 2 -optimal fault detection problems Andras Varga German Aerospace Center (DLR) DAE 10 Workshop Banff, Canada, October 25-29, 2010 Outline approximate fault detection
More informationEL 625 Lecture 10. Pole Placement and Observer Design. ẋ = Ax (1)
EL 625 Lecture 0 EL 625 Lecture 0 Pole Placement and Observer Design Pole Placement Consider the system ẋ Ax () The solution to this system is x(t) e At x(0) (2) If the eigenvalues of A all lie in the
More informationCourse Outline. FRTN10 Multivariable Control, Lecture 13. General idea for Lectures Lecture 13 Outline. Example 1 (Doyle Stein, 1979)
Course Outline FRTN Multivariable Control, Lecture Automatic Control LTH, 6 L-L Specifications, models and loop-shaping by hand L6-L8 Limitations on achievable performance L9-L Controller optimization:
More informationCONSTRUCTION OF A COMPLEX JACOBI MATRIX FROM TWO-SPECTRA
Hacettepe Journal of Mathematics and Statistics Volume 40(2) (2011), 297 303 CONSTRUCTION OF A COMPLEX JACOBI MATRIX FROM TWO-SPECTRA Gusein Sh Guseinov Received 21:06 :2010 : Accepted 25 :04 :2011 Abstract
More informationConvex Optimization 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Convex Optimization 1 Many optimization objectives generated
More informationAssignment 1: From the Definition of Convexity to Helley Theorem
Assignment 1: From the Definition of Convexity to Helley Theorem Exercise 1 Mark in the following list the sets which are convex: 1. {x R 2 : x 1 + i 2 x 2 1, i = 1,..., 10} 2. {x R 2 : x 2 1 + 2ix 1x
More informationLecture 2 INF-MAT : , LU, symmetric LU, Positve (semi)definite, Cholesky, Semi-Cholesky
Lecture 2 INF-MAT 4350 2009: 7.1-7.6, LU, symmetric LU, Positve (semi)definite, Cholesky, Semi-Cholesky Tom Lyche and Michael Floater Centre of Mathematics for Applications, Department of Informatics,
More informationChapter Two Elements of Linear Algebra
Chapter Two Elements of Linear Algebra Previously, in chapter one, we have considered single first order differential equations involving a single unknown function. In the next chapter we will begin to
More informationAn Infeasible Interior-Point Algorithm with full-newton Step for Linear Optimization
An Infeasible Interior-Point Algorithm with full-newton Step for Linear Optimization H. Mansouri M. Zangiabadi Y. Bai C. Roos Department of Mathematical Science, Shahrekord University, P.O. Box 115, Shahrekord,
More information20D - Homework Assignment 5
Brian Bowers TA for Hui Sun MATH D Homework Assignment 5 November 8, 3 D - Homework Assignment 5 First, I present the list of all matrix row operations. We use combinations of these steps to row reduce
More informationEigenvalues and Eigenvectors
/88 Chia-Ping Chen Department of Computer Science and Engineering National Sun Yat-sen University Linear Algebra Eigenvalue Problem /88 Eigenvalue Equation By definition, the eigenvalue equation for matrix
More informationA CONIC DANTZIG-WOLFE DECOMPOSITION APPROACH FOR LARGE SCALE SEMIDEFINITE PROGRAMMING
A CONIC DANTZIG-WOLFE DECOMPOSITION APPROACH FOR LARGE SCALE SEMIDEFINITE PROGRAMMING Kartik Krishnan Advanced Optimization Laboratory McMaster University Joint work with Gema Plaza Martinez and Tamás
More informationRobust Multi-Objective Control for Linear Systems
Robust Multi-Objective Control for Linear Systems Elements of theory and ROMULOC toolbox Dimitri PEAUCELLE & Denis ARZELIER LAAS-CNRS, Toulouse, FRANCE Part of the OLOCEP project (includes GloptiPoly)
More informationAA242B: MECHANICAL VIBRATIONS
AA242B: MECHANICAL VIBRATIONS 1 / 50 AA242B: MECHANICAL VIBRATIONS Undamped Vibrations of n-dof Systems These slides are based on the recommended textbook: M. Géradin and D. Rixen, Mechanical Vibrations:
More informationSymmetric Tridiagonal Inverse Quadratic Eigenvalue Problems with Partial Eigendata
Symmetric Tridiagonal Inverse Quadratic Eigenvalue Problems with Partial Eigendata Zheng-Jian Bai Revised: October 18, 2007 Abstract In this paper we concern the inverse problem of constructing the n-by-n
More informationIntroduction to Applied Linear Algebra with MATLAB
Sigam Series in Applied Mathematics Volume 7 Rizwan Butt Introduction to Applied Linear Algebra with MATLAB Heldermann Verlag Contents Number Systems and Errors 1 1.1 Introduction 1 1.2 Number Representation
More informationTheory and Computations of Some Inverse Eigenvalue Problems for the Quadratic Pencil
Contemporary Mathematics Theory and Computations of Some Inverse Eigenvalue Problems for the Quadratic Pencil Biswa N. Datta and Daniil R. Sarkissian Abstract. The paper reviews, both from theoretical
More informationMatrix stabilization using differential equations.
Matrix stabilization using differential equations. Nicola Guglielmi Universitá dell Aquila and Gran Sasso Science Institute, Italia NUMOC-2017 Roma, 19 23 June, 2017 Inspired by a joint work with Christian
More informationGlobal Optimization of Polynomials
Semidefinite Programming Lecture 9 OR 637 Spring 2008 April 9, 2008 Scribe: Dennis Leventhal Global Optimization of Polynomials Recall we were considering the problem min z R n p(z) where p(z) is a degree
More informationSelected Examples of CONIC DUALITY AT WORK Robust Linear Optimization Synthesis of Linear Controllers Matrix Cube Theorem A.
. Selected Examples of CONIC DUALITY AT WORK Robust Linear Optimization Synthesis of Linear Controllers Matrix Cube Theorem A. Nemirovski Arkadi.Nemirovski@isye.gatech.edu Linear Optimization Problem,
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 informationRecent Advances in the Numerical Solution of Quadratic Eigenvalue Problems
Recent Advances in the Numerical Solution of Quadratic Eigenvalue Problems Françoise Tisseur School of Mathematics The University of Manchester ftisseur@ma.man.ac.uk http://www.ma.man.ac.uk/~ftisseur/
More informationConvexity of the Joint Numerical Range
Convexity of the Joint Numerical Range Chi-Kwong Li and Yiu-Tung Poon October 26, 2004 Dedicated to Professor Yik-Hoi Au-Yeung on the occasion of his retirement. Abstract Let A = (A 1,..., A m ) be an
More informationSemidefinite Programming
Semidefinite Programming Basics and SOS Fernando Mário de Oliveira Filho Campos do Jordão, 2 November 23 Available at: www.ime.usp.br/~fmario under talks Conic programming V is a real vector space h, i
More informationResearch Reports on Mathematical and Computing Sciences
ISSN 1342-284 Research Reports on Mathematical and Computing Sciences Exploiting Sparsity in Linear and Nonlinear Matrix Inequalities via Positive Semidefinite Matrix Completion Sunyoung Kim, Masakazu
More informationRecurrence Relations and Fast Algorithms
Recurrence Relations and Fast Algorithms Mark Tygert Research Report YALEU/DCS/RR-343 December 29, 2005 Abstract We construct fast algorithms for decomposing into and reconstructing from linear combinations
More informationControl Systems Design, SC4026. SC4026 Fall 2010, dr. A. Abate, DCSC, TU Delft
Control Systems Design, SC4026 SC4026 Fall 2010, dr. A. Abate, DCSC, TU Delft Lecture 4 Controllability (a.k.a. Reachability) and Observability Algebraic Tests (Kalman rank condition & Hautus test) A few
More informationTutorials in Optimization. Richard Socher
Tutorials in Optimization Richard Socher July 20, 2008 CONTENTS 1 Contents 1 Linear Algebra: Bilinear Form - A Simple Optimization Problem 2 1.1 Definitions........................................ 2 1.2
More informationCertifying the Global Optimality of Graph Cuts via Semidefinite Programming: A Theoretic Guarantee for Spectral Clustering
Certifying the Global Optimality of Graph Cuts via Semidefinite Programming: A Theoretic Guarantee for Spectral Clustering Shuyang Ling Courant Institute of Mathematical Sciences, NYU Aug 13, 2018 Joint
More informationApplied Linear Algebra
Applied Linear Algebra Peter J. Olver School of Mathematics University of Minnesota Minneapolis, MN 55455 olver@math.umn.edu http://www.math.umn.edu/ olver Chehrzad Shakiban Department of Mathematics University
More informationABSTRACT. Name: Daniil R. Sarkissian Department: Mathematical Sciences
ABSTRACT Name: Daniil R. Sarkissian Department: Mathematical Sciences Title: Theory and Computations of Partial Eigenvalue and Eigenstructure Assignment Problems in Matrix Second-Order and Distributed-Parameter
More informationLinear Matrix Inequalities in Robust Control. Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University MTNS 2002
Linear Matrix Inequalities in Robust Control Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University MTNS 2002 Objective A brief introduction to LMI techniques for Robust Control Emphasis on
More information21 Linear State-Space Representations
ME 132, Spring 25, UC Berkeley, A Packard 187 21 Linear State-Space Representations First, let s describe the most general type of dynamic system that we will consider/encounter in this class Systems may
More informationLecture 4 Eigenvalue problems
Lecture 4 Eigenvalue problems Weinan E 1,2 and Tiejun Li 2 1 Department of Mathematics, Princeton University, weinan@princeton.edu 2 School of Mathematical Sciences, Peking University, tieli@pku.edu.cn
More informationModeling and Analysis of Dynamic Systems
Modeling and Analysis of Dynamic Systems by Dr. Guillaume Ducard Fall 2016 Institute for Dynamic Systems and Control ETH Zurich, Switzerland based on script from: Prof. Dr. Lino Guzzella 1/33 Outline 1
More informationConvex sets, conic matrix factorizations and conic rank lower bounds
Convex sets, conic matrix factorizations and conic rank lower bounds Pablo A. Parrilo Laboratory for Information and Decision Systems Electrical Engineering and Computer Science Massachusetts Institute
More informationAn Algorithm for Solving the Convex Feasibility Problem With Linear Matrix Inequality Constraints and an Implementation for Second-Order Cones
An Algorithm for Solving the Convex Feasibility Problem With Linear Matrix Inequality Constraints and an Implementation for Second-Order Cones Bryan Karlovitz July 19, 2012 West Chester University of Pennsylvania
More informationDELFT UNIVERSITY OF TECHNOLOGY
DELFT UNIVERSITY OF TECHNOLOGY REPORT -09 Computational and Sensitivity Aspects of Eigenvalue-Based Methods for the Large-Scale Trust-Region Subproblem Marielba Rojas, Bjørn H. Fotland, and Trond Steihaug
More informationLecture 1. 1 Conic programming. MA 796S: Convex Optimization and Interior Point Methods October 8, Consider the conic program. min.
MA 796S: Convex Optimization and Interior Point Methods October 8, 2007 Lecture 1 Lecturer: Kartik Sivaramakrishnan Scribe: Kartik Sivaramakrishnan 1 Conic programming Consider the conic program min s.t.
More informationProblem structure in semidefinite programs arising in control and signal processing
Problem structure in semidefinite programs arising in control and signal processing Lieven Vandenberghe Electrical Engineering Department, UCLA Joint work with: Mehrdad Nouralishahi, Tae Roh Semidefinite
More informationEventually reducible matrix, eventually nonnegative matrix, eventually r-cyclic
December 15, 2012 EVENUAL PROPERIES OF MARICES LESLIE HOGBEN AND ULRICA WILSON Abstract. An eventual property of a matrix M C n n is a property that holds for all powers M k, k k 0, for some positive integer
More informationACM/CMS 107 Linear Analysis & Applications Fall 2016 Assignment 4: Linear ODEs and Control Theory Due: 5th December 2016
ACM/CMS 17 Linear Analysis & Applications Fall 216 Assignment 4: Linear ODEs and Control Theory Due: 5th December 216 Introduction Systems of ordinary differential equations (ODEs) can be used to describe
More informationThe quadratic eigenvalue problem (QEP) is to find scalars λ and nonzero vectors u satisfying
I.2 Quadratic Eigenvalue Problems 1 Introduction The quadratic eigenvalue problem QEP is to find scalars λ and nonzero vectors u satisfying where Qλx = 0, 1.1 Qλ = λ 2 M + λd + K, M, D and K are given
More informationMulti Degrees of Freedom Systems
Multi Degrees of Freedom Systems MDOF s http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano March 9, 07 Outline, a System
More informationLinear Programming: Simplex
Linear Programming: Simplex Stephen J. Wright 1 2 Computer Sciences Department, University of Wisconsin-Madison. IMA, August 2016 Stephen Wright (UW-Madison) Linear Programming: Simplex IMA, August 2016
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 informationLecture: Cone programming. Approximating the Lorentz cone.
Strong relaxations for discrete optimization problems 10/05/16 Lecture: Cone programming. Approximating the Lorentz cone. Lecturer: Yuri Faenza Scribes: Igor Malinović 1 Introduction Cone programming is
More information