piecewise 函数随机取样的问题

Question

大家好，我正在对一个piecewise 函数进行随机取样，但是最后取出的数值总是来自于一边，另外一边总是取不到，不知道是怎么回事，下面是我的函数的具体细节：

ProbabilityDistribution[Piecewise[{{ (1/806.573)*Exp[-(5713.73 - x)/806.573], x < 5713.73}, {(1/806.573)*Exp[(5713.73 - x)/806.573], x >= 5713.73}}], {x, -\[Infinity], \[Infinity]}]

换句话说，我只能取到小于 5713.73的数值，大于5713.73的样品根本没有。

请大侠帮忙看看，小弟出学。

Answer 1

这是由于你定义的概率密度函数有问题，其在实数域上的积分等于2，所以除以2就可以了：

pdf = ProbabilityDistribution[
   Piecewise[{{(1/806.573)*Exp[-(5713.73 - x)/806.573]/2, 
      x < 5713.73}, {(1/806.573)*Exp[(5713.73 - x)/806.573]/2, 
      x >= 5713.73}}], {x, -\[Infinity], \[Infinity]}];
RandomVariate[pdf, 100]

Plot[PDF[pdf, x], {x, 5713.73 - 5004, 5713.73 + 5004}, 
 Filling -> Axis]