python生成正态分布随机数

        使用python基础库random生成具备正态分布特征的随机数,并对生成的随机数的统计性质进行了计算。

        经过多次实验,生成的随机数序列统计特征随着样本数量的增加逐渐逼近均值,在百万数量级的数据测试中,均值和方差都向下接近于设定参数。

import random

def randomTestCase(mu, sigma, loop=5000000):
    data = []
    for _ in range(loop):
        data.append(random.normalvariate(mu=mu, sigma=sigma))
    # print("data:\n", data)

    s = float(sum(data))/len(data)
    print("均值:{}".format(s))
    data1 = [(x - s) ** 2 for x in data]
    var = sum(data1) / len(data)
    print("方差:{}".format(var))


if __name__ == "__main__":
    randomTestCase(mu=50, sigma=4)

  

 

  • 3
    点赞
  • 31
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
正态分布随机数产生方法有多种,其中一种常用的方法是Box-Muller方法。该方法的推导过程较为复杂,但得到的结果却是很令人满意的。使用Box-Muller方法,我们可以通过抽取两个相互独立的均匀分布的随机数来得到正态分布随机数。具体步骤如下: 1. 随机抽取两个从均匀分布\[0,1\]中抽取的数字u和v。 2. 计算z = sqrt(-2 * log(u)) * cos(2 * pi * v)。 3. 最后,将z乘以标准差sigma,并加上均值mu,即可得到符合正态分布随机数。 以下是一个使用Python实现的Box-Muller方法的示例代码: ```python import numpy as np def boxmullersampling(mu=0, sigma=1, size=1): u = np.random.uniform(size=size) v = np.random.uniform(size=size) z = np.sqrt(-2 * np.log(u)) * np.cos(2 * np.pi * v) return mu + z * sigma ``` 另外,根据中心极限定理,我们也可以通过生成n个独立同分布的均匀分布来近似生成正态分布。具体步骤如下: 1. 生成n个独立同分布的均匀分布随机数。 2. 将这些随机数相加,得到一个和。 3. 重复上述步骤多次,得到多个和。 4. 最后,将这些和进行归一化处理,即可得到符合正态分布随机数。 以下是一个使用MATLAB实现的生成正态分布的示例代码: ```matlab clear all n = unifrnd(0, 1, 200000, 1); N = 50; w = zeros(1, 4000); w(1) = 0; for t = 1:4000 for j = 1:N w(t) = w(t) + n((j-1)*4000 + t); end end figure(1); hist(w, 400); ``` 这些方法都可以用来生成符合正态分布随机数。具体选择哪种方法取决于你的需求和编程环境。 #### 引用[.reference_title] - *1* *2* [正态分布随机数生成算法](https://blog.csdn.net/u010899985/article/details/80784810)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^control_2,239^v3^insert_chatgpt"}} ] [.reference_item] - *3* [正态分布随机数产生方法](https://blog.csdn.net/fengying2016/article/details/80601991)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^control_2,239^v3^insert_chatgpt"}} ] [.reference_item] [ .reference_list ]

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值