Los polinomios de Hermite fueron definidos por Pierre-Simon Laplace en 1810, [1] [2] aunque en una forma apenas reconocible, y estudiados en detalle por Pafnuty Chebyshev en 1859. [3] El trabajo de Chebyshev fue pasado por alto, y fueron nombrados más tarde en honor a Charles Hermite , quien escribió sobre los polinomios en 1864, describiéndolos como nuevos. [4] En consecuencia, no eran nuevos, aunque Hermite fue el primero en definir los polinomios multidimensionales.
Definición
Al igual que los demás polinomios ortogonales clásicos , los polinomios de Hermite se pueden definir a partir de varios puntos de partida diferentes. Teniendo en cuenta desde el principio que existen dos estandarizaciones diferentes de uso común, un método conveniente es el siguiente:
Los "polinomios de Hermite del probabilista" están dados por
Mientras que los "polinomios de Hermite del físico" se dan por
Estas ecuaciones tienen la forma de una fórmula de Rodrigues y también se pueden escribir como,
Las dos definiciones no son exactamente idénticas; cada una es una reescalación de la otra:
Éstas son secuencias polinomiales de Hermite de diferentes varianzas; consulte el material sobre varianzas a continuación.
Los primeros once polinomios de Hermite del probabilista son:
Los primeros once polinomios de Hermite del físico son:
Propiedades
El polinomio de Hermite de orden n es un polinomio de grado n . La versión probabilística He n tiene un coeficiente principal de 1, mientras que la versión física H n tiene un coeficiente principal de 2 n .
Simetría
De las fórmulas de Rodrigues dadas anteriormente, podemos ver que H n ( x ) y He n ( x ) son funciones pares o impares dependiendo de n :
Ortogonalidad
H n ( x ) y He n ( x ) sonpolinomios de grado n para n = 0, 1, 2, 3,... . Estos polinomios son ortogonales con respecto a la función de peso ( medida )
o
sea, tenemos
Los polinomios probabilistas son entonces ortogonales con respecto a la función de densidad de probabilidad normal estándar.
Lo completo
Los polinomios de Hermite (de probabilista o de físico) forman una base ortogonal del espacio de Hilbert de funciones que satisfacen
en el que el producto interno está dado por la integral
que incluye la función de peso gaussiana w ( x ) definida en la sección anterior.
Una base ortogonal para L 2 ( R , w ( x ) dx ) es un sistema ortogonal completo . Para un sistema ortogonal, la completitud es equivalente al hecho de que la función 0 es la única función f ∈ L 2 ( R , w ( x ) dx ) ortogonal a todas las funciones del sistema.
Dado que el espacio lineal de los polinomios de Hermite es el espacio de todos los polinomios, uno tiene que demostrar (en el caso de la física) que si f satisface
para cada n ≥ 0 , entonces f = 0 .
Una forma posible de hacer esto es apreciar que toda la función
se anula de manera idéntica. El hecho entonces de que F ( it ) = 0 para cada t real significa que la transformada de Fourier de f ( x ) e − x 2 es 0, por lo tanto f es 0 casi en todas partes. Las variantes de la prueba de completitud anterior se aplican a otros pesos con decaimiento exponencial.
En el caso de Hermite, también es posible probar una identidad explícita que implica completitud (ver sección sobre la relación de completitud a continuación).
Una formulación equivalente del hecho de que los polinomios de Hermite son una base ortogonal para L 2 ( R , w ( x ) dx ) consiste en introducir funciones de Hermite (ver más abajo), y decir que las funciones de Hermite son una base ortonormal para L 2 ( R ) .
Ecuación diferencial de Hermite
Los polinomios de Hermite del probabilista son soluciones de la ecuación diferencial
donde λ es una constante. Si se impone la condición de contorno de que u debe estar polinómicamente acotado en el infinito, la ecuación tiene soluciones solo si λ es un entero no negativo, y la solución está dada únicamente por , donde denota una constante.
Al reescribir la ecuación diferencial como un problema de valores propios ,
los polinomios de Hermite pueden entenderse como funciones propias del operador diferencial . Este problema de valores propios se denomina ecuación de Hermite , aunque el término también se utiliza para la ecuación estrechamente relacionada
cuya solución se da únicamente en términos de polinomios de Hermite del físico en la forma , donde denota una constante, después de imponer la condición de contorno de que u debe estar polinomialmente acotado en el infinito.
Las soluciones generales de las ecuaciones diferenciales de segundo orden anteriores son, de hecho, combinaciones lineales de polinomios de Hermite y funciones hipergeométricas confluentes de primera clase. Por ejemplo, para la ecuación de Hermite del físico,
la solución general toma la forma
donde y son constantes, son polinomios de Hermite del físico (de primera clase) y son funciones de Hermite del físico (de segunda clase). Las últimas funciones se representan de forma compacta como donde son Funciones hipergeométricas confluentes de primera clase . Los polinomios de Hermite convencionales también se pueden expresar en términos de funciones hipergeométricas confluentes, véase más abajo.
La secuencia de polinomios de Hermite del probabilista también satisface la relación de recurrencia
. Los coeficientes individuales están relacionados por la siguiente fórmula de recursión:
y a 0,0 = 1 , a 1,0 = 0 , a 1,1 = 1 .
Para los polinomios del físico, suponiendo que
tenemos
Los coeficientes individuales están relacionados por la siguiente fórmula de recursión:
y a 0,0 = 1 , a 1,0 = 0 , a 1,1 = 2 .
Los polinomios de Hermite constituyen una secuencia de Appell , es decir, son una secuencia polinómica que satisface la identidad
Una recurrencia integral que se deduce y demuestra en [6] es la siguiente:
De manera equivalente, mediante la expansión de Taylor ,
estas identidades umbrales son evidentes y se incluyen en la representación del operador diferencial que se detalla a continuación,
En consecuencia, para las derivadas m se cumplen las siguientes relaciones:
Estas últimas relaciones, junto con los polinomios iniciales H 0 ( x ) y H 1 ( x ) , se pueden utilizar en la práctica para calcular los polinomios rápidamente.
Los polinomios de Hermite del físico se pueden escribir explícitamente como
Estas dos ecuaciones se pueden combinar en una usando la función floor :
Los polinomios de Hermite del probabilista tienen fórmulas similares, que pueden obtenerse a partir de estas reemplazando la potencia de 2 x con la potencia correspondiente de √ 2 x y multiplicando la suma entera por 2 − norte/2:
Expresión explícita inversa
La inversa de las expresiones explícitas anteriores, es decir, aquellas para monomios en términos de polinomios de Hermite del probabilista son
Las expresiones correspondientes para los polinomios de Hermite del físico H se deducen directamente escalando adecuadamente esto: [7]
Esta igualdad es válida para todos los valores complejos de x y t , y se puede obtener escribiendo la expansión de Taylor en x de toda la función z → e − z 2 (en el caso del físico). También se puede derivar la función generadora (del físico) utilizando la fórmula integral de Cauchy para escribir los polinomios de Hermite como
Usando esto en la suma
se puede evaluar la integral restante usando el cálculo de residuos y llegar a la función generadora deseada.
Los momentos de la normal estándar (con valor esperado cero) se pueden leer directamente de la relación para índices pares:
donde (2 n − 1)!! es el factorial doble . Nótese que la expresión anterior es un caso especial de la representación de los polinomios de Hermite del probabilista como momentos:
Expansión asintótica
Asintóticamente, cuando n → ∞ , la expansión [8]
es válida. Para ciertos casos que involucran un rango más amplio de evaluación, es necesario incluir un factor para cambiar la amplitud:
que, utilizando la aproximación de Stirling , puede simplificarse aún más, en el límite, a
Una mejor aproximación, que tiene en cuenta la variación en la frecuencia, viene dada por
Una aproximación más fina, [9] que tiene en cuenta el espaciamiento desigual de los ceros cerca de los bordes, hace uso de la sustitución
con la que se tiene la aproximación uniforme
Se aplican aproximaciones similares para las regiones monótonas y de transición. En concreto, si
entonces
mientras que para con t complejo y acotado, la aproximación es
donde Ai es la función de Airy de primera clase.
Valores especiales
Los polinomios de Hermite del físico evaluados con argumento cero H n (0) se denominan números de Hermite .
que satisfacen la relación de recursión H n (0) = −2( n − 1) H n − 2 (0) .
En términos de los polinomios probabilistas esto se traduce en
Relaciones con otras funciones
Polinomios de Laguerre
Los polinomios de Hermite se pueden expresar como un caso especial de los polinomios de Laguerre :
Relación con funciones hipergeométricas confluentes
De manera similar a la expansión de Taylor, algunas funciones se pueden expresar como una suma infinita de polinomios de Hermite. Específicamente, si , entonces tiene una expansión en los polinomios de Hermite del físico. [10]
Dado tal , las sumas parciales de la expansión de Hermite de convergen a en la norma si y sólo si . [11]
Representación mediante operadores diferenciales
Los polinomios de Hermite del probabilista satisfacen la identidad donde D representa la diferenciación con respecto a x y la exponencial se interpreta desarrollándola como una serie de potencias . No hay cuestiones delicadas de convergencia de esta serie cuando opera sobre polinomios, ya que todos los términos, salvo un número finito, se anulan.
Dado que los coeficientes de la serie de potencias de la exponencial son bien conocidos y las derivadas de orden superior del monomio x n se pueden escribir explícitamente, esta representación del operador diferencial da lugar a una fórmula concreta para los coeficientes de H n que se puede utilizar para calcular rápidamente estos polinomios.
Como la expresión formal para la transformada de Weierstrass W es e D 2 , vemos que la transformada de Weierstrass de ( √ 2 ) n He n ( incógnita/√ 2 ) es x n . Esencialmente, la transformada de Weierstrass convierte una serie de polinomios de Hermite en una serie de Maclaurin correspondiente .
La existencia de alguna serie de potencias formal g ( D ) con coeficiente constante distinto de cero, tal que He n ( x ) = g ( D ) x n , es otro equivalente a la afirmación de que estos polinomios forman una sucesión de Appell . Puesto que son una sucesión de Appell, son a fortiori una sucesión de Sheffer .
Representación integral de contorno
De la representación de la función generadora anterior, vemos que los polinomios de Hermite tienen una representación en términos de una integral de contorno , como
con el contorno que rodea el origen.
Generalizaciones
Los polinomios de Hermite del probabilista definidos anteriormente son ortogonales con respecto a la distribución de probabilidad normal estándar, cuya función de densidad es
que tiene un valor esperado 0 y una varianza 1.
Escalando, se puede hablar análogamente de polinomios de Hermite generalizados [12]
de varianza α , donde α es cualquier número positivo. Éstos son entonces ortogonales con respecto a la distribución de probabilidad normal cuya función de densidad es
Vienen dados por
Ahora bien, si
la secuencia polinómica cuyo término n es
se denomina composición umbral de las dos secuencias polinómicas, se puede demostrar que satisface las identidades
y
La última identidad se expresa diciendo que esta familia parametrizada de secuencias polinómicas se conoce como secuencia cruzada. (Véase la sección anterior sobre secuencias de Appell y sobre la representación del operador diferencial, que conduce a una fácil derivación de la misma. Esta identidad de tipo binomial , para α = β = 1/2, ya se ha encontrado en la sección anterior sobre relaciones de #Recursión).
"Varianza negativa"
Como las secuencias polinómicas forman un grupo bajo la operación de composición umbral , se puede denotar por
la secuencia que es inversa a la denotada de manera similar, pero sin el signo menos, y así hablar de polinomios de Hermite de varianza negativa. Para α > 0 , los coeficientes de son simplemente los valores absolutos de los coeficientes correspondientes de .
Estos surgen como momentos de distribuciones de probabilidad normales: El momento n de la distribución normal con valor esperado μ y varianza σ 2 es
donde X es una variable aleatoria con la distribución normal especificada. Un caso especial de la identidad entre secuencias dice entonces que
Funciones de Hermite
Definición
Se pueden definir las funciones de Hermite (a menudo llamadas funciones Hermite-Gaussianas) a partir de los polinomios del físico:
Así,
Dado que estas funciones contienen la raíz cuadrada de la función de peso y se han escalado adecuadamente , son ortonormales
y forman una base ortonormal de L 2 ( R ) . Este hecho es equivalente al enunciado correspondiente para los polinomios de Hermite (ver arriba).
Las funciones de Hermite satisfacen la ecuación diferencial
Esta ecuación es equivalente a la ecuación de Schrödinger para un oscilador armónico en mecánica cuántica, por lo que estas funciones son las funciones propias .
Relación de recursión
Siguiendo las relaciones de recursión de los polinomios de Hermite, las funciones de Hermite obedecen
y
Extender la primera relación a las derivadas m arbitrarias para cualquier entero positivo m conduce a
Esta fórmula se puede utilizar en conexión con las relaciones de recurrencia para He n y ψ n para calcular cualquier derivada de las funciones de Hermite de manera eficiente.
Desigualdad de Cramer
Para x real , las funciones de Hermite satisfacen el siguiente límite debido a Harald Cramér [13] [14] y Jack Indritz: [15]
Funciones de Hermite como funciones propias de la transformada de Fourier
Las funciones de Hermite ψ n ( x ) son un conjunto de funciones propias de la transformada continua de Fourier F . Para ver esto, tome la versión física de la función generadora y multiplíquela por e − 1/2 x 2 . Esto da
La transformada de Fourier del lado izquierdo está dada por
La transformada de Fourier del lado derecho está dada por
Igualando potencias iguales de t en las versiones transformadas de los lados izquierdo y derecho finalmente obtenemos
Las funciones de Hermite ψ n ( x ) son, por tanto, una base ortonormal de L 2 ( R ) , que diagonaliza el operador de transformada de Fourier . [16]
Distribuciones de Wigner de funciones de Hermite
La función de distribución de Wigner de la función de Hermite de orden n está relacionada con el polinomio de Laguerre de orden n . Los polinomios de Laguerre
conducen a las funciones de Laguerre del oscilador
Para todos los enteros naturales n , es sencillo ver [17] que
donde la distribución de Wigner de una función x ∈ L 2 ( R , C ) se define como
Este es un resultado fundamental para el oscilador armónico cuántico descubierto por Hip Groenewold en 1946 en su tesis doctoral. [18] Es el paradigma estándar de la mecánica cuántica en el espacio de fases .
Existen más relaciones entre las dos familias de polinomios.
Interpretación combinatoria de coeficientes
En el polinomio de Hermite He n ( x ) de varianza 1, el valor absoluto del coeficiente de x k es el número de particiones (no ordenadas) de un conjunto de n elementos en k singletons yn - k/2 pares (desordenados). Equivalentemente, es el número de involuciones de unconjunto de n elementos con precisamente k puntos fijos, o en otras palabras, el número de emparejamientos en el grafo completo en n vértices que dejan k vértices sin cubrir (de hecho, los polinomios de Hermite son los polinomios de emparejamiento de estos grafos). La suma de los valores absolutos de los coeficientes da el número total de particiones en singletons y pares, los llamados números de teléfono .
1, 1, 2, 4, 10, 26, 76, 232, 764, 2620, 9496,... (secuencia A000085 en la OEIS ).
Esta interpretación combinatoria se puede relacionar con polinomios de Bell exponenciales completos como
donde x i = 0 para todo i > 2 .
Estos números también pueden expresarse como un valor especial de los polinomios de Hermite: [19]
Además, la siguiente identidad de completitud para las funciones de Hermite anteriores se cumple en el sentido de las distribuciones :
donde δ es la función delta de Dirac , ψ n las funciones de Hermite, y δ ( x − y ) representa la medida de Lebesgue en la línea y = x en R 2 , normalizada de modo que su proyección en el eje horizontal es la medida de Lebesgue habitual.
Esta identidad distribucional sigue a Wiener (1958) al tomar u → 1 en la fórmula de Mehler , válida cuando −1 < u < 1 :
que a menudo se enuncia de manera equivalente como un núcleo separable, [20] [21]
La función ( x , y ) → E ( x , y ; u ) es la densidad de probabilidad gaussiana bivariada en R 2 , que está, cuando u está cerca de 1, muy concentrada alrededor de la línea y = x , y muy dispersa en esa línea. De ello se deduce que
cuando f y g son continuas y están soportadas de forma compacta.
Esto da como resultado que f se puede expresar en funciones de Hermite como la suma de una serie de vectores en L 2 ( R ) , es decir,
El polinomio de Hermite se representa entonces como
Con esta representación para H n ( x ) y H n ( y ) , es evidente que
y esto produce la resolución deseada del resultado de identidad, utilizando nuevamente la transformada de Fourier de núcleos gaussianos bajo la sustitución
^ Laplace (1811). "Mémoire sur les intégrales définies et leur application aux probabilités, et spécialement a la recherche du milieu qu'il faut choisir entre les resultats des observes" [Memoria sobre integrales definidas y su aplicación a las probabilidades, y especialmente a la búsqueda de la media que debe ser elegido entre los resultados de las observaciones]. Mémoires de la Classe des Sciences Mathématiques et Physiques de l'Institut Impérial de France (en francés). 11 : 297–347.
^ Laplace, P.-S. (1812), Théorie analytique des probabilités [ Teoría analítica de la probabilidad ], vol. 2, págs. 194-203Recogido en Œuvres complètes VII.
^ Tchébychef, P. (1860). "Sur le développement des fonctions à une seule variable" [Sobre el desarrollo de funciones de una sola variable]. Bulletin de l'Académie impériale des sciences de St.-Pétersbourg (en francés). 1 : 193–200.Recopilado en Obras I, 501–508.
^ Hermita, C. (1864). "Sur un nouveau développement en série de fonctions" [Sobre un nuevo desarrollo en serie de funciones]. CR Acad. Ciencia. París (en francés). 58 : 93–100, 266–273.Recopilado en Œuvres II , 293–308.
^ Tom H. Koornwinder, Roderick SC Wong y Roelof Koekoek et al. (2010) y Abramowitz y Stegun .
^ Hurtado Benavides, Miguel Ángel. (2020). De las sumas de potencias a las sucesiones de Appell y su caracterización a través de funcionales. [Tesis de maestría]. Universidad Sergio Arboleda.
^ "18. Polinomios ortogonales, polinomios ortogonales clásicos, sumas". Biblioteca digital de funciones matemáticas . Instituto Nacional de Estándares y Tecnología . Consultado el 30 de enero de 2015 .
^ Abramowitz y Stegun 1983, pág. 508–510, 13.6.38 y 13.5.16.
^ Szegő 1955, pág. 201
^ "Tutorial de MATHEMATICA, parte 2.5: Expansión de Hermite". www.cfm.brown.edu . Consultado el 24 de diciembre de 2023 .
^ Askey, Richard; Wainger, Stephen (1965). "Convergencia media de expansiones en series de Laguerre y Hermite". American Journal of Mathematics . 87 (3): 695–708. doi :10.2307/2373069. ISSN 0002-9327.
^ Roman, Steven (1984), El cálculo umbral , Matemáticas puras y aplicadas, vol. 111 (1.ª ed.), Academic Press, págs. 87-93, ISBN978-0-12-594380-2
^ Erdélyi y col. 1955, pág. 207.
^ Szegő 1955.
^ Indritz, Jack (1961), "Una desigualdad para polinomios de Hermite", Actas de la American Mathematical Society , 12 (6): 981–983, doi : 10.1090/S0002-9939-1961-0132852-2 , MR 0132852
^ En este caso, utilizamos la versión unitaria de la transformada de Fourier, por lo que los valores propios son (− i ) n . La resolución resultante de la identidad sirve entonces para definir potencias, incluidas las fraccionarias, de la transformada de Fourier, es decir, una generalización de la transformada de Fourier fraccionaria , en efecto un núcleo de Mehler .
^ Folland, GB (1989), Análisis armónico en el espacio de fases , Anales de estudios matemáticos, vol. 122, Princeton University Press, ISBN978-0-691-08528-9
^ Groenewold, HJ (1946). "Sobre los principios de la mecánica cuántica elemental". Physica . 12 (7): 405–460. Bibcode :1946Phy....12..405G. doi :10.1016/S0031-8914(46)80059-4.
^ Banderier, Cirilo; Bousquet-Mélou, Mireille ; Denise, Alain; Flajolet, Philippe ; Gardy, Danièle; Gouyou-Beauchamps, Dominique (2002), "Generación de funciones para generar árboles", Matemáticas discretas , 246 (1–3): 29–55, arXiv : math/0411250 , doi :10.1016/S0012-365X(01)00250-3 , SEÑOR 1884885, S2CID 14804110
^ Mehler, FG (1866), "Ueber die Entwicklung einer Function von beliebig vielen Variabeln nach Laplaceschen Functionen höherer Ordnung" [Sobre el desarrollo de una función de muchas variables arbitrarias según funciones de Laplace de orden superior], Journal für die Reine und Angewandte Mathematik (en alemán) (66): 161–176, ISSN 0075-4102, ERAM 066.1720cj. Véase pág. 174, ecuación (18) y pág. 173, ecuación (13).
^ Erdélyi y col. 1955, pág. 194, 10.13 (22).
Referencias
Abramowitz, Milton ; Stegun, Irene Ann , eds. (1983) [junio de 1964]. "Capítulo 22". Manual de funciones matemáticas con fórmulas, gráficos y tablas matemáticas . Serie de Matemáticas Aplicadas. Vol. 55 (Novena reimpresión con correcciones adicionales de la décima impresión original con correcciones (diciembre de 1972); primera ed.). Washington DC; Nueva York: Departamento de Comercio de los Estados Unidos, Oficina Nacional de Normas; Dover Publications. pág. 773. ISBN 978-0-486-61272-0. LCCN 64-60036. MR 0167642. LCCN 65-12253.
Laplace, PS (1810), "Mémoire sur les intégrales définies et leur application aux probabilités, et spécialement a la recherche du milieu qu'il faut choisir entre les résultats des observes", Mémoires de l'Académie des Sciences : 279–347Obras completas 12, pp.357-412, traducción al inglés Archivado el 4 de marzo de 2016 en Wayback Machine .
Shohat, JA; Hille, Einar; Walsh, Joseph L. (1940), Una bibliografía sobre polinomios ortogonales , Boletín del Consejo Nacional de Investigación, Washington DC: Academia Nacional de Ciencias- 2000 referencias de Bibliografía sobre polinomios de Hermite.
Szegő, Gábor (1955) [1939], Polinomios ortogonales , Colloquium Publications, vol. 23 (4.ª ed.), American Mathematical Society, ISBN 978-0-8218-1023-1
Temme, Nico (1996), Funciones especiales: una introducción a las funciones clásicas de la física matemática , Nueva York: Wiley, ISBN 978-0-471-11313-3
Wiener, Norbert (1958) [1933], La integral de Fourier y algunas de sus aplicaciones (edición revisada), Nueva York: Dover Publications, ISBN 0-486-60272-9
Biblioteca científica GNU: incluye la versión C de polinomios de Hermite, funciones, sus derivadas y ceros (consulte también Biblioteca científica GNU )