问题

我计算了一个取决于两个变量的概率密度函数。我想使用此多元分布生成一些随机数，这些随机数的出现概率与 PDF 成正比。

看起来，SciPy 目前只支持 univariate distributions .是否有允许二维分布的简单方法或易于使用的软件包？

作为解决方法，我可能会尝试在感兴趣的域上创建随机数并将它们丢弃或保留它们与我的 PDF 相关的机会，但仍然可能有其他选择。随机数生成不必很快。

感谢您的帮助!

这是一个可能的解决方案

根据答案(非常感谢!)，我破解了一些您可能会在 this gist 中找到的代码.如果您使用 sin^2*Gauss PDF 运行此示例，则将在 PDF 上绘制满足给定条件(在圆圈内)的 2000 个随机随机变量。也许这对其他人也有帮助。

请您参考如下方法:

所以你有一个 PDF F(x,y) 并且你想生成根据此分布的 x 和 y 对PDF？

我会说，除非您可以使用多变量版本的反演技术 (wiki)，否则拒绝抽样是可行的方法。