72 votos

¿Cómo enderezar una parábola?

Considere la función $f(x)=a_0x^2$ para algunos $a_0\in \mathbb{R}^+$ . Toma $x_0\in\mathbb{R}^+$ de modo que la longitud de arco $L$ entre $(0,0)$ y $(x_0,f(x_0))$ es fijo. Dada una $a_1$ ¿cómo se encuentra el punto $(x_1,y_1)$ para que la longitud del arco sea la misma?

Esquemáticamente,

enter image description here

En otras palabras, estoy buscando una función $g:\mathbb{R}^3\to\mathbb{R}$ , $g(a_0,a_1,x_0)$ que toma un coeficiente cuadrático fijo inicial $a_0$ y punto y devuelve el punto correspondiente después de "enderezarlo" mediante el nuevo coeficiente $a_1$ manteniendo la longitud del arco con respecto a $(0,0)$ . Tenga en cuenta que el $y$ las coordenadas vienen dadas simplemente por $y_0=f(x_0)$ y $y_1=a_1x_1^2$ . ¿Alguna idea?

Mi enfoque: Sabiendo que el la longitud del arco viene dada por $$ L=\int_0^{x_0}\sqrt{1+(f'(x))^2}\,dx=\int_0^{x_0}\sqrt{1+(2a_0x)^2}\,dx $$ podemos utilizar la conservación de $L$ escribir $$ \int_0^{x_0}\sqrt{1+(2a_0x)^2}\,dx=\int_0^{x_1}\sqrt{1+(2a_1x)^2}\,dx $$ que resolvemos para $x_1$ . Esto funciona, pero no es muy rápido computacionalmente y sólo se puede hacer numéricamente (creo), ya que $$ \int_0^{x_1}\sqrt{1+(2a_1x)^2}\,dx=\frac{1}{4a_1}\left(2a_1x_1\sqrt{1+(a_1x_1)^2}+\arcsin{(2a_1x_1)}\right) $$ ¿Alguna idea sobre cómo hacerlo de forma más eficaz? ¿Quizás utilizando las líneas tangentes de la parábola?

En general , para longitudes de arco fijas, supongo que mi pregunta realmente es cuáles son las expresiones de las siguientes curvas rojas para longitudes de arco fijas:

enter image description here

Además, ¿podría determinarse para cualquier $f$ ?

Edita: Curiosamente, encontré este clip de 3Blue1Brown. El punto de origen no es fijo como en mi caso, pero me pregunto cómo se hizo la animación (no pude encontrar el vídeo original, sólo un clip, pero aquí está el enlace )

enter image description here

Para cualquier Mathematica entusiastas, también se está debatiendo una implementación computacional del efecto de enderezamiento aquí con algunas aplicaciones.

33voto

Ninad Munshi Puntos 801

Dicho de otro modo, lo que queremos son las curvas de nivel de la función

$$\frac{1}{2}f(x,y) = \int_0^x\sqrt{1+\frac{4y^2t^2}{x^4}}\:dt = \frac{1}{2}\int_0^2 \sqrt{x^2+y^2t^2}\:dt$$

que siempre será perpendicular al gradiente en ese punto

$$\nabla f = \int_0^2 dt\left(\frac{x}{\sqrt{x^2+y^2t^2}},\frac{yt^2}{\sqrt{x^2+y^2t^2}}\right)$$

Ahora es el momento de reintroducir de forma natural $a$ como parámetro para estas curvas. Por lo tanto lo que queremos es resolver la ecuación diferencial

$$x'(a) = \int_0^2 \frac{-axt^2}{\sqrt{1+a^2x^2t^2}}dt \hspace{20 pt} x(0) = L$$

donde sustituimos $y(a) = a\cdot x^2(a)$ Por lo tanto, si resolvemos un componente, obtendremos automáticamente el otro.


EDIT: Una investigación más profunda me ha llevado a algunas conclusiones interesantes. Parece que si $y=f_a(x)$ es una familia de funciones continuas estrictamente crecientes monotónicamente y $$\lim_{a\to0^+}f_a(x) = \lim_{a\to\infty}f_a^{-1}(y) = 0$$

Entonces las curvas de arclongitud constante comenzarán y terminarán en los puntos $(0,L)$ y $(L,0)$ . Tomemos como ejemplo la familia de curvas de aspecto similar

$$y = \frac{\cosh(ax)-1}{a}\implies L = \frac{\sinh(ax)}{a}$$

Las curvas de arclitud constante son de la forma

$$\vec{r}(a) = \left(\frac{\sinh^{-1}(aL)}{a},\frac{\sqrt{1+a^2L^2}-1}{a}\right)$$

A continuación se muestra un gráfico (lateral) de la curva de arclength $L=1$ (junto con la familia de curvas evaluadas en $a=\frac{1}{2},1,2,4,$ y $10$ ), que tiene una ecuación explícita de la forma

$$x = \frac{\tanh^{-1}y}{y}\cdot(1-y^2)$$ enter image description here

Estas curvas y la familia original de parábolas en cuestión tienen ambas esta propiedad, así como los círculos perfectos obtenidos de la familia $f_a(x) = ax$ . La razón por la que la pregunta original era difícil de resolver era la fórmula arclength no analíticamente invertible

23voto

Claude Leibovici Puntos 54392

$L$ siendo la longitud de arco conocida, sea $x_1=\frac t{2a}$ y $k=4a_1L$ entonces hay que resolver $t$ la ecuación $$k=t\sqrt{t^2+1} +\sinh ^{-1}(t)$$ Una buena aproximación viene dada por $t_0=\sqrt k$ .

Ahora, utilizando una serie de Taylor alrededor de $t=t_0$ y luego la reversión en serie da $$t_1=\sqrt{k}+z-\frac{\sqrt{k} }{2 (k+1)}z^2+\frac{(3 k-1) }{6 (k+1)^2}z^3+\frac{(13-15 k) \sqrt{k} }{24 (k+1)^3}z^4+\cdots$$ donde $$z=-\frac{\sqrt{k(k+1)} +\sinh ^{-1}\left(\sqrt{k}\right)-k}{2 \sqrt{k+1}}$$

Intentemos $k=10^n$ $$\left( \begin{array}{ccc} n & \text{estimate} & \text{solution} \\ 0 & 0.4810185 & 0.4819447 \\ 1 & 2.7868504 & 2.7868171 \\ 2 & 9.8244940 & 9.8244940 \\ 3 & 31.549250 & 31.549250 \\ 4 & 99.971006 & 99.971006 \\ 5 & 316.21678 & 316.21678 \\ 6 & 999.99595 & 999.99595 \end{array} \right)$$ Parece bastante decente.

11voto

David Puntos 11

Esto es lo más eficiente que puedo ver:

Tome su antiderivada (sustituyendo el sin por un sinh) y defina $a_1 x \equiv y$ para que

$$f(y) = a_1 * (\textrm{arc length}),$$

$$f(y) \equiv \frac{1}{2}(y \sqrt{1+y^2}+ \sinh^{-1} y).$$

$f(y)$ es monótona y tiene algunas buenas aproximaciones cuando $y \ll 1, y \gg 1$ . En esos casos podría ser posible obtenerla analíticamente. En general, inviértela numéricamente. Entonces,

$$x = a_1^{-1} f^{-1}(a_1 * (\textrm{arc length}).$$

La utilidad de hacerlo de esta manera es que no tienes que seguir invirtiendo una nueva función para cada nuevo $a_1$ sólo tienes que hacerlo una vez para poder aplanar cualquier parábola que quieras.

9voto

Tim Pederick Puntos 141

TLDR:

Quiere resolver para $v$ en esta ecuación: $$v\sqrt{v^2+1} + \sinh^{-1} v = 2a_1x_0\sqrt{u^2+1} + \frac{\sinh^{-1} u}{2a_0}$$

Y entonces $x_1=\frac{v}{2a_1}$ es su solución.

Llevo horas dándole vueltas a esto, no porque crea que pueda ayudarte -mucho de esto es nuevo para mí- sino porque me ha parecido interesante. Todavía no he encontrado una solución (aproximada) para $v$ Sin embargo. Y parece que Respuesta de Claude se me ha adelantado, de todos modos.

Pero aquí está mi trabajo, sólo porque no puedo soportar golpear Descartar en todo esto.


Creo que nunca había hecho la longitud de una curva parabólica. Así que mientras esto probablemente iría más suavemente si se hace de la definición integral de longitud de arco Voy a tomar el camino más fácil: Wikipedia tiene una fórmula de longitud de arco parabólico . Intentémoslo.


Entonces tenemos una parábola $f\left(x\right)=a_0 x^2$ . Wikipedia nos dice que podemos encontrar la longitud del arco, desde el vértice en $\left(0,0\right)$ a cualquier punto $\left(x,f\left(x\right)\right)$ en la parábola, utilizando estos valores:

  • La distancia focal $l$ de la parábola; en este caso, $l=\frac{1}{4a_0}$
  • La distancia perpendicular $p$ entre el punto y el eje de simetría; en este caso es simplemente $p=x$

Entonces, dado $h=\frac{p}{2}$ y $q=\sqrt{l^2+h^2}$ la longitud del arco es:

$$s=\frac{hq}{l}+l\ln\frac{h+q}{l}$$

Simplifiquemos. Dado que $h=\frac{x}{2}$

$$\begin{align} q &= \sqrt{\frac{1}{16a_0^2}+\frac{x^2}{4}} \\ &= \sqrt{\frac{4a_0^2x^2+1}{16a_0^2}} \\ &= l\sqrt{4a_0^2x^2+1} \end{align}$$

Así:

$$s=\frac{x}{2}\sqrt{4a_0^2x^2+1}+\frac{1}{4a_0}\ln\left(2a_0x+\sqrt{4a_0^2x^2+1}\right)$$


Ahora, tenemos otra parábola $g(x)=a_1 x^2$ tal que las longitudes de arco de $f\left(x_0\right)$ y $g\left(x_1\right)$ son iguales, es decir

$$\begin{align} \frac{x_0}{2}\sqrt{4a_0^2x_0^2+1}+\frac{1}{4a_0}\ln\left(2a_0x_0+\sqrt{4a_0^2x_0^2+1}\right) &= \frac{x_1}{2}\sqrt{4a_1^2x_1^2+1}+\frac{1}{4a_1}\ln\left(2a_1x_1+\sqrt{4a_1^2x_1^2+1}\right) \\ \therefore x_0\sqrt{4a_0^2x_0^2+1}+\frac{1}{2a_0}\ln\left(2a_0x_0+\sqrt{4a_0^2x_0^2+1}\right) &= x_1\sqrt{4a_1^2x_1^2+1}+\frac{1}{2a_1}\ln\left(2a_1x_1+\sqrt{4a_1^2x_1^2+1}\right) \end{align}$$

Y queremos resolver para $x_1$ en términos de $a_0$ , $a_1$ y $x_0$ . Nada más sencillo. </sarc>


Definamos $u=2a_0x_0$ (así $u^2=4a_0^2x_0^2$ ) y $v=2a_1x_1$ (así $v^2=4a_1^2x_1^2$ ). Eso acorta las cosas a:

$$x_0\sqrt{u^2+1} + \frac{1}{2a_0}\ln\left(u+\sqrt{u^2+1}\right) = x_1\sqrt{v^2+1} + \frac{1}{2a_1}\ln\left(v+\sqrt{v^2+1}\right)$$

Ahora veo dónde $\sinh$ ¡entra en las otras respuestas! Es porque $\sinh^{-1} x = \ln\left(x+\sqrt{x^2+1}\right)$ así que tenemos:

$$\begin{align} x_0\sqrt{u^2+1} + \frac{\sinh^{-1} u}{2a_0} &= x_1\sqrt{v^2+1} + \frac{\sinh^{-1} v}{2a_1} \\ \therefore 2a_1x_0\sqrt{u^2+1} + \frac{\sinh^{-1} u}{2a_0} &= 2a_1x_1\sqrt{v^2+1} + \sinh^{-1} v \\ &= v\sqrt{v^2+1} + \sinh^{-1} v \end{align}$$

Ese lado derecho hace no parecen fáciles de invertir. Lo admito, me evadí y le pedí a Wolfram Alpha que lo hiciera. Y, por supuesto, me dice: "no se ha encontrado ningún resultado en términos de funciones matemáticas estándar". [ suspiro ]


Basado en Respuesta de Tyma Gaidash Fui a buscar en el Teorema de inversión de Lagrange . Mi educación orientada a la ingeniería nunca cubrió esto, pero yo piense en He comprendido lo básico. Según tengo entendido, para resolver $y=f(x)$ para $x$ elegimos unos $z$ tal que $f(z)$ se define y $f'(z)\ne 0$ .

Acortemos todo el lado izquierdo de la ecuación a $w$ y defina $w=g\left(v\right)=v\sqrt{v^2+1} + \sinh^{-1} v$ . Primero, encontremos la derivada haciendo trampas y usando Wolfram Alpha: $g'(v)=2\sqrt{v^2+1}$ .

Necesitamos un valor $z$ donde $g'(z)\ne 0$ . Convenientemente, esta derivada no es cero en ningún lugar de los reales, así que es trivial. Creo (pero no estoy seguro) que $z$ debe aproximarse a $v$ así que supongamos aleatoriamente que $x_1\approx 1$ y así $v\approx 2a_1 = z$ .

Ahora, por el teorema de inversión, la función inversa $v=g^{-1}\left(w\right)$ es:

$$\begin{align} v &= z+\sum_{n=1}^\infty \left[\frac{\left(w-g\left(z\right)\right)^n}{n!} \lim_{t\to z} \frac{d^{n-1}}{dt^{n-1}}\left(\frac{t-z}{g(t)-g(z)}\right)^n\right] \\ &= z+\left(w-g\left(z\right)\right)\lim_{t\to z}\frac{t-z}{g(t)-g(z)} + \frac{\left(w - g\left(z\right)\right)^2}{2} \lim_{t\to z}\frac{d}{dt} \left(\frac{t-z}{g(t)-g(z)}\right)^2 + \cdots \end{align}$$

¿No parece que el primer límite igual que ¿el recíproco de la derivada? Así que el segundo término de la serie se convierte en $\frac{w-g\left(z\right)}{2\sqrt{z^2+1}}$ .

Y he estado trabajando en esto durante demasiado tiempo, así que aquí es donde voy a parar por esta noche.

3voto

Narasimham Puntos 7596

Para una parábola con parametrización

$$ x= at,y= a t^2/2 ;\;x^2= 2 a y \;; \text{slope} \;t=\tan \phi; \tag1$$

Diferencie $x^2$ , primed wrt arc length

$$ 2 x \cos \phi = 2 a \sin \phi ;\; x= a \tan \phi ;\;x'= \cos \phi= a \sec^2 \phi \;\phi' \tag 2 $$

de donde proviene la curvatura

$$ a \phi'=a \kappa=\cos^3 \phi \tag 3$$

Una salida más fácil/directa es la integración numérica directa de la oda en(3). Se puede establecer una fracción k de la longitud máxima del arco como parámetro para las integradas requeridas ( k= 2/3 en este caso particular), haciendo las parábolas del subconjunto más planas o más profundas ajustando $a$ .

enter image description here

La longitud total de un lado es una constante dada $L$

$$ L = \int _0^{\phi_m} \frac{ d \phi}{\kappa} ; \text{ now plug in curvature from(3) and integrate }$$

$$ \frac{L}{a}=\int _0^{\phi_m} \sec^3 \phi\; d\phi = \frac12\bigg[\log\bigg(\tan\bigg(\frac{\pi}{4} + \frac{\phi_m}{2} \bigg)\bigg) +\sec \phi_m \tan \phi_m \bigg]\tag 4 $$

Enchufe de (2) $x_{m}=a \tan \phi_m $ $$ \frac{2L}{a}= \log\left(\tan\bigg(\frac{\pi}{4} + \frac{\tan^{-1}(x_m/a)}{2}\right)\bigg)+ (x_m/a) \sqrt{(x_m/a)^2-1} \tag 5 $$

que es una función implícita $ f(a,x_m,L) $ , trazado suponiendo que un brazo de parábola tiene una longitud de arco dada $=1.8,$ en Mathematica lo que permite trazar parábolas más planas o profundas.

enter image description here

Ha quedado claro que hay dos criterios para arcos parabólicos de igual longitud que conectan $x_{max} $ a $ \text { a = 2* focal-length }$ .

EDITAR 1/2:

Para salir del aparente dilema he calculado/trazado cuidadosamente casos especiales para $( x_{max},a)$ combinaciones:

$$(0.4,0.0458),(0.8,0.1948), (1.2,0.4704),(1.6,0.8874)(2.0, 1.42264),(2.4,2.0101),(2.8,2.5825)$$

Todos los arcos tienen la misma longitud pero para $( x_{max} =2.0,2.4,2.8) $ la elección de $a$ parece que se ha pasado a la segunda crítica. Relación $(x_{max},a) $ no es único por la primera trama... ahora se examina más a fondo.

enter image description here

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