从python中的矩阵正态分布中随机抽取样本

Question

有没有办法使用python从矩阵正态分布（https://en.wikipedia.org/wiki/Matrix_normal_distribution#Drawing_values_from_the_distribution）中随机抽取样本？ Numpy具有从1-D正态分布和多元正态分布中绘制的功能，但我在矩阵正态分布上找不到任何东西。谢谢。

Answer 1

从您链接的维基百科页面：

矩阵法线通过以下方式与多元正态分布相关：

当且仅当

转换为numpy，应该是

numpy.random.multivariate_normal(M.ravel(), np.kron(V, U)).reshape(M.shape)

Answer 2

numpy.random.normal()可以绘制2d矩阵，只要您指定大小如下：

numpy.random.normal(loc = mean, scale = standard deviation, size = (m,n))