Ma2213 lecture 8
1 / 35

MA2213 Lecture 8 - PowerPoint PPT Presentation

  • Uploaded on

MA2213 Lecture 8. Eigenvectors. Application of Eigenvectors. Vufoil 18, lecture 7 : The Fibonacci sequence satisfies. Fibonacci Ratio Sequence. Fibonacci Ratio Sequence. Another Biomathematics Application.

I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.
Download Presentation

PowerPoint Slideshow about ' MA2213 Lecture 8' - hewitt

An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript
Ma2213 lecture 8

MA2213 Lecture 8


Application of eigenvectors
Application of Eigenvectors

Vufoil 18, lecture 7 : The Fibonacci sequence satisfies

Another biomathematics application
Another Biomathematics Application

Leonardo da Pisa, better known as Fibonacci, invented his famous sequence to compute the reproductive success of rabbits* Similar sequences describe frequencies in males, females of a sex-linked gene.

For genes (2 alleles) carried in the X chromosome**

The solution has the form


*page i, ** pages 10-12 in The Theory of Evolution and Dynamical Systems ,J. Hofbauer and K. Sigmund, 1984.

Eigenvector problem pages 333 351
Eigenvector Problem (pages 333-351)

is a square matrix then a nonzero

Recall that if


is an eigenvector corresponding to the



Eigenvectors and eigenvalues arise in biomathematics where they describe growth and population genetics

They arise in numerical solution of linear equations because they determine convergence properties

They arise in physical problems, especially those that involve vibrations in which eigenvalues are related to vibration frequencies

Example 7 2 1 pages 333 334
Example 7.2.1 pages 333-334


the eigenvalue-eigenvector pairs are


We observe that every (column) vector


Example 7 2 1 pages 333 3341
Example 7.2.1 pages 333-334

Therefore, since x  Ax is a linear transformation

and since

are eigenvectors

We can repeat this process to obtain


Question What happens as

Example 7 2 1 pages 333 3342
Example 7.2.1 pages 333-334

General Principle : If a vector v can be expressed as a linear combination of eigenvectors of a matrix A, then it is very easy to compute Av

It is possible to express every vector as a linear combination of eigenvectors of an n by n matrix A iff either of the following equivalent conditions is satisfied :

(i) there exists a basis consisting of eigenvectors of A

(ii) the sum of dimensions of eigenspaces of A = n

Question Does this condition hold for


Question What special form does this matrix have ?

Example 7 2 1 pages 333 3343
Example 7.2.1 pages 333-334

The characteristic polynomial of


2 is the (only) eigenvalue, it has algebraicmultiplicity 2

so the eigenspace

for eigenvalue 5

has dimension 1

the eigenvalue 5 is said to have geometricmultiplicity 1

Question What are alg.&geom. mult. in Example 7.2.7 ?

Characteristic polynomials pp 335 337
Characteristic Polynomials pp. 335-337

Example 7.22 (p. 335) The eigenvalue-eigenvector pairs

in Example 7.2.1 are

of the matrix

corresponding eigenvectors

Question What is the equation for


Eigenvalues of symmetric matrices
Eigenvalues of Symmetric Matrices

The following real symmetric matrices that we studied

have real eigenvalues and eigenvectors correspondingto distinct eigenvectors are orthogonal.

Question What are the eigenvalues of these matrices ?

Question What are the corresponding eigenvectors ?

Question Compute their scalar products

Eigenvalues of symmetric matrices1
Eigenvalues of Symmetric Matrices

Theorem 1. All eigenvalues of real symmetric matrices

are real valued.

Proof For a matrix

with complex (or real) entries


denote the matrix whose entries are the

complex conjugates of the entries of

Question Prove

is real (all entries are real) iff

Question Prove

Assume that

and observe that



Eigenvalues of symmetric matrices2
Eigenvalues of Symmetric Matrices

Theorem 2. Eigenvectors of a real symmetric matrix that

correspond to distinct eigenvalues are orthogonal.

Proof Assume that

Then compute

and observe that

Orthogonal matrices
Orthogonal Matrices

Definition A matrix

is orthogonal if


is orthogonal then

therefore either



is nonsingular and has an inverse




Permutation matrices
Permutation Matrices

Definition A matrix

is called a permutation

matrix if there exists a function (called a permutation)

that is 1-to-1 (and therefore onto) such that


Question Why is every permutation matrix orthogonal ?

Eigenvalues of symmetric matrices3
Eigenvalues of Symmetric Matrices

Theorem 7.2.4 pages 337-338 If

is symmetric


then there exists a set

eigenvalue-eigenvector pairs

Proof Uses Theorems 1 and 2 and a little linear algebra.

Choose eigenvectors so that

construct matrices

and observe that

Matlab eig command

>> help eig

EIG Eigenvalues and eigenvectors.

E = EIG(X) is a vector containing the eigenvalues of a square

matrix X.

[V,D] = EIG(X) produces a diagonal matrix D of eigenvalues and a

full matrix V whose columns are the corresponding eigenvectors so

that X*V = V*D.

[V,D] = EIG(X,'nobalance') performs the computation with balancing

disabled, which sometimes gives more accurate results for certain

problems with unusual scaling. If X is symmetric, EIG(X,'nobalance')

is ignored since X is already balanced.

E = EIG(A,B) is a vector containing the generalized eigenvalues

of square matrices A and B.

[V,D] = EIG(A,B) produces a diagonal matrix D of generalized

eigenvalues and a full matrix V whose columns are the

corresponding eigenvectors so that A*V = B*V*D.

EIG(A,B,'chol') is the same as EIG(A,B) for symmetric A and symmetric

positive definite B. It computes the generalized eigenvalues of A and B

using the Cholesky factorization of B.

EIG(A,B,'qz') ignores the symmetry of A and B and uses the QZ algorithm.

In general, the two algorithms return the same result, however using the

QZ algorithm may be more stable for certain problems.

The flag is ignored when A and B are not symmetric.


Matlab eig command1

Example 7.2.3 page 336

>> A = [-7 13 -16;13 -10 13;-16 13 -7]

A =

-7 13 -16

13 -10 13

-16 13 -7

>> [U,D] = eig(A);

>> U

U =

-0.5774 0.4082 0.7071

0.5774 0.8165 -0.0000

-0.5774 0.4082 -0.7071

>> D

D =

-36.0000 0 0

0 3.0000 0

0 0 9.0000

>> A*U

ans =

20.7846 1.2247 6.3640

-20.7846 2.4495 -0.0000

20.7846 1.2247 -6.3640

>> U*D

ans =

20.7846 1.2247 6.3640

-20.7846 2.4495 -0.0000

20.7846 1.2247 -6.3640

Positive definite symmetric matrices
Positive Definite Symmetric Matrices

is [lec4,slide24]

Theorem 4 A symmetric matrix

(semi) positive definite iff all of its eigenvalues

Proof Let

be the orthogonal, diagonal

matrices on the previous page that satisfy

Then for every



is nonsingular

is (semi) positive definite iff


Clearly this condition holds iff

Singular value decomposition
Singular Value Decomposition

then there


Theorem 3 If

exist orthogogonal matrices

has the form

such that


Singular Values

= sqrt eig

Proof Outline Choose



are diagonal, then


try to finish

Matlab svd command

>> help svd

SVD Singular value decomposition.

[U,S,V] = SVD(X) produces a diagonal matrix S, of the same dimension as X and with nonnegative diagonal elements in decreasing order, and unitary matrices U and V so that X =U*S*V'.

S = SVD(X) returns a vector containing the singular values.

[U,S,V] = SVD(X,0) produces the "economy size“ decomposition. If X is m-by-n with m > n, then only the first n columns of U are computed and S is n-by-n.

See also SVDS, GSVD.

Matlab svd command1

>> M = [0 1; 0.5 0.5]

M =

0 1.0000

0.5000 0.5000

>> [U,S,V] = svd(M)

U =

-0.8507 -0.5257

-0.5257 0.8507

S =

1.1441 0

0 0.4370

V =

-0.2298 0.9732

-0.9732 -0.2298

>> U*S*V'

ans =

0.0000 1.0000

0.5000 0.5000

Square roots
Square Roots

Theorem 5 A symmetric positive definite matrix

has a symmetric positive definite ‘square root’.

Proof Let

be the orthogonal, diagonal

matrices on the previous page that satisfy

Then construct the matrices

and observe that

is symmetric positive definite

and satisfies

Polar decomposition
Polar Decomposition

Theorem 6 Every nonsingular matrix

can be factored as


is symmetric

and positive definite and

is orthogonal.

Proof Construct

and observe that

is symmetric and positive definite. Let

be symmetric

and construct

positive definite and satisfy


and clearly

L wdin orthonormalization
Löwdin Orthonormalization

(1) Per-Olov Löwdin, On the Non-Orthogonality Problem Connected with the use of Atomic Wave Functions in the Theory of Molecules and Crystals, J. Chem. Phys. 18, 367-370 (1950).

in an inner product space

Proof Start with

(assumed to be linearly independent), compute the

Gramm matrix


is symmetric and positive definite, Theorem 5

gives (and provides a method to compute) a matrix

that is symmetric and positive definite and

are orthonormal.


The power method pages 340 345
The Power Method pages 340-345

Finds the eigenvalue with largest absolute value of a

whose eigenvalues satisfy


Step 1 Compute a vector with random entries

Step 2 Compute



Step 3 Compute

( recall that


Step 4 Compute






The inverse power method
The Inverse Power Method

Result If

is an eigevector of


corresponding to eigenvalue


is an eigenvector of

corresponding to


Furthermore, if


corresponding to

is an eigenvector of


Definition The inverse power method is the power

method applied to the matrix

It can find the eigenvalue-eigenvector pair if there

is one eigenvalue that has smallest absolute value.

Inverse power method with shifts
Inverse Power Method With Shifts

Computes eigenvalue


closest to

and a corresponding eigenvector

Step 1 Apply 1 or more interations of the power method

to estimate an eigenvalue

using the matrix

- eigenvector pair

- better estimate of

Step 2 Compute

Step 3 Apply 1 or more interations of the power method

using the matrix

to estimate an eigenvalue

and iterate. Then

- eigenvector pair

with cubic rate of convergence !

Unitary and hermitian matrices
Unitary and Hermitian Matrices

Definition The adjoint of a matrix

is the matrix


Definition A matrix

is unitary if

Definition A matrix

is hermitian if

is (semi) positive definite

Definition A matrix


(or self-adjoint)

Super Theorem : All previous theorems true for complex

matrices if orthogonal is replaced by unitary, symmetric by hermitian, and old with new (semi) positive definite.

Homework Due Tutorial 5 (Week 11, 29 Oct – 2 Nov)

1. Do Problem 1 on page 348.

2. Read Convergence of the Power Method (pages 342-346)

and do Problem 16 on page 350.

3. Do problem 19 on pages 350-351.

4. Estimate eigenvalue-eigenvector pairs of the

matrix M using the power and inverse power

methods – use 4 iterations and compute errors

5. Compute the eigenvalue-eigenvector

pairs of the orthogonal matrix O

6. Prove that the vectors

defined at the bottom of

slide 29 are orthonormal by computing their inner products

Extra Fun and Adventure

We have discussed several matrix decompositions :



Singular Value


Find out about other matrix decompositions. How are

they derived / computed ? What are their applications ?