Scientific Computing, Lecture Notes - Computer Science - 2, Study notes of Computer Numerical Control

Eigenvalue problem, eigenvector, determinant, hermitian matrices, eigen decomposition, rayleigh quotient scaling

Typology: Study notes

2010/2011

Uploaded on 09/09/2011

jennyfromtheblock
jennyfromtheblock 🇬🇧

2.3

(3)

225 documents

1 / 35

Toggle sidebar

This page cannot be seen from the preview

Don't miss anything!

bg1
Eigenvalue Problems
Solving linear systems Ax= bis one part
of numerical linear algebra, and involves
manipulating the rows of a matrix.
The second main part of numerical linear
algebra is about transforming a matrix to
leave its eigenvalues unchanged.
Ax= lx
where lis an eigenvalue of A and non-
zero x is the corresponding eigenvector.
pf3
pf4
pf5
pf8
pf9
pfa
pfd
pfe
pff
pf12
pf13
pf14
pf15
pf16
pf17
pf18
pf19
pf1a
pf1b
pf1c
pf1d
pf1e
pf1f
pf20
pf21
pf22
pf23

Partial preview of the text

Download Scientific Computing, Lecture Notes - Computer Science - 2 and more Study notes Computer Numerical Control in PDF only on Docsity!

Eigenvalue Problems

  • Solving linear systems A

x^ =

b

is one part

of numerical linear algebra, and involvesmanipulating the rows of a matrix. • The second main part of numerical linearalgebra is about transforming a matrix toleave its eigenvalues unchanged.

A x

l x

where

l^

is an

eigenvalue

of A and non-

zero

x^

is the corresponding

eigenvector

What are Eigenvalues?

  • Eigenvalues are important in physical,biological, and financial problems (andothers) that can be represented byordinary differential equations. • Eigenvalues often represent things likenatural frequencies of vibration, energylevels in quantum mechanical systems,stability parameters.

Determinants

-^ Suppose we if eliminate the components of

x

from A

x =

0 using Gaussian Elimination without

pivoting. • We do the kth forward elimination step bysubtracting a

times row k from ajk

times row j,kk^

for j=k,k+1,…n. • Then at the end of forward elimination we haveU x

=^0

, and U

nn^

is the determinant of A, det(A).

-^ For nonzero solutions to A

x =

0 we must have

det(A) = 0. • Determinants are defined only for squarematrices.

Determinant Example

-^ Suppose we have a^3

μ3 matrix.

⎞⎟ ⎟ ⎟⎠

⎛⎜ =⎜ ⎜⎝

33 32 31

23 22 21

13 12 11

a a a

a a a

a a a A

-^ So A

x =

0 is the same as: a^11

x^ +a^1

x 12

+a 2

x 13

= 0 3

a^21

x^ +a^1

x 22

+a 2

x 23

= 0 3

a^31

x^ +a^1

x 32

+a 2

x 33

= 0 3

Determinant Example (continued)

-^ Step k=2:

) times equation 3. and so:

  • a – subtract (a
  • -a
  • a
  • ) times equation 2 from (a
  • a
    • a
    • a
      • [(a • So we have:
      • a-a
      • a
      • )(a
      • a-a
      • a
      • a )- (a
      • -a
      • a)(a
      • a
      • -a
      • a)]x
      • =
        • [a which becomes:
        • (a
        • a
        • –a
        • a) – a
        • (a
        • a
        • -a
        • a) + a
        • (a
        • a
        • -a
        • a)]x
        • =
          • a det(A) =
          • (a
          • a
          • –a
          • a) – a
          • (a
          • a
          • -a
          • a) + a
          • (a
          • a
          • -a

Definitions

-^ Minor M

of matrix A is the determinant of theij^

matrix obtained by removing row i and column jfrom A. • Cofactor C

= (-1)ij^

i+j M

ij

-^ If A is a 1

μ1 matrix then det(A)=a

. 11 -^ In general,

ij n j

Caij

A^

=^ ∑=^1

det(

where i can be any value i=1,…n.

Characteristic Equation

  • A

x^ =

l x

can be written as^ (A-

lI)

x^ =

^0

which holds for x

∫0, so (A-

lI) is singular

and

det(A-

lI) = 0

  • This is called the characteristic polynomial.If A is n

μn the polynomial is of degree n

and so A has n eigenvalues, countingmultiplicities.

Example

⎞ ⎟ ⎟ ⎠ ⎛ ⎜ = ⎜ ⎝

2 1

3 4

a
a

⎞ ⎟⎟ ⎠

⎛ ⎜⎜ ⎝

− − = −

λ λ

λ^

2 1

3

4 I A

0 ) 3 )( (^1) ( ) (^2) )( (^4) (

0 )

det(

= − − − ⇒ = −

λ I

A

2

λ • Hence the two eigenvalues are 1 and 5.

Positive Definite Matrices

  • A complex matrix A is

positive definite

if

for every nonzero complex vector

x^

the

quadratic form

H x

A x

is real and:

H x

A x

where

H x

denotes the conjugate transpose

of^

x^ (i.e., change the sign of the imaginary

part of each component of

x^

and then

transpose).

Eigenvalues of Positive Definite

Matrices

  • If A is positive definite and

l^

and

x^

are an

eigenvalue/eigenvector pair, then:

A x

l x

fl

H x

A x

l x

H x

  • Since

H x

A x

and

H x

x^ are both real and

positive it follows that

l^

is real and

positive.

Hermitian Matrices

-^ A square matrix for which A = A

H^ is said to be an

Hermitian

matrix.

-^ If A is real and Hermitian it is said to be^ symmetric

, and A = A

T^.

-^ Every Hermitian matrix is positive definite. •^ Every eigenvalue of an Hermitian matrix is real. •^ Different eigenvectors of an Hermitian matrix areorthogonal to each other, i.e., their scalarproduct is zero.

Eigen Decomposition

-^ Let

l^1

,l^2

,…,

ln^

be the eigenvalues of the n

μn

matrix A and

x^1

, x^2

,…,

x n^

the corresponding

eigenvectors. • Let

L be the diagonal matrix with

l^1

,l^2

,…,

ln^

on

the main diagonal. • Let X be the n

μn matrix whose jth column is

x .j

-^ Then AX = X

L, and so we have the

eigen

decomposition

of A:^ A = X

LX

  • -^ This requires X to be invertible, thus theeigenvectors of A must be linearly independent.

The Power Method

  • Label the eigenvalues in order ofdecreasing absolute value so |

l^1

l^2

|ln

  • Consider the iteration formula:

y k+

= A

y k

where we start with some initial

y^0

, so that:

y k

= A

k y

0

  • Then

y k

converges to the eigenvector

x^1

corresponding the eigenvalue

l^1

Proof

-^ We know that A

k^ = X

kL

-1X , so:

y k^

= A

k y^0

= X

kL -1X y^0

-^ Now we have:

⎞ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎠

⎛ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎝ ⎞⎟ ⎟^ =⎟ ⎟ ⎟ ⎠

⎛ ⎜ ⎜ =⎜ ⎜ ⎜ ⎝ Λ

k n^ k

k k k k n

k k k

1

2 1 1

2 1

1

λ^ λ

λ λ λ λ

λ λ

O

O

-^ The terms on the diagonal get smaller inabsolute value as k increases, since

l^1

is the

dominant eigenvalue.