我想知道如何使用特定的歸一化分布函數生成隨機數

Question

基本上，如果我的分布 function 是 f(v)= NormalDistribution(-u,sigma)+ NormalDistribution(u,sigma)

如何將 f 定義為 PDF，對其進行規范化，然后將一些隨機變量命令應用於我的 PDF？

Answer 1

好吧，你有兩個歸一化高斯的總和，

f(x) = N(x|μ,σ) + N(x|-μ,σ)

∫ f(x) dx = 1 + 1 = 2

PDF(x|μ,σ) = N(x|μ,σ)/2 + N(x|-μ,σ)/2

因為 PDF 是對稱的，所以采樣很簡單： select 一個高斯或另一個有 50% 的概率

沿線（未經測試）

import random

def sample(μ, σ):
    if random.random() < 0.5:
        return random.gauss(μ,σ)
    return random.gauss(-μ,σ)