Pseudospectra and Nonnormal Dynamical Systems
|
|
- Anabel Newman
- 5 years ago
- Views:
Transcription
1 Pseudospectra and Nonnormal Dynamical Systems Mark Embree and Russell Carden Computational and Applied Mathematics Rice University Houston, Texas ELGERSBURG MARCH 1
2 Overview of the Course These lectures describe modern tools for the spectral analysis of dynamical systems We shall cover a mix of theory, computation, and applications Lecture 1: Introduction to Nonnormality and Pseudospectra Lecture : Functions of Matrices Lecture : Toeplitz Matrices and Model Reduction Lecture : Model Reduction, Numerical Algorithms, Differential Operators Lecture : Discretization, Extensions, Applications
3 Outline for Today Lecture : Functions of Matrices and Model Reduction Toeplitz matrices Model Reduction: Balanced Truncation Nonnormality and Lyapunov Equations Model Reduction: Moment Matching
4 (a) Toeplitz Matrices
5 Pseudospectra Recall the example that began our investigation of pseudospectra yesterday Example Compute eigenvalues of three similar 1 1 matrices using MATLAB s eig / 1/ 1/ 1/
6 Toeplitz Matrices Consider the pseudospectra of the 1 1 matrix in the middle of the last slide, A = tridiag(,, 1/) A is diagonalizable (it has distinct eigenvalues), but Bauer Fike is useless here: κ(v) =
7 Jordan Blocks We ve already analyzed pseudospectra of Jordan blocks near λ for small ε > Here we want to investigate the entire pseudospectrum for larger ε S n = 1 1 Near the eigenvalue, the resolvent norm grows with dimension n; outside the unit disk, the resolvent norm does not seem to get big We would like to prove this
8 Jordan Blocks We ve already analyzed pseudospectra of Jordan blocks near λ for small ε > Here we want to investigate the entire pseudospectrum for larger ε S n = 1 1 Near the eigenvalue, the resolvent norm grows with dimension n; outside the unit disk, the resolvent norm does not seem to get big We would like to prove this
9 Jordan Blocks We ve already analyzed pseudospectra of Jordan blocks near λ for small ε > Here we want to investigate the entire pseudospectrum for larger ε S n = 1 1 Near the eigenvalue, the resolvent norm grows with dimension n; outside the unit disk, the resolvent norm does not seem to get big We would like to prove this
10 Jordan Blocks/Shift Operator Consider the generalization of the Jordan block to the domain l (N) = {(x 1, x, ) : X x j < } j=1 The shift operator S on l (N) is defined as S(x 1, x, ) = (x, x, )
11 Jordan Blocks/Shift Operator Consider the generalization of the Jordan block to the domain l (N) = {(x 1, x, ) : X x j < } j=1 The shift operator S on l (N) is defined as S(x 1, x, ) = (x, x, ) In particular, S(1, z, z, ) = (z, z, z, ) = z(1, z, z, ) So if (1, z, z, ) l (N), then z σ(s)
12 Jordan Blocks/Shift Operator Consider the generalization of the Jordan block to the domain l (N) = {(x 1, x, ) : X x j < } j=1 The shift operator S on l (N) is defined as S(x 1, x, ) = (x, x, ) In particular, S(1, z, z, ) = (z, z, z, ) = z(1, z, z, ) So if (1, z, z, ) l (N), then z σ(s) If z < 1, then So, X z j 1 1 = 1 z < j=1 {z C : z < 1} σ(s)
13 Jordan Blocks/Shift Operator S(x 1, x, ) = (x, x, ) We have seen that {z C : z < 1} σ(s) Observe that and so The spectrum is closed, so S = sup Sx = 1, x =1 σ(s) {z C : z 1} σ(a) = {z C : z 1} For any finite dimensional n n Jordan block S n, σ(s n) = {}
14 Jordan Blocks/Shift Operator S(x 1, x, ) = (x, x, ) We have seen that {z C : z < 1} σ(s) Observe that and so The spectrum is closed, so S = sup Sx = 1, x =1 σ(s) {z C : z 1} σ(a) = {z C : z 1} For any finite dimensional n n Jordan block S n, σ(s n) = {} So the S n S strongly, but there is a discontinuity in the spectrum: σ(s n) σ(s)
15 Pseudospectra of Jordan Blocks Pseudospectra resolve this unpleasant discontinuity Recall the eigenvectors (1, z, z, ) for S Truncate this vector to length n, and apply it to S n: z z n z n 1 = z z z n 1
16 Pseudospectra of Jordan Blocks Pseudospectra resolve this unpleasant discontinuity Recall the eigenvectors (1, z, z, ) for S Truncate this vector to length n, and apply it to S n: z z n z n 1 = z z z n 1 = z 1 z z n z n 1 z n
17 Pseudospectra of Jordan Blocks Pseudospectra resolve this unpleasant discontinuity Recall the eigenvectors (1, z, z, ) for S Truncate this vector to length n, and apply it to S n: z z n z n 1 = z z z n 1 = z 1 z z n z n 1 z n Hence, S nx z x = z n, so for all ε > z n x = z n p 1 z n p 1 z, z σ ε(s n) We conclude that for fixed z < 1, the resolvent norm (z S n) 1 grows exponentially with n
18 Upper Triangular Toeplitz Matrices Consider an upper triangular Toeplitz matrix giving the matrix with constant diagonals containing the Laurent coefficients: a a 1 a a A n = a1 a C n n a a 1 a Definition (Symbol, Symbol Curve) Toeplitz matrices are described by their symbol a with Taylor expansion X a(z) = a k z k Call the image of the unit circle T under a the symbol curve, a(t) k=
19 Pseudospectra of Upper Triangular Toeplitz Matrices Apply the same approximate eigenvector we used for the Jordan block: a a 1 a a n 1 a a1 a a a 1 a 1 z z n z n 1 = P n 1 k= a k z k P n 1 k= a kz k+1 a z n + a 1z n 1 a z n 1 If the matrix has fixed bandwidth b n, (ie, a k = for k > b), then a a b a ab a a 1 z z n z n 1 = P b k= a k z k P b k= a kz k+1 P b 1 k= a kz k+n b a z n 1
20 Pseudospectra of Upper Triangular Toeplitz Matrices a a b a ab a a 1 z z n z n 1 = P b k= a k z k P b k= a kz k+1 P b 1 k= a kz k+n b a z n 1 = bx k= a k z k 1 z z n z n 1 a b z n P b k=1 a kz n 1+k
21 Pseudospectra of Upper Triangular Toeplitz Matrices a a b a ab a a 1 z z n z n 1 = P b k= a k z k P b k= a kz k+1 P b 1 k= a kz k+n b a z n 1 = bx k= a k z k 1 z z n z n 1 a b z n P b k=1 a kz n 1+k Hence A nx a(z)x gets very small in size as n In fact, this reveals the spectrum of the infinite Toeplitz operator on l (N)
22 Spectrum of Toeplitz Operators on l (Z) For the Toeplitz operator (semi-infinite matrix) A with the same symbol: a a b a ab a 1 z z = P b k= a k z k P b k= a kz k+1 P b k= a kz k+ = a(z) 1 z z
23 Spectrum of Toeplitz Operators on l (Z) For the Toeplitz operator (semi-infinite matrix) A with the same symbol: a a b a ab a 1 z z = P b k= a k z k P b k= a kz k+1 P b k= a kz k+ = a(z) 1 z z Thus a(z) σ(a) for all z < 1 In fact, one can show that for this banded upper triangular symbol, σ(a ) = {a(z) : z 1}
24 Spectrum of Toeplitz Operators on l (Z) For the Toeplitz operator (semi-infinite matrix) A with the same symbol: a a b a ab a 1 z z = P b k= a k z k P b k= a kz k+1 P b k= a kz k+ = a(z) 1 z z Thus a(z) σ(a) for all z < 1 In fact, one can show that for this banded upper triangular symbol, σ(a ) = {a(z) : z 1} The calculation on the last slide guarantees that for any ε >, there exists N > such that if n > N, σ(a ) σ ε(a n)
25 Symbols of Upper Triangular Toeplitz Matrices Symbol curves for banded upper triangular Toeplitz matrices a(z) = z + z a(z) = z + z + z + z a(z) = z + z z a(z) = z + z a(z) = z + z
26 General Toeplitz Matrices More generally, the dense Toeplitz matrix a a 1 a a 1 a A n = a a1 a a 1 a a 1 a a 1 a C n n follows the same terminology Definition (Symbol, Symbol Curve) Toeplitz matrices are described by their symbol a with Laurent expansion X a(z) = a k z k k= Call the image of the unit circle T under a the symbol curve, a(t)
27 Spectrum of a General Toeplitz Operator Theorem (Spectrum of a Toeplitz Operator) Suppose the Toeplitz operator A : l (N) l (N) has a symbol that is continuous Then σ(a ) = a(t ) {all points a(t ) encloses with nonzero winding number} Due variously to: Wintner; Gohberg; Krein; Calderón, Spitzer, & Widom See Albrecht Böttcher and colleagues for many more details
28 Spectrum of a General Toeplitz Operator Theorem (Spectrum of a Toeplitz Operator) Suppose the Toeplitz operator A : l (N) l (N) has a symbol that is continuous Then σ(a ) = a(t ) {all points a(t ) encloses with nonzero winding number} Due variously to: Wintner; Gohberg; Krein; Calderón, Spitzer, & Widom See Albrecht Böttcher and colleagues for many more details a(t)
29 Spectrum of a General Toeplitz Operator Theorem (Spectrum of a Toeplitz Operator) Suppose the Toeplitz operator A : l (N) l (N) has a symbol that is continuous Then σ(a ) = a(t ) {all points a(t ) encloses with nonzero winding number} Due variously to: Wintner; Gohberg; Krein; Calderón, Spitzer, & Widom See Albrecht Böttcher and colleagues for many more details a(t) σ(a )
30 Spectrum of a General Toeplitz Matrix What can be said of the eigenvalues of a finite-dimensional Toeplitz matrix? Theorem (Limiting Spectrum of Finite Toeplitz Matrices) Consider the family of banded Toeplitz matrices {A n} n N with upper bandwidth b and lower bandwidth d For any fixed λ C, label the roots ζ 1,, ζ b+d of the polynomial z d (a(z) λ) by increasing modulus If ζ d = ζ d+1, then λ lim n σ(an) This result, proved by [Schmidt & Spitzer, 19], shows that in general: lim σ(an) σ(a ) n
31 Spectrum of a General Toeplitz Matrix What can be said of the eigenvalues of a finite-dimensional Toeplitz matrix? Theorem (Limiting Spectrum of Finite Toeplitz Matrices) Consider the family of banded Toeplitz matrices {A n} n N with upper bandwidth b and lower bandwidth d For any fixed λ C, label the roots ζ 1,, ζ b+d of the polynomial z d (a(z) λ) by increasing modulus If ζ d = ζ d+1, then λ lim n σ(an) This result, proved by [Schmidt & Spitzer, 19], shows that in general: lim σ(an) σ(a ) n a(t) lim σ(an) n
32 Spectrum of a General Toeplitz Matrix: Example a(t) lim σ(an) n σ(a )
33 Spectrum of a General Toeplitz Matrix: Example a(t) lim σ(an) n σ(a 1)
34 Spectrum of a General Toeplitz Matrix: Example a(t) lim σ(an) n σ(a )
35 Spectrum of a General Toeplitz Matrix: Example a(t) lim σ(an) n σ(a )
36 Spectrum of a General Toeplitz Matrix: Example a(t) lim σ(an) n σ(a 8)
37 Pseudospectra of Toeplitz Matrices Theorem (Landau; Reichel and Trefethen; Böttcher) Let a(z) = P M k= M a kz k be the symbol of a banded Toeplitz operator The pseudospectra of A n converge to the pseudospectra of the Toeplitz operator on l (N) as n Let z C have nonzero winding number wrt the symbol curve a(t ) (z A n) 1 grows exponentially in n For all ε >, z σ ε(a n) for all n sufficiently large n = 1 n = n = n = 18
38 Pseudospectra of Toeplitz Matrices symbol curve σ ε(a 1)
39 Pseudospectra of Toeplitz Matrices symbol curve σ ε(a )
40 Hermitian Toeplitz Matrices We have seen large pseudospectra arise for generic Toeplitz matrices But what about Hermitian Toeplitz matrices? For example, the matrix has symbol a(z) = z 1 + z, so a(e i θ ) = e i θ + e i θ = cos(θ) [, ] R
41 Hermitian Toeplitz Matrices We have seen large pseudospectra arise for generic Toeplitz matrices But what about Hermitian Toeplitz matrices? For example, the matrix has symbol a(z) = z 1 + z, so a(e i θ ) = e i θ + e i θ = cos(θ) [, ] R In general, Hermitian symbols have a k = a k, so a(e i θ ) = a + X 1 k= k=1 a k e i k θ + X a k e i k θ k=1 X X X = a + a k e i k θ + a k e i k θ = a + Re(a k e i k θ ) R As n, eigenvalues of A n distribute according to Szegő s theorem k=1 k=1
42 Tridiagonal Toeplitz Matrices Consider the case of a tridiagonal Toeplitz matrix, β γ α β A n = γ Cn n α β with symbol a(z) = α z + β + γz The symbol curve is the ellipse a(t) = {z T : α z + β + γz}
43 Eigenvectors of Tridiagonal Toeplitz Matrices A = tridiag(1,, 1), n = Normal matrix: orthogonal eigenvectors
44 Eigenvectors of Tridiagonal Toeplitz Matrices A = tridiag(,, 1/), n = Nonnormal matrix: non-orthogonal eigenvectors
45 Circulant Matrices and Laurent Operators In contrast, consider the circulant matrix a a 1 a a 1 a 1 a a C n = a1 a a 1 a a 1 a 1 a a 1 a C n n
46 Circulant Matrices and Laurent Operators In contrast, consider the circulant matrix a a 1 a a 1 a 1 a a C n = a1 a a 1 a a 1 a 1 a a 1 a C n n C n is normal for all symbols C n is diagonalized by the Discrete Fourier Transform matrix σ(c n) = {a(z) : z T n}, where T n := {e kπi/n, k =,, n 1}: ie, σ(c n) comprises the image of the nth roots of unity under the symbol Infinite dimensional generalization: Laurent operators C on l (Z) (doubly-infinite matrices) with spectrum σ(c ) = a(t)
47 Circulant Matrix and Laurent Operators σ(c ) = a(t) σ(c 1) σ(c ) σ(c )
48 Piecewise Continuous Symbols It is possible for the symbol a(z) = X k= a k z k to, eg, have a jump discontinuity This compromises the exponential growth of the norm of the resolvent [Böttcher, E, Trefethen, ] For example, take a(e i θ ) = θe i θ, a symbol studied by [Basor & Morrison, 199] a(t) σ(a ) σ ε(a ) σ ε(a 1)
49 Twisted Toeplitz Matrices A twisted Toeplitz matrix is a Toeplitz-like matrix with varying coefficients [Trefethen & Chapman, ] For example, for x j = πj/n, set A n = x 1 1 x1 x 1 xn 1 x n The symbol now depends on two variables: a(x, z) = x + 1 xz The pseudospectra resemble those of standard Toeplitz matrices, but the (pseudo)-eigenvectors have an entirely different character
50 Twisted Toeplitz Matrices For x j = πj/n, set A n = x 1 1 x1 x 1 xn 1 x n Pseudospectra of σ ε(a 1)
51 Eigenvectors of a Twisted Toeplitz Matrix Eigenvectors form wave packets for twisted Toeplitz matrices Pseudoeigenvectors for z σ ε(a) have a similar form Eigenvectors of A for a(x, z) = x + 1 xz
52 Looking Forward to Differential Operators Where do Toeplitz (and twisted Toeplitz) matrices come from? Numerous applications be we will highlight just one of them: discretization of differential operators Consider the steady state convection diffusion equation Lu = f, where Lu = u + c u posed over for x [, 1] with u() = u(1) = Fix a discretization parameter n Approximate the problem on a simple grid {x j } n+1 j= h = 1/(n + 1): x j = j h with spacing Replace the first and second derivatives with second-order accurate formulas on the grid: u (x j ) = u(x j+1) u(x j+1 ) h + O(h ) u (x j ) = u(x j 1) u(x j ) + u(x j+1 ) h + O(h )
53 Finite Difference Discretizations Now let u j u(x j ) with u = u n+1 =, so that Approximate Lu = f, with u (x j ) u j+1 u j+1 h u (x j ) u j 1 u j + u j+1 h Lu = u + c u, u() = u(1) =, via 1 h 1 + ch/ 1 ch/ 1 + ch/ 1 ch/ u 1 u u n = f (x 1) f (x ) f (x n) Label these components: A nu = f A n is Toeplitz for this constant-coefficient differential operator
54 Finite Difference Discretizations A n := 1 h 1 + ch/ 1 ch/ 1 + ch/ 1 ch/ Notice that the symbol of A n depends on n (recall h = 1/(n + 1)): 1 a n(z) = h c ««z h h h + c «z h
55 Finite Difference Discretizations A n := 1 h 1 + ch/ 1 ch/ 1 + ch/ 1 ch/ Notice that the symbol of A n depends on n (recall h = 1/(n + 1)): 1 a n(z) = h c ««z h h h + c «z h For all n, the symbol curve a n(t) is an ellipse
56 Finite Difference Discretizations A n := 1 h 1 + ch/ 1 ch/ 1 + ch/ 1 ch/ Notice that the symbol of A n depends on n (recall h = 1/(n + 1)): 1 a n(z) = h c ««z h h h + c «z h For all n, the symbol curve a n(t) is an ellipse If c = (no convection), then a n(t) is a real line segment: A n and L are both self-adjoint, hence normal If c, the eigenvalues of A n will be real, but the pseudospectra can be far from these eigenvalues
57 Finite Difference Discretizations A n := 1 h 1 + ch/ 1 ch/ 1 + ch/ 1 ch/ Notice that the symbol of A n depends on n (recall h = 1/(n + 1)): 1 a n(z) = h c ««z h h h + c «z h For all n, the symbol curve a n(t) is an ellipse If c = (no convection), then a n(t) is a real line segment: A n and L are both self-adjoint, hence normal If c, the eigenvalues of A n will be real, but the pseudospectra can be far from these eigenvalues Rightmost part of σ ε(a n) approximates the corresponding part of σ ε(l)
58 Finite Differences: Symbol Curves Symbol curves a n(t) for n = 1,,
59 Finite Differences: Symbol Curves Symbol curves a n(t) for n = 1,,
60 Finite Differences: Approximate Pseudospectra Rightmost part of σ ε(a n) for n =
61 Finite Differences: Approximate Pseudospectra Rightmost part of σ ε(a n) for n = 18
62 Finite Differences: Approximate Pseudospectra Rightmost part of σ ε(a n) for n =
63 Finite Differences: Approximate Pseudospectra Rightmost part of σ ε(a n) for n = 1
64 (b) Balanced Truncation Model Reduction
65 Balanced Truncation Model Reduction Consider the single-input, single-output (SISO) linear dynamical system: ẋ(t) = Ax(t) + bu(t) y(t) = cx(t), A C n n, b, c T C n We assume that A is stable: α(a) < We wish to reduce the dimension of the dynamical system by projecting onto well-chosen subspaces Balanced truncation: Change basis to match states that are easy to reach and easy to observe, then project onto that prominent subspace
66 Controllability and Observability Gramians To guage the observability of an initial state x = bx, measure the energy in its output (when there is no input, u = ): Then Z y(t) dt = y(t) = ce ta bx Z = bx Z bx e ta c ce ta bx dt e ta c ce ta dt bx =: bx Qbx
67 Controllability and Observability Gramians To guage the observability of an initial state x = bx, measure the energy in its output (when there is no input, u = ): Then Z y(t) dt = y(t) = ce ta bx Z = bx Z bx e ta c ce ta bx dt e ta c ce ta dt bx =: bx Qbx Similarly, we measure the controllability by the total input energy required to steer x = to a target state bx as t The special form of u that drives the system to bx with minimal energy satisfies Z u(t) dt = bx Z 1bx e ta bb e dt ta =: bx P 1 bx
68 Controllability and Observability Gramians To guage the observability of an initial state x = bx, measure the energy in its output (when there is no input, u = ): Then Z y(t) dt = y(t) = ce ta bx Z = bx Z bx e ta c ce ta bx dt e ta c ce ta dt bx =: bx Qbx Similarly, we measure the controllability by the total input energy required to steer x = to a target state bx as t The special form of u that drives the system to bx with minimal energy satisfies Z u(t) dt = bx Z 1bx e ta bb e dt ta =: bx P 1 bx Thus we have the infinite controllability and observability gramians P and Q: P := Z See, eg, [Antoulas, ] e ta bb e ta dt, Q := Z e ta c ce ta dt
69 Balanced Truncation Model Reduction The gramians P := Z e ta bb e ta dt, Q := Z e ta c ce ta dt (Hermitian positive definite, for a controllable and observable stable system) can be determined by solving the Lyapunov equations see the next lecture If x =, the minimum energy of u required to drive x to state bx is bx P 1 bx Starting from x = bx with u(t), the energy of output y is bx Qbx bx P 1 bx: bx is hard to reach if it is rich in the lowest modes of P bx Q bx: bx is hard to observe if it is rich in the lowest modes of Q Balanced truncation transforms the state space coordinate system to make these two gramians the same, then it truncates the lowest modes
70 Balanced Truncation Model Reduction Consider a generic coordinate transformation, for S invertible: (Sx) (t) = (SAS 1 )(Sx(t)) + (Sb)u(t) y(t) = (cs 1 )(Sx(t)) + du(t), (Sx)() = Sx With this transformation, the controllability and observability gramians are bp = SPS, b Q = S QS 1 For balancing, we seek S so that b P = b Q are diagonal
71 Balanced Truncation Model Reduction Consider a generic coordinate transformation, for S invertible: (Sx) (t) = (SAS 1 )(Sx(t)) + (Sb)u(t) y(t) = (cs 1 )(Sx(t)) + du(t), (Sx)() = Sx With this transformation, the controllability and observability gramians are bp = SPS, b Q = S QS 1 For balancing, we seek S so that b P = b Q are diagonal Observation (How does nonnormality affect balancing?) σ ε/κ(s) (SAS 1 ) σ ε(a) σ εκ(s) (SAS 1 ) The choice of internal coordinates will affect P, Q,
72 Balanced Truncation Model Reduction Consider a generic coordinate transformation, for S invertible: (Sx) (t) = (SAS 1 )(Sx(t)) + (Sb)u(t) y(t) = (cs 1 )(Sx(t)) + du(t), (Sx)() = Sx With this transformation, the controllability and observability gramians are bp = SPS, b Q = S QS 1 For balancing, we seek S so that b P = b Q are diagonal Observation (How does nonnormality affect balancing?) σ ε/κ(s) (SAS 1 ) σ ε(a) σ εκ(s) (SAS 1 ) The choice of internal coordinates will affect P, Q, but not the Hankel singular values: b P b Q = SPQS 1, and not the transfer function: (cs 1 )(z SAS 1 ) 1 (Sb) = d + c(z A) 1 b, and not the system moments: (cs 1 )(Sb) = cb, (cs 1 )(SAS 1 )(Sb) = cab,
MATH 5524 MATRIX THEORY Problem Set 5
MATH 554 MATRIX THEORY Problem Set 5 Posted Tuesday April 07. Due Tuesday 8 April 07. [Late work is due on Wednesday 9 April 07.] Complete any four problems, 5 points each. Recall the definitions of the
More informationFEL3210 Multivariable Feedback Control
FEL3210 Multivariable Feedback Control Lecture 8: Youla parametrization, LMIs, Model Reduction and Summary [Ch. 11-12] Elling W. Jacobsen, Automatic Control Lab, KTH Lecture 8: Youla, LMIs, Model Reduction
More informationGrammians. Matthew M. Peet. Lecture 20: Grammians. Illinois Institute of Technology
Grammians Matthew M. Peet Illinois Institute of Technology Lecture 2: Grammians Lyapunov Equations Proposition 1. Suppose A is Hurwitz and Q is a square matrix. Then X = e AT s Qe As ds is the unique solution
More informationNon-Self-Adjoint Operators and Pseudospectra
Non-Self-Adjoint Operators and Pseudospectra E.B. Davies King s College London 1 The Relevance of Pseudospectra The theory of pseudospectra is a way of imposing some order on the chaos of non-self-adjoint
More informationHOMEWORK PROBLEMS FROM STRANG S LINEAR ALGEBRA AND ITS APPLICATIONS (4TH EDITION)
HOMEWORK PROBLEMS FROM STRANG S LINEAR ALGEBRA AND ITS APPLICATIONS (4TH EDITION) PROFESSOR STEVEN MILLER: BROWN UNIVERSITY: SPRING 2007 1. CHAPTER 1: MATRICES AND GAUSSIAN ELIMINATION Page 9, # 3: Describe
More informationLecture Note 1: Background
ECE5463: Introduction to Robotics Lecture Note 1: Background Prof. Wei Zhang Department of Electrical and Computer Engineering Ohio State University Columbus, Ohio, USA Spring 2018 Lecture 1 (ECE5463 Sp18)
More informationLecture 1 INF-MAT3350/ : Some Tridiagonal Matrix Problems
Lecture 1 INF-MAT3350/4350 2007: Some Tridiagonal Matrix Problems Tom Lyche University of Oslo Norway Lecture 1 INF-MAT3350/4350 2007: Some Tridiagonal Matrix Problems p.1/33 Plan for the day 1. Notation
More informationChapter 7. Canonical Forms. 7.1 Eigenvalues and Eigenvectors
Chapter 7 Canonical Forms 7.1 Eigenvalues and Eigenvectors Definition 7.1.1. Let V be a vector space over the field F and let T be a linear operator on V. An eigenvalue of T is a scalar λ F such that there
More information1 Math 241A-B Homework Problem List for F2015 and W2016
1 Math 241A-B Homework Problem List for F2015 W2016 1.1 Homework 1. Due Wednesday, October 7, 2015 Notation 1.1 Let U be any set, g be a positive function on U, Y be a normed space. For any f : U Y let
More informationProblem Set (T) If A is an m n matrix, B is an n p matrix and D is a p s matrix, then show
MTH 0: Linear Algebra Department of Mathematics and Statistics Indian Institute of Technology - Kanpur Problem Set Problems marked (T) are for discussions in Tutorial sessions (T) If A is an m n matrix,
More informationNonnormality in Lyapunov Equations
RICE UNIVERSITY Nonnormality in Lyapunov Equations by Jonathan Baker A Thesis Submitted in Partial Fulfillment of the Requirements for the Degree Master of Arts Approved, Thesis Committee: Danny Sorensen,
More informationMath 108b: Notes on the Spectral Theorem
Math 108b: Notes on the Spectral Theorem From section 6.3, we know that every linear operator T on a finite dimensional inner product space V has an adjoint. (T is defined as the unique linear operator
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 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 informationIntroduction to Modern Control MT 2016
CDT Autonomous and Intelligent Machines & Systems Introduction to Modern Control MT 2016 Alessandro Abate Lecture 2 First-order ordinary differential equations (ODE) Solution of a linear ODE Hints to nonlinear
More informationforms Christopher Engström November 14, 2014 MAA704: Matrix factorization and canonical forms Matrix properties Matrix factorization Canonical forms
Christopher Engström November 14, 2014 Hermitian LU QR echelon Contents of todays lecture Some interesting / useful / important of matrices Hermitian LU QR echelon Rewriting a as a product of several matrices.
More informationThe goal of this chapter is to study linear systems of ordinary differential equations: dt,..., dx ) T
1 1 Linear Systems The goal of this chapter is to study linear systems of ordinary differential equations: ẋ = Ax, x(0) = x 0, (1) where x R n, A is an n n matrix and ẋ = dx ( dt = dx1 dt,..., dx ) T n.
More informationME 234, Lyapunov and Riccati Problems. 1. This problem is to recall some facts and formulae you already know. e Aτ BB e A τ dτ
ME 234, Lyapunov and Riccati Problems. This problem is to recall some facts and formulae you already know. (a) Let A and B be matrices of appropriate dimension. Show that (A, B) is controllable if and
More informationECE504: Lecture 8. D. Richard Brown III. Worcester Polytechnic Institute. 28-Oct-2008
ECE504: Lecture 8 D. Richard Brown III Worcester Polytechnic Institute 28-Oct-2008 Worcester Polytechnic Institute D. Richard Brown III 28-Oct-2008 1 / 30 Lecture 8 Major Topics ECE504: Lecture 8 We are
More informationECEN 605 LINEAR SYSTEMS. Lecture 7 Solution of State Equations 1/77
1/77 ECEN 605 LINEAR SYSTEMS Lecture 7 Solution of State Equations Solution of State Space Equations Recall from the previous Lecture note, for a system: ẋ(t) = A x(t) + B u(t) y(t) = C x(t) + D u(t),
More informationFinal A. Problem Points Score Total 100. Math115A Nadja Hempel 03/23/2017
Final A Math115A Nadja Hempel 03/23/2017 nadja@math.ucla.edu Name: UID: Problem Points Score 1 10 2 20 3 5 4 5 5 9 6 5 7 7 8 13 9 16 10 10 Total 100 1 2 Exercise 1. (10pt) Let T : V V be a linear transformation.
More informationNumerical Methods I Eigenvalue Problems
Numerical Methods I Eigenvalue Problems Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 October 2nd, 2014 A. Donev (Courant Institute) Lecture
More information8. Diagonalization.
8. Diagonalization 8.1. Matrix Representations of Linear Transformations Matrix of A Linear Operator with Respect to A Basis We know that every linear transformation T: R n R m has an associated standard
More informationAnalysis Preliminary Exam Workshop: Hilbert Spaces
Analysis Preliminary Exam Workshop: Hilbert Spaces 1. Hilbert spaces A Hilbert space H is a complete real or complex inner product space. Consider complex Hilbert spaces for definiteness. If (, ) : H H
More informationExamples include: (a) the Lorenz system for climate and weather modeling (b) the Hodgkin-Huxley system for neuron modeling
1 Introduction Many natural processes can be viewed as dynamical systems, where the system is represented by a set of state variables and its evolution governed by a set of differential equations. Examples
More informationEXERCISES ON DETERMINANTS, EIGENVALUES AND EIGENVECTORS. 1. Determinants
EXERCISES ON DETERMINANTS, EIGENVALUES AND EIGENVECTORS. Determinants Ex... Let A = 0 4 4 2 0 and B = 0 3 0. (a) Compute 0 0 0 0 A. (b) Compute det(2a 2 B), det(4a + B), det(2(a 3 B 2 )). 0 t Ex..2. For
More informationPseudospectra and Nonnormal Dynamical Systems
Pseudospectra and Nonnormal Dynamical Systems Mark Embree and Russell Carden Computational and Applied Mathematics Rice University Houston, Texas ELGERSBURG MARCH 2012 Overview of the Course These lectures
More informationMath 18, Linear Algebra, Lecture C00, Spring 2017 Review and Practice Problems for Final Exam
Math 8, Linear Algebra, Lecture C, Spring 7 Review and Practice Problems for Final Exam. The augmentedmatrix of a linear system has been transformed by row operations into 5 4 8. Determine if the system
More informationModule 09 From s-domain to time-domain From ODEs, TFs to State-Space Modern Control
Module 09 From s-domain to time-domain From ODEs, TFs to State-Space Modern Control Ahmad F. Taha EE 3413: Analysis and Desgin of Control Systems Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/
More informationELEMENTARY LINEAR ALGEBRA
ELEMENTARY LINEAR ALGEBRA K. R. MATTHEWS DEPARTMENT OF MATHEMATICS UNIVERSITY OF QUEENSLAND Second Online Version, December 1998 Comments to the author at krm@maths.uq.edu.au Contents 1 LINEAR EQUATIONS
More informationThis is a closed book exam. No notes or calculators are permitted. We will drop your lowest scoring question for you.
Math 54 Fall 2017 Practice Final Exam Exam date: 12/14/17 Time Limit: 170 Minutes Name: Student ID: GSI or Section: This exam contains 9 pages (including this cover page) and 10 problems. Problems are
More informationBalanced Truncation 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 2004: MODEL REDUCTION Balanced Truncation This lecture introduces balanced truncation for LTI
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 informationUniform individual asymptotics for the eigenvalues and eigenvectors of large Toeplitz matrices
Uniform individual asymptotics for the eigenvalues and eigenvectors of large Toeplitz matrices Sergei Grudsky CINVESTAV, Mexico City, Mexico The International Workshop WIENER-HOPF METHOD, TOEPLITZ OPERATORS,
More informationDiagonalization. P. Danziger. u B = A 1. B u S.
7., 8., 8.2 Diagonalization P. Danziger Change of Basis Given a basis of R n, B {v,..., v n }, we have seen that the matrix whose columns consist of these vectors can be thought of as a change of basis
More informationAdvanced Control Theory
State Space Solution and Realization chibum@seoultech.ac.kr Outline State space solution 2 Solution of state-space equations x t = Ax t + Bu t First, recall results for scalar equation: x t = a x t + b
More informationNovember 18, 2013 ANALYTIC FUNCTIONAL CALCULUS
November 8, 203 ANALYTIC FUNCTIONAL CALCULUS RODICA D. COSTIN Contents. The spectral projection theorem. Functional calculus 2.. The spectral projection theorem for self-adjoint matrices 2.2. The spectral
More informationToeplitz matrices. Niranjan U N. May 12, NITK, Surathkal. Definition Toeplitz theory Computational aspects References
Toeplitz matrices Niranjan U N NITK, Surathkal May 12, 2010 Niranjan U N (NITK, Surathkal) Linear Algebra May 12, 2010 1 / 15 1 Definition Toeplitz matrix Circulant matrix 2 Toeplitz theory Boundedness
More informationControl Systems. Dynamic response in the time domain. L. Lanari
Control Systems Dynamic response in the time domain L. Lanari outline A diagonalizable - real eigenvalues (aperiodic natural modes) - complex conjugate eigenvalues (pseudoperiodic natural modes) - phase
More informationDiagonalizing Matrices
Diagonalizing Matrices Massoud Malek A A Let A = A k be an n n non-singular matrix and let B = A = [B, B,, B k,, B n ] Then A n A B = A A 0 0 A k [B, B,, B k,, B n ] = 0 0 = I n 0 A n Notice that A i B
More informationTOEPLITZ OPERATORS. Toeplitz studied infinite matrices with NW-SE diagonals constant. f e C :
TOEPLITZ OPERATORS EFTON PARK 1. Introduction to Toeplitz Operators Otto Toeplitz lived from 1881-1940 in Goettingen, and it was pretty rough there, so he eventually went to Palestine and eventually contracted
More informationSPRING 2006 PRELIMINARY EXAMINATION SOLUTIONS
SPRING 006 PRELIMINARY EXAMINATION SOLUTIONS 1A. Let G be the subgroup of the free abelian group Z 4 consisting of all integer vectors (x, y, z, w) such that x + 3y + 5z + 7w = 0. (a) Determine a linearly
More informationControl Systems Design
ELEC4410 Control Systems Design Lecture 14: Controllability Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science Lecture 14: Controllability p.1/23 Outline
More informationModel reduction of large-scale dynamical systems
Model reduction of large-scale dynamical systems Lecture III: Krylov approximation and rational interpolation Thanos Antoulas Rice University and Jacobs University email: aca@rice.edu URL: www.ece.rice.edu/
More informationNATIONAL BOARD FOR HIGHER MATHEMATICS. Research Scholarships Screening Test. Saturday, January 20, Time Allowed: 150 Minutes Maximum Marks: 40
NATIONAL BOARD FOR HIGHER MATHEMATICS Research Scholarships Screening Test Saturday, January 2, 218 Time Allowed: 15 Minutes Maximum Marks: 4 Please read, carefully, the instructions that follow. INSTRUCTIONS
More informationAn Arnoldi Gram-Schmidt process and Hessenberg matrices for Orthonormal Polynomials
[ 1 ] University of Cyprus An Arnoldi Gram-Schmidt process and Hessenberg matrices for Orthonormal Polynomials Nikos Stylianopoulos, University of Cyprus New Perspectives in Univariate and Multivariate
More information33AH, WINTER 2018: STUDY GUIDE FOR FINAL EXAM
33AH, WINTER 2018: STUDY GUIDE FOR FINAL EXAM (UPDATED MARCH 17, 2018) The final exam will be cumulative, with a bit more weight on more recent material. This outline covers the what we ve done since the
More informationAnalysis II: Basic knowledge of real analysis: Part V, Power Series, Differentiation, and Taylor Series
.... Analysis II: Basic knowledge of real analysis: Part V, Power Series, Differentiation, and Taylor Series Kenichi Maruno Department of Mathematics, The University of Texas - Pan American March 4, 20
More informationThe Eigenvalue Problem: Perturbation Theory
Jim Lambers MAT 610 Summer Session 2009-10 Lecture 13 Notes These notes correspond to Sections 7.2 and 8.1 in the text. The Eigenvalue Problem: Perturbation Theory The Unsymmetric Eigenvalue Problem Just
More informationMATH 1120 (LINEAR ALGEBRA 1), FINAL EXAM FALL 2011 SOLUTIONS TO PRACTICE VERSION
MATH (LINEAR ALGEBRA ) FINAL EXAM FALL SOLUTIONS TO PRACTICE VERSION Problem (a) For each matrix below (i) find a basis for its column space (ii) find a basis for its row space (iii) determine whether
More information(a) If A is a 3 by 4 matrix, what does this tell us about its nullspace? Solution: dim N(A) 1, since rank(a) 3. Ax =
. (5 points) (a) If A is a 3 by 4 matrix, what does this tell us about its nullspace? dim N(A), since rank(a) 3. (b) If we also know that Ax = has no solution, what do we know about the rank of A? C(A)
More informationModule 03 Linear Systems Theory: Necessary Background
Module 03 Linear Systems Theory: Necessary Background Ahmad F. Taha EE 5243: Introduction to Cyber-Physical Systems Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ taha/index.html September
More informationOn the transient behaviour of stable linear systems
On the transient behaviour of stable linear systems D. Hinrichsen Institut für Dynamische Systeme Universität Bremen D-28334 Bremen Germany dh@math.uni-bremen.de A. J. Pritchard Mathematics Institute University
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 informationLINEAR ALGEBRA MICHAEL PENKAVA
LINEAR ALGEBRA MICHAEL PENKAVA 1. Linear Maps Definition 1.1. If V and W are vector spaces over the same field K, then a map λ : V W is called a linear map if it satisfies the two conditions below: (1)
More informationEvaluation of Chebyshev pseudospectral methods for third order differential equations
Numerical Algorithms 16 (1997) 255 281 255 Evaluation of Chebyshev pseudospectral methods for third order differential equations Rosemary Renaut a and Yi Su b a Department of Mathematics, Arizona State
More informationSPECTRAL THEOREM FOR COMPACT SELF-ADJOINT OPERATORS
SPECTRAL THEOREM FOR COMPACT SELF-ADJOINT OPERATORS G. RAMESH Contents Introduction 1 1. Bounded Operators 1 1.3. Examples 3 2. Compact Operators 5 2.1. Properties 6 3. The Spectral Theorem 9 3.3. Self-adjoint
More informationSolving large scale eigenvalue problems
arge scale eigenvalue problems, Lecture 5, March 23, 2016 1/30 Lecture 5, March 23, 2016: The QR algorithm II http://people.inf.ethz.ch/arbenz/ewp/ Peter Arbenz Computer Science Department, ETH Zürich
More informationEigenvalues, Eigenvectors and the Jordan Form
EE/ME 701: Advanced Linear Systems Eigenvalues, Eigenvectors and the Jordan Form Contents 1 Introduction 3 1.1 Review of basic facts about eigenvectors and eigenvalues..... 3 1.1.1 Looking at eigenvalues
More informationControllability, Observability, Full State Feedback, Observer Based Control
Multivariable Control Lecture 4 Controllability, Observability, Full State Feedback, Observer Based Control John T. Wen September 13, 24 Ref: 3.2-3.4 of Text Controllability ẋ = Ax + Bu; x() = x. At time
More informationIdentification Methods for Structural Systems
Prof. Dr. Eleni Chatzi System Stability Fundamentals Overview System Stability Assume given a dynamic system with input u(t) and output x(t). The stability property of a dynamic system can be defined from
More informationMATH 304 Linear Algebra Lecture 20: The Gram-Schmidt process (continued). Eigenvalues and eigenvectors.
MATH 304 Linear Algebra Lecture 20: The Gram-Schmidt process (continued). Eigenvalues and eigenvectors. Orthogonal sets Let V be a vector space with an inner product. Definition. Nonzero vectors v 1,v
More informationCayley-Hamilton Theorem
Cayley-Hamilton Theorem Massoud Malek 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 Let A be an n n matrix Although det (λ I n A
More informationThe following definition is fundamental.
1. Some Basics from Linear Algebra With these notes, I will try and clarify certain topics that I only quickly mention in class. First and foremost, I will assume that you are familiar with many basic
More informationELEMENTARY LINEAR ALGEBRA
ELEMENTARY LINEAR ALGEBRA K R MATTHEWS DEPARTMENT OF MATHEMATICS UNIVERSITY OF QUEENSLAND Second Online Version, December 998 Comments to the author at krm@mathsuqeduau All contents copyright c 99 Keith
More informationLecture notes: Applied linear algebra Part 1. Version 2
Lecture notes: Applied linear algebra Part 1. Version 2 Michael Karow Berlin University of Technology karow@math.tu-berlin.de October 2, 2008 1 Notation, basic notions and facts 1.1 Subspaces, range and
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 informationName: Final Exam MATH 3320
Name: Final Exam MATH 3320 Directions: Make sure to show all necessary work to receive full credit. If you need extra space please use the back of the sheet with appropriate labeling. (1) State the following
More information2 Two-Point Boundary Value Problems
2 Two-Point Boundary Value Problems Another fundamental equation, in addition to the heat eq. and the wave eq., is Poisson s equation: n j=1 2 u x 2 j The unknown is the function u = u(x 1, x 2,..., x
More information7. Symmetric Matrices and Quadratic Forms
Linear Algebra 7. Symmetric Matrices and Quadratic Forms CSIE NCU 1 7. Symmetric Matrices and Quadratic Forms 7.1 Diagonalization of symmetric matrices 2 7.2 Quadratic forms.. 9 7.4 The singular value
More informationNumerical Methods in Matrix Computations
Ake Bjorck Numerical Methods in Matrix Computations Springer Contents 1 Direct Methods for Linear Systems 1 1.1 Elements of Matrix Theory 1 1.1.1 Matrix Algebra 2 1.1.2 Vector Spaces 6 1.1.3 Submatrices
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 informationx 3y 2z = 6 1.2) 2x 4y 3z = 8 3x + 6y + 8z = 5 x + 3y 2z + 5t = 4 1.5) 2x + 8y z + 9t = 9 3x + 5y 12z + 17t = 7
Linear Algebra and its Applications-Lab 1 1) Use Gaussian elimination to solve the following systems x 1 + x 2 2x 3 + 4x 4 = 5 1.1) 2x 1 + 2x 2 3x 3 + x 4 = 3 3x 1 + 3x 2 4x 3 2x 4 = 1 x + y + 2z = 4 1.4)
More informationA fast and well-conditioned spectral method: The US method
A fast and well-conditioned spectral method: The US method Alex Townsend University of Oxford Leslie Fox Prize, 24th of June 23 Partially based on: S. Olver & T., A fast and well-conditioned spectral method,
More informationMath 21b Final Exam Thursday, May 15, 2003 Solutions
Math 2b Final Exam Thursday, May 5, 2003 Solutions. (20 points) True or False. No justification is necessary, simply circle T or F for each statement. T F (a) If W is a subspace of R n and x is not in
More informationA Smooth Operator, Operated Correctly
Clark Department of Mathematics Bard College at Simon s Rock March 6, 2014 Abstract By looking at familiar smooth functions in new ways, we can make sense of matrix-valued infinite series and operator-valued
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 information6 Inner Product Spaces
Lectures 16,17,18 6 Inner Product Spaces 6.1 Basic Definition Parallelogram law, the ability to measure angle between two vectors and in particular, the concept of perpendicularity make the euclidean space
More informationNotes on Eigenvalues, Singular Values and QR
Notes on Eigenvalues, Singular Values and QR Michael Overton, Numerical Computing, Spring 2017 March 30, 2017 1 Eigenvalues Everyone who has studied linear algebra knows the definition: given a square
More informationMATH 5524 MATRIX THEORY Pledged Problem Set 2
MATH 5524 MATRIX THEORY Pledged Problem Set 2 Posted Tuesday 25 April 207 Due by 5pm on Wednesday 3 May 207 Complete any four problems, 25 points each You are welcome to complete more problems if you like,
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 informationLinear Algebra March 16, 2019
Linear Algebra March 16, 2019 2 Contents 0.1 Notation................................ 4 1 Systems of linear equations, and matrices 5 1.1 Systems of linear equations..................... 5 1.2 Augmented
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 informationControl Systems. Laplace domain analysis
Control Systems Laplace domain analysis L. Lanari outline introduce the Laplace unilateral transform define its properties show its advantages in turning ODEs to algebraic equations define an Input/Output
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 informationNATIONAL BOARD FOR HIGHER MATHEMATICS. M. A. and M.Sc. Scholarship Test. September 25, Time Allowed: 150 Minutes Maximum Marks: 30
NATIONAL BOARD FOR HIGHER MATHEMATICS M. A. and M.Sc. Scholarship Test September 25, 2010 Time Allowed: 150 Minutes Maximum Marks: 30 Please read, carefully, the instructions on the following page 1 INSTRUCTIONS
More informationCompound matrices and some classical inequalities
Compound matrices and some classical inequalities Tin-Yau Tam Mathematics & Statistics Auburn University Dec. 3, 04 We discuss some elegant proofs of several classical inequalities of matrices by using
More informationLINEAR ALGEBRA 1, 2012-I PARTIAL EXAM 3 SOLUTIONS TO PRACTICE PROBLEMS
LINEAR ALGEBRA, -I PARTIAL EXAM SOLUTIONS TO PRACTICE PROBLEMS Problem (a) For each of the two matrices below, (i) determine whether it is diagonalizable, (ii) determine whether it is orthogonally diagonalizable,
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 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 informationProblems Session. Nikos Stylianopoulos University of Cyprus
[ 1 ] University of Cyprus Problems Session Nikos Stylianopoulos University of Cyprus Hausdorff Geometry Of Polynomials And Polynomial Sequences Institut Mittag-Leffler Djursholm, Sweden May-June 2018
More informationNumerical Linear Algebra Homework Assignment - Week 2
Numerical Linear Algebra Homework Assignment - Week 2 Đoàn Trần Nguyên Tùng Student ID: 1411352 8th October 2016 Exercise 2.1: Show that if a matrix A is both triangular and unitary, then it is diagonal.
More informationA Brief Outline of Math 355
A Brief Outline of Math 355 Lecture 1 The geometry of linear equations; elimination with matrices A system of m linear equations with n unknowns can be thought of geometrically as m hyperplanes intersecting
More informationANSWERS. E k E 2 E 1 A = B
MATH 7- Final Exam Spring ANSWERS Essay Questions points Define an Elementary Matrix Display the fundamental matrix multiply equation which summarizes a sequence of swap, combination and multiply operations,
More informationBergman shift operators for Jordan domains
[ 1 ] University of Cyprus Bergman shift operators for Jordan domains Nikos Stylianopoulos University of Cyprus Num An 2012 Fifth Conference in Numerical Analysis University of Ioannina September 2012
More informationMATH 31 - ADDITIONAL PRACTICE PROBLEMS FOR FINAL
MATH 3 - ADDITIONAL PRACTICE PROBLEMS FOR FINAL MAIN TOPICS FOR THE FINAL EXAM:. Vectors. Dot product. Cross product. Geometric applications. 2. Row reduction. Null space, column space, row space, left
More informationFINITE-DIMENSIONAL LINEAR ALGEBRA
DISCRETE MATHEMATICS AND ITS APPLICATIONS Series Editor KENNETH H ROSEN FINITE-DIMENSIONAL LINEAR ALGEBRA Mark S Gockenbach Michigan Technological University Houghton, USA CRC Press Taylor & Francis Croup
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationAnalysis of Transient Behavior in Linear Differential Algebraic Equations
Analysis of Transient Behavior in Linear Differential Algebraic Equations Dr. Mark Embree and Blake Keeler May 6, 5 Abstract The linearized Navier-Stokes equations form a type of system called a differential
More informationNew York Journal of Mathematics. Spectral Approximation of Multiplication Operators
New York Journal of Mathematics New York J. Math. 1 (1995) 75 96. Spectral Approximation of Multiplication Operators Kent E. Morrison Abstract. A linear operator on a Hilbert space may be approximated
More information