06.09.2021 Views

Linear Algebra, Theory And Applications, 2012a

Linear Algebra, Theory And Applications, 2012a

Linear Algebra, Theory And Applications, 2012a

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

15.1. THE POWER METHOD FOR EIGENVALUES 385<br />

Now pick the eigenvalue λ q which is closest to q. Then<br />

|Ax − qx| 2 =<br />

n∑<br />

∑<br />

n<br />

|a k | 2 (λ k − q) 2 ≥ (λ q − q) 2 |a k | 2 =(λ q − q) 2 |x| 2<br />

k=1<br />

which implies (15.6). <br />

Example 15.1.10 Consider the symmetric matrix A =<br />

⎛<br />

⎝ 1 2 2 2 3<br />

1<br />

⎞<br />

⎠ . Let x =(1, 1, 1) T .<br />

3 1 4<br />

How close is the Rayleigh quotient to some eigenvalue of A? Find the eigenvector and eigenvalue<br />

to several decimal places.<br />

Everything is real and so there is no need to worry about taking conjugates. Therefore,<br />

the Rayleigh quotient is<br />

⎛<br />

( ) 1 1 1 ⎝ 1 2 3 ⎞ ⎛<br />

2 2 1 ⎠ ⎝ 1 ⎞<br />

1 ⎠<br />

3 1 4 1<br />

= 19<br />

3<br />

3<br />

According to the above theorem, there is some eigenvalue of this matrix λ q such that<br />

∣⎛<br />

⎞ ⎛ ⎞ ⎛ ⎞∣<br />

∣ λ q − 19<br />

3 ∣<br />

⎝ 1 2 3<br />

2 2 1 ⎠ ⎝ 1 1<br />

∣ 3 1 4 1<br />

≤<br />

√<br />

3<br />

⎛ ⎞<br />

= √ 1 − 1 3<br />

⎝ − 4 ⎠<br />

3<br />

3 5<br />

3<br />

=<br />

√<br />

1<br />

) 2<br />

k=1<br />

⎠ − 19<br />

3<br />

9 + ( 4 2 (<br />

3)<br />

+<br />

5<br />

3<br />

√ =1. 247 2<br />

3<br />

⎝ 1 1<br />

1<br />

Could you find this eigenvalue and associated eigenvector? Of course you could. This is<br />

what the shifted inverse power method is all about.<br />

Solve ⎛⎛<br />

⎝⎝ 1 2 3<br />

⎞ ⎛<br />

2 2 1 ⎠ − 19 ⎝ 1 0 0<br />

⎞⎞<br />

⎛<br />

0 1 0 ⎠⎠<br />

⎝ x ⎞ ⎛<br />

y ⎠ = ⎝ 1 ⎞<br />

1 ⎠<br />

3<br />

3 1 4<br />

0 0 1 z 1<br />

In other words solve ⎛<br />

⎞ ⎛<br />

− 16<br />

3<br />

2 3<br />

⎝ 2 − 13 3<br />

1 ⎠ ⎝ x ⎞<br />

y ⎠ =<br />

3 1 − 7 z<br />

3<br />

and divide by the entry which is largest, 3. 870 7, to get<br />

⎛<br />

u 2 = ⎝<br />

. 699 25<br />

. 493 89<br />

1.0<br />

⎞<br />

⎠<br />

⎛<br />

⎝ 1 1<br />

1<br />

⎞<br />

⎠<br />

⎠<br />

∣<br />

Now solve<br />

⎛<br />

⎝<br />

− 16 3<br />

2 3<br />

2 − 13<br />

3<br />

1<br />

3 1 − 7 3<br />

⎞ ⎛<br />

⎠<br />

⎝ x y<br />

z<br />

⎞ ⎛<br />

⎠ = ⎝<br />

. 699 25<br />

. 493 89<br />

1.0<br />

⎞<br />

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!