Python隨機(jī)生成均勻分布在單位圓內(nèi)的點(diǎn)代碼示例
Python有一隨機(jī)函數(shù)可以產(chǎn)生[0,1)區(qū)間內(nèi)的隨機(jī)數(shù)，但是如果我們想生成隨機(jī)分布在單位圓上的，那么我們可以首先生成隨機(jī)分布在單位圓邊上的點(diǎn)，然后隨機(jī)調(diào)整每個(gè)點(diǎn)距離原點(diǎn)的距離，但是我們發(fā)現(xiàn)這個(gè)距離不是均勻分布于[0,1]的，而是與扇形的面積相關(guān)的
我們使用另外的隨機(jī)函數(shù)生成從[0,1)的隨機(jī)數(shù)r，我們發(fā)現(xiàn)r<s0的概率為s0，顯而易見，如果r為0，那么對應(yīng)的距離應(yīng)該為0，如果是1，對應(yīng)的距離自然也應(yīng)該是1，假設(shè)我們產(chǎn)生了m個(gè)隨機(jī)數(shù)，那么小于s0的隨機(jī)數(shù)應(yīng)該為s0*m左右，而且這些應(yīng)該對應(yīng)于扇形面積的s0倍處即圖2的小扇形區(qū)域，落在這一區(qū)域的點(diǎn)應(yīng)該為s0*m，此時(shí)扇形邊長為s0^0.5,因此s0對應(yīng)的距離應(yīng)該為s0^0.5，因此我們得到的映射函數(shù)為y=x^0.5(圖1)
圖1
圖2
因此我們對于每個(gè)頂點(diǎn)的邊長便是產(chǎn)生隨機(jī)數(shù)的算術(shù)平方根的大小
附代碼如下：
# -*- 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
總結(jié)
以上就是本文關(guān)于Python隨機(jī)生成均勻分布在單位圓內(nèi)的點(diǎn)代碼示例的全部內(nèi)容，希望對大家有所幫助。感興趣的朋友可以繼續(xù)參閱本站：python實(shí)現(xiàn)圖片處理和特征提取詳解、python中實(shí)現(xiàn)k-means聚類算法詳解、Python內(nèi)存管理方式和垃圾回收算法解析等，有什么問題可以隨時(shí)留言，小編會及時(shí)回復(fù)大家的。感謝朋友們對本站的支持！
相關(guān)文章
在本篇文章里小編給大家整理的是關(guān)于python實(shí)例化對象的具體方法，有興趣的朋友們可以學(xué)習(xí)下。2020-06-06
在 Python 中,邊遍歷邊刪除列表元素通常是一個(gè)不推薦的操作,因?yàn)樗鼤淖兞斜淼拇笮?可能會導(dǎo)致一些意料之外的行為,例如,元素被刪除后,列表的索引會發(fā)生變化,可能導(dǎo)致漏掉某些元素或者遍歷到錯誤的位置,所以本文介紹了Python邊遍歷邊刪除列表元素的幾種方法2024-12-12
今天小編就為大家分享一篇Anaconda 離線安裝 python 包的操作方法，具有很好的參考價(jià)值，希望對大家有所幫助。一起跟隨小編過來看看吧2018-06-06
這篇文章主要介紹了Python 數(shù)字轉(zhuǎn)化成列表，主要以代碼實(shí)現(xiàn)了將輸入的數(shù)字轉(zhuǎn)化成一個(gè)列表，輸入數(shù)字中的每一位按照從左到右的順序成為列表中的一項(xiàng)。,需要的朋友可以參考下2021-11-11
這篇文章主要介紹了python基于windows平臺鎖定鍵盤輸入的方法,涉及Python中調(diào)用user32.dll中BlockInput方法的技巧,具有一定參考借鑒價(jià)值,需要的朋友可以參考下2015-03-03
pip更新問題的解決:'python -m pip install -
這篇文章主要介紹了pip更新問題的解決:'python -m pip install --upgrade pip' 報(bào)錯問題，本文給大家介紹的非常詳細(xì)，對大家的學(xué)習(xí)或工作具有一定的參考借鑒價(jià)值,需要的朋友可以參考下2023-04-04
本文主要介紹了Python隨機(jī)生成迷宮游戲的代碼示例，Python可以通過生成二維數(shù)組模擬迷宮的結(jié)構(gòu)，使用深度優(yōu)先搜索和廣度優(yōu)先搜索等算法找到通路，文中通過示例代碼介紹的非常詳細(xì)，需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧2023-07-07
這篇文章主要介紹了對于Python裝飾器使用的一些建議,裝飾器是Python學(xué)習(xí)進(jìn)階中的重要知識,需要的朋友可以參考下2015-06-06
這篇文章主要介紹了Python yield 使用淺析,本文給出了多個(gè)使用實(shí)例來分析yield的使用方法,需要的朋友可以參考下2015-05-05
配置文件管理已經(jīng)成為了一個(gè)不可或缺的環(huán)節(jié)。它們包含了程序運(yùn)行所需的各種參數(shù)，本文主要介紹了Python中使用YAML文件進(jìn)行配置文件管理，感興趣的可以了解一下2023-12-12