La elección de la columna no importa: la distribución resultante en las matrices ortogonales especiales, SO ( n ) , sigue siendo uniforme.
Explicaré esto usando un argumento que se extiende, de manera obvia, a muchas preguntas relacionadas sobre la generación uniforme de elementos de grupos. Cada paso de este argumento es trivial, no requiere más que una referencia a definiciones adecuadas o un cálculo simple (como señalar que la matriz es ortogonal y autoinversa).yo1
El argumento es una generalización de una situación familiar. Considere la tarea de elaborar positivos números reales de acuerdo a una distribución continua especificada . Esto se puede hacer extrayendo cualquier número real de una distribución continua y negando el resultado, si es necesario, para garantizar un valor positivo (casi seguro). Para que este proceso tenga la distribución , debe tener la propiedad queG F GFsolFsol
G ( x ) - G ( - x ) = F( X ) .
La forma más sencilla de lograr esto es cuando es simétrico alrededor de modo que , lo que implica : toda probabilidad positiva las densidades simplemente se duplican y se eliminan todos los resultados negativos. La relación familiar entre la distribución medio normal ( ) y la distribución normal ( ) es de este tipo.0 G ( x ) - 1 / 2 = 1 / 2 - G ( - x ) F ( x ) = 2 G ( x ) - 1 F Gsol0 0G ( x ) - 1 / 2 = 1 / 2 - G ( - x )F( x ) = 2 G ( x ) - 1Fsol
A continuación, el grupo desempeña el papel de los números reales distintos de cero (considerados como un grupo multiplicativo ) y su subgrupo desempeña el papel de los números reales positivos . La medida de Haar es invariable bajo negación, por lo que cuando se "dobla" de a , la distribución de los valores positivos no cambia . (Esta medida, desafortunadamente, no se puede normalizar a una medida de probabilidad, pero esa es la única forma en que la analogía se rompe).S O ( n ) R + d x / x R - { 0 } R +O ( n )SO ( n )R+rex / xR -{0}R+
Negar una columna específica de una matriz ortogonal (cuando su determinante es negativo) es el análogo de negar un número real negativo para plegarlo en el subgrupo positivo. En términos más generales, puede elegir de antemano cualquier matriz ortogonal de determinante negativo y usarla en lugar de : los resultados serían los mismos.I 1Jyo1
Aunque la pregunta se formula en términos de generar variables aleatorias, realmente pregunta sobre las distribuciones de probabilidad en los grupos de matriz y . La conexión entre estos grupos se describe en términos de la matriz ortogonal.S O ( n , R ) = S O ( n )O ( n , R ) = O ( n )SO ( n , R ) = SO ( n )
yo1= ⎛⎝⎜⎜⎜⎜- 10 0⋮0 00 01⋮0 00 00 0⋮0 0............0 00 00 01⎞⎠⎟⎟⎟⎟
porque negar la primera columna de una matriz ortogonal significa multiplicar a la derecha por . Observe que y es la unión disjuntaX I 1 S O ( n ) ⊂ O ( n ) O ( n )XXyo1SO ( n ) ⊂ O ( n )O ( n )
O ( n ) = SO ( n ) ∪ SO ( n )yo- 11.
Dado un espacio de probabilidad definido en , el proceso descrito en la pregunta define un mapaO ( n )(O(n),S,P)O(n)
f:O(n)→SO(n)
configurando
f(X)=X
cuando yX∈SO(n)
f(X)=XI1
para .X∈SO(n)I1−1
La pregunta se refiere a generar elementos aleatorios en obteniendo elementos aleatorios : es decir, "empujándolos hacia adelante" a través de para producir . El avance crea un espacio de probabilidad conω ∈ O ( n ) f f ∗ ω = f ( ω ) ∈ S O ( n ) ( S O ( n ) , S ′ , P ′ )SO(n)ω∈O(n)ff∗ω=f(ω)∈SO(n)(SO(n),S′,P′)
S′=f∗S={f(E)|E⊂S}
y
P′(E)=(f∗P)(E)=P(f−1(E))=P(E∪EI1)
para todo .E⊂S′
Suponiendo que la multiplicación correcta por preserva la medida y observa que, en cualquier caso, , se inmediatamente que para todos , E∩EI1E∈ S ′E∩EI1=∅E∈S′
P′(E)=P(E∪EI−11)=P(E)+P(EI−11)=2P(E).
En particular, cuando es invariante bajo la multiplicación a la derecha en (que es lo que normalmente significa "uniforme"), el hecho obvio de que y su inverso (que es igual a sí) son ambas ortogonales significa lo anterior, lo que demuestra que es uniforme. Por lo tanto, no es necesario seleccionar una columna aleatoria para la negación. O ( n ) I 1 I 1 P ′PO(n)I1I1P′