Following is the example from this book.
My question is, what value of λ did he put in? He did not explain that, can anybody explain how did he get 1.24?
Asked
Active
Viewed 2,188 times
0
user4904589
- 365
1 Answers
1
For the matrix $T_j$, he calculates the eigenvalues. We have the characteristic polynomial:
$$\tag 1 \det(T_j - \lambda I) = -\lambda(\lambda^2 -0.625) \implies \lambda_{1,2,3} = 0, \pm \sqrt{0.625}$$
The spectral radius is given by the maximum magnitude eigenvalue:
$$\rho(A) = \displaystyle \max ~ |\lambda_i| $$
From $(1)$, we have:
$$\rho(T_j) = \displaystyle \max ~ |\lambda_i| = \sqrt{0.625}$$
Now, when $A$ is positive definite and tridiagonal, we have:
$$\omega = \dfrac{2}{1 + \sqrt{1 - [\rho(T_j)]^2}}$$
The calculation follows.
There is a write-up of the theory of this choice in Section 3.
Moo
- 11,311
-
But what if the spectral radius is greater than 1? The expression inside the square root will become negative and hence $\omega$ becomes imaginary? – Aditi Narware Nov 16 '17 at 13:39
-
@AditiNarware: See https://math.stackexchange.com/questions/2390456/matrix-with-spectral-radius-greater-or-equal-1-such-that-fixed-point-iteration-c – Moo Nov 16 '17 at 15:23
-
@AditiNarware: Also see page $7$ of http://www.it.uom.gr/teaching/linearalgebra/ExamplesToIterativeMethods.pdf – Moo Nov 16 '17 at 15:31

