Polinomios de Hermite

Sucesión polinómica

En matemáticas , los polinomios de Hermite son una secuencia polinomial ortogonal clásica .

Los polinomios surgen en:

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 Él norte ( incógnita ) = ( 1 ) norte mi incógnita 2 2 d norte d incógnita norte mi incógnita 2 2 , {\displaystyle \operatorname {He} _{n}(x)=(-1)^{n}e^{\frac {x^{2}}{2}}{\frac {d^{n}}{dx^{n}}}e^{-{\frac {x^{2}}{2}}},}
  • Mientras que los "polinomios de Hermite del físico" se dan por yo norte ( incógnita ) = ( 1 ) norte mi incógnita 2 d norte d incógnita norte mi incógnita 2 . {\displaystyle H_{n}(x)=(-1)^{n}e^{x^{2}}{\frac {d^{n}}{dx^{n}}}e^{-x^{2}}.}

Estas ecuaciones tienen la forma de una fórmula de Rodrigues y también se pueden escribir como, Él norte ( incógnita ) = ( incógnita d d incógnita ) norte 1 , yo norte ( incógnita ) = ( 2 incógnita d d incógnita ) norte 1. {\displaystyle \operatorname {He} _{n}(x)=\left(x-{\frac {d}{dx}}\right)^{n}\cdot 1,\quad H_{n}(x)=\left(2x-{\frac {d}{dx}}\right)^{n}\cdot 1.}

Las dos definiciones no son exactamente idénticas; cada una es una reescalación de la otra: yo norte ( incógnita ) = 2 norte 2 Él norte ( 2 incógnita ) , Él norte ( incógnita ) = 2 norte 2 yo norte ( incógnita 2 ) . {\displaystyle H_{n}(x)=2^{\frac {n}{2}}\operatorname {He} _{n}\left({\sqrt {2}}\,x\right),\quad \operatorname {He} _{n}(x)=2^{-{\frac {n}{2}}}H_{n}\left({\frac {x}{\sqrt {2}}}\right).}

Éstas son secuencias polinomiales de Hermite de diferentes varianzas; consulte el material sobre varianzas a continuación.

La notación He y H es la utilizada en las referencias estándar. [5] Los polinomios He n a veces se denotan por H n , especialmente en teoría de probabilidad, porque es la función de densidad de probabilidad para la distribución normal con valor esperado 0 y desviación estándar 1. 1 2 π mi incógnita 2 2 {\displaystyle {\frac {1}{\sqrt {2\pi }}}e^{-{\frac {x^{2}}{2}}}}

Los primeros seis polinomios de Hermite del probabilista He n ( x )
Los primeros seis polinomios de Hermite (del físico) H n ( x )
  • Los primeros once polinomios de Hermite del probabilista son: Él 0 ( incógnita ) = 1 , Él 1 ( incógnita ) = incógnita , Él 2 ( incógnita ) = incógnita 2 1 , Él 3 ( incógnita ) = incógnita 3 3 incógnita , Él 4 ( incógnita ) = incógnita 4 6 incógnita 2 + 3 , Él 5 ( incógnita ) = incógnita 5 10 incógnita 3 + 15 incógnita , Él 6 ( incógnita ) = incógnita 6 15 incógnita 4 + 45 incógnita 2 15 , Él 7 ( incógnita ) = incógnita 7 21 incógnita 5 + 105 incógnita 3 105 incógnita , Él 8 ( incógnita ) = incógnita 8 28 incógnita 6 + 210 incógnita 4 420 incógnita 2 + 105 , Él 9 ( incógnita ) = incógnita 9 36 incógnita 7 + 378 incógnita 5 1260 incógnita 3 + 945 incógnita , Él 10 ( incógnita ) = incógnita 10 45 incógnita 8 + 630 incógnita 6 3150 incógnita 4 + 4725 incógnita 2 945. {\displaystyle {\begin{aligned}\operatorname {He} _{0}(x)&=1,\\\operatorname {He} _{1}(x)&=x,\\\operatorname {He} _{2}(x)&=x^{2}-1,\\\operatorname {He} _{3}(x)&=x^{3}-3x,\\\operatorname {He} _{4}(x)&=x^{4}-6x^{2}+3,\\\operatorname {He} _{5}(x)&=x^{5}-10x^{3}+15x,\\\operatorname {He} _{6}(x)&=x^{6}-15x^{4}+45x^{2}-15,\\\operatorname {He} _{7}(x)&=x^{7}-21x^{5}+105x^{3}-105x,\\\nombre del operador {He} _{8}(x)&=x^{8}-28x^{6}+210x^{4}-420x^{2}+105,\\\nombre del operador {He} _{9}(x)&=x^{9}-36x^{7}+378x^{5}-1260x^{3}+945x,\\\nombre del operador {He} _{10}(x)&=x^{10}-45x^{8}+630x^{6}-3150x^{4}+4725x^{2}-945.\end{aligned}}}
  • Los primeros once polinomios de Hermite del físico son: yo 0 ( incógnita ) = 1 , yo 1 ( incógnita ) = 2 incógnita , yo 2 ( incógnita ) = 4 incógnita 2 2 , yo 3 ( incógnita ) = 8 incógnita 3 12 incógnita , yo 4 ( incógnita ) = 16 incógnita 4 48 incógnita 2 + 12 , yo 5 ( incógnita ) = 32 incógnita 5 160 incógnita 3 + 120 incógnita , yo 6 ( incógnita ) = 64 incógnita 6 480 incógnita 4 + 720 incógnita 2 120 , yo 7 ( incógnita ) = 128 incógnita 7 1344 incógnita 5 + 3360 incógnita 3 1680 incógnita , yo 8 ( incógnita ) = 256 incógnita 8 3584 incógnita 6 + 13440 incógnita 4 13440 incógnita 2 + 1680 , yo 9 ( incógnita ) = 512 incógnita 9 9216 incógnita 7 + 48384 incógnita 5 80640 incógnita 3 + 30240 incógnita , yo 10 ( incógnita ) = 1024 incógnita 10 23040 incógnita 8 + 161280 incógnita 6 403200 incógnita 4 + 302400 incógnita 2 30240. {\displaystyle {\begin{aligned}H_{0}(x)&=1,\\H_{1}(x)&=2x,\\H_{2}(x)&=4x^{2}-2,\\H_{3}(x)&=8x^{3}-12x,\\H_{4}(x)&=16x^{4}-48x^{2}+12,\\H_{5}(x)&=32x^{5}-160x^{3}+120x,\\H_{6}(x)&=64x^{6}-480x^{4}+720x^{2}-120,\\H_{7}(x)&=128x^{7}-1344x^{5}+3360x^{3}-1680x,\\H_{8}(x)&=256x^{8}-3584x^{6}+13440x^{4}-13440x^{2}+1680,\\H_{9}(x)&=512x^{9}-9216x^{7}+48384x^{5}-80640x^{3}+30240x,\\H_{10}(x)&=1024x^{10}-23040x^{8}+161280x^{6}-403200x^{4}+302400x^{2}-30240.\end{aligned}}}

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 : H n ( x ) = ( 1 ) n H n ( x ) , He n ( x ) = ( 1 ) n He n ( x ) . {\displaystyle H_{n}(-x)=(-1)^{n}H_{n}(x),\quad \operatorname {He} _{n}(-x)=(-1)^{n}\operatorname {He} _{n}(x).}

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 w ( x ) = e x 2 2 ( for  He ) {\displaystyle w(x)=e^{-{\frac {x^{2}}{2}}}\quad ({\text{for }}\operatorname {He} )} w ( x ) = e x 2 ( for  H ) , {\displaystyle w(x)=e^{-x^{2}}\quad ({\text{for }}H),} H m ( x ) H n ( x ) w ( x ) d x = 0 for all  m n . {\displaystyle \int _{-\infty }^{\infty }H_{m}(x)H_{n}(x)\,w(x)\,dx=0\quad {\text{for all }}m\neq n.}

Además, ¿ dónde está el delta de Kronecker ? H m ( x ) H n ( x ) e x 2 d x = π 2 n n ! δ n m , {\displaystyle \int _{-\infty }^{\infty }H_{m}(x)H_{n}(x)\,e^{-x^{2}}\,dx={\sqrt {\pi }}\,2^{n}n!\,\delta _{nm},} He m ( x ) He n ( x ) e x 2 2 d x = 2 π n ! δ n m , {\displaystyle \int _{-\infty }^{\infty }\operatorname {He} _{m}(x)\operatorname {He} _{n}(x)\,e^{-{\frac {x^{2}}{2}}}\,dx={\sqrt {2\pi }}\,n!\,\delta _{nm},} δ n m {\displaystyle \delta _{nm}}

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. | f ( x ) | 2 w ( x ) d x < , {\displaystyle \int _{-\infty }^{\infty }{\bigl |}f(x){\bigr |}^{2}\,w(x)\,dx<\infty ,} f , g = f ( x ) g ( x ) ¯ w ( x ) d x {\displaystyle \langle f,g\rangle =\int _{-\infty }^{\infty }f(x){\overline {g(x)}}\,w(x)\,dx}

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 fL 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 . f ( x ) x n e x 2 d x = 0 {\displaystyle \int _{-\infty }^{\infty }f(x)x^{n}e^{-x^{2}}\,dx=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. F ( z ) = f ( x ) e z x x 2 d x = n = 0 z n n ! f ( x ) x n e x 2 d x = 0 {\displaystyle F(z)=\int _{-\infty }^{\infty }f(x)e^{zx-x^{2}}\,dx=\sum _{n=0}^{\infty }{\frac {z^{n}}{n!}}\int f(x)x^{n}e^{-x^{2}}\,dx=0}

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. ( e 1 2 x 2 u ) + λ e 1 2 x 2 u = 0 , {\displaystyle \left(e^{-{\frac {1}{2}}x^{2}}u'\right)'+\lambda e^{-{\frac {1}{2}}x^{2}}u=0,} u ( x ) = C 1 He λ ( x ) {\displaystyle u(x)=C_{1}\operatorname {He} _{\lambda }(x)} C 1 {\displaystyle C_{1}}

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. L [ u ] = u x u = λ u , {\displaystyle L[u]=u''-xu'=-\lambda u,} He λ ( x ) {\displaystyle \operatorname {He} _{\lambda }(x)} L [ u ] {\displaystyle L[u]} u 2 x u = 2 λ u . {\displaystyle u''-2xu'=-2\lambda u.} u ( x ) = C 1 H λ ( x ) {\displaystyle u(x)=C_{1}H_{\lambda }(x)} C 1 {\displaystyle C_{1}}

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. u 2 x u + 2 λ u = 0 , {\displaystyle u''-2xu'+2\lambda u=0,} u ( x ) = C 1 H λ ( x ) + C 2 h λ ( x ) , {\displaystyle u(x)=C_{1}H_{\lambda }(x)+C_{2}h_{\lambda }(x),} C 1 {\displaystyle C_{1}} C 2 {\displaystyle C_{2}} H λ ( x ) {\displaystyle H_{\lambda }(x)} h λ ( x ) {\displaystyle h_{\lambda }(x)} h λ ( x ) = 1 F 1 ( λ 2 ; 1 2 ; x 2 ) {\displaystyle h_{\lambda }(x)={}_{1}F_{1}(-{\tfrac {\lambda }{2}};{\tfrac {1}{2}};x^{2})} 1 F 1 ( a ; b ; z ) {\displaystyle {}_{1}F_{1}(a;b;z)}

Con condiciones de contorno más generales , los polinomios de Hermite se pueden generalizar para obtener funciones analíticas más generales para λ de valor complejo . También es posible una fórmula explícita de polinomios de Hermite en términos de integrales de contorno (Courant y Hilbert 1989).

Relación de recurrencia

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 . He n + 1 ( x ) = x He n ( x ) He n ( x ) . {\displaystyle \operatorname {He} _{n+1}(x)=x\operatorname {He} _{n}(x)-\operatorname {He} _{n}'(x).} a n + 1 , k = { ( k + 1 ) a n , k + 1 k = 0 , a n , k 1 ( k + 1 ) a n , k + 1 k > 0 , {\displaystyle a_{n+1,k}={\begin{cases}-(k+1)a_{n,k+1}&k=0,\\a_{n,k-1}-(k+1)a_{n,k+1}&k>0,\end{cases}}}

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 . H n ( x ) = k = 0 n a n , k x k , {\displaystyle H_{n}(x)=\sum _{k=0}^{n}a_{n,k}x^{k},} H n + 1 ( x ) = 2 x H n ( x ) H n ( x ) . {\displaystyle H_{n+1}(x)=2xH_{n}(x)-H_{n}'(x).} a n + 1 , k = { a n , k + 1 k = 0 , 2 a n , k 1 ( k + 1 ) a n , k + 1 k > 0 , {\displaystyle a_{n+1,k}={\begin{cases}-a_{n,k+1}&k=0,\\2a_{n,k-1}-(k+1)a_{n,k+1}&k>0,\end{cases}}}

Los polinomios de Hermite constituyen una secuencia de Appell , es decir, son una secuencia polinómica que satisface la identidad He n ( x ) = n He n 1 ( x ) , H n ( x ) = 2 n H n 1 ( x ) . {\displaystyle {\begin{aligned}\operatorname {He} _{n}'(x)&=n\operatorname {He} _{n-1}(x),\\H_{n}'(x)&=2nH_{n-1}(x).\end{aligned}}}

Una recurrencia integral que se deduce y demuestra en [6] es la siguiente: He n + 1 ( x ) = ( n + 1 ) 0 x He n ( t ) d t H e n ( 0 ) , {\displaystyle \operatorname {He} _{n+1}(x)=(n+1)\int _{0}^{x}\operatorname {He} _{n}(t)dt-He'_{n}(0),}

H n + 1 ( x ) = 2 ( n + 1 ) 0 x H n ( t ) d t H n ( 0 ) . {\displaystyle H_{n+1}(x)=2(n+1)\int _{0}^{x}H_{n}(t)dt-H'_{n}(0).}

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, He n ( x + y ) = k = 0 n ( n k ) x n k He k ( y ) = 2 n 2 k = 0 n ( n k ) He n k ( x 2 ) He k ( y 2 ) , H n ( x + y ) = k = 0 n ( n k ) H k ( x ) ( 2 y ) n k = 2 n 2 k = 0 n ( n k ) H n k ( x 2 ) H k ( y 2 ) . {\displaystyle {\begin{aligned}\operatorname {He} _{n}(x+y)&=\sum _{k=0}^{n}{\binom {n}{k}}x^{n-k}\operatorname {He} _{k}(y)&&=2^{-{\frac {n}{2}}}\sum _{k=0}^{n}{\binom {n}{k}}\operatorname {He} _{n-k}\left(x{\sqrt {2}}\right)\operatorname {He} _{k}\left(y{\sqrt {2}}\right),\\H_{n}(x+y)&=\sum _{k=0}^{n}{\binom {n}{k}}H_{k}(x)(2y)^{n-k}&&=2^{-{\frac {n}{2}}}\cdot \sum _{k=0}^{n}{\binom {n}{k}}H_{n-k}\left(x{\sqrt {2}}\right)H_{k}\left(y{\sqrt {2}}\right).\end{aligned}}} He n ( x ) = e D 2 2 x n , H n ( x ) = 2 n e D 2 4 x n . {\displaystyle {\begin{aligned}\operatorname {He} _{n}(x)&=e^{-{\frac {D^{2}}{2}}}x^{n},\\H_{n}(x)&=2^{n}e^{-{\frac {D^{2}}{4}}}x^{n}.\end{aligned}}}

En consecuencia, para las derivadas m se cumplen las siguientes relaciones: He n ( m ) ( x ) = n ! ( n m ) ! He n m ( x ) = m ! ( n m ) He n m ( x ) , H n ( m ) ( x ) = 2 m n ! ( n m ) ! H n m ( x ) = 2 m m ! ( n m ) H n m ( x ) . {\displaystyle {\begin{aligned}\operatorname {He} _{n}^{(m)}(x)&={\frac {n!}{(n-m)!}}\operatorname {He} _{n-m}(x)&&=m!{\binom {n}{m}}\operatorname {He} _{n-m}(x),\\H_{n}^{(m)}(x)&=2^{m}{\frac {n!}{(n-m)!}}H_{n-m}(x)&&=2^{m}m!{\binom {n}{m}}H_{n-m}(x).\end{aligned}}}

De ello se deduce que los polinomios de Hermite también satisfacen la relación de recurrencia. He n + 1 ( x ) = x He n ( x ) n He n 1 ( x ) , H n + 1 ( x ) = 2 x H n ( x ) 2 n H n 1 ( x ) . {\displaystyle {\begin{aligned}\operatorname {He} _{n+1}(x)&=x\operatorname {He} _{n}(x)-n\operatorname {He} _{n-1}(x),\\H_{n+1}(x)&=2xH_{n}(x)-2nH_{n-1}(x).\end{aligned}}}

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.

Las desigualdades de Turán son H n ( x ) 2 H n 1 ( x ) H n + 1 ( x ) = ( n 1 ) ! i = 0 n 1 2 n i i ! H i ( x ) 2 > 0. {\displaystyle {\mathit {H}}_{n}(x)^{2}-{\mathit {H}}_{n-1}(x){\mathit {H}}_{n+1}(x)=(n-1)!\sum _{i=0}^{n-1}{\frac {2^{n-i}}{i!}}{\mathit {H}}_{i}(x)^{2}>0.}

Además, se cumple el siguiente teorema de multiplicación : H n ( γ x ) = i = 0 n 2 γ n 2 i ( γ 2 1 ) i ( n 2 i ) ( 2 i ) ! i ! H n 2 i ( x ) , He n ( γ x ) = i = 0 n 2 γ n 2 i ( γ 2 1 ) i ( n 2 i ) ( 2 i ) ! i ! 2 i He n 2 i ( x ) . {\displaystyle {\begin{aligned}H_{n}(\gamma x)&=\sum _{i=0}^{\left\lfloor {\tfrac {n}{2}}\right\rfloor }\gamma ^{n-2i}(\gamma ^{2}-1)^{i}{\binom {n}{2i}}{\frac {(2i)!}{i!}}H_{n-2i}(x),\\\operatorname {He} _{n}(\gamma x)&=\sum _{i=0}^{\left\lfloor {\tfrac {n}{2}}\right\rfloor }\gamma ^{n-2i}(\gamma ^{2}-1)^{i}{\binom {n}{2i}}{\frac {(2i)!}{i!}}2^{-i}\operatorname {He} _{n-2i}(x).\end{aligned}}}

Expresión explícita

Los polinomios de Hermite del físico se pueden escribir explícitamente como H n ( x ) = { n ! l = 0 n 2 ( 1 ) n 2 l ( 2 l ) ! ( n 2 l ) ! ( 2 x ) 2 l for even  n , n ! l = 0 n 1 2 ( 1 ) n 1 2 l ( 2 l + 1 ) ! ( n 1 2 l ) ! ( 2 x ) 2 l + 1 for odd  n . {\displaystyle H_{n}(x)={\begin{cases}\displaystyle n!\sum _{l=0}^{\frac {n}{2}}{\frac {(-1)^{{\tfrac {n}{2}}-l}}{(2l)!\left({\tfrac {n}{2}}-l\right)!}}(2x)^{2l}&{\text{for even }}n,\\\displaystyle n!\sum _{l=0}^{\frac {n-1}{2}}{\frac {(-1)^{{\frac {n-1}{2}}-l}}{(2l+1)!\left({\frac {n-1}{2}}-l\right)!}}(2x)^{2l+1}&{\text{for odd }}n.\end{cases}}}

Estas dos ecuaciones se pueden combinar en una usando la función floor : H n ( x ) = n ! m = 0 n 2 ( 1 ) m m ! ( n 2 m ) ! ( 2 x ) n 2 m . {\displaystyle H_{n}(x)=n!\sum _{m=0}^{\left\lfloor {\tfrac {n}{2}}\right\rfloor }{\frac {(-1)^{m}}{m!(n-2m)!}}(2x)^{n-2m}.}

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: He n ( x ) = n ! m = 0 n 2 ( 1 ) m m ! ( n 2 m ) ! x n 2 m 2 m . {\displaystyle \operatorname {He} _{n}(x)=n!\sum _{m=0}^{\left\lfloor {\tfrac {n}{2}}\right\rfloor }{\frac {(-1)^{m}}{m!(n-2m)!}}{\frac {x^{n-2m}}{2^{m}}}.}

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 x n = n ! m = 0 n 2 1 2 m m ! ( n 2 m ) ! He n 2 m ( x ) . {\displaystyle x^{n}=n!\sum _{m=0}^{\left\lfloor {\tfrac {n}{2}}\right\rfloor }{\frac {1}{2^{m}m!(n-2m)!}}\operatorname {He} _{n-2m}(x).}

Las expresiones correspondientes para los polinomios de Hermite del físico H se deducen directamente escalando adecuadamente esto: [7] x n = n ! 2 n m = 0 n 2 1 m ! ( n 2 m ) ! H n 2 m ( x ) . {\displaystyle x^{n}={\frac {n!}{2^{n}}}\sum _{m=0}^{\left\lfloor {\tfrac {n}{2}}\right\rfloor }{\frac {1}{m!(n-2m)!}}H_{n-2m}(x).}

Función generadora

Los polinomios de Hermite están dados por la función generadora exponencial e x t 1 2 t 2 = n = 0 He n ( x ) t n n ! , e 2 x t t 2 = n = 0 H n ( x ) t n n ! . {\displaystyle {\begin{aligned}e^{xt-{\frac {1}{2}}t^{2}}&=\sum _{n=0}^{\infty }\operatorname {He} _{n}(x){\frac {t^{n}}{n!}},\\e^{2xt-t^{2}}&=\sum _{n=0}^{\infty }H_{n}(x){\frac {t^{n}}{n!}}.\end{aligned}}}

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 ze 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 H n ( x ) = ( 1 ) n e x 2 d n d x n e x 2 = ( 1 ) n e x 2 n ! 2 π i γ e z 2 ( z x ) n + 1 d z . {\displaystyle H_{n}(x)=(-1)^{n}e^{x^{2}}{\frac {d^{n}}{dx^{n}}}e^{-x^{2}}=(-1)^{n}e^{x^{2}}{\frac {n!}{2\pi i}}\oint _{\gamma }{\frac {e^{-z^{2}}}{(z-x)^{n+1}}}\,dz.}

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. n = 0 H n ( x ) t n n ! , {\displaystyle \sum _{n=0}^{\infty }H_{n}(x){\frac {t^{n}}{n!}},}

Valores esperados

Si X es una variable aleatoria con una distribución normal con desviación estándar 1 y valor esperado μ , entonces E [ He n ( X ) ] = μ n . {\displaystyle \operatorname {\mathbb {E} } \left[\operatorname {He} _{n}(X)\right]=\mu ^{n}.}

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: E [ X 2 n ] = ( 1 ) n He 2 n ( 0 ) = ( 2 n 1 ) ! ! , {\displaystyle \operatorname {\mathbb {E} } \left[X^{2n}\right]=(-1)^{n}\operatorname {He} _{2n}(0)=(2n-1)!!,} He n ( x ) = 1 2 π ( x + i y ) n e y 2 2 d y . {\displaystyle \operatorname {He} _{n}(x)={\frac {1}{\sqrt {2\pi }}}\int _{-\infty }^{\infty }(x+iy)^{n}e^{-{\frac {y^{2}}{2}}}\,dy.}

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 e x 2 2 H n ( x ) 2 n π Γ ( n + 1 2 ) cos ( x 2 n n π 2 ) {\displaystyle e^{-{\frac {x^{2}}{2}}}\cdot H_{n}(x)\sim {\frac {2^{n}}{\sqrt {\pi }}}\Gamma \left({\frac {n+1}{2}}\right)\cos \left(x{\sqrt {2n}}-{\frac {n\pi }{2}}\right)} e x 2 2 H n ( x ) 2 n π Γ ( n + 1 2 ) cos ( x 2 n n π 2 ) ( 1 x 2 2 n + 1 ) 1 4 = 2 Γ ( n ) Γ ( n 2 ) cos ( x 2 n n π 2 ) ( 1 x 2 2 n + 1 ) 1 4 , {\displaystyle e^{-{\frac {x^{2}}{2}}}\cdot H_{n}(x)\sim {\frac {2^{n}}{\sqrt {\pi }}}\Gamma \left({\frac {n+1}{2}}\right)\cos \left(x{\sqrt {2n}}-{\frac {n\pi }{2}}\right)\left(1-{\frac {x^{2}}{2n+1}}\right)^{-{\frac {1}{4}}}={\frac {2\Gamma (n)}{\Gamma \left({\frac {n}{2}}\right)}}\cos \left(x{\sqrt {2n}}-{\frac {n\pi }{2}}\right)\left(1-{\frac {x^{2}}{2n+1}}\right)^{-{\frac {1}{4}}},} e x 2 2 H n ( x ) ( 2 n e ) n 2 2 cos ( x 2 n n π 2 ) ( 1 x 2 2 n + 1 ) 1 4 . {\displaystyle e^{-{\frac {x^{2}}{2}}}\cdot H_{n}(x)\sim \left({\frac {2n}{e}}\right)^{\frac {n}{2}}{\sqrt {2}}\cos \left(x{\sqrt {2n}}-{\frac {n\pi }{2}}\right)\left(1-{\frac {x^{2}}{2n+1}}\right)^{-{\frac {1}{4}}}.}

Esta expansión es necesaria para resolver la función de onda de un oscilador armónico cuántico de manera que concuerde con la aproximación clásica en el límite del principio de correspondencia .

Una mejor aproximación, que tiene en cuenta la variación en la frecuencia, viene dada por e x 2 2 H n ( x ) ( 2 n e ) n 2 2 cos ( x 2 n + 1 x 2 3 n π 2 ) ( 1 x 2 2 n + 1 ) 1 4 . {\displaystyle e^{-{\frac {x^{2}}{2}}}\cdot H_{n}(x)\sim \left({\frac {2n}{e}}\right)^{\frac {n}{2}}{\sqrt {2}}\cos \left(x{\sqrt {2n+1-{\frac {x^{2}}{3}}}}-{\frac {n\pi }{2}}\right)\left(1-{\frac {x^{2}}{2n+1}}\right)^{-{\frac {1}{4}}}.}

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 x = 2 n + 1 cos ( φ ) , 0 < ε φ π ε , {\displaystyle x={\sqrt {2n+1}}\cos(\varphi ),\quad 0<\varepsilon \leq \varphi \leq \pi -\varepsilon ,} e x 2 2 H n ( x ) = 2 n 2 + 1 4 n ! ( π n ) 1 4 ( sin φ ) 1 2 ( sin ( 3 π 4 + ( n 2 + 1 4 ) ( sin 2 φ 2 φ ) ) + O ( n 1 ) ) . {\displaystyle e^{-{\frac {x^{2}}{2}}}\cdot H_{n}(x)=2^{{\frac {n}{2}}+{\frac {1}{4}}}{\sqrt {n!}}(\pi n)^{-{\frac {1}{4}}}(\sin \varphi )^{-{\frac {1}{2}}}\cdot \left(\sin \left({\frac {3\pi }{4}}+\left({\frac {n}{2}}+{\frac {1}{4}}\right)\left(\sin 2\varphi -2\varphi \right)\right)+O\left(n^{-1}\right)\right).}

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. x = 2 n + 1 cosh ( φ ) , 0 < ε φ ω < , {\displaystyle x={\sqrt {2n+1}}\cosh(\varphi ),\quad 0<\varepsilon \leq \varphi \leq \omega <\infty ,} e x 2 2 H n ( x ) = 2 n 2 3 4 n ! ( π n ) 1 4 ( sinh φ ) 1 2 e ( n 2 + 1 4 ) ( 2 φ sinh 2 φ ) ( 1 + O ( n 1 ) ) , {\displaystyle e^{-{\frac {x^{2}}{2}}}\cdot H_{n}(x)=2^{{\frac {n}{2}}-{\frac {3}{4}}}{\sqrt {n!}}(\pi n)^{-{\frac {1}{4}}}(\sinh \varphi )^{-{\frac {1}{2}}}\cdot e^{\left({\frac {n}{2}}+{\frac {1}{4}}\right)\left(2\varphi -\sinh 2\varphi \right)}\left(1+O\left(n^{-1}\right)\right),} x = 2 n + 1 + t {\displaystyle x={\sqrt {2n+1}}+t} e x 2 2 H n ( x ) = π 1 4 2 n 2 + 1 4 n ! n 1 12 ( Ai ( 2 1 2 n 1 6 t ) + O ( n 2 3 ) ) , {\displaystyle e^{-{\frac {x^{2}}{2}}}\cdot H_{n}(x)=\pi ^{\frac {1}{4}}2^{{\frac {n}{2}}+{\frac {1}{4}}}{\sqrt {n!}}\,n^{-{\frac {1}{12}}}\left(\operatorname {Ai} \left(2^{\frac {1}{2}}n^{\frac {1}{6}}t\right)+O\left(n^{-{\frac {2}{3}}}\right)\right),}

Valores especiales

Los polinomios de Hermite del físico evaluados con argumento cero H n (0) se denominan números de Hermite .

H n ( 0 ) = { 0 for odd  n , ( 2 ) n 2 ( n 1 ) ! ! for even  n , {\displaystyle H_{n}(0)={\begin{cases}0&{\text{for odd }}n,\\(-2)^{\frac {n}{2}}(n-1)!!&{\text{for even }}n,\end{cases}}} 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 He n ( 0 ) = { 0 for odd  n , ( 1 ) n 2 ( n 1 ) ! ! for even  n . {\displaystyle \operatorname {He} _{n}(0)={\begin{cases}0&{\text{for odd }}n,\\(-1)^{\frac {n}{2}}(n-1)!!&{\text{for even }}n.\end{cases}}}

Relaciones con otras funciones

Polinomios de Laguerre

Los polinomios de Hermite se pueden expresar como un caso especial de los polinomios de Laguerre : H 2 n ( x ) = ( 4 ) n n ! L n ( 1 2 ) ( x 2 ) = 4 n n ! k = 0 n ( 1 ) n k ( n 1 2 n k ) x 2 k k ! , H 2 n + 1 ( x ) = 2 ( 4 ) n n ! x L n ( 1 2 ) ( x 2 ) = 2 4 n n ! k = 0 n ( 1 ) n k ( n + 1 2 n k ) x 2 k + 1 k ! . {\displaystyle {\begin{aligned}H_{2n}(x)&=(-4)^{n}n!L_{n}^{\left(-{\frac {1}{2}}\right)}(x^{2})&&=4^{n}n!\sum _{k=0}^{n}(-1)^{n-k}{\binom {n-{\frac {1}{2}}}{n-k}}{\frac {x^{2k}}{k!}},\\H_{2n+1}(x)&=2(-4)^{n}n!xL_{n}^{\left({\frac {1}{2}}\right)}(x^{2})&&=2\cdot 4^{n}n!\sum _{k=0}^{n}(-1)^{n-k}{\binom {n+{\frac {1}{2}}}{n-k}}{\frac {x^{2k+1}}{k!}}.\end{aligned}}}

Relación con funciones hipergeométricas confluentes

Los polinomios de Hermite del físico se pueden expresar como un caso especial de las funciones de cilindro parabólico : en el semiplano derecho , donde U ( a , b , z ) es la función hipergeométrica confluente de Tricomi . De manera similar, donde 1 F 1 ( a , b ; z ) = M ( a , b ; z ) es la función hipergeométrica confluente de Kummer . H n ( x ) = 2 n U ( 1 2 n , 1 2 , x 2 ) {\displaystyle H_{n}(x)=2^{n}U\left(-{\tfrac {1}{2}}n,{\tfrac {1}{2}},x^{2}\right)} H 2 n ( x ) = ( 1 ) n ( 2 n ) ! n ! 1 F 1 ( n , 1 2 ; x 2 ) , H 2 n + 1 ( x ) = ( 1 ) n ( 2 n + 1 ) ! n ! 2 x 1 F 1 ( n , 3 2 ; x 2 ) , {\displaystyle {\begin{aligned}H_{2n}(x)&=(-1)^{n}{\frac {(2n)!}{n!}}\,_{1}F_{1}{\big (}-n,{\tfrac {1}{2}};x^{2}{\big )},\\H_{2n+1}(x)&=(-1)^{n}{\frac {(2n+1)!}{n!}}\,2x\,_{1}F_{1}{\big (}-n,{\tfrac {3}{2}};x^{2}{\big )},\end{aligned}}}

Expansión del polinomio de Hermite

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] e x 2 f ( x ) 2 d x < {\displaystyle \int e^{-x^{2}}f(x)^{2}dx<\infty }

Dado tal , las sumas parciales de la expansión de Hermite de convergen a en la norma si y sólo si . [11] f {\displaystyle f} f {\displaystyle f} L p {\displaystyle L^{p}} 4 / 3 < p < 4 {\displaystyle 4/3<p<4} x n = n ! 2 n k = 0 n / 2 1 k ! ( n 2 k ) ! H n 2 k ( x ) = n ! k = 0 n / 2 1 k ! 2 k ( n 2 k ) ! He n 2 k ( x ) , n Z + . {\displaystyle x^{n}={\frac {n!}{2^{n}}}\,\sum _{k=0}^{\left\lfloor n/2\right\rfloor }{\frac {1}{k!\,(n-2k)!}}\,H_{n-2k}(x)=n!\sum _{k=0}^{\left\lfloor n/2\right\rfloor }{\frac {1}{k!\,2^{k}\,(n-2k)!}}\,\operatorname {He} _{n-2k}(x),\qquad n\in \mathbb {Z} _{+}.} e a x = e a 2 / 4 n 0 a n n ! 2 n H n ( x ) , a C , x R . {\displaystyle e^{ax}=e^{a^{2}/4}\sum _{n\geq 0}{\frac {a^{n}}{n!\,2^{n}}}\,H_{n}(x),\qquad a\in \mathbb {C} ,\quad x\in \mathbb {R} .} e a 2 x 2 = n 0 ( 1 ) n a 2 n n ! ( 1 + a 2 ) n + 1 / 2 2 2 n H 2 n ( x ) . {\displaystyle e^{-a^{2}x^{2}}=\sum _{n\geq 0}{\frac {(-1)^{n}a^{2n}}{n!\left(1+a^{2}\right)^{n+1/2}2^{2n}}}\,H_{2n}(x).} erf ( x ) = 2 π 0 x e t 2   d t = 1 2 π k 0 ( 1 ) k k ! ( 2 k + 1 ) 2 3 k H 2 k ( x ) . {\displaystyle \operatorname {erf} (x)={\frac {2}{\sqrt {\pi }}}\int _{0}^{x}e^{-t^{2}}~dt={\frac {1}{\sqrt {2\pi }}}\sum _{k\geq 0}{\frac {(-1)^{k}}{k!(2k+1)2^{3k}}}H_{2k}(x).} cosh ( 2 x ) = e k 0 1 ( 2 k ) ! H 2 k ( x ) , sinh ( 2 x ) = e k 0 1 ( 2 k + 1 ) ! H 2 k + 1 ( x ) . {\displaystyle \cosh(2x)=e\sum _{k\geq 0}{\frac {1}{(2k)!}}\,H_{2k}(x),\qquad \sinh(2x)=e\sum _{k\geq 0}{\frac {1}{(2k+1)!}}\,H_{2k+1}(x).} cos ( x ) = e 1 / 4 k 0 ( 1 ) k 2 2 k ( 2 k ) ! H 2 k ( x ) sin ( x ) = e 1 / 4 k 0 ( 1 ) k 2 2 k + 1 ( 2 k + 1 ) ! H 2 k + 1 ( x ) {\displaystyle \cos(x)=e^{-1/4}\,\sum _{k\geq 0}{\frac {(-1)^{k}}{2^{2k}\,(2k)!}}\,H_{2k}(x)\quad \sin(x)=e^{-1/4}\,\sum _{k\geq 0}{\frac {(-1)^{k}}{2^{2k+1}\,(2k+1)!}}\,H_{2k+1}(x)}

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. He n ( x ) = e D 2 2 x n , {\displaystyle \operatorname {He} _{n}(x)=e^{-{\frac {D^{2}}{2}}}x^{n},}

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. He n ( x ) = n ! 2 π i C e t x t 2 2 t n + 1 d t , H n ( x ) = n ! 2 π i C e 2 t x t 2 t n + 1 d t , {\displaystyle {\begin{aligned}\operatorname {He} _{n}(x)&={\frac {n!}{2\pi i}}\oint _{C}{\frac {e^{tx-{\frac {t^{2}}{2}}}}{t^{n+1}}}\,dt,\\H_{n}(x)&={\frac {n!}{2\pi i}}\oint _{C}{\frac {e^{2tx-t^{2}}}{t^{n+1}}}\,dt,\end{aligned}}}

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. 1 2 π e x 2 2 , {\displaystyle {\frac {1}{\sqrt {2\pi }}}e^{-{\frac {x^{2}}{2}}},}

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 He n [ α ] ( x ) {\displaystyle \operatorname {He} _{n}^{[\alpha ]}(x)} ( 2 π α ) 1 2 e x 2 2 α . {\displaystyle (2\pi \alpha )^{-{\frac {1}{2}}}e^{-{\frac {x^{2}}{2\alpha }}}.} He n [ α ] ( x ) = α n 2 He n ( x α ) = ( α 2 ) n 2 H n ( x 2 α ) = e α D 2 2 ( x n ) . {\displaystyle \operatorname {He} _{n}^{[\alpha ]}(x)=\alpha ^{\frac {n}{2}}\operatorname {He} _{n}\left({\frac {x}{\sqrt {\alpha }}}\right)=\left({\frac {\alpha }{2}}\right)^{\frac {n}{2}}H_{n}\left({\frac {x}{\sqrt {2\alpha }}}\right)=e^{-{\frac {\alpha D^{2}}{2}}}\left(x^{n}\right).}

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 α = β = He n [ α ] ( x ) = k = 0 n h n , k [ α ] x k , {\displaystyle \operatorname {He} _{n}^{[\alpha ]}(x)=\sum _{k=0}^{n}h_{n,k}^{[\alpha ]}x^{k},} ( He n [ α ] He [ β ] ) ( x ) k = 0 n h n , k [ α ] He k [ β ] ( x ) {\displaystyle \left(\operatorname {He} _{n}^{[\alpha ]}\circ \operatorname {He} ^{[\beta ]}\right)(x)\equiv \sum _{k=0}^{n}h_{n,k}^{[\alpha ]}\,\operatorname {He} _{k}^{[\beta ]}(x)} ( He n [ α ] He [ β ] ) ( x ) = He n [ α + β ] ( x ) {\displaystyle \left(\operatorname {He} _{n}^{[\alpha ]}\circ \operatorname {He} ^{[\beta ]}\right)(x)=\operatorname {He} _{n}^{[\alpha +\beta ]}(x)} He n [ α + β ] ( x + y ) = k = 0 n ( n k ) He k [ α ] ( x ) He n k [ β ] ( y ) . {\displaystyle \operatorname {He} _{n}^{[\alpha +\beta ]}(x+y)=\sum _{k=0}^{n}{\binom {n}{k}}\operatorname {He} _{k}^{[\alpha ]}(x)\operatorname {He} _{n-k}^{[\beta ]}(y).} 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 . He n [ α ] ( x ) {\displaystyle \operatorname {He} _{n}^{[-\alpha ]}(x)} He n [ α ] ( x ) {\displaystyle \operatorname {He} _{n}^{[-\alpha ]}(x)} He n [ α ] ( x ) {\displaystyle \operatorname {He} _{n}^{[\alpha ]}(x)}

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 E [ X n ] = He n [ σ 2 ] ( μ ) , {\displaystyle E[X^{n}]=\operatorname {He} _{n}^{[-\sigma ^{2}]}(\mu ),} k = 0 n ( n k ) He k [ α ] ( x ) He n k [ α ] ( y ) = He n [ 0 ] ( x + y ) = ( x + y ) n . {\displaystyle \sum _{k=0}^{n}{\binom {n}{k}}\operatorname {He} _{k}^{[\alpha ]}(x)\operatorname {He} _{n-k}^{[-\alpha ]}(y)=\operatorname {He} _{n}^{[0]}(x+y)=(x+y)^{n}.}

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í, ψ n ( x ) = ( 2 n n ! π ) 1 2 e x 2 2 H n ( x ) = ( 1 ) n ( 2 n n ! π ) 1 2 e x 2 2 d n d x n e x 2 . {\displaystyle \psi _{n}(x)=\left(2^{n}n!{\sqrt {\pi }}\right)^{-{\frac {1}{2}}}e^{-{\frac {x^{2}}{2}}}H_{n}(x)=(-1)^{n}\left(2^{n}n!{\sqrt {\pi }}\right)^{-{\frac {1}{2}}}e^{\frac {x^{2}}{2}}{\frac {d^{n}}{dx^{n}}}e^{-x^{2}}.} 2 ( n + 1 )     ψ n + 1 ( x ) = ( x d d x ) ψ n ( x ) . {\displaystyle {\sqrt {2(n+1)}}~~\psi _{n+1}(x)=\left(x-{d \over dx}\right)\psi _{n}(x).}

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). ψ n ( x ) ψ m ( x ) d x = δ n m , {\displaystyle \int _{-\infty }^{\infty }\psi _{n}(x)\psi _{m}(x)\,dx=\delta _{nm},}

Las funciones de Hermite están estrechamente relacionadas con la función de Whittaker (Whittaker y Watson 1996) D n ( z ) : y por tanto con otras funciones de cilindros parabólicos . D n ( z ) = ( n ! π ) 1 2 ψ n ( z 2 ) = ( 1 ) n e z 2 4 d n d z n e z 2 2 {\displaystyle D_{n}(z)=\left(n!{\sqrt {\pi }}\right)^{\frac {1}{2}}\psi _{n}\left({\frac {z}{\sqrt {2}}}\right)=(-1)^{n}e^{\frac {z^{2}}{4}}{\frac {d^{n}}{dz^{n}}}e^{\frac {-z^{2}}{2}}}

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 . ψ n ( x ) + ( 2 n + 1 x 2 ) ψ n ( x ) = 0. {\displaystyle \psi _{n}''(x)+\left(2n+1-x^{2}\right)\psi _{n}(x)=0.}

Funciones de Hermite: 0 (azul, sólido), 1 (naranja, discontinuo), 2 (verde, discontinuo), 3 (rojo, discontinuo), 4 (violeta, sólido) y 5 (marrón, discontinuo)

ψ 0 ( x ) = π 1 4 e 1 2 x 2 , ψ 1 ( x ) = 2 π 1 4 x e 1 2 x 2 , ψ 2 ( x ) = ( 2 π 1 4 ) 1 ( 2 x 2 1 ) e 1 2 x 2 , ψ 3 ( x ) = ( 3 π 1 4 ) 1 ( 2 x 3 3 x ) e 1 2 x 2 , ψ 4 ( x ) = ( 2 6 π 1 4 ) 1 ( 4 x 4 12 x 2 + 3 ) e 1 2 x 2 , ψ 5 ( x ) = ( 2 15 π 1 4 ) 1 ( 4 x 5 20 x 3 + 15 x ) e 1 2 x 2 . {\displaystyle {\begin{aligned}\psi _{0}(x)&=\pi ^{-{\frac {1}{4}}}\,e^{-{\frac {1}{2}}x^{2}},\\\psi _{1}(x)&={\sqrt {2}}\,\pi ^{-{\frac {1}{4}}}\,x\,e^{-{\frac {1}{2}}x^{2}},\\\psi _{2}(x)&=\left({\sqrt {2}}\,\pi ^{\frac {1}{4}}\right)^{-1}\,\left(2x^{2}-1\right)\,e^{-{\frac {1}{2}}x^{2}},\\\psi _{3}(x)&=\left({\sqrt {3}}\,\pi ^{\frac {1}{4}}\right)^{-1}\,\left(2x^{3}-3x\right)\,e^{-{\frac {1}{2}}x^{2}},\\\psi _{4}(x)&=\left(2{\sqrt {6}}\,\pi ^{\frac {1}{4}}\right)^{-1}\,\left(4x^{4}-12x^{2}+3\right)\,e^{-{\frac {1}{2}}x^{2}},\\\psi _{5}(x)&=\left(2{\sqrt {15}}\,\pi ^{\frac {1}{4}}\right)^{-1}\,\left(4x^{5}-20x^{3}+15x\right)\,e^{-{\frac {1}{2}}x^{2}}.\end{aligned}}}

Funciones de Hermite: 0 (azul, sólido), 2 (naranja, discontinuo), 4 (verde, discontinuo) y 50 (rojo, sólido)

Relación de recursión

Siguiendo las relaciones de recursión de los polinomios de Hermite, las funciones de Hermite obedecen y ψ n ( x ) = n 2 ψ n 1 ( x ) n + 1 2 ψ n + 1 ( x ) {\displaystyle \psi _{n}'(x)={\sqrt {\frac {n}{2}}}\,\psi _{n-1}(x)-{\sqrt {\frac {n+1}{2}}}\psi _{n+1}(x)} x ψ n ( x ) = n 2 ψ n 1 ( x ) + n + 1 2 ψ n + 1 ( x ) . {\displaystyle x\psi _{n}(x)={\sqrt {\frac {n}{2}}}\,\psi _{n-1}(x)+{\sqrt {\frac {n+1}{2}}}\psi _{n+1}(x).}

Extender la primera relación a las derivadas m arbitrarias para cualquier entero positivo m conduce a ψ n ( m ) ( x ) = k = 0 m ( m k ) ( 1 ) k 2 m k 2 n ! ( n m + k ) ! ψ n m + k ( x ) He k ( x ) . {\displaystyle \psi _{n}^{(m)}(x)=\sum _{k=0}^{m}{\binom {m}{k}}(-1)^{k}2^{\frac {m-k}{2}}{\sqrt {\frac {n!}{(n-m+k)!}}}\psi _{n-m+k}(x)\operatorname {He} _{k}(x).}

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] | ψ n ( x ) | π 1 4 . {\displaystyle {\bigl |}\psi _{n}(x){\bigr |}\leq \pi ^{-{\frac {1}{4}}}.}

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/2x 2 . Esto da e 1 2 x 2 + 2 x t t 2 = n = 0 e 1 2 x 2 H n ( x ) t n n ! . {\displaystyle e^{-{\frac {1}{2}}x^{2}+2xt-t^{2}}=\sum _{n=0}^{\infty }e^{-{\frac {1}{2}}x^{2}}H_{n}(x){\frac {t^{n}}{n!}}.}

La transformada de Fourier del lado izquierdo está dada por F { e 1 2 x 2 + 2 x t t 2 } ( k ) = 1 2 π e i x k e 1 2 x 2 + 2 x t t 2 d x = e 1 2 k 2 2 k i t + t 2 = n = 0 e 1 2 k 2 H n ( k ) ( i t ) n n ! . {\displaystyle {\begin{aligned}{\mathcal {F}}\left\{e^{-{\frac {1}{2}}x^{2}+2xt-t^{2}}\right\}(k)&={\frac {1}{\sqrt {2\pi }}}\int _{-\infty }^{\infty }e^{-ixk}e^{-{\frac {1}{2}}x^{2}+2xt-t^{2}}\,dx\\&=e^{-{\frac {1}{2}}k^{2}-2kit+t^{2}}\\&=\sum _{n=0}^{\infty }e^{-{\frac {1}{2}}k^{2}}H_{n}(k){\frac {(-it)^{n}}{n!}}.\end{aligned}}}

La transformada de Fourier del lado derecho está dada por F { n = 0 e 1 2 x 2 H n ( x ) t n n ! } = n = 0 F { e 1 2 x 2 H n ( x ) } t n n ! . {\displaystyle {\mathcal {F}}\left\{\sum _{n=0}^{\infty }e^{-{\frac {1}{2}}x^{2}}H_{n}(x){\frac {t^{n}}{n!}}\right\}=\sum _{n=0}^{\infty }{\mathcal {F}}\left\{e^{-{\frac {1}{2}}x^{2}}H_{n}(x)\right\}{\frac {t^{n}}{n!}}.}

Igualando potencias iguales de t en las versiones transformadas de los lados izquierdo y derecho finalmente obtenemos F { e 1 2 x 2 H n ( x ) } = ( i ) n e 1 2 k 2 H n ( k ) . {\displaystyle {\mathcal {F}}\left\{e^{-{\frac {1}{2}}x^{2}}H_{n}(x)\right\}=(-i)^{n}e^{-{\frac {1}{2}}k^{2}}H_{n}(k).}

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 xL 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 . L n ( x ) := k = 0 n ( n k ) ( 1 ) k k ! x k , {\displaystyle L_{n}(x):=\sum _{k=0}^{n}{\binom {n}{k}}{\frac {(-1)^{k}}{k!}}x^{k},} l n ( x ) := e x 2 L n ( x ) . {\displaystyle l_{n}(x):=e^{-{\frac {x}{2}}}L_{n}(x).} W ψ n ( t , f ) = ( 1 ) n l n ( 4 π ( t 2 + f 2 ) ) , {\displaystyle W_{\psi _{n}}(t,f)=(-1)^{n}l_{n}{\big (}4\pi (t^{2}+f^{2}){\big )},} W x ( t , f ) = x ( t + τ 2 ) x ( t τ 2 ) e 2 π i τ f d τ . {\displaystyle W_{x}(t,f)=\int _{-\infty }^{\infty }x\left(t+{\frac {\tau }{2}}\right)\,x\left(t-{\frac {\tau }{2}}\right)^{*}\,e^{-2\pi i\tau f}\,d\tau .}

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 . He n ( x ) = B n ( x , 1 , 0 , , 0 ) , {\displaystyle \operatorname {He} _{n}(x)=B_{n}(x,-1,0,\ldots ,0),}

Estos números también pueden expresarse como un valor especial de los polinomios de Hermite: [19] T ( n ) = He n ( i ) i n . {\displaystyle T(n)={\frac {\operatorname {He} _{n}(i)}{i^{n}}}.}

Relación de completitud

La fórmula de Christoffel-Darboux para polinomios de Hermite se lee k = 0 n H k ( x ) H k ( y ) k ! 2 k = 1 n ! 2 n + 1 H n ( y ) H n + 1 ( x ) H n ( x ) H n + 1 ( y ) x y . {\displaystyle \sum _{k=0}^{n}{\frac {H_{k}(x)H_{k}(y)}{k!2^{k}}}={\frac {1}{n!2^{n+1}}}\,{\frac {H_{n}(y)H_{n+1}(x)-H_{n}(x)H_{n+1}(y)}{x-y}}.}

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 δ ( xy ) 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. n = 0 ψ n ( x ) ψ n ( y ) = δ ( x y ) , {\displaystyle \sum _{n=0}^{\infty }\psi _{n}(x)\psi _{n}(y)=\delta (x-y),}

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] E ( x , y ; u ) := n = 0 u n ψ n ( x ) ψ n ( y ) = 1 π ( 1 u 2 ) exp ( 1 u 1 + u ( x + y ) 2 4 1 + u 1 u ( x y ) 2 4 ) , {\displaystyle E(x,y;u):=\sum _{n=0}^{\infty }u^{n}\,\psi _{n}(x)\,\psi _{n}(y)={\frac {1}{\sqrt {\pi (1-u^{2})}}}\,\exp \left(-{\frac {1-u}{1+u}}\,{\frac {(x+y)^{2}}{4}}-{\frac {1+u}{1-u}}\,{\frac {(x-y)^{2}}{4}}\right),} n = 0 H n ( x ) H n ( y ) n ! ( u 2 ) n = 1 1 u 2 e 2 u 1 + u x y u 2 1 u 2 ( x y ) 2 . {\displaystyle \sum _{n=0}^{\infty }{\frac {H_{n}(x)H_{n}(y)}{n!}}\left({\frac {u}{2}}\right)^{n}={\frac {1}{\sqrt {1-u^{2}}}}e^{{\frac {2u}{1+u}}xy-{\frac {u^{2}}{1-u^{2}}}(x-y)^{2}}.}

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. n = 0 u n f , ψ n ψ n , g = E ( x , y ; u ) f ( x ) g ( y ) ¯ d x d y f ( x ) g ( x ) ¯ d x = f , g {\displaystyle \sum _{n=0}^{\infty }u^{n}\langle f,\psi _{n}\rangle \langle \psi _{n},g\rangle =\iint E(x,y;u)f(x){\overline {g(y)}}\,dx\,dy\to \int f(x){\overline {g(x)}}\,dx=\langle f,g\rangle }

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, f = n = 0 f , ψ n ψ n . {\displaystyle f=\sum _{n=0}^{\infty }\langle f,\psi _{n}\rangle \psi _{n}.}

Para demostrar la igualdad anterior para E ( x , y ; u ) , se utiliza repetidamente la transformada de Fourier de funciones gaussianas : ρ π e ρ 2 x 2 4 = e i s x s 2 ρ 2 d s for  ρ > 0. {\displaystyle \rho {\sqrt {\pi }}e^{-{\frac {\rho ^{2}x^{2}}{4}}}=\int e^{isx-{\frac {s^{2}}{\rho ^{2}}}}\,ds\quad {\text{for }}\rho >0.}

El polinomio de Hermite se representa entonces como H n ( x ) = ( 1 ) n e x 2 d n d x n ( 1 2 π e i s x s 2 4 d s ) = ( 1 ) n e x 2 1 2 π ( i s ) n e i s x s 2 4 d s . {\displaystyle H_{n}(x)=(-1)^{n}e^{x^{2}}{\frac {d^{n}}{dx^{n}}}\left({\frac {1}{2{\sqrt {\pi }}}}\int e^{isx-{\frac {s^{2}}{4}}}\,ds\right)=(-1)^{n}e^{x^{2}}{\frac {1}{2{\sqrt {\pi }}}}\int (is)^{n}e^{isx-{\frac {s^{2}}{4}}}\,ds.}

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 E ( x , y ; u ) = n = 0 u n 2 n n ! π H n ( x ) H n ( y ) e x 2 + y 2 2 = e x 2 + y 2 2 4 π π ( n = 0 1 2 n n ! ( u s t ) n ) e i s x + i t y s 2 4 t 2 4 d s d t = e x 2 + y 2 2 4 π π e u s t 2 e i s x + i t y s 2 4 t 2 4 d s d t , {\displaystyle {\begin{aligned}E(x,y;u)&=\sum _{n=0}^{\infty }{\frac {u^{n}}{2^{n}n!{\sqrt {\pi }}}}\,H_{n}(x)H_{n}(y)e^{-{\frac {x^{2}+y^{2}}{2}}}\\&={\frac {e^{\frac {x^{2}+y^{2}}{2}}}{4\pi {\sqrt {\pi }}}}\iint \left(\sum _{n=0}^{\infty }{\frac {1}{2^{n}n!}}(-ust)^{n}\right)e^{isx+ity-{\frac {s^{2}}{4}}-{\frac {t^{2}}{4}}}\,ds\,dt\\&={\frac {e^{\frac {x^{2}+y^{2}}{2}}}{4\pi {\sqrt {\pi }}}}\iint e^{-{\frac {ust}{2}}}\,e^{isx+ity-{\frac {s^{2}}{4}}-{\frac {t^{2}}{4}}}\,ds\,dt,\end{aligned}}} s = σ + τ 2 , t = σ τ 2 . {\displaystyle s={\frac {\sigma +\tau }{\sqrt {2}}},\quad t={\frac {\sigma -\tau }{\sqrt {2}}}.}

Véase también

Notas

  1. ^ 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.
  2. ^ 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.
  3. ^ 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.
  4. ^ 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.
  5. ^ Tom H. Koornwinder, Roderick SC Wong y Roelof Koekoek et al. (2010) y Abramowitz y Stegun .
  6. ^ 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.
  7. ^ "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 .
  8. ^ Abramowitz y Stegun 1983, pág. 508–510, 13.6.38 y 13.5.16.
  9. ^ Szegő 1955, pág. 201
  10. ^ "Tutorial de MATHEMATICA, parte 2.5: Expansión de Hermite". www.cfm.brown.edu . Consultado el 24 de diciembre de 2023 .
  11. ^ 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.
  12. ^ Roman, Steven (1984), El cálculo umbral , Matemáticas puras y aplicadas, vol. 111 (1.ª ed.), Academic Press, págs. 87-93, ISBN 978-0-12-594380-2
  13. ^ Erdélyi y col. 1955, pág. 207.
  14. ^ Szegő 1955.
  15. ^ 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
  16. ^ 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 .
  17. ^ Folland, GB (1989), Análisis armónico en el espacio de fases , Anales de estudios matemáticos, vol. 122, Princeton University Press, ISBN 978-0-691-08528-9
  18. ^ 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.
  19. ^ 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
  20. ^ 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).
  21. ^ Erdélyi y col. 1955, pág. 194, 10.13 (22).

Referencias

Retrieved from "https://en.wikipedia.org/w/index.php?title=Hermite_polynomials&oldid=1238599588"