Articulo de referencia

Radio espectral

En matemáticas , el radio espectral de una matriz cuadrada es el máximo de los valores absolutos de sus autovalores . [ 1 ] De forma más general, el radio espectral de un operad...

En matemáticas , el radio espectral de una matriz cuadrada es el máximo de los valores absolutos de sus autovalores . [ 1 ] De forma más general, el radio espectral de un operador lineal acotado es el supremo de los valores absolutos de los elementos de su espectro . El radio espectral se suele denotar porρ(){\displaystyle \rho (\cdot )}.

Definición

Matrices

Sean λ 1 , ..., λ n los valores propios de una matriz AC n × n . El radio espectral de A se define como

ρ(A)=máximo{|λ1|,,|λnorte|}.{\displaystyle \rho (A)=\max \left\{|\lambda _{1}|,\dotsc ,|\lambda _{n}|\right\}.}

El radio espectral puede pensarse como un ínfimo de todas las normas de una matriz. De hecho, por un lado,ρ(A)A{\displaystyle \rho (A)\leqslant \|A\|}para cada norma de matriz natural{\displaystyle \|\cdot \|}; y por otro lado, la fórmula de Gelfand establece queρ(A)=límitekAk1/k{\displaystyle \rho (A)=\lim _{k\to \infty }\|A^{k}\|^{1/k}}Ambos resultados se muestran a continuación.

Sin embargo, el radio espectral no necesariamente satisfaceAvρ(A)v{\displaystyle \|A\mathbf {v} \|\leqslant \rho (A)\|\mathbf {v} \|}para vectores arbitrariosvdonorte{\displaystyle \mathbf {v} \in \mathbb {C} ^{n}}Para ver por qué, dejemos que...r>1{\displaystyle r>1}ser arbitrario y considerar la matriz

dor=(0r1r0){\displaystyle C_{r}={\begin{pmatrix}0&r^{-1}\\r&0\end{pmatrix}}}.

El polinomio característico dedor{\displaystyle C_{r}}esλ21{\displaystyle \lambda ^{2}-1}, por lo que sus valores propios son{1,1}{\displaystyle \{-1,1\}}y por lo tantoρ(dor)=1{\displaystyle \rho (C_{r})=1}. Sin embargo,dormi1=rmi2{\ Displaystyle C_ {r} \ mathbf {e} _ {1} = r \ mathbf {e} _ {2}}. Como resultado,

dormi1=r>1=ρ(dor)mi1.{\displaystyle \|C_{r}\mathbf {e} _{1}\|=r>1=\rho (C_{r})\|\mathbf {e} _{1}\|.}

Como ilustración de la fórmula de Gelfand, observe quedork1/k1{\displaystyle \|C_{r}^{k}\|^{1/k}\to 1}comok{\displaystyle k\to \infty }, desdedork=I{\displaystyle C_{r}^{k}=I}sik{\displaystyle k}es par ydork=dor{\displaystyle C_{r}^{k}=C_{r}}sik{\displaystyle k}es extraño.

Un caso especial en el queAvρ(A)v{\displaystyle \|A\mathbf {v} \|\leqslant \rho (A)\|\mathbf {v} \|}a pesar devdonorte{\displaystyle \mathbf {v} \in \mathbb {C} ^{n}}es cuandoA{\displaystyle A}es una matriz hermitiana y{\displaystyle \|\cdot \|}es la norma euclidiana . Esto se debe a que cualquier matriz hermitiana es diagonalizable por una matriz unitaria , y las matrices unitarias preservan la longitud del vector. Como resultado,

Av=UDUv=DUvρ(A)Uv=ρ(A)v.{\displaystyle \|A\mathbf {v} \|=\|U^{*}DU\mathbf {v} \|=\|DU\mathbf {v} \|\leqslant \rho (A)\|U\mathbf {v} \|=\rho (A)\|\mathbf {v} \|.}

Operadores lineales acotados

En el contexto de un operador lineal acotado A en un espacio de Banach , los autovalores deben reemplazarse con los elementos del espectro del operador , es decir, los valores.λ{\displaystyle \lambda }para quéAλI{\displaystyle A-\lambda I}no es biyectiva. Denotamos el espectro por

σ(A)={λdo:AλIno es biyectivo}{\displaystyle \sigma (A)=\left\{\lambda \in \mathbb {C} :A-\lambda I\;{\text{no es biyectivo}}\right\}}

El radio espectral se define entonces como el supremo de las magnitudes de los elementos del espectro:

ρ(A)=sorberλσ(A)|λ|{\displaystyle \rho (A)=\sup _{\lambda \in \sigma (A)}|\lambda |}

La fórmula de Gelfand, también conocida como fórmula del radio espectral, también es válida para operadores lineales acotados: dejando{\displaystyle \|\cdot \|}denotemos la norma del operador , tenemos

ρ(A)=límitekAk1k=infknorteAk1k.{\displaystyle \rho (A)=\lim _{k\to \infty }\|A^{k}\|^{\frac {1}{k}}=\inf _{k\in \mathbb {N} ^{*}}\|A^{k}\|^{\frac {1}{k}}.}

Un operador acotado (en un espacio de Hilbert complejo) se denomina operador espectroide si su radio espectral coincide con su radio numérico . Un ejemplo de dicho operador es un operador normal .

Gráficos

El radio espectral de un grafo finito se define como el radio espectral de su matriz de adyacencia .

Esta definición se extiende al caso de grafos infinitos con grados de vértices acotados (es decir, existe algún número real C tal que el grado de cada vértice del grafo es menor que C ). En este caso, para el grafo G definimos:

2(GRAMO)={F:V(GRAMO)R : vV(GRAMO)F(v)2<}.{\displaystyle \ell ^{2}(G)=\left\{f:V(G)\to \mathbf {R} \ :\ \sum \nolimits _{v\in V(G)}\left\|f(v)^{2}\right\|<\infty \right\}.}

Sea γ el operador de adyacencia de G :

{γ:2(GRAMO)2(GRAMO)(γF)(v)=(,v)mi(GRAMO)F(){\displaystyle {\begin{cases}\gamma :\ell ^{2}(G)\to \ell ^{2}(G)\\(\gamma f)(v)=\sum _{(u,v)\in E(G)}f(u)\end{cases}}}

El radio espectral de G se define como el radio espectral del operador lineal acotado γ .

límites superiores

Límites superiores del radio espectral de una matriz

La siguiente proposición proporciona límites superiores sencillos pero útiles para el radio espectral de una matriz.

Proposición. Sea AC n × n con radio espectral ρ ( A ) y una norma matricial submultiplicativa ||⋅|| . Entonces, para cada enterok1{\displaystyle k\geqslant 1}:

ρ(A)Ak1k.{\displaystyle \rho (A)\leq \|A^{k}\|^{\frac {1}{k}}.}

Prueba

Sea ( v , λ ) un par vector propio - valor propio para una matriz A. Por la submultiplicatividad de la norma de la matriz, obtenemos:

|λ|kv=λkv=AkvAkv.{\displaystyle |\lambda |^{k}\|\mathbf {v} \|=\|\lambda ^{k}\mathbf {v} \|=\|A^{k}\mathbf {v} \|\leq \|A^{k}\|\cdot \|\mathbf {v} \|.}

Dado que v ≠ 0 , tenemos

|λ|kAk{\displaystyle |\lambda |^{k}\leq \|A^{k}\|}

y por lo tanto

ρ(A)Ak1k.{\displaystyle \rho (A)\leq \|A^{k}\|^{\frac {1}{k}}.}

Concluyendo la demostración.

Límites superiores para el radio espectral de un grafo

Existen muchos límites superiores para el radio espectral de un grafo en función de su número n de vértices y su número m de aristas. Por ejemplo, si

(k2)(k3)2metronortek(k3)2{\displaystyle {\frac {(k-2)(k-3)}{2}}\leq m-n\leq {\frac {k(k-3)}{2}}}

dónde3knorte{\displaystyle 3\leq k\leq n}es un número entero, entonces [ 2 ]

ρ(GRAMO)2metronortek+52+2metro2norte+94{\displaystyle \rho (G)\leq {\sqrt {2m-n-k+{\frac {5}{2}}+{\sqrt {2m-2n+{\frac {9}{4}}}}}}}

Matrices simétricas

Para matrices de valor realA{\displaystyle A}la desigualdadρ(A)A2{\displaystyle \rho (A)\leq {\|A\|}_{2}}sostiene en particular, donde2{\displaystyle {\|\cdot \|}_{2}}denota la norma espectral . En el caso dondeA{\displaystyle A}es simétrica , esta desigualdad es ajustada:

Teorema. SeaARnorte×norte{\displaystyle A\in \mathbb {R} ^{n\times n}}ser simétrico, es decir,A=AT.{\displaystyle A=A^{T}.}Entonces sostiene queρ(A)=A2.{\displaystyle \rho (A)={\|A\|}_{2}.}

Prueba

Dejar(vi,λi)i=1norte{\displaystyle (v_{i},\lambda _{i})_{i=1}^{n}}sean los pares propios de A. Debido a la simetría de A , todosvi{\displaystyle v_{i}}yλi{\displaystyle \lambda _{i}}son de valor real y los autovectoresvi{\displaystyle v_{i}}son ortonormales . Por definición de la norma espectral, existe unaincógnitaRnorte{\displaystyle x\in \mathbb {R} ^{n}}con incógnita2=1{\displaystyle {\|x\|}_{2}=1}de tal manera queA2=Aincógnita2.{\displaystyle {\|A\|}_{2}={\|Ax\|}_{2}.}Dado que los autovectoresvi{\displaystyle v_{i}}formar una base deRnorte,{\displaystyle \mathbb {R} ^{n},}existen factoresδ1,,δnorteRnorte{\displaystyle \delta _{1},\ldots ,\delta _{n}\in \mathbb {R} ^{n}}de tal manera queincógnita=i=1norteδivi{\displaystyle \textstyle x=\sum _{i=1}^{n}\delta _{i}v_{i}}lo cual implica que

Aincógnita=i=1norteδiAvi=i=1norteδiλivi.{\displaystyle Ax=\sum _{i=1}^{n}\delta _{i}Av_{i}=\sum _{i=1}^{n}\delta _{i}\lambda _{i}v_{i}.}

A partir de la ortonormalidad de los autovectoresvi{\displaystyle v_{i}}resulta que

Aincógnita2=i=1norteδiλivi2=i=1norte|δi||λi|vi2=i=1norte|δi||λi|{\displaystyle {\|Ax\|}_{2}=\|\sum _{i=1}^{n}\delta _{i}\lambda _{i}v_{i}\|_{2}=\sum _{i=1}^{n}{|\delta _{i}|}\cdot {|\lambda _{i}|}\cdot {\|v_{i}\|}_{2}=\sum _{i=1}^{n}{|\delta _{i}|}\cdot {|\lambda _{i}|}}

y

incógnita2=i=1norteδivi2=i=1norte|δi|vi2=i=1norte|δi|.{\displaystyle {\|x\|}_{2}=\|\sum _{i=1}^{n}\delta _{i}v_{i}\|_{2}=\sum _{i=1}^{n}{|\delta _{i}|}\cdot {\|v_{i}\|}_{2}=\sum _{i=1}^{n}{|\delta _{i}|}.}

Desdeincógnita{\displaystyle x}se elige de tal manera que maximiceAincógnita2{\displaystyle {\|Ax\|}_{2}}mientras satisfaceincógnita2=1,{\displaystyle {\|x\|}_{2}=1,}los valores deδi{\displaystyle \delta _{i}}deben ser tales que maximiceni=1norte|δi||λi|{\displaystyle \textstyle \sum _{i=1}^{n}{|\delta _{i}|}\cdot {|\lambda _{i}|}}mientras satisfacei=1norte|δi|=1.{\displaystyle \textstyle \sum _{i=1}^{n}{|\delta _{i}|}=1.}Esto se logra configurandoδk=1{\displaystyle \delta _{k}=1}parak=argramometroaincógnitai=1norte|λi|{\displaystyle k=\mathrm {arg\,max} _{i=1}^{n}{|\lambda _{i}|}}yδi=0{\displaystyle \delta _{i}=0}de lo contrario, obteniendo un valor deAincógnita2=|λk|=ρ(A).{\displaystyle {\|Ax\|}_{2}={|\lambda _{k}|}=\rho (A).}

secuencia de potencia

El radio espectral está estrechamente relacionado con el comportamiento de convergencia de la secuencia de potencias de una matriz; concretamente como lo demuestra el siguiente teorema.

Teorema. Sea AC n × n con radio espectral ρ ( A ) . Entonces ρ ( A ) < 1 si y solo si

límitekAk=0.{\displaystyle \lim _{k\to \infty }A^{k}=0.}

Por otro lado, si ρ ( A ) > 1 ,límitekAk={\displaystyle \lim _{k\to \infty }\|A^{k}\|=\infty }. La afirmación es válida para cualquier elección de norma matricial en C n × n .

Prueba

Supongamos queAk{\displaystyle A^{k}}va a cero comok{\displaystyle k}tiende a infinito. Demostraremos que ρ ( A ) < 1 . Sea ( v , λ ) un par vector propio - valor propio para A . Dado que A k v = λ k v , tenemos

0=(límitekAk)v=límitek(Akv)=límitekλkv=vlímitekλk{\displaystyle {\begin{aligned}0&=\left(\lim _{k\to \infty }A^{k}\right)\mathbf {v} \\&=\lim _{k\to \infty }\left(A^{k}\mathbf {v} \right)\\&=\lim _{k\to \infty }\lambda ^{k}\mathbf {v} \\&=\mathbf {v} \lim _{k\to \infty }\lambda ^{k}\end{aligned}}}

Dado que v ≠ 0 por hipótesis, debemos tener

límitekλk=0,{\displaystyle \lim _{k\to \infty }\lambda ^{k}=0,}

lo cual implica|λ|<1{\displaystyle |\lambda |<1}Dado que esto debe ser cierto para cualquier valor propioλ{\displaystyle \lambda }, podemos concluir que ρ ( A ) < 1 .

Ahora, supongamos que el radio de A es menor que 1. Del teorema de la forma normal de Jordan , sabemos que para todo AC n × n , existen V , JC n × n con V no singular y J diagonal por bloques tales que:

A=VJV1{\displaystyle A=VJV^{-1}}

con

J=[Jmetro1(λ1)0000Jmetro2(λ2)0000Jmetros1(λs1)000Jmetros(λs)]{\displaystyle J={\begin{bmatrix}J_{m_{1}}(\lambda _{1})&0&0&\cdots &0\\0&J_{m_{2}}(\lambda _{2})&0&\cdots &0\\\vdots &\cdots &\ddots &\cdots &\vdots \\0&\cdots &0&J_{m_{s-1}}(\lambda _{s-1})&0\\0&\cdots &\cdots &0&J_{m_{s}}(\lambda _{s})\end{bmatrix}}}

dónde

Jmetroi(λi)=[λi1000λi1000λi1000λi]dometroi×metroi,1is.{\displaystyle J_{m_{i}}(\lambda _{i})={\begin{bmatrix}\lambda _{i}&1&0&\cdots &0\\0&\lambda _{i}&1&\cdots &0\\\vdots &\vdots &\ddots &\ddots &\vdots \\0&0&\cdots &\lambda _{i}&1\\0&0&\cdots &0&\lambda _{i}\end{bmatrix}}\in \mathbf {C} ^{m_{i}\times m_{i}},1\leq i\leq s.}

Es fácil ver que

Ak=VJkV1{\displaystyle A^{k}=VJ^{k}V^{-1}}

y, puesto que J es diagonal por bloques,

Jk=[Jmetro1k(λ1)0000Jmetro2k(λ2)0000Jmetros1k(λs1)000Jmetrosk(λs)]{\displaystyle J^{k}={\begin{bmatrix}J_{m_{1}}^{k}(\lambda _{1})&0&0&\cdots &0\\0&J_{m_{2}}^{k}(\lambda _{2})&0&\cdots &0\\\vdots &\cdots &\ddots &\cdots &\vdots \\0&\cdots &0&J_{m_{s-1}}^{k}(\lambda _{s-1})&0\\0&\cdots &\cdots &0&J_{m_{s}}^{k}(\lambda _{s})\end{bmatrix}}}

Ahora, un resultado estándar sobre la potencia k de unmetroi×metroi{\displaystyle m_{i}\times m_{i}}El bloque de Jordania afirma que, parakmetroi1{\displaystyle k\geq m_{i}-1}:

Jmetroik(λi)=[λik(k1)λik1(k2)λik2(kmetroi1)λikmetroi+10λik(k1)λik1(kmetroi2)λikmetroi+200λik(k1)λik1000λik]{\displaystyle J_{m_{i}}^{k}(\lambda _{i})={\begin{bmatrix}\lambda _{i}^{k}&{k \choose 1}\lambda _{i}^{k-1}&{k \choose 2}\lambda _{i}^{k-2}&\cdots &{k \choose m_{i}-1}\lambda _{i}^{k-m_{i}+1}\\0&\lambda _{i}^{k}&{k \choose 1}\lambda _{i}^{k-1}&\cdots &{k \choose m_{i}-2}\lambda _{i}^{k-m_{i}+2}\\\vdots &\vdots &\ddots &\ddots &\vdots \\0&0&\cdots &\lambda _{i}^{k}&{k \choose 1}\lambda _{i}^{k-1}\\0&0&\cdots &0&\lambda _{i}^{k}\end{bmatrix}}}

Por lo tanto, siρ(A)<1{\displaystyle \rho (A)<1}entonces por todo lo que yo|λi|<1{\displaystyle |\lambda _{i}|<1}Por lo tanto, para todo i tenemos:

límitekJmetroik=0{\displaystyle \lim _{k\to \infty }J_{m_{i}}^{k}=0}

lo cual implica

límitekJk=0.{\displaystyle \lim _{k\to \infty }J^{k}=0.}

Por lo tanto,

límitekAk=límitekVJkV1=V(límitekJk)V1=0{\displaystyle \lim _{k\to \infty }A^{k}=\lim _{k\to \infty }VJ^{k}V^{-1}=V\left(\lim _{k\to \infty }J^{k}\right)V^{-1}=0}

Por otro lado, siρ(A)>1{\displaystyle \rho (A)>1}, hay al menos un elemento en J que no permanece acotado a medida que k aumenta, lo que demuestra la segunda parte de la afirmación.

La fórmula de Gelfand

La fórmula de Gelfand, que recibe su nombre de Israel Gelfand , proporciona el radio espectral como límite de normas matriciales.

Teorema

Para cualquier norma matricial ||⋅||, tenemos [ 3 ]

ρ(A)=límitekAk1k{\displaystyle \rho (A)=\lim _{k\to \infty }\left\|A^{k}\right\|^{\frac {1}{k}}}.

Además, en el caso de una norma matricial consistentelímitekAk1k{\displaystyle \lim _{k\to \infty }\left\|A^{k}\right\|^{\frac {1}{k}}}aprochesρ(A){\displaystyle \rho (A)}desde arriba (de hecho, en ese caso)ρ(A)Ak1k{\displaystyle \rho (A)\leq \left\|A^{k}\right\|^{\frac {1}{k}}}a pesar dek{\displaystyle k}).

Prueba

Para cualquier ε > 0 , definamos las dos matrices siguientes:

A±=1ρ(A)±εA.{\displaystyle A_{\pm }={\frac {1}{\rho (A)\pm \varepsilon }}A.}

De este modo,

ρ(A±)=ρ(A)ρ(A)±ε,ρ(A+)<1<ρ(A).{\displaystyle \rho \left(A_{\pm }\right)={\frac {\rho (A)}{\rho (A)\pm \varepsilon }},\qquad \rho (A_{+})<1<\rho (A_{-}).}

Comenzamos aplicando el teorema anterior sobre límites de secuencias de potencias a A + :

límitekA+k=0.{\displaystyle \lim _{k\to \infty }A_{+}^{k}=0.}

Esto demuestra la existencia de N +N tal que, para todo kN + ,

A+k<1.{\displaystyle \left\|A_{+}^{k}\right\|<1.}

Por lo tanto,

Ak1k<ρ(A)+ε.{\displaystyle \left\|A^{k}\right\|^{\frac {1}{k}}<\rho (A)+\varepsilon .}

De manera similar, el teorema sobre secuencias de potencias implica queAk{\displaystyle \|A_{-}^{k}\|}no está acotada y existe N N tal que, para todo kN ,

Ak>1.{\displaystyle \left\|A_{-}^{k}\right\|>1.}

Por lo tanto,

Ak1k>ρ(A)ε.{\displaystyle \left\|A^{k}\right\|^{\frac {1}{k}}>\rho (A)-\varepsilon .}

Sea N = max{ N + , N }. Entonces,

ε>0nortenorteknorteρ(A)ε<Ak1k<ρ(A)+ε,{\displaystyle \forall \varepsilon >0\quad \exists N\in \mathbf {N} \quad \forall k\geq N\quad \rho (A)-\varepsilon <\left\|A^{k}\right\|^{\frac {1}{k}}<\rho (A)+\varepsilon ,}

eso es,

límitekAk1k=ρ(A).{\displaystyle \lim _{k\to \infty }\left\|A^{k}\right\|^{\frac {1}{k}}=\rho (A).}

Con esto concluye la demostración.

Corolario

La fórmula de Gelfand proporciona una cota para el radio espectral de un producto de matrices conmutativas: siA1,,Anorte{\displaystyle A_{1},\ldots ,A_{n}}son matrices que conmutan entre sí, entonces

ρ(A1Anorte)ρ(A1)ρ(Anorte).{\displaystyle \rho (A_{1}\cdots A_{n})\leq \rho (A_{1})\cdots \rho (A_{n}).}

Ejemplo numérico

La convergencia de las 3 normas matriciales al radio espectral.

Consideremos la matriz

A=[912284118]{\displaystyle A={\begin{bmatrix}9&-1&2\\-2&8&4\\1&1&8\end{bmatrix}}}

cuyos autovalores son 5, 10, 10 ; por definición, ρ ( A ) = 10. En la siguiente tabla, los valores deAk1k{\displaystyle \|A^{k}\|^{\frac {1}{k}}}para las cuatro normas más utilizadas se enumeran frente a varios valores crecientes de k (tenga en cuenta que, debido a la forma particular de esta matriz,.1=.{\displaystyle \|.\|_{1}=\|.\|_{\infty }}):

Notas y referencias

  1. Gradshteĭn, IS (1980). Tabla de integrales, series y productos . IM Ryzhik, Alan Jeffrey (Ed. corregida y ampliada  ). Nueva York: Academic Press. ISBN 0-12-294760-6OCLC 5892996 
  2. Guo, Ji-Ming; Wang, Zhi-Wen; Li, Xin (2019). "Límites superiores precisos del radio espectral de un grafo" . Matemáticas Discretas . 342 (9): 2559– 2563. doi : 10.1016/j.disc.2019.05.017 . S2CID 198169497 . 
  3. La fórmula es válida para cualquier álgebra de Banach ; véase el lema IX.1.8 en Dunford y Schwartz 1963 y Lax 2002 , págs. 195-197. 

Bibliografía

  • Dunford, Nelson; Schwartz, Jacob (1963), Operadores lineales II. Teoría espectral: operadores autoadjuntos en el espacio de Hilbert , Interscience Publishers, Inc.
  • Lax, Peter D. (2002), Análisis funcional , Wiley-Interscience, ISBN 0-471-55604-1

Véase también