4 votos

Distribución de valores propios límite de $(I-A)^T(I-A)$

Dejemos que $A\in\mathbb R^{n\times n}$ sea una matriz gaussiana aleatoria con entradas i.i.d de $\mathcal N (0, \frac{a}{\sqrt{n}})$ . Por Marchenko-Pastur conocemos la distribución límite del valor propio de $A^TA$ . Ahora quiero entender el comportamiento de los valores propios de $(I-A)^{T}(I-A)$ . Según algunas simulaciones, si $a$ es grande, entonces la distribución converge a la ley del cuarto de círculo (lo cual es intuitivo ya que el efecto de $I$ se encoge), y un pequeño $a$ conduce a un semicírculo deformado. ¿Existe una representación analítica de esta distribución? Gracias.

a=10

a=0.01

0 votos

El semicírculo se explica por la Ley del semicírculo de Wigner .

0voto

Matthias Benkard Puntos 11264

Sé que la distribución de valores propios límite satisface un principio variacional (véase, por ejemplo, el trabajo conjunto con A. Kuijlaars Grandes desviaciones para una matriz Wishart no centrada ) a partir de la cual se puede derivar una ecuación algebraica de grado tres para la transformada de Sieltjes, y eventualmente recuperar una fórmula analítica para la densidad límite utilizando la fórmula de inversión de Cauchy-Plemelj.

De forma más general, la distribución límite puede expresarse como la convolución libre rectangular (introducida por Benaych-Georges en Matrices aleatorias rectangulares, convolución relacionada ) de la distribución de Marchenko-Pastur y una masa de Dirac; a continuación se puede recuperar esta ecuación algebraica utilizando la transformada R rectangular asociada.

Espero que esto te ayude y que puedas resolver los detalles por ti mismo. Por lo demás, yo partiría de la prueba original de Marchenko-Pastur y la adaptaría para calcular la asintótica de la transformada de Stieljes de tu modelo matricial, lo que eventualmente llevaría también, tras tediosos cálculos, a la misma ecuación algebraica de grado tres.

0 votos

Estas referencias son muy útiles. Para volver a comprobar lo que he entendido: 1) la convolución libre de Marchenko-Pastur con una masa de Dirac describe la distribución de valores propios de $A^TA+I$ ¿cierto? Así que la distribución que quiero es ligeramente diferente (que creo que es $(\sqrt{\mu_{MP}}\oplus \sqrt{\delta_1})^2$ ); y 2) el procedimiento estándar después de aplicar la transformada R sigue siendo encontrar la transformada de Stieljes y aplicar la relación de Plemelj para encontrar el pdf. ¿Es esto correcto o estoy pasando algo por alto? Gracias.

0 votos

1) Existen diferentes circunvalaciones libres. La que mencionas describe efectivamente la distribución límite para A+B para A,B asintóticamente libres (suele llamarse convolución libre aditiva); también existe la convolución libre mutiplicativa para A^{1/2} B A^{1/2} y la rectangular para (A+B)^T(A+B). No creo que se recupere la ley que mencionas. 2) Efectivamente, esta es la aproximación que sugiero.

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X