scipy.stats.unitary_group

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

一個矩陣值的 U(N) 隨機變數。

返回一個隨機么正矩陣。

dim 關鍵字指定維度 N。

參數:

dim純量

矩陣的維度，必須大於 1。

seed{None, int, np.random.RandomState, np.random.Generator}, optional

用於繪製隨機變數。如果 seed 為 None，則使用 RandomState 單例。如果 seed 是一個整數，則使用一個新的 RandomState 實例，並以 seed 作為種子。如果 seed 已經是一個 RandomState 或 Generator 實例，則使用該物件。預設值為 None。

另請參閱

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)

返回在新分頁中開啟

方法

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

從 U(N) 中抽取隨機樣本。