Verteilung von Skalarprodukten zweier zufälliger Einheitsvektoren in Dimensionen

27

Wenn und zwei unabhängige zufällige Einheitsvektoren in (gleichmäßig auf einer Einheitskugel verteilt), wie lautet die Verteilung ihres Skalarprodukts (Skalarprodukt) ?xyRDxy

Ich vermute, als wächst, wird die Verteilung schnell normal (?), Wobei der Mittelwert Null und die Varianz in höheren Dimensionen abnehmen. aber gibt es eine explizite Formel für \ Sigma ^ 2 (D) ?D σ 2 (D)

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

Aktualisieren

Ich habe einige schnelle Simulationen durchgeführt. Wenn man 10000 Paare zufälliger Einheitsvektoren für D=1000 , ist leicht zu erkennen, dass die Verteilung ihrer Punktprodukte perfekt Gauß ist (tatsächlich ist sie bereits für D = 100 ziemlich Gauß D=100), siehe die Nebenkurve auf der linken Seite. Zweitens habe ich für jedes D Bereich von 1 bis 10000 (mit zunehmenden Schritten) 1000 Paare generiert und die Varianz berechnet. Log-Log-Plot wird rechts gezeigt, und es ist klar, dass die Formel sehr gut durch 1 / D angenähert ist 1/D. Beachten Sie, dass diese Formel für D=1 und D=2 sogar genaue Ergebnisse liefert (aber ich bin nicht sicher, was später passiert).

Punktprodukte zwischen zufälligen Einheitsvektoren

Amöbe sagt Reinstate Monica
quelle
@KarlOskar: Danke, dieser Link ist sehr relevant, und in der Tat macht meine Frage fast ein Duplikat, aber nicht ganz. Es gibt also eine explizite Formel für die eine kumulative Verteilungsfunktion der Punktprodukte ist. Man kann eine Ableitung nehmen, um das PDF zu erhalten und dann das Limit studieren . Die Formel wird jedoch in Form von Betafunktionen und unvollständigen Betafunktionen angegeben, sodass die Berechnungen wahrscheinlich unangenehm sind. D P{(x,y)>ϵ}D
Amöbe sagt Reinstate Monica
@KarlOskar: von der gleichmäßigen Verteilung auf einem Einheitskugel in . Um einen Zufallsvektor aus dieser Verteilung zu erzeugen, kann man einen Zufallsvektor aus einem Gaußschen mit einer Einheitsvarianz erzeugen und ihn dann normalisieren. RD
Amöbe sagt Reinstate Monica

Antworten:

30

Da ( wie allgemein bekannt ) eine gleichmäßige Verteilung auf der Einheitskugel durch Normalisieren einer Variantennormalverteilung erhalten wird und das Punktprodukt der normalisierten Vektoren ihr Korrelationskoeffizient ist, sind die Antworten auf die drei Fragen sind: DSD1Dt

  1. ( ( D - 1 ) / 2 , ( D - 1 ) / 2 )u=(t+1)/2 hat eine Beta -Verteilung.((D1)/2,(D1)/2)

  2. Die Varianz von ist (wie in der Frage spekuliert).1 / Dt1/D

  3. Die standardisierte Verteilung von nähert sich der Normalität mit einer Rate vonO ( 1tO(1D).


Methode

Die genaue Verteilung des Punktprodukts von Einheitsvektoren ist geometrisch leicht zu erhalten, da dies die Komponente des zweiten Vektors in Richtung des ersten ist. Da der zweite Vektor unabhängig vom ersten ist und gleichmäßig auf der Einheitskugel verteilt ist, ist seine Komponente in der ersten Richtung genauso verteilt wie jede Koordinate der Kugel. (Beachten Sie, dass die Verteilung des ersten Vektors keine Rolle spielt.)

Die Dichte finden

Wenn diese Koordinate die letzte ist, ist die Dichte bei daher proportional zur Oberfläche, die in einer Höhe zwischen und auf der Einheitskugel liegt. Dieser Anteil tritt innerhalb eines Riemens mit der Höhe und dem Radius der im Wesentlichen ein Kegelstumpf ist , der aus einem mit dem Radius der Höhe und Steigung . Woher ist die Wahrscheinlichkeit proportional zut t + d t d t t[1,1]tt+dtdtS D - 2 1t2,SD2dt1/1t2,dt1/1t2

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

Wenn ist, ist . Wenn man das in das Vorhergehende einfügt, ergibt sich für das Wahrscheinlichkeitselement eine normalisierende Konstante:T = 2 u - 1u=(t+1)/2[0,1]t=2u1

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

Es ist unmittelbar, dass eine Beta -Verteilung hat, weil (per Definition) auch seine Dichte proportional zu ist( ( D - 1 ) / 2 , ( D - 1 ) / 2 )u=(t+1)/2((D1)/2,(D1)/2)

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

Festlegen des Begrenzungsverhaltens

Daraus ergeben sich mit elementaren Techniken leicht Informationen über das Begrenzungsverhalten: kann integriert werden, um die Proportionalitätskonstante ; kann integriert werden (z. B. unter Verwendung von Eigenschaften von Beta-Funktionen), um Momente zu erhalten, die zeigen, dass die Varianz und auf schrumpft (ausgehend von Chebyshevs Theorem wird die Wahrscheinlichkeit in der Nähe von konzentriert) ); und die Grenzverteilung wird dann gefunden, indem Werte der Dichte der standardisierten Verteilung, proportional zu für kleine Werte von berücksichtigt werdenΓ ( nfDtkfD(t)1/D0t=0fD(t/Γ(n2)πΓ(D12)tkfD(t)1/D0t=0tfD(t/D),t :

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

wobei die 's (log) Integrationskonstanten darstellen. Offensichtlich ist die Rate, mit der sich dies der Normalität nähert (für die die logarithmische Dichte gleich ),- 1CO(112t2O(1D).

Zahl

Dieses Diagramm zeigt die Dichten des Punktprodukts für , wie auf Einheitsvarianz standardisiert, und ihre Grenzdichte. Die Werte bei erhöhen sich mit (von blau über rot, gold und dann grün für die normale Standarddichte). Die Dichte für wäre bei dieser Auflösung nicht von der normalen Dichte zu unterscheiden.0 D D = 1000D=4,6,100DD=1000

whuber
quelle
4
(+1) Vielen Dank, @whuber, das ist eine großartige Antwort! Besonderer Dank für die Erwähnung des Wortes "Kegelstumpf". Es ist so, dass ich eine weitere Antwort nur wenige Minuten vor Ihrem Posting angenommen habe und ich sie jetzt nicht mehr ablehnen möchte. Ich hoffe du verstehst. Schade, dass es nicht möglich ist, beides zu akzeptieren! Übrigens, beachten Sie einen sehr einfachen Beweis des Ausdrucks für die Abweichung von dieser Antwort: Man kann ihn direkt sehen, ohne mit Betafunktionen herumzuspielen! Varianz des Punktprodukts zu Varianz jeder Kugel koordinieren gleich (wie Sie geschrieben haben ) und eine Summe aller von ihnen sein sollte , QEDD 11/DD1
Amöbe sagt wieder einzusetzen Monica
1
Das ist eine schöne Beobachtung über die Abweichungen.
whuber
2
@amoeba, die jüngste Aktivität hat auch hier wieder meine Aufmerksamkeit erregt, und so sehr ich es schätze, dass Sie meine Antwort akzeptiert haben, so viel voller ist diese. Es würde mir nichts ausmachen, wenn Sie sich ändern würden.
Freitag,
1
@ Student001: das ist ein fairer und großzügiger Kommentar. Ich habe die akzeptierte Antwort gewechselt. Ich habe auch ein Q und ein A von Ihnen gefunden, um das auszugleichen :)
Amöbe sagt Reinstate Monica
1
@mat Die Verteilung von ist die von . Das macht es zu einer Beta-Verteilung, die skaliert und vom Intervall in das Intervall verschoben wurde . 2 U - 1 [ 0 , 1 ] [ - 1 , 1 ]t2U1[0,1][1,1]
whuber
11

Lassen Sie uns anhand der Standardergebnisse die Verteilung und dann die Varianz bestimmen. Betrachten Sie das Vektorprodukt und schreiben Sie es in seine Kosinusform, dh beachten Sie, dass wir wobei der Winkel zwischen und . Im letzten Schritt habe ich das für alle Ereignisse undBetrachte nun den Term . Es ist klar, dass es keine Rolle spielt, was ist, da Bezug auf die Kugeloberfläche gleichmäßig gewählt wirdθ x y A B E P ( A B ) : = E [ E [ χ

P(xyt)=P(|x||y|cosθt)=P(cosθt)=EP(cosθty),
θxyAB P ( cos & thgr; t | y ) x y x y y y = [ 1 , 0 , 0 , ... ] " . P ( x ' y t ) = P ( x 1t ) . x 1
EP(AB):=E[E[χAB]]=EχA=P(A).
P(cosθty)xyeigentlich ist nur der winkel zwischen und bedeutung. Somit ist der Term innerhalb der Erwartung als Funktion von tatsächlich konstant und wir können annehmen, dassDann erhalten wiraber da die erste Koordinate eines normalisierten in wir, wir das asymptotische Ergebnis dieser Arbeit aufrufen dass Gauß mit der Varianz .xyyy=[1,0,0,].
P(xyt)=P(x1t).
x1 x ' y1 / nRn,xy1/n

Um ein explizites Ergebnis der Varianz zu erhalten, verwenden Sie die Tatsache, dass das Skalarprodukt Null als Unabhängigkeit bedeutet und, wie oben gezeigt, wie die erste Koordinate von . Nach diesen Ergebnissen entspricht das Finden von dem Finden von . Beachten Sie nun, dass pro Konstruktion und somit wobei die letzte Gleichheit daraus folgt, dass die Koordinaten von identisch verteilt sind. Beim Zusammenfügen haben wir festgestellt, dassVar ( x ' y ) E x 2 1 x ' x = 1 1 = E x ' x = E n i = 1 x 2 i = n i = 1 E x 2 i = n E x 2 1 , x Var ( x ' y ) = E x 2 1xVar(xy)Ex12xx=1

1=Exx=Ei=1nxi2=i=1nExi2=nEx12,
xVar(xy)=Ex12=1/n
ekvall
quelle
Danke, aber ich bin verwirrt: Was genau ist "das gewünschte Ergebnis" und wie folgt es aus der letzten Gleichung? Die endgültige Wahrscheinlichkeitsverteilung sollte von abhängen . D
Amöbe sagt Reinstate Monica
Wie das Ergebnis aus Ihrer letzten Gleichung hervorgeht, ist genau das, was Sie im math.SE-Thread gefunden haben. Es handelt sich um Betaverteilungen usw., und das einschränkende Verhalten ist (für mich) alles andere als offensichtlich. Ich denke, es sollte einen einfacheren direkten Weg geben, um dieses . σ2(D)1/D
Amöbe sagt Reinstate Monica
Dies hängt von der Dimension ab, da , wobei der generierte Gaußsche Vektor ist. Ich werde die Antwort später heute oder morgen aktualisieren. x1=z1|z|1z
ekvall
Wow, toll, Ihr letztes Glied stellt die Grenze dieses Ausdrucks , den inversen Beta - Funktionen (was ich habe Angst , zu berechnen) in der dritten Gleichung auf Seite 1. So die Argumentation zu vervollständigen: Wenn die Kugel des Radius , dann wird (asymptotisch) als . Das bedeutet , dass für den Bereich des Einheitsradius Varianz - mal kleiner, dh . Ich habe jedoch immer noch Bedenken: Ich habe nach von 1 bis 4 gesucht, und scheint eine exakte Varianz zu ergeben , obwohl die Verteilungen für D = 1 oder D = 2 weit vom Normalen entfernt sind. Es sollte einen tieferen Grund dafür geben. Dx1N(0,1)D1/DD1/D
Amöbe sagt Reinstate Monica
@amoeba Ja, aktualisiert mit einem Beweis dafür.
Freitag,
2

Um den ersten Teil Ihrer Frage zu beantworten, bezeichnen Sie . Definiere Das Produkt der -Elemente von und die hier als werden, werden gemäß der gemeinsamen Verteilung von und . dann seit , 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

Für den zweiten Teil denke ich, dass Sie, wenn Sie etwas Interessantes über das asymptotische Verhalten von sagen möchten, mindestens die Unabhängigkeit von und annehmen und dann eine CLT anwenden müssen.σXY

Wenn Sie beispielsweise annehmen , dass die mit und Sie dies Sagen Sie, dass und .{Z1,,ZD}E[Zi]=μV[Zi]=σ2σ2(D)=σ2DlimDσ2(D)=0

tom
quelle
Danke, aber ich bin verwirrt über den zweiten Teil. und sollen natürlich unabhängig sein, das werde ich der Frage hinzufügen. Sie sagen, dass , und das klingt vernünftig, aber wie ist das asymptotische Verhalten von ? Ich denke, der gesuchte Ausdruck sollte nur von abhängen . Übrigens in 2D wenn ich mich nicht irre, frage ich mich, ob dies in höheren Dimensionen wahr bleibt ...Y & sgr; 2 ( D ) = V a r ( z i ) / D V a r ( z i ) D V a r ( Z i ) = 1 / 2XYσ2(D)=Var(zi)/DVar(zi)DVar(zi)=1/2
Amöbe sagt Reinstate Monica
Ist es wirklich möglich, dass unabhängig ist, , dass und eine Längeneinheit haben? X YziXY
ekvall
@ Tom: Übrigens, ich war verwechselt: in 2D 1 ist, ist es , die gleich 1/2 ist. Ich habe meine Frage mit einigen Simulationsergebnissen aktualisiert. Die richtige Formel scheint . V a r ( z ) , 1 / DVar(zi)Var(z)1/D
Amöbe sagt Reinstate Monica