stringtranslate.com

Distribución de Tracy-Widom

Densidades de distribuciones de Tracy-Widom para β  = 1, 2, 4

La distribución de Tracy-Widom es una distribución de probabilidad de la teoría de matrices aleatorias introducida por Craig Tracy y Harold Widom  (1993, 1994). Es la distribución del valor propio normalizado más grande de una matriz hermítica aleatoria . La distribución se define como un determinante de Fredholm .

En términos prácticos, Tracy-Widom es la función de cruce entre las dos fases de componentes débilmente acoplados versus fuertemente acoplados en un sistema. [1] También aparece en la distribución de la longitud de la subsecuencia creciente más larga de permutaciones aleatorias , [2] como estadísticas a gran escala en la ecuación de Kardar-Parisi-Zhang , [3] en fluctuaciones actuales del proceso de exclusión simple asimétrico (ASEP) con condición inicial escalonada, [4] y en modelos matemáticos simplificados del comportamiento del problema de la subsecuencia común más larga en entradas aleatorias. [5] Consulte Takeuchi y Sano (2010) y Takeuchi et al. (2011) para pruebas experimentales (y verificación) de que las fluctuaciones de la interfaz de una gota en crecimiento (o sustrato) se describen mediante la distribución TW (o ) como lo predijeron Prähofer y Spohn (2000).

La distribución es de particular interés en las estadísticas multivariadas . [6] Para una discusión de la universalidad de , , véase Deift (2007). Para una aplicación de para inferir la estructura de la población a partir de datos genéticos, véase Patterson, Price y Reich (2006). En 2017 se demostró que la distribución F no es infinitamente divisible. [7]

Definición como ley de grandes números

Sea la función de distribución acumulativa de la distribución de Tracy-Widom con . Puede definirse como una ley de los grandes números, similar al teorema del límite central .

Normalmente hay tres distribuciones de Tracy-Widom, , con . Corresponden a los tres conjuntos gaussianos : ortogonal ( ), unitario ( ) y simpléctico ( ).

En general, considere un conjunto gaussiano con valor beta , con sus entradas diagonales que tienen varianza 1, y entradas fuera de la diagonal que tienen varianza , y sea probabilidad de que una matriz muestreada del conjunto tenga valor propio máximo , luego defina [8] donde denota el valor propio más grande de la matriz aleatoria. El desplazamiento por centra la distribución, ya que en el límite, la distribución de valores propios converge a la distribución semicircular con radio . La multiplicación por se utiliza porque la desviación estándar de la distribución escala como (derivada por primera vez en [9] ).

Por ejemplo: [10]

donde la matriz se muestrea del conjunto unitario gaussiano con varianza fuera de la diagonal .

La definición de las distribuciones de Tracy-Widom puede extenderse a todas (Diapositiva 56 en Edelman (2003), Ramírez, Rider y Virág (2006)).

Naturalmente, se puede preguntar por la distribución límite de los segundos y terceros valores propios más grandes, etc. Son conocidos. [11] [8]

Formas funcionales

Determinante de Fredholm

puede darse como el determinante de Fredholm

del núcleo ("núcleo de Airy") sobre funciones integrables cuadradas en la semirrecta , dadas en términos de funciones de Airy Ai por

Trascendentes de Painlevé

También se puede dar como una integral

en términos de una solución [nota 1] de una ecuación de Painlevé de tipo II

con condición de contorno Esta función es una trascendente de Painlevé .

Otras distribuciones también se pueden expresar en términos de lo mismo : [10]

Ecuaciones funcionales

Definir entonces [8]

Ocurrencias

Además de en la teoría de matrices aleatorias, las distribuciones de Tracy-Widom aparecen en muchos otros problemas de probabilidad. [12]

Sea la longitud de la subsecuencia creciente más larga en una permutación aleatoria muestreada uniformemente de , el grupo de permutación en n elementos. Entonces la función de distribución acumulativa de converge a . [13]

Asintóticos

Función de densidad de probabilidad

Sea la función de densidad de probabilidad para la distribución, entonces [12] En particular, vemos que está severamente sesgada hacia la derecha: es mucho más probable que sea mucho mayor que que mucho menor. Esto se puede intuir al ver que la distribución límite es la ley del semicírculo, por lo que hay "repulsión" desde la mayor parte de la distribución, lo que obliga a que no sea mucho menor que .

En el límite, una expresión más precisa es (ecuación 49 [12] ) para algún número positivo que depende de .

Función de distribución acumulativa

En el límite, [14] y en el límite, donde es la función zeta de Riemann , y .

Esto permite derivar el comportamiento de . Por ejemplo,

Painlevé trascendente

La trascendente de Painlevé tiene expansión asintótica en (ecuación 4.1 de [15] ). Esto es necesario para los cálculos numéricos, ya que la solución es inestable: cualquier desviación de ella tiende a dejarla caer en la rama. [16]

Numéricos

Las técnicas numéricas para obtener soluciones numéricas a las ecuaciones de Painlevé de los tipos II y V, y evaluar numéricamente distribuciones de valores propios de matrices aleatorias en los conjuntos beta fueron presentadas por primera vez por Edelman & Persson (2005) usando MATLAB . Estas técnicas de aproximación fueron justificadas analíticamente en Bejan (2005) y utilizadas para proporcionar una evaluación numérica de las distribuciones de Painlevé II y Tracy–Widom (para ) en S-PLUS . Estas distribuciones han sido tabuladas en Bejan (2005) a cuatro dígitos significativos para valores del argumento en incrementos de 0,01; una tabla estadística para valores p también fue proporcionada en este trabajo. Bornemann (2010) proporcionó algoritmos precisos y rápidos para la evaluación numérica de y las funciones de densidad para . Estos algoritmos pueden usarse para calcular numéricamente la media , la varianza , la asimetría y el exceso de curtosis de las distribuciones . [17]

Las funciones para trabajar con las leyes de Tracy-Widom también se presentan en el paquete R 'RMTstat' de Johnstone et al. (2009) y en el paquete MATLAB 'RMLab' de Dieng (2006).

Para una aproximación simple basada en una distribución gamma desplazada, consulte Chiani (2014).

Shen y Serkh (2022) desarrollaron un algoritmo espectral para la descomposición propia del operador integral , que se puede utilizar para evaluar rápidamente las distribuciones de Tracy-Widom o, de forma más general, las distribuciones del nivel más grande en el límite de escala del borde suave de los conjuntos gaussianos, con precisión de máquina.

Tracy-Widom y la universalidad de KPZ

La distribución de Tracy-Widom aparece como una distribución límite en la clase de universalidad de la ecuación KPZ . Por ejemplo, aparece bajo el escalamiento de la ecuación KPZ unidimensional con tiempo fijo. [18]

Véase también

Notas al pie

  1. ^ Una misteriosa ley estadística podría finalmente tener una explicación, wired.com 2014-10-27
  2. ^ Baik, Deift y Johansson (1999).
  3. ^ Sasamoto y Spohn (2010)
  4. ^ Johansson (2000); Tracy y Widom (2009)).
  5. ^ Majumdar y Nechaev (2005).
  6. ^ Johnstone (2007, 2008, 2009).
  7. ^ Domínguez-Molina (2017).
  8. ^ abc Tracy, Craig A.; Widom, Harold (2009b). "Las distribuciones de la teoría de matrices aleatorias y sus aplicaciones". En Sidoravičius, Vladas (ed.). Nuevas tendencias en física matemática . Dordrecht: Springer Netherlands. págs. 753–765. doi :10.1007/978-90-481-2810-5_48. ISBN 978-90-481-2810-5.
  9. ^ Forrester, PJ (9 de agosto de 1993). "El borde espectral de conjuntos de matrices aleatorias". Física nuclear B . 402 (3): 709–728. Código Bibliográfico :1993NuPhB.402..709F. doi :10.1016/0550-3213(93)90126-A. ISSN  0550-3213.
  10. ^ por Tracy y Widom (1996).
  11. ^ Dieng, Momar (2005). "Funciones de distribución para valores propios de aristas en conjuntos ortogonales y simplécticos: representaciones de Painlevé". International Mathematics Research Notices . 2005 (37): 2263–2287. doi : 10.1155/IMRN.2005.2263 . ISSN  1687-0247.
  12. ^ abc Majumdar, Satya N; Schehr, Grégory (31 de enero de 2014). "Valor propio superior de una matriz aleatoria: grandes desviaciones y transición de fase de tercer orden". Journal of Statistical Mechanics: Theory and Experiment . 2014 (1): 01012. arXiv : 1311.0580 . Bibcode :2014JSMTE..01..012M. doi :10.1088/1742-5468/2014/01/p01012. ISSN  1742-5468. S2CID  119122520.
  13. ^ Baik, Deift y Johansson 1999
  14. ^ Baik, Jinho; Buckingham, Robert; DiFranco, Jeffery (26 de febrero de 2008). "Asintótica de distribuciones de Tracy-Widom y la integral total de una función de Painlevé II". Communications in Mathematical Physics . 280 (2): 463–497. arXiv : 0704.3636 . Código Bibliográfico :2008CMaPh.280..463B. doi :10.1007/s00220-008-0433-5. ISSN  0010-3616. S2CID  16324715.
  15. ^ Tracy, Craig A.; Widom, Harold (mayo de 1993). "Distribuciones de espaciado de niveles y el núcleo de Airy". Physics Letters B . 305 (1–2): 115–118. arXiv : hep-th/9210074 . Código Bibliográfico :1993PhLB..305..115T. doi :10.1016/0370-2693(93)91114-3. ISSN  0370-2693. S2CID  13912236.
  16. ^ Bender, Carl M.; Orszag, Steven A. (29 de octubre de 1999). Métodos matemáticos avanzados para científicos e ingenieros I: métodos asintóticos y teoría de perturbaciones. Springer Science & Business Media. págs. 163–165. ISBN 978-0-387-98931-0.
  17. ^ Su, Zhong-gen; Lei, Yu-huan; Shen, Tian (1 de marzo de 2021). "Distribución Tracy-Widom, proceso Airy2 y sus propiedades de ruta de muestra". Matemáticas Aplicadas: una revista de universidades chinas . 36 (1): 128–158. doi : 10.1007/s11766-021-4251-2 . ISSN  1993-0445. S2CID  237903590.
  18. ^ Amir, Gideon; Corwin, Ivan; Quastel, Jeremy (2010). "Distribución de probabilidad de la energía libre del polímero aleatorio dirigido continuo en 1 + 1 dimensiones". Communications on Pure and Applied Mathematics . 64 (4). Wiley: 466–537. arXiv : 1003.0443 . doi :10.1002/cpa.20347.
  1. ^ llamada "solución Hastings–McLeod". Publicada por Hastings, SP, McLeod, JB: Un problema de valor de contorno asociado con el segundo trascendente de Painlevé y la ecuación de Korteweg-de Vries. Arch. Ration. Mech. Anal. 73 , 31–51 (1980)

Referencias

Lectura adicional

Enlaces externos