3 Diagonalización

Objetivo: Dada una matriz A𝔐n(𝕂)A\in\mathfrak{{M}}_{n}(\mathbb{K}), saber si existen PP invertible y DD diagonal tal que P1AP=DP^{-1}AP=D (si existen y cuales son).

Equivalentemente, dada una aplicación lineal f:VVf\colon V\to V, buscar una base adecuada de vectores de forma que la matriz de ff respecto de dicha base sea diagonal:

B={v1,,vn} tal que 𝔐B,B(f)=D=(α100αn)B=\left\{v_{1},\ldots,v_{n}\right\}\text{ tal que }\mathfrak{{M}}_{B,B}(f)=D=% \begin{pmatrix}\alpha_{1}&&0\\ &\ddots&\\ 0&&\alpha_{n}\\ \end{pmatrix}

¿Por qué equivalente a lo anterior?

Suponer que ff respecto de la base GG tiene matriz igual a AA.

VG{{V_{G}}}VG{{V_{G}}}D=P1AP{{D=P^{-1}AP}}VB{{V_{B}}}VB{{V_{B}}}A\scriptstyle{A}Cambio de baseP\scriptstyle{\begin{array}[]{c}\text{Cambio de base}\\ P\end{array}}D\scriptstyle{D}Cambio de baseP\scriptstyle{\begin{array}[]{c}\text{Cambio de base}\\ P\end{array}}

Aplicaciones: Entre otras, sirve para calcular facilmente AmA^{m}. En caso de que AA sea diagonalizable (P1AP=D)(P^{-1}AP=D) entonces:

Am=(PDP1)m=PDP1PDP1PDP1=PDmP1A^{m}=(PDP^{-1})^{m}=PD\cancel{P^{-1}}\cdot\cancel{P}D\cancel{P^{-1}}\cdots% \cancel{P}DP^{-1}=P\cdot D^{m}\cdot P^{-1}

donde Dm=(α1m00αnm)D^{m}=\begin{pmatrix}\alpha^{m}_{1}&&0\\ &\ddots&\\ 0&&\alpha^{m}_{n}\\ \end{pmatrix} (fácil).

Y también a sistemas dinámicos, sucesiones recursivas, etc.

Definition 3.1.

Dada una matriz A𝔐n(𝕂)A\in\mathfrak{{M}}_{n}(\mathbb{K}), llamamos polinomio característico de AA al siguiente polinomio de Kn[x]K_{n}[x]:

pA(x)=det(xInA)p_{A}(x)=\det(xI_{n}-A)
Remark 3.2.

El coeficiente de xnx^{n} es 1.

Lemma 3.3.

Si AA y BB son matrices semejantes (ie matrices de la misma aplicación lineal, equivalentemente B=P1APB=P^{-1}AP para PP invertible), entonces

pA(x)=pB(x)p_{A}(x)=p_{B}(x)
Proof 3.4.
pB(x)\displaystyle p_{B}(x)
=det(xInB)=B=P1APdet(xInP1AP)=In=P1Pdet(xP1PP1AP)=\displaystyle{}=det(xI_{n}-B)\overset{B=P^{-1}AP}{=}\det(xI_{n}-P^{-1}AP)% \overset{I_{n}=P^{-1}P}{=}\det(xP^{-1}P-P^{-1}AP)=
=det(P1(xInA)P)=Prop.determinantesdet(P1)det(xInA)det(P)=pA(x)\displaystyle{}=\det(P^{-1}(xI_{n}-A)P)\overset{\begin{subarray}{c}\text{Prop.% }\\ \text{determinantes}\end{subarray}}{=}\cancel{\det(P^{-1})}\cdot\det(xI_{n}-A)% \cdot\cancel{\det(P)}=p_{A}(x)
Definition 3.5.

Se llama polinomio caracteristico de la aplicacion lineal f:VVf\colon V\to V al polinomio caracteristico de cualquiera de las matrices de ff. Notacion: pf(x)p_{f}(x).

Definition 3.6.

Dada una aplicacion lineal f:VVf\colon V\to V decimos que vVv\in V es un vector propio de ff asociado al valor propio λ𝕂\lambda\in\mathbb{K} si v0v\neq 0 y f(v)=λvf(v)=\lambda v.

Remark 3.7.

La base que estamos buscando esta formada por vectores propios y la diagonal de DD por valores propios.

Example 3.8.

Sea la aplicacion lineal

f:2\displaystyle f\colon\mathbb{R}^{2}
2\displaystyle{}\longrightarrow\mathbb{R}^{2}
(x,y)\displaystyle(x,y)
(2x+2y,2x+2y)\displaystyle{}\longmapsto(2x+2y,2x+2y)
  • v1=(1,1)v_{1}=(1,-1) cumple f(v1)=(0,0)=0v1v1f(v_{1})=(0,0)=0\cdot v_{1}\Rightarrow v_{1} es un vector propio asociado al valor propio 0.

  • v2=(1,1)v_{2}=(1,1) cumple f(v2)=(4,4)=4v2v2f(v_{2})=(4,4)=4\cdot v_{2}\Rightarrow v_{2} es un vector propio asociado al valor propio 44.

Además, ff si B={v1,v2}B=\left\{v_{1},v_{2}\right\}, M(f,B,B)=(0004)M(f,B,B)=\left(\begin{array}[]{c|c}0&0\\ 0&4\end{array}\right) (diagonal!)

Definition 3.9 (Subespacio propio).

Dada f:VVf\colon V\to V lineal y λ𝕂\lambda\in\mathbb{K}, llamamos subespacio propio asociado a λ𝕂\lambda\in\mathbb{K} a

Sλ=V(λ)={vVf(v)=λv}.S_{\lambda}=V(\lambda)=\left\{v\in V\mid f(v)=\lambda v\right\}.

Se cumple que SλS_{\lambda} es un subespacio vectorial.

Theorem 3.10.

Sea f:VVf\colon V\to V lineal. Los valores propios de ff son exactamente las raices del polinomio pf(x)=0p_{f}(x)=0.

Proof 3.11.

Lo probaremos escribiendo todo respecto de una base B={v1,,vn}B=\left\{v_{1},\ldots,v_{n}\right\}: A=M(f,B,B)A=M(f,B,B), (v)B=(x1,,xn)(v)_{B}=(x_{1},\cdots,x_{n})

λ𝕂 valor propio\displaystyle\lambda\in\mathbb{K}\text{ valor propio }
def0vV con f(v)=λv\displaystyle{}\overset{def}{\iff}\exists 0\neq v\in V\text{ con }f(v)=\lambda v\iff
0(x1,,xn)𝕂n con A(x1xn)=λ(x1xn)\displaystyle{}\iff\exists 0\neq(x_{1},\ldots,x_{n})\in\mathbb{K}^{n}\text{ % con }A\begin{pmatrix}x_{1}\\ \vdots\\ x_{n}\\ \end{pmatrix}=\lambda\begin{pmatrix}x_{1}\\ \vdots\\ x_{n}\\ \end{pmatrix}\iff
0(x1,,xn)𝕂n con (λInA)(x1xn)=(00)\displaystyle{}\iff\exists 0\neq(x_{1},\ldots,x_{n})\in\mathbb{K}^{n}\text{ % con }(\lambda I_{n}-A)\begin{pmatrix}x_{1}\\ \vdots\\ x_{n}\\ \end{pmatrix}=\begin{pmatrix}0\\ \vdots\\ 0\\ \end{pmatrix}\iff
 el sistema homogeneo (λInA)(x1xn)=(00) es SCI \displaystyle{}\iff\text{ el sistema homogeneo }(\lambda I_{n}-A)\begin{% pmatrix}x_{1}\\ \vdots\\ x_{n}\\ \end{pmatrix}=\begin{pmatrix}0\\ \vdots\\ 0\\ \end{pmatrix}\text{ es SCI }\iff
det(λInA)=0.\displaystyle{}\iff\det(\lambda I_{n}-A)=0.
Example 3.12.

Sea f:22f\colon\mathbb{R}^{2}\to\mathbb{R}^{2}, (x,y)(2x+2y,2x+2y)(x,y)\mapsto(2x+2y,2x+2y).

Tomemos como base la base canonica BC={(1,0),(0,1)}BC=\left\{(1,0),(0,1)\right\}. Entonces, M(f,BC)=(2222)=AM(f,BC)=\begin{pmatrix}2&2\\ 2&2\\ \end{pmatrix}=A.

pA(x)=|x(1001)(2222)|=|x222x2|=x24x=x(x4)p_{A}(x)=\left|x\begin{pmatrix}1&0\\ 0&1\\ \end{pmatrix}-\begin{pmatrix}2&2\\ 2&2\\ \end{pmatrix}\right|=\left|\begin{array}[]{cc }x-2&-2\\ -2&x-2\end{array}\right|=x^{2}-4x=x(x-4)

con lo que pA(x)=0x{0,4}p_{A}(x)=0\iff x\in\left\{0,4\right\} (valores propios).

Lemma 3.13.

Si λ1,,λk\lambda_{1},\ldots,\lambda_{k} son valores propios distintos de ff entonces la suma siguiente es suma directa:

Sλ1SλkS_{\lambda_{1}}\oplus\cdots\oplus S_{\lambda_{k}}

es decir, las bases de vectores de Sλ1,,SλkS_{\lambda_{1}},\ldots,S_{\lambda_{k}} son independientes entre sí.

Proof 3.14.

Sean λ1,,λk𝕂\lambda_{1},\ldots,\lambda_{k}\in\mathbb{K} valores propios distintos.

Veamos que Sλ1SλkS_{\lambda_{1}}\oplus\cdots\oplus S_{\lambda_{k}}. Por ello veremos que si a1v1++akvk=0a_{1}v_{1}+\cdots+a_{k}v_{k}=0 entonces a1==ak=0a_{1}=\cdots=a_{k}=0 (linealmente independientes).

Por inducción en kk:

  • Si k=1k=1: a1v1=0a_{1}\cdot v_{1}=0. Como v10v_{1}\neq 0, a1=0a_{1}=0 y se cumple.

  • Suponer cierto para k1k-1 y veamos el caso kk:

    a1v1++akvk=0f(v1++vk)=a1λ1v1++akλkvk=0a1λkv1++akλkvk=0}a1(λ1λk)v1++ak1(λk1λk)vk1+(λkλk)vk=0{}\begin{rcases}a_{1}v_{1}+\cdots+a_{k}v_{k}=0\Rightarrow f(v_{1}+\cdots+v_{k}% )=a_{1}\lambda_{1}v_{1}+\cdots+a_{k}\lambda_{k}v_{k}=0\\ a_{1}\lambda_{k}v_{1}+\cdots+a_{k}\lambda_{k}v_{k}=0\end{rcases}\Rightarrow\\ {}\Rightarrow a_{1}(\lambda_{1}-\lambda_{k})v_{1}+\cdots+a_{k-1}(\lambda_{k-1}% -\lambda_{k})v_{k-1}+\cancel{(\lambda_{k}-\lambda_{k})v_{k}}=0

    Puesto que por hipótesis de inducción los vectores v1,,vk1v_{1},\ldots,v_{k-1} son linealmente independientes, se tiene que

    a1(λ1λk)==ak1(λk1λk)=0a_{1}(\lambda_{1}-\lambda_{k})=\cdots=a_{k-1}(\lambda_{k-1}-\lambda_{k})=0

    y como los autovalores son distintos entre si,

    a1==ak1=0.a_{1}=\cdots=a_{k-1}=0.

    Sustituyendo en la ecuacion anterior, nos queda akvk=0a_{k}v_{k}=0, con lo que ak=0a_{k}=0 y hemos llegado a que se cumple la tesis de inducción.

Lemma 3.15.

Si λ\lambda es un valor propio de ff, entonces dimSλma(λ)\dim S_{\lambda}\leq m_{a}(\lambda), donde ma(λ)m_{a}(\lambda), llamada multiplicidad algebraica, es el numero de veces que aparece λ\lambda como raiz de pf(x)p_{f}(x).

Proof 3.16.

Sea {v1,,vs}\left\{v_{1},\ldots,v_{s}\right\} una base de SλS_{\lambda}, donde s=dimSλs=\dim S_{\lambda}, y ampliemos a una base de VV: B={v1,,vs,vs+1,vn}B=\left\{v_{1},\ldots,v_{s},v_{s+1},\cdots v_{n}\right\}. Entonces, siendo u1,,udu_{1},\ldots,u_{d} autovectores asociados al autovalor λ\lambda, la matriz asociada a ff respecto de la base BB sera de la forma

B=M(f,B,B)=(λ000λ000λ000)B=M(f,B,B)=\begin{pmatrix}\lambda&0&\cdots&0&*&\cdots&*\\ 0&\lambda&\cdots&0&*&\cdots&*\\ \vdots&\vdots&\ddots&\vdots&\vdots&\ddots&\vdots\\ 0&0&\cdots&\lambda&*&\cdots&*\\ \vdots&\vdots&\ddots&\vdots&\vdots&\ddots&\vdots\\ 0&0&\cdots&0&*&\cdots&*\\ \end{pmatrix}

y el polinomio caracteristico de ff es, desarrollando,

pf(x)=(λx)sq(x)p_{f}(x)=(\lambda-x)^{s}\cdot q(x)

con lo que ma(λ)s=dimSλm_{a}(\lambda)\geq s=\dim S_{\lambda}.

Definition 3.17.

Llamamos multiplicidad geometrica del valor propio λ\lambda, mg(λ)m_{g}(\lambda), a la dimension del subespacio propio VλV_{\lambda}, es decir,

mg(λ)=dimSλm_{g}(\lambda)=\dim S_{\lambda}

Por tanto, el lema 3.15 afirma que mg(λ)ma(λ)m_{g}(\lambda)\leq m_{a}(\lambda).

Definition 3.18.

Decimos que f:VVf\colon V\to V lineal es diagonalizable si existe una base de VV formada por vectores propios. Respecto de dicha base, la matriz de ff es diagonal.

Theorem 3.19 (Criterio de diagonalizacion).

Sea f:VVf\colon V\to V aplicación lineal. Se cumple que

f es diagonalizable {1. la suma de las multiplicidades algebraicas detodos los valores propios es igual a n=dimV2. para cada valor propio λ,ma(λ)=mg(λ)f\text{ es diagonalizable }\iff\begin{cases}\text{1. la suma de las % multiplicidades algebraicas de}\\ \;\;\;\;\text{todos los valores propios es igual a }n=\dim V\\ \text{2. para cada valor propio }\lambda,m_{a}(\lambda)=m_{g}(\lambda)\end{cases}
Proof 3.20.
\Rightarrow

] Sea B={v1,,vn}B=\left\{v_{1},\ldots,v_{n}\right\} la base de vectores propios de ff.

M(f,B)=(λ100λn)pf(x)=(xλ1)(xλn)M(f,B)=\begin{pmatrix}\lambda_{1}&&0\\ &\ddots&\\ 0&&\lambda_{n}\\ \end{pmatrix}\qquad p_{f}(x)=(x-\lambda_{1})\cdots(x-\lambda_{n})

Por un lado, se tiene que el número de λi\lambda_{i}’s es igual a n=dimVn=\dim V (1).

Si λ\lambda es valor propio, por el lema 3.15 ma(λ)mg(λ)m_{a}(\lambda)\geq m_{g}(\lambda). Ademas, si λ\lambda aparece ma(λ)m_{a}(\lambda) veces en pf(x)p_{f}(x) eso quiere decir que aparece al menos ese nº de veces en la diagonal de M(f,B)M(f,B), que provienen de vectores de una base de VV (independientes), por lo que ma(λ)dimSλ=mg(λ)m_{a}(\lambda)\leq\dim S_{\lambda}=m_{g}(\lambda).

\Leftarrow

] Por (1)(1) pf(x)=(xλ1)t1(xλk)tkp_{f}(x)=(x-\lambda_{1})^{t_{1}}\cdots(x-\lambda_{k})^{t_{k}} donde t1=ma(λ1),,tk=ma(λk)t_{1}=m_{a}(\lambda_{1}),\ldots,t_{k}=m_{a}(\lambda_{k}).

Ademas, para cada λi\lambda_{i}, ma(λi)=mg(λi)m_{a}(\lambda_{i})=m_{g}(\lambda_{i}) (2), luego Sλ1SλkS_{\lambda_{1}}\oplus\cdots\oplus S_{\lambda_{k}} es un subespacio de VV y su dimension es dim(Sλ1Sλk)=dimSλ1++dimSλk=mg(λ1)++mg(λk)=2)ma(λ1)++ma(λk)=grado de pf(x)n\dim(S_{\lambda_{1}}\oplus\cdots\oplus S_{\lambda_{k}})=\dim S_{\lambda_{1}}+% \cdots+\dim S_{\lambda_{k}}=m_{g}(\lambda_{1})+\cdots+m_{g}(\lambda_{k})% \overset{2)}{=}m_{a}(\lambda_{1})+\cdots+m_{a}(\lambda_{k})\overset{\text{% grado de }p_{f}(x)}{=}n.

Como consecuencia, Sλ1Sλk=VS_{\lambda_{1}}\oplus\cdots\oplus S_{\lambda_{k}}=V y la base buscada de VV es BSλ1BSλ2BSλkB_{S_{\lambda_{1}}}\cup B_{S_{\lambda_{2}}}\cup\cdots\cup B_{S_{\lambda_{k}}}.

Corollary 3.21.

Si dimV=n\dim V=n y ff tiene nn valores propios distintos entonces ff siempre es diagonalizable.

Proof 3.22.

Si AA tiene nn autovalores distintos, entonces cada autovalor tendrá multiplicidad algebraica ma(λi)=1m_{a}(\lambda_{i})=1. Además, para cada λi\lambda_{i}, mg(λi)1m_{g}(\lambda_{i})\geq 1 ya que λi\lambda_{i} es valor propio si y sólo si Sλi0S_{\lambda_{i}}\neq 0.

Aplicando el lema 3.15,

1=ma(λi)mg(λi)1mg(λi)=ma(λi)=1.1=m_{a}(\lambda_{i})\geq m_{g}(\lambda_{i})\geq 1\Rightarrow m_{g}(\lambda_{i}% )=m_{a}(\lambda_{i})=1.
{metodo}
  1. 1.

    Calcular pf(x)p_{f}(x) y los valores propios con sus multiplicidades algebraicas.

  2. 2.

    Comprobar la primera condición del criterio de diagonalización (que todas las raíces del polinomio característico estén en 𝕂\mathbb{K}).

  3. 3.

    Si ma(λ)=1m_{a}(\lambda)=1 siempre mg(λ)=1=ma(λ)m_{g}(\lambda)=1=m_{a}(\lambda). Comprobar solo los que ma(λ)2m_{a}(\lambda)\geq 2. Para ello hallar SλS_{\lambda}, su dimensión y una base

    Sλ={vf(v)=λv}={v(λIA)v=0}S_{\lambda}=\left\{v\mid f(v)=\lambda v\right\}=\left\{v\mid(\lambda I-A)v=0\right\}

    Si se cumple que ma(λ)=mg(λ)m_{a}(\lambda)=m_{g}(\lambda) en todos los casos, la base buscada es la unión de las bases de SλS_{\lambda}’s.

    DD es la matriz de los vectores propios en la diagonal y PP es la matriz de los vectores propios con sus coordenadas por columnas.

Remark 3.23.

El que una aplicacion lineal sea diagonalizable depende fuertemente del cuerpo sobre el que se trabaje. Por ejemplo, si A=(0110)A=\begin{pmatrix}0&1\\ -1&0\\ \end{pmatrix}:

pa(x)=|x11x|=x2+1p_{a}(x)=\begin{vmatrix}x&-1\\ 1&x\\ \end{vmatrix}=x^{2}+1

  • Si 𝕂=\mathbb{K}=\mathbb{R}, x2+1x^{2}+1 no tiene raices y por tanto no se cumple la primera condicion del criterio de diagonalizacion.

  • Si 𝕂=\mathbb{K}=\mathbb{C}, x2+1=(x+i)(xi){λ=ima(i)=1=mg(i)λ=ima(i)=1=mg(i)x^{2}+1=(x+i)(x-i)\Rightarrow\begin{cases}\lambda=i\;m_{a}(i)=1=m_{g}(i)\\ \lambda=-i\;m_{a}(-i)=1=m_{g}(-i)\end{cases}

    Es decir, AA si es diagonalizable en \mathbb{C}.

Aplicaciones

  • Potencias de matrices: si AA es diagonalizable, P1AP=DAk=(PDP1)k=PDkP1P^{-1}AP=D\Rightarrow A^{k}=(PDP^{-1})^{k}=PD^{k}P^{-1}

  • Sucesiones recursivas:

    {xn=axn1+bxn2x2=αx1=β{(xnxn1)=(ab10)(xn1xn2)(x2x1)=(αβ)\displaystyle\begin{cases}x_{n}=ax_{n-1}+bx_{n-2}\\ x_{2}=\alpha\\ x_{1}=\beta\end{cases}\Rightarrow\begin{cases}\binom{x_{n}}{x_{n-1}}=\begin{% pmatrix}a&b\\ 1&0\\ \end{pmatrix}\binom{x_{n-1}}{x_{n-2}}\\ \binom{x_{2}}{x_{1}}=\binom{\alpha}{\beta}\end{cases}\Rightarrow
    (xnxn1)=A(xn1xn2)=A2(xn2xn3)==An2potencias de A(x2x1)\displaystyle\Rightarrow\binom{x_{n}}{x_{n-1}}=A\binom{x_{n-1}}{x_{n-2}}=A^{2}% \binom{x_{n-2}}{x_{n-3}}=\cdots=\underbrace{A^{n-2}}_{\text{potencias de }A}% \binom{x_{2}}{x_{1}}
  • Sistemas dinamicos …