Articulo de referencia

Integral de Fresnel

Gráficas de S ( x ) y C ( x ) . El máximo de C ( x ) es aproximadamente 0,977 451 424 . Si los integrandos de S y C se definieran usando 2 "}},"i":0}}]}"> ⁠ π / 2 ⁠ t 2 en lugar...

Gráficas de S ( x ) y C ( x ) . El máximo de C ( x ) es aproximadamente0,977 451 424 . Si los integrandos de S y C se definieran usando π / 2 t 2 en lugar de t 2 , entonces la imagen se escalaría vertical y horizontalmente (ver más abajo).

Las integrales de Fresnel S ( x ) y C ( x ) , y sus funciones auxiliares F ( x ) y G ( x ) son funciones trascendentales que reciben su nombre de Augustin-Jean Fresnel y que se utilizan en óptica . Están estrechamente relacionadas con la función de error ( erf ). Surgen en la descripción de los fenómenos de difracción de Fresnel de campo cercano y se definen mediante las siguientes representaciones integrales :

S(incógnita)=0incógnitapecado(t2)dt,do(incógnita)=0incógnitaporque(t2)dt,F(incógnita)=(12π2S(incógnita))porque(incógnita2)(12π2do(incógnita))pecado(incógnita2),GRAMO(incógnita)=(12π2S(incógnita))pecado(incógnita2)+(12π2do(incógnita))porque(incógnita2).{\displaystyle {\begin{aligned}S(x)&=\int _{0}^{x}\sin \left(t^{2}\right)\,dt,\\C(x)&=\int _{0}^{x}\cos \left(t^{2}\right)\,dt,\\F(x)&=\left({\frac {1}{2}}{\sqrt {\frac {\pi }{2}}}-S\left(x\right)\right)\cos \left(x^{2}\right)-\left({\frac {1}{2}}{\sqrt {\frac {\pi }{2}}}-C\left(x\right)\right)\sin \left(x^{2}\right),\\G(x)&=\left({\frac {1}{2}}{\sqrt {\frac {\pi }{2}}}-S\left(x\right)\right)\sin \left(x^{2}\right)+\left({\frac {1}{2}}{\sqrt {\frac {\pi }{2}}}-C\left(x\right)\right)\cos \left(x^{2}\right).\end{aligned}}}

La curva paramétrica(S(t),do(t)){\displaystyle {\bigl (}S(t),C(t){\bigr )}} es la espiral de Euler o clotoide, una curva cuya curvatura varía linealmente con la longitud del arco.

El término integral de Fresnel también puede referirse a la integral definida compleja.

mi±iaincógnita2dincógnita=πami±iπ/4{\displaystyle \int _{-\infty }^{\infty }e^{\pm iax^{2}}dx={\sqrt {\frac {\pi }{a}}}e^{\pm i\pi /4}}

donde a es real y positivo; esto se puede evaluar cerrando un contorno en el plano complejo y aplicando el teorema integral de Cauchy .

Definición

Las integrales de Fresnel con argumentos π / 2 t 2 en lugar de t 2 convergen a 1 / 2 en lugar de 1 / 2 · π2 .

Las integrales de Fresnel admiten las siguientes series de Maclaurin que convergen para todo x : S(incógnita)=0incógnitapecado(t2)dt=norte=0(1)norteincógnita4norte+3(2norte+1)¡(4norte+3),do(incógnita)=0incógnitaporque(t2)dt=norte=0(1)norteincógnita4norte+1(2norte)¡(4norte+1).{\displaystyle {\begin{aligned}S(x)&=\int _{0}^{x}\sin \left(t^{2}\right)\,dt=\sum _{n=0}^{\infty }(-1)^{n}{\frac {x^{4n+3}}{(2n+1)!(4n+3)}},\\C(x)&=\int _{0}^{x}\cos \left(t^{2}\right)\,dt=\sum _{n=0}^{\infty }(-1)^{n}{\frac {x^{4n+1}}{(2n)!(4n+1)}}.\end{aligned}}}

Algunas tablas de uso común [ 1 ] [ 2 ] usan π / 2 t 2 en lugar de t 2 para el argumento de las integrales que definen S ( x ) y C ( x ) . Esto cambia sus límites en el infinito de 1 / 2 · π / 2 a 1 / 2[ 3 ] y la longitud de arco para la primera vuelta de la espiral de 2 π a 2 (en t = 2 ). Estas funciones alternativas se conocen generalmente como integrales de Fresnel normalizadas .

Las funciones auxiliares F ( x ) y G ( x ) proporcionan cotas monótonas para las integrales de Fresnel: [ 4 ]12π2F(incógnita)GRAMO(incógnita)do(incógnita)12π2+F(incógnita)+GRAMO(incógnita),12π2F(incógnita)GRAMO(incógnita)S(incógnita)12π2+F(incógnita)+GRAMO(incógnita).{\displaystyle {\begin{aligned}{\frac {1}{2}}{\sqrt {\frac {\pi }{2}}}-F(x)-G(x)\leq C(x)\leq {\frac {1}{2}}{\sqrt {\frac {\pi }{2}}}+F(x)+G(x),\\{\frac {1}{2}}{\sqrt {\frac {\pi }{2}}}-F(x)-G(x)\leq S(x)\leq {\frac {1}{2}}{\sqrt {\frac {\pi }{2}}}+F(x)+G(x).\end{aligned}}}

espiral de Euler

Espiral de Euler ( x , y ) = ( C ( t ), S ( t )) . La espiral converge al centro de los agujeros en la imagen cuando t tiende a infinito positivo o negativo.
Animación que representa la evolución de una espiral de Cornu con el círculo tangente con el mismo radio de curvatura que en su extremo, también conocido como círculo osculador .

La espiral de Euler, también conocida como espiral de Cornu o clotoide, es la curva generada por una gráfica paramétrica de S ( t ) frente a C ( t ) . La espiral de Euler fue estudiada por primera vez a mediados del siglo XVIII por Leonhard Euler en el contexto de la teoría de vigas de Euler-Bernoulli . Un siglo después, Marie Alfred Cornu construyó la misma espiral como nomograma para cálculos de difracción.

A partir de las definiciones de las integrales de Fresnel, los infinitesimales dx y dy son, por lo tanto: dincógnita=do(t)dt=porque(t2)dt,dy=S(t)dt=pecado(t2)dt.{\displaystyle {\begin{aligned}dx&=C'(t)\,dt=\cos \left(t^{2}\right)\,dt,\\dy&=S'(t)\,dt=\sin \left(t^{2}\right)\,dt.\end{aligned}}}

Por lo tanto, la longitud de la espiral medida desde el origen se puede expresar como L=0t0dincógnita2+dy2=0t0dt=t0.{\displaystyle L=\int _{0}^{t_{0}}{\sqrt {dx^{2}+dy^{2}}}=\int _{0}^{t_{0}}dt=t_{0}.}

Es decir, el parámetro t es la longitud de la curva medida desde el origen (0, 0) , y la espiral de Euler tiene longitud infinita . El vector (cos( ​​t 2 ), sin( t 2 )) , donde θ = t 2 , también expresa el vector tangente unitario a lo largo de la espiral. Dado que t es la longitud de la curva, la curvatura κ se puede expresar como κ=1R=dθdt=2t.{\displaystyle \kappa ={\frac {1}{R}}={\frac {d\theta }{dt}}=2t.}

Por lo tanto, la tasa de cambio de la curvatura con respecto a la longitud de la curva es dκdt=d2θdt2=2.{\displaystyle {\frac {d\kappa }{dt}}={\frac {d^{2}\theta }{dt^{2}}}=2.}

Una espiral de Euler tiene la propiedad de que su curvatura en cualquier punto es proporcional a la distancia a lo largo de la espiral, medida desde el origen. Esta propiedad la hace útil como curva de transición en ingeniería de carreteras y ferrocarriles: si un vehículo sigue la espiral a velocidad unitaria, el parámetro t en las derivadas anteriores también representa el tiempo. En consecuencia, un vehículo que sigue la espiral a velocidad constante tendrá una tasa constante de aceleración angular .

Es común incorporar secciones de espirales de Euler en la forma de los bucles de las montañas rusas para crear lo que se conoce como bucles clotoides .

Propiedades

C ( x ) y S ( x ) son funciones impares de x ,

do(incógnita)=do(incógnita),S(incógnita)=S(incógnita).{\displaystyle C(-x)=-C(x),\quad S(-x)=-S(x).}

lo cual se puede apreciar fácilmente por el hecho de que sus expansiones en series de potencias solo tienen términos de grado impar, o bien porque son antiderivadas de funciones pares que también son cero en el origen.

Las asintóticas de las integrales de Fresnel cuando x → ∞ vienen dadas por las fórmulas:

S(incógnita)=18πsgnincógnita[1+O(incógnita4)](porque(incógnita2)2incógnita+pecado(incógnita2)4incógnita3),do(incógnita)=18πsgnincógnita+[1+O(incógnita4)](pecado(incógnita2)2incógnitaporque(incógnita2)4incógnita3).{\displaystyle {\begin{aligned}S(x)&={\sqrt {{\tfrac {1}{8}}\pi }}\operatorname {sgn} x-\left[1+O\left(x^{-4}\right)\right]\left({\frac {\cos \left(x^{2}\right)}{2x}}+{\frac {\sin \left(x^{2}\right)}{4x^{3}}}\right),\\[6px]C(x)&={\sqrt {{\tfrac {1}{8}}\pi }}\operatorname {sgn} x+\left[1+O\left(x^{-4}\right)\right]\left({\frac {\sin \left(x^{2}\right)}{2x}}-{\frac {\cos \left(x^{2}\right)}{4x^{3}}}\right).\end{aligned}}}

Integral de Fresnel compleja S ( z )

Utilizando las expansiones en series de potencias anteriores, las integrales de Fresnel se pueden extender al dominio de los números complejos , donde se convierten en funciones enteras de la variable compleja z .

Las integrales de Fresnel se pueden expresar utilizando la función de error de la siguiente manera: [ 5 ]

Integral de Fresnel compleja C ( z )

S(z)=π21+i4[terreno(1+i2z)iterreno(1i2z)],do(z)=π21i4[terreno(1+i2z)+iterreno(1i2z)].{\displaystyle {\begin{aligned}S(z)&={\sqrt {\frac {\pi }{2}}}\cdot {\frac {1+i}{4}}\left[\operatorname {erf} \left({\frac {1+i}{\sqrt {2}}}z\right)-i\operatorname {erf} \left({\frac {1-i}{\sqrt {2}}}z\right)\right],\\[6px]C(z)&={\sqrt {\frac {\pi }{2}}}\cdot {\frac {1-i}{4}}\left[\operatorname {erf} \left({\frac {1+i}{\sqrt {2}}}z\right)+i\operatorname {erf} \left({\frac {1-i}{\sqrt {2}}}z\right)\right].\end{aligned}}}

o

do(z)+iS(z)=π21+i2terreno(1i2z),S(z)+ido(z)=π21+i2terreno(1+i2z).{\displaystyle {\begin{aligned}C(z)+iS(z)&={\sqrt {\frac {\pi }{2}}}\cdot {\frac {1+i}{2}}\operatorname {erf} \left({\frac {1-i}{\sqrt {2}}}z\right),\\[6px]S(z)+iC(z)&={\sqrt {\frac {\pi }{2}}}\cdot {\frac {1+i}{2}}\operatorname {erf} \left({\frac {1+i}{\sqrt {2}}}z\right).\end{aligned}}}

Límites cuando x tiende a infinito

Las integrales que definen C ( x ) y S ( x ) no pueden evaluarse en forma cerrada en términos de funciones elementales , excepto en casos especiales. Se conocen los límites de estas funciones cuando x tiende a infinito:0porque(t2)dt=0pecado(t2)dt=2π4=π80,6267.{\displaystyle \int _{0}^{\infty }\cos \left(t^{2}\right)\,dt=\int _{0}^{\infty }\sin \left(t^{2}\right)\,dt={\frac {\sqrt {2\pi }}{4}}={\sqrt {\frac {\pi }{8}}}\approx 0.6267.}

Generalización

El integral incógnitametromiiincógnitanortedincógnita=k=0ikincógnitametro+nortekk¡dincógnita=k=0ik(metro+nortek+1)incógnitametro+nortek+1k¡{\displaystyle \int x^{m}e^{ix^{n}}\,dx=\int \sum _{k=0}^{\infty }{\frac {i^{k}x^{m+nk}}{k!}}\,dx=\sum _{k=0}^{\infty }{\frac {i^{k}}{(m+nk+1)}}{\frac {x^{m+nk+1}}{k!}}} is a confluent hypergeometric function and also an incomplete gamma function[7]xmeixndx=xm+1m+11F1(m+1n1+m+1nixn)=1nim+1nγ(m+1n,ixn),{\displaystyle {\begin{aligned}\int x^{m}e^{ix^{n}}\,dx&={\frac {x^{m+1}}{m+1}}\,_{1}F_{1}\left({\begin{array}{c}{\frac {m+1}{n}}\\1+{\frac {m+1}{n}}\end{array}}\mid ix^{n}\right)\\[6px]&={\frac {1}{n}}i^{\frac {m+1}{n}}\gamma \left({\frac {m+1}{n}},-ix^{n}\right),\end{aligned}}} which reduces to Fresnel integrals if real or imaginary parts are taken: xmsin(xn)dx=xm+n+1m+n+11F2(12+m+12n32+m+12n,32x2n4).{\displaystyle \int x^{m}\sin(x^{n})\,dx={\frac {x^{m+n+1}}{m+n+1}}\,_{1}F_{2}\left({\begin{array}{c}{\frac {1}{2}}+{\frac {m+1}{2n}}\\{\frac {3}{2}}+{\frac {m+1}{2n}},{\frac {3}{2}}\end{array}}\mid -{\frac {x^{2n}}{4}}\right).} The leading term in the asymptotic expansion is 1F1(m+1n1+m+1nixn)m+1nΓ(m+1n)eiπm+12nxm1,{\displaystyle _{1}F_{1}\left({\begin{array}{c}{\frac {m+1}{n}}\\1+{\frac {m+1}{n}}\end{array}}\mid ix^{n}\right)\sim {\frac {m+1}{n}}\,\Gamma \left({\frac {m+1}{n}}\right)e^{i\pi {\frac {m+1}{2n}}}x^{-m-1},} and therefore 0xmeixndx=1nΓ(m+1n)eiπm+12n.{\displaystyle \int _{0}^{\infty }x^{m}e^{ix^{n}}\,dx={\frac {1}{n}}\,\Gamma \left({\frac {m+1}{n}}\right)e^{i\pi {\frac {m+1}{2n}}}.}

For m = 0, the imaginary part of this equation in particular is 0sin(xa)dx=Γ(1+1a)sin(π2a),{\displaystyle \int _{0}^{\infty }\sin \left(x^{a}\right)\,dx=\Gamma \left(1+{\frac {1}{a}}\right)\sin \left({\frac {\pi }{2a}}\right),} with the left-hand side converging for |a| > 1 and the right-hand side being its analytical extension to the whole plane less where lie the poles of Γ(a−1).

The Kummer transformation of the confluent hypergeometric function is xmeixndx=Vn,m(x)eixn,{\displaystyle \int x^{m}e^{ix^{n}}\,dx=V_{n,m}(x)e^{ix^{n}},} with Vn,m:=xm+1m+11F1(11+m+1nixn).{\displaystyle V_{n,m}:={\frac {x^{m+1}}{m+1}}\,_{1}F_{1}\left({\begin{array}{c}1\\1+{\frac {m+1}{n}}\end{array}}\mid -ix^{n}\right).}

Numerical approximation

For computation to arbitrary precision, the power series is suitable for small argument. For large argument, asymptotic expansions converge faster.[8] Continued fraction methods may also be used.[9]

For computation to particular target precision, other approximations have been developed. Cody[10] developed a set of efficient approximations based on rational functions that give relative errors down to 2×10−19. A FORTRAN implementation of the Cody approximation that includes the values of the coefficients needed for implementation in other languages was published by van Snyder.[11] Boersma developed an approximation with error less than 1.6×10−9.[12]

Applications

The Fresnel integrals were originally used in the calculation of the electromagnetic field intensity in an environment where light bends around opaque objects.[13] More recently, they have been used in the design of highways and railways, specifically their curvature transition zones, see track transition curve.[14] Other applications are rollercoasters[13] or calculating the transitions on a velodrome track to allow rapid entry to the bends and gradual exit.

Véase también

Notas

  1. ^ Abramowitz y Stegun 1983 , ecuación 7.3.1–7.3.2.
  2. Temme 2010 .
  3. ^ Abramowitz y Stegun 1983 , ecuación 7.3.20.
  4. Oldham, Keith B.; Myland, Jan C.; Spanier, Jerome; Myland, Jan (2009). Un atlas de funciones: con equator, la calculadora de funciones del atlas . Nueva York, NY: Springer US Springer e-books. ISBN 978-0-387-48807-3.
  5. functions.wolfram.com, Integral de Fresnel S: Representaciones mediante funciones equivalentes e Integral de Fresnel C: Representaciones mediante funciones equivalentes . Nota: Wolfram utiliza la convención de Abramowitz y Stegun, que difiere de la de este artículo por factores deπ2 .
  6. Otro método basado en la integración paramétrica se describe, por ejemplo, en Zajta y Goel 1989 .
  7. Mathar 2012 .
  8. Temme 2010 , §7.12(ii).
  9. Prensa et al. 2007 .
  10. Cody 1968 .
  11. van Snyder 1993 .
  12. Boersma 1960 .
  13. 1 2 Beatty 2013 .
  14. Stewart 2008 , pág. 383.

Referencias

  • Abramowitz, Milton ; Stegun, Irene Ann , eds. (1983) [junio de 1964]. «Capítulo 7». Manual de funciones matemáticas con fórmulas, gráficas y tablas matemáticas . Serie de Matemáticas Aplicadas. Vol.  55 (novena reimpresión con correcciones adicionales de la décima edición original con correcciones (diciembre de 1972); primera  ed.). Washington D. C.; Nueva York: Departamento de Comercio de los Estados Unidos, Oficina Nacional de Normas; Dover Publications. ISBN 978-0-486-61272-0. LCCN 64-60036 . MR 0167642 . LCCN 65-12253 .   
  • Alazah, Mohammad (2012). "Cálculo de integrales de Fresnel mediante reglas del trapecio modificadas". Matemática numérica . 128 (4): 635–661 . arXiv : 1209.3451 . Código Bib : 2012arXiv1209.3451A . doi : 10.1007/s00211-014-0627-z . S2CID 13934493 . 
  • Beatty, Thomas (2013). "Cómo evaluar integrales de Fresnel" (PDF) . FGCU Math - Verano 2013. Archivado del original (PDF) el 9 de abril de 2023. Recuperado el 27 de julio de 2013 .
  • Boersma, J. (1960). "Cálculo de integrales de Fresnel" . Math. Comp . 14 (72): 380. doi : 10.1090/S0025-5718-1960-0121973-3 . MR 0121973 . 
  • Bulirsch, Roland (1967). "Cálculo numérico de las integrales de seno, coseno y Fresnel". Numer. Math . 9 (5): 380– 385. doi : 10.1007/BF02162153 . S2CID 121794086 . 
  • Cody, William J. (1968). "Aproximaciones de Chebyshev para las integrales de Fresnel" (PDF) . Math. Comp . 22 (102): 450– 453. doi : 10.1090/S0025-5718-68-99871-2 .
  • Hangelbroek, RJ (1967). "Aproximación numérica de integrales de Fresnel mediante polinomios de Chebyshev". J. Eng. Math . 1 (1): 37– 50. Bibcode : 1967JEnMa...1...37H . doi : 10.1007/BF01793638 . S2CID 122271446 . 
  • Mathar, RJ (2012). "Expansión en serie de integrales de Fresnel generalizadas". arXiv : 1211.3963 [ math.CA ].
  • Nave, R. (2002). «La espiral de Cornu» .(Utiliza π / 2 t 2 en lugar de t 2 .)
  • Press, WH; Teukolsky, SA; Vetterling, WT; Flannery, BP (2007). «Sección 6.8.1. Integrales de Fresnel» . Numerical Recipes: The Art of Scientific Computing (3.ª  ed.). Nueva York: Cambridge University Press. ISBN 978-0-521-88068-8Archivado del original el 11 de agosto de 2011. Consultado el 9 de agosto de 2011 .
  • van Snyder, W. (1993). "Algoritmo 723: Integrales de Fresnel" . ACM Trans. Math. Softw . 19 (4): 452– 456. doi : 10.1145/168173.168193 . S2CID 12346795 . 
  • Stewart, James (2008). Cálculo: Primeros trascendentales . Cengage Learning EMEA. ISBN 978-0-495-38273-7.
  • Temme, NM (2010), "Funciones de error, integrales de Dawson y Fresnel" , en Olver, Frank WJ ; Lozier, Daniel M.; Boisvert, Ronald F.; Clark, Charles W. (eds.), NIST Handbook of Mathematical Functions , Cambridge University Press, ISBN 978-0-521-19225-5, MR 2723248 .
  • van Wijngaarden, A.; Scheen, WL (1949). Tabla de Integrales de Fresnel . Verhandl. Konink. Ned. Akád. Wetenschapen. vol.  19.
  • Zajta, Aurel J.; Goel, Sudhir K. (1989). "Técnicas de integración paramétrica". Mathematics Magazine . 62 (5): 318– 322. doi : 10.1080/0025570X.1989.11977462 .