Transformada_de_Hankel Knowpia

En matemáticas, la transformada de Hankel es una transformada integral, desarrollada por primera vez por el matemático Hermann Hankel, que expresa una función $f(r)$ como suma ponderada de un número infinito de funciones de Bessel de primer tipo $J_{\nu }(kr)$ . También se conoce como transformada de Fourier-Bessel. Las funciones de Bessel del núcleo de la integral son todas del mismo orden. $\nu$ , pero difieren en el factor de escala $k$ a lo largo del eje $r$ . El coeficiente $F_{\nu }$ de cada función de Bessel, vista como una función del factor de escala $k$ , constituye la transformada de Hankel. La transformada de Hankel está estrechamente relacionada con la serie de Fourier-Bessel, de la misma manera que la transformada de Fourier para un intervalo infinito está relacionada con la serie de Fourier en un intervalo finito.

Definición

La transformada de Hankel de orden $\nu$ de una función $f(r)$ es:

F_{\nu }(k)=\int _{0}^{\infty }f(r)J_{\nu }(kr)\,r\,\mathrm {d} r,

dónde $J_{\nu }$ es la función de Bessel del primer tipo de orden $\nu$ , con $\nu \geq -1/2$ . La transformada inversa de Hankel de $F_{\nu }(k)$ Se define como

f(r)=\int _{0}^{\infty }F_{\nu }(k)J_{\nu }(kr)\,k\,\mathrm {d} k,

por lo que puede verificarse una relación de ortogonalidad entre las funciones de Bessel.

Dominio de definición

La inversión de la transformada de Hankel de una función. $f(r)$ es válido en todos los puntos donde $f(r)$ es continuo, siempre que esté definido y sea continuo por partes $(0,\infty )$ , con variación limitada en cada subintervalo finito de $(0,\infty )$ y

\int _{0}^{\infty }|f(r)|\,r^{\frac {1}{2}}\,\mathrm {d} r<\infty .

Sin embargo, en analogía con la transformada de Fourier, se puede ampliar el dominio mediante el razonamiento de densidad, incluyendo algunas funciones para las cuales la integral anterior no es finita, como $f(r)=(1+r)^{-3/2}$ .

Definición alternativa

Una definición alternativa establece que la transformada de Hankel de $g(r)$ es ^[1]

h_{\nu }(k)=\int _{0}^{\infty }g(r)J_{\nu }(kr)\,{\sqrt {kr}}\,\mathrm {d} r.

Las dos definiciones están relacionadas:

Si

g(r)=f(r){\sqrt {r}}

, entonces

h_{\nu }(k)=F_{\nu }(k){\sqrt {k}}.

Esto significa que, al igual que en la definición anterior, la transformada de Hankel definida de esta manera es su propia inversa:

g(r)=\int _{0}^{\infty }h_{\nu }(k)J_{\nu }(kr)\,{\sqrt {kr}}\,\mathrm {d} k.

El dominio ahora tiene la condición

\int _{0}^{\infty }|g(r)|\,\mathrm {d} r<\infty ,

pero se puede ampliar. Según de Branges, se puede tomar la integral como el límite con el límite superior tendiendo al infinito (una integral impropia en lugar de una integral de Lebesgue ), y de esta manera la transformada de Hankel y su inversa se definen para cada función en L ² ( 0, ∞).

Ortogonalidad

Las funciones de Bessel forman una base ortogonal cuando se ponderan con la función $r$ : ^[2]

\int _{0}^{\infty }J_{\nu }(kr)J_{\nu }(k'r)\,r\,\mathrm {d} r={\frac {\delta (k-k')}{k}},\quad k,k'>0.

Teorema de Plancherel y de Parseval

Si las funciones $f(r)$ y $g(r)$ poseen transformaciones de Hankel $F_{\nu }(k)$ y $G_{\nu }(k)$ bien definidas, entonces el teorema de Plancherel establece que

\int _{0}^{\infty }f(r)g(r)\,r\,\mathrm {d} r=\int _{0}^{\infty }F_{\nu }(k)G_{\nu }(k)\,k\,\mathrm {d} k.

El teorema de Parseval, que establece

\int _{0}^{\infty }|f(r)|^{2}\,r\,\mathrm {d} r=\int _{0}^{\infty }|F_{\nu }(k)|^{2}\,k\,\mathrm {d} k,

es un caso especial del teorema de Plancherel. Estos teoremas se pueden demostrar utilizando la propiedad de ortogonalidad.

Relaciones transformadas con los demás.

Relación con la transformada de Fourier (simetría circular)

La transformada de Hankel de orden cero es esencialmente la transformada de Fourier bidimensional de una función circularmente simétrica.

Si se considera una función bidimensional $f(\mathbf {r} )$ del radio vectorial $r$ . Su transformada de Fourier es

F(\mathbf {k} )={\frac {1}{2\pi }}\iint f(\mathbf {r} )e^{-i\mathbf {k} \cdot \mathbf {r} }\,\mathrm {d} \mathbf {r} .

Sin pérdida de generalidad, se puede elegir un sistema de coordenadas polares. $(r,\theta )$ para que el vector $\mathbf {k}$ acostarse en el eje $\theta =0$ (en el espacio K). La transformada de Fourier ahora se escribe en estas coordenadas como

F(\mathbf {k} )={\frac {1}{2\pi }}\int _{r=0}^{\infty }\int _{\theta =0}^{2\pi }f(r,\theta )e^{-ikr\cos(\theta )}\,r\,\mathrm {d} \theta \,\mathrm {d} r,

dónde $\theta$ es el ángulo entre los vectores $\mathbf {k}$ y $\mathbf {r}$ . Si la función $f$ es circularmente simétrico y no depende de la variable angular $\theta$ y se puede escribir como $f(r)$ . Por lo tanto, puede quedar fuera de la integración en $\theta$ , y en este caso la transformada de Fourier se convierte en

F(\mathbf {k} )=F(k)=\int _{0}^{\infty }f(r)J_{0}(kr)\,r\,\mathrm {d} r,

que es exactamente la transformada de Hankel de orden cero de $f(r)$ . De manera similar para la transformada inversa,

f(\mathbf {r} )={\frac {1}{2\pi }}\iint F(\mathbf {k} )e^{i\mathbf {k} \cdot \mathbf {r} }\,\mathrm {d} \mathbf {k} =\int _{0}^{\infty }F(k)J_{0}(kr)\,k\,\mathrm {d} k,

por lo tanto $f(r)$ es la transformada de Hankel de orden cero de $F(k)$ .

Relación con la transformada de Fourier (simetría radial en n dimensiones)

Para una transformada de Fourier de n dimensiones,

F(\mathbf {k} )={\frac {1}{(2\pi )^{n/2}}}\int f(\mathbf {r} )e^{-i\mathbf {k} \cdot \mathbf {r} }\,d^{n}\mathbf {r} ,

Si la función $f$ es radialmente simétrico, entonces ^[3]

k^{\frac {n-2}{2}}F(k)=\int _{0}^{\infty }r^{\frac {n-2}{2}}f(r)J_{\frac {n-2}{2}}(kr)\,r\,dr.

Relación con la transformada de Fourier (caso general en dos dimensiones)

Para generalizar, si $f$ se puede expandir en una serie de multipolos ,

f(r,\theta )=\sum _{m=-\infty }^{\infty }f_{m}(r)e^{im\theta },

y si $\theta _{k}$ es el ángulo entre la dirección de $\mathbf {k}$ y el eje $\theta =0$ , En ese tiempo

{\begin{aligned}F(\mathbf {k} )&={\frac {1}{2\pi }}\int _{0}^{\infty }r\,\mathrm {d} r\,\int _{0}^{2\pi }\mathrm {d} \theta \,f(r,\theta )e^{-ikr\cos(\theta -\theta _{k})}\\&={\frac {1}{2\pi }}\sum _{m}\int _{0}^{\infty }r\,\mathrm {d} r\,\int _{0}^{2\pi }\mathrm {d} \theta \,f_{m}(r)e^{im\theta }e^{-ikr\cos(\theta -\theta _{k})}\\&={\frac {1}{2\pi }}\sum _{m}e^{im\theta _{k}}\int _{0}^{\infty }r\,\mathrm {d} r\,f_{m}(r)\int _{0}^{2\pi }\mathrm {d} \varphi \,e^{im\varphi }e^{-ikr\cos \varphi }&&(\varphi =\theta -\theta _{k})\\&=\sum _{m}e^{im\theta _{k}}\int _{0}^{\infty }r\,\mathrm {d} r\,f_{m}(r)i^{-m}J_{m}(kr)\\&=\sum _{m}i^{-m}e^{im\theta _{k}}\int _{0}^{\infty }f_{m}(r)J_{m}(kr)\,r\,\mathrm {d} r\\&=\sum _{m}i^{-m}e^{im\theta _{k}}F_{m}(k),\end{aligned}}

dónde $F_{m}(k)$ es la transformada de orden de Hankel $m$ de $f_{m}(r)$ .

Funciones dentro de un radio limitado

Además, si $f_{m}$ es suficientemente suave cerca del origen y es cero fuera de una bola de radio $R$ , luego se puede ampliar a la serie Chebyshev :

f_{m}(r)=r^{m}\sum _{t\geq 0}f_{mt}\left(1-\left({\tfrac {r}{R}}\right)^{2}\right)^{t},\quad 0\leq r\leq R.

Sustituyéndolo en la última ecuación de la sección anterior se obtiene

{\begin{aligned}F(\mathbf {k} )&=\sum _{m}i^{-m}e^{im\theta _{k}}\sum _{t}f_{mt}\int _{0}^{R}r^{m}\left(1-\left({\tfrac {r}{R}}\right)^{2}\right)^{t}J_{m}(kr)r\,\mathrm {d} r&&\\&=\sum _{m}i^{-m}e^{im\theta _{k}}R^{m+2}\sum _{t}f_{mt}\int _{0}^{1}x^{m+1}(1-x^{2})^{t}J_{m}(kxR)\,\mathrm {d} x&&(x={\tfrac {r}{R}})\\&=\sum _{m}i^{-m}e^{im\theta _{k}}R^{m+2}\sum _{t}f_{mt}{\frac {t!2^{t}}{(kR)^{1+t}}}J_{m+t+1}(kR),\end{aligned}}

donde la última igualdad se deriva del §6.567.1 de ^[4] . Se trata de un caso mucho más general que el abordado en el apartado anterior. El aspecto numérico importante es que los coeficientes $f_{mt}$ se puede obtener utilizando técnicas de la transformada discreta de Fourier .

Este es un vistazo de la rápida transformada de Hankel.

Relación con las transformadas de Fourier y Abel

En dos dimensiones, si lo defines $A$ como el operador de transformación de Abel, $F$ como operador de transformada de Fourier, e $H$ como la transformada de Hankel de orden cero, entonces el caso especial del teorema de corte de proyección para funciones circularmente simétricas establece que

FA=H.

En otras palabras, aplicar la transformada de Abel a una función en una dimensión y luego realizar la transformada de Fourier equivale a aplicar la transformada de Hankel a la función. Este concepto se puede extender a todas las dimensiones.

Transformada de algunas funciones particulares^[5]

$f(r)$	$F_{0}(k)$
$1$	${\frac {\delta (k)}{k}}$
${\frac {1}{r}}$	${\frac {1}{k}}$
$r$	$-{\frac {1}{k^{3}}}$
$r^{3}$	${\frac {9}{k^{5}}}$
$r^{m}$	${\frac {2^{m+1}\Gamma \left({\tfrac {m}{2}}+1\right)}{k^{m+2}\Gamma \left(-{\tfrac {m}{2}}\right)}},\quad -2<\Re (m)<-{\tfrac {1}{2}}$
${\frac {1}{\sqrt {r^{2}+z^{2}}}}$	${\frac {e^{-k\|z\|}}{k}}$
${\frac {1}{z^{2}+r^{2}}}$	$K_{0}(kz),\quad z\in \mathbf {C}$
${\frac {e^{iar}}{r}}$	${\frac {i}{\sqrt {a^{2}-k^{2}}}},\quad a>0,k<a$
${\frac {e^{iar}}{r}}$	${\frac {1}{\sqrt {k^{2}-a^{2}}}},\quad a>0,k>a$
$e^{-{\frac {1}{2}}a^{2}r^{2}}$	${\frac {1}{a^{2}}}e^{-{\tfrac {k^{2}}{2a^{2}}}}$
${\frac {1}{r}}J_{0}(lr)e^{-sr}$	${\frac {2}{\pi {\sqrt {(k+l)^{2}+s^{2}}}}}K{\bigg (}{\sqrt {\frac {4kl}{(k+l)^{2}+s^{2}}}}{\bigg )}$
$-r^{2}f(r)$	${\frac {\mathrm {d} ^{2}F_{0}}{\mathrm {d} k^{2}}}+{\frac {1}{k}}{\frac {\mathrm {d} F_{0}}{\mathrm {d} k}}$

$f(r)$	$F_{\nu }(k)$
$r^{s}$	${\frac {2^{s+1}}{k^{s+2}}}{\frac {\Gamma \left({\tfrac {1}{2}}(2+\nu +s)\right)}{\Gamma ({\tfrac {1}{2}}(\nu -s))}}$
$r^{\nu -2s}\Gamma (s,r^{2}h)$	${\tfrac {1}{2}}\left({\tfrac {k}{2}}\right)^{2s-\nu -2}\gamma \left(1-s+\nu ,{\tfrac {k^{2}}{4h}}\right)$
$e^{-r^{2}}r^{\nu }U(a,b,r^{2})$	${\frac {\Gamma (2+\nu -b)}{2\Gamma (2+\nu -b+a)}}\left({\tfrac {k}{2}}\right)^{\nu }e^{-{\frac {k^{2}}{4}}}\,_{1}F_{1}\left(a,2+a-b+\nu ,{\tfrac {k^{2}}{4}}\right)$
$r^{n}J_{\mu }(lr)e^{-sr}$	expresable en términos de integrales elípticas . ^[6]
$-r^{2}f(r)$	${\frac {\mathrm {d} ^{2}F_{\nu }}{\mathrm {d} k^{2}}}+{\frac {1}{k}}{\frac {\mathrm {d} F_{\nu }}{\mathrm {d} k}}-{\frac {\nu ^{2}}{k^{2}}}F_{\nu }$

$K_{n}(z)$ es la función de Bessel modificada de segundo tipo . $K(z)$ es la integral elíptica completa de primer tipo .

La expresión

{\frac {\mathrm {d} ^{2}F_{0}}{\mathrm {d} k^{2}}}+{\frac {1}{k}}{\frac {\mathrm {d} F_{0}}{\mathrm {d} k}}

coincide con la expresión del operador de Laplace en coordenadas polares $(k,\theta )$ aplicado a la función esféricamente simétrica $F_{0}(k)$ .

La transformada de Hankel de los polinomios de Zernike son esencialmente funciones de Bessel (Noll 1976):

R_{n}^{m}(r)=(-1)^{\frac {n-m}{2}}\int _{0}^{\infty }J_{n+1}(k)J_{m}(kr)\,\mathrm {d} k

para $n-m\geq 0$ incluso.

Referencias

↑ Louis de Branges (1968). Hilbert spaces of entire functions. London: Prentice-Hall. p. 189. ISBN 978-0-13-388900-0.
↑ Ponce de Leon, J. (2015). «Revisiting the orthogonality of Bessel functions of the first kind on an infinite interval». European Journal of Physics 36 (1). doi:10.1088/0143-0807/36/1/015016.
↑ Faris, William G. «Radial functions and the Fourier transform: Notes for Math 583A, Fall 2008». University of Arizona, Department of Mathematics.
↑ Gradshteyn; Ryzhik, I. M. (2015). Table of Integrals, Series, and Products (Eighth edición). Academic Press. p. 687. ISBN 978-0-12-384933-5.
↑ Papoulis, Athanasios (1981). Systems and Transforms with Applications to Optics. Florida USA: Krieger Publishing Company. pp. 140-175. ISBN 978-0-89874-358-6.
↑ «Laplace transform of products of Bessel functions: A visitation of earlier formulas». Quarterly of Applied Mathematics 70: 77-97. 2012. doi:10.1090/s0033-569x-2011-01239-2.

Bibliografía

Gaskill, Jack D. (1978). Linear Systems, Fourier Transforms, and Optics. New York: John Wiley & Sons. ISBN 978-0-471-29288-3.
Polyanin; Manzhirov, A. V. (1998). Handbook of Integral Equations. Boca Raton: CRC Press. ISBN 978-0-8493-2876-3.
Smythe, William R. (1968). Static and Dynamic Electricity (3rd edición). New York: McGraw-Hill. pp. 179-223.
«On Hankel transforms». Proceedings of the London Mathematical Society 39 (2): 49-67. 1935. doi:10.1112/plms/s2-39.1.49.
«On certain integrals of Lipschitz-Hankel type involving products of Bessel Functions». Philosophical Transactions of the Royal Society A 247 (935): 529-551. 1955. Bibcode:1955RSPTA.247..529E. JSTOR 91565. doi:10.1098/rsta.1955.0005.
«Calculation of integrals of products of Bessel functions». Mathematics of Computation 21 (99): 407-412. 1967. doi:10.1090/S0025-5718-67-99149-1.
«The asymptotic expansions of Hankel transforms and related integrals». Mathematics of Computation 26 (118): 515-527. 1972. JSTOR 2003243. doi:10.1090/S0025-5718-1972-0308695-9.
«A note on the computation of integrals involving products of trigonometric and Bessel functions». Mathematics of Computation 27 (124): 871-872. 1973. JSTOR 2005522. doi:10.2307/2005522.
Noll, Robert J (1976). «Zernike polynomials and atmospheric turbulence». Journal of the Optical Society of America 66 (3): 207-211. Bibcode:1976JOSA...66..207N. doi:10.1364/JOSA.66.000207.
«Quasi-fast Hankel transform». Opt. Lett. 1 (1): 13-15. 1977. Bibcode:1977OptL....1...13S. doi:10.1364/OL.1.000013.
«High-accuracy fast Hankel transform for optical beam propagation». J. Opt. Soc. Am. A 9 (11): 2031-2033. 1992. Bibcode:1992JOSAA...9.2031M. doi:10.1364/JOSAA.9.002031.
«Numerical evaluation of the Hankel transform: remarks». Journal of the Optical Society of America A 10 (9): 1872. 1993. Bibcode:1993JOSAA..10.1872A. doi:10.1364/JOSAA.10.001872.
«Numerical evaluation of the zero-order Hankel transform using Filon quadrature philosophy». Applied Mathematics Letters 9 (5): 21-26. 1996. doi:10.1016/0893-9659(96)00067-5.
«Fast Hankel transform of nth order». J. Opt. Soc. Am. A 16 (10): 2581-2582. 1999. Bibcode:1999JOSAA..16.2581F. doi:10.1364/JOSAA.16.002581.
«Numerical evaluation of the Hankel transform». Comp. Phys. Comm. 116 (2–3): 278-294. 1999. Bibcode:1999CoPhC.116..278S. doi:10.1016/S0010-4655(98)00108-8.
«Algorithm 794: Numerical Hankel transform by the Fortran program HANKEL». ACM Trans. Math. Softw. 25 (2): 240-250. 1999. doi:10.1145/317275.317284.
«Fast Hankel transform by fast sine and cosine transforms: the Mellin connection». IEEE Trans. Signal Process. 48 (6): 1695-1701. 2000. Bibcode:2000ITSP...48.1695K. doi:10.1109/78.845927. Consultado el sì. (enlace roto disponible en Internet Archive; véase el historial, la primera versión y la última).
«Fast Hankel transform and its application for studying the propagation of cylindrical electromagnetic fields». Opt. Express 10 (12): 521-525. 2002. Bibcode:2002OExpr..10..521Z. doi:10.1364/oe.10.000521.
«Numerical evaluation of Hankel transforms for oscillating functions». J. Opt. Soc. Am. A 20 (4): 621-630. 2003. Bibcode:2003JOSAA..20..621M. doi:10.1364/JOSAA.20.000621.
«Fast Hankel transform of nth order with improved performance». J. Opt. Soc. Am. A 21 (9): 1811. 2004. Bibcode:2004JOSAA..21.1811P. doi:10.1364/JOSAA.21.001811.
«Computation of quasi-discrete Hankel transform of integer order for propagating optical wave fields». J. Opt. Soc. Am. A 21 (1): 53-58. 2004. Bibcode:2004JOSAA..21...53G. doi:10.1364/JOSAA.21.000053.
«The Zernike-Bessel representation and its application to Hankel transforms». J. Opt. Soc. Am. A 24 (6): 1609-1616. 2007. Bibcode:2007JOSAA..24.1609C. doi:10.1364/JOSAA.24.001609.