兰道随机数发生器

Question

我正在尝试生成 1000 个具有 MPV（最可能值）为 25 的兰道分布的值，但在 scipy 或 numpy 中找不到兰道随机数生成器。 我尝试了 pylandau { pip install pylandau } 但这似乎只适合兰道斯而不生成随机数。 python 或 pyroot 中欢迎任何这样做的方式。

Answer 1

package pylandau只是在计算 Landau function 值。 它可用于拟合、随机数生成等。您似乎想知道如何从分布中生成随机数，正如此处的另一个问题所回答的那样。 以下示例从朗道分布中创建 100 个随机值：

import numpy as np
import pylandau

x = np.arange(-1, 1, 0.01)
y = pylandau.landau(x)
random_values = np.random.choice(x, size=100, p=y / y.sum())

Answer 2

Landaupy实现了朗道（和朗高斯）分布的纯 Python 实现，并允许生成随机样本。 代码将是

from landaupy import landau

samples = landau.sample(x_mpv=25, xi=1, n_samples=1000)

你可以安装这个

pip install git+https://github.com/SengerM/landaupy