11.08.2013 Views

Apunte

Apunte

Apunte

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

2. MÉTODOS ITERATIVOS 57<br />

Demostración. Los λ son autovalores de BGS y por lo anterior verifican que<br />

det(λ(L + D) + U) = 0<br />

pero λ(L + D) + U es tridiagonal y por el lema anterior<br />

Si λ = 0 sea α tal que α2 = 1<br />

λ . Entonces,<br />

det(λD + αλL + α −1 U) = 0<br />

0 = α −N det(λαD + L + U)<br />

y como los autovalores de BJ no nulos son las raíces µ de det(µD + L + U) = 0 resulta que<br />

Pero como α2 1<br />

λ<br />

se tiene que<br />

µ = λα<br />

µ 2 = λ.<br />

Ahora observemos que en lo anterior, dado λ = 0 autovalor de BGS encontramos µ autovalor de<br />

BJ con µ 2 = λ, pero en realidad es un si y sólo si. Es decir, dado µ autovalor de BJ, λ = µ 2<br />

resulta ser un autovalor de BGS,<br />

si y sólo si<br />

si y sólo si (por el lema previo)<br />

y tomando α = µ se tiene<br />

Entonces µ 2 es autovalor de BGS.<br />

det(µD + L + U) = 0<br />

det(µ 2 D + µ(L + U)) = 0<br />

det(µ 2 D + αµL + α −1 µU) = 0<br />

det(µ 2 (D + L) + U)) = 0.<br />

Convergencia del método de Gauss-Seidel para A ∈ IR N×N simétrica<br />

Como los autovalores de A ∈ IR N×N pueden ser complejos, trabajaremos directamente con<br />

A ∈ C N×N .<br />

Recordemos algunas definiciones<br />

Definición 3.20. Si A ∈ CN×N se define A∗ = AT o sea A∗ es la matriz que en el lugar i, j<br />

tiene al elemento a∗ ij = aji.<br />

Definición 3.21. A ∈ C N×N es Hermitiana si<br />

A ∗ = A.

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

Saved successfully!

Ooh no, something went wrong!