Генерация случайных переменных с определенной функцией плотности вероятности с помощью Python ⇐ Python

[Гость]

Я хочу сгенерировать случайную величину с энергией Максвелла (например, f(ene)) или произвольным распределением.
Функция плотности вероятности имеет вид:

Код: Выделить всё

def f(ene):
le=3
return 2(ene/pi/le**3)*np.exp(-ene/le)

Я хотел бы сгенерировать 10000 образцов, например

Код: Выделить всё

f.rvs(scale,size)

который генерирует список, включающий 10000 элементов, функция плотности которых соответствует f(ene)

Как я мог сделать это?

Пример:
http://docs.scipy.org/doc/scipy/referen ... .htmlscipy обеспечивает распределение Максвелла:
с методом: rvs(loc=0, Scale=1, size=1) Случайные переменные.

Код: Выделить всё

print maxwell.rvs(1,10000)

сгенерирует 10000 выборок с распределением Максвелла. Но scipy предоставляет только некоторые функции. Теперь у меня есть еще одна функция, которой нет в scipy. Как я мог это сделать?

Подробнее здесь: https://stackoverflow.com/questions/283 ... -by-python