Python SciPy stats.unitary_group用法及代碼示例

本文簡要介紹 python 語言中 scipy.stats.unitary_group 的用法。

用法:

scipy.stats.unitary_group = <scipy.stats._multivariate.unitary_group_gen object>#

matrix-valued U(N) 隨機變量。

返回一個隨機酉矩陣。

dim 關鍵字指定維度 N。

參數 ：：

dim： 標量

矩陣的維數

seed： {無，int，np.random.RandomState，np.random.Generator}，可選

用於繪製隨機變量。如果種子是None， 這RandomState使用單例。如果種子是一個 int，一個新的RandomState使用實例，用種子播種。如果種子已經是一個RandomState或者Generator實例，然後使用該對象。默認為None.

注意：

此類類似於 ortho_group 。

參考：

[1]

F. Mezzadri，“如何從經典緊群生成隨機矩陣”，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]])

這會從 U(3) 生成一個隨機矩陣。點積確認它是單一的，直到機器精度。

或者，可以調用該對象(作為函數)來修複 dim 參數，返回 “frozen” unitary_group 隨機變量：

>>> rv = unitary_group(5)