Python随机生成均匀分布在单位圆内的点代码示例
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 np
import matplotlib.pyplot as plt
if __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内存管理方式和垃圾回收算法解析等，有什么问题可以随时留言，小编会及时回复大家的。感谢朋友们对本站的支持！
相关文章
网站被入侵，担心被挂马，因此就想自己写个脚本来查找那些被挂马的文件2014-06-06
这篇文章主要介绍了解密Python中的描述符（descriptor）,本文详细讲解了描述符（descriptor）的作用、访问描述符、对描述符赋值、删除描述符等内容,需要的朋友可以参考下2015-06-06
这篇文章主要介绍了python之拟合的实现，文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值，需要的朋友们下面随着小编来一起学习学习吧2019-07-07
这篇文章主要给大家介绍了关于numpy版本兼容问题解决方案的相关资料,文中通过图文介绍的非常详细,对大家学习或者使用numpy具有一定的参考借鉴价值,需要的朋友可以参考下2023-10-10
这篇文章主要介绍了如何使用Python脚本实现文件拷贝,文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下2019-11-11
这篇文章主要介绍了Python从文件中读取数据的方法步骤，文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值，需要的朋友们下面随着小编来一起学习学习吧2020-11-11
上一篇我们已经知道怎么简单使用selenium了,那么我们就从这篇博客来动手爬取网站吧,文中有非常详细的代码示例,需要的朋友可以参考下2021-06-06
今天小编就为大家分享一篇Python3 串口接收与发送16进制数据包的实例，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧2019-06-06
这篇文章主要介绍了Python脚本判断 Linux 是否运行在虚拟机上,本文讲解了判断 OpenVZ/Xen PV/UML、判断 VMware/Xen HVM/KVM、判断 VirtualBox/Virtual PC等方法,需要的朋友可以参考下2015-04-04
这篇文章主要是和大家分享一个Python实现PDF文字识别与提取并写入 CSV文件的脚本。文中的示例代码讲解详细，感兴趣的小伙伴可以了解一下2022-03-03
最新评论