鉴于生成 normal variates 是多么容易，我会这样做，然后直接转换为球坐标。

如果你需要，其中是维度，那么在球坐标中，，并且角度都是相互独立的和长度，在它们各自的范围内均匀分布用于第一个（仅与余弦一起使用的那个），用于其余的。如果是偶数，，您可以生成 ,。如果是奇数，，那么您需要为此添加另一个法线平方（参见上文如何生成它们），$N_p(0,I_p)$$p$$r\sim \chi^2_p = \Gamma(p/2,1/2)$$[0,\pi)$$[0,2\pi)$$p$$p=2m$$r_p = \sum_{k=1}^{m} (-\frac12) \ln U_k$$U_k \sim \mbox{i.i.d. } U[0,1]$$p$$p=2m+1$$r_p = r_{p-1} + z^2$ ,。$z\sim N(0,1)$

这本书在这里有这本书的参考。不是很便宜，但我认为它是您正在寻找的那个。;)