Wolfram Language

Matrices aleatorias

Densidad espectral de una matriz

La densidad espectral conjunta de muchas distribuciones de matriz tiene una forma correspondiente universal de límite. Para las matrices aleatorias hermitianas con entradas independientes, esta es la ley de semicírculo de Wigner.

Para los conjuntos de Gauss, la densidad espectral de las matrices tiene una forma cerrada para la dimensión de matriz finita, relacionada con las funciones propias del oscilador armónico cuántico.

Use MatrixPropertyDistribution para representar el espectro escalado del conjunto unitario de Gauss y proporcione la expresión de forma cerrada de su densidad espectral conjunta.

In[1]:=
Click for copyable input
scaledSpectrum\[ScriptCapitalD][n_] := MatrixPropertyDistribution[ Eigenvalues[\[Lambda]]/(2 Sqrt[n]), \[Lambda] \[Distributed] GaussianUnitaryMatrixDistribution[n]];
In[2]:=
Click for copyable input
spectralPDF[n_Integer, \[Lambda]_] := Sqrt[2/(\[Pi] n)] Exp[-2 n \[Lambda]^2] \!\( \*UnderoverscriptBox[\(\[Sum]\), \(j = 0\), \(n - 1\)] \*FractionBox[ SuperscriptBox[\(HermiteH[j, \*SqrtBox[\(2\ n\)]\ \[Lambda]]\), \(2\)], \( \*SuperscriptBox[\(2\), \(j\)]\ \(j!\)\)]\)

Para la distribución de matriz pequeña, hay un patrón oscilatorio de características, cuyo número de densidad máxima es igual al tamaño de la matriz.

In[3]:=
Click for copyable input
scaledSpectra = Flatten[RandomVariate[scaledSpectrum\[ScriptCapitalD][#], 10^5]] & /@ {3, 4, 5};
muestre la entrada completa de Wolfram Language
In[4]:=
Click for copyable input
Row@MapThread[ Show[Plot[spectralPDF[#1, x], {x, -1.5, 1.5}, PlotTheme -> "Scientific", PlotStyle -> ColorData[97, 1], PlotLegends -> Placed["n = " <> ToString@#1, Above]], Histogram[#2, {0.05}, PDF]] &, {{3, 4, 5}, scaledSpectra}]
Out[4]=

En el límite de la dimensión grande, la densidad converge con WignerSemicircleDistribution.

In[5]:=
Click for copyable input
n = 250; scaledSpectrum = Flatten[RandomVariate[scaledSpectrum\[ScriptCapitalD][n], 10^2]];
muestre la entrada completa de Wolfram Language
In[6]:=
Click for copyable input
Show[Histogram[scaledSpectrum, {0.05}, PDF, PlotTheme -> "Detailed"], Plot[PDF[WignerSemicircleDistribution[1], x], {x, -1.5, 1.5}, Exclusions -> None, PlotTheme -> "Scientific", PlotLegends -> None, PlotStyle -> ColorData[97, 1]], ImageSize -> Medium]
Out[6]=

Ejemplos relacionados

de en fr ja ko pt-br ru zh