Python使用numpy产生正态分布随机数的向量或矩阵操作示例
2018-08-22 10:25
1251 查看
本文实例讲述了Python使用numpy产生正态分布随机数的向量或矩阵操作。分享给大家供大家参考,具体如下:
简单来说,正态分布(Normal distribution)又名高斯分布(Gaussian distribution),是一个在数学、物理及工程等领域都非常重要的概率分布,在统计学的许多方面有着重大的影响力。一般的正态分布可以通过标准正态分布配合数学期望向量和协方差矩阵得到。如下代码,可以得到满足一维和二维正态分布的样本。
示例1(一维正态分布):
# coding=utf-8 ''' 作者:采石工 来源:知乎 ''' import numpy as np from numpy.linalg import cholesky import matplotlib.pyplot as plt sampleNo = 1000; # 一维正态分布 # 下面三种方式是等效的 mu = 3 sigma = 0.1 np.random.seed(0) s = np.random.normal(mu, sigma, sampleNo ) plt.subplot(141) plt.hist(s, 30, normed=True) np.random.seed(0) s = sigma * np.random.randn(sampleNo ) + mu plt.subplot(142) plt.hist(s, 30, normed=True) np.random.seed(0) s = sigma * np.random.standard_normal(sampleNo ) + mu plt.subplot(143) plt.hist(s, 30, normed=True) # 二维正态分布 mu = np.array([[1, 5]]) Sigma = np.array([[1, 0.5], [1.5, 3]]) R = cholesky(Sigma) s = np.dot(np.random.randn(sampleNo, 2), R) + mu plt.subplot(144) # 注意绘制的是散点图,而不是直方图 plt.plot(s[:,0],s[:,1],'+') plt.show()
运行结果:
示例2(正态分布):
#-*- coding:utf-8 -*- # Python实现正态分布 # 绘制正态分布概率密度函数 import numpy as np import matplotlib.pyplot as plt import math u = 0 # 均值μ u01 = -2 sig = math.sqrt(0.2) # 标准差δ x = np.linspace(u - 3*sig, u + 3*sig, 50) y_sig = np.exp(-(x - u) ** 2 /(2* sig **2))/(math.sqrt(2*math.pi)*sig) print x print "="*20 print y_sig plt.plot(x, y_sig, "r-", linewidth=2) plt.grid(True) plt.show()
运行结果:
更多关于Python相关内容感兴趣的读者可查看本站专题:《Python数据结构与算法教程》、《Python编码操作技巧总结》、《Python函数使用技巧总结》、《Python字符串操作技巧汇总》及《Python入门与进阶经典教程》
希望本文所述对大家Python程序设计有所帮助。
您可能感兴趣的文章:
相关文章推荐
- python numpy 使用笔记 矩阵操作
- python 正态分布随机数 numpy.random.randn 使用小技
- python 正态分布随机数 numpy.random.randn 使用小技
- Python使用numpy模块创建数组操作示例
- 两个使用Python脚本操作文件的小示例分享
- 关于用mapreduce做kmeans聚类以及python的numpy和list做矩阵、向量乘法的速度对比
- [Python]使用Numpy对矩阵进行转置
- 两个使用Python脚本操作文件的小示例分享
- 使用python删除nginx缓存文件示例(python文件操作)
- python使用cx_Oracle操作oracle示例
- Python-Numpy(3)矩阵基本操作
- Python使用当前时间、随机数产生一个唯一数字的方法
- matlab函数randn:产生正态分布的随机数或矩阵的函数
- Python列表list解析操作示例【整数操作、字符操作、矩阵操作】
- python的Numpy之矩阵操作
- python列表操作使用示例分享
- python文件读写并使用mysql批量插入示例分享(python操作mysql)
- python使用numpy向量运算注意点
- Python使用scipy和numpy操作处理图像
- Python中矩阵库Numpy基本操作