Existe uma expressão canônica ou analítica para a distribuição de probabilidade da variável aleatória complexa simétrica circular :
que ?Zθ∼U(0,2π)
Notas laterais:
Sabe-se que as partes reais e imaginárias, ou seja:
têm densidades marginais dadas por :
mas porque não são independentes, calculando seu PDF conjunto não é trivial.ℜ(Z)=cosθℑ(Z)=sinθfℜ(Z)(z)=fℑ(Z)(z)=1π√1−z2,−1<z<1,
EDIT: é diferente de um complexo normal, aqui, a amplitudeé determinístico e identicamente 1, enquanto que se fosse complexo normal,seria Rayleigh distribuído.Z
complex-random-variable
Robert L.
fonte
fonte
Respostas:
Como as partes reais e imaginárias são muito dependentes uma da outra (se você tem o valor de uma, você sabe exatamente o valor da outra), parece que você poderia aplicar o pdf marginal da parte real , dado um valor da parte imaginária :rr ii
fri(r,i)=fr|i(r | i)fi(i)fri(r,i)=fr|i(r | i)fi(i)
Você anotou o pdf das partes reais e imaginárias individualmente:
fr(z)=fi(z)=1π√1−z2fr(z)=fi(z)=1π1−z2−−−−−√
Isso deixa o pdf marginal . Lembre-se de que, para uma determinada realização da variável aleatória , os dois componentes estão relacionados de forma determinística:fr|i(r | i)fr|i(r | i) ZZ
r2+i2=cos2(θ)+sin2(θ)=1r2+i2=cos2(θ)+sin2(θ)=1
Dada essa relação, podemos resolver em termos de :rr ii
r2=1−i2r2=1−i2
r=±√1−i2r=±1−i2−−−−−√
Portanto, o pdf marginal de dado um valor de é um par de impulsos:rr ii
fr|i(r | i)=12δ(r−√1−i2)+12δ(r+√1−i2)fr|i(r | i)=12δ(r−1−i2−−−−−√)+12δ(r+1−i2−−−−−√)
Colocá-los juntos produziria:
fri(r,i)=δ(r−√1−i2)+δ(r+√1−i2)2π√1−i2fri(r,i)=δ(r−1−i2−−−−−√)+δ(r+1−i2−−−−−√)2π1−i2−−−−−√
Pensando nisso geometricamente, para cada linha horizontal (para ) no plano , existem apenas dois pontos que são diferentes de zero , e o pdf tem uma altura infinita nesses pontos. Como poderíamos esperar, esses pontos de interseção (ou seja, pontos em que o pdf é diferente de zero) são onde a linha horizontal se cruza com o círculo unitário!i=i0i=i0 i0∈[−1,1]i0∈[−1,1] riri r0=±√1−i20r0=±1−i20−−−−−√
Isso significa que o pdf da junta é de valor zero, exceto ao longo do círculo unitário, onde assume altura infinita. Isso está alinhado com a intuição, pois a definição da variável aleatória garante que eles possam assumir apenas valores que estão no círculo unitário.ZZ
Não há nada de especial na maneira específica como descrevi isso; você também pode transpor o problema e observar as linhas verticais no plano da forma e encontrará a mesma relação devido ao acoplamento próximo das duas variáveis aleatórias.riri r=r0r=r0
Acredito que essa formulação seja equivalente à da resposta de AlexTP , mas sua derivação é provavelmente mais intuitiva.
fonte
Evitar cálculos complicados, deixar e ser iid normais padrão variáveis aleatórias, a variável aleatória tem a mesma distribuição de (fácil de ver e o ângulo de é equivalente ao ângulo de um normal circularmente simétrico, portanto uniforme) .XX YY ZZ VV
V≜(X√X2+Y2,Y√X2+Y2)V≜(XX2+Y2−−−−−−−√,YX2+Y2−−−−−−−√) ‖V‖=1∥V∥=1 VV
Esse tipo de é uma das construções de um ponto uniformemente distribuído no círculo (que pode ser generalizado para a esfera , consulte Seleção de ponto de esfera e, por exemplo, esta resposta ).VV (n−1)(n−1)
Assim, o PDF de é simplesmente o inverso da circunferência do círculo unitário. Para com fixo e uniforme ,ZZ Zρ=ρejΘZρ=ρejΘ ρρ ΘΘ
em coordenadas polares (onde a área infinitesimal é ), rdrdθrdrdθ fR,Θ(r,θ)=12πδ(r−ρ)fR,Θ(r,θ)=12πδ(r−ρ)
fonte
Com base nas respostas existentes, que abriram meus olhos para o que está acontecendo aqui, gostaria de apresentar mais uma expressão muito simples para a solução, que é apenas ligeiramente diferente daquela na resposta da AlexTP (e que acabou sendo equivalente a o indicado na resposta de Jason R , como mostrado abaixo na parte EDIT).
[EDIT: agora que AlexTP editou sua resposta, nossas expressões para o PDF são idênticas; então todas as três respostas finalmente concordam uma com a outra].
Seja a variável aleatória complexa definida comoZ=X+jYZ=X+jY
Z=ρejθZ=ρejθ(1)
onde o raio é determinístico e dado, enquanto o ângulo é aleatório e uniformemente distribuído em . Afirmo, sem mais provas, que é circular simétrico, do qual se conclui que sua função de densidade de probabilidade (PDF) deve satisfazerρρ θθ [0,2π)[0,2π) ZZ
fZ(z)=fZ(x+jy)=fZ(r),withr=√x2+y2fZ(z)=fZ(x+jy)=fZ(r),withr=x2+y2−−−−−−√(2)
isto é, pode ser escrito como uma função do raio (magnitude) .rr
Como o PDF deve ser zero em qualquer lugar, exceto , e como deve ser integrado à unidade (quando integrado no plano bidimensional), o único PDF possível ér=ρr=ρ
Pode ser mostrado que conduz às densidades marginais correcção para as variáveis aleatórias e .(3)(3) XX YY
EDITAR:
Após uma discussão muito útil nos comentários, parece que conseguimos chegar a um acordo sobre uma solução para o problema. Vou mostrar a seguir que a fórmula despretensioso é realmente equivalente a fórmula olhando mais envolvido na resposta de Jason R . Note que eu uso para a magnitude (raio) do complexo RV , enquanto que na resposta de Jason denota a parte real . Usarei e para as partes reais e imaginárias, respectivamente. Aqui vamos nós:(3)(3) rr ZZ rr ZZ xx yy
fZ(r)=12πδ(r-ρ)=12πδ(√x2+y2-ρ)fZ(r)=12πδ(r−ρ)=12πδ(x2+y2−−−−−−√−ρ)(4)
Nós sabemos isso δ(g(x))δ(g(x)) É dado por
δ(g(x))=∑Euδ(x-xEu)|g′(xEu)|δ(g(x))=∑iδ(x−xi)|g′(xi)|(5)
Onde xEuxi são as raízes (simples) de g(x)g(x) . Nós temos
g(x)=√x2+y2-ρeg′(x)=x√x2+y2=xrg(x)=x2+y2−−−−−−√−ρandg′(x)=xx2+y2−−−−−−√=xr(6)
As duas raízes xEuxi estamos
x1,2=±√ρ2-y2x1,2=±ρ2−y2−−−−−−√(7)
Consequentemente,
|g′(x1)|=|g′(x2)|=√ρ2-y2ρ=√1-(yρ)2|g′(x1)|=|g′(x2)|=ρ2−y2−−−−−−√ρ=1−(yρ)2−−−−−−−−√(8)
Com (5)(5) -(8)(8) , Eq. (4)(4) pode ser escrito como
fX,Y(x,y)=12π√1-(yρ)2[δ(x-√ρ2-y2)+δ(x+√ρ2-y2)]fX,Y(x,y)=12π1−(yρ)2−−−−−−−−√[δ(x−ρ2−y2−−−−−−√)+δ(x+ρ2−y2−−−−−−√)](9)
Para ρ=1ρ=1 , Eq. (9)(9) é idêntica à expressão dada na resposta de Jason R .
Acho que agora podemos concordar que a Eq. (3)(3) é uma expressão correta (e muito simples) para o PDF do complexo RV Z=ρejθZ=ρejθ com determinista ρρ e distribuído uniformemente θθ .
fonte