On some properties of the Lyapunov equation for damped systems
|
|
- Prudence Johns
- 5 years ago
- Views:
Transcription
1 On some properties of the Lyapuno equation for damped systems Ninosla Truhar, Uniersity of Osijek, Department of Mathematics, 31 Osijek, Croatia 1 ntruhar@mathos.hr Krešimir Veselić Lehrgebiet Mathematische Physik, Fernuniersität, 5884 Hagen, Germany kresimir.eselic@fernuni-hagen.de Abstract We consider a damped linear ibrational system whose dampers depend linearly on the iscosity parameter. We show that the trace of the corresponding Lyapuno solution can be represented as a rational function of whose poles are the eigenalues of a certain skew symmetric matrix. This makes it possible to derie an asymptotic expansion of the solution in the neighborhood of zero (small damping). 1 Introduction We consider a damped linear ibrational system Mẍ + Cẋ + Kx = (1.1) the matrices M, C, K (mass, damping, stiffness) are symmetric, M, K are positie definite and C is positie semidefinite. If internal damping is neglected C has often small rank as it describes a few dampers built in to calm down dangerous ibrations. Often C has the form C = C is a ariable iscosity and C describes the geometry of a damper. C will hae rank 1,2 or 3 according to whether the damper can exhibit linear, planar or spatial displacements. 1 A part of this work was written while the author was isiting researcher on the Lehrgebiet Mathematische Physik, Fernuniersität, Hagen 1
2 An example is the so-called n-mass oscillator or oscillator ladder (Fig. 1) M = diag(m 1,m 2,...,m n ) k 1 k 1 k 1 k 1 + k 2 k 2 K = , k n 2 k n 2 + k n 1 k n 1 k n 1 k n 1 + k n C = e 1 e T 1 + (e 3 e 2 )(e 3 e 2 ) T. (1.2) Figure 1: The n-mass oscillator with two dampers Here m i > are the masses, k i > the spring constants or stiffnesses, e i is the i-th canonical basis ector, and is the iscosity of the damper applied on the i-th mass. After the substitution y 1 = ΩΦ 1 x, y 2 = Φ 1 ẋ, (1.3) KΦ = MΦΩ 2, Φ T MΦ = I (1.4) Ω = diag(ω 1,...,ω n ), ω 1 <... < ω n (1.5) is the eigenreduction of the symmetric positie definite matrix pair K,M, the system (1.1) goes oer into y1 ẏ = Ay, y =, (1.6) Then E = 1 2 y 2 Ω A =, D = Φ T CΦ. (1.7) Ω D (ẋt Mẋ + x T Kx ) dt = 2 y 2 dt = e A t y 2 dt, (1.8)
3 y is the initial data. Thus, soles the Lyapuno equation E E(y ) = y T Xy, X = e AT t e At dt (1.9) A T X + XA = I. (1.1) Our penalty function is obtained by aeraging E oer all initial data with the equal energy, that is, we form Ē = y T Xy dσ y =1 dσ is a probability measure on the unit sphere in R 2n. Since by the map X y T Xy dσ y =1 is gien a linear functional on the space of the symmetric matrices, by Riesz theorem there exists a symmetric matrix Z such that X y T Xy dσ = Tr(ZX), for all symmetric matrices X. y =1 Let y R 2n be arbitrary. Set X = yy T. Then y T Xy dσ = Tr(ZX) = Tr(Zyy T ) = Tr(y T Zy), y =1 hence Z is always positie semi-definite. For any gien measure there is a unique positie semidefinite matrix Z such that Ē = Tr(ZX). (1.11) For the measure σ generated by the Lebesgue measure (i.e. the usual surface measure) on R 2n, we obtain Z = 1 2nI. For the conenience of the reader, we gie a sketch of the proof: Recall, Z ij = y i y j σ(dy). S One can easily see using Minkowski formula (see [2]) that Z ij = y i y j σ(dy) = 1 2ε lim y i y j σ(dy), ε S d(y,s) ε 3
4 here S denotes the unit sphere in R 2n and d(y,s) is a corresponding distance. Obiously, Z ij = for i j and Z ii = Z jj, for i,j {1,2,...,2n}. Since 1 Z 11 + Z Z 2n2n = lim ε 2ε ol ( y R 2n : d(y,s) ε ) = 1, it follows Z = 1 2nI. The more details about the structure of the matrix Z one can find in [4]. We hae shown that y T Xy dσ = min is equialent to y =1 Tr(ZX) = min. (1.12) Z is a symmetric positie semidefinite matrix which may be normalized to hae a unit trace. If one is interested in damping a certain part of the spectrum of the matrix A (which is ery important in applications) then the matrix Z will hae a special structure. For example, let σ = σ 1 σ 2 σ 1 σ 2, σ 1 is a measure on the frequency subspace determined by ω ω max ω s generated by Lebesgue measure, that is σ 1 is a measure on the frequency subspace which corresponds to the eigenfrequencies (defined by (1.5)) ω 1,...ω s and σ 2 is Dirac measure on the complement. Then we obtain that the corresponding matrix Z has the form I s Z s = Z = 1 2s I s, (1.13) I s is the identity matrix of the dimension s which is defined by ω max = ω s. Here ω max = ω s is critical frequency with the property that the eigenfrequencies from (1.5) greater then ω s are not dangerous. Hence, we damp first s eigenfrequencies. The construction of Z from (1.13) is a generalization of the simplest case Z = 1 2n I. In [7] a simple solution of the problem (1.11) has been presented for A from (1.7) and rank(c) = 1. In particular, Tr(ZX()) = const + a + b, a,b >, (1.14) which made it possible to find the minimum by a simple formula explicitly. The case rank(c) > 1 seems to be essentially more difficult to handle. The main result of this paper is the explicit formula: X() = Ψ 1 Ψ Ψ 1 + s i=1 λ i (λ i Φ i Υ i ) λ 2 i + 2. (1.15) 4
5 Ψ 1, Ψ, Ψ 1, Φ i and Υ i are m m matrices and λ i are eigenalues of the pencil (A, D) A and D are matrices which correspond to the linear operators X A X + XA, X DX + XD, respectiely. Thus, technically, we hae turned the iscosity into the spectral parameter. Further, we hae obtained a simple formula for the Ψ 1 in (1.15). This matrix is responsible for the behaior of the solution X() in the neighborhood of zero (small damping): D Ψ 1 =, D D = (diag (D)) 1. We will use the following notation: matrices written in the simple Roman fonts, M, D or K for example will hae n 2 entries. Matrices written in the mathematical bold fonts, A, B will hae m 2 entries, m = 2n (that is A, B are matrices defined on the 2n-dimensional phase space). Finally, matrices written in the Blackboard bold fonts A, or D will hae more than m 2 entries. 2 The main result As we hae said in the Introduction, our aim is to obtain the solution X of the Lyapuno equation A() T X + XA() = I, (2.1) Z is defined by (1.13) and I is m m identity matrix. From (1.7) it follows that A() can be written as Ω 1 Ω 2 A() A D, A =, Ω i =... Ωn [ ] ωi ω i (2.2) and D = D D T,... d 11 d d 1r... D = d 21 d d 2r d n1 d n2... d nr, 5
6 d ij are entries of the matrix D = P, C = L C L T C, L C and P is the perfect shuffling permutation. Now, we proceed with soling equation (2.1). As it is well known, Lyapuno equation (2.1) is equialent to ([3, Theorem ]) ( I (A D) T + (A D) T I ) ec(x) = ec(i), (2.3) L T denotes the Kronecker product of L and T, and ec(i) is the ector formed by stacking the columns of I into one long ector. Further, we will need the following two m 2 m 2 matrices defined by A = I A T + A T I, D = I D D T + D D T I. (2.4) It is easy to show that D = D F D T F, D F = [ I D D I ]. (2.5) Now, using (2.5) and (2.4) it follows that solution ec(x) of equation (2.3) can be written as ec(x) = ( A D F D T F) 1 ec(i). (2.6) Obiously, there exists a unitary matrix U such that [ U T A U =, (2.7) Â] A is the skew-symmetric matrix corresponds with linear operator defined in (2.4) and  is a non-singular block diagonal matrix defined by µi  = diag(ξ 1,...,Ξ m2 ) Ξ i = (2.8) µ i ±iµ i, i = 1,...,m 2 are non-zero eigenalues of matrix A, that is, µ i = (±ω i ) (±ω j ) for i,j = 1,...,n (see [3, Corollary ]). Note that m 2 = (m 2 m)/2. Set D1 = U T D F, D 2 D F is defined in (2.5). Now, ([ ) A D F D T D1 D F = U Â] T 1 D 1 D T 2 D 2 D T 1 D 2 D T U T. (2.9) 2 Taking [ ] I G = D 2 D T 1 (D 1 D T 1 ) 1 I (2.1) 6
7 we obtain Note that we can write A D F D T F = U G 1 [ D1 D T 1 ] G T U T, (2.11) à à =  D 2 (I D T 1 (D 1 D T 1 ) 1 D 1 )D T 2. (2.12) D 2 (I D T 1 (D 1 D T 1 ) 1 D 1 )D T 2 = FF T. Further, we hae to find the inerse of à =  FF T. This is obtained by using the Sherman-Morrison-Woodbury formula ([1, (2.1.4), pg.51.]), that is, ( ) 1 ( FF T ) 1 =  1 +  1 F I F T  1 F F T  1. (2.13) The inerse of I F T  1 Λ = diag F remains to be found. Let ] λ2,,..., λ 2 ([ λ1 λ 1 [ λs λ s ]), ±iλ 1, ±iλ 2,..., ±iλ s are non-anishing finite eigenalues of the problem (A λd)ec(y) =. Since F T  1 F is skew-symmetric, then there exists an orthogonal matrix U S of order 2(r 1)m such that U T S F T  1 FU S = Γ = Λ 1. Using (2.11), (2.13), (2.14) and (2.6) it follows, (2.14) Γ ec(x) = U G T [ 1 2 ] GU T ec(i), (2.15) 1 = (D 1D T 1 ) 1, 2 =  1 +  1 FU S Since Γ is block diagonal we hae ( [ (I Γ) 1 1 λ 2 = diag 1 λ 1 λ λ 1 λ 2 1 Using (2.16) and (2.15) it follows ( V 1 s ec(x) = V V 1 + i=1 7 ],..., I (I Γ) 1 Us T F T  1. 1 λ 2 s + 2 (2.16) ) λ 2 s λ s λ s λ 2. s (2.17) ) λ i (λ i W i Z i ) λ 2 i + ec(i), (2.18) 2
8 matrices V 1, V, V 1, W i and Z i are constructed using (2.1), (2.16) and (2.15). By reshaping ectors in (2.18) back into m m matrices we obtain the solution of equation (2.1) X = Ψ 1 Ψ Ψ 1 + s i=1 λ i (λ i Φ i Υ i ) λ 2 i + 2. (2.19) If one is interested in deriing an optimal damping, then according to (1.11) one has to minimize the function f() = tr(zx()), This gies tr(zx()) = ec(z) T ec(x) = ec(z) T ( A D F D T F) 1ec(I). Tr(ZX) = X 1 X X 1 + s i=1 X 1 = ec(z) T V 1 ec(i), λ i (λ i X i Y i ) λ 2 i + 2, (2.2) X = ec(z) T V ec(i), (2.21) X 1 = ec(z) T V 1 ec(i), and X i = ec(z) T W i ec(i), Y i = ec(z) T Z i ec(i). (2.22) Note that the function Tr(ZX) from (2.2) is a generalization of the function Tr(X) defined in (1.14). Remark 2.1 In the case when Z is a diagonal matrix, the function Tr(ZX) from (2.2) has the following simpler form: Tr(ZX) = X 1 X 1 + s i=1 λ 2 i X i λ 2 i + 2. Finally we derie an explicit formula for the matrix Ψ 1. Let U be that part of U corresponding to the null-space of A. From (2.18) it follows that ec(ψ 1 ) = V 1 ec(i) U (D 1 D T 1 ) 1 U T ec(i). (2.23) It is easily seen that i-th column of the matrix U can be written as ec(o i ), O i is a block-diagonal matrix and O 2i 1 = diag(,...,,o ii,,...,) O 2i = diag(,...,,ôii,,...,) 8
9 O ii and Ôii are orthogonal solutions of Ω i O ij + O ij Ω j = i,j = 1,...,n, that is O ii = 1 2, and Ôii = 1. 2 Let A,B = Tr(A T B) be the usual Frobenius scalar product. Then, we can write the (p,q)-th element of U T D F D T F U as (U ) T (:,p) D F D T F(U ) (:,q) = D pq O p + O p D pq,o q p,q = 1,...,m, D pq = (D pq ) and D pq = D (p,:) D T. (q,:) The orthonormality property O p,o q = δ pq, implies ( U T D F D T F U ) 1 = diag(1/(d)22,1/(d) 22,...,1/(D) mm,1/(d) mm ). Using the fact that U U T ec(i) = ec(i), from (2.23) it follows that ( ) Ψ 1 = diag,,...,,. (2.24) D 22 D 22 D mm D mm After applying a perfect shuffle permutation we hae D Ψ 1 =, (2.25) D D = (diag (D)) 1. The explicitness of the obtained formulas is attractie for possible numerical computation. Our first attempts to perform this task did not succeed due to unexpected complexity problems. We will come back to this issue in our future research. References [1] G.H. Golub and Ch.F. an Loan, Matrix computations, J. Hopkins Uniersity Press, Baltimore,
10 [2] Herbert Federer, Geometric Measure Theory. Springer-Verlag, New York Inc., New York, [3] P. Lancaster and M. Tismenetsky, The Theory of Matrices, Academic Press, New York [4] I. Nakić, Optimal damping of ibrational systems, Ph. D. thesis Fernuniersität, Hagen, 22. [5] N. Truhar, An efficient algorithm for damper optimization for linear ibrating systems using Lyapuno equation, J. Comput. Appl. Math (24), [6] K. Veselić, K. Brabender and K. Delinić, Passie control of linear systems, in: Applied Mathematics and Computation, (M. Rogina et al. Eds.) Dept. of Math. Uni. Zagreb, 21, [7] K. Veselić, On linear ibrational systems with one dimensional damping II, Integral Eq. Operator Th., 13(199),
Journal of Computational and Applied Mathematics. New matrix iterative methods for constraint solutions of the matrix
Journal of Computational and Applied Mathematics 35 (010 76 735 Contents lists aailable at ScienceDirect Journal of Computational and Applied Mathematics journal homepage: www.elseier.com/locate/cam New
More informationDepartment of Mathematics J. J. Strossmayer University of Osijek
Department of Mathematics J. J. Strossmayer University of Osijek Series: Technical Reports Optimal Damping of Selected Eigenfrequencies Using Dimension Reduction Peter Benner, Zoran Tomljanović Ninoslav
More informationAN EFFICIENT APPROXIMATION FOR OPTIMAL DAMPING IN MECHANICAL SYSTEMS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 14, Number 2, Pages 201 217 c 2017 Institute for Scientific Computing and Information AN EFFICIENT APPROXIMATION FOR OPTIMAL DAMPING IN MECHANICAL
More informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016
Prof. Dr. Eleni Chatzi Lecture 4-09. March, 2016 Fundamentals Overview Multiple DOF Systems State-space Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural
More informationReview of Matrices and Vectors 1/45
Reiew of Matrices and Vectors /45 /45 Definition of Vector: A collection of comple or real numbers, generally put in a column [ ] T "! Transpose + + + b a b a b b a a " " " b a b a Definition of Vector
More informationDifferential Geometry of Surfaces
Differential Geometry of urfaces Jordan mith and Carlo équin C Diision, UC Berkeley Introduction These are notes on differential geometry of surfaces ased on reading Greiner et al. n. d.. Differential
More informationPassive Control of Linear Systems
Proceedings of the 1. Conference on Applied Mathematics and Computation Dubrovnik, Croatia, September 13 18, 1999 pp. 39 68 Passive Control of Linear Systems Krešimir Veselić, Katrin Brabender, and Krešimir
More informationCollocated versus non-collocated control [H04Q7]
Collocated versus non-collocated control [H04Q7] Jan Swevers September 2008 0-0 Contents Some concepts of structural dynamics Collocated versus non-collocated control Summary This lecture is based on parts
More information22.2. Applications of Eigenvalues and Eigenvectors. Introduction. Prerequisites. Learning Outcomes
Applications of Eigenvalues and Eigenvectors 22.2 Introduction Many applications of matrices in both engineering and science utilize eigenvalues and, sometimes, eigenvectors. Control theory, vibration
More informationA Geometric Review of Linear Algebra
A Geometric Reiew of Linear Algebra The following is a compact reiew of the primary concepts of linear algebra. The order of presentation is unconentional, with emphasis on geometric intuition rather than
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 information1 :: Mathematical notation
1 :: Mathematical notation x A means x is a member of the set A. A B means the set A is contained in the set B. {a 1,..., a n } means the set hose elements are a 1,..., a n. {x A : P } means the set of
More informationLecture #9: Harmonic Oscillator: Creation and Annihilation Operators
5.61 Fall, 013 Lecture #9 Page 1 Last time Lecture #9: Harmonic Oscillator: Creation and Annihilation Operators 1/ Simplified Schrödinger equation: ξ=α x, α =(kμ) 1/ n E +ξ ψ=0 (dimensionless) ξ nω reduced
More information0 a 3 a 2 a 3 0 a 1 a 2 a 1 0
Chapter Flow kinematics Vector and tensor formulae This introductory section presents a brief account of different definitions of ector and tensor analysis that will be used in the following chapters.
More informationON THE HÖLDER CONTINUITY OF MATRIX FUNCTIONS FOR NORMAL MATRICES
Volume 10 (2009), Issue 4, Article 91, 5 pp. ON THE HÖLDER CONTINUITY O MATRIX UNCTIONS OR NORMAL MATRICES THOMAS P. WIHLER MATHEMATICS INSTITUTE UNIVERSITY O BERN SIDLERSTRASSE 5, CH-3012 BERN SWITZERLAND.
More informationMA 266 Review Topics - Exam # 2 (updated)
MA 66 Reiew Topics - Exam # updated Spring First Order Differential Equations Separable, st Order Linear, Homogeneous, Exact Second Order Linear Homogeneous with Equations Constant Coefficients The differential
More informationEIGENVALUES AND EIGENVECTORS
CHAPTER 6 EIGENVALUES AND EIGENVECTORS SECTION 6. INTRODUCTION TO EIGENVALUES In each of Problems we first list the characteristic polynomial p( λ) A λi of the gien matrix A and then the roots of p( λ
More informationTwo Results About The Matrix Exponential
Two Results About The Matrix Exponential Hongguo Xu Abstract Two results about the matrix exponential are given. One is to characterize the matrices A which satisfy e A e AH = e AH e A, another is about
More informationMATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators.
MATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators. Adjoint operator and adjoint matrix Given a linear operator L on an inner product space V, the adjoint of L is a transformation
More informationA Geometric Review of Linear Algebra
A Geometric Reiew of Linear Algebra The following is a compact reiew of the primary concepts of linear algebra. I assume the reader is familiar with basic (i.e., high school) algebra and trigonometry.
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 Riesz-Fischer sequences and lower frame bounds
On Riesz-Fischer sequences and lower frame bounds P. Casazza, O. Christensen, S. Li, A. Lindner Abstract We investigate the consequences of the lower frame condition and the lower Riesz basis condition
More informationNOTES ON THE REGULAR E-OPTIMAL SPRING BALANCE WEIGHING DESIGNS WITH CORRE- LATED ERRORS
REVSTAT Statistical Journal Volume 3, Number 2, June 205, 9 29 NOTES ON THE REGULAR E-OPTIMAL SPRING BALANCE WEIGHING DESIGNS WITH CORRE- LATED ERRORS Authors: Bronis law Ceranka Department of Mathematical
More informationOctober 25, 2013 INNER PRODUCT SPACES
October 25, 2013 INNER PRODUCT SPACES RODICA D. COSTIN Contents 1. Inner product 2 1.1. Inner product 2 1.2. Inner product spaces 4 2. Orthogonal bases 5 2.1. Existence of an orthogonal basis 7 2.2. Orthogonal
More informationReview of similarity transformation and Singular Value Decomposition
Review of similarity transformation and Singular Value Decomposition Nasser M Abbasi Applied Mathematics Department, California State University, Fullerton July 8 7 page compiled on June 9, 5 at 9:5pm
More informationReview of Some Concepts from Linear Algebra: Part 2
Review of Some Concepts from Linear Algebra: Part 2 Department of Mathematics Boise State University January 16, 2019 Math 566 Linear Algebra Review: Part 2 January 16, 2019 1 / 22 Vector spaces A set
More informationTo string together six theorems of physics by Pythagoras theorem
To string together six theorems of physics by Pythagoras theorem H. Y. Cui Department of Applied Physics Beijing Uniersity of Aeronautics and Astronautics Beijing, 00083, China ( May, 8, 2002 ) Abstract
More informationNotes on Linear Minimum Mean Square Error Estimators
Notes on Linear Minimum Mean Square Error Estimators Ça gatay Candan January, 0 Abstract Some connections between linear minimum mean square error estimators, maximum output SNR filters and the least square
More informationWould you risk your live driving drunk? Intro
Martha Casquete Would you risk your lie driing drunk? Intro Motion Position and displacement Aerage elocity and aerage speed Instantaneous elocity and speed Acceleration Constant acceleration: A special
More informationWeek Quadratic forms. Principal axes theorem. Text reference: this material corresponds to parts of sections 5.5, 8.2,
Math 051 W008 Margo Kondratieva Week 10-11 Quadratic forms Principal axes theorem Text reference: this material corresponds to parts of sections 55, 8, 83 89 Section 41 Motivation and introduction Consider
More informationInner product spaces. Layers of structure:
Inner product spaces Layers of structure: vector space normed linear space inner product space The abstract definition of an inner product, which we will see very shortly, is simple (and by itself is pretty
More information1 Quasi-definite matrix
1 Quasi-definite matrix The matrix H is a quasi-definite matrix, if there exists a permutation matrix P such that H qd P T H11 H HP = 1 H1, 1) H where H 11 and H + H1H 11 H 1 are positive definite. This
More informationMath 520 Exam 2 Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008
Math 520 Exam 2 Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008 Exam 2 will be held on Tuesday, April 8, 7-8pm in 117 MacMillan What will be covered The exam will cover material from the lectures
More informationResidual migration in VTI media using anisotropy continuation
Stanford Exploration Project, Report SERGEY, Noember 9, 2000, pages 671?? Residual migration in VTI media using anisotropy continuation Tariq Alkhalifah Sergey Fomel 1 ABSTRACT We introduce anisotropy
More informationDynamics of Structures
Dynamics of Structures Elements of structural dynamics Roberto Tomasi 11.05.2017 Roberto Tomasi Dynamics of Structures 11.05.2017 1 / 22 Overview 1 SDOF system SDOF system Equation of motion Response spectrum
More informationarxiv: v1 [math.gt] 2 Nov 2010
CONSTRUCTING UNIVERSAL ABELIAN COVERS OF GRAPH MANIFOLDS HELGE MØLLER PEDERSEN arxi:101105551 [mathgt] 2 No 2010 Abstract To a rational homology sphere graph manifold one can associate a weighted tree
More informationThroughout these notes we assume V, W are finite dimensional inner product spaces over C.
Math 342 - Linear Algebra II Notes Throughout these notes we assume V, W are finite dimensional inner product spaces over C 1 Upper Triangular Representation Proposition: Let T L(V ) There exists an orthonormal
More informationProblem Set 1. Homeworks will graded based on content and clarity. Please show your work clearly for full credit.
CSE 151: Introduction to Machine Learning Winter 2017 Problem Set 1 Instructor: Kamalika Chaudhuri Due on: Jan 28 Instructions This is a 40 point homework Homeworks will graded based on content and clarity
More information2. Review of Linear Algebra
2. Review of Linear Algebra ECE 83, Spring 217 In this course we will represent signals as vectors and operators (e.g., filters, transforms, etc) as matrices. This lecture reviews basic concepts from linear
More informationCases of integrability corresponding to the motion of a pendulum on the two-dimensional plane
Cases of integrability corresponding to the motion of a pendulum on the two-dimensional plane MAXIM V. SHAMOLIN Lomonoso Moscow State Uniersity Institute of Mechanics Michurinskii Ae.,, 99 Moscow RUSSIAN
More informationMath 144 Activity #9 Introduction to Vectors
144 p 1 Math 144 ctiity #9 Introduction to Vectors Often times you hear people use the words speed and elocity. Is there a difference between the two? If so, what is the difference? Discuss this with your
More informationLINEAR ALGEBRA BOOT CAMP WEEK 4: THE SPECTRAL THEOREM
LINEAR ALGEBRA BOOT CAMP WEEK 4: THE SPECTRAL THEOREM Unless otherwise stated, all vector spaces in this worksheet are finite dimensional and the scalar field F is R or C. Definition 1. A linear operator
More informationMATH 425-Spring 2010 HOMEWORK ASSIGNMENTS
MATH 425-Spring 2010 HOMEWORK ASSIGNMENTS Instructor: Shmuel Friedland Department of Mathematics, Statistics and Computer Science email: friedlan@uic.edu Last update April 18, 2010 1 HOMEWORK ASSIGNMENT
More informationElements of Positive Definite Kernel and Reproducing Kernel Hilbert Space
Elements of Positive Definite Kernel and Reproducing Kernel Hilbert Space Statistical Inference with Reproducing Kernel Hilbert Space Kenji Fukumizu Institute of Statistical Mathematics, ROIS Department
More informationRecursion Systems and Recursion Operators for the Soliton Equations Related to Rational Linear Problem with Reductions
GMV The s Systems and for the Soliton Equations Related to Rational Linear Problem with Reductions Department of Mathematics & Applied Mathematics University of Cape Town XIV th International Conference
More informationCHAPTER 3 : VECTORS. Definition 3.1 A vector is a quantity that has both magnitude and direction.
EQT 101-Engineering Mathematics I Teaching Module CHAPTER 3 : VECTORS 3.1 Introduction Definition 3.1 A ector is a quantity that has both magnitude and direction. A ector is often represented by an arrow
More informationQUASI-UNIFORMLY POSITIVE OPERATORS IN KREIN SPACE. Denitizable operators in Krein spaces have spectral properties similar to those
QUASI-UNIFORMLY POSITIVE OPERATORS IN KREIN SPACE BRANKO CURGUS and BRANKO NAJMAN Denitizable operators in Krein spaces have spectral properties similar to those of selfadjoint operators in Hilbert spaces.
More informationStructural Consensus Controllability of Singular Multi-agent Linear Dynamic Systems
Structural Consensus Controllability of Singular Multi-agent Linear Dynamic Systems M. ISAL GARCÍA-PLANAS Universitat Politècnica de Catalunya Departament de Matèmatiques Minería 1, sc. C, 1-3, 08038 arcelona
More informationLecture 4. CP and KSVD Representations. Charles F. Van Loan
Structured Matrix Computations from Structured Tensors Lecture 4. CP and KSVD Representations Charles F. Van Loan Cornell University CIME-EMS Summer School June 22-26, 2015 Cetraro, Italy Structured Matrix
More informationON THE SINGULAR DECOMPOSITION OF MATRICES
An. Şt. Univ. Ovidius Constanţa Vol. 8, 00, 55 6 ON THE SINGULAR DECOMPOSITION OF MATRICES Alina PETRESCU-NIŢǍ Abstract This paper is an original presentation of the algorithm of the singular decomposition
More informationLinear Algebra Formulas. Ben Lee
Linear Algebra Formulas Ben Lee January 27, 2016 Definitions and Terms Diagonal: Diagonal of matrix A is a collection of entries A ij where i = j. Diagonal Matrix: A matrix (usually square), where entries
More informationDynamics ( 동역학 ) Ch.2 Motion of Translating Bodies (2.1 & 2.2)
Dynamics ( 동역학 ) Ch. Motion of Translating Bodies (. &.) Motion of Translating Bodies This chapter is usually referred to as Kinematics of Particles. Particles: In dynamics, a particle is a body without
More informationSection 6: PRISMATIC BEAMS. Beam Theory
Beam Theory There are two types of beam theory aailable to craft beam element formulations from. They are Bernoulli-Euler beam theory Timoshenko beam theory One learns the details of Bernoulli-Euler beam
More informationECE 275A Homework #3 Solutions
ECE 75A Homework #3 Solutions. Proof of (a). Obviously Ax = 0 y, Ax = 0 for all y. To show sufficiency, note that if y, Ax = 0 for all y, then it must certainly be true for the particular value of y =
More informationReview problems for MA 54, Fall 2004.
Review problems for MA 54, Fall 2004. Below are the review problems for the final. They are mostly homework problems, or very similar. If you are comfortable doing these problems, you should be fine on
More informationLecture 2: Linear Algebra Review
EE 227A: Convex Optimization and Applications January 19 Lecture 2: Linear Algebra Review Lecturer: Mert Pilanci Reading assignment: Appendix C of BV. Sections 2-6 of the web textbook 1 2.1 Vectors 2.1.1
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 informationSeminar 6: COUPLED HARMONIC OSCILLATORS
Seminar 6: COUPLED HARMONIC OSCILLATORS 1. Lagrangian Equations of Motion Let consider a system consisting of two harmonic oscillators that are coupled together. As a model, we will use two particles attached
More informationDuke University, Department of Electrical and Computer Engineering Optimization for Scientists and Engineers c Alex Bronstein, 2014
Duke University, Department of Electrical and Computer Engineering Optimization for Scientists and Engineers c Alex Bronstein, 2014 Linear Algebra A Brief Reminder Purpose. The purpose of this document
More informationNumerical Linear Algebra
University of Alabama at Birmingham Department of Mathematics Numerical Linear Algebra Lecture Notes for MA 660 (1997 2014) Dr Nikolai Chernov April 2014 Chapter 0 Review of Linear Algebra 0.1 Matrices
More informationCSL361 Problem set 4: Basic linear algebra
CSL361 Problem set 4: Basic linear algebra February 21, 2017 [Note:] If the numerical matrix computations turn out to be tedious, you may use the function rref in Matlab. 1 Row-reduced echelon matrices
More informationMatrices and Vectors. Definition of Matrix. An MxN matrix A is a two-dimensional array of numbers A =
30 MATHEMATICS REVIEW G A.1.1 Matrices and Vectors Definition of Matrix. An MxN matrix A is a two-dimensional array of numbers A = a 11 a 12... a 1N a 21 a 22... a 2N...... a M1 a M2... a MN A matrix can
More information1 Linear Algebra Problems
Linear Algebra Problems. Let A be the conjugate transpose of the complex matrix A; i.e., A = A t : A is said to be Hermitian if A = A; real symmetric if A is real and A t = A; skew-hermitian if A = A and
More informationTHE CAUCHY PROBLEM FOR ONE-DIMENSIONAL FLOW OF A COMPRESSIBLE VISCOUS FLUID: STABILIZATION OF THE SOLUTION
GLASNIK MATEMATIČKI Vol. 4666, 5 3 THE CAUCHY POBLEM FO ONE-DIMENSIONAL FLOW OF A COMPESSIBLE VISCOUS FLUID: STABILIZATION OF THE SOLUTION Nermina Mujakoić and Ian Dražić Uniersity of ijeka, Croatia Abstract.
More informationUniqueness of the Solutions of Some Completion Problems
Uniqueness of the Solutions of Some Completion Problems Chi-Kwong Li and Tom Milligan Abstract We determine the conditions for uniqueness of the solutions of several completion problems including the positive
More information11 a 12 a 21 a 11 a 22 a 12 a 21. (C.11) A = The determinant of a product of two matrices is given by AB = A B 1 1 = (C.13) and similarly.
C PROPERTIES OF MATRICES 697 to whether the permutation i 1 i 2 i N is even or odd, respectively Note that I =1 Thus, for a 2 2 matrix, the determinant takes the form A = a 11 a 12 = a a 21 a 11 a 22 a
More informationLECTURE 2: CROSS PRODUCTS, MULTILINEARITY, AND AREAS OF PARALLELOGRAMS
LECTURE : CROSS PRODUCTS, MULTILINEARITY, AND AREAS OF PARALLELOGRAMS MA1111: LINEAR ALGEBRA I, MICHAELMAS 016 1. Finishing up dot products Last time we stated the following theorem, for which I owe you
More informationMotion in Two and Three Dimensions
PH 1-1D Spring 013 Motion in Two and Three Dimensions Lectures 5,6,7 Chapter 4 (Halliday/Resnick/Walker, Fundamentals of Physics 9 th edition) 1 Chapter 4 Motion in Two and Three Dimensions In this chapter
More information6.1.1 Angle between Two Lines Intersection of Two lines Shortest Distance from a Point to a Line
CHAPTER 6 : VECTORS 6. Lines in Space 6.. Angle between Two Lines 6.. Intersection of Two lines 6..3 Shortest Distance from a Point to a Line 6. Planes in Space 6.. Intersection of Two Planes 6.. Angle
More informationSQUARE ROOTS OF 2x2 MATRICES 1. Sam Northshield SUNY-Plattsburgh
SQUARE ROOTS OF x MATRICES Sam Northshield SUNY-Plattsburgh INTRODUCTION A B What is the square root of a matrix such as? It is not, in general, A B C D C D This is easy to see since the upper left entry
More informationSTAT200C: Review of Linear Algebra
Stat200C Instructor: Zhaoxia Yu STAT200C: Review of Linear Algebra 1 Review of Linear Algebra 1.1 Vector Spaces, Rank, Trace, and Linear Equations 1.1.1 Rank and Vector Spaces Definition A vector whose
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 informationAN INVERSE EIGENVALUE PROBLEM AND AN ASSOCIATED APPROXIMATION PROBLEM FOR GENERALIZED K-CENTROHERMITIAN MATRICES
AN INVERSE EIGENVALUE PROBLEM AND AN ASSOCIATED APPROXIMATION PROBLEM FOR GENERALIZED K-CENTROHERMITIAN MATRICES ZHONGYUN LIU AND HEIKE FAßBENDER Abstract: A partially described inverse eigenvalue problem
More informationFinding eigenvalues for matrices acting on subspaces
Finding eigenvalues for matrices acting on subspaces Jakeniah Christiansen Department of Mathematics and Statistics Calvin College Grand Rapids, MI 49546 Faculty advisor: Prof Todd Kapitula Department
More informationQ1 Q2 Q3 Q4 Tot Letr Xtra
Mathematics 54.1 Final Exam, 12 May 2011 180 minutes, 90 points NAME: ID: GSI: INSTRUCTIONS: You must justify your answers, except when told otherwise. All the work for a question should be on the respective
More informationMath 443 Differential Geometry Spring Handout 3: Bilinear and Quadratic Forms This handout should be read just before Chapter 4 of the textbook.
Math 443 Differential Geometry Spring 2013 Handout 3: Bilinear and Quadratic Forms This handout should be read just before Chapter 4 of the textbook. Endomorphisms of a Vector Space This handout discusses
More informationEconometrics II - EXAM Outline Solutions All questions have 25pts Answer each question in separate sheets
Econometrics II - EXAM Outline Solutions All questions hae 5pts Answer each question in separate sheets. Consider the two linear simultaneous equations G with two exogeneous ariables K, y γ + y γ + x δ
More informationMathematical Methods wk 2: Linear Operators
John Magorrian, magog@thphysoxacuk These are work-in-progress notes for the second-year course on mathematical methods The most up-to-date version is available from http://www-thphysphysicsoxacuk/people/johnmagorrian/mm
More informationMIT Final Exam Solutions, Spring 2017
MIT 8.6 Final Exam Solutions, Spring 7 Problem : For some real matrix A, the following vectors form a basis for its column space and null space: C(A) = span,, N(A) = span,,. (a) What is the size m n of
More informationMATH 431: FIRST MIDTERM. Thursday, October 3, 2013.
MATH 431: FIRST MIDTERM Thursday, October 3, 213. (1) An inner product on the space of matrices. Let V be the vector space of 2 2 real matrices (that is, the algebra Mat 2 (R), but without the mulitiplicative
More informationKernel Method: Data Analysis with Positive Definite Kernels
Kernel Method: Data Analysis with Positive Definite Kernels 2. Positive Definite Kernel and Reproducing Kernel Hilbert Space Kenji Fukumizu The Institute of Statistical Mathematics. Graduate University
More informationA spectral Turán theorem
A spectral Turán theorem Fan Chung Abstract If all nonzero eigenalues of the (normalized) Laplacian of a graph G are close to, then G is t-turán in the sense that any subgraph of G containing no K t+ contains
More informationDamped oscillations of linear systems a mathematical introduction
Krešimir Veselić Damped oscillations of linear systems a mathematical introduction Wednesday 30 th March, 2011 Springer To my wife. viii Foreword The theory of linear damped oscillations has been studied
More information1 The space of linear transformations from R n to R m :
Math 540 Spring 20 Notes #4 Higher deriaties, Taylor s theorem The space of linear transformations from R n to R m We hae discssed linear transformations mapping R n to R m We can add sch linear transformations
More information2 trees T 1,T 2,...,T
Volume 3, pp 95-305, March 0 http://mathtechnionacil/iic/ela ORDERING TREES BY THE MINIMAL ENTRIES OF THEIR DOUBLY STOCHASTIC GRAPH MATRICES SHUCHAO LI AND QIN ZHAO Abstract Gien an n-ertex graph G, the
More informationbe ye transformed by the renewing of your mind Romans 12:2
Lecture 12: Coordinate Free Formulas for Affine and rojectie Transformations be ye transformed by the reing of your mind Romans 12:2 1. Transformations for 3-Dimensional Computer Graphics Computer Graphics
More informationarxiv: v1 [math.na] 1 Sep 2018
On the perturbation of an L -orthogonal projection Xuefeng Xu arxiv:18090000v1 [mathna] 1 Sep 018 September 5 018 Abstract The L -orthogonal projection is an important mathematical tool in scientific computing
More informationStat 206: Linear algebra
Stat 206: Linear algebra James Johndrow (adapted from Iain Johnstone s notes) 2016-11-02 Vectors We have already been working with vectors, but let s review a few more concepts. The inner product of two
More informationMath 425 Lecture 1: Vectors in R 3, R n
Math 425 Lecture 1: Vectors in R 3, R n Motiating Questions, Problems 1. Find the coordinates of a regular tetrahedron with center at the origin and sides of length 1. 2. What is the angle between the
More informationMATH 583A REVIEW SESSION #1
MATH 583A REVIEW SESSION #1 BOJAN DURICKOVIC 1. Vector Spaces Very quick review of the basic linear algebra concepts (see any linear algebra textbook): (finite dimensional) vector space (or linear space),
More informationOVERVIEW OF TATE S THESIS
OVERVIEW O TATE S THESIS ALEX MINE Abstract. This paper gies an oeriew of the main idea of John Tate s 1950 PhD thesis. I will explain the methods he used without going into too much technical detail.
More informationMath 102, Winter Final Exam Review. Chapter 1. Matrices and Gaussian Elimination
Math 0, Winter 07 Final Exam Review Chapter. Matrices and Gaussian Elimination { x + x =,. Different forms of a system of linear equations. Example: The x + 4x = 4. [ ] [ ] [ ] vector form (or the column
More informationThe perturbed Riemann problem for the chromatography system of Langmuir isotherm with one inert component
Aailable online at www.tjnsa.com J. Nonlinear Sci. Appl. 9 2016, 5382 5397 Research Article The perturbed Riemann problem for the chromatography system of Langmuir isotherm with one inert component Pengpeng
More informationMath 1060 Linear Algebra Homework Exercises 1 1. Find the complete solutions (if any!) to each of the following systems of simultaneous equations:
Homework Exercises 1 1 Find the complete solutions (if any!) to each of the following systems of simultaneous equations: (i) x 4y + 3z = 2 3x 11y + 13z = 3 2x 9y + 2z = 7 x 2y + 6z = 2 (ii) x 4y + 3z =
More informationApproximations - the method of least squares (1)
Approximations - the method of least squares () In many applications, we have to consider the following problem: Suppose that for some y, the equation Ax = y has no solutions It could be that this is an
More informationNumerically Stable Cointegration Analysis
Numerically Stable Cointegration Analysis Jurgen A. Doornik Nuffield College, University of Oxford, Oxford OX1 1NF R.J. O Brien Department of Economics University of Southampton November 3, 2001 Abstract
More informationClifford Algebras and Spin Groups
Clifford Algebras and Spin Groups Math G4344, Spring 2012 We ll now turn from the general theory to examine a specific class class of groups: the orthogonal groups. Recall that O(n, R) is the group of
More information2.3. VECTOR SPACES 25
2.3. VECTOR SPACES 25 2.3 Vector Spaces MATH 294 FALL 982 PRELIM # 3a 2.3. Let C[, ] denote the space of continuous functions defined on the interval [,] (i.e. f(x) is a member of C[, ] if f(x) is continuous
More informationMotion in Two and Three Dimensions
PH 1-A Fall 014 Motion in Two and Three Dimensions Lectures 4,5 Chapter 4 (Halliday/Resnick/Walker, Fundamentals of Physics 9 th edition) 1 Chapter 4 Motion in Two and Three Dimensions In this chapter
More information6.3 Vectors in a Plane
6.3 Vectors in a Plane Plan: Represent ectors as directed line segments. Write the component form of ectors. Perform basic ector operations and represent ectors graphically. Find the direction angles of
More information