19 votos

Valor propio más pequeño de una matriz aleatoria complicada

Mientras experimentaba con funciones definidas positivas, llegué a la siguiente conclusión:

Sea $n$ sea un número entero positivo y $x_1,\ldots,x_n$ ser muestreada a partir de una gaussiana de media cero y varianza unitaria. Consideremos la matriz (positiva-definida) $$M_{ij}=\frac{1}{1+|x_i-x_j|}.$$ Ahora deseo saber:

¿Cómo puedo obtener una estimación para el valor propio más pequeño $\lambda_n$ de $M$ ?

Experimentos preliminares (véase el gráfico; eje x: $n$ eje y: $\lambda_n$ ) sugieren que $\lambda_n \approx 1/n^2$ ¿pero cómo puedo probarlo o tal vez un resultado más exacto?

enter image description here

0 votos

¿Quieres decir $M_{ij}=\frac{1}{1+|x_i-x_j|}$ ? Porque $M$ como lo has definido no parece ser una matriz.

1 votos

Aquí tienes un enlace: mathoverflow.net/preguntas/24287/

0 votos

@Gilead: Gracias por detectar la errata; además, no busco un método numérico, sino un límite o aproximación "de forma cerrada".

19voto

steevc Puntos 211

Creo que puedo obtener un límite superior de $O(1/n^2)$ exhibiendo un vector $v$ de magnitud comparable a $1$ que se asigna a un vector de magnitud $O(1/n^2)$ . La idea básica es explotar la paradoja del cumpleaños para encontrar (con alta probabilidad) dos índices $i \neq i'$ tal que $x_i-x_{i'} = O(1/n^2)$ . También debería ser posible encontrar otro índice adicional $i''$ tal que $x_{i''} = x_i + O(1/n)$ .

Ahora mira el $i^{th}$ y $(i')^{th}$ filas, que tienen componentes $1/(1+|x_i-x_j|)$ y $1/(1+|x_{i'}-x_j|)$ . Estas filas se diferencian por $O(n^{-2})$ en cada coeficiente. Esto ya da un límite superior de $O(n^{-3/2})$ para el valor propio más pequeño, pero se puede hacer mejor utilizando la expansión de Taylor para observar que la diferencia entre los dos componentes es $(x_i-x_{i'}) \hbox{sgn}( x_i-x_j ) / (1 + |x_i-x_j|)^2 + O(n^{-4})$ excepto cuando $x_j$ está muy cerca de $x_i$ en cuyo punto sólo tenemos el límite bruto de $O(n^{-2})$ . Del mismo modo, la diferencia entre el $i''$ y $i$ filas es algo así como $(x_i-x_{i''}) \hbox{sgn}( x_i-x_j ) / (1 + |x_i-x_j|)^2 + O(n^{-4})$ excepto cuando $x_j$ está demasiado cerca de $x_i$ . Por tanto, podemos utilizar un múltiplo de la segunda diferencia para anular en su mayor parte la primera diferencia, y acabar con una combinación lineal de tres filas en la que la mayoría de las entradas tienen tamaño $O(n^{-4})$ y sólo unos $O(1)$ entradas tienen tamaño $O(n^{-2})$ . Esto parece dar un límite superior de $O(n^{-2})$ para el menor valor propio (o menor valor singular), aunque no he comprobado todos los detalles.

Obtener un límite inferior equivalente es más complicado. Puede que haya que pasar a una representación de Fourier de la matriz, ya que así se captaría más fácilmente la definición positiva de la matriz (como sugiere el teorema de Bochner).

13 votos

Es incluso más fácil que eso: una vez que hayas encontrado un par de índices que den orden de distancia $1/n^2$ definen un $2 \times 2$ con un valor propio de orden $1/n^2$ y luego por entrelazamiento de valores propios para matrices hermitianas (o simétricas reales) toda la matriz tiene un valor propio no mayor que ese.

0 votos

@Terry: gracias por el buen argumento para un límite superior (y @Tracy, buen punto el de invocar el teorema del entrelazamiento de Cauchy para argumentar sobre el límite superior de los valores propios). En cuanto a la conexión de Fourier, efectivamente esta matriz está generada por la función posdef $f(x)=1/(1+|x|)$ . De hecho, tal vez se pueda demostrar aún más, porque en realidad dicha matriz es infinitamente divisible.

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