时间:2021-05-22
Python有一随机函数可以产生[0,1)区间内的随机数,但是如果我们想生成随机分布在单位圆上的,那么我们可以首先生成随机分布在单位圆边上的点,然后随机调整每个点距离原点的距离,但是我们发现这个距离不是均匀分布于[0,1]的,而是与扇形的面积相关的
我们使用另外的随机函数生成从[0,1)的随机数r,我们发现r<s0的概率为s0,显而易见,如果r为0,那么对应的距离应该为0,如果是1,对应的距离自然也应该是1,假设我们产生了m个随机数,那么小于s0的随机数应该为s0*m左右,而且这些应该对应于扇形面积的s0倍处即图2的小扇形区域,落在这一区域的点应该为s0*m,此时扇形边长为s0^0.5,因此s0对应的距离应该为s0^0.5,因此我们得到的映射函数为y=x^0.5(图1)
图1
图2
因此我们对于每个顶点的边长便是产生随机数的算术平方根的大小
附代码如下:
# -*- coding:utf-8 -*-import numpy as npimport matplotlib.pyplot as pltif __name__=='__main__': samples_num = 800 t = np.random.random(size=samples_num) * 2 * np.pi - np.pi x = np.cos(t) y = np.sin(t) i_set = np.arange(0,samples_num,1) for i in i_set: len = np.sqrt(np.random.random()) x[i] = x[i] * len y[i] = y[i] * len plt.figure(figsize=(10,10.1),dpi=125) plt.plot(x,y,'ro') _t = np.arange(0,7,0.1) _x = np.cos(_t) _y = np.sin(_t) plt.plot(_x,_y,'g-') plt.xlim(-1.1,1.1) plt.ylim(-1.1,1.1) plt.xlabel('x') plt.ylabel('y') plt.title('Random Scatter') plt.grid(True) plt.savefig('imag.png') plt.show()图3
总结
以上就是本文关于Python随机生成均匀分布在单位圆内的点代码示例的全部内容,希望对大家有所帮助。感兴趣的朋友可以继续参阅本站:python实现图片处理和特征提取详解、python中实现k-means聚类算法详解、Python内存管理方式和垃圾回收算法解析等,有什么问题可以随时留言,小编会及时回复大家的。感谢朋友们对本站的支持!
声明:本页内容来源网络,仅供用户参考;我单位不保证亦不表示资料全面及准确无误,也不保证亦不表示这些资料为最新信息,如因任何原因,本网内容或者用户因倚赖本网内容造成任何损失或损害,我单位将不会负任何法律责任。如涉及版权问题,请提交至online#300.cn邮箱联系删除。
1.均匀分布函数:runif(n,min=0,max=1),n表示生成的随机数数量,min表示均匀分布的下限,max表示均匀分布的上限,若省略参min、max,
今天学习了用python生成仿真数据的一些基本方法和技巧,写成博客和大家分享一下。本篇博客主要讲解如何从给定参数的的正态分布/均匀分布中生成随机数以及如何以给定
在用tensorflow构建神经网络的时候,有很多随机的因素,比如参数的随机初始化:正态分布随机变量tf.random_normal([m,n]),均匀分布的随
本文实例讲述了C++抽奖程序实现方法。分享给大家供大家参考。具体实现方法如下:一、intrand()可以生成从[0,65536)之间均匀分布的随机数。现要求实现
前言大家应该都知道,在Oracle数据库中,CBO会默认目标列的数据在其最小值low_value和最大值high_value之间均匀分布,并按照均匀分布原则,来