Usuari:Freutci/semicercle

Distribucó del semicercle de Wigner

La distribució del semicercle de Wigner, o senzillament distribució del semicercle, va ser introduïda pel físic hongarès-americà Eugene Wigner (1902-1995) en relació amb la distribució asimptòtica dels valors propis de certes matrius aleatòries simètriques.

Funció de densitat i moments

La funció de densitat de la distribució del semicercle és

$$f(x)={\begin{cases}{\frac {1}{2\pi }}{\sqrt {4-x^{2}}},&{\text{si}}\ x\in [-2,2],\\0,&{\text{en cas contrari.}}\\\end{cases}}$$

Correspon a una distribució beta amb quatre paràmetres amb $\alpha =\beta =3/2$, i $a=-2\ {\text{i}}\ b=2$ . El gràfic d'aquesta densitat és una semi-el·lipse, vegeu la la Figura 1.

Figura 1. Funció de densitat de la distribució del semicercle

Aquesta distribució té moments de tots els ordres. Si designem per $m_{k}=\int _{-\infty }^{\infty }x^{k}\,f(x)\,dx$ el moment d'ordre $k$, tenim que

$$m_{k}={\begin{cases}0,&{\text{si}}\ k\ {\text{és senar}},\\C_{n},&{\text{si}}\ k=2n,\end{cases}}$$

on

$$C_{n}={\frac {(2n)!}{(n+1)!\,n!}}$$

és l'$n$-èsim nombre de Catalan. Compareu aquests moments amb els d'una distribució normal estàndard. En particular, si $X$ té distribució del semicercle, llavors

$$E[X]=0\quad {\text{i}}\quad {\text{Var}}(X)=1.$$

Prova

Els moments d'ordre senar són zero degut a la simetria respecte 0 de la funció de densitat. Pel cas parell, es redueix la integral a una funció beta. En primer lloc, per simetria,

$$m_{2n}={\frac {1}{2\pi }}\int _{-2}^{2}x^{2n}{\sqrt {4-x^{2}}}\,dx={\frac {1}{\pi }}\int _{0}^{2}x^{2n}{\sqrt {4-x^{2}}}\,dx.$$

Ara es fa en canvi de variables $x=2{\sqrt {y}}$. Aleshores, $m_{2n}={\frac {2^{2n+1}}{\pi }}\int _{0}^{1}y^{m-1/2}(1-y)^{1/2}\,dy={\frac {2^{2n+1}}{\pi }}\,{\text{B}}(m+1/2,3/2)={\frac {2^{2n+1}}{\pi }}\,{\frac {\Gamma (m+1/2)\,\Gamma (3/2)}{\Gamma (m+2)}}=C_{n}.$

Funció generatriu de moments, funció característica i transformada de Stieltjes

La distribució del semicercle té funció generatriu de moments en tot $\mathbb {R}$,

$$M(t)={\frac {1}{t}}I_{1}(2t),\quad t\in \mathbb {R} ,$$

on $I_{1}$ és la funció de Bessel modificada amb $\nu =1$.

La funció característica val

$$\varphi (t)={\frac {1}{t}}J_{1}(2t),\quad t\in \mathbb {R} ,$$

on $J_{1}$ és la funció de Bessel amb $\nu =1$.

Prova

La funció generatriu de moments és

$$M(t)={\frac {1}{2\pi }}\int _{-2}^{2}e^{tx}{\sqrt {4-x^{2}}}\,dx={\frac {2}{\pi }}\int _{-1}^{1}e^{2ux}{\sqrt {1-u^{2}}}\,du={\frac {1}{t}}I_{1}(2t),$$

on hem aplicat el canvi $x=2u$ i després la relació amb la funció modificada de Bessel amb $\nu =1$
La funció característica es calcula de manera similar:

$$\varphi (t)={\frac {1}{2\pi }}\int _{-2}^{2}e^{itx}{\sqrt {4-x^{2}}}\,dx={\frac {1}{2\pi }}\int _{-2}^{2}\cos(tx){\sqrt {4-x^{2}}}\,dx+{\frac {1}{2\pi }}i\int _{-2}^{2}\sin(tx){\sqrt {4-x^{2}}}\,dx={\frac {1}{t}}J_{1}(2t),$$

ja que la $\int _{-2}^{2}\sin(tx){\sqrt {4-x^{2}}}\,dx=0$ per raons de simetria, i fent el mateix canvi de variables que abans a la integral amb el cosinus, s'obté una integral que es redueix a una expressió en termes de la funció $J_{1}$ de Bessel amb $\nu =1$ .

Per estudiar les propietats de les matrius aleatòries, s'utlitza la transformada de Stieltjes. La transformada de Stieltjes d'una probabilitat $p$ a $\mathbb {R}$ es defineix per

$$s(z)=\int _{\mathbb {R} }{\frac {1}{x-z}}\,dp(x),\quad z\in \mathbb {C} \setminus {\text{sup}}(p),$$

on ${\text{sup}}(p)$ és el suport de $p$; en particular, $s$ està ben definida en $\mathbb {C} \setminus \mathbb {R}$ . Si la probabilitat $p$ té densitat $g$, aleshores

$$s(z)=\int _{\mathbb {R} }{\frac {g(x)}{x-z}}\,dx.$$

La transformada de Stieltjes determina unívocament una distribució de probabilitat i té molt bones propietats respecte la convergència feble de probabilitats .
La distribució del semicercle té transformada de Stieltjes donada per

$$s(z)=-{\frac {1}{2}}{\big (}z-{\sqrt {z^{2}-4}}),\quad z\in \mathbb {C} \setminus [-2,2].$$

Per a una demostració mitjançant càlcul de residus, vegeu Bai and Silverstein .

La distribució del semicercle i probabilitats lliures

La distribució dels semicercle té, en probabilitats lliures (free probability), en molts aspectes un paper anàleg al de la distribució normal. Per exemple, en el teorema central del límit per a variables lliures el límit és una distribució del semicercle .

Distribució del semicercle amb paràmetres

Diversos autors introdueixen una distribució del semicercle amb un o dos paràmetres. Cal dir que les notacions no són universals. La distribució del semicercle amb paràmetre $a\in \mathbb {R}$ i $R>0$ ve donada per la densitat

$$f(x)={\begin{cases}{\dfrac {2}{\pi R^{2}}}{\sqrt {R^{2}-(x-a)^{2}}},&{\text{si}}\ x\in [a-R,a+R],\\\\0,&{\text{en cas contrari.}}\\\end{cases}}$$

La definició inicial que hem donat a la primera secció correspon als paràmetres $a=0$ i $R=2$. Designarem aquesta distribució per $w_{a,R}$. Quan $a=0$, llavors la denotarem per $w_{R}$.

Si $X$ és una variable aleatòria amb la distribució del semicercle, $X\sim w_{2}$ , aleshores $(R/2)X+a\sim w_{a,R}$ , la qual cosa permet deduir els moments, la funció generatriu de moments, etc. de la distribució $w_{a,R}$.

Estudiem amb més detall la distribució $w_{R}$: la seva densitat és

$$f(x)={\begin{cases}{\dfrac {2}{\pi R^{2}}}{\sqrt {R^{2}-x^{2}}},&{\text{si}}\ x\in [-R,R],\\\\0,&{\text{en cas contrari.}}\\\end{cases}}$$

Vegeu a la Figura 2 diverses densitats segons el paràmetre $R$.

Figura 2. Funció de densitat de la distribució del semicercle de $w_{R}$ amb paràmetre $R$

Aplicant que si $X\sim w_{2}$ , aleshores $(R/2)X\sim w_{R}$ , deduïm que el moment d'ordre $k$ de $w_{R}$ és

$$m_{k}^{(R)}={\begin{cases}0,&{\text{si}}\ k\ {\text{és senar}},\\\\{\dfrac {R^{2n}}{2^{2n}}}\,C_{n},&{\text{si}}\ k=2n,\end{cases}}$$

En particular, si considerem una variable aleatòria $Y\sim w_{R}$ llavors,

$$E[Y]=0\quad {\text{i}}\quad {\text{Var}}(Y)={\frac {R^{2}}{2}}.$$

La funció generatriu de moments és

$$M_{R}(t)={\frac {2}{Rt}}I_{1}(Rt).$$

La funció característica val

$$\varphi _{R}(t)={\frac {2}{Rt}}J_{1}(Rt).$$

El resultat de Wigner

Ens limitarem al cas de matrius aleatòries reals; per a resultats més generals, vegeu, per exemple, Anderson et al. . Considerem un espai de probabilitat $(\Omega ,{\cal {A,P)}}$. Una matriu aleatòria (real) és una matriu

$${\boldsymbol {M}}={\begin{pmatrix}X_{11}&X_{12}&\cdots &X_{1n}\\X_{21}&X_{22}&\cdots &X_{2n}\\\vdots &\vdots &&\vdots \\X_{n1}&X_{n2}&\cdots &X_{nn}\\\end{pmatrix}}$$

on cada component $X_{ij}$ és una variable aleatòria ordinària.
Considerarem només matrius simètriques i suposarem certes condicions d'independència i sobre els moments de les variables. Concretament, suposarem:

1. Les variables de la diagonal $X_{11},\dots ,X_{nn}$ són independents i tenen la mateixa distribució (iid), amb variància finita.

2. Les variables $X_{ij},\ i<j$ són independents i tenen la mateixa distribució (iid), amb esperança $E[X_{ij}]=0$ i variància ${\text{Var}}(X_{ij})=1$.

3. Les variables del triangle inferior esquerra s'obtenen per simetria respecte de la diagonal: $X_{ji}=X_{ij},\ j>i$, amb la qual cosa la matriu serà simètrica.

Així, la matriu serà de la forma

$${\boldsymbol {M}}={\begin{pmatrix}X_{11}&X_{12}&\cdots &X_{1n}\\X_{12}&X_{22}&\cdots &X_{2n}\\\vdots &\vdots &&\vdots \\X_{1n}&X_{2n}&\cdots &X_{nn}\\\end{pmatrix}}$$

Les matrius amb aquestes característiques o similars (depenent dels autors, vegeu les referències citades) s'anomenen matrius de Wigner.

Atès que la matriu és simètrica, els seus valors propis seran tots reals; els designarem per $\lambda _{1},\dots ,\lambda _{n}$, i degut al caràcter aleatori de la matriu, seran variables aleatòries. Per a $A\subset \mathbb {R}$ (en rigor, $A$ és un subconjunt de Borel de $\mathbb {R}$) definim

$$\mu (A)={\frac {{\text{Card}}\{j:\ \lambda _{j}\in A\}}{n}},$$

on ${\text{Card}}(C)$ és el cardinal d'un conjunt $C$. Es tracta d'una probabilitat aleatòria, ja que per a cada realització de l'experiment aleatori, $\omega \in \Omega$,

$$\mu (A)(\omega )={\frac {{\text{Card}}\{j:\ \lambda _{j}(\omega )\in A\}}{n}},$$

serà una probabilitat ordinària sobre $\mathbb {R}$. L funció $\mu$ s'anomena mesura espectral empírica de la matriu ${\boldsymbol {M}}$.

Teorema. Considerem una successió de matrius aleatòries , $\{{\boldsymbol {M}}_{n},\ n\geq 1\}$, ${\boldsymbol {M}}_{n}$ de dimensions $n\times n$, amb les propietats 1,2 i 3 que hem enunciat abans. Definim $${\boldsymbol {W}}_{n}={\frac {1}{\sqrt {n}}}\,{\boldsymbol {M}}_{n}$$ i designem per $\mu _{n}$ la mesura espectral empírica de ${\boldsymbol {W}}_{n}$. Aleshores, amb probabilitat 1, $$\lim _{n\to \infty }\mu _{n}=S,\quad {\text{en distribució}},$$ on $S$ és la distribució del semicercle.

Figura 3. Il·Iustració del teorema de Wigner. En blau histograma dels valors propis d'una matriu $M_{1000}$ i en vermell la densitat de la distribució de Wigner.

Per fer una il·lustració empírica d'aquest teorema, s'han calculat els valors propis d'una matriu ${\boldsymbol {M}}_{n}$ amb $n=1000$, amb totes les variables amb distribució normal estàndard ${\mathcal {N}}(0,1)$. A la Figura 3 hi ha l'histograma dels 1000 valors propis (en blau), on s'ha superposat el gràfic de la funció de densitat de la llei del semicercle (en vermell).

Figura 3. Les variables aleatòries de la part triangular dreta són independents

que les variables aleatòries de la part triangular dreta són independents, vegeu la Figura 3. Per tenir la simetria de la matriu prendrem $X_{ji}=X_{ij},\ j>i$. Així, la matriu serà de la forma

$${\boldsymbol {M}}={\begin{pmatrix}X_{11}&X_{12}&\cdots &X_{1n}\\X_{12}&X_{22}&\cdots &X_{2n}\\\vdots &\vdots &&\vdots \\X_{1n}&X_{2n}&\cdots &X_{nn}\\\end{pmatrix}}$$
A més suposarem certes condicions sobre les distribucions de les variables. Concretament, suposarem:

• Les variables de la diagonal $X_{11},\dots ,X_{nn}$ són independents i amb la mateixa llei, amb variància finita.
• Les variables $X_{ij},\ i<j$ són independents i tenen la mateixa llei, amb esperança $E[X_{ij}]=0$ i variància ${\text{Var}}(X_{ij})=1$ .
• Les variables del triangle inferior esquerra s'obtenen per simetria respecte de la diagonal: $X_{ji}=X_{ij},\ j>i$, amb la qual cosa la matriu és simètrica.
• Les variables $X_{ij},\ 1\leq i\leq j\leq n$, són independents.
• Les variables de la diagonal $X_{11},\dots ,X_{nn}$ tenen la mateixa llei, amb variància finita.
• Les variables $X_{ij},\ i\neq j$ tenen la mateixa llei, amb esperança $E[X_{ij}]=0$ i variància ${\text{Var}}(X_{ij})=1$ .

Bibliografia

• Anderson, Greg W.; Guionnet, Alice; Zeituni, Ofer. An introduction to random matrices. Cambridge: Cambridge University Press, 2010. ISBN 978-0-521-19452-5. 

• Bai, Zhidong; Silverstein, Jack W.. Spectral analysis of large dimensional random matrices. 2nd ed. New York: Springer, 2010. ISBN 978-1-4419-0661-8. 

• Hiai, Fumio; Petz, Dénes. The semicircle law, free random variables, and entropy. Providence, RI: American Mathematical Society, 2000. ISBN 0-8218-2081-8. 

• Mingo, James A.; Speicher, Roland. Free probability and random matrices, 2017. ISBN 978-1-4939-6942-5. 

• Tao, Terence. Topics in random matrix theory. Providence, R.I.: American Mathematical Society, 2012. ISBN 978-0-8218-7430-1. 

• Voiculescu, D. V.; Dykema, K.J.; Nica, A. Free random variables : a noncommutative probability approach to free products with applications to random matrices, operator algebras, and harmonic analysis on free groups. Providence, R.I., USA: American Mathematical Society, 1992. ISBN 0-8218-6999-X. 

• Wigner, Eugene P. «Characteristic Vectors of Bordered Matrices With Infinite Dimensions». The Annals of Mathematics, 62, 3, 1955, pàg. 548-564. DOI: 10.2307/1970079.

• Wigner, Eugene P. «On the Distribution of the Roots of Certain Symmetric Matrices». Annals of Mathematics, 67, 2, 1958, pàg. 325–327. DOI: 10.2307/1970008. ISSN: 0003-486X.