Wolfram Language

Matrices aléatoires

Propriétés de distributions de matrice

Les statistiques de plus faibles dimensions dérivées des matrices aléatoires jouent un rôle important dans la caractérisation des ensembles de matrice. Dans diverses situations de limites, les distributions de ces statistiques tombent dans différentes classes d'universalité. MatrixPropertyDistribution fournit un accès pratique à l'échantillon et le calcul de l'approximation numérique de ces propriétés dérivées.

Prenez un échantillon des deux plus grandes valeurs propres à partir d'un ensemble unitaire de Gauss.

In[1]:=
Click for copyable input
dist = MatrixPropertyDistribution[With[{spectrum = Eigenvalues[x]}, {Max[spectrum], RankedMax[spectrum, 2]}], x \[Distributed] GaussianUnitaryMatrixDistribution[50]];
In[2]:=
Click for copyable input
RandomVariate[dist]
Out[2]=

Visualisez la distribution conjointe des deux plus grandes valeurs propres sur la base du résultat échantillonné.

In[3]:=
Click for copyable input
sample = RandomVariate[dist, 10^4];
Montrer l'entrée complète de Wolfram Language
In[4]:=
Click for copyable input
SmoothHistogram3D[sample, "Scott", "PDF", PlotTheme -> "Detailed", ImageSize -> Medium]
Out[4]=

Prenez léchantillonnage du déterminant de matrices à partir d'un ensemble orthogonal gaussien et comparez la distribution empirique avec l'expression de forme fermée.

In[5]:=
Click for copyable input
dist = MatrixPropertyDistribution[Det[x], x \[Distributed] GaussianOrthogonalMatrixDistribution[2]]; dets = RandomVariate[dist, 10^6];
Montrer l'entrée complète de Wolfram Language
In[6]:=
Click for copyable input
empirical = Histogram[dets, {0.1}, PDF, PlotTheme -> "Detailed", ImageSize -> Medium]; detpdf[y_] := 1/Sqrt[2] Exp[ y] Piecewise[{{Erfc[Sqrt[2 y]], y >= 0}, {1, y < 0}}]; analytical = Plot[detpdf[y], {y, -5, 4}, Exclusions -> None, PlotLegends -> None, PlotTheme -> "Detailed", ImageSize -> Medium, Filling -> Axis]; Show[empirical, analytical]
Out[6]=

Approximez la moyenne du déterminant par au moyen de la méthode de Monte Carlo et comparez-la avec la valeur actuelle.

In[7]:=
Click for copyable input
{N@Mean[dist], Integrate[x detpdf[x], {x, -Infinity, Infinity}]}
Out[7]=

Exemples connexes

de en es ja ko pt-br ru zh