使用BoxMuller算法生成符合正态分布随机数

正态分布随机数.py

# 从这抄的
# https://blog.csdn.net/Yonggie/article/details/97404027

import math
from random import random
import matplotlib.pyplot as plt

# mu:均值/期望
# sigma:标准差
def Generate(mu, sigma):
    data = []
    for i in range(10000):
        data.append(
            mu
            + sigma
            * (math.sqrt(-2 * math.log(random())) * math.sin(2 * math.pi * random()))
        )
    plt.figure()
    plt.hist(data, bins=50)
    plt.show()
    exit(0)



if __name__ == "__main__":
    # Generate(0,1)
    # 标准差0.25表示生成的随机数大多数落在[-1,1]区间内
    Generate(0,0.25)

.