Polinomios de Jacobi





En matemáticas, los polinomios de Jacobi (ocasionalmente llamados polinomios hipergeométricos) P(α, β)
n
(x)
son una clase de polinomios ortogonales clásicos. Son ortogonales con respecto al peso
(1 − x)α(1 + x)β en el intervalo [−1, 1]. Los polinomios de Gegenbauer, y por lo tanto también los de Legendre, de Zernike y de Chebyshov, son casos especiales de los polinomios de Jacobi.[1]


Los polinomios de Jacobi fueron introducidos por el matemático alemán Carl Gustav Jakob Jacobi (1804-1851).




Índice






  • 1 Definiciones


    • 1.1 A través de la función hipergeométrica


    • 1.2 La fórmula de Rodrigues


    • 1.3 Expresión alternativa para el argumento real




  • 2 Propiedades básicas


    • 2.1 Ortogonalidad


    • 2.2 Relación de simetría


    • 2.3 Derivadas


    • 2.4 Ecuación diferencial


    • 2.5 Relaciones de recurrencia


    • 2.6 Función de generación




  • 3 Polinomios de Jacobi asintóticos


  • 4 Aplicaciones


    • 4.1 Matriz D de Wigner




  • 5 Véase también


  • 6 Referencias


  • 7 Lecturas relacionadas


  • 8 Enlaces externos





Definiciones



A través de la función hipergeométrica


Los polinomios de Jacobi se definen a través de la función hipergeométrica de la siguiente manera:[2]


Pn(α)(z)=(α+1)nn!2F1(−n,1+α+n;α+1;12(1−z)),{displaystyle P_{n}^{(alpha ,beta )}(z)={frac {(alpha +1)_{n}}{n!}},{}_{2}F_{1}left(-n,1+alpha +beta +n;alpha +1;{tfrac {1}{2}}(1-z)right),}{displaystyle P_{n}^{(alpha ,beta )}(z)={frac {(alpha +1)_{n}}{n!}},{}_{2}F_{1}left(-n,1+alpha +beta +n;alpha +1;{tfrac {1}{2}}(1-z)right),}

donde +1)n{displaystyle (alpha +1)_{n}}{displaystyle (alpha +1)_{n}} es un símbolo de Pochhammer (para el factorial ascendente). En este caso, la serie para la función hipergeométrica es finita, por lo tanto, se obtiene la siguiente expresión equivalente:


Pn(α)(z)=Γ+n+1)n!Γ+n+1)∑m=0n(nm)Γ+n+m+1)Γ+m+1)(z−12)m.{displaystyle P_{n}^{(alpha ,beta )}(z)={frac {Gamma (alpha +n+1)}{n!,Gamma (alpha +beta +n+1)}}sum _{m=0}^{n}{n choose m}{frac {Gamma (alpha +beta +n+m+1)}{Gamma (alpha +m+1)}}left({frac {z-1}{2}}right)^{m}.}{displaystyle P_{n}^{(alpha ,beta )}(z)={frac {Gamma (alpha +n+1)}{n!,Gamma (alpha +beta +n+1)}}sum _{m=0}^{n}{n choose m}{frac {Gamma (alpha +beta +n+m+1)}{Gamma (alpha +m+1)}}left({frac {z-1}{2}}right)^{m}.}


La fórmula de Rodrigues


La fórmula de Rodrigues da una definición equivalente:[1][3]


Pn(α)(z)=(−1)n2nn!(1−z)−α(1+z)−βdndzn{(1−z)α(1+z)β(1−z2)n}.{displaystyle P_{n}^{(alpha ,beta )}(z)={frac {(-1)^{n}}{2^{n}n!}}(1-z)^{-alpha }(1+z)^{-beta }{frac {d^{n}}{dz^{n}}}left{(1-z)^{alpha }(1+z)^{beta }left(1-z^{2}right)^{n}right}.}{displaystyle P_{n}^{(alpha ,beta )}(z)={frac {(-1)^{n}}{2^{n}n!}}(1-z)^{-alpha }(1+z)^{-beta }{frac {d^{n}}{dz^{n}}}left{(1-z)^{alpha }(1+z)^{beta }left(1-z^{2}right)^{n}right}.}

Si α=0{displaystyle alpha =beta =0}{displaystyle alpha =beta =0}, entonces se reduce a los polinomios de Legendre:


Pn(z)=12nn!dndzn(z2−1)n.{displaystyle P_{n}(z)={frac {1}{2^{n}n!}}{frac {d^{n}}{dz^{n}}}(z^{2}-1)^{n};.}{displaystyle P_{n}(z)={frac {1}{2^{n}n!}}{frac {d^{n}}{dz^{n}}}(z^{2}-1)^{n};.}


Expresión alternativa para el argumento real


Para x real, el polinomio de Jacobi puede escribirse alternativamente como


Pn(α)(x)=∑s=0n(n+αn−s)(n+βs)(x−12)s(x+12)n−s.{displaystyle P_{n}^{(alpha ,beta )}(x)=sum _{s=0}^{n}{n+alpha choose n-s}{n+beta choose s}left({frac {x-1}{2}}right)^{s}left({frac {x+1}{2}}right)^{n-s}.}{displaystyle P_{n}^{(alpha ,beta )}(x)=sum _{s=0}^{n}{n+alpha  choose n-s}{n+beta  choose s}left({frac {x-1}{2}}right)^{s}left({frac {x+1}{2}}right)^{n-s}.}

y para un número entero n


(zn)={Γ(z+1)Γ(n+1)Γ(z−n+1)n≥00n<0{displaystyle {z choose n}={begin{cases}{frac {Gamma (z+1)}{Gamma (n+1)Gamma (z-n+1)}}&ngeq 0\0&n<0end{cases}}}{displaystyle {z choose n}={begin{cases}{frac {Gamma (z+1)}{Gamma (n+1)Gamma (z-n+1)}}&ngeq 0\0&n<0end{cases}}}

donde Γ(z) es la Función gamma.


En el caso especial de que las cuatro cantidades n, n + α, n + β y n + α + β son enteros no negativos, el polinomio de Jacobi se puede escribir como









Pn(α)(x)=(n+α)!(n+β)!∑s1s!(n+αs)!(β+s)!(n−s)!(x−12)n−s(x+12)s.{displaystyle P_{n}^{(alpha ,beta )}(x)=(n+alpha )!(n+beta )!sum _{s}{frac {1}{s!(n+alpha -s)!(beta +s)!(n-s)!}}left({frac {x-1}{2}}right)^{n-s}left({frac {x+1}{2}}right)^{s}.}{displaystyle P_{n}^{(alpha ,beta )}(x)=(n+alpha )!(n+beta )!sum _{s}{frac {1}{s!(n+alpha -s)!(beta +s)!(n-s)!}}left({frac {x-1}{2}}right)^{n-s}left({frac {x+1}{2}}right)^{s}.}













 



 



 



 





(1)



La suma se extiende sobre todos los valores enteros de s para los cuales los argumentos de los factoriales no son negativos.



Propiedades básicas



Ortogonalidad


Los polinomios de Jacobi satisfacen la condición de ortogonalidad


11(1−x)α(1+x)βPm(α)(x)Pn(α)(x)dx=2α+12n+α+1Γ(n+α+1)Γ(n+β+1)Γ(n+α+1)n!δnm,α, βMAYORQUE−1.{displaystyle int _{-1}^{1}(1-x)^{alpha }(1+x)^{beta }P_{m}^{(alpha ,beta )}(x)P_{n}^{(alpha ,beta )}(x),dx={frac {2^{alpha +beta +1}}{2n+alpha +beta +1}}{frac {Gamma (n+alpha +1)Gamma (n+beta +1)}{Gamma (n+alpha +beta +1)n!}}delta _{nm},qquad alpha , beta MAYORQUE-1.}{displaystyle int _{-1}^{1}(1-x)^{alpha }(1+x)^{beta }P_{m}^{(alpha ,beta )}(x)P_{n}^{(alpha ,beta )}(x),dx={frac {2^{alpha +beta +1}}{2n+alpha +beta +1}}{frac {Gamma (n+alpha +1)Gamma (n+beta +1)}{Gamma (n+alpha +beta +1)n!}}delta _{nm},qquad alpha , beta MAYORQUE-1.}

Como se define, no tienen una norma unitaria con respecto al peso. Esto se puede corregir dividiendo por la raíz cuadrada del lado derecho de la ecuación anterior, cuando n=m{displaystyle n=m}{displaystyle n=m}.


Aunque no proporciona una base ortonormal, a veces se prefiere una normalización alternativa debido a su simplicidad:


Pn(α)(1)=(n+αn).{displaystyle P_{n}^{(alpha ,beta )}(1)={n+alpha choose n}.}{displaystyle P_{n}^{(alpha ,beta )}(1)={n+alpha  choose n}.}


Relación de simetría


Los polinomios tienen la relación de simetría


Pn(α)(−z)=(−1)nPn(β)(z);{displaystyle P_{n}^{(alpha ,beta )}(-z)=(-1)^{n}P_{n}^{(beta ,alpha )}(z);}{displaystyle P_{n}^{(alpha ,beta )}(-z)=(-1)^{n}P_{n}^{(beta ,alpha )}(z);}

por lo tanto, el otro valor terminal es


Pn(α)(−1)=(−1)n(n+βn).{displaystyle P_{n}^{(alpha ,beta )}(-1)=(-1)^{n}{n+beta choose n}.}{displaystyle P_{n}^{(alpha ,beta )}(-1)=(-1)^{n}{n+beta  choose n}.}


Derivadas


La k-ésima derivada de la expresión explícita conduce a


dkdzkPn(α)(z)=Γ+n+1+k)2kΓ+n+1)Pn−k(α+k,β+k)(z).{displaystyle {frac {d^{k}}{dz^{k}}}P_{n}^{(alpha ,beta )}(z)={frac {Gamma (alpha +beta +n+1+k)}{2^{k}Gamma (alpha +beta +n+1)}}P_{n-k}^{(alpha +k,beta +k)}(z).}{displaystyle {frac {d^{k}}{dz^{k}}}P_{n}^{(alpha ,beta )}(z)={frac {Gamma (alpha +beta +n+1+k)}{2^{k}Gamma (alpha +beta +n+1)}}P_{n-k}^{(alpha +k,beta +k)}(z).}


Ecuación diferencial


El polinomio de Jacobi P(α, β)
n
es una solución de una ecuación diferencial lineal de segundo orden[1]


(1−x2)y″+(βα+2)x)y′+n(n+α+1)y=0.{displaystyle left(1-x^{2}right)y''+(beta -alpha -(alpha +beta +2)x)y'+n(n+alpha +beta +1)y=0.}{displaystyle left(1-x^{2}right)y''+(beta -alpha -(alpha +beta +2)x)y'+n(n+alpha +beta +1)y=0.}


Relaciones de recurrencia


La relación de recurrencia para los polinomios de Jacobi de α, β fija es:[1]


2n(n+α)(2n+α2)Pn(α)(z)=(2n+α1){(2n+α)(2n+α2)z+α2−β2}Pn−1(α)(z)−2(n+α1)(n+β1)(2n+α)Pn−2(α)(z),{displaystyle {begin{aligned}&2n(n+alpha +beta )(2n+alpha +beta -2)P_{n}^{(alpha ,beta )}(z)\&qquad =(2n+alpha +beta -1){Big {}(2n+alpha +beta )(2n+alpha +beta -2)z+alpha ^{2}-beta ^{2}{Big }}P_{n-1}^{(alpha ,beta )}(z)-2(n+alpha -1)(n+beta -1)(2n+alpha +beta )P_{n-2}^{(alpha ,beta )}(z),end{aligned}}}{displaystyle {begin{aligned}&2n(n+alpha +beta )(2n+alpha +beta -2)P_{n}^{(alpha ,beta )}(z)\&qquad =(2n+alpha +beta -1){Big {}(2n+alpha +beta )(2n+alpha +beta -2)z+alpha ^{2}-beta ^{2}{Big }}P_{n-1}^{(alpha ,beta )}(z)-2(n+alpha -1)(n+beta -1)(2n+alpha +beta )P_{n-2}^{(alpha ,beta )}(z),end{aligned}}}

para n = 2, 3, ....


Dado que los polinomios de Jacobi se pueden describir en términos de la función hipergeométrica, las recurrencias de la función hipergeométrica dan recidivas equivalentes de los polinomios de Jacobi. En particular, las relaciones contiguas de Gauss corresponden a las identidades


(z−1)ddzPn(α)(z)=12(z−1)(1+α+n)Pn−1(α+1,β+1)=nPn(α)−+n)Pn−1(α+1)=(1+α+n)(Pn(α+1)−Pn(α))=(α+n)Pn(α1,β+1)−αPn(α)=2(n+1)Pn+1(α1)−(z(1+α+n)+α+1+n−β)Pn(α)1+z=(2β+n+nz)Pn(α)−2(β+n)Pn(α1)1+z=1−z1+z(βPn(α)−+n)Pn(α+1,β1)).{displaystyle {begin{aligned}(z-1){frac {d}{dz}}P_{n}^{(alpha ,beta )}(z)&={frac {1}{2}}(z-1)(1+alpha +beta +n)P_{n-1}^{(alpha +1,beta +1)}\&=nP_{n}^{(alpha ,beta )}-(alpha +n)P_{n-1}^{(alpha ,beta +1)}\&=(1+alpha +beta +n)left(P_{n}^{(alpha ,beta +1)}-P_{n}^{(alpha ,beta )}right)\&=(alpha +n)P_{n}^{(alpha -1,beta +1)}-alpha P_{n}^{(alpha ,beta )}\&={frac {2(n+1)P_{n+1}^{(alpha ,beta -1)}-left(z(1+alpha +beta +n)+alpha +1+n-beta right)P_{n}^{(alpha ,beta )}}{1+z}}\&={frac {(2beta +n+nz)P_{n}^{(alpha ,beta )}-2(beta +n)P_{n}^{(alpha ,beta -1)}}{1+z}}\&={frac {1-z}{1+z}}left(beta P_{n}^{(alpha ,beta )}-(beta +n)P_{n}^{(alpha +1,beta -1)}right),.end{aligned}}}{displaystyle {begin{aligned}(z-1){frac {d}{dz}}P_{n}^{(alpha ,beta )}(z)&={frac {1}{2}}(z-1)(1+alpha +beta +n)P_{n-1}^{(alpha +1,beta +1)}\&=nP_{n}^{(alpha ,beta )}-(alpha +n)P_{n-1}^{(alpha ,beta +1)}\&=(1+alpha +beta +n)left(P_{n}^{(alpha ,beta +1)}-P_{n}^{(alpha ,beta )}right)\&=(alpha +n)P_{n}^{(alpha -1,beta +1)}-alpha P_{n}^{(alpha ,beta )}\&={frac {2(n+1)P_{n+1}^{(alpha ,beta -1)}-left(z(1+alpha +beta +n)+alpha +1+n-beta right)P_{n}^{(alpha ,beta )}}{1+z}}\&={frac {(2beta +n+nz)P_{n}^{(alpha ,beta )}-2(beta +n)P_{n}^{(alpha ,beta -1)}}{1+z}}\&={frac {1-z}{1+z}}left(beta P_{n}^{(alpha ,beta )}-(beta +n)P_{n}^{(alpha +1,beta -1)}right),.end{aligned}}}


Función de generación


La función generadora de los polinomios de Jacobi está dada por


n=0∞Pn(α)(z)tn=2αR−1(1−t+R)−α(1+t+R)−β,{displaystyle sum _{n=0}^{infty }P_{n}^{(alpha ,beta )}(z)t^{n}=2^{alpha +beta }R^{-1}(1-t+R)^{-alpha }(1+t+R)^{-beta },}{displaystyle sum _{n=0}^{infty }P_{n}^{(alpha ,beta )}(z)t^{n}=2^{alpha +beta }R^{-1}(1-t+R)^{-alpha }(1+t+R)^{-beta },}

de donde


R=R(z,t)=(1−2zt+t2)12 ,{displaystyle R=R(z,t)=left(1-2zt+t^{2}right)^{frac {1}{2}}~,}{displaystyle R=R(z,t)=left(1-2zt+t^{2}right)^{frac {1}{2}}~,}

y la rama de la raíz cuadrada se elige para que R (z, 0) = 1.[1]



Polinomios de Jacobi asintóticos


Para x en el interior de [−1, 1], el término asintótico de P(α, β)
n
para n grande viene dado por la fórmula de Darboux[1]


Pn(α)(cos⁡θ)=n−12k(θ)cos⁡(Nθ)+O(n−32),{displaystyle P_{n}^{(alpha ,beta )}(cos theta )=n^{-{frac {1}{2}}}k(theta )cos(Ntheta +gamma )+Oleft(n^{-{frac {3}{2}}}right),}{displaystyle P_{n}^{(alpha ,beta )}(cos theta )=n^{-{frac {1}{2}}}k(theta )cos(Ntheta +gamma )+Oleft(n^{-{frac {3}{2}}}right),}

donde


k(θ)=π12sin−α12⁡θ2cos−β12⁡θ2,N=n+12(α+1),γ=−π2(α+12),{displaystyle {begin{aligned}k(theta )&=pi ^{-{frac {1}{2}}}sin ^{-alpha -{frac {1}{2}}}{tfrac {theta }{2}}cos ^{-beta -{frac {1}{2}}}{tfrac {theta }{2}},\N&=n+{tfrac {1}{2}}(alpha +beta +1),\gamma &=-{tfrac {pi }{2}}left(alpha +{tfrac {1}{2}}right),end{aligned}}}{displaystyle {begin{aligned}k(theta )&=pi ^{-{frac {1}{2}}}sin ^{-alpha -{frac {1}{2}}}{tfrac {theta }{2}}cos ^{-beta -{frac {1}{2}}}{tfrac {theta }{2}},\N&=n+{tfrac {1}{2}}(alpha +beta +1),\gamma &=-{tfrac {pi }{2}}left(alpha +{tfrac {1}{2}}right),end{aligned}}}

y el término "O" es uniforme en el intervalo [ε, π-ε] para cada ε > 0.


Los polinomios asintóticos de Jacobi cerca de los puntos ±1 vienen dados por la formula de Mehler-Heine


limn→n−αPn(α)(cos⁡(zn))=(z2)−α(z)limn→n−βPn(α)(cos⁡zn))=(z2)−β(z){displaystyle {begin{aligned}lim _{nto infty }n^{-alpha }P_{n}^{(alpha ,beta )}left(cos left({tfrac {z}{n}}right)right)&=left({tfrac {z}{2}}right)^{-alpha }J_{alpha }(z)\lim _{nto infty }n^{-beta }P_{n}^{(alpha ,beta )}left(cos left(pi -{tfrac {z}{n}}right)right)&=left({tfrac {z}{2}}right)^{-beta }J_{beta }(z)end{aligned}}}{displaystyle {begin{aligned}lim _{nto infty }n^{-alpha }P_{n}^{(alpha ,beta )}left(cos left({tfrac {z}{n}}right)right)&=left({tfrac {z}{2}}right)^{-alpha }J_{alpha }(z)\lim _{nto infty }n^{-beta }P_{n}^{(alpha ,beta )}left(cos left(pi -{tfrac {z}{n}}right)right)&=left({tfrac {z}{2}}right)^{-beta }J_{beta }(z)end{aligned}}}

donde los límites son uniformes para z en un dominio delimitado.


Los polinomios asintóticos fuera de [−1, 1] son menos explícitos.



Aplicaciones



Matriz D de Wigner


La expresión (1) permite la expresión de la matriz D de Wigner djm’,m(φ) (para 0 ≤ φ ≤ 4Plantilla:Pi) en términos de polinomios de Jacobi:[4]


dm′mj(ϕ)=[(j+m)!(j−m)!(j+m′)!(j−m′)!]12(sin⁡ϕ2)m−m′(cos⁡ϕ2)m+m′Pj−m(m−m′,m+m′)(cos⁡ϕ).{displaystyle d_{m'm}^{j}(phi )=left[{frac {(j+m)!(j-m)!}{(j+m')!(j-m')!}}right]^{frac {1}{2}}left(sin {tfrac {phi }{2}}right)^{m-m'}left(cos {tfrac {phi }{2}}right)^{m+m'}P_{j-m}^{(m-m',m+m')}(cos phi ).}{displaystyle d_{m'm}^{j}(phi )=left[{frac {(j+m)!(j-m)!}{(j+m')!(j-m')!}}right]^{frac {1}{2}}left(sin {tfrac {phi }{2}}right)^{m-m'}left(cos {tfrac {phi }{2}}right)^{m+m'}P_{j-m}^{(m-m',m+m')}(cos phi ).}


Véase también



  • Inecuación de Askey-Gasper

  • Grandes polinomios q de Jacobi

  • Polinomios continuos q de Jacobi

  • Pequeños polinomios q de Jacobi

  • Pseudo polinomios de Jacobi

  • Proceso de Jacobi

  • Polinomios de Gegenbauer

  • Polinomios de Romanovski



Referencias




  1. abcdef Szegő, Gábor (1939). «IV. Jacobi polynomials.». Orthogonal Polynomials. Colloquium Publications. XXIII. American Mathematical Society. ISBN 978-0-8218-1023-1. MR 0372517.  La definición figura en IV.1; la ecuación diferencial en IV.2; la fórmula de Rodrigues formula aparece en IV.3; la función generadora IV.4; y la relación recurrente está en IV.5.


  2. Abramowitz, Milton; Stegun, Irene Ann, eds. (1983) [June 1964]. "Chapter 22". Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables. Applied Mathematics Series. 55 (Ninth reprint with additional corrections of tenth original printing with corrections (December 1972); first ed.). Washington D.C.; New York: United States Department of Commerce, National Bureau of Standards; Dover Publications. p. 561. ISBN 978-0-486-61272-0. LCCN 64-60036. MR 0167642. LCCN 65-12253.


  3. P.K. Suetin (2001) [1994], "Jacobi_polynomials", in Hazewinkel, Michiel, Encyclopedia of Mathematics, Springer Science+Business Media B.V. / Kluwer Academic Publishers, ISBN 978-1-55608-010-4


  4. Biedenharn, L.C.; Louck, J.D. (1981). Angular Momentum in Quantum Physics. Reading: Addison-Wesley. 



Lecturas relacionadas




  • Andrews, George E.; Askey, Richard; Roy, Ranjan (1999), Special functions, Encyclopedia of Mathematics and its Applications 71, Cambridge University Press, ISBN 978-0-521-62321-6, MR 1688958, ISBN 978-0-521-78988-2 


  • Koornwinder, Tom H.; Wong, Roderick S. C.; Koekoek, Roelof; Swarttouw, René F. (2010), Cambridge University Press, ed., Orthogonal Polynomials, ISBN 978-0521192255, MR 2723248 



Enlaces externos



  • Weisstein, Eric W. «Jacobi Polynomial». En Weisstein, Eric W. MathWorld (en inglés). Wolfram Research. 



Popular posts from this blog

Plaza Victoria

In PowerPoint, is there a keyboard shortcut for bulleted / numbered list?

How to put 3 figures in Latex with 2 figures side by side and 1 below these side by side images but in...