Создание равномерно распределенных точек по кругу

Question

'' 'Я пытаюсь создать точки в круге, которые должны быть равномерно распределены, но получаю странный узор. Если я увеличиваю радиус R до очень большого значения, распределение кажется нормальным, но при меньших значениях R оно генерирует спирали. Есть предложения по улучшению кода? '' '

from numpy.random import uniform
#from scipy.stats import uniform
import matplotlib.pyplot as plt
import numpy as np
import math  

R = 5
# Generate uniformly distributed random numbers.
rand_num = []
for i in range(30000):
    rand_num.append(np.random.uniform(0,1))

# Use these generated numbers to obtain the CDF of the radius which is the true radius i.e. r = R*sqrt(random()).
radius = []
for n,data in enumerate(rand_num):
    radius.append(R*math.sqrt(data))

# Generate the angle using the same uniformly distributed random numbers.
theta2 = []
for n, k in enumerate(radius):
    theta2.append(2*math.pi*radius[n])

# Calculate the corresponding x-coordinate.
x = []
for j,v in enumerate(radius):
    x.append(radius[j]*math.cos(theta2[j]))
x = np.array(x)   

# Calculate the correspoding y-coordinate.
y = []    
for j,v in enumerate(radius):
    y.append(radius[j]*math.sin(theta2[j]))
y = np.array(y)

# Finally plot the coordinates.
plt.figure(figsize=(10,10))
plt.scatter(x, y, marker='o')

Answer 1

Да, код должен дать вам спираль, поскольку и theta, и radius пропорциональны rand_num. Вместо этого вы должны генерировать theta и radius независимо. Также используйте векторизованные операторы numpy вместо математических

R = 5
num_points = 10000

np.random.seed(1)
theta = np.random.uniform(0,2*np.pi, num_points)
radius = np.random.uniform(0,R, num_points) ** 0.5

x = radius * np.cos(theta)
y = radius * np.sin(theta)

# visualize the points:
plt.scatter(x,y, s=1)

. Вывод: