Diagonalización

Plantilla:Referencias En matemáticas y, en particular, en álgebra lineal, la diagonalización es un proceso que permite simplificar la descripción de ciertos endomorfismos de un espacio vectorial. En particular, identificando el endomorfismo con su matriz asociada en cierta base, se puede hablar de diagonalización de matrices. Consiste en encontrar una base del espacio vectorial formada por vectores propios, si existe alguna. Esto se refleja en obtener una base tal que la matriz asociada al endomorfismo en la misma es una matriz diagonal.

El proceso se reduce, pues, a una reducción máxima del endomorfismo, es decir, a una descomposición del espacio vectorial en suma directa de subespacios vectoriales invariantes por el endomorfismo. Restringido sobre cada uno de ellos, el endomorfismo se reduce a una homotecia. Por tanto, la diagonalización permite una mejor visualización geométrico de la actuación del endomorfismo sobre el espacio vectorial. Además, la diagonalización permite un cálculo rápido y simple de potencias y exponenciales de matrices (entendidas como matrices de un endomorfismo), lo que permite expresar numéricamente ciertos sistemas dinámicos lineales, obtenidos por iteración o por ecuaciones diferenciales.

Introducción y deducción del método

Supongamos que tenemos un endomorfismo $f:{ℝ}^3\to {ℝ}^3$ del espacio vectorial ${ℝ}^3$ que en base canónica tiene por matriz asociada la matriz $A=\left(\begin{array}{ccc}2& 1& 1\\ 0& 4& 1\\ 0& 3& 2\end{array}\right)$.

Supongamos que queremos estudiar, por ejemplo, el endomorfismo $f^{10}$, que tiene por matriz asociada $A^{10}$. Esta matriz es difícil de calcular. Sin embargo, si encontráramos otra base ${ℬ}^{\prime }$ tal que la matriz de $f$ en esa base fuera una matriz diagonal $D=\left(\begin{array}{ccc}{d}_1& 0& 0\\ 0& d_2& 0\\ 0& 0& d_3\end{array}\right)$, con $P$ la matriz de cambio de base de ${ℬ}^{\prime }$ a la base canónica, es decir, $A=PD{P}^{-1}$, tendríamos que $A^{10}=(PD{P}^{-1})(PD{P}^{-1})\stackrel{10\text{ veces}}{\dots }(PD{P}^{-1})(PD{P}^{-1})=PD(P^{-1}P)D(P^{-1}P)\dots (P^{-1}P)D{P}^{-1}=P{D}^{10}{P}^{-1}=P\left(\begin{array}{ccc}{d}_1^{10}& 0& 0\\ 0& d_2^{10}& 0\\ 0& 0& d_3^{10}\end{array}\right)P^{-1}$.

Así, el problema se vería reducido a calcular potencias de números reales, algo mucho más sencillo.

Veamos ahora cómo podemos construir la base ${ℬ}^{\prime }=\{u_1,u_2,u_3\}$ para que la matriz de $f$ en esa base sea diagonal. Por construcción de la matriz de una aplicación lineal en una cierta base, que la matriz de $f$ en base ${ℬ}^{\prime }$ sea $D$ significa que $\{\begin{array}{c}f(u_1)=d_1{u}_1\\ f(u_2)=d_2{u}_2\\ f(u_3)=d_3{u}_3\end{array}$. Tomemos por ejemplo la primera condición y denotemos $ℬ$ la base canónica: $f(u_1)=d_1{u}_1\iff A(u_1{)}_{ℬ}=d_1(u_1{)}_{ℬ}\iff (A-d_{1}I)(u_1{)}_{ℬ}=0$. Como $u_1$ forma parte de una base, no puede ser nulo, de forma que esto quiere decir que el sistema homogéneo $(A-d_{1}I)x=0$ tiene soluciones no triviales (distintas de 0). Por el teorema de Rouché–Frobenius, esto quiere decir que $\mathrm{rg}(A-d_{1}I)<3\iff \det (A-d_{1}I)=0$. Si consideramos el polinomio $Q_f(t)=\det (A-tI)$ esto quiere decir que $d_1$ tiene que ser raíz de $Q_f$. Simétricamente para $u_2$ y $u_3$, obtenemos que $d_2,d_3$ también son raíces de $Q_f$. Es decir, $\{d_i\}$ es el conjunto de raíces del polinomio $Q_f$, al que llamaremos polinomio característico de $f$. Por tanto, si no hay ninguna raíz en el cuerpo donde estamos trabajando, podemos afirmar que el endomorfismo no diagonaliza.

Por tanto, el primer paso es encontrar las raíces de $Q_f$ a las que llamaremos valores propios (o VAPS) de $f$. Una vez encontradas, encontramos los vectores $u_i$ como soluciones no triviales de los sistemas $(A-d_{i}I)x=0$ que sabemos que existen porque hemos impuesto que el rango de la matriz $A-d_{i}I$ baje para los $d_i$ encontrados. A los vectores $u_i$ los llamaremos vectores propios (o VEPS) de $f$. Una vez encontrados, hay que comprobar que formen entre ellos una base. Si la forman, ya tenemos la base ${ℬ}^{\prime }$ que buscábamos. Si no, diremos que $f$ era un endomorfismo no diagonalizable. Este método es válido en general para cualquier dimensión. Es decir, para diagonalizar una matriz los pasos a seguir son:

(1) Encontrar los valores propios: las raíces de $Q_f(t)=\det (A-tI)$. Si no tiene, $f$ no diagonaliza.

(2) Encontrar los vectores propios: para cada valor propio $d_i$, las soluciones de $(A-d_{i}I)x=0$. Tomar, de entre los vectores propios de cada valor propio, tantos linealmente independientes como sea posible.

(3) Comprobar que el conjunto de vectores obtenidos es, efectivamente una base. Si lo es, ya estamos. Si no, $f$ no es diagonalizable.

En el ejemplo anterior haríamos lo siguiente:

$(1)Q_f(t)=\det (A-tI)=\left|\begin{array}{ccc}2-t& 1& 1\\ 0& 4-t& 1\\ 0& 3& 2-t\end{array}\right|=(2-t)\left|\begin{array}{cc}4-t& 1\\ 3& 2-t\end{array}\right|=$

$=(2-t)((4-t)(2-t)-3)=(2-t)(5-6t+t^2)=(2-t)(5-t)(1-t)$. Por lo que los valores propios son $1,2$ y $5$.

$(2)$Encontramos los vectores propios:

$t=1$: $A-I=\left(\begin{array}{ccc}1& 1& 1\\ 0& 3& 1\\ 0& 3& 1\end{array}\right)\sim \left(\begin{array}{ccc}1& 1& 1\\ 0& 3& 1\\ 0& 0& 0\end{array}\right)\Rightarrow \ker (A-I)=\mathrm{span}(-2,1,3)\Rightarrow u_1=(-2,1,3)$

$t=2$: $A-2I=\left(\begin{array}{ccc}0& 1& 1\\ 0& 2& 1\\ 0& 3& 0\end{array}\right)\sim \left(\begin{array}{ccc}0& 1& 0\\ 0& 0& 1\\ 0& 0& 0\end{array}\right)\Rightarrow \ker (A-2I)=\mathrm{span}(1,0,0)\Rightarrow u_2=(1,0,0)$

$t=5$: $A-5I=\left(\begin{array}{ccc}-3& 1& 1\\ 0& -1& 1\\ 0& 3& -3\end{array}\right)\sim \left(\begin{array}{ccc}-3& 0& 2\\ 0& -1& 1\\ 0& 0& 0\end{array}\right)\Rightarrow \ker (A-5I)=\mathrm{span}(2,3,3)\Rightarrow u_3=(2,3,3)$

$(3)$ Nuestro candidato a base es ${ℬ}^{\prime }=\{(-2,1,3),(1,0,0),(2,3,3)\}$, que efectivamente, es una base. Por tanto,

$A=\left(\begin{array}{ccc}-2& 1& 2\\ 1& 0& 3\\ 3& 0& 3\end{array}\right)\left(\begin{array}{ccc}1& 0& 0\\ 0& 2& 0\\ 0& 0& 5\end{array}\right){\left(\begin{array}{ccc}-2& 1& 2\\ 1& 0& 3\\ 3& 0& 3\end{array}\right)}^{-1}$, y podemos calcular $A^k=\left(\begin{array}{ccc}-2& 1& 2\\ 1& 0& 3\\ 3& 0& 3\end{array}\right)\left(\begin{array}{ccc}{1}^k& 0& 0\\ 0& 2^k& 0\\ 0& 0& 5^k\end{array}\right){\left(\begin{array}{ccc}-2& 1& 2\\ 1& 0& 3\\ 3& 0& 3\end{array}\right)}^{-1}$.

Sin embargo, no todos los endomorfismos son diagonalizables, pero podemos caracterizar aquellos que sí que lo son.

Caracterización de endomorfismos diagonalizables

Definimos la multiplicidad algebraica y geométrica de un valor propio $\lambda$. La multiplicidad algebraica $m_a(\lambda )$ es su multiplicidad como raíz del polinomio característico. La multiplicidad geométrica $m_g(\lambda )$ de $\lambda$ es la dimensión del espacio de vectores propios de valor propio $\lambda$ (que denotaremos $E_{\lambda }$). Así, $m_g(\lambda )=\dim (E_{\lambda })$.

Antes de la caracterización, hace falta demostrar dos lemas: Plantilla:DemostraciónPlantilla:Demostración Veamos ahora el teorema que caracteriza a los endomorfismos diagonalizables y nos permitirá asegurar que no podemos diagonalizar $f$ antes de acabar el algoritmo presentado en el primer apartado del artículo. Plantilla:Demostración

Por tanto, en el algoritmo anterior podemos afirmar que $f$ no es diagonalizable sin acabar en dos casos: si el polinomio característico no descompone o si al calcular los subespacios de vectores propios, encontramos alguno cuya dimensión no coincida con la multiplicidad algebraica del VAP correspondiente.

Además, si no hemos parado en ninguno de los dos casos anteriores, el teorema afirma que $f$ es diagonalizable directamente, es decir, no hace falta el tercer paso, donde comprobábamos que el conjunto de VEPs obtenidos eran efectivamente una base. Así, podríamos actualizar el algoritmo como sigue:

(1) Encontrar los valores propios: las raíces de $Q_f(t)=\det (A-tI)$. Si no descompone completamente, $f$ no diagonaliza. Fin.

(2) Encontrar los vectores propios: para cada valor propio $d_i$, $\ker (A-d_{i}I)$. Si la dimensión de este núcleo no es igual a la multiplidad del VAP correspondiente como raíz del polinomio característico, $f$ no diagonaliza. Fin.

(3) Llegado a este punto, $f$ diagonaliza, y la unión de las bases de los espacios de VEPs es la base que buscamos para diagonalizar $f$.

Referencias

• Plantilla:Ref-libro

Bibliografía

Plantilla:En Richard S. Varga, Matrix Iterative Analysis, Springer, 2010