如何生成正行列式的均匀随机正交矩阵?


9

p11det=1

我们可以通过将正交矩阵的的正负号更改为负号来更改它的任何一列(或更一般地,其任意奇数个列)。det

我的问题是:鉴于我们会反复生成这样的随机矩阵,如果每次我们选择只恢复特定列的符号(例如,始终为第一个或始终为最后一个),是否会在它们的统一随机性上引入一些偏差?还是我们必须随机选择列以使矩阵表示随机均匀分布的集合?

Answers:


7

列的选择无关紧要:在特殊正交矩阵上所得的分布仍然是均匀的。SO(n)

我将通过使用一种显而易见的论点来解释这一点,该论点扩展到有关组元素统一生成的许多相关问题。此参数的每个步骤都很简单,只需要引用适当的定义或简单的计算即可(例如,注意矩阵I1是正交的并且是自反的)。

该论点是对熟悉情况的概括。 考虑根据指定的连续分布F绘制实数的任务。这可以通过从连续分布G提取任何实数并在必要时取反结果来保证正值(几乎可以肯定)来完成。为了使该过程具有分布FG必须具有以下特性:FGFG

G(x)G(x)=F(x).

做到这一点的最简单的方法是,当是对称的周围0,使得G ^ X - 1 / 2 = 1 / 2 - g ^ - X ,将会导致˚F X = 2 g ^ X - 1:所有阳性简单地将概率密度加倍,并消除所有负面结果。半正态分布(F)和正态分布(GG0G(x)1/2=1/2G(x)F(x)=2G(x)1FG)就是这种类型。

在下文中,组充当非零实数(被视为乘法组)的角色,其子组S O n 充当正实数R +的角色。Haar度量d x / x在否定下是不变的,因此当它从R{ 0 }折叠到R +时O(n)SO(n)R+dx/xR{0}R+,正值的分布不变。(不幸的是,该量度无法归一化为概率量度,但这是类推失败的唯一方法。)

否定正交矩阵的特定列(当行列式为负时)类似于否定负实数以将其折叠为正子组。更一般而言,您可以提前选择任何一个负行列式正交矩阵并使用它代替I 1JI1:结果将是相同的。


尽管此问题是根据生成随机变量来表达的,但它实际上是询问矩阵组S O n R= S O n 上的概率分布。这些组之间的连接用正交矩阵描述O(n,R)=O(n)SO(n,R)=SO(n)

I1=(1000010000001)

因为否定的正交矩阵的第一列手段右乘以X 1。注意,小号ø Ñ Ô Ñ ø Ñ 是不相交并XXI1SO(n)O(n)O(n)

O(n)=SO(n)SO(n)I11.

给定的概率空间上定义ø Ñ ,在问题中描述的方法定义了一个地图(O(n),S,P)O(n)

f:O(n)SO(n)

通过设置

f(X)=X

XSO(n)

f(X)=XI1

对于XSO(n)I11

问题是关心生成随机元素通过获得随机元素ω ø Ñ :,通过由“推动它们向前” ˚F以产生˚F * ω = ˚F ω 小号Ô n 。的前推创建概率空间小号ø Ñ 小号 'P 'SO(n)ωO(n)ffω=f(ω)SO(n)(SO(n),S,P)

S=fS={f(E)|ES}

P(E)=(fP)(E)=P(f1(E))=P(EEI1)

所有ES

通过假定右乘法是保测,并注意到在任何情况下ë I1,它会立即跟随所有Ë小号EEI1=ES

P(E)=P(EEI11)=P(E)+P(EI11)=2P(E).

尤其是,当O n )的右乘下不变(这通常是“均匀”的意思)时,很明显的事实是I 1及其逆数(恰好等于I 1本身)都是正交的,这意味着前述成立,证明P也是均匀的。因此,没有必要选择一个随机列进行求反。PO(n)I1I1P


+1。这是一个非常不错的文章,感谢您发布此答案。
变形虫

1
一个了不起的答案。但是从The question is concerned about generating我开始,我很难通过象征主义来推动我前进。您能用言语来概括推理的原因吗,请早点外行?
ttnphns
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.