Box–Muller变换最初由 George E. P. Box 与 Mervin E. Muller 在1958年提出

George E. P. Box 是统计学的一代大师,统计学中的很多名词术语都以他的名字命名

统计学中的名言 “All models are wrong, but some are useful”(所有模型都是错的,但其中一些是有用的)出自 Box 之口

程序检验

import random
import numpy as np
import matplotlib.pyplot as plt
%matplotlib inline
import seaborn as snsu1 = np.random.random(10000)
u2 = np.random.random(10000)
z1 = np.sqrt(-2*np.log(u1))*np.cos(2*np.pi*u2)
z2 = np.sqrt(-2*np.log(u1))*np.sin(2*np.pi*u2)

核密度估计

sns.distplot(z1,        bins = 100, hist = True,kde = True, norm_hist = False, rug = True, vertical = False,color = 'b', label = 'distplot', axlabel = 'x'
)
plt.legend()

正态检验

import scipy.statsscipy.stats.normaltest(z1)
"""
NormaltestResult(statistic=0.6940204433592375, pvalue=0.7067981035570555)
"""

两变量独立

import pandas as pd
df = pd.DataFrame({'A':z1, 'B':z2})#两个维度数据生成曲线密度图,以颜色作为密度的衰减显示
sns.kdeplot(df['A'],df['B'],cbar = True,          #是否显示颜色的图例shade = True,         #是否填充cmap = 'Reds',        #设置调色盘shade_lowest = False, #最外围颜色是否显示n_levels = 10        #曲线个数(如果非常多,则会越平滑))#注意设置x,y轴
plt.scatter(df['A'],df['B'],s = 5, alpha = 0.5, color = 'k')
sns.rugplot(df['A'], color = 'g', axis = 'x', alpha = 0.5)
sns.rugplot(df['B'], color = 'r', axis = 'y', alpha = 0.5)

Box-Muller 变换相关推荐

  1. R语言构建回归模型并进行模型诊断(线性关系不满足时)、进行变量变换(Transforming variables)、使用car包中的boxTidwell函数对预测变量进行Box–Tidwell变换

    R语言构建回归模型并进行模型诊断(线性关系不满足时).进行变量变换(Transforming variables).使用car包中的boxTidwell函数对预测变量进行Box–Tidwell变换 目 ...

  2. 计算机视觉资源汇总 - Part I

    转载自:http://blog.csdn.net/zhubenfulovepoem/article/details/7191794 以下是computer vision:algorithm and a ...

  3. opencv图像处理学习(二十六)——噪声

    1.噪声量化 图像噪声是图像在获取或传输的过程中受到随机信号干扰,妨碍人们对图像处理及分析处理的信号.很多时候将图像噪声看作多位随机过程,因而描述噪声的方法完全可以借用随机过程的描述,即使用其概率分布 ...

  4. 蒙特卡洛采样与Gibbs采样

    蒙特卡洛采样与Gibbs采样 采样 采样问题指的是给定一个特定的概率分布p(z),得到一批符合这个概率分布的样本点. 采样的方法有很多,MCMC是其中的一类方法,意思是利用Mento Carlo和Ma ...

  5. 23.代码简单实现模拟噪声(图像噪声/一、二阶矩/功率谱密度/at函数/rand函数)-- OpenCV从零开始到图像(人脸 + 物体)识别系列

    本文作者:小嗷 微信公众号:aoxiaoji 吹比QQ群:736854977 简书链接:https://www.jianshu.com/u/45da1fbce7d0 本文你会找到以下问题的答案: 图像 ...

  6. 史上最全采样方法详细解读与代码实现

    项目github地址:bitcarmanlee easy-algorithm-interview-and-practice 欢迎大家star,留言,一起学习进步 1.什么是采样 在信号系统.数字信号处 ...

  7. Computer Vision: Algorithms and Application书籍章节介绍及其源码

    ComputerVision: Algorithms and Applications Richard Szeliski 本书网址:书的最好附录中,我总结了一些对学生,教授和研究者有用的附加材料.这本 ...

  8. 随机数生成方法及其应用

    目录 一.前言 二.伪随机数发生器介绍 三.均匀随机数发生器 四.任意概率分布的随机数发生器 五.产生特定要求的随机数 六.参考资料 一.前言 最近公司的一个项目,需要对一个功能模块进行性能验证,具体 ...

  9. 根据概率密度函数生成随机数的代码

    我这里并不是要讲"伪随机"."真随机"这样的问题,而是关于如何生成服从某个概率分布的随机数(或者说 sample)的问题.比如,你想要从一个服从正态分布的随机变 ...

最新文章

  1. 微信小程序商城demo
  2. 循环数组函数c语言,C语言练习题2(分支结构循环结构数组函数2009-2012二级真题)..doc...
  3. Centos 6安装python3.5
  4. c语言函数打印零星,C语言零星笔记--(const、typedef)
  5. 另我模糊的URL 重写
  6. ustc小道消息20211223
  7. 给文物安全信息化支个招
  8. 我做的小学二年级公开课《生命,生态和安全》的教案
  9. jdk8永久代从方法区移除的验证
  10. Android 5.0新特性
  11. 释放数据价值:DAYU数据运营新能力解读
  12. ubuntu虚拟机进程被杀死_Linux进程被杀掉(OOM killer),查看系统日志
  13. 2 --> PCIE 协议栈的框架详解
  14. win10安装Visual Studio 2019失败
  15. 抖音上好看的小姐姐,Python给你都下载了
  16. Android 字体颜色渐变效果 Span实现
  17. 什么百度霸屏?百度霸屏是什么意思?
  18. 计算机课导入语,信息技术课的提问和导入技巧
  19. 项目实训(十六)——总结
  20. SeEnableDelegationPrivilege权限的利用

热门文章

  1. 初级入门华为交换机学习
  2. JS读取字幕srt文件内容
  3. 自学实前后端践项目4 MMall商城 1
  4. 如何使用iPhone生存僵尸启示录
  5. 二进制转十进制在c51单片机中的问题
  6. 出货量第一的三星为何在手机AI芯片竞争中落后华为和苹果?
  7. 如何实现身份证输入框的输入效果
  8. 如何判断某一个IP地址是否被占用
  9. 导出excel,sheetName名字乱码
  10. 数据中的知识发现KDD和机器学习概念