我需要生成一个随机变量列表受约束，这些约束可以用的形式表示，其中是一个矩阵，如果有个条目。在我处理的所有情况下，，例如大约为 14,000，为 50。我不确定我将使用什么方法进行随机抽样，无论是正常的还是统一的，它不是清楚哪个最适合我要解决的问题，但我需要从具有相同均值和范围/方差的分布中采样所有变量。$\bf{x}$$\bf{E}x=b$$\bf{E}$$m \times n$$\bf{x}$$n$$n >> m$$n$$m$

我一直在做的解决这个问题是将减少为行梯形形式，将与最后一个枢轴右侧的列对应的所有变量设置为随机值，然后解决剩余的方阵相等。$\bf{E}$

但是有一个问题，为了解决方阵相等，我从右侧减去已经设置的值。不幸的是，方差也增加了，所以我最后的 50 个值往往变化很大，不幸的是，这在这个问题中是不可接受的。

有一个更好的方法吗？我想不出一种方法来修复我正在使用的当前方法。我用R。