www.wikidata.de-de.nina.az
Eine Zufallsmatrix bezeichnet in der Stochastik eine matrixwertige Zufallsvariable englisch Random Matrix Ihre Verteilung nennt man zur Abgrenzung von den multivariaten Verteilungen eine matrixvariate Wahrscheinlichkeitsverteilung Zufallsmatrizen spielen eine wichtige Rolle in der statistischen sowie mathematischen Physik insbesondere in der statistischen Mechanik Aus historischer Sicht hat sich die Theorie aus dem Versuch entwickelt Systeme mit vielen stochastischen aber miteinander agierenden Teilchen zu beschreiben Viele der Grundlagen der Theorie stammen deshalb von mathematischen Physikern und viele Modelle haben eine physikalische Interpretation Die Theorie der Zufallsmatrizen ist auch in der multivariaten Statistik relevant wo man sie zur Analyse von Kovarianzmatrizen benotigt Insbesondere im Zusammenhang mit hoch dimensionalen Daten und spektralstatistischen Verfahren wie der Hauptkomponentenanalyse PCA Zufallsmatrizen sind zu unterscheiden von der stochastischen Matrix Inhaltsverzeichnis 1 Haarsches Mass und Weylsche Integralformel 2 Definition 2 1 Zentrale Begriffe 2 1 1 Partitionsfunktion 2 1 2 Wignersche Matrix 2 1 3 Invariante Ensembles 2 1 4 Gausssche Ensembles 2 1 4 1 Herleitung des GUE durch Ornstein Uhlenbeck Prozesse 2 1 5 Zirkulare Ensembles 2 1 6 b Ensembles und Dysons Threefolded Way 2 2 Theorie der Zufallsmatrizen 3 Spektraltheorie der Zufallsmatrizen 3 1 Globale Situation 3 1 1 Das Variationsproblem der Verteilung der Eigenwerte 3 1 1 1 Beispiel Wignersche Halbkreis 3 2 Lokale Situation 4 Geschichte 5 Anwendungen 5 1 Statistik 5 2 Physik 5 3 Weitere Anwendungen 6 Literatur 7 Weblinks 8 Einzelnachweise und AnmerkungenHaarsches Mass und Weylsche Integralformel BearbeitenAuf jeder Lie Gruppe G displaystyle G nbsp existiert ein eindeutiges links invariantes Mass m L displaystyle mu L nbsp d h fur jedes g G displaystyle g in G nbsp und jede Borel messbare Menge S G displaystyle S subseteq G nbsp gilt m L S g m L g displaystyle mu L Sg mu L g nbsp Dieses Mass nennt man linkes Haarsches Mass und es ist eindeutig bis auf Multiplikation mit einer Konstanten Betrachtet man nun eine kompakte Lie Gruppe G displaystyle G nbsp so existiert ein eindeutiges linkes Haarsches Mass m H displaystyle mu H nbsp welches zu gleich auch rechts invariant und normalisiert ist genannt das Haarsche Wahrscheinlichkeitsmass auf G displaystyle G nbsp Das heisst fur jedes g G displaystyle g in G nbsp und jede Borel messbare Menge B G displaystyle B subseteq G nbsp gilt m H B g m H g B m H g displaystyle mu H Bg mu H gB mu H g nbsp Fur kompakte Lie Gruppen lasst sich mit Hilfe der Integralformel von Weyl eine Formel fur die Wahrscheinlichkeitsdichte bezuglich der Eigenwerte finden Als Beispiel sei G U n displaystyle G mathbb U n nbsp die unitare Gruppe die Eigenwerte sind von der Form e i ϕ 1 e i ϕ n displaystyle e i phi 1 dots e i phi n nbsp mit ϕ 1 ϕ n R displaystyle phi 1 ldots phi n in mathbb R nbsp Weiter sei f displaystyle f nbsp eine Klassenfunktion und T displaystyle T nbsp der maximale Torus aller Diagonalmatrizen t diag e i ϕ 1 e i ϕ n displaystyle t operatorname diag e i phi 1 ldots e i phi n nbsp von U n displaystyle mathbb U n nbsp Ad displaystyle operatorname Ad nbsp bezeichnet die adjungierte Darstellung und R displaystyle R nbsp bezeichne das Wurzelsystem dann gilt 1 det Id Ad G T t 1 a R 1 e a j lt k e i ϕ k e i ϕ j 2 displaystyle det operatorname Id operatorname Ad G T t 1 prod limits alpha in R 1 e alpha prod limits j lt k e i phi k e i phi j 2 nbsp und somit kriegt man mit Hilfe von Weyl s Integralformel ein Integral uber den maximalen Torus T displaystyle T nbsp G f d g 1 n T f t j lt k e i ϕ k e i ϕ j 2 d t displaystyle int G f mathrm d g frac 1 n int T f t prod limits j lt k e i phi k e i phi j 2 mathrm d t nbsp Definition BearbeitenEine formale mathematische Definition lautet 2 Sei M n displaystyle mathcal M n nbsp der Raum der n n displaystyle n times n nbsp Matrizen uber dem Korper K displaystyle K nbsp mit einer s Algebra A displaystyle mathcal A nbsp und W F P displaystyle Omega mathcal F mathbb P nbsp ein Wahrscheinlichkeitsraum Eine F A displaystyle mathcal F mathcal A nbsp messbare Funktion A W M n displaystyle A colon Omega to mathcal M n nbsp heisst Zufallsmatrix Als s displaystyle sigma nbsp Algebra kann die borelsche s Algebra des euklidischen Umgebungsraumes der Mannigfaltigkeit M n displaystyle mathcal M n nbsp verwendet werden Eine Zufallsmatrix ist somit das matrixwertige Analogon zu einer Zufallsvariablen Zentrale Begriffe Bearbeiten Partitionsfunktion Bearbeiten Sei E displaystyle mathcal E nbsp ein Matrix Raum z B der hermiteschen n n displaystyle n times n nbsp Matrizen H n C X M n C X X displaystyle mathcal H n mathbb C X in mathcal M n mathbb C X X nbsp und sei m d M displaystyle mu mathrm d M nbsp ein komplexes Mass auf diesem Raum welches in der Regel nicht normalisiert ist Dann nennt man das Integral Z E m d M displaystyle Z int mathcal E mu mathrm d M nbsp Partitionsfunktion und man erhalt einen Erwartungswert zur Funktion f E C displaystyle f mathcal E to mathbb C nbsp E f M 1 Z E f M m d M displaystyle mathbb E f M frac 1 Z int mathcal E f M mu mathrm d M nbsp Wignersche Matrix Bearbeiten Seien X i j 1 i lt j C displaystyle X i j 1 leq i lt j in mathbb C nbsp und Y i i 1 i R displaystyle Y i i 1 leq i in mathbb R nbsp i i d Zufallsvariablen mit gleichem Erwartungswert 0 displaystyle 0 nbsp sowie E X i j 2 0 displaystyle mathbb E X i j 2 0 nbsp und E X i j 2 1 displaystyle mathbb E X i j 2 1 nbsp Man nennt eine Zufallsmatrix H n h i j 1 i j n displaystyle H n h i j 1 leq i j leq n nbsp eine komplexe Wignersche Matrix wenn sie hermitesch ist und folgendes gilt h i j X i j n i lt j Y i i n i j displaystyle h i j begin cases X i j sqrt n amp i lt j Y i i sqrt n amp i j end cases nbsp Die Matrix wird oft mit n 1 displaystyle sqrt n 1 nbsp skaliert Manche Autoren definieren sie aber auch ohne Skalierung Sie ist ein wichtiger Typ von Zufallsmatrizen und benannt nach Eugene Wigner Wignersche Matrizen mit einer zugrundeliegenden Normalverteilung fuhren zu dem Begriff der gaussschen invarianten Ensembles Allgemeine Wignersche Matrizen sind nicht invariant Das GUE erhalt man wenn zusatzlich E Y i i 2 1 displaystyle mathbb E Y i i 2 1 nbsp gilt und die Eintrage normalverteilt sind Das GOE erhalt man wenn alle Eintrage reell und normalverteilt sind und zusatzlich E Y i i 2 2 displaystyle mathbb E Y i i 2 2 nbsp gilt Invariante Ensembles Bearbeiten Zentrale Studienobjekte sind die invarianten Ensembles welche durch die folgenden Masse auf dem entsprechenden Raum der Matrizen induziert werden P b n d A 1 Z b n e b 2 n Tr Q A d A displaystyle P beta n mathrm d A frac 1 Z beta n e frac beta 2 n operatorname Tr Q A mathrm d A nbsp wobei b displaystyle beta nbsp der Dyson Index ist und Q A displaystyle Q A nbsp das Potential Man setzt an Q displaystyle Q nbsp voraus dass Q x displaystyle Q x to infty nbsp genugend schnell wenn x displaystyle x to pm infty nbsp damit alle Momente existieren In der Regel ist Q displaystyle Q nbsp ein Polynom Man erhalt fur b 1 displaystyle beta 1 nbsp das orthogonale Ensemble OE auf dem Raum der n n displaystyle n times n nbsp reellen symmetrischen Matrizen b 2 displaystyle beta 2 nbsp das unitare Ensemble UE auf dem Raum der n n displaystyle n times n nbsp hermiteschen Matrizen b 4 displaystyle beta 4 nbsp das symplektische Ensemble SE auf dem Raum der n n displaystyle n times n nbsp hermiteschen quaternionen Matrizen Die freie Energie der unitaren Ensembles ist 3 F n 0 n 2 ln Z 2 n n 2 ln H n e n Tr Q M d M displaystyle F n 0 n 2 ln Z 2 n n 2 ln int mathcal H n e n operatorname Tr Q M mathrm d M nbsp wobei H n displaystyle mathcal H n nbsp den Raum der hermiteschen Matrizen bezeichnet Mit Hilfe der weylschen Integralformel lasst sich zeigen dass das kanonische unnormalisierte Haarsche Mass d U h displaystyle mathrm d U h nbsp auf der entsprechenden kompakten Lie Gruppe O n U n T n displaystyle mathbb O n mathbb U n mathbb T n nbsp oder U S p n T n displaystyle mathbb USp n mathbb T n nbsp folgende Darstellung zulasst d A i lt j l j l i b d L d U h displaystyle mathrm d A prod i lt j left lambda j lambda i right beta mathrm d Lambda mathrm d U h nbsp wobei d L displaystyle mathrm d Lambda nbsp das Lebesgue Mass der Eigenwerte ist Fur skalierte Eintrage des Gaussschen Ensembles erhalt man eine geschlossene Form des Wahrscheinlichkeitsmasses uber der Weyl Kammer mit l n l 1 displaystyle lambda n geq dots geq lambda 1 nbsp P b n d l 1 d l n 1 Z b n i lt j l j l i b k 1 n e b n 2 Q l k d l k displaystyle P beta n mathrm d lambda 1 dots mathrm d lambda n frac 1 widetilde Z beta n prod i lt j left lambda j lambda i right beta prod k 1 n e frac beta n 2 Q lambda k mathrm d lambda k nbsp Das Wahrscheinlichkeitsmass enthalt den Boltzmann Faktor e b n E l 1 l n displaystyle e beta n mathcal E lambda 1 dots lambda n nbsp wobei E displaystyle mathcal E nbsp die totale potentielle Energie bezeichnet E l 1 l n 1 2 j 1 n Q l k 1 n i lt j log l j l i displaystyle mathcal E lambda 1 dots lambda n frac 1 2 sum limits j 1 n Q lambda k frac 1 n sum i lt j log lambda j lambda i nbsp Die Konstante Z b n n i lt j l j l i b i 1 n e b n 2 Q l i d l i displaystyle widetilde Z beta n n left int infty infty cdots int infty infty prod i lt j left lambda j lambda i right beta prod limits i 1 n e frac beta n 2 Q lambda i mathrm d lambda i right nbsp lasst sich mit Hilfe des Selberg Integrals berechnen Gausssche Ensembles Bearbeiten Wichtige Spezialfalle der invariante Ensembles sind die Gaussschen Ensembles welche durch das Potential Q A 1 2 A 2 displaystyle Q A tfrac 1 2 A 2 nbsp mit Q l k 1 2 l k 2 displaystyle Q lambda k tfrac 1 2 lambda k 2 nbsp und die folgenden Gaussschen Masse erzeugt werden P b n d A 1 Z b n e b 4 n Tr A 2 d A displaystyle P beta n mathrm d A frac 1 Z beta n e frac beta 4 n operatorname Tr A 2 mathrm d A nbsp Man erhalt fur b 1 displaystyle beta 1 nbsp das Gausssches Orthogonales Ensemble GOE auf dem Raum der n n displaystyle n times n nbsp reellen symmetrischen Matrizen b 2 displaystyle beta 2 nbsp das Gausssches Unitares Ensemble GUE auf dem Raum der n n displaystyle n times n nbsp hermiteschen Matrizen b 4 displaystyle beta 4 nbsp das Gausssches Symplektisches Ensemble GSE auf dem Raum der n n displaystyle n times n nbsp hermiteschen quaternionen Matrizen Die Bezeichnung orthogonal unitar symplektisch bezeichnet unter welcher Matrix Konjugation die Verteilung invariant ist Beispielsweise gilt fur eine Matrix M n displaystyle M n nbsp aus dem GOE und einer Matrix O displaystyle O nbsp aus der orthogonalen Gruppe O n displaystyle mathbb O n nbsp dass M n O T M n O displaystyle M n sim O T M n O nbsp In der Quantenmechanik werden sie verwendet um Hamiltonoperatoren zu modellieren Herleitung des GUE durch Ornstein Uhlenbeck Prozesse Bearbeiten Man betrachte das System stochastischer Differentialgleichungen von Ornstein Uhlenbeck Prozesse d X i t X i t d t d W i t 1 i n d X i j t X i j t d t d W i j 1 t 1 i lt j n d Y i j t Y i j t d t d W i j 2 t 1 i lt j n displaystyle begin cases mathrm d X i t X i t mathrm d t mathrm d W i t amp quad 1 leq i leq n mathrm d X ij t X ij t mathrm d t mathrm d W ij 1 t amp quad 1 leq i lt j leq n mathrm d Y ij t Y ij t mathrm d t mathrm d W ij 2 t amp quad 1 leq i lt j leq n end cases nbsp wobei W i W i j 1 W i j 2 displaystyle W i W ij 1 W ij 2 nbsp unabhangige brownsche Bewegungen sind mit E W i t 1 W j t 2 2 D d i j d t 1 t 2 E W i j s 1 t 1 W s t s 2 t 2 D d s 1 s 2 d s t d i j d t 1 t 2 displaystyle mathbb E dot W i t 1 dot W j t 2 2D delta ij delta t 1 t 2 qquad mathbb E dot W ij sigma 1 t 1 dot W st sigma 2 t 2 D delta sigma 1 sigma 2 delta st delta ij delta t 1 t 2 nbsp und die Initialwerte X i 0 X i j 0 Y i j 0 displaystyle X i 0 X ij 0 Y ij 0 nbsp beliebig sind Definiert man nun eine hermitesche Zufallsmatrix H h i j t n 1 i j n displaystyle H h ij t sqrt n 1 leq i j leq n nbsp fur t 0 displaystyle t geq 0 nbsp mit h i j t X i t i j X i j t i Y i j t i lt j displaystyle h ij t begin cases X i t amp quad i j X ij t mathrm i Y ij t amp quad i lt j end cases nbsp und bezeichnet mit f t d H displaystyle f t mathrm d H nbsp das zugehorige Wahrscheinlichkeitsmass dann gilt fur D 2 displaystyle D 2 nbsp und t displaystyle t to infty nbsp f t d H P 2 n d H displaystyle f t mathrm d H propto P 2 n mathrm d H nbsp wobei P 2 n d H displaystyle P 2 n mathrm d H nbsp das GUE bezeichnet 4 Zirkulare Ensembles Bearbeiten Man erhalt das Zirkulare Unitare Ensemble ZUE durch das haarsche Mass auf dem Raum der unitaren Matrizen Das Zirkulare Orthogonale Ensemble ZOE erhalt man durch das haarsche Mass auf dem Raum der symmetrischen unitaren Matrizen Das Zirkulare Symplektische Ensemble ZSE erhalt man durch das haarsche Mass auf dem Raum der selbst dualen unitaren Quaternionen Matrizen Die Dichte der Eigenwerte l j e i 8 j displaystyle lambda j e i theta j nbsp der Zirkularen Ensembles ist p b l 1 l n 1 Z b n c i r c i lt j e i 8 j e i 8 i b p lt 8 l p displaystyle p beta lambda 1 dots lambda n frac 1 Z beta n circ prod i lt j left e i theta j e i theta i right beta quad pi lt theta l leq pi nbsp wobei b 1 displaystyle beta 1 nbsp fur das ZOE b 2 displaystyle beta 2 nbsp fur das ZUE und b 4 displaystyle beta 4 nbsp fur das ZSE gilt 5 b Ensembles und Dysons Threefolded Way Bearbeiten Man spricht von Dysons b displaystyle beta nbsp Ensemble da Freeman Dyson in seiner wissenschaftlichen Schrift The Threefolded Way 6 diese b 1 2 4 displaystyle beta 1 2 4 nbsp Klassifizierungen der Zufallsmatrizen herleitete basierend auf physikalisch moglichen Zeitumkehr Eigenschaften der Quantenmechanik orthogonal unitar symplektisch Der Fall b 3 displaystyle beta 3 nbsp ist aufgrund des Satzes von Frobenius nicht moglich Neben den Gaussschen Ensembles spielen auch die b displaystyle beta nbsp Wishart Laguerre Ensembles und die b displaystyle beta nbsp Jacobi Manova Ensembles eine zentrale Rolle in der Theorie der Zufallsmatrizen Es ist ublich nur von Laguerre Ensembles bzw Jacobi Ensembles zu sprechen statt von Wishart bzw Manova EnsemblesAllgemeine b displaystyle beta nbsp spielen in der klassischen Theorie der Zufallsmatrizen eine untergeordnete Rolle Theorie der Zufallsmatrizen Bearbeiten Die Theorie der Zufallsmatrizen befasst sich weniger mit einer konkreten Zufallsmatrix sondern mit dem Matrizenraum dahinter Konkret geht es um Wahrscheinlichkeitsmasse auf Matrixraumen und Lie Gruppen dies erklart den Begriff Ensembles Ein klassisches Problem der Theorie der Zufallsmatrizen ist das Finden einer multivariaten Wahrscheinlichkeitsdichte fur die Eigenwerte unterschiedlicher Matrix Ensembles Eine der fruhsten Arbeiten stammt von Dyson welcher eine geschlossene Form fur eine grosse Menge von Matrizen fand abhangig von der zugrundeliegenden Symmetrie der Matrizen und Wahrscheinlichkeitsverteilung Die Spektraleigenschaften von grossen Zufallsmatrizen haben universelle Eigenschaften und man kann beim Studium komplizierter deterministischer Operatoren wie zum Beispiel dem Dirac Operator aus der Physik diese Operatoren mit Zufallsmatrizen ersetzen und die Theorie der Zufallsmatrizen anwenden Beim Studium von Integralen uber Matrix Raume verwendet man zum Teil Resultate aus der Theorie der Lie Gruppen und Lie Algebren Auch die freie Wahrscheinlichkeitstheorie von Voiculescu ist von Relevanz fur grosse Zufallsmatrizen Generell untersucht man Matrizen mit bestimmten Symmetrie Eigenschaften z B hermitesche und hat bestimmte stochastische Anforderungen an die Wahrscheinlichkeitsverteilung auf den Raum jener Matrizen z B obere Dreiecksmatrix unabhangig Des Weiteren interessiert man sich vor allem fur die Spektraltheorie und dessen asymptotisches Verhalten wenn die Dimension N displaystyle N to infty nbsp Die Spektraltheorie ist engverbunden mit der Theorie der Punktprozesse da die Eigenwerte einen zufalligen Punktprozess formen Bei vielen Ensembles taucht in der gleichen Region derselbe Punktprozess in unendlicher Dimension auf Universalitat Matrix wertige Funktionen wie die Determinante oder die Spur konnen nicht einfach auf unendlich dimensionale Matrizen ubertragen werden Fur bestimmte Operatoren lasst sich aber mit der abstrakten Fredholmtheorie eine Erweiterung auf unendlich dimensionale separable Hilbertraume uber die aussere Algebra finden Es lassen sich Determinanten fur Operatoren aus den Schatten von Neumann Klassen definieren Definiert man die Eintrage der Matrix als Brownsche Bewegungen so lasst sich auch das matrixwertige Analogon eines stochastischen Prozess bilden und die Theorie der stochastischen Analysis und die Martingal Theorie ist anwendbar siehe Dysons brownsche Bewegung und Wishart Prozess Spektraltheorie der Zufallsmatrizen BearbeitenSind die Eintrage einer hermiteschen Zufallsmatrix M n displaystyle M n nbsp von der Grosse O n displaystyle mathcal O sqrt n nbsp so konvergiert das empirische Spektralmass n n l M n n 1 j 0 n d l l j M n displaystyle nu n lambda M n n 1 sum limits j 0 n delta lambda lambda j M n nbsp wobei d displaystyle delta nbsp das Dirac Delta bezeichnet Da die zufalligen Ensembles Punktprozesse sind kann man die n displaystyle n nbsp Punkt Korrelationsfunktion fur die Eigenwerte l 1 l N displaystyle lambda 1 dots lambda N nbsp herleiten Sei f displaystyle varphi nbsp eine Testfunktion und definiere das Funktional E N f n 0 N 1 n N n E f l 1 f l n displaystyle E N varphi sum limits n 0 N 1 n binom N n mathbb E left varphi lambda 1 cdots varphi lambda n right nbsp Dann ist die n displaystyle n nbsp Punkt Korrelationsfunktion folgende ausgewertete Funktionalableitung 7 1 n d n d f l 1 d f l n E N f f 0 R N n l 1 l n displaystyle 1 n frac delta n delta varphi lambda 1 cdots delta varphi lambda n E N varphi Biggr varphi 0 R N n lambda 1 dots lambda n nbsp Mit dem Darstellungssatz von Frechet Riesz lasst sich Konvergenz im Erwartungswert fur f C c R displaystyle varphi in C c mathbb R nbsp definieren f d E n n l M n E f d n n l M n displaystyle int varphi mathrm d mathbb E nu n lambda M n mathbb E int varphi mathrm d nu n lambda M n nbsp Globale Situation Bearbeiten Eines der wichtigsten Ergebnisse ist das sogenannte Wignersche Halbkreisgesetz siehe Eugen Wigner Es besagt dass das skalierte empirische Spektralmass n n displaystyle nu n nbsp einer Wignerischen Zufallsmatrix in der Physik bekannt als die sogenannte Zustandsdichte einer charakteristischen Halbkreis Verteilung genugt Das Variationsproblem der Verteilung der Eigenwerte Bearbeiten Allgemeiner handelt es sich bei der Grenzwertverteilung der Eigenwerte um die Losung eines Variationsproblem Definiere den Raum der Masse M 1 R n n 0 R d n 1 displaystyle M 1 mathbb R left nu nu geq 0 int mathbb R mathrm d nu 1 right nbsp und betrachte das Funktional E Q inf n M 1 R x y ln x y d n x d n y Q x d n x displaystyle E Q inf limits nu in M 1 mathbb R int int x neq y ln x y mathrm d nu x mathrm d nu y int Q x mathrm d nu x nbsp Das Funktional erklart sich durch die Integralschreibweise der totalen potentiellen Energie x y ln x y 1 n i 1 n Q x displaystyle sum x neq y ln x y 1 n sum i 1 n Q x nbsp bezuglich des empirischen Spektralmasses n n displaystyle nu n nbsp Fur E Q displaystyle E Q nbsp wird ein eindeutiges Equilibriummass n Q displaystyle nu Q nbsp durch die Euler Lagrange Variationsbedingung fur eine reelle Konstante l displaystyle l nbsp 3 2 R log x y d n y Q x l x J displaystyle 2 int mathbb R log x y mathrm d nu y Q x l quad x in J nbsp 2 R log x y d n y Q x l x R J displaystyle 2 int mathbb R log x y mathrm d nu y Q x leq l quad x in mathbb R setminus J nbsp definiert wobei J j 1 q a j b j displaystyle J bigcup limits j 1 q a j b j nbsp der Trager des Masses ist und q Q x 2 2 Q x Q y x y d n Q y displaystyle q left frac Q x 2 right 2 int frac Q x Q y x y mathrm d nu Q y nbsp Das Equilibirummass n Q displaystyle nu Q nbsp besitzt folgende Radon Nikodym Dichte d n Q x d x 1 p q x displaystyle frac mathrm d nu Q x mathrm d x frac 1 pi sqrt q x nbsp Beispiel Wignersche Halbkreis Bearbeiten Im Fall des GUE konvergiert das zufallige Mass schwach in Wahrscheinlichkeit gegen die deterministische Verteilung s d x 1 2 p 4 x 2 d x displaystyle sigma mathrm d x frac 1 2 pi sqrt 4 x 2 mathrm d x nbsp Es gilt fur eine Funktion f C b R displaystyle f in C b mathbb R nbsp und e gt 0 displaystyle varepsilon gt 0 nbsp lim n P n n f s f gt e 0 displaystyle lim limits n to infty P langle nu n f rangle langle sigma f rangle gt varepsilon 0 nbsp Der Satz kann mit Mitteln der Kombinatorik und der Momentmethode bewiesen werden Fur eine Zufallsvariable X s d x displaystyle X sim sigma mathrm d x nbsp gilt dass E X k d k 2 Z C k 2 displaystyle mathbb E left X k right delta k 2 in mathbb Z C k 2 nbsp wobei C n displaystyle C n nbsp die Catalan Zahlen sind Durch die oben erwahnte Equilibriummass Methode der statistischen Mechanik gibt es eine Verbindung zur Theorie der grossen Abweichungen Einen analytischen konstruktiven Beweis ergibt sich uber die Stieltjes Transformation Fur Wishart bzw Laguerre Matrizen konvergiert das empirische Spektralmass n n displaystyle nu n nbsp gegen die Martschenko Pastur Verteilung und fur MANOVA bzw Jacobi Matrizen gegen die Kesten Mckey Verteilung Fur quadratische Zufallsmatrizen A n displaystyle A n nbsp mit i i d komplexen Eintragen a i j displaystyle a ij nbsp mit E a i j 0 displaystyle mathbb E a ij 0 nbsp und Var a i j 1 displaystyle operatorname Var a ij 1 nbsp gilt das Kreisgesetz Tao Vu 8 welches besagt dass n n displaystyle nu n nbsp gegen 8 d x 1 p 1 x 2 y 2 1 d x d y displaystyle theta mathrm d x frac 1 pi 1 x 2 y 2 leq 1 mathrm d x mathrm d y nbsp konvergiert Man spricht von Universalitat weil die Satze unabhangig von der zugrundeliegenden Verteilung sind Lokale Situation Bearbeiten LimitverhaltenLokal ergibt sich bei Skalierung ein Punktprozess fur die Eigenwerte Der Fall b 2 displaystyle beta 2 nbsp von hermiteschen Matrizen ist signifikant einfacher Man kann mittels der Theorie der orthogonale Polynome eine determinantale Form fur die Korrelationsfunktion finden welche dann zu Fredholm Determinanten von Integraloperatoren fuhren Die Falle b 1 displaystyle beta 1 nbsp und b 4 displaystyle beta 4 nbsp lassen sich mit Quaternionen Determinanten und schief orthogonalen Polynome losen 9 Es gilt fur die n displaystyle n nbsp Punkt Korrelationsfunktion R N n l 1 l n N N n R l n 1 R l N p b l 1 l N j n 1 N d l j displaystyle R N n lambda 1 dots lambda n frac N N n int mathbb R lambda n 1 cdots int mathbb R lambda N p beta lambda 1 dots lambda N prod limits j n 1 N mathrm d lambda j nbsp wobei p b l 1 l n displaystyle p beta lambda 1 dots lambda n nbsp die multivariate Dichte der Eigenwerte ist Fur das GUE erhalt man einen determinantal point process ein einfacher Punktprozess mit Kern bezuglich eines Masses m displaystyle mu nbsp dessen R N n displaystyle R N n nbsp existiert so dass fur alle n 1 displaystyle n geq 1 nbsp gilt R N n d x 1 d x n det K x i x j 1 i j n m d x 1 d x n displaystyle R N n mathrm d x 1 dots mathrm d x n operatorname det K x i x j 1 leq i j leq n mu mathrm d x 1 dots mathrm d x n nbsp Skaliert man den Integralkern konvergiert dieser entweder zu dem Sinus oder Airy Kern Die benotigten asymptotischen Entwicklungen konnen mittels der nicht trivialen Methode des steilsten Anstiegs gezeigt werden asymptotische Entwicklungen vom Plancherel Rotach Typ K Sine x y sin p x y p x y displaystyle K text Sine x y tfrac sin pi x y pi x y nbsp Die Wahrscheinlichkeit dass eine kompakte Menge V R displaystyle V subset mathbb R nbsp keine unskalierte Eigenwerte l i i n displaystyle lambda i i leq n nbsp enthalt lasst sich als Fredholm Determinante formulieren Gaudin Mehta lim n P n l 1 n l n V 1 k 1 1 k k V V det K Sine x i x j i j 1 k j 1 k d x j displaystyle lim limits n to infty P sqrt n lambda 1 dots sqrt n lambda n not in V 1 sum k 1 infty frac 1 k k int V cdots int V operatorname det K text Sine x i x j i j 1 k prod limits j 1 k mathrm d x j nbsp Universalitat im Hauptteil2010 zeigten Erdos Ramirez Schlein Tao Vu Yau fur wignerische Matrizen mit subexponentialer Abnahme Universalitat des Sinus Kern 10 RandBetrachtet man den Rand des Spektrums so erhalt man einen Airy Prozess und bekommt die Tracy Widom Verteilung mit Kern K Airy x y Ai x Ai y Ai x Ai y x y displaystyle K text Airy x y tfrac operatorname Ai x operatorname Ai y operatorname Ai x operatorname Ai y x y nbsp wobei Ai displaystyle operatorname Ai nbsp die Airy Funktion bezeichnet Fur das GSE und GOE erhalt man eine Verallgemeinerung ein sogenannter pfaffian point processes Im Falle des Laguerre Ensembles ergibt sich bei dem hard edge harten Rand ein Bessel Prozess und bei dem soft edge weichen Rand ein Airy Prozess Geschichte BearbeitenBereits 1928 untersuchte John Wishart als einer der ersten die Zufallsmatrizen die bei einer standard multivariaten normalverteilten Stichprobe entstehen die Kovarianzmatrix Dies fuhrte zu der Wishart Verteilung die matrixvariate Verallgemeinerung der x2 Verteilung bzw Gamma Verteilung In den 1950er untersuchte Eugene Wigner die Verteilung zwischen benachbarten Energieniveaus von schweren Atomkernen Das Energieniveau wird durch die Eigenwerte des Hamiltonian der zeitunabhangigen Schrodingergleichung beschrieben H PS E PS displaystyle operatorname hat H Psi rangle E Psi rangle nbsp Fur schwere Atomkerne ist dieses Problem zu komplex um es theoretisch zu losen deshalb kam Wigner auf die Idee dieses Problem als statistisches Problem zu losen und stattdessen die Spektraldichte von grossen endlichen Zufallsmatrizen zu untersuchen Empirische Daten aus Experimenten zeigten dass die Verteilung von der Form w d x C b x b e k b x 2 d x displaystyle omega mathrm d x C beta x beta e k beta x 2 mathrm d x nbsp sein musste und somit das Energieniveau korreliert ist da sonst eine Poisson Verteilung zugrunde liegen sollte und es erklarte auch das Phanomen dass sich die Energieniveaus gegenseitig abstiessen Dieses Resultat wird als Wigners Vermutung englisch Wigner s surmise bezeichnet Die Konstanten C b k b displaystyle C beta k beta nbsp sind von b 1 2 4 displaystyle beta in 1 2 4 nbsp abhangig und b displaystyle beta nbsp beschreibt die zugrundeliegende Symmetrie der Atomkerne unter Zeitumkehr und Spinrotation Wigner postulierte dass die Abstande zwischen den Linien des Spektrums den Abstanden der Eigenwerte einer Zufallsmatrix gleichen Aus den 1960ern stammen bedeutende Arbeiten zur mathematischen Theorie der Zufallsmatrizen von Gaudin Mehta und Dyson Parallel dazu entwickelte sich auch wichtige Arbeiten zu den Kovarianzmatrizen Die traditionelle Ausgangslage der Statistik hat eine kleine fixe Anzahl p displaystyle p nbsp von Parametern und n displaystyle n to infty nbsp Observationen Die Theorie der Zufallsmatrizen hat sich aus der Situation entwickelt wenn p displaystyle p nbsp sehr gross ist und man interessiert sich auch fur die Falle wenn n p displaystyle n p to infty nbsp In den 1970ern entdeckte Montgomery und Dyson eine Verbindung zwischen den Zufallsmatrizen und der Zahlentheorie respektive zwischen schweren Atomkernen und den kritischen Nullstellen der riemannschen Zeta Funktion Anwendungen BearbeitenStatistik Bearbeiten In der Multivariaten Statistik insbesondere die Wishart Verteilung ist zentraler Untersuchungsgegenstand Bei der Analyse von hoch dimensionalen Daten bei der Auswahl von statistischen Modellen der Hauptkomponentenanalyse PCA Untersuchungsgegenstand sind somit insbesondere Probleme der multivariaten Statistik bei grossen Daten mit vielen Parametern Im Maschinellen Lernen und Deep Learning gibt es auch Anwendungen 11 Neue Forschungsergebnisse haben Equilibriumsmasse fur einige kunstliche neuronale Netze gezeigt welche fur die Beschleunigung von neuronalen Netzwerken genutzt werden konnen Physik Bearbeiten Im Fall eines ungeordneten physikalischen Systems z B bei sog amorphem Material sind die betreffenden Matrix Elemente Zufallsgrossen Die Physik dieser Systeme kann im Wesentlichen durch die Kenngrossen der jeweiligen Matrizen erfasst werden z B durch Mittelwert und Schwankung der jeweiligen Grosse Beispielsweise kann die Warmeleitfahigkeit eines kristallinen Festkorpers direkt aus der sogenannten dynamischen Matrix der Teilchen Teilchen Wechselwirkung des Kristallgitters berechnet werden Anwendungen u a bei magnetischen Systemen z B bei Multilagensystemen magnetischer Dunnschicht Systeme 12 dem Quanten Hall Effekt 13 14 sogenannter Quanten Dots 15 und bei Supraleitern 16 Anwendungen in der Kernphysik betreffen u a das oben erwahnte Gauss sche orthogonale das unitare und das symplektische Ensemble Energiespektrum und Wirkungsquerschnitt eines Atomkerns sind zwar extrem komplex aber gerade deshalb der Theorie des sogenannten chaotischen Verhaltens zuganglich sowie das sogenannte Quantenchaos und die mesoskopische Physik 17 Ferner gibt es Anwendungen in der sogenannten Quantengravitation bei zweidimensionalen Systemen 18 Weitere Anwendungen Bearbeiten In der Mathematik die L Funktionen von Dirichlet und andere Ferner gibt es zahlreiche Anwendungen in der analytischen Zahlentheorie in der Operatoralgebra 19 und in der sogenannten freien Wahrscheinlichkeitstheorie 20 21 Insbesondere liefert es durch Montgomerys Paar Korrelation Vermutung einen neuen Zugang zur Riemannschen Vermutung Des Weiteren gibt es auch Anwendungen in der Graphentheorie und der Knotentheorie In der Portfoliotheorie bei der Analyse der Kovarianzmatrix von Portfolio Renditen 22 In der Signalverarbeitung und fur drahtlose Netzwerke 23 In der Genetik zur Klassifizierung von RNA Strukturen wie den Pseudoknoten Aus aktuellen Untersuchungen ergibt sich als Vermutung dass die Theorie der Zufallsmatrizen zu Verbesserungen bei Suchmaschinen im Web fuhren konnte 24 Literatur BearbeitenGreg Anderson Alice Guionnet Ofer Zeitouni An Introduction to Random Matrices Cambridge University Press 2010 Alice Guionnet Large Random Matrices Lectures on Macroscopic Asymptotics Springer Verlag 2009 Persi Diaconis Patterns in Eigenvalues The 70th Josiah Willard Gibbs Lecture In American Mathematical Society Bulletin New Series 2003 S 155 178 doi 10 1090 S0273 0979 03 00975 3 Persi Diaconis What Is a Random Matrix In Notices of the American Mathematical Society 2005 S 1348 1349 ISSN 0002 9920 M L Mehta Random matrices 3 Auflage In Pure and Applied Mathematics 142 Elsevier Academic Press Amsterdam 2004 xviii 688 S Guhr Muller Groening Weidenmuller Random Matrix Theories in Quantum Physics Common Concepts In Physics Reports Band 299 1998 S 189 425 arxiv cond mat 9707301 Alan Edelman N Raj Rao Random Matrix Theory PDF In Acta Numerica Band 14 2005 S 233 297 Terence Tao Topics in Random Matrix Theory American Mathematical Society 2012Weblinks BearbeitenRandom Matrix at MathWorld RMTool A MATLAB based Random Matrix Calculator Vorlesung von Terence Tao Kriecherbauer Eine kurze und selektive Einfuhrung in die Theorie der Zufallsmatrizen PDF 402 kB Uni Bochum 2008 Einzelnachweise und Anmerkungen Bearbeiten Elizabeth Meckes The Random Matrix Theory of the Classical Compact Groups 1 Auflage Cambridge University Press 2019 ISBN 978 1 108 41952 9 doi 10 1017 9781108303453 Greg W Anderson Alice Guionnet Ofer Zeitouni An Introduction to Random Matrices Cambridge University Press 2009 a b John Harnad Random Matrices Random Processes and Integrable Systems Hrsg Springer ISBN 978 1 4614 2877 0 Yan V Fyodorov Introduction to the Random Matrix Theory Gaussian Unitary Ensemble and Beyond doi 10 48550 ARXIV MATH PH 0412017 Peter J Forrester Log Gases and Random Matrices In Princeton University Press Hrsg London Mathematical Society Monographs LMS 34 Band 34 Freeman Dyson The Threefold Way Algebraic Structure of Symmetry Groups and Ensembles in Quantum Mechanics Abgerufen am 23 Mai 2021 Leonid Pastur und Mariya Shcherbina Eigenvalue distribution of large random matrices Hrsg American Mathematical Society 2011 2010 gezeigt durch Terence Tao und Van H Vu G Mahoux und M L Mehta A method of integration over matrix variables IV In EDP Sciences Hrsg Journal de Physique I Nr 1 1991 S 1093 1108 Laszlo Erdos Jose Ramirez and Benjamin Schlein and Terence Tao and Van Vu and Horng Tzer Yau Bulk universality for Wigner hermitian matrices with subexponential decay arxiv 0906 4400 math https tel archives ouvertes fr tel 02397287 document VS Rychkov S Borlenghi H Jaffres A Fert X Waintal Spin Torque and Waviness in Magnetic Multilayers A Bridge Between Valet Fert Theory and Quantum Approaches In Phys Rev Lett 103 Jahrgang Nr 6 August 2009 S 066602 doi 10 1103 PhysRevLett 103 066602 PMID 19792592 DJE Callaway Random Matrices Fractional Statistics and the Quantum Hall Effect In Phys Rev B Condens Matter 43 Jahrgang Nr 10 April 1991 S 8641 8643 doi 10 1103 PhysRevB 43 8641 PMID 9996505 M Janssen K Pracz Correlated Random Band Matrices Localization delocalization Transitions In Phys Rev E Stat Phys Plasmas Fluids Relat Interdiscip Topics 61 Jahrgang 6 Pt A Juni 2000 S 6278 86 doi 10 1103 PhysRevE 61 6278 PMID 11088301 DM Zumbuhl JB Miller CM Marcus K Campman AC Gossard Spin orbit Coupling Antilocalization and Parallel Magnetic Fields in Quantum Dots In Phys Rev Lett 89 Jahrgang Nr 27 Dezember 2002 S 276803 doi 10 1103 PhysRevLett 89 276803 PMID 12513231 SR Bahcall Random Matrix Model for Superconductors in a Magnetic Field In Phys Rev Lett 77 Jahrgang Nr 26 Dezember 1996 S 5276 5279 doi 10 1103 PhysRevLett 77 5276 PMID 10062760 D Sanchez M Buttiker Magnetic field Asymmetry of Nonlinear Mesoscopic Transport In Phys Rev Lett 93 Jahrgang Nr 10 September 2004 S 106802 doi 10 1103 PhysRevLett 93 106802 PMID 15447435 F Franchini VE Kravtsov Horizon in Random Matrix Theory the Hawking Radiation and Flow of Cold Atoms In Phys Rev Lett 103 Jahrgang Nr 16 Oktober 2009 S 166401 doi 10 1103 PhysRevLett 103 166401 PMID 19905710 arxiv math OA 0412545 James A Mingo Roland Speicher Free Probability and Random Matrices Fields Institute Monographs Bd 35 Springer Verlag New York 2017 arxiv math OA 0501238 http web mit edu sea06 agenda talks Harding pdf Antonia M Tulino Sergio Verdu Random Matrix Theory and Wireless Communications Now 2004 newscientist comSpezielle Matrizen in der Statistik Datenmatrix Produktsummenmatrix Pradiktionsmatrix residuenerzeugende Matrix zentrierende Matrix Kovarianzmatrix Korrelationsmatrix Prazisionsmatrix Gewichtsmatrix Restriktionsmatrix Fisher Informationsmatrix Bernoulli Matrix Leslie Matrix Zufallsmatrix Ubergangsmatrix Abgerufen von https de wikipedia org w index php title Zufallsmatrix amp oldid 239487594