Articulo de referencia

Distribución normal compleja

\\mathbf{\\mu} \\in \\mathbb{C}^n — [[location parameter|location]] \n \\Gamma \\in \\mathbb{C}^{n \\times n} — [[Complex covariance matrix|covariance matrix]] ([[positive semi-...

En teoría de la probabilidad , la familia de distribuciones normales complejas , denotadadonorte{\displaystyle {\mathcal {CN}}}onortedo{\displaystyle {\mathcal {N}}_{\mathcal {C}}}, caracteriza variables aleatorias complejas cuyas partes real e imaginaria son conjuntamente normales . [ 1 ] La familia normal compleja tiene tres parámetros: parámetro de localización μ , matriz de covarianzaΓ{\displaystyle \Gamma }y la matriz de relacióndo{\displaystyle C}. La distribución normal compleja estándar es la distribución univariada conμ=0{\displaystyle \mu =0},Γ=1{\displaystyle \Gamma =1}, ydo=0{\displaystyle C=0}.

Una subclase importante de la familia normal compleja se denomina normal compleja circularmente simétrica (central) y corresponde al caso de matriz de relación cero y media cero:μ=0{\displaystyle \mu =0} ydo=0{\displaystyle C=0}. [ 2 ] Este caso se utiliza ampliamente en el procesamiento de señales , donde a veces se le denomina simplemente normal complejo en la literatura.

Definiciones

variable aleatoria normal estándar compleja

La variable aleatoria normal compleja estándar o variable aleatoria gaussiana compleja estándar es una variable aleatoria compleja.Z{\displaystyle Z}cuyas partes real e imaginaria son variables aleatorias independientes distribuidas normalmente con media cero y varianza1/2{\displaystyle 1/2}. [ 3 ] : pág. 494 [ 4 ] : págs. 501 Formalmente,

dónde{\displaystyle \perp \!\!\!\perp }denota independencia yZdonorte(0,1){\displaystyle Z\sim {\mathcal {CN}}(0,1)}indica queZ{\displaystyle Z}es una variable aleatoria normal compleja estándar.

variable aleatoria normal compleja

Suponerincógnita{\displaystyle X}yY{\displaystyle Y}son variables aleatorias reales tales que(incógnita,Y)T{\displaystyle (X,Y)^{\mathrm {T} }}es un vector aleatorio normal bidimensional . Entonces, la variable aleatoria complejaZ=incógnita+iY{\displaystyle Z=X+iY}Se denomina variable aleatoria normal compleja o variable aleatoria gaussiana compleja . [ 3 ] : pág. 500

vector aleatorio normal estándar complejo

Un vector aleatorio complejo n-dimensionalZ=(Z1,,Znorte)T{\displaystyle \mathbf {Z} =(Z_{1},\ldots,Z_{n})^{\mathrm {T} }}es un vector aleatorio normal estándar complejo o un vector aleatorio gaussiano estándar complejo si sus componentes son independientes y todas ellas son variables aleatorias normales complejas estándar como se definió anteriormente. [ 3 ] : pág. 502 [ 4 ] : págs. 501 QueZ{\displaystyle \mathbf {Z} }es un vector aleatorio normal complejo estándar se denotaZdonorte(0,Inorte){\displaystyle \mathbf {Z} \sim {\mathcal {CN}}(0,{\boldsymbol {I}}_{n})}.

vector aleatorio normal complejo

Siincógnita=(incógnita1,,incógnitanorte)T{\displaystyle \mathbf {X} =(X_{1},\ldots ,X_{n})^{\mathrm {T} }}yY=(Y1,,Ynorte)T{\displaystyle \mathbf {Y} =(Y_{1},\ldots ,Y_{n})^{\mathrm {T} }}son vectores aleatorios enRnorte{\displaystyle \mathbb {R} ^{n}}de tal manera que[incógnita,Y]{\displaystyle [\mathbf {X} ,\mathbf {Y} ]}es un vector aleatorio normal con2norte{\displaystyle 2n}componentes. Entonces decimos que el vector aleatorio complejo

Z=incógnita+iY{\displaystyle \mathbf {Z} =\mathbf {X} +i\mathbf {Y} \,}

es un vector aleatorio normal complejo o un vector aleatorio gaussiano complejo .

Media, covarianza y relación

La distribución gaussiana compleja se puede describir con 3 parámetros: [ 5 ]

μ=mi[Z],Γ=mi[(Zμ)(Zμ)H],do=mi[(Zμ)(Zμ)T],{\displaystyle \mu =\operatorname {E} [\mathbf {Z} ],\quad \Gamma =\operatorname {E} [(\mathbf {Z} -\mu )({\mathbf {Z} }-\mu )^{\mathrm {H} }],\quad C=\operatorname {E} [(\mathbf {Z} -\mu )(\mathbf {Z} -\mu )^{\mathrm {T} }],}

dóndeZT{\displaystyle \mathbf {Z} ^{\mathrm {T} }}denota la transpuesta de la matriz deZ{\displaystyle \mathbf {Z} }, yZH{\displaystyle \mathbf {Z} ^{\mathrm {H} }}denota transpuesta conjugada . [ 3 ] : pág. 504 [ 4 ] : págs. 500

Aquí el parámetro de ubicaciónμ{\displaystyle \mu }es un vector complejo n-dimensional; la matriz de covarianzaΓ{\displaystyle \Gamma }es hermitiana y semidefinida positiva ; y, la matriz de relación o matriz de pseudocovarianzado{\displaystyle C}es simétrico . El vector aleatorio normal complejoZ{\displaystyle \mathbf {Z} }ahora se puede denotar comoZ  donorte(μ, Γ, do).{\displaystyle \mathbf {Z} \ \sim \ {\mathcal {CN}}(\mu,\ \Gamma,\ C).}Además, matricesΓ{\displaystyle \Gamma }ydo{\displaystyle C}son tales que la matriz

PAG=Γ¯doHΓ1do{\displaystyle P={\overline {\Gamma }}-{C}^{\mathrm {H} }\Gamma ^{-1}C}

también es definida no negativa dondeΓ¯{\displaystyle {\overline {\Gamma }}}denota el conjugado complejo deΓ{\displaystyle \Gamma }. [ 5 ]

Relaciones entre matrices de covarianza

Como ocurre con cualquier vector aleatorio complejo, las matricesΓ{\displaystyle \Gamma }ydo{\displaystyle C}puede estar relacionado con las matrices de covarianza deincógnita=(Z){\displaystyle \mathbf {X} =\Re (\mathbf {Z} )}yY=(Z){\displaystyle \mathbf {Y} =\Soy (\mathbf {Z} )}a través de expresiones

Vincógnitaincógnitami[(incógnitaμincógnita)(incógnitaμincógnita)T]=12Re[Γ+do],VincógnitaYmi[(incógnitaμincógnita)(YμY)T]=12Soy[Γ+do],VYincógnitami[(YμY)(incógnitaμincógnita)T]=12Soy[Γ+do],VYYmi[(YμY)(YμY)T]=12Re[Γdo],{\displaystyle {\begin{aligned}&V_{XX}\equiv \operatorname {E} [(\mathbf {X} -\mu _{X})(\mathbf {X} -\mu _{X})^{\mathrm {T} }]={\tfrac {1}{2}}\operatorname {Re} [\Gamma +C],\quad V_{XY}\equiv \operatorname {E} [(\mathbf {X} -\mu _{X})(\mathbf {Y} -\mu _{Y})^{\mathrm {T} }]={\tfrac {1}{2}}\operatorname {Im} [-\Gamma +C],\\&V_{YX}\equiv \operatorname {E} [(\mathbf {Y} -\mu _{Y})(\mathbf {X} -\mu _{X})^{\mathrm {T} }]={\tfrac {1}{2}}\operatorname {Im} [\Gamma +C],\quad \,V_{YY}\equiv \operatorname {E} [(\mathbf {Y} -\mu _{Y})(\mathbf {Y} -\mu _{Y})^{\mathrm {T} }]={\tfrac {1}{2}}\operatorname {Re} [\Gamma -C],\end{aligned}}}

y viceversa

Γ=Vincógnitaincógnita+VYY+i(VYincógnitaVincógnitaY),do=VincógnitaincógnitaVYY+i(VYincógnita+VincógnitaY).{\displaystyle {\begin{aligned}&\Gamma =V_{XX}+V_{YY}+i(V_{YX}-V_{XY}),\\&C=V_{XX}-V_{YY}+i(V_{YX}+V_{XY}).\end{aligned}}}

Función de densidad

La función de densidad de probabilidad para la distribución normal compleja se puede calcular como

F(z)=1πnortedet(Γ)det(PAG)exp{12[zμz¯μ¯]H[Γdodo¯Γ¯]1[zμz¯μ¯]}=det(PAG1¯RPAG1R)det(PAG1)πnortemi(zμ)PAG1¯(zμ)+Re((zμ)RPAG1¯(zμ)),{\displaystyle {\begin{aligned}f(z)&={\frac {1}{\pi ^{n}{\sqrt {\det(\Gamma )\det(P)}}}}\,\exp \!\left\{-{\frac {1}{2}}{\begin{bmatrix}z-\mu \\{\overline {z}}-{\overline {\mu }}\end{bmatrix}}^{\mathrm {H} }{\begin{bmatrix}\Gamma &C\\{\overline {C}}&{\overline {\Gamma }}\end{bmatrix}}^{\!\!-1}\!{\begin{bmatrix}z-\mu \\{\overline {z}}-{\overline {\mu }}\end{bmatrix}}\right\}\\[8pt]&={\tfrac {\sqrt {\det \left({\overline {P^{-1}}}-R^{\ast }P^{-1}R\right)\det(P^{-1})}}{\pi ^{n}}}\,e^{-(z-\mu )^{\ast }{\overline {P^{-1}}}(z-\mu )+\operatorname {Re} \left((z-\mu )^{\intercal }R^{\intercal }{\overline {P^{-1}}}(z-\mu )\right)},\end{aligned}}}

dóndeR=doHΓ1{\displaystyle R=C^{\mathrm {H} }\Gamma ^{-1}}yPAG=Γ¯Rdo{\displaystyle P={\overline {\Gamma }}-RC}.

Función característica

La función característica de la distribución normal compleja viene dada por [ 5 ].

φ(w)=exp{iRe(w¯μ)14(w¯Γw+Re(w¯dow¯))},{\displaystyle \varphi (w)=\exp \!{\big \{}i\operatorname {Re} ({\overline {w}}'\mu )-{\tfrac {1}{4}}{\big (}{\overline {w}}'\Gamma w+\operatorname {Re} ({\overline {w}}'C{\overline {w}}){\big )}{\big \}},}

donde el argumentow{\displaystyle w}es un vector complejo n -dimensional.

Propiedades

  • SiZ{\displaystyle \mathbf {Z} }es un vector normal complejo de n dimensiones,A{\displaystyle {\boldsymbol {A}}}una matriz m×n , yb{\displaystyle b}un vector m constante , luego la transformación linealAZ+b{\displaystyle {\boldsymbol {A}}\mathbf {Z} +b}también se distribuirá de forma compleja-normal:
Z  donorte(μ,Γ,do)AZ+b  donorte(Aμ+b,AΓAH,AdoAT){\displaystyle Z\ \sim \ {\mathcal {CN}}(\mu ,\,\Gamma ,\,C)\quad \Rightarrow \quad AZ+b\ \sim \ {\mathcal {CN}}(A\mu +b,\,A\Gamma A^{\mathrm {H} },\,ACA^{\mathrm {T} })}
  • SiZ{\displaystyle \mathbf {Z} }es un vector normal complejo de n dimensiones , entonces
2[(Zμ)HPAG1¯(Zμ)Re((Zμ)TRTPAG1¯(Zμ))]  χ2(2norte){\displaystyle 2{\Big [}(\mathbf {Z} -\mu )^{\mathrm {H} }{\overline {P^{-1}}}(\mathbf {Z} -\mu )-\operatorname {Re} {\big (}(\mathbf {Z} -\mu )^{\mathrm {T} }R^{\mathrm {T} }{\overline {P^{-1}}}(\mathbf {Z} -\mu ){\big )}{\Big ]}\ \sim \ \chi ^{2}(2n)}
  • Teorema del límite central . SiZ1,,ZT{\displaystyle Z_{1},\ldots ,Z_{T}}son variables aleatorias complejas independientes e idénticamente distribuidas, entonces
T(1Tt=1TZtmi[Zt]) d donorte(0,Γ,do),{\displaystyle {\sqrt {T}}{\Big (}{\tfrac {1}{T}}\textstyle \sum _{t=1}^{T}Z_{t}-\operatorname {E} [Z_{t}]{\Big )}\ {\xrightarrow {d}}\ {\mathcal {CN}}(0,\,\Gamma ,\,C),}
dóndeΓ=mi[ZZH]{\displaystyle \Gamma =\operatorname {E} [ZZ^{\mathrm {H} }]}ydo=mi[ZZT]{\displaystyle C=\operatorname {E} [ZZ^{\mathrm {T} }]}.

Caso central de simetría circular

Definición

Un vector aleatorio complejoZ{\displaystyle \mathbf {Z} }se denomina simétrico circularmente si para cada deterministaφ[π,π){\displaystyle \varphi \in [-\pi ,\pi )}la distribución demiiφZ{\displaystyle e^{\mathrm {i} \varphi }\mathbf {Z} }es igual a la distribución deZ{\displaystyle \mathbf {Z} }[ 4 ] : págs . 500–501

Los vectores aleatorios complejos normales centrales que son circularmente simétricos son de particular interés porque están completamente especificados por la matriz de covarianza.Γ{\displaystyle \Gamma }.

La distribución normal compleja circularmente simétrica (central) corresponde al caso de media cero y matriz de relación cero, es decirμ=0{\displaystyle \mu =0}ydo=0{\displaystyle C=0}. [ 3 ] : pág. 507 [ 7 ] Esto se suele denotar

Zdonorte(0,Γ){\displaystyle \mathbf {Z} \sim {\mathcal {CN}}(0,\,\Gamma )}

Distribución de las partes real e imaginaria

SiZ=incógnita+iY{\displaystyle \mathbf {Z} =\mathbf {X} +i\mathbf {Y} }es circularmente simétrico (central) normal complejo, entonces el vector[incógnita,Y]{\displaystyle [\mathbf {X} ,\mathbf {Y} ]}es una distribución normal multivariada con estructura de covarianza

(incógnitaY)  norte([00], 12[ReΓSoyΓSoyΓReΓ]){\displaystyle {\begin{pmatrix}\mathbf {X} \\\mathbf {Y} \end{pmatrix}}\ \sim \ {\mathcal {N}}{\Big (}{\begin{bmatrix}0\\0\end{bmatrix}},\ {\tfrac {1}{2}}{\begin{bmatrix}\operatorname {Re} \,\Gamma &-\operatorname {Im} \,\Gamma \\\operatorname {Im} \,\Gamma &\operatorname {Re} \,\Gamma \end{bmatrix}}{\Big )}}

dóndeΓ=mi[ZZH]{\displaystyle \Gamma =\operatorname {E} [\mathbf {Z} \mathbf {Z} ^{\mathrm {H} }]}.

Función de densidad de probabilidad

Para una matriz de covarianza no singularΓ{\displaystyle \Gamma }, su distribución también puede simplificarse como [ 3 ] : pág. 508

FZ(z)=1πnortedet(Γ)mi(zμ)HΓ1(zμ){\displaystyle f_{\mathbf {Z} }(\mathbf {z} )={\tfrac {1}{\pi ^{n}\det(\Gamma )}}\,e^{-(\mathbf {z} -\mathbf {\mu } )^{\mathrm {H} }\Gamma ^{-1}(\mathbf {z} -\mathbf {\mu } )}}.

Por lo tanto, si la media no es ceroμ{\displaystyle \mu }y matriz de covarianzaΓ{\displaystyle \Gamma }son desconocidos, una función de verosimilitud logarítmica adecuada para un único vector de observaciónz{\displaystyle z}sería

ln(L(μ,Γ))=ln(det(Γ))(zμ)¯Γ1(zμ)norteln(π).{\displaystyle \ln(L(\mu ,\Gamma ))=-\ln(\det(\Gamma ))-{\overline {(z-\mu )}}'\Gamma ^{-1}(z-\mu )-n\ln(\pi ).}

La distribución normal compleja estándar (definida en la ecuación 1 ) corresponde a la distribución de una variable aleatoria escalar conμ=0{\displaystyle \mu =0},do=0{\displaystyle C=0}yΓ=1{\displaystyle \Gamma =1}Por lo tanto, la distribución normal compleja estándar tiene densidad

FZ(z)=1πmiz¯z=1πmi|z|2.{\displaystyle f_{Z}(z)={\tfrac {1}{\pi }}e^{-{\overline {z}}z}={\tfrac {1}{\pi }}e^{-|z|^{2}}.}

Propiedades

La expresión anterior demuestra por qué el casodo=0{\displaystyle C=0},μ=0{\displaystyle \mu =0}se denomina “circularmente simétrica”. La función de densidad depende únicamente de la magnitud dez{\displaystyle z}pero no en su argumento . Como tal, la magnitud|z|{\displaystyle |z|}de una variable aleatoria normal compleja estándar tendrá la distribución de Rayleigh y la magnitud al cuadrado|z|2{\displaystyle |z|^{2}}tendrá la distribución exponencial , mientras que el argumento se distribuirá uniformemente en[π,π]{\displaystyle [-\pi ,\pi ]}.

Si{Z1,,Zk}{\displaystyle \left\{\mathbf {Z} _{1},\ldots ,\mathbf {Z} _{k}\right\}}son vectores aleatorios normales complejos circulares n -dimensionales independientes e idénticamente distribuidos conμ=0{\displaystyle \mu =0}, entonces la norma cuadrática aleatoria

Q=j=1kZjHZj=j=1kZj2{\displaystyle Q=\sum _{j=1}^{k}\mathbf {Z} _{j}^{\mathrm {H} }\mathbf {Z} _{j}=\sum _{j=1}^{k}\|\mathbf {Z} _{j}\|^{2}}

tiene la distribución chi-cuadrado generalizada y la matriz aleatoria

W=j=1kZjZjH{\displaystyle W=\sum _{j=1}^{k}\mathbf {Z} _{j}\mathbf {Z} _{j}^{\mathrm {H} }}

tiene la distribución compleja de Wishart conk{\displaystyle k}grados de libertad. Esta distribución se puede describir mediante la función de densidad.

F(w)=det(Γ1)kdet(w)knorteπnorte(norte1)/2j=1k(kj)¡ mitr(Γ1w){\displaystyle f(w)={\frac {\det(\Gamma ^{-1})^{k}\det(w)^{k-n}}{\pi ^{n(n-1)/2}\prod _{j=1}^{k}(k-j)!}}\ e^{-\operatorname {tr} (\Gamma ^{-1}w)}}

dóndeknorte{\displaystyle k\geq n}, yw{\displaystyle w}es unnorte×norte{\displaystyle n\times n}matriz semidefinida positiva.

Véase también

Referencias

  1. Goodman, NR (1963). "Análisis estadístico basado en una cierta distribución gaussiana compleja multivariada (una introducción)" . The Annals of Mathematical Statistics . 34 (1): 152– 177. doi : 10.1214/aoms/1177704250 . JSTOR 2991290 . 
  2. capítulo del libro, Gallager.R , pág. 9.
  3. 1 2 3 4 5 6 Lapidoth, A. (2009). Fundamentos de la comunicación digital . Cambridge University Press. ISBN 9780521193955.
  4. 1 2 3 4 Tse, David (2005). Fundamentos de la comunicación inalámbrica . Cambridge University Press. ISBN 9781139444668.
  5. 1 2 3 Picinbono, Bernard (1996). "Vectores aleatorios complejos de segundo orden y distribuciones normales" . IEEE Transactions on Signal Processing . 44 (10): 2637– ​​2640. Bibcode : 1996ITSP...44.2637P . doi : 10.1109/78.539051 .
  6. Daniel Wollschlaeger. "La distribución de Hoyt (Documentación para el paquete de R 'shotGroups' versión 0.6.2)" .
  7. capítulo del libro, Gallager.R