Distribución de

7

¿Existe una expresión canónica o analítica para la distribución de probabilidad de la variable aleatoria compleja circularmente simétrica Z:

Z=ejθ,
dónde θU(0,2π)?

Notas al margen:

Se sabe que las partes reales e imaginarias, es decir:

(Z)=cosθ(Z)=sinθ
tienen densidades marginales dadas por :
f(Z)(z)=f(Z)(z)=1π1z2,1<z<1,
pero como no son independientes, calcular su PDF conjunto no es trivial.

EDITAR: Z es diferente de una normal compleja en que aquí, la amplitud |Z|es determinista e idénticamente 1, mientras que si Z eran complejos normales, |Z| Sería distribuido por Rayleigh.

Robert L.
fuente
¿Cómo es esto diferente del complejo circular simétrico normal ?
Maxtron
1
@Maxtron la amplitud aquí es 1, mientras que la normal compleja tiene una amplitud que en sí misma es una variable aleatoria con una distribución de Rayleigh.
Robert L.
@OlliNiemitalo tenía una respuesta para algo similar a lo que recuerdo ...
Fat32
@ Fat32 Sí, esto suena una campana. La otra pregunta es: ¿Cuál es su distribución?
Olli Niemitalo el
@OlliNiemitalo ¡Sí, esa era la pregunta!
Fat32

Respuestas:

6

Dado que las partes real e imaginaria dependen mucho una de la otra (si tiene el valor de una, conoce exactamente el valor de la otra), parece que podría aplicar el pdf marginal de la parte real r, dado un valor de la parte imaginaria i:

fri(r,i)=fr|i(r | i)fi(i)

Usted notó el pdf de las partes reales e imaginarias individualmente:

fr(z)=fi(z)=1π1z2

Eso deja el pdf marginal fr|i(r | i). Recuerde que para una realización dada de la variable aleatoriaZ, los dos componentes están relacionados determinísticamente:

r2+i2=cos2(θ)+sin2(θ)=1

Dada esta relación, podemos resolver por r en términos de i:

r2=1i2
r=±1i2

Por lo tanto, el pdf marginal de r dado un valor de i es un par de impulsos:

fr|i(r | i)=12δ(r1i2)+12δ(r+1i2)

Poner estos juntos produciría:

fri(r,i)=δ(r1i2)+δ(r+1i2)2π1i2

Pensando en esto geométricamente, para cada línea horizontal i=i0 (para i0[1,1]) en el ri avión, solo hay dos puntos r0=±1i02que no son cero, y el pdf tiene una altura infinita en esos puntos. ¡Como podríamos esperar, esos puntos de intersección (es decir, puntos donde el pdf no es cero) están donde la línea horizontal se cruza con el círculo unitario!

Esto significa que el pdf conjunto tiene un valor cero, excepto a lo largo del círculo unitario, donde toma una altura infinita. Eso se alinea con la intuición, como la definición de la variable aleatoriaZ asegura que solo pueda tomar valores que están en el círculo unitario.

No hay nada especial en la forma específica en que lo expuse; También podría transponer el problema y mirar las líneas verticales en elri plano de la forma r=r0 y encontraría la misma relación debido al acoplamiento estrecho de las dos variables aleatorias.

Creo que esta formulación es equivalente a la de la respuesta de AlexTP , pero su derivación es probablemente más intuitiva.

Jason R
fuente
1
¿hay un error tipográfico en r=1i2; posiblemente te refieresr2=1i2? También podemos mostrar que el pdf conjunto se integra a 1 en el plano xy (o eqv ri). Nuevamente también el pdf conjuntofri(r,i)no es circularmente simétrico? (el denominador?) ¿No podemos simplemente usar la lógica que daz=ejθ=cos(θ)+jsin(θ)=x+jy con θ uniforme en [0,2π] y la restricción x2+y2=1 simplemente ofrece una junta circular simétrica pdf del tipo fxy(x,y)=Kδ2(x2+y21) o fR,θ=Kδ2(R1); ¿Un impulso de anillo en el plano xy? donde K es posiblemente1/2π?
Fat32
1
Corregí los errores tipográficos a los que te referías, gracias.
Jason R
@ Fat32: Hice algunas correcciones que me llevaron a lo que suena a lo que te referías en tu comentario anterior. Sin embargo, la respuesta de AlexTP es probablemente más intuitiva.
Jason R
1
Hay un error tipográfico en la copia fr|i(r|i) a fr,i(r,i). Y observo quefr,i(0,1)= es diferente de fr,i(1,0)=12π. Creo que deberían tener el mismo valor, ¿no?
AlexTP
1
todos, tal vez una simulación rápida sería útil mostrando sus puntos :)
AlexTP
9

Evite cálculos complicados, deje X y Yser iid variables aleatorias normales estándar , su variable aleatoriaZ tiene la misma distribución de V

V(XX2+Y2,YX2+Y2)
(fácil de ver V=1 y el ángulo de V es equivalente al ángulo de una Normal circular simétrica, por lo tanto uniforme)

Este tipo de V es una de las construcciones de un punto uniformemente distribuido en círculo (que puede generalizarse para (n1)-sphere, ver Esfera Punto de Recogida y por ejemplo esta respuesta ).

Así, el PDF de Zes simplemente el recíproco de la circunferencia del círculo unitario. porZρ=ρejΘ con fijo ρ y uniforme Θ,

en coordenadas polares (donde el área infinitesimal es rdrdθ),

fR,Θ(r,θ)=12πδ(rρ)

AlexTP
fuente
2
i think this is correct.
robert bristow-johnson
The intuition seems to be so simple behind the proof. All you need to do is pick points from the surface of a unit sphere :)
Maxtron
Can't we just remove the ||z|| in the denominator? Otherwise, how would you evaluate the expression for ||z||0? Of course, we would like fZ(0)=0 to hold.
Matt L.
@MattL. you are right, it is wrong having z in the denominator, not due to the evaluation at limits but because it is simply wrong. I have fixed the final result.
AlexTP
1
@AlexTP: I still think that there shouldn't be any variable in the denominator. Have a look at my answer and please tell me if/where I go wrong.
Matt L.
4

Based on the existing answers, which opened my eyes for what's going on here, I would like to present yet another very simple expression for the solution, which is just slightly different from the one in AlexTP's answer (and which turned out to be equivalent to the one given in Jason R's answer, as shown below in the EDIT-part).

[EDIT: now that AlexTP has edited his answer, our expressions for the PDF are identical; so all three answers finally agree with each other].

Let the complex random variable Z=X+jY be defined as

(1)Z=ρejθ

where the radius ρ is deterministic and given, whereas the angle θ is random and uniformly distributed on [0,2π). I state without further proof that Z is circularly symmetrical, from which it follows that its probability density function (PDF) must satisfy

(2)fZ(z)=fZ(x+jy)=fZ(r),withr=x2+y2

i.e., it can be written as a function of the radius (magnitude) r.

Since the PDF must be zero everywhere except for r=ρ, and since it must integrate to unity (when integrated over the 2-dimensional plane), the only possible PDF is

(3)fZ(r)=12πδ(rρ)

It can be shown that (3) leads to the correct marginal densities for the random variables X and Y.


EDIT:

After some very useful discussion in the comments it appears that we've managed to agree on one solution to the problem. I will show in the following that the unassuming formula (3) is actually equivalent to the more involved looking formula in Jason R's answer. Note that I use r for the magnitude (radius) of the complex RV Z, whereas in Jason's answer r denotes the real part of Z. I will use x and y for the real and imaginary parts, respectively. Here we go:

(4)fZ(r)=12πδ(rρ)=12πδ(x2+y2ρ)

We know that δ(g(x)) is given by

(5)δ(g(x))=iδ(xxi)|g(xi)|

where xi are the (simple) roots of g(x). We have

(6)g(x)=x2+y2ρandg(x)=xx2+y2=xr

The two roots xi are

(7)x1,2=±ρ2y2

Consequently,

(8)|g(x1)|=|g(x2)|=ρ2y2ρ=1(yρ)2

With (5)-(8), Eq. (4) can be written as

(9)fX,Y(x,y)=12π1(yρ)2[δ(xρ2y2)+δ(x+ρ2y2)]

For ρ=1, Eq. (9) is identical to the expression given in Jason R's answer.

I think we can now agree that Eq. (3) is a correct (and very simple) expression for the PDF of the complex RV Z=ρejθ with deterministic ρ and uniformly distributed θ.

Matt L.
fuente
I think the determinant of Jacobian transform should be r then
θ=02πr=0+fZ(r)rdrdθ1
AlexTP
@AlexTP: Take as an extreme example ρ=0, which makes the RV Z deterministic, it's always zero. It's PDF should then be a Dirac impulse at r=z=y=0, which is what I get with the formula I suggested. With the formula in your answer we get an undefined term in that case.
Matt L.
it's true, but I think the value of PDF can be infinity and only the CDF should be bounded. In other words, the Dirac impulse matters only inside integrals.
AlexTP
Yes, but if you integrate your pdf (for ρ=0), what do you get (and how would you do it)?
Matt L.
I would do like this
limρ0+θr=0+12πρδ(rρ)rdrdθ=limρ0+θ12πρρdθ=θ12πdθ=1
The integral is finite while the PDF is not thanks to the Jacobian determinant. My point is that you need that Jacobian determinant r whenever you integrate the PDF in polar coordinate. Tell me if you find it unconvinced.
AlexTP