Teorema de Gerschgorin
Teorema de Gerschgorin
Teorema de Gerschgorin
Departamento de Matematicas
Humacao, Puerto Rico 00791
MATE 4061
Analisis Numerico
Dada una matriz A de tama~no n n y entradas complejas, un numero 2 C (conjunto de los numeros
complejos) es un valor propio de A si existe un vector x 6= 0 tal que Ax = x. Usando propiedades de
determinantes, es facil ver que el problema de hallar los valores propios de una matriz es equivalente al de
hallar las ra
ces de un polinomio con coecientes complejos. Por el Teorema Fundamental del Algebra,
tenemos ahora que toda matriz n n tiene exactamente n valores propios contando multiplicidades.
Tambien tenemos, por un teorema famoso de Galois, que el calculo de los valores propios de una matriz
no se puede hacer utilizando unicamente formulas algebraicas. Por tal razon el calculo de los valores
propios de matrices es un problema computacionalmente complejo y el poder estimar los mismos es
de vital importancia. El estimado mas crudo de los valores propios de una matriz esta dado por la
desigualdad (A) kAk, donde
(A) = maxfjj : 2 (A)g
(1)
se conoce como el radio espectral de A. Este estimado aunque util en muchas ocasiones, no es muy
preciso en cuanto a la localizacion de los valores propios de A. El Teorema de Gerschgorin va mas alla
en este sentido. Para A = (a ) denimos los radios
ij
r =
i
X
j
=i
ja j 1 i n
(2)
ij
y los discos
D = fz 2 C : jz ; a j r g:
(3)
El Teorema de Gerschgorin establece que cada valor propio de A pertenece al menos a uno de los D y
que si k de los discos de Gerschgorin se intersecan entre si y a la vez estan aislados de los otros discos,
entonces su union contiene exactamente k de los valores propios de A.
Teorema 1. Sea A una matriz n n y dena los discos D por (2), (3). Entonces
i
ii
2
n
=1
D
(4)
para cualquier valor propio de A. Adem as si S es la uni on de m discos los cuales son disjuntos de los
restantes n ; m, entonces S contiene exactamente m valores propios de A.
Como A y At tienen los mismos valores propios, el Teorema de Gerschgorin tambien es cierto para
los discos columnas :
9
8
<
R = :z 2 C : jz ; a j
i
ii
X
k
=i
6
8 ;2
4
1@;
;1 6 2
A = 16
2 2 ;10
ja j :
ki
(5)
1
A:
Note que kAk = (1=16) maxf14 9 14g = 7=8, de modo que los valores propios de A cumplen con
jj 7=8. Podemos mejorar este estimado con el Teorema de Gerschgorin. Tenemos aplicando (2) que:
1
0.3
D1
0.2
D2
D3
0.1
0.1
0.2
0.3
0.4
1
0.8
0.6
0.4
0.2
0.2
0.4
0.6
Recuerde que un circulo de radio a y centro en (h k) se puede especicar por las ecuaciones parametricas:
ii
en el plano complejo. Usando las formulas de arriba, podemos calcular los discos en MATLAB como
sigue:
t= 0:0.1:6.3]'
X=zeros(length(t),n) Y=X
for i=1:n
X(:,i)=r(i)*cos(t)+real(A(i,i))
Y(:,i)=r(i)*sin(t)+imag(A(i,i))
end
Podemos ahora trazar los valores propios y los discos de Gerschgorin mediante:
e=eig(A)
plot(X,Y,real(e),imag(e),'x')
Combinamos ahora todas estas instrucciones en una funcion que dada la matriz A, esta calcula y dibuja
los discos de Gerschgorin:
function gersch(A)
n=length(A)
diagonal=abs(diag(A))
r=sum(abs(A),2)-diagonal
t= 0:0.1:6.3]'
X=zeros(length(t),n) Y=X
for i=1:n
X(:,i)=r(i)*cos(t)+real(A(i,i))
Y(:,i)=r(i)*sin(t)+imag(A(i,i))
end
e=eig(A)
plot(X,Y,real(e),imag(e),'x')
Puede probar esta funcion trazando los discos de Gerschgorin para una matriz 5 5 arbitraria de entradas
complejas:
A=round(10*rand(5))+sqrt(-1)*round(10*rand(5))
gersch(A)
Ejemplo 2. Dena
B = 3,0.1,20.1,7,22,2,50]
Utilizando la funcion gersch tenemos que los discos Gerschgorin de B quedan como en la Figura (2).
Como B es simetrica, sus valores propios son todos reales y por esto estan en el eje real. Note que
4
10
20
30
40
50
60
lo cual resulta en
C =
3.0000
0.1000
20.0000
0.1000
7.0000
20.0000
0.2000
0.2000
50.0000
La nueva matriz C tiene los mismos valores propios que B ya que son similares. Utilizando C podemos
hallar intervalos que contengan los otros dos valores propios de B . De hecho los discos de C son todos
disjuntos lo que nos da que A tiene un valor propio en cada uno de los intervalos 2:7 3:3], 6:7 7:3] y
con lo que hab
amos calculado antes, otro valor propio en 46 54]. Los discos de C se muestran en la
Figura (3).
40
30
20
10
10
20
30
40
10
20
30
40
50
60
70
80
90