scipy.stats.unitary_group#

scipy.stats.unitary_group = <scipy.stats._multivariate.unitary_group_gen object>[源代码]#

一个矩阵值的 U(N) 随机变量。

返回一个随机酉矩阵。

dim 关键字指定维度 N。

参数:
dim标量

矩阵的维度。

seed{None, int, np.random.RandomState, np.random.Generator}, 可选

用于绘制随机变量。如果 seedNone,则使用 RandomState 单例。如果 seed 是一个整数,则使用一个新的 RandomState 实例,并使用 seed 进行播种。如果 seed 已经是一个 RandomStateGenerator 实例,则使用该对象。默认为 None

方法

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

从 U(N) 抽取随机样本。

参见

ortho_group

注释

此类类似于 ortho_group

参考

[1]

F. Mezzadri, “How to generate random matrices from the classical compact groups”, arXiv:math-ph/0609050v2.

示例

>>> import numpy as np
>>> from scipy.stats import unitary_group
>>> x = unitary_group.rvs(3)
>>> np.dot(x, x.conj().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]])  # may vary

这将从 U(3) 生成一个随机矩阵。点积确认它是酉矩阵,精度达到机器精度。

或者,可以调用该对象(作为函数)来固定 dim 参数,并返回一个“冻结的” unitary_group 随机变量

>>> rv = unitary_group(5)