Diagonalización de Matrices Simétricas
Diagonalización de Matrices Simétricas
Diagonalización de Matrices Simétricas
2011
1
0 0 0 n donde 1, 2, .. ,n son los valore propios de A. Si C es una matriz cuyas columnas son vectores propios linealmente independientes de A, entonces D = C-1AC Una matriz diremos que es ortogonal si su transpuesta coincide con su inversa. P orotgonal <=> P-1 = Pt Si P=(u1|u2||un) resulta que decir que P es ortogonal, es equivalente a decir que los vectores {u1,u2,,un} son ortonormales (respecto al producto escalar habitual) Para las matrices reales y
2
simtricas podemos dar una diagonalizacin donde la matriz de paso es ortogonal. Esto es lo que se entiende por diagonalizacin ortogonal.
MATRIZ SIMTRICA
Una matriz de elementos:
es simtrica, si es una matriz cuadrada (m = n) y aij = aji para todo i distinto de j con i, j =1,2,3,4,...,n. Ntese que la simetra es respecto a la diagonal principal. Ejemplo para n = 3:
A es tambin la matriz traspuesta de s misma: At = A. Esta ltima igualdad es una definicin alternativa de matriz simtrica. Las matrices simtricas son un caso particular de las matrices hermticas.
PROPIEDADES
Uno de los teoremas bsicos que concierne este tipo de matrices es el teorema espectral de dimensin finita, que dice que toda matriz simtrica cuyos elementos sean reales es diagonalizable. En particular, es semejante a una matriz ortogonal.
AUTOVALORES
Como las matrices simtricas son un caso particular de las matrices hermticas, todos sus auto valores son reales. Con base en las propiedades de los auto valores de una matriz simtrica, se pueden clasificar en los siguientes tipos:
3
definida positiva: Una matriz simtrica es definida positiva si y solo si todos sus auto valores son estrictamente positivos. definida negativa: Una matriz simtrica es definida negativa si y solo si todos sus auto valores son estrictamente negativos. semidefinida positiva: Una matriz simtrica es semidefinida positiva si y solo si todos sus auto valores son mayores o iguales a cero. semidefinida negativa: Una matriz simtrica es semidefinida negativa si y solo si todos sus auto valores son menores o iguales a cero.
1 0 0 0 0 2 0 0 0 0 3 0
D=
0 n
donde 1, 2, .. ,n son los valore propios de A. Si C es una matriz cuyas columnas son vectores propios linealmente independientes de A, entonces
D = C-1AC
Demostracin Primero se supone que A tiene n vectores propios linealmente independientes V1, V2, , Vn que corresponden a los valores no propios (no necesariamente diferentes) 1, 2, , n.
V1 =
V2 =
, Vn =
Cn1
Cn2
Cnn
Y sea
C11 C21
C12 C22
C1n C2n
Entonces C es invertible ya que sus columnas son linealmente independientes. Ahora bien
A11 A21
AC =
An1
An2
Ann
C1
Y se ve que la columna AC es A =
C2i
Cni
columna i es ivi y
pero
1 0
0 0 2 0
Entonces AC = CD (1)
Y como C es invertible, se pueden multiplicar ambos lados de (1) por la izquierda por C-1 para obtener D = C-1 AC (2)
Esto prueba si A tiene n vectores propios linealmente independientes, entonces A es diagonalizable. Inversamente suponga que A es diagonalizable; esto es, suponga que (2) se cumple para alguna matriz invertible C. Sean v1, v2, vn las columnas de C. entonces AC = 7
CD, e invirtiendo los argumentos anteriores, se ve de inmediato que Avi = ivi para i = 1, 2,, n. Entonces v1, v2, , vn son los vectores propios de A y son linealmente independientes porque C es invertible. Notacin: para indicar que D es la matriz diagonal con componentes diagonales 1, 2, , n, se escribir D = diag (1, 2, , n). Corolario: Si la matriz A de n x n tiene n valores propios diferentes, entonces A es diagonalizable. Observacin: Si se selecciona al azar los coeficientes reales de un polinomio de grado n entonces, con probabilidad 1, el polinomio tendr n races diferentes. No es difcil ver, intuitivamente, por que se cumple esto. Si n = 2, por ejemplo, entonces la ecuacin 2 + a + b = 0 tiene races reales si y solo si a2 = 4b --- un evento muy improbable se a y b se eligen aleatoria mente. Por supuesto, se pueden escribir polinomios que tienen races de multiplicidad algebraica mayor que 1, pero son excepcionales. Por lo tanto, sin pretender precisin matemtica, es posible decir que la mayora de los polinomios tienen races distintas. As, la mayora de las matrices tienen valores propios diferentes y como se estableci al principio, la mayor parte de las matrices son diagonalizables.
donde
CLASIFICACIN:
Las cudricas se clasifican de acuerdo a su signatura, es decir, el mdulo de la diferencia entre el nmero de autovalores positivos y negativos de A00 . Sin embargo, para calcular la signatura de la cudrica no es necesario diagonalizar la matriz. Ello es debido a la existencia de unas cantidades invariantes asociadas a A00 que permiten es. Vemoslo: los autovalores son las races del polinomio caracterstico, es decir, las soluciones de la ecuacin . Ahora bien,
con
Cuando los tres autovalores de A00 son no nulos , es decir det A00 , si escribimos la sucesin K, J, I, 1 y denotamos P y V al nmero de permanencias y variaciones de signo que hay en ella, respectivamente, entonces |P-V| = I, J, K se conocen como invariantes de la cudrica. De esta forma se tiene: 1. Si = 3 : 1. det A > 0 ---> elipsoide real 2. det A < 0 ---> elipsoide imaginario (no existen puntos reales que verifican la ecuacin) 3. det A = 0 ---> cono imaginario 2. Si = 1 : 1. det A > 0 ---> hiperboloide hiperblico (de una hoja) 2. det A < 0 ---> hiperboloide elptico (de dos hojas) 3. det A = 0 ---> cono real Si alguno de los autovalores es nulo (det A00 = 0) pero el determinante de A es distinto de cero, entonces; 1. Si J > 0 ---> paraboloide elptico 2. Si J < 0 ---> paraboloide hiperblico
9
Si det A = det A00 = 0 hay que introducir nuevos invariantes para completar la clasificacin
donde Aii representa la matriz adjunta del elemento aii en A para i=1,2,3. Con estos nuevos invariantes se tiene 1. J > 0 1. K' y signo K' = signo I ---> cilindro elptico imaginario 2. K' y signo K' ---> cilindro elptico real 3. K' = 0 ---> par de planos imaginarios secantes 1. J < 0 1. K' 0 ----> cilindro hiperblico 2. K' = 0 ----> par de planos reales secantes 1. J = 0 y I 1. K' 0 ----> cilindro parablico 2. K' = 0 y J' > 0 -- --> par de planos imaginarios paralelos distintos 3. K' = 0 y J' < 0 -----> par de planos reales paralelos distintos 4. K' = 0 y J' = 0 ----> par de planos coincidentes
Centro:
Plano polar: Dado un punto P = (x0,y0,z0 de matriz A como el plano de ecuacin
3
10
Si P pertenece a la cudrica, entonces el plano polar de P coincide con el plano tangente a dicha superficie en P. No todos los puntos poseen plano polar. La condicin para que un punto (x,y,z) no lo tenga es que verifique el sistema de ecuaciones
que geomtricamente se interpreta como la interseccin de tres planos. Si det A00 entonces el sistema es compatible y tiene solucin nica. El punto solucin se conoce como CENTROde la cudrica. Si det A00 cuando det A=0 y los rangos de ambas matrices son iguales a 2, entonces se dice que la cudrica tiene una recta de centros. Cuando el rango de ambas matrices es igual a 1 hay un plano de soluciones: la cudrica tiene un plano de centros. Finalmente el sistema no tiene solucin si los rangos difieren o det A cudrica carece de centro, recta o plano de centros. As se tiene: 1 Cudricas con centro: elipsoides, hiperboloides y conos. 2 Cudricas con eje de centros: cilindros elpticos e hiperblicos y pares de planos secantes. 3 Cudricas con plano de centros: pares de planos paralelos o coincidentes. 4 El resto de las cudricas no posee centro (lo tiene en el infinito): paraboloides y cilindros parablicos. El centro es un punto de simetra de la cudrica, el eje y el plano de centros son a su vez eje y plano de simetra. Ejemplo: Consideremos la cudrica de ecuacin
Esta cudrica es un elipsoide (vase la tabla de clasificacin). El plano polar por el punto (2,1,3) es el plano de ecuacin
que corta a la superficie (ntese que (2,1,3) es exterior a la superficie (vase la figura).
11
12
En las figuras siguientes vemos los planos polares en los puntos (0,1,1/2) y (0,2,0):
En el primer caso el punto es interior a la superficie y el plano polar es exterior a la misma, mientras que en el segundo caso el punto est sobre el elipsoide y el plano polar coincide con el plano tangente a la superficie en dicho punto.
13
Ecuacin reducida:
La ecuacin reducida de una cudrica es aquella ecuacin simplificada que representa la superficie con su centro (si lo tiene) situado en el origen de coordenadas mientras que los ejes coordenados tienen relaciones particulares con la cudrica. Partiendo de la ecuacin general de una cudrica se puede llegar a su ecuacin reducida aplicandole consecutivamente un giro y una translacin de forma adecuada aunque en algunos casos especiales es necesario aplicar despus de esta ltima un giro plano. A continuacin recogemos los tipos de ecuaciones reducidas y que cudricas representan as como la forma de obtenerlas a partir de los invariantes.
Denotemos por
las races de
entonces:
elipsoide
hiperboloide hiperblico
14
Paraboloides:
donde
donde
15
cilindro elptico
cilindro hiperblico
Cilindro parablico:
donde
16
cilindro parablico
donde
Cudricas no degeneradas:
17
Elipsoide
Ecuacin reducida:
La cudrica tiene signatura 3 y los autovalores de la matriz A00 son los tres positivos. Los cortes del elipsoide por planos paralelos a los coordenados son curvas cnicas de tipo elipse (en lo siguiente se supone que el elipsoide esta centrado en el origen de coordenadas y tiene la ecuacin reducida que se da arriba): 1. por planos z
18
Si c no hay interseccin real, mientras que si c la interseccin se reduce a un punto siendo el plano tangente a la superficie elptica.
Para los planos de la forma y = o x= el resultado es anlogo al anterior intercambiando el papel de las variables de forma adecuada.
19
20
21
La cudrica tiene signatura 1 y los autovalores de la matriz A00 son dos positivos y uno negativo. Los cortes del hiperboloide por planos paralelos a los coordenados son curvas cnicas: (en lo siguiente se supone que el hiperboloide esta centrado en el origen de coordenadas y tiene la pimera de las ecuaciones reducidas dadas arriba):
donde
22
23
(elipse de garganta )
por planos x=
24
por planos y=
El corte es una hiprbola como la del caso anterior donde los papeles de x e y se han intercambiado
El hiperboloide hiperblico puede ser visto tambin como una superficie de revolucin engendrada al girar una hiperbola alrededor del eje de la cudrica ( en el caso de la ecuacin reducida que estamos utilizando, el eje z) describiendo una elipse. Adems el hiperboloide hiperblico es una superficie doblemente reglada puesto que contiene a las dos familias de rectas. Veamoslo. La ecuacin del hiperboloide se puede escribir como
Entonces cualquier punto que satisface la ecuacin del hiperboloide satisface el siguiente conjunto de ecuaciones para algun valor del parametro.
Cada una de las ecuaciones anteriores representa un plano luego finalmente tenemos un par de rectas contenidas en el hiperboloide.
25
Hiperboloide elptico
Ecuacin reducida:
La cudrica tiene signatura 1 y los autovalores de la matriz A00 son dos negativos y uno positivo. Los cortes del hiperboloide por planos paralelos a los coordenados son curvas cnicas: (El desarrollo que sigue se ha hecho utilizando la primera de las ecuaciones reducidas)
26
por planos x= si |a entonces la curva interseccin resulta ser una elipse de ecuacin
donde
27
( los planos: x=a y x=- -a si |a no hay interseccin real mientras que si |=a entonces la interseccin se reduce a un punto y el plano en cuestin es tangente a la superficie.
28
( a= 0 ) A continuacin incluimos otros dibujos en los cuales el hiperboloide tiene la segunda ecuacin reducida y los parmetros a,b y c son distintos
29
30
31
Paraboloide elptico
Ecuacin reducida:
32
Los cortes del paraboloide por planos paralelos a los coordenados son curvas cnicas (en lo siguiente se supone que el parabolide tiene la ecuacin reducida que se da arriba):
por planos z si entonces la curva interseccin resulta ser una elipse de semiejes a y b y ecuacin
33
Si entonces no existe interseccin mientras que para la interseccin se reduce a un punto siendo la superficie cudrica tangente al plano en dicho punto.
(a<0)
por planos y = o por planos x= las curvas interseccin son las parbolas
34
Paraboloide hiperblico
35
Ecuacin reducida:
36
En lo que sigue utilizaremos la primera de las ecuaciones reducidas. El paraboloide hiperblico es una superficie doblemente reglada por las familias de rectas:
Los cortes del paraboloide por planos paralelos a los coordenados son curvas cnicas:
37
( > 0)
( < 0)
38
39
parbolas
y respectivamente.
( y = 0)
40
(x = 0)
A continuacin presentamos figuras donde se ha cortado el paraboloide hiperblico por plano oblicuos no paralelos a los coordenados
41
IDENTIFICACIN DE CNICAS
La circunferencia, la elipse, la parbola o la hiprbola son curvas planas de todos conocidas. Estas curvas aparecan ya en la geometra griega y fueron denominadas secciones cnicas, ya que los griegos de la poca de Platn consideraban que tales curvas procedan de la interseccin de un cono con un plano.
42
Finalmente se estableci una teora algebraica general que engloba todas estas curvas y las describe como curvas cuadrticas. Es esta teora la que presentamos a continuacin.
43
donde
Una cnica queda pues definida por una matriz simtrica En lo que sigue denotaremos por Aii a la matriz adjunta en A del elemento aii i=0,1,2 .
44
Ejemplo:
Las figuras que represetan las ecuaciones cuadrticas pueden ser, adems de elipses, hiprbolas y parbolas, pares de rectas tanto secantes como paralelas y estas ltimas pueden ser distintas o coincidentes. Tambin puede darse el caso de que la ecuacin sea verificada por un nico punto o por ninguno. Alguna de estas ltimas tambin se pueden obtener como secciones cnicas como se ve en las imgenes siguientes:
45
A continuacin estudiamos como se puede determinar que tipo de curva que define una ecuacin cuadrtica dada.
Si y son las matrices asociadas a la cnica despus de que sta ha sufrido un giro y una traslacin, respectivamente, entonces 1) det A=det A'=det A'', 2) a11 + a22 = a'11+ a'22 = a''11 + a''22, 3) det A00 = det A'00 = det A''00. 46
Tabla de Clasificacin
det A 0
det A00 0
Elipse
Elipse real
Hiprbola
det A00 = 0
Parbola
det A11 + det A22 > 0 imaginarias det A11 + det A22 < 0 reales
Rectas paralelas
Rectas paralelas
Rectas coincidentes
Centro:
Polar Dado un punto P=(x0,y0) se llama polar de P respecto de una cnica C de matriz A a la recta de
ecuacin
47
Si el punto P est en la cnica C entonces la recta polar de P respecto a C es precisamente la recta tangente a la cnica en dicho punto P.
Ejemplo:
Consideremos la cnica de ecuacin
Utilizando la tabla de clasificacin vemos que se trata de una elipse real puesto que
Observamos que el punto (1,2) pertenece a la cnica y por lo tanto la polar coincide con la recta tangente en dicho punto.
48
49
Es posible que un punto P=(x0,y0) no tenga polar respecto a una cnica C. Las coordenadas de los puntos que no tienen polar deben verificar el sistema de ecuaciones
que impone que los coeficientes de x e y en la recta polar sean nulos (es decir que no exista dicha recta). Para que este sistema tenga alguna solucin se ha de verificar
50
Adems si det A00 es no nulo, entonces la solucin del sistema es nica y por lo tanto habra un nico punto que no poseer recta polar. Este punto se denomina centro de la cnica. No todas las cnicas tienen centro. El centro de la cnica tiene la particularidad de ser su centro de simetra. Si C es una elipse o una hiprbola entonces det A00 0 y el sistema es compatible determinado lo que indica que estas cnicas tienen centro y que ste es nico.
Ejemplo: En la elipse del ejemplo anterior el centro ser el punto (2,2) nica solucin del sistema de ecuaciones
Sin embargo si la cnica es una parbola, todos sus puntos tienen polar. La parbola es por tanto una cnica sin centro.
Polo Dada una recta r diremos que un punto P es un polo de r respecto a una cnica C si r es la polar
de P respecto a C
51
Por supuesto hay rectas que no tienen polos: en la elipse y la hiprbola son todas las que pasan por el centro y en la parbola son las rectas de direccin (-a12, a11) que son adems perpendiculares al vector (a01 a12 - a02 a11, a01 a22 - a02 a12) y las nicas que cortan a la parbola en un solo punto.
Diremos que dos dimetros son conjugados si no son asntotas (en el caso de la hiprbola) y uno de ellos coincide con el lugar geomtrico de los puntos medios de las cuerdas determinadas por la cnica en las rectas paralelas al otro.
52
Ecuacin reducida
La ecuacin reducida de una cnica es aquella ecuacin simplificada de la curva que sita el centro (si lo tiene) de la cnica como origen de coordenadas mientras que los ejes presentan unas relaciones particulares con la cnica. Partiendo de la ecuacin general de una cnica se puede llegar a su ecuacin reducida aplicndole consecutivamente un giro y una traslacin de forma adecuada. Clasificaremos en tres tipos las ecuaciones reducidas de las cnicas:
53
Parbola
54
con
con
55
56
57
Si F es un punto fijo del plano y D una recta, el lugar geomtrico de los puntos del plano cuyas distancias al punto F y a la recta D estn en proporcin constante es una cnica no degenerada (elipse, hiprbola, parbola). Al punto F se le denomina foco de la cnica y a la recta D directriz asociada al foco F.
Ecuacin focal
Si C es una cnica propia (no degenerada), en un sistema de referencia determinado, su ecuacin ser
donde F=(x0, y0) es el foco de la cnica y la recta ax+by+d=0 es la directriz asociada al foco F. La ecuacin anterior se puede transformar fcilmente en
Renombrando los coeficientes y haciendo a2+b2=e2 ya que a2+b2 >0 la ecuacin queda finalmente como
Esta ecuacin se conoce como ecuacin focal de la cnica y a'x+b'y+d'=0 es la ecuacin de Euler para la directriz. Ms adelante veremos como se puede obtener la ecuacin focal partiendo de la ecuacin reducida de la cnica. La cantidad e se denomina excentricidad de la cnica y es un invariante puesto que, como se comprueba fcilmente, det A00=1-e2 En trminos de distancia, la ecuacin focal queda como
Es decir, en una cnica no degenerada el cociente entre las distancias de cualquiera de sus puntos al foco y a la directriz es constante (dada por la excentricidad).
58
Se puede comprobar que en una cnica la directriz asociada a un foco F es precisamente su recta polar. El foco de la cnica es pues cualquier punto del plano para el que la razn de la distancia de un punto cualquiera P de la cnica a F y la distancia de P a la polar de F es constante
59
60
Cnicas no degeneradas
Cnicas con centro
Como hemos visto las cnicas propias (o no degeneradas) con centro son la elipse, dentro de la cual se incluye la circunferencia como caso particular (cuando a11= a22), y la hiprbola. Antes de estudiar cada una de ellas por separado veamos algunas de las caractersticas conjuntas que presentan.
Para estas cnicas se definen los ejes de la cnica como un par de dimetros conjugados perpendiculares que son adems ejes de simetra de la cnica. Adems slo existe un par de ejes salvo en el caso de la circunferencia en el que hay infinitos.
Aviso: esta seccin requiere conocimientos de lgebra lineal de la carrera. Si quieres saltarla pincha aqu
A continuacin vemos como se pueden calcular los ejes. La matriz A00 asociada a la cnica C define un endomorfismo de IR2. Los valores propios de este endomorfismo, 1 y 2, tendrn vectores propios asociados, v1 y v2, ortonormales. Como 1 y 2 son las races del polinomio caracterstico de A00, si el discriminante de la ecuacin en l det(A00-I) = 0 es no nulo entonces hay dos soluciones distintas y v1 y v2 quedan determinados de forma nica (salvo cambio de orden y signo) y reciben el nombre de direcciones principales de la cnica. Si el discriminante es nulo entonces l1=l2 y cualquier vector es propio. Se dice entonces que toda direccin es principal. Como v1 y v2 se pueden elegir ahora cualquier par de vectores ortonormales. Si C es el centro de la cnica y v1 y v2 sus direcciones principales entonces las rectas que pasan por el centro y tienen como vectores directores a v1 y v2 son los ejes de la cnica.
61
Su matriz asociada es
Como det A=-1/4 0 y det A00 = -9/4 < 0, la cnica es una hiprbola. El centro de la hiprbola es la solucin del sistema
es decir el punto (2/9,-10/9). Para calcular los ejes de la hiprbola hallamos en primer lugar los autovalores de A00 o lo que es lo mismo las races del polinomio caracterstico de A00 es decir necesitamos resolver la ecuacin en l det(A00-I)=0, que en este caso queda
que sern las direcciones de los ejes de la hiprbola. Teniendo en cuenta que los ejes pasan por el centro ya calculado, sus ecuaciones sern
62
63
64
Supongamos que los focos son F1=(c,0) y F2=(-c,0) y llamemos 2a a la suma de las distancias, entonces los puntos (x,y) de la elipse verifican
Esta es la ecuacin reducida de la elipse en la que los ejes coordenados son los ejes de simetra de la elipse y el origen de coordenadas es su centro. Encontremos la ecuacin focal de la elipse
65
En esta ecuacin focal tenemos que el foco es el punto (c,0) y la directriz es la recta paralela al eje X: x=a2/c. La excentricidad es e=c/a que es estrictamente menor que 1 puesto que c=(a2-b2)1/2< a. Una cnica propia es una elipse si la excentricidad es menor que 1: e< 1. Cuando e=0 la elipse es una circunferencia: la excentricidad en la elipse mide, por tanto, lo que sta se aleja de la circularidad. En la circunferencia los dos focos se confunden y son a su vez el centro de la cnica. En la animacin siguiente se ve como vara la elipse al ir disminuyendo su excentricidad.
66
67
hecho, ni la Tierra ni el cometa seran visibles en la animacin en el caso de que se representasen a escala.
Supongamos que los focos son F1=(c,0) y F2=(-c,0) y llamemos 2a a la diferencia de las distancias, entonces los puntos (x,y) de la hiprbola verifican (c>a)
68
De igual forma a como se hizo para la elipse, se consigue la ecuacin focal de la hiprbola
pero ahora c2=a2+b2 La excentricidad en la hiprbola es e=c/a >1 Luego las cnicas no degeneradas de excentricidad mayor que 1 son hiprbolas.
69
Supongamos que el foco es F=(p/2,0) con p>0 y la recta directriz es x=-p/2, entonces los puntos (x,y) de la parbola verifican la ecuacin
La excentricidad en la parbola es 1. 71
72