5 votos

integral definida

En la relatividad general, null geodesics (en el unbounded caso) puede escribirse en la siguiente forma : $$\frac{d\varphi}{dr}=\frac{1}{r^2\sqrt{\frac{R_0-R_S}{R_0^3}-\frac{1}{r^2}\left(1-\frac{R_{s}}{r}\right)}}$$ con:

  • $\left(r, \varphi\right)$ las coordenadas polares de los fotones
  • $R_S$ el Schwarschild radio de la central de objeto ($R_S\in\mathbb{R^{+}_{*}}$)
  • $R_0$ de la distancia más cercana de enfoque ($R_0 > \frac{3\sqrt{3}}{2}R_S$)

En consecuencia, para calcular la trayectoria exacta de hasta un radio de $R$ ( $R > R_0$ ), se puede evaluar: $$\varphi\left(R\right)=\int_{R_0}^{R}\frac{dr}{r^2\sqrt{\frac{R_0-R_S}{R_0^3}-\frac{1}{r^2}\left(1-\frac{R_{s}}{r}\right)}}$$

Y ahora viene mi pregunta: ¿existe una fórmula analítica (en términos de funciones especiales, por ejemplo) correspondiente a esta integral ? (Me gustaría para calcular esta integral numéricamente y una expresión en términos de funciones especiales ayudaría mucho).

1voto

Joe Gauterin Puntos 9526

La integral se puede expresar en términos de la inversa de Weierstrass de la elíptica funciones y, por tanto, en sentido contrario a la Jacobi funciones elípticas. La otra respuesta ha mencionado el incompleta integral elíptica de primera especie. Que la integral es esencialmente uno de los inversos Jacobi elíptica funciones en el disfraz.

Voy a darle a la expresión en términos de la inversa de Weierstrass elíptica funciones sólo. A diferencia de la expresión de Jacobi elíptica funciones de los/las integrales elípticas, esta parte es mucho más fácil de derivar.

Vamos $\alpha = \frac{R_s}{R_0} = \beta + \frac13$, $u = \frac{R_s}{r} = p + \frac13$, la integral se puede escribir como

$$\varphi(R) = \int_{\frac{R_s}{R}}^\alpha \frac{du}{\sqrt{\alpha^2(1-\alpha) - u^2(1-u)}} = \int_{\frac{R_s}{R}-\frac13}^{\beta}\frac{dp}{\sqrt{p^3 - \frac{p}{3} - ( \beta^3 - \frac{\beta}{3} )}} $$ Deje $\wp(z)$ ser el de Weierstrass elíptica función asociada con la educación a distancia

$$\wp'(z)^2 = 4\wp(z)^3 - g_2 \wp(z) - g_3\quad\text{ donde }\quad \begin{cases}g_2 &= \frac43\\g_3 &= 4(\beta^3 - \frac{\beta}{3})\end{casos}$$

and substitute $p$ by $\wp(z)$, we have

$$\varphi(R) = 2 \int_{\frac{R_s}{R}-\frac13}^\beta \frac{d\wp (z)}{\sqrt{4\wp(z)^3 - g_2 \wp(z) - g_3}}$$

If one change variable to $z$, the factor $\frac{d\wp}{\sqrt{4\wp^3 - g_2\wp - g_3}}$ in above integrand reduces to either $+dz$ or $-dz$.

Let $\omega$ and $\omega'$ be the half periods of $\wp(z)$ with $\omega \en (0,\infty)$ and $\Im\omega' > 0$. Let $\omega_c = \omega + \omega'$. We will adopt the convention as $\rho$ decreases from $+\infty$ to $-\infty$ along the real axis, $\wp^{-1}(\rho)$ will move along the path

$$0 \to \omega \to \omega_c = \omega + \omega' \to \omega' \to 0$$

For the parameter range $R > R_0 > \frac{3\sqrt{3}}{2} R_s$ we are interested in, $\wp^{-1}(\beta) = \omega_c$ and in general $\wp^{-1}(p)$ will lie on the line segment joining $\omega'$ and $\omega + \omega'$. On this line segment, above factor reduces to $dz$ y llegamos

$$\varphi(R) = 2 \int_{\wp^{-1}\left(\frac{R_s}{R}-\frac13\right)}^{\wp^{-1}(\beta)} dz = 2 \left[ \omega_c - \wp^{-1}\left(\frac{R_s}{R} - \frac13\right) \right] $$ Computacionalmente, esto no es conveniente para el uso. Esto es debido a que para el rango de $r$ nos importa, $\wp^{-1}(p)= \wp^{-1}\left(\frac{R_s}{r} - \frac13\right)$ es un número complejo! Utilizando la siguiente formula de sumar de a $\wp(z)$:

$$\wp(z+\omega_c) + \wp(z) + \wp(\omega_c) = \frac14 \left(\frac{\wp'(z) - \wp'(\omega_c)}{\wp(z) - \wp(\omega_c)}\right)^2$$ Nos encontramos

$$ \wp(\omega_c-z) = \wp(z+\omega_c) = \frac{\left(p^3-\frac{p}{3}\right) - \left(\beta^3 - \frac{\beta}{3}\right)}{( (p-\beta)^2} - (p + \beta) = \beta + \frac{9\beta^2-1}{3(p-\beta)} $$ y por lo tanto $$\varphi(R) = 2\wp^{-1}\left( \beta + \frac{9\beta^2-1}{3(p-\beta)} \right) = 2\wp^{-1}\left(\alpha\frac13 + \frac{\alpha(3\alpha-2)}{\frac{R_s}{R} - \alpha}\right) \etiqueta{*1}$$ Para demostrar cómo este trabajo, vamos a considerar el caso de $(R_s, R_0, R ) = (1,3,4)$ como un ejemplo. Tenemos $\alpha = \frac13$ $\varphi(R)$ se convierte en

$$\begin{align} \varphi(R) = & \int_3^4 \frac{dr}{r^2\sqrt{\frac{2}{27} - \frac{1}{r^2}\left(1 - \frac{1}{r}\right)}}\\ \approx & 1.00210163826249415497545279509980330387176483490377 \end{align} $$ El valor numérico se evalúa utilizando el siguiente comando en WA (Wolfram Alpha).

$\quad(a)\quad\quad$ N[Integrate[1/(r^2*Sqrt[2/27-1/r^2*(1-1/r)]),{r,3,4}],50]

Para este ejemplo, $(g_2, g_3) = (4/3,0)$. RHS de $(*1)$ hace $\wp^{-1}(4)$ y una vez más, podemos evaluar en WA usar otro comando

$\quad(b)\quad\quad$ N[-2*InverseWeierstrassP[ 4, {4/3,0}],50]

y obtener $$2\wp^{-1}(4) \approx 1.0021016382624941549754898784512602241681634415061$$

To those with the sharp eyes, you will notice two things.

  1. the numerical output of these two commands are very close but not the same. The only thing I can said is WA label the output for $(a)$ como aproximación decimal, mientras que el salida para $(b)$ como resultado.

  2. Hay un misterioso signo menos en el comando $(b)$. Resulta WA utiliza una convención diferente para $\wp^{-1}(p)$ de mí. Para grandes y posible $p$, InverseWeierstrass[p,{g_2,g_3}] devuelve un negativo en vez de un número positivo y de ahí la necesidad de que extra signo menos.

Actualización

Acerca de la expresión de todo en términos de la inversa de Jacobi funciones elípticas. Finalmente me desenterrar mis libros de referencia y encontrar algo útil. Deje $e_1, e_2, e_3$ serían las 3 raíces de los polinomios de $4p^3 - g_2 p - g_3 = 0$. Para la concreción, nos deja elegir

$$e_1 = \wp(\omega),\quad e_2 = \wp(\omega_c)\quad\text{ and }\quad e_3 = \wp(\omega')$$

We have

$$\wp'(z)^2 = 4 \wp(z)^3 - g_2 \wp(z) - g_3 = 4(\wp(z) - e_1)(\wp(z) - e_2)(\wp(z) - e_3)$$

When $0 < \alfa < \frac23$, uno puede mostrar que estos 3 raíces son reales y distintas. Además, satisfacen $e_1 > e_2 = \beta > e_3$. Si uno define una función $\text{sn}(\cdot)$ por la siguiente relación

$$\wp(z) = e_3 + \frac{e_1-e_3}{\text{sn}(z\sqrt{e_1-e_3})^2} \quad\ffi\quad \text{sn}(z) = \sqrt{\frac{e_1 - e_3}{\wp\left(\frac{z}{\sqrt{e_1-e_3}}\right) - e_3}}$$

One can show that $\text{sn}(z)$ satisfy the ODE for the Legendre form of Jacobi elliptic sine function:

$$\text{sn}'(z)^2 = (1 - \text{sn}(z)^2)(1 - m\,\text{sn}(z)^2)\quad\text{ where }\quad m = \frac{e_2 - e_3}{e_1-e_3}$$

Uno puede utilizar esta relación para convertir la expresión a partir de la inversa de Weierstrass elíptica funciones a la inversa Jacobi funciones elípticas.

0voto

Integrals Puntos 2841

No hay una analítica muy cerrados forma a este, pero todavía podemos intentar aprender un poco más sobre él. Sin embargo, es relativa a las integrales Elípticas de primera especie. Nota, Jacobi Elíptica funciones son inversas para las integrales elípticas.

Si desea calcular numéricamente aún así, hágamelo saber, así lo hice y tienen códigos en mathematica del uso de NIntegrate. Yo también se calcula utilizando sólo Integrar en mathematica como bien y tienen este código.

Espero que esta ayuda, la escritura de la integral tenemos $$\varphi\left(R\right)=\int_{R_0}^{R}\frac{dr}{r^2\sqrt{\frac{R_0-R_S}{R_0^3}-\frac{1}{r^2}\left(1-\frac{R_{s}}{r}\right)}}$$ Definimos $a\equiv (R_0-R_s)/R^3_0$, $b\equiv R_s$ y escribir $$ \varphi\left(R\right)=\int_{R_0}^{R}\frac{dr}{r^2\sqrt {- \frac{1}{r^2}\left(1-\frac{b}{r}\right)}}$$

Ahora nos factor de algunas de nuestras constantes en virtud de la raíz cuadrada y hacer el álgebra para obtener $$ \varphi(r)=a^{-1/2}\int_{R_0}^R \frac{dr}{r^2\sqrt{1-\frac{1}{r^2}\left(a^{-1}-\frac{b}{ar}\right)}}=a^{-1/2}\int_{R_0}^R \frac{dr}{r \sqrt{r^2+\frac{b}{ar}-a^{-1}}}. $$ Ahora podemos analizar brevemente la integral indefinida de este $$ \varphi(r)=a^{-1/2}\int \frac{dr}{r \sqrt{r^2+\frac{b}{ar}-a^{-1}}} $$ Esta integral no tiene una forma cerrada, sin embargo, puede ser escrito en un no tan bonito Integral Elíptica de Primera Especie (Elíptica-F) de la Función. Así que para responder a su pregunta, es expresable en términos de integrales elípticas. Es algo parecido a $I\propto r F(\arcsin(...)) $. Mathematica no podía calcular mucho más a menos que desee utilizar NIntegrate.

Las integrales elípticas de primera clase en Jacobi de la forma está dada por $$ F(x,k)=\int_0^x \frac{dt}{\sqrt{(1-t^2)(1-k^2t^2)}}, $$ donde la forma usual, usted puede saber que es $$ F(\zeta,k)=\int_0^\zeta \frac{d\theta}{\sqrt{1-k^2\sin^2\theta}}, $$ están relacionados con un cambio de variables $t=\sin \theta$ y k es la relativa a la excentricidad de la órbita. Lo siento, no se puede encontrar una forma cerrada para su post, pero espero elíptica F funciones va a hacer. Definitivamente, esta es también la razón por la que no ponga la solución en ningún libro, porque no es tan limpio.

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