Top Banner
Lecture 16 Cramer’s Rule, Eigenvalue and Eigenvector Shang-Hua Teng
28

Lecture 16 Cramer’s Rule, Eigenvalue and Eigenvector Shang-Hua Teng.

Dec 22, 2015

Download

Documents

Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Page 1: Lecture 16 Cramer’s Rule, Eigenvalue and Eigenvector Shang-Hua Teng.

Lecture 16Cramer’s Rule, Eigenvalue and

Eigenvector

Shang-Hua Teng

Page 2: Lecture 16 Cramer’s Rule, Eigenvalue and Eigenvector Shang-Hua Teng.

Determinants and Linear SystemCramer’s Rule

333231

232221

131211

33323

23222

13121

1

33323

23222

13121

3

2

1

333231

232221

131211

3

2

1

3

2

1

333231

232221

131211

det

det

10

01

00

aaa

aaa

aaa

aab

aab

aab

x

aab

aab

aab

x

x

x

aaa

aaa

aaa

b

b

b

x

x

x

aaa

aaa

aaa

Page 3: Lecture 16 Cramer’s Rule, Eigenvalue and Eigenvector Shang-Hua Teng.

Cramer’s Rule

• If det A is not zero, then Ax = b has the unique solution

ni

A

aabaax niii

,...,2,1

det

],,,,,det[ 111

Page 4: Lecture 16 Cramer’s Rule, Eigenvalue and Eigenvector Shang-Hua Teng.

Cramer’s Rule for Inverse

nji

A

CA ji

ij

,...,2,1,

det1

A

aaeaaA niji

ij det

],,,,,det[ 1111 Proof:

Page 5: Lecture 16 Cramer’s Rule, Eigenvalue and Eigenvector Shang-Hua Teng.

Where Does Matrices Come From?

Page 6: Lecture 16 Cramer’s Rule, Eigenvalue and Eigenvector Shang-Hua Teng.

Computer Science

• Graphs: G = (V,E)

Page 10: Lecture 16 Cramer’s Rule, Eigenvalue and Eigenvector Shang-Hua Teng.

Graphs in Scientific Computing

Page 11: Lecture 16 Cramer’s Rule, Eigenvalue and Eigenvector Shang-Hua Teng.

Resource Allocation Graph

Page 12: Lecture 16 Cramer’s Rule, Eigenvalue and Eigenvector Shang-Hua Teng.

Road Map

Page 13: Lecture 16 Cramer’s Rule, Eigenvalue and Eigenvector Shang-Hua Teng.

Matrices Representation of graphs

Adjacency matrix: ( ) , # edgesij ijA a a ij

Page 14: Lecture 16 Cramer’s Rule, Eigenvalue and Eigenvector Shang-Hua Teng.

Adjacency Matrix:

01001

10100

01011

00101

10110

A

1

2

34

5

edgean not is j)(i, if 0

edgean is j)(i, if 1ijA

Page 15: Lecture 16 Cramer’s Rule, Eigenvalue and Eigenvector Shang-Hua Teng.

Matrix of GraphsAdjacency Matrix:• If A(i, j) = 1: edge exists

Else A(i, j) = 0.

0101

0000

1100

00101 2

34

1

-3

3

2 4

Page 16: Lecture 16 Cramer’s Rule, Eigenvalue and Eigenvector Shang-Hua Teng.

01001

12100

01311

00121

10113

L

1

2

34

5

Laplacian of Graphs

Page 17: Lecture 16 Cramer’s Rule, Eigenvalue and Eigenvector Shang-Hua Teng.

Matrix of Weighted GraphsWeighted Matrix:• If A(i, j) = w(i,j): edge exists

Else A(i, j) = infty.

032

0

430

101 2

34

1

-3

3

2 4

Page 18: Lecture 16 Cramer’s Rule, Eigenvalue and Eigenvector Shang-Hua Teng.

Random walks

How long does it take to get completely lost?

Page 19: Lecture 16 Cramer’s Rule, Eigenvalue and Eigenvector Shang-Hua Teng.

Random walks Transition Matrix

1

2

345

6

02

1

4

100

2

13

10

4

1000

3

1

2

10

2

1

3

10

004

10

3

10

004

1

2

10

2

13

1000

3

10

P

Page 20: Lecture 16 Cramer’s Rule, Eigenvalue and Eigenvector Shang-Hua Teng.

Markov Matrix

• Every entry is non-negative

• Every column adds to 1

• A Markov matrix defines a Markov chain

Page 21: Lecture 16 Cramer’s Rule, Eigenvalue and Eigenvector Shang-Hua Teng.

Other Matrices

• Projections

• Rotations

• Permutations

• Reflections

Page 22: Lecture 16 Cramer’s Rule, Eigenvalue and Eigenvector Shang-Hua Teng.

Term-Document Matrix• Index each document (by human or by

computer)– fij counts, frequencies, weights, etc

m term

2 term

1 term

n docdoc21 doc

21

22221

1 1211

mnmm

n

n

fff

fff

fff

• Each document can be regarded as a point in m dimensions

Page 23: Lecture 16 Cramer’s Rule, Eigenvalue and Eigenvector Shang-Hua Teng.

Document-Term Matrix• Index each document (by human or by

computer)– fij counts, frequencies, weights, etc

m doc

2 doc

1 doc

n term2 term1 term

21

22221

1 1211

mnmm

n

n

fff

fff

fff

• Each document can be regarded as a point in n dimensions

Page 24: Lecture 16 Cramer’s Rule, Eigenvalue and Eigenvector Shang-Hua Teng.

Term Occurrence Matrix

Page 25: Lecture 16 Cramer’s Rule, Eigenvalue and Eigenvector Shang-Hua Teng.

c1 c2 c3 c4 c5 m1 m2 m3 m4 human 1 0 0 1 0 0 0 0 0 interface 1 0 1 0 0 0 0 0 0 computer 1 1 0 0 0 0 0 0 0 user 0 1 1 0 1 0 0 0 0 system 0 1 1 2 0 0 0 0 0 response 0 1 0 0 1 0 0 0 0 time 0 1 0 0 1 0 0 0 0 EPS 0 0 1 1 0 0 0 0 0 survey 0 1 0 0 0 0 0 0 1 trees 0 0 0 0 0 1 1 1 0 graph 0 0 0 0 0 0 1 1 1 minors 0 0 0 0 0 0 0 1 1

Page 26: Lecture 16 Cramer’s Rule, Eigenvalue and Eigenvector Shang-Hua Teng.

Matrix in Image Processing

Page 27: Lecture 16 Cramer’s Rule, Eigenvalue and Eigenvector Shang-Hua Teng.

Random walks

How long does it take to get completely lost?

0

0

0

0

0

1

Page 28: Lecture 16 Cramer’s Rule, Eigenvalue and Eigenvector Shang-Hua Teng.

Random walks Transition Matrix

1

2

345

6

0

0

0

0

0

1

02

1

4

100

2

13

10

4

1000

3

1

2

10

2

1

3

10

004

10

3

10

004

1

2

10

2

13

1000

3

10

100

P