Loading [MathJax]/jax/element/mml/optable/BasicLatin.js

2 votos

Un conejo intenta escapar de su depredador zigzagueando de la siguiente manera.

Este problema trata de un conejo que intenta escapar de su depredador zigzagueando. Inicialmente se fija un parámetro α(0,π2) (en radianes), y en cada segundo n1 escoge Un Uniforme (0,α) ángulo. Entre los tiempos n1 y n (segundos) corre a velocidad unitaria en dirección Un en relación con el Este, entonces escoge Un+1 y corre en esa dirección con respecto al Este, etc. Por ejemplo, si Un=π4 , entonces a partir del tiempo n1 al tiempo n exactamente corre hacia el noreste. Por lo tanto, la posición del conejo al final de la n El segundo es (Sn,Tn)=(ni=1cosUi,ni=1sinUi). Supongamos que el Un son independientes.

  • Demuestre que el conejo tiene una velocidad asintótica en el sentido de que su distancia Dn=S2n+T2n desde su punto de partida en el momento n dividido por n converge a.s. a un límite.
  • Calcula esta velocidad asintótica limnDnn en función del parámetro α . ¿Qué pasa si α está cerca de 0 ?

Hasta ahora he hecho algunos progresos. Tenemos que Dnn=(Snn)2+(Tnn)2a.s.(EcosUi)2+(EsinUi)2. Así que ahora todo lo que tengo que hacer es calcular este límite. Así que en primer lugar EcosUi=α01αcosxdx=1α[sinx]α0=sinαα y EsinUi=α01αsinxdx=1α[cosx]α0=1α(1cosα). Así que ahora tenemos que lim Luego, a la hora de analizar el comportamiento de este límite en \alpha=0 Empecé a sospechar que había hecho algo mal en mi trabajo, porque habría pensado que esto sería un límite obvio de 1 ya que el conejo iría casi seguro en la misma dirección a velocidad unitaria. Sin embargo, cuando he intentado evaluar este límite, no he podido averiguar cómo, lo que me hace sospechar que he hecho algo mal, ya que esperaba que esta parte de la pregunta fuera fácil. Así que, ¿alguien puede ver en qué me he equivocado o ayudarme a calcular este límite? Gracias.

2voto

LoveTooNap29 Puntos 96

Su argumento y sus cálculos son correctos hasta su punto de preocupación, por lo que he comprobado. Aquí hay una simulación de D_n/n como n\to 9000 , para \alpha = \pi/4 . Converge al valor f(\pi/4) \approx 0.9744954 dado por su fórmula f(\alpha)=\sqrt{2-2\cos(\alpha)}/\alpha . enter image description here

Aquí está el código para producir esto en R,

n <- 9000
f <- function(a)  sqrt(2-2*cos(a))/a
alpha <- pi/4
u <- runif(n, 0, alpha)
x <- cumsum(cos(u))
y <- cumsum(sin(u))
d <- sqrt(x^2+y^2)
v <- d/(1:n)
plot(1:n, v, type = "l")
abline(h = f(alpha), col = "blue", lty = "dashed")
w <- c(exact = f(alpha), mc = sqrt(mean(cos(u))^2+mean(sin(u))^2))
print(w)

Por último, sobre el comportamiento como \alpha \to 0 . Afirmamos que \lim_{\alpha \to 0} \frac{\sqrt{2-2\cos(\alpha)}}{\alpha}=1. Ya que, si se conocen las series de Taylor/expansiones de Little-O, tenemos \cos(\alpha)=1-\alpha^2/2+o(\alpha^4) como \alpha \to 0 . Así, obtenemos \frac{\sqrt{2}}{\alpha} \sqrt{1-\left(1- \frac 12 \alpha^2+o(\alpha^4)\right)} =\sqrt{2} \sqrt{\frac 12 + o(\alpha^2)}\to 1, como \alpha \to 0 . Numéricamente, esto se comprueba también si elegimos algunos pequeños \alpha cerca de cero en el código anterior.

Por favor, comente para aclaraciones o correcciones.

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X