¿Por qué dibujar círculos repetidamente alrededor de un punto crea un patrón de "anillo"?

Definir un punto PAG 0 Al origen. Defina una "iteración" poniendo un círculo de radio uno alrededor PAG y eligiendo un punto uniformemente aleatorio en el círculo. Deja que ese nuevo punto se convierta PAG 1 . Luego repita ese proceso con PAG 1 , digamos 15 veces. Luego repite los pasos antes de un número grande, digamos 250,000 veces (ver esto aquí ). ¿Por qué la imagen proporcionada produce anillos de color? (Por anillos quiero decir que hay áreas de color)

El color se define a partir del número de iteraciones para llegar a ese punto, y cuanto más oscuro es el color, más iteraciones se necesitan. Los puntos que ocurren varias veces tienen su color promediado, y el negro significa que ningún punto llega a esa ubicación. El círculo blanco es la primera iteración. (Además, los puntos fuera del círculo NO impiden la iteración, simplemente no se dibujan)

Además, solo algunas preguntas curiosas:

  1. ¿Hay alguna manera de calcular la distancia promedio para norte iteraciones (sin simularlo)?
  2. ¿Hay alguna función que pueda "suavizar" cualquiera de las imágenes (no tiene que ser perfecta)?
  3. ¿Qué pasaría si pudieras elegir un radio aleatorio? (250,000 simulaciones aquí , 35px es un radio)
Este es un ejemplo de un paseo aleatorio.
Gracias, agregué esa etiqueta.
¿No es el segundo párrafo de su pregunta una respuesta perfectamente buena a la pregunta que hace en el primer párrafo?
Más o menos, pero hubiera esperado un degradado y no un anillo.
(Había generalizado demasiado ansiosamente en mi respuesta; ahora actualicé con un análisis más cuidadoso. Muestra que en realidad solo obtenemos bordes afilados en el radio 1, 2, 3, 4; no es solo que los límites en radios más grandes son más débiles, también son más suaves ).

Respuestas (1)

Primero considere solo los segundos puntos en cada recorrido aleatorio. Su dirección desde el centro está (por simetría) uniformemente distribuida, pero la distancia D desde el centro se distribuye como 2 | pecado ( X / 2 ) | , dónde X es la variable aleatoria que denota la diferencia de ángulo entre los dos saltos. Esta distribución de probabilidad no cae suavemente; de ​​hecho, la densidad de probabilidad a distancia D < 2 debería ser algo como

2 π d d D arcsen D 2 = 2 π 4 D 2
que aumenta hacia el infinito cuando D enfoques 2 , pero luego de repente cae a cero probabilidad para D > 2 .

Esto crea el borde visible a distancia. 2 en tu parcela.

La densidad de probabilidad por área debe tener en cuenta además que hay más área a una distancia de 1,8 a 1,9 que a una distancia de 0,1 a 0,2, por lo que se convierte en 1 π 2 D 4 D 2 . Esto va al infinito tanto en D 0 y D 2 , y de hecho vemos un pico en la densidad justo en el centro.


Para más de dos saltos, las cosas se vuelven menos dramáticas y también más complejas de calcular. Consideraré solo el comportamiento asintótico para puntos finales muy cercanos a las distancias máximas.

Para norte + 1 salta, deja variables aleatorias X 1 , X 2 , , X norte ser los ángulos entre cada uno de los saltos no primeros y el primer salto. Estos son entonces independientes y cada uno distribuido uniformemente en [ π , π ] . Como solo nos interesan los casos en los que terminamos muy cerca de la distancia máxima, podemos aproximar la distancia total como

D 1 + porque X 1 + + porque X norte norte + 1 X 1 2 + + X norte 2 2
Así, la probabilidad de aterrizar más cerca del máximo que ε es el volumen de un norte -bola de radio 2 ε , dividido por ( 2 π ) norte :

PAG ( D > norte + 1 ε ) ε norte / 2 ( 2 π ) norte / 2 Γ ( norte / 2 + 1 )

Esta probabilidad se distribuye sobre un anillo de área 2 π ( norte + 1 ) ε , por lo que podemos estimar la densidad de probabilidad por unidad de área justo dentro del círculo máximo como:

ε norte / 2 1 ( norte + 1 ) ( 2 π ) norte / 2 + 1 Γ ( norte / 2 + 1 ) como  ε 0

El exponente de ε es norte / 2 1 , por lo que el comportamiento es cualitativamente diferente entre diferentes norte .

  • Para dos saltos ( norte = 1 ), como hemos visto anteriormente, la densidad de puntos llega al infinito cerca del borde.

  • Para tres saltos ( norte = 2 ), la densidad se aproxima a un límite positivo finito cerca del borde, y luego de repente salta a 0 . Podemos ver un límite ligeramente más débil en el radio 3 en la trama

  • Para cuatro saltos ( norte = 3 ), la densidad queda como ε , por lo que en realidad alcanza 0 continuamente. Pero justo en el borde, todavía varía lo suficientemente rápido como para que el límite se pueda discernir en la trama.

  • En más de cuatro saltos, la densidad de probabilidad cae lo suficientemente suave hacia el límite que no se registra visiblemente como un borde .