Distribuição de


7

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 ?ZZZ=ejθ,

Z=ejθ,
θU(0,2π)θ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θ

R(Z)=cosθI(Z)=sinθ
f(Z)(z)=f(Z)(z)=1π1z2,1<z<1,
fR(Z)(z)=fI(Z)(z)=1π1z2,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.ZZ|Z||Z| ZZ|Z||Z|


Como isso é diferente do complexo circular-simétrico normal ?
Maxtron

1
@ Maxtron, a amplitude aqui é 1, enquanto o complexo normal tem uma amplitude que é uma variável aleatória com uma distribuição Rayleigh.
Robert L.

@OlliNiemitalo tinha uma resposta para algo semelhante como eu me lembro ...
Fat32

@ Fat32 Sim, isso toca uma campainha. A outra pergunta é: qual é a distribuição dele?
Olli Niemitalo 16/11

@OlliNiemitalo Sim, essa era a pergunta!
precisa

Respostas:


6

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 :rrii

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π1z2

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

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(θ)=1

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

Dada essa relação, podemos resolver em termos de :rrii

r2=1i2

r2=1i2
r=±1i2
r=±1i2

Portanto, o pdf marginal de dado um valor de é um par de impulsos:rrii

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

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

Colocá-los juntos produziria:

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

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

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=i0i0[1,1]i0[1,1]ririr0=±1i20r0=±1i20

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.ririr=r0r=r0

Acredito que essa formulação seja equivalente à da resposta de AlexTP , mas sua derivação é provavelmente mais intuitiva.


1
existe um erro de digitação em ; possivelmente você quer dizer ? Também podemos mostrar que o pdf conjunto se integra ao plano 1 em xy (ou eqv ri). Novamente, também a articulação pdf não é circularmente simétrica? (o denominador?). Não podemos simplesmente usar a lógica que forneceu com uniforme em [ ] e o a restrição simplesmente oferece uma junta circular simétrica do tipo ou ; um impulso de anel no plano xy? onde K é possivelmente ?r=1i2r=1i2r2=1i2r2=1i2fri(r,i)fri(r,i)z=ejθ=cos(θ)+jsin(θ)=x+jyz=ejθ=cos(θ)+jsin(θ)=x+jyθθ0,2π0,2πx2+y2=1x2+y2=1fxy(x,y)=Kδ2(x2+y21)fxy(x,y)=Kδ2(x2+y21)fR,θ=Kδ2(R1)fR,θ=Kδ2(R1)1/2π1/2π
precisa

1
Corrigi os erros de digitação a que você se referiu, obrigado.
Jason R

@ Fat32: Fiz algumas correções que me levaram ao que parece ser o que você estava recebendo no seu comentário acima. A resposta da AlexTP ainda é provavelmente mais intuitivamente agradável.
Jason R

1
Há um erro de digitação ao copiar para . E observo que é diferente de . Eu acho que eles deveriam ter o mesmo valor, não deveriam? fr|i(r|i)fr|i(r|i)fr,i(r,i)fr,i(r,i)fr,i(0,1)=fr,i(0,1)=fr,i(1,0)=12πfr,i(1,0)=12π
AlexTP #

1
pessoal, talvez uma simulação rápida seja útil para mostrar seus pontos :) #
AlexTP

9

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) .XXYYZZVV V(XX2+Y2,YX2+Y2)

V(XX2+Y2,YX2+Y2)
V=1V=1VV

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(n1)(n1)

Assim, o PDF de é simplesmente o inverso da circunferência do círculo unitário. Para com fixo e uniforme ,ZZZρ=ρejΘZρ=ρejΘρρΘΘ

em coordenadas polares (onde a área infinitesimal é ), rdrdθrdrdθfR,Θ(r,θ)=12πδ(rρ)

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


2
Eu acho que está correto.
precisa

A intuição parece ser tão simples por trás da prova. Tudo que você precisa fazer é escolher pontos da superfície de uma esfera unidade :)
Maxtron

Não podemos simplesmente remover ono denominador? Caso contrário, como você avaliaria a expressão para ? Obviamente, gostaríamos que permanecesse. ||z||||z||||z||0||z||0fZ(0)=0fZ(0)=0
Matt L.

@MattL. você está certo, é errado ter no denominador, não devido à avaliação nos limites, mas porque é simplesmente errado. Eu consertei o resultado final. zz
AlexTP #

1
@ AlexTP: Eu ainda acho que não deveria haver nenhuma variável no denominador. Dê uma olhada na minha resposta e me diga se / onde eu errei.
Matt L.

4

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+y2

fZ(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=ρ

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

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

Pode ser mostrado que conduz às densidades marginais correcção para as variáveis aleatórias e .(3)(3)XXYY


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)rrZZrrZZxxyy

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δ(xxi)|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)=xx2+y2=xr

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

As duas raízes xEuxi estamos

x1,2=±ρ2-y2

x1,2=±ρ2y2(7)

Consequentemente,

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

|g(x1)|=|g(x2)|=ρ2y2ρ=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ρ2y2)+δ(x+ρ2y2)](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 θθ.


Eu acho que o determinante da transformação jacobiana deve ser rr então 2πθ=0 0+r=0 0fZ(r)rdrdθ1
2πθ=0+r=0fZ(r)rdrdθ1
AlexTP

@AlexTP: Tome como exemplo extremo ρ=0 0ρ=0, o que torna o RV ZZdeterminístico, é sempre zero. Seu PDF deve ser um impulso Dirac emr=z=y=0 0r=z=y=0, que é o que recebo com a fórmula que sugeri. Com a fórmula na sua resposta, obtemos um termo indefinido nesse caso.
Matt L.

é verdade, mas acho que o valor do PDF pode ser infinito e apenas o CDF deve ser limitado. Em outras palavras, o impulso Dirac importa apenas dentro de integrais.
AlexTP #

Sim, mas se você integrar seu pdf (por ρ=0 0ρ=0), o que você recebe (e como você o faria)?
Matt L.

Eu faria assim limρ0 0+θ+r=0 012πρδ(r-ρ)rdrdθ=limρ0 0+θ12πρρdθ=θ12πdθ=1
limρ0+θ+r=012πρδ(rρ)rdrdθ=limρ0+θ12πρρdθ=θ12πdθ=1
A integral é finita, enquanto o PDF não é graças ao determinante jacobiano. O que quero dizer é que você precisa desse determinante jacobianorrsempre que você integrar o PDF na coordenada polar. Diga-me se você não está convencido.
AlexTP #
Ao utilizar nosso site, você reconhece que leu e compreendeu nossa Política de Cookies e nossa Política de Privacidade.
Licensed under cc by-sa 3.0 with attribution required.