Die Verteilung für das Kronecker-Produkt zweier einheitlicher Zufallsvektoren in der Einheitskugel?


7

Angenommen, zwei Zufallsvektoren und sind gleichmäßig auf der Einheitskugel . Kann gezeigt werden, dass das Kronecker-Produkt von und gleichmäßig auf einer Teilmenge der Einheitskugel ?xySn1xySn1Sn1

Antworten:


4

Ja. Dies wird beim Durcharbeiten der Definitionen deutlich.

Das Kronecker-Produkt

Die Einheitskugel ist die Menge der Einheitsvektoren im euklidischen Raum wobei ist die euklidische Norm.Sn1En=(Rn,||||)

||(x1,x2,,xn)||=x12+x22++xn2

Das "Kronecker-Produkt" ist das übliche Tensorprodukt. Es gibt verschiedene Möglichkeiten, darüber nachzudenken und damit zu rechnen. Eine besteht darin, es als Matrix zu definierenn×n

xy=xy=(x1y1x1y2x1ynx2y1x2y2x2ynxny1xny2xnyn)Mat(Rn,Rn).

Eine andere äquivalente Methode löst die Komponenten dieser Matrix in einen Vektor mit Komponenten, sodass wir als ein Element von . Beachten Sie, dass die euklidische Metrik für geschrieben werden kannn2xyRn2Rn2

(1)||Z||2=Tr(ZZ).

Beides sind Möglichkeiten, die Summe der Quadrate aller Komponenten zu schreiben .n2

Das Kronecker-Produkt ist in dem Sinne mit den euklidischen Metriken für und kompatibelEnEn2

||xy||2=||x||2||y||2.

Dies ist leicht zu demonstrieren, da die linke Seite als die Summe der Quadrate aller während die rechte Seite das Produkt ihrer Quadratsummen ist. Erweitern Sie einfach dieses Produkt:xichyj

||xy||2=ichj(xichyj)2=(ichxich2)(jyj2)=||x||2||y||2.

Insbesondere wenn sowohl als auch eine Einheitslänge haben, hat eine Einheitslänge. Deshalbxyxy

S.n- -1S.n- -1S.n2- -1.

Gleichmäßige Verteilung

Die euklidischen Kugeln erben ein Maß von dem üblichen Maß für die euklidischen Räume (das Lebesgue-Maß, das letztendlich durch den euklidischen Abstand bestimmt wird). Dieses Maß wird durch jede Isometrie einer Kugel erhalten, da (per Definition) eine Isometrie Entfernungen beibehält und das Maß letztendlich durch Entfernungen bestimmt wird. Die Gruppe der Isometrien der Einheitskugel in wird mit , der orthogonalen Gruppe. Es ist ein klassisches Ergebnis und einfach zu zeigen, dass es aus den linearen Transformationen besteht, die durch alle Matrizen für die .R.mÖ(m)m×mP.P.P.'=P.'P.=ichm

Die orthogonale Gruppe wirkt transitiv auf . (Hier ist ein Beweis, den Euklid möglicherweise erbracht hat: Wählen Sie zwei unterschiedliche Punkte und auf der Kugel aus. Zeichnen Sie das Liniensegment zwischen sich. Es bestimmt eine eindeutige Hyperebene senkrecht zu die durch den Mittelpunkt von . Die Reflexion in dieser Hyperebene wird abgebildet für sich und behält alle Abstände bei, von wo es in . Diese Reflexion tauscht und und zeigt, dass es eine orthogonale Transformation gibt, die an sendet .)Ö(n)S.n- -1xyxyxyxyS.n- -1Ö(n)xyxy

Zu sagen, dass "Vektoren gleichmäßig auf " bedeutet, dass die Verteilung unter einer transitiven Gruppe von Isometrien wie invariant ist .S.n- -1Ö(n)

Hier ist die Pointe: Der "Hypertorus" genießt eine transitive Gruppe von Isometrien, die isomorph zu einem Quotienten von . In der Tat, wenn ein beliebiges und ein anderes , wählen Sie für das und . Sei eine beliebige Matrix und definiereS.n- -1S.n- -1S.n2- -1Ö(n)×Ö(n)xyS.n- -1x2y2S.n- -1P.,Q.Ö(n)P.x=x2Q.y=y2Z.=(zichj)n×n

(2)(PQ)(Z)=PZQ.

Dies ist eine Isometrie weil, unter Verwendung der Formel ,(1)

||(PQ)(Z)||2=Tr((PZQ)(PZQ))=Tr(PZZP)=Tr(PPZZ)=Tr(ZZ)=||Z||2.

Diese Schritte nutzten die Orthogonalität von und und die Tatsache, dass für alle quadratischen Matrizen und .PQTr(AB)=Tr(BA)AB

Da die Isometrien von über eine transitive Gruppe von Isometrien von induzieren , ist die gleichmäßige (dh gruppeninvariante) Verteilung on wird auf eine gleichmäßige Verteilung auf abgebildet , QED .Sn1Sn1Sn1(2)Sn1Sn1Sn1


@ whuber In Ihrem Beweis gibt es einige Verwirrungen wie folgt: 1. Die rechte Hand von Gleichung (2) ist eine Matrix von by- . Wie können wir die linke Hand der Gleichung (2) verstehen? nn
Mao-lin Che

2. Warum müssen Sie und wählen ? Weil Sie im Rest des Prozesses diese beiden Vektoren nicht verwenden. x1y1x2y2
Mao-lin Che

@ Mao-linChe Gleichung definiert die linke Seite. Ich weiß nicht, worauf Sie sich mit " " beziehen, da dies in meiner Antwort nicht erscheint. (2) x1y1
whuber

@whuber ich bin soory. Ich meine, warum Sie und . Hängt die Matrix von diesen Vektoren ab? xyx1y1Z
Mao-lin Che

@ Mao-linChe Diese vier Vektoren bestimmen teilweise die orthogonalen Matrizen und , wie unmittelbar nach Einführung dieser Vektoren erläutert wird. hängt nicht von ihnen ab: Das ist es, was "any" in der Phrase "sei eine beliebige Matrix ..." bedeutet. PQZZ n×n
whuber

0

Solange und unabhängig sind, ist das, was Sie sagen, wahr.xy

Es ist am einfachsten, über Winkel und nicht über Punkte im Raum nachzudenken.

Beispielsweise in der 2-D - Bereich in , wir Punkte gleichmäßig auf der Oberfläche dieser Kugel durch Abtasten erzeugen können und unabhängig von sowie unter jenen wie Zenith und Azimut jeweils.R3θ1θ2U[0,2π]

Wir können uns also als und als . Ihr KP ist dann und aus der Konstruktion geht hervor, dass eine gleichmäßige Verteilung folgt ebenfalls.x(θ1,...,θn1)U[0,2π]n1y(ϕ1,...,ϕn1)U[0,2π]n1(θ1,...,θn1,ϕ1,...,ϕn1)U[0,2π]2n2


1
Dieses Argument ist grundlegend fehlerhaft. Der Zenit verläuft nur von bis und eine gleichmäßige Verteilung auf der Kugel induziert eine ungleichmäßige Verteilung des Zenits. (In hat der Sinus eine gleichmäßige Verteilung.) Weitere Informationen finden Sie unter stats.stackexchange.com/questions/7977 . π/2π/2R3
whuber

Beschäftigen wir uns mit der Oberfläche der Kugel oder dem gesamten Volumen? Ich nahm das erstere an.
JDL

Das ist richtig: wird implizit mit der Menge der Einheitsvektoren in identifiziert . Sn1Rn
whuber

@whuber; Du hast recht! D'oh. Ich denke immer noch, dass die Frage mehr oder weniger trivial wahr ist. Das Kronecker-Produkt zweier unabhängiger Gleichverteilungen kann immer nur auf dem Produktsatz einheitlich sein.
JDL

1
@ JDL hast du versucht hist(runif(5000) %x% runif(5000))..?
Tim
Durch die Nutzung unserer Website bestätigen Sie, dass Sie unsere Cookie-Richtlinie und Datenschutzrichtlinie gelesen und verstanden haben.
Licensed under cc by-sa 3.0 with attribution required.