构造随机正交矩阵

Posted 2021-12-28 weiyinfu

import numpy as np
from scipy.stats import ortho_group

d = 3
seed = 1
size = 2
a, b = np.float32(ortho_group.rvs(size=size, dim=d, random_state=seed))
print(a)
print(np.linalg.det(a))
print(a @ a.T)

参考资料

https://docs.scipy.org/doc/scipy/reference/generated/scipy.stats.ortho_group.html