def random_orthogonal_matrix(n):
"""Returns a random, orthogonal matrix of n by n."""
a = np.random.randn(n, n)
U, _, _ = linalg.svd(a)
assert U.shape == (n, n)
return U
微信
公众号
扫码关注公众号