Distribución de productos escalares de dos vectores unitarios aleatorios en dimensiones

27

Si x e y son dos vectores unitarios aleatorios independientes en RD (distribuidos uniformemente en una esfera unitaria), ¿cuál es la distribución de su producto escalar (producto de puntos) xy ?

D

limDσ2(D)0,
σ2(D)

Actualizar

Ejecuté algunas simulaciones rápidas. Primero, al generar 10000 pares de vectores unitarios aleatorios para , es fácil ver que la distribución de sus productos de puntos es perfectamente gaussiana (de hecho, ya es bastante gaussiana para ), vea la subtrama a la izquierda. En segundo lugar, para cada varía de 1 a 10000 (con pasos crecientes), generé 1000 pares y calculé la varianza. Gráfico log-log se muestra a la derecha, y es evidente que la fórmula está muy bien aproximada por . Tenga en cuenta que para y esta fórmula incluso da resultados exactos (pero no estoy seguro de lo que sucede más adelante).D=1000D=100D1/DD=1D=2

productos de punto entre vectores unitarios aleatorios

ameba dice Reinstate Monica
fuente
@KarlOskar: gracias, este enlace es muy relevante y, de hecho, hace que mi pregunta sea casi un duplicado, pero no del todo. Entonces, hay una fórmula explícita para que es una función de distribución acumulativa de los productos de punto. Uno puede tomar una derivada para obtener el PDF y luego estudiar el límite . Sin embargo, la fórmula se da en términos de funciones beta y funciones beta incompletas, por lo que es probable que los cálculos sean desagradables. P{(x,y)>ϵ}D
ameba dice Reinstate Monica
@KarlOskar: de la distribución uniforme en una esfera unidad en . Para generar un vector aleatorio a partir de esta distribución, se puede generar un vector aleatorio a partir de un gaussiano con una varianza unitaria, y luego normalizarlo. RD
ameba dice Reinstate Monica

Respuestas:

30

Debido a que ( como es bien sabido ) se obtiene una distribución uniforme en la esfera de la unidad normalizando una distribución normal de variación y el producto de punto de los vectores normalizados es su coeficiente de correlación, las respuestas a los tres las preguntas son:SD1Dt

  1. u=(t+1)/2 tiene una distribución Beta .((D1)/2,(D1)/2)

  2. La varianza de es igual a (como se especula en la pregunta).t1/D

  3. La distribución estandarizada de aproxima a la normalidad a una tasa detO(1D).


Método

La distribución exacta del producto escalar de los vectores unitarios se obtiene fácilmente geométricamente, porque este es el componente del segundo vector en la dirección del primero. Dado que el segundo vector es independiente del primero y se distribuye uniformemente en la esfera de la unidad, su componente en la primera dirección se distribuye de la misma manera que cualquier coordenada de la esfera. (Tenga en cuenta que la distribución del primer vector no importa).

Encontrar la densidad

Dejando que esa coordenada sea la última, la densidad en es, por lo tanto, proporcional al área de superficie que se encuentra a una altura entre y en la esfera de la unidad. Esa proporción ocurre dentro de un cinturón de altura y radio que es esencialmente un tronco cónico construido a partir de un de radio de altura y pendiente . De donde la probabilidad es proporcional at[1,1]tt+dtdt1t2,SD21t2,dt1/1t2

(1t2)D21t2dt=(1t2)(D3)/2dt.

Dejar implica . Sustituyendo eso en el precedente da el elemento de probabilidad hasta una constante de normalización:u=(t+1)/2[0,1]t=2u1

fD(u)du(1(2u1)2)(D3)/2d(2u1)=2D2(uu2)(D3)/2du.

Es inmediato que tiene una distribución Beta , porque (por definición) su densidad también es proporcional au=(t+1)/2((D1)/2,(D1)/2)

u(D1)/21(1u)(D1)/21=(uu2)(D3)/2fD(u).

Determinando el comportamiento limitante

La información sobre el comportamiento limitante se deduce fácilmente de esto mediante técnicas elementales: se puede integrar para obtener la constante de proporcionalidad ; puede integrarse (usando las propiedades de las funciones Beta, por ejemplo) para obtener momentos, mostrando que la varianza es y se reduce a (de donde, por el Teorema de Chebyshev, la probabilidad se concentra cerca de ); y la distribución limitante se encuentra luego considerando los valores de la densidad de la distribución estandarizada, proporcional a para valores pequeños defDΓ(n2)πΓ(D12)tkfD(t)1/D0t=0fD(t/D),t :

log(fD(t/D))=C(D)+D32log(1t2D)=C(D)(1/2+32D)t2+O(t4D)C12t2

donde las representan constantes (log) de integración. Evidentemente, la velocidad a la que esto se aproxima a la normalidad (para la cual la densidad logarítmica es igual a ) esC12t2O(1D).

Figura

Este gráfico muestra las densidades del producto escalar para , según la varianza estandarizada, y su densidad limitante. Los valores en aumentan con (de azul a rojo, dorado y luego verde para la densidad normal estándar). La densidad para sería indistinguible de la densidad normal en esta resolución.D=4,6,100DD=1000

whuber
fuente
44
(+1) Muchas gracias, @whuber, ¡esta es una gran respuesta! Un agradecimiento especial por mencionar la palabra "frustum". Sucede que acepté otra respuesta unos minutos antes de que publicaras la tuya, y no me gustaría cancelarla ahora; espero que entiendas. Lástima que no sea posible aceptar ambos! Por cierto, observe una prueba muy simple de la expresión para la varianza de esa respuesta: ¡uno puede verla directamente sin perder el tiempo con las funciones beta! La varianza del producto punto es igual a la varianza de cualquier coordenada de esfera (como usted escribió), y una suma de todos los de ellos debería ser , QED1/DD1
ameba dice Reinstate Monica
1
Esa es una buena observación sobre las variaciones.
whuber
2
@amoeba, la actividad reciente también me llamó la atención aquí, y aunque aprecio que hayas aceptado mi respuesta, esta es mucho más completa. No me importaría en absoluto si cambiaras.
ekvall
1
@ Student001: este es un comentario justo y generoso. Cambié la respuesta aceptada. También he encontrado una Q y una A para votar para compensarlo :)
ameba dice Reinstate Monica
1
@mat La distribución de es la de . Eso lo convierte en una distribución Beta que ha sido escalada y desplazada del intervalo al intervalo . t2U1[0,1][1,1]
whuber
11

Encontremos la distribución y luego la varianza sigue los resultados estándar. Considere el producto vectorial y escríbalo en su forma coseno, es decir, tenga en cuenta que donde es el ángulo entre e . En el último paso lo he usado para cualquier evento yAhora considere el término . Está claro que dado que se elige de manera uniforme con respecto a la superficie de la esfera, no importa qué

P(xyt)=P(|x||y|cosθt)=P(cosθt)=EP(cosθty),
θxyAB
EP(AB):=E[E[χAB]]=EχA=P(A).
P(cosθty)xyen realidad es, solo importa el ángulo entre e . Por lo tanto, el término dentro de la expectativa es en realidad constante en función de y podemos asumir queLuego obtenemos quepero dado que es la primera coordenada de un vector gaussiano normalizado en tenemos que es gaussiano con varianza invocando el resultado asintótico de este artículo .xyyy=[1,0,0,].
P(xyt)=P(x1t).
x1Rn,xy1/n

Para obtener un resultado explícito de la varianza, use el hecho de que el producto de punto es cero por independencia y, como se muestra arriba, se distribuye como la primera coordenada de . Según estos resultados, encontrar equivale a encontrar . Ahora, tenga en cuenta que por construcción y así podemos escribir donde la última igualdad se deduce de que las coordenadas de están idénticamente distribuidas. Al poner las cosas juntas, hemos encontrado quexVar(xy)Ex12xx=1

1=Exx=Ei=1nxi2=i=1nExi2=nEx12,
xVar(xy)=Ex12=1/n
ekvall
fuente
Gracias, pero estoy confundido: ¿cuál es exactamente "el resultado deseado" y cómo se deduce de la última ecuación? La distribución final de probabilidad debe depender de . D
ameba dice Reinstate Monica
En realidad, el resultado de su última ecuación es exactamente lo que se discutió en matemáticas . Implica distribuciones beta, etc., y el comportamiento limitante está (para mí) lejos de ser obvio. Creo que debe haber una manera directa más fácil de ver que . σ2(D)1/D
ameba dice Reinstate Monica
Depende de la dimensión ya que , donde es el vector gaussiano generado. Actualizaré la respuesta más tarde hoy o mañana. x1=z1|z|1z
ekvall
Wow, genial, tu último enlace proporciona el límite de esa expresión que involucra funciones beta inversas (que tenía miedo de calcular) en la tercera ecuación en la página 1. Entonces, para completar el razonamiento: si la esfera tiene radio , entonces se distribuye (asintóticamente) como . Lo que significa que para la esfera de la varianza radio unidad es veces más pequeño, es decir, . Sin embargo, todavía tengo una preocupación: verifiqué para de 1 a 4, y parece dar una variación exacta , a pesar de que las distribuciones para D = 1 o D = 2 están muy lejos de lo normal. Debería haber una razón más profunda detrás de eso. Dx1N(0,1)D1/DD1/D
ameba dice Reinstate Monica
@amoeba Sí, actualizado con una prueba de eso.
ekvall
2

Para responder la primera parte de su pregunta, denote . Definir El producto de los elementos de e indicados aquí ya que se distribuirá de acuerdo con la distribución conjunta de e . entonces desde , Z=X,Y=XiYi

fZi(zi)=fZ1,,ZD(z1,,zD)dzi
ithXYZiXiYi
fZi(zi)=fXi,Yi(x,zix)1|x|dx
Z=Zi
fZ(z)=fZ1,,ZD(z1,,zd)δ(zzi)dz1dzd

Para la segunda parte, creo que si quieres decir algo interesante sobre el comportamiento asintótico de necesitas al menos asumir la independencia de e , y luego aplicar un CLT.σXY

Por ejemplo, si estaba dispuesto a asumir que los se combinan con y podría digamos que y .{Z1,,ZD}E[Zi]=μV[Zi]=σ2σ2(D)=σ2DlimDσ2(D)=0

tom
fuente
Gracias, pero estoy confundido acerca de la segunda parte. supuesto, se supone que e son independientes, agregaré esto a la pregunta. Usted dice que , y eso suena razonable, pero ¿cuál es el comportamiento asintótico de ? Creo que la expresión Estoy buscando debería depender sólo de . Por cierto en 2D si no me equivoco, me pregunto si esto sigue siendo cierto en las dimensiones superiores ...XYσ2(D)=Var(zi)/DVar(zi)DVar(zi)=1/2
Amoeba dice Reinstate Monica
¿Es realmente posible que sea ​​independiente dado el requisito de que e son de longitud unitaria? ziXY
ekvall
@ Tom: Por cierto, yo estaba equivocada: en 2D es 1, es que es igual 1/2. He actualizado mi pregunta con algunos resultados de simulación. Parece que la fórmula correcta es . Var(zi)Var(z)1/D
ameba dice Reinstate Monica