A distribuição para o produto Kronecker de dois vetores aleatórios uniformes na esfera unitária?


7

Suponha que dois vetores aleatórios e sejam distribuídos uniformemente na esfera unitária . É possível demonstrar que o produto de Kronecker de e é uniformemente distribuído sobre um subconjunto da esfera unidade ?xySn−1xySn−1⊗Sn−1

Respostas:


4

Sim. Isso fica claro ao trabalhar com as definições.

O produto Kronecker

A esfera unitária é o conjunto de vetores de unidades no espaço euclidiano onde é a norma euclidiana.Sn−1En=(Rn,||⋅||)

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

O "produto Kronecker" é o produto tensorial usual. Existem várias maneiras de pensar e calcular com ele. Um é defini-lo como uma matrizn×n

x⊗y=xy′=(x1y1x1y2⋯x1ynx2y1x2y2⋯x2yn⋮⋮⋯⋮xny1xny2⋯xnyn)∈Mat⁡(Rn,Rn).

Outra maneira equivalente desvenda os componentes dessa matriz em um vetor com componentes, permitindo visualizar como um elemento de . Observe que a métrica euclidiana para pode ser escritan2x⊗yRn2Rn2

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

Ambas são maneiras de escrever a soma dos quadrados de todos os componentes.n2

O produto Kronecker é compatível com as métricas euclidianas de e no sentido de queEnEn2

||x⊗y||2=||x||2||y||2.

Isso é facilmente demonstrado, porque o lado esquerdo é definido como a soma dos quadrados de todos os enquanto o lado direito é o produto de suas somas de quadrados. Apenas expanda esse produto:xiyj

||x⊗y||2=∑i∑j(xiyj)2=(∑ixi2)(∑jyj2)=||x||2||y||2.

Em particular, quando e têm comprimento unitário, possui comprimento unitário. Portantoxyx⊗y

Sn−1⊗Sn−1⊂Sn2−1.

Distribuições uniformes

As esferas euclidianas herdam uma medida da medida usual nos espaços euclidianos (a medida de Lebesgue, determinada em última análise pela distância euclidiana). Essa medida é preservada por qualquer isometria de uma esfera, porque (por definição) uma isometria preserva distâncias e a medida é finalmente determinada por distâncias. O grupo de isometrias da esfera unitária em é denotado , o grupo ortogonal. É um resultado clássico, e direto para mostrar, que consiste nas transformações lineares representadas por todas as matrizes para as quais .RmO(m)m×mPPP′=P′P=Im

O grupo ortogonal atua transitivamente em . (Aqui está uma prova de Euclides poderia ter feito: pegar quaisquer dois pontos distintos e na esfera Desenhe o segmento de linha. entre eles Ela determina uma perpendicular hyperplane única para. que passa pelo ponto médio de A reflexão em que os mapas hiperplano. para si e preserva todas as distâncias, onde é em . isso reflexão trocas de e , que mostra que existe uma transformação ortogonal envio de para ).O(n)Sn−1xyxyxyxySn−1O(n)xyxy

Dizer que "os vetores são distribuídos uniformemente em " significa que a distribuição é invariável sob um grupo transitivo de isometrias como .Sn−1O(n)

Aqui está o argumento: o "hypertorus" desfruta de um grupo transitivo de isometrias isomórficas para um quociente de . De fato, dado um arbitrários e outro , escolha para o qual e . Seja qualquer matriz e definaSn−1⊗Sn−1⊂Sn2−1O(n)×O(n)x⊗y∈Sn−1x2⊗y2∈Sn−1P,Q∈O(n)Px=x2Qy=y2Z=(zij)n×n

(2)(P⊗Q)(Z)=PZQ′.

Esta é uma isometria porque, usando a fórmula ,(1)

||(P⊗Q)(Z)||2=Tr⁡((PZQ′)(PZQ′)′)=Tr⁡(PZZ′P′)=Tr⁡(PP′ZZ′)=Tr⁡(ZZ′)=||Z||2.

Estes passos explorada a ortogonalidade de e e o facto de para quaisquer matrizes quadrados e .PQTr⁡(AB)=Tr⁡(BA)AB

Como as isometrias de induzem um grupo transitivo de isometrias de via , a distribuição uniforme (ou seja, invariável ao grupo) em é mapeado para uma distribuição uniforme em , QED .Sn−1Sn−1⊗Sn−1(2)Sn−1Sn−1⊗Sn−1


@ whuber Na sua prova, existem algumas confusões como segue: 1. a mão direita da equação (2) é uma matriz de -by- , como podemos entender a mão esquerda da equação (2)? nn
— Mao-lin Che

2. Por que você precisa escolher e ? Porque no restante do processo, você não usa esses dois vetores. x1⊗y1x2⊗y2
— Mao-lin Che

@ A Equação Mao-linChe define seu lado esquerdo. Não sei a que você se refere por " " porque isso não aparece na minha resposta. (2) x1⊗y1
— whuber

@ Whuber eu sou tão sombrio. Quero dizer por que você precisa para definir e . A matriz depende desses vetores? x⊗yx1⊗y1Z
— Mao-lin Che

@ Mao-linChe Esses quatro vetores determinam parcialmente as matrizes ortogonais e , conforme explicado imediatamente após a introdução desses vetores. não depende deles: isso é o que "nenhum" significa na frase "vamos ser qualquer matriz ...". PQZZ n×n
— whuber

0

Enquanto e são independentes, então o que você diz é verdade.xy

É mais fácil pensar em ângulos em vez de pontos no espaço.

Por exemplo, na esfera 2-D em , podemos gerar pontos uniformemente na superfície dessa esfera, amostrando e independentemente de e considerando-os como zênite e azimute, respectivamente.R3θ1θ2U[0,2π]

Assim, podemos pensar em como e como . O KP deles é então e fica claro a partir da construção que segue uma distribuição uniforme também.x(θ1,...,θn−1)∼U[0,2π]n−1y(ϕ1,...,ϕn−1)∼U[0,2π]n−1(θ1,...,θn−1,ϕ1,...,ϕn−1)U[0,2π]2n−2


11
Este argumento é fundamentalmente falho. O zênite é executado apenas de a e uma distribuição uniforme na esfera induz uma distribuição não uniforme do zênite. (Em , seu seno tem uma distribuição uniforme.) Consulte stats.stackexchange.com/questions/7977 para obter detalhes. −π/2π/2R3
— whuber

Estamos lidando com a superfície da esfera ou com todo o volume? Eu estava assumindo o primeiro.
— JDL

Isso está correto: é implicitamente identificado com o conjunto de vetores de unidades em . Sn−1Rn
— whuber

@whuber; você está certo! D'oh. Ainda acho que a pergunta é mais ou menos trivialmente verdadeira. O produto kronecker de duas distribuições uniformes independentes só pode ser uniforme no conjunto de produtos.
— JDL

11
@JDL você já tentou hist(runif(5000) %x% runif(5000))..?
— Tim
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.