computing the eigenvalues of a
play

Computing the eigenvalues of a Joint work with Raf Vandebril and - PowerPoint PPT Presentation

Mar 19, 2024 •469 likes •1.12k views

Computing the eigenvalues of a companion matrix Marc Van Barel Computing the eigenvalues of a Joint work with Raf Vandebril and Paul Van Dooren companion matrix Structured Linear Algebra Problems: Analysis, Algorithms, and Applications

  1. Computing the eigenvalues of a companion matrix Marc Van Barel Computing the eigenvalues of a Joint work with Raf Vandebril and Paul Van Dooren companion matrix Structured Linear Algebra Problems: Analysis, Algorithms, and Applications The problem Cortona, Italy Companion matrix Working with Givens September 15-19, 2008 transformations Representation Fusion and shift-through operations Unitary plus rank one matrices Structure under a QR -step A representation for the unitary matrix Representation of the unitary plus rank one matrix Implicit QR -algorithm with single shift Initialization Marc Van Barel The chasing Joint work with Raf Vandebril and Paul Van Dooren The last Givens transformation Dept. Computer Science, K.U.Leuven, Belgium Numerical Experiments Dept. of Math. Eng., Catholic University of Louvain, Belgium Scaling Comparison 1 / 64

  2. Contents Computing the eigenvalues of a companion matrix Marc Van Barel Joint work with Raf The problem Vandebril and Paul Van Dooren Companion matrix Working with Givens transformations Representation Fusion and shift-through operations The problem Companion matrix Unitary plus rank one matrices Working with Givens transformations Structure under a QR -step Representation A representation for the unitary matrix Fusion and shift-through operations Representation of the unitary plus rank one matrix Unitary plus rank one matrices Implicit QR -algorithm with single shift Structure under a QR -step A representation for the Initialization unitary matrix Representation of the The chasing unitary plus rank one matrix Implicit QR -algorithm The last Givens transformation with single shift Initialization Numerical Experiments The chasing The last Givens Scaling transformation Numerical Comparison Experiments Scaling Comparison 2 / 64

  3. Outline Computing the eigenvalues of a companion matrix Marc Van Barel Joint work with Raf The problem Vandebril and Paul Van Dooren Companion matrix Working with Givens transformations Representation Fusion and shift-through operations The problem Companion matrix Unitary plus rank one matrices Working with Givens transformations Structure under a QR -step Representation A representation for the unitary matrix Fusion and shift-through operations Representation of the unitary plus rank one matrix Unitary plus rank one matrices Implicit QR -algorithm with single shift Structure under a QR -step A representation for the Initialization unitary matrix Representation of the The chasing unitary plus rank one matrix Implicit QR -algorithm The last Givens transformation with single shift Initialization Numerical Experiments The chasing The last Givens Scaling transformation Numerical Comparison Experiments Scaling Comparison 3 / 64

  4. Outline Computing the eigenvalues of a companion matrix Marc Van Barel Joint work with Raf The problem Vandebril and Paul Van Dooren Companion matrix Working with Givens transformations Representation Fusion and shift-through operations The problem Companion matrix Unitary plus rank one matrices Working with Givens transformations Structure under a QR -step Representation A representation for the unitary matrix Fusion and shift-through operations Representation of the unitary plus rank one matrix Unitary plus rank one matrices Implicit QR -algorithm with single shift Structure under a QR -step A representation for the Initialization unitary matrix Representation of the The chasing unitary plus rank one matrix Implicit QR -algorithm The last Givens transformation with single shift Initialization Numerical Experiments The chasing The last Givens Scaling transformation Numerical Comparison Experiments Scaling Comparison 4 / 64

  5. Companion matrix Computing the eigenvalues of a companion matrix Marc Van Barel Definition Joint work with Raf Vandebril and Paul Van Given a monic polynomial Dooren p ( z ) = p 0 + p 1 z + p 2 z 2 + ... + p n − 1 z n − 1 + z n with the coefficients p i ∈ R or C The problem the associated companion matrix C p is defined as: Companion matrix Working with Givens transformations   0 − p 0 Representation Fusion and shift-through ... operations   1 − p 1   C p = . Unitary plus rank one  .  ... matrices .   . Structure under a QR -step   A representation for the 1 − p n − 1 unitary matrix Representation of the unitary plus rank one matrix Implicit QR -algorithm with single shift The eigenvalues of the companion matrix coincide with the Initialization The chasing zeros of the associated polynomial p ( z ) , because The last Givens transformation p ( z ) = det ( zI − C p ) . Numerical Experiments Scaling Comparison 5 / 64

  6. Unitary plus rank one matrix Computing the eigenvalues of a companion matrix Marc Van Barel Joint work with Raf Vandebril and Paul Van C p = H upper Hessenberg Dooren U + uv H = unitary plus rank one with   0 0 ··· 0 ± 1 The problem 1 0 ··· 0 0 Companion matrix   Working with Givens   0 1 ··· 0 0 U = transformations   . . . .  ...  Representation . . . .   . . . . Fusion and shift-through   operations 0 0 ··· 1 0 Unitary plus rank one matrices v T Structure under a QR -step = ( 0 , 0 ,..., 0 , 1 ) A representation for the unitary matrix u T = ( − p 0 ∓ 1 , − p 1 ,..., − p n − 1 ) Representation of the unitary plus rank one matrix Hence, computing the zeros of the polynomial p ( z ) is Implicit QR -algorithm with single shift equivalent to computing the eigenvalues of the upper Initialization The chasing Hessenberg, unitary plus rank one matrix The last Givens transformation C p = H = U + uv H . Numerical Experiments Scaling Comparison 6 / 64

  7. Outline Computing the eigenvalues of a companion matrix Marc Van Barel Joint work with Raf The problem Vandebril and Paul Van Dooren Companion matrix Working with Givens transformations Representation Fusion and shift-through operations The problem Companion matrix Unitary plus rank one matrices Working with Givens transformations Structure under a QR -step Representation A representation for the unitary matrix Fusion and shift-through operations Representation of the unitary plus rank one matrix Unitary plus rank one matrices Implicit QR -algorithm with single shift Structure under a QR -step A representation for the Initialization unitary matrix Representation of the The chasing unitary plus rank one matrix Implicit QR -algorithm The last Givens transformation with single shift Initialization Numerical Experiments The chasing The last Givens Scaling transformation Numerical Comparison Experiments Scaling Comparison 7 / 64

  8. Outline Computing the eigenvalues of a companion matrix Marc Van Barel Joint work with Raf The problem Vandebril and Paul Van Dooren Companion matrix Working with Givens transformations Representation Fusion and shift-through operations The problem Companion matrix Unitary plus rank one matrices Working with Givens transformations Structure under a QR -step Representation A representation for the unitary matrix Fusion and shift-through operations Representation of the unitary plus rank one matrix Unitary plus rank one matrices Implicit QR -algorithm with single shift Structure under a QR -step A representation for the Initialization unitary matrix Representation of the The chasing unitary plus rank one matrix Implicit QR -algorithm The last Givens transformation with single shift Initialization Numerical Experiments The chasing The last Givens Scaling transformation Numerical Comparison Experiments Scaling Comparison 8 / 64

  9. Representation Computing the eigenvalues of a companion matrix ◮ Givens transformations are a powerful tool for working with Marc Van Barel Joint work with Raf structured matrices. Vandebril and Paul Van Dooren ◮ representation of sequences of Givens transformations by a graphical scheme. ◮ example: the QR -factorization of a 6 × 6 Hessenberg matrix The problem ➊ × × × × × × Companion matrix � � ➋ × × × × × Working with Givens � � transformations ➌ × × × × � Representation � ➍ × × × � Fusion and shift-through � × × operations ➎ � � × ➏ Unitary plus rank one matrices 5 4 3 2 1 Structure under a QR -step A representation for the unitary matrix The figure corresponds to G 5 G 4 ... G 1 R . Representation of the Rewriting the formula we have G H 1 ... G H unitary plus rank one matrix 5 H = R . Implicit QR -algorithm Hence G H 5 annihilates the first subdiagonal element, G H with single shift 4 Initialization the second and so forth. The chasing The last Givens transformation Numerical Experiments Scaling Comparison 9 / 64

Recommend

Matrix Calculations: Eigenvalues and Eigenvectors H. Geuvers (and A. Kissinger) Institute for

Matrix Calculations: Eigenvalues and Eigenvectors H. Geuvers (and A. Kissinger) Institute for

Eigenvalues and Eigenvectors Radboud University Nijmegen Applications of Eigenvalues and Eigenvectors Matrix Calculations: Eigenvalues and Eigenvectors H. Geuvers (and A. Kissinger) Institute for Computing and Information Sciences Radboud

923 views • 35 slides
A refined Lanczos method for computing eigenvalues and eigenvectors of unsymmetric matrices Jean

A refined Lanczos method for computing eigenvalues and eigenvectors of unsymmetric matrices Jean

A refined Lanczos method for computing eigenvalues and eigenvectors of unsymmetric matrices Jean Christophe Tremblay and Tucker Carrington Chemistry Department Queens University 23 ao ut 2007 We want to solve the eigenvalue problem GX (

1.03k views • 34 slides
Chapter 5 Eigenvalues and Eigenvectors Section 5.1 Eigenvectors and Eigenvalues Motivation:

Chapter 5 Eigenvalues and Eigenvectors Section 5.1 Eigenvectors and Eigenvalues Motivation:

Chapter 5 Eigenvalues and Eigenvectors Section 5.1 Eigenvectors and Eigenvalues Motivation: Difference equations A Biology Question How to predict a population of rabbits with given dynamics : 1. half of the newborn rabbits survive their first

353 views • 22 slides
Eigenvalues, Eigenvectors, Matrix Factoring, and Principal Components The eigenvalues and

Eigenvalues, Eigenvectors, Matrix Factoring, and Principal Components The eigenvalues and

Eigenvalues, Eigenvectors, Matrix Factoring, and Principal Components The eigenvalues and eigenvectors of a square matrix play a key role in some important operations in statistics. In particular, they are intimately connected with the

376 views • 25 slides
Section 5.5 Complex Eigenvalues (Part II) Motivation: Complex Versus Two Real Eigenvalues

Section 5.5 Complex Eigenvalues (Part II) Motivation: Complex Versus Two Real Eigenvalues

Section 5.5 Complex Eigenvalues (Part II) Motivation: Complex Versus Two Real Eigenvalues Todays decomposition is very analogous to diagonalization. Theorem Let A be a 2 2 matrix with linearly independent eigenvectors v 1 , v 2 and

700 views • 16 slides
For matrix A ( p p ) with real eigenvalues, define F A , the empirical distribution function of

For matrix A ( p p ) with real eigenvalues, define F A , the empirical distribution function of

For matrix A ( p p ) with real eigenvalues, define F A , the empirical distribution function of the eigenvalues of A , to be F A ( x ) (1 /p ) (number of eigenvalues of A x ) . For and p.d.f. G the Stieltjes transform of G is defined

447 views • 19 slides
For matrix A ( p p ) with real eigenvalues, define F A , the empirical distribution function of

For matrix A ( p p ) with real eigenvalues, define F A , the empirical distribution function of

For matrix A ( p p ) with real eigenvalues, define F A , the empirical distribution function of the eigenvalues of A , to be F A ( x ) (1 /p ) (number of eigenvalues of A x ) . For and p.d.f. G the Stieltjes transform of G is defined

451 views • 24 slides
Eigenvalues, Eigenvectors, and Diagonalization Diagonalization Math 240 Calculus III Summer

Eigenvalues, Eigenvectors, and Diagonalization Diagonalization Math 240 Calculus III Summer

Eigenvalues, Eigenvectors, and Diagonal- ization Math 240 Eigenvalues and Eigenvectors Eigenvalues, Eigenvectors, and Diagonalization Diagonalization Math 240 Calculus III Summer 2015, Session II Thursday, July 16, 2015 Eigenvalues,

171 views • 15 slides
Homogeneous Linear Systems Three Cases: Distinct Real Eigenvalues Repeated Eigenvalues

Homogeneous Linear Systems Three Cases: Distinct Real Eigenvalues Repeated Eigenvalues

Homogeneous Linear Systems Three Cases: Distinct Real Eigenvalues Repeated Eigenvalues Complex Eigenvalues MSW MTH 291 Distinct Eigenvalues General Solution-Homogeneous Systems Let 1 , 2 , . . . , n be n distinct real

264 views • 4 slides
Incremental Classification with Generalized Eigenvalues Mario Rosario Guarracino September 17,

Incremental Classification with Generalized Eigenvalues Mario Rosario Guarracino September 17,

High Performance Computing and Networking Institute National Research Council, Italy The Data Reference Model: Incremental Classification with Generalized Eigenvalues Mario Rosario Guarracino September 17, 2007 11/30/2007 8:53 AM

1.26k views • 49 slides
( A I ) = 0 However, there is no exact formula for polynomials of degree five

( A I ) = 0 However, there is no exact formula for polynomials of degree five

Lesson 21 C OMPUTATION OF EIGENVALUES We consider the problem of computation of eigenvalues To make life (much!) simpler, let's assume we have a symmetric matrix A R n n , which has decomposition A = Q Q Unlike previous

503 views • 26 slides
Fredholm Determinants: 2 , 000 A Robust 0 Approach to 2 , 000 Computing Stokes

Fredholm Determinants: 2 , 000 A Robust 0 Approach to 2 , 000 Computing Stokes

Fredholm Determinants: 2 , 000 A Robust 0 Approach to 2 , 000 Computing Stokes Eigenvalues 0 2 4 6 8 10 k Travis Askham (New Jersey Institute of Technology) SIAM CSE 2019. Spokane, WA, USA. Joint work with Manas Rachh (Flatiron

1.39k views • 105 slides
Linear Algebra Chapter 5: Eigenvalues and Eigenvectors Section 5.1. Eigenvalues and

Linear Algebra Chapter 5: Eigenvalues and Eigenvectors Section 5.1. Eigenvalues and

Linear Algebra Chapter 5: Eigenvalues and Eigenvectors Section 5.1. Eigenvalues and EigenvectorsProofs of Theorems April 10, 2020 () Linear Algebra April 10, 2020 1 / 23 Table of contents Page 300 Number 8 1 Page 300 Number 14 2

738 views • 60 slides
Eigenvalues Definition 1 Given an n n matrix A , a scalar C , and a non zero vector v

Eigenvalues Definition 1 Given an n n matrix A , a scalar C , and a non zero vector v

7.1 Eigenvalues & Eigenvectors P. Danziger Eigenvalues Definition 1 Given an n n matrix A , a scalar C , and a non zero vector v R n we say that is an eigenvalue of A , with corresponding eigenvalue v if A v = v Notes

623 views • 13 slides
Overview The previous lecture introduced eigenvalues and eigenvectors. Well review these

Overview The previous lecture introduced eigenvalues and eigenvectors. Well review these

Overview The previous lecture introduced eigenvalues and eigenvectors. Well review these definitions before considering the following question: Question Given a square matrix A, how can you find the eigenvalues of A? Well discuss an

685 views • 42 slides
Eigenvalues and Eigenvectors Let A R n n be a matrix. If R and v R n , v = 0,

Eigenvalues and Eigenvectors Let A R n n be a matrix. If R and v R n , v = 0,

Eigenvalues and Eigenvectors Let A R n n be a matrix. If R and v R n , v = 0, with Av = v , then we call 1. an eigenvalue of A , 2. v an eigenvector of A , 3. and ( , v ) an eigenpair of A Eigen, adjective: own,

730 views • 51 slides
Eigenvalues and Eigenvectors Artem Los (arteml@kth.se) February 9th, 2017 Artem Los

Eigenvalues and Eigenvectors Artem Los (arteml@kth.se) February 9th, 2017 Artem Los

Eigenvalues and Eigenvectors Artem Los (arteml@kth.se) February 9th, 2017 Artem Los (arteml@kth.se) Eigenvalues and Eigenvectors February 9th, 2017 1 / 16 Overview What is an Eigenvalue? 1 Finding eigenvalues and eigenvectors 2 Exam

780 views • 45 slides
Transmission Eigenvalues in Inverse Scattering Theory David Colton Department of Mathematical

Transmission Eigenvalues in Inverse Scattering Theory David Colton Department of Mathematical

TE and Scattering Theory Spherically Stratified Media Transmission Eigenvalues Open Problem Transmission Eigenvalues in Inverse Scattering Theory David Colton Department of Mathematical Sciences University of Delaware Newark, DE 19716, USA

495 views • 28 slides
Linear Algebra Chapter 5: Eigenvalues and Eigenvectors Section 5.2. DiagonalizationProofs of

Linear Algebra Chapter 5: Eigenvalues and Eigenvectors Section 5.2. DiagonalizationProofs of

Linear Algebra Chapter 5: Eigenvalues and Eigenvectors Section 5.2. DiagonalizationProofs of Theorems April 13, 2020 () Linear Algebra April 13, 2020 1 / 34 Table of contents Theorem 5.2. Matrix Summary of Eigenvalues of A 1 Corollary

1.11k views • 79 slides
Overview Before the break, we began to study eigenvectors and eigenvalues, introducing the

Overview Before the break, we began to study eigenvectors and eigenvalues, introducing the

Overview Before the break, we began to study eigenvectors and eigenvalues, introducing the characteristic equation as a tool for finding the eigenvalues of a matrix: det( A I ) = 0 . The roots of the characteristic equation are the

423 views • 15 slides
Eigenvalues and Eigenvectors Raibatak Sen Gupta 2019 Eigenvalues Characteristic Equation and

Eigenvalues and Eigenvectors Raibatak Sen Gupta 2019 Eigenvalues Characteristic Equation and

Eigenvalues and Eigenvectors Raibatak Sen Gupta Eigenvalues and Eigenvectors Raibatak Sen Gupta 2019 Eigenvalues Characteristic Equation and Eigenvectors Raibatak Sen Gupta Let A be an n n matrix. Then det ( A xI n ) gives a

355 views • 24 slides
Eigenvalues and Eigenvectors Eigenvalue problem Let ! be an "" matrix: $ & is

Eigenvalues and Eigenvectors Eigenvalue problem Let ! be an "" matrix: $ & is

Eigenvalues and Eigenvectors Eigenvalue problem Let ! be an "" matrix: $ & is an eigenvector of ! if there exists a scalar ' such that ! $ = ' $ I l eigenvectors where ' is called an eigenvalue . eigenvalues eigeipairs If $ is

723 views • 13 slides
Fast Stable Computation of the Eigenvalues of Companion Matrices David S. Watkins Department of

Fast Stable Computation of the Eigenvalues of Companion Matrices David S. Watkins Department of

Fast Stable Computation of the Eigenvalues of Companion Matrices David S. Watkins Department of Mathematics Washington State University Spa, Belgium, 2014 David S. Watkins Eigenvalues of Companion Matrices Collaborators This is joint work

1.57k views • 111 slides
Matrices with Integer Eigenvalues 4 1 1 0 1 1 4 0 1 1

Matrices with Integer Eigenvalues 4 1 1 0 1 1 4 0 1 1

Matrices with Integer Eigenvalues 4 1 1 0 1 1 4 0 1 1 1 0 3 0 0 0 1 0 3 0 1 1 0 0 3 Ron Adin Bar-Ilan University

758 views • 27 slides

More recommend