如何在Python中生成截断正态分布的相关随机数？

Question

我正在尝试使用numpy.random.multivariate_normal()生成三个变量的相关随机数，其中均值和协方差矩阵（根据数据计算）作为输入。

正态分布在0和1之间被截断，因此，生成的随机数（对于所有三个变量）应该在0和1之间。但是，一些生成的随机数超出了界限。

如何控制为每个变量生成正态分布随机数的界限？

编辑：我可以使用scipy.stats.truncnorm独立地从三个截断的正态分布生成不相关的随机数。但是，我在这里寻找能产生相关随机数的东西。