scipy.stats.ortho_group

scipy.stats.ortho_group = <scipy.stats._multivariate.ortho_group_gen object>[源代码]

矩阵值O(N)随机变量。

返回一个随机正交矩阵,取自O(N)Haar分布(O(N)上唯一的均匀分布)。

这个 dim 关键字指定尺寸N。

参数
dim标量

矩阵的维数

注意事项

这门课与以下内容密切相关 special_ortho_group

根据Mezzadri的论文,采取了一些谨慎的措施来避免数字错误。

参考文献

1

F.Mezzadri,“如何从经典紧群生成随机矩阵”, arXiv:math-ph/0609050v2

示例

>>> from scipy.stats import ortho_group
>>> x = ortho_group.rvs(3)
>>> np.dot(x, x.T)
array([[  1.00000000e+00,   1.13231364e-17,  -2.86852790e-16],
       [  1.13231364e-17,   1.00000000e+00,  -1.46845020e-16],
       [ -2.86852790e-16,  -1.46845020e-16,   1.00000000e+00]])
>>> import scipy.linalg
>>> np.fabs(scipy.linalg.det(x))
1.0

这将从O(3)生成一个随机矩阵。它是正交的,行列式为+1或-1。

方法:

``rvs(dim=None, size=1, random_state=None)``

从O(N)中随机抽取样本。