Articulo de referencia

funciones elípticas de Jacobi

En matemáticas , las funciones elípticas de Jacobi son un conjunto de funciones elípticas básicas . Se encuentran en la descripción del movimiento de un péndulo , así como en el...

En matemáticas , las funciones elípticas de Jacobi son un conjunto de funciones elípticas básicas . Se encuentran en la descripción del movimiento de un péndulo , así como en el diseño de filtros elípticos electrónicos . Mientras que las funciones trigonométricas se definen con referencia a un círculo, las funciones elípticas de Jacobi son una generalización que se refiere a otras secciones cónicas , en particular la elipse. La relación con las funciones trigonométricas está contenida en la notación, por ejemplo, mediante la notación correspondiente.sn{\displaystyle \operatorname {sn} }parapecado{\displaystyle \sin }Las funciones elípticas de Jacobi se utilizan con mayor frecuencia en problemas prácticos que las funciones elípticas de Weierstrass, ya que no requieren la definición ni la comprensión de conceptos de análisis complejo . Fueron introducidas por Carl Gustav Jakob Jacobi ( 1829 ) . Carl Friedrich Gauss ya había estudiado funciones elípticas de Jacobi especiales en 1797, en particular las funciones elípticas lemniscatas [ 1 ] , pero su trabajo se publicó mucho más tarde. 

Descripción general

El rectángulo fundamental en el plano complejo de{\displaystyle u}

Existen doce funciones elípticas de Jacobi denotadas porpq(,metro){\displaystyle \operatorname {pq} (u,m)}, dóndepag{\displaystyle \mathrm {p} }yq{\displaystyle \mathrm {q} }¿Alguna de las letras?do{\displaystyle \mathrm {c} },s{\displaystyle \mathrm {s} },norte{\displaystyle \mathrm {n} }, yd{\displaystyle \mathrm {d} }. (Funciones de la formapáginas(,metro){\displaystyle \operatorname {pp} (u,m)}se establecen trivialmente a la unidad para mayor completitud notacional.){\displaystyle u}ese es el argumento, ymetro{\displaystyle m}es el parámetro, ambos pueden ser complejos. De hecho, las funciones elípticas de Jacobi son meromorfas en ambos{\displaystyle u}ymetro{\displaystyle m}. [ 2 ] La distribución de los ceros y polos en el{\displaystyle u}El plano -es bien conocido. Sin embargo, existen cuestiones sobre la distribución de los ceros y polos en elmetro{\displaystyle m}-plane aún por investigar. [ 2 ]

En el plano complejo del argumento{\displaystyle u}, las doce funciones forman una red repetitiva de polos y ceros simples . [ 3 ] Dependiendo de la función, un paralelogramo repetitivo , o celda unitaria, tendrá lados de longitud2K{\displaystyle 2K}o4K{\displaystyle 4K}en el eje real, y2K{\displaystyle 2K'}o4K{\displaystyle 4K'}en el eje imaginario, dondeK=K(metro){\displaystyle K=K(m)}yK=K(1metro){\displaystyle K'=K(1-m)}se conocen como los períodos trimestrales conK(){\displaystyle K(\cdot )}siendo la integral elíptica de primera especie. La naturaleza de la celda unitaria se puede determinar inspeccionando el "rectángulo auxiliar" (generalmente un paralelogramo), que es un rectángulo formado por el origen.(0,0){\displaystyle (0,0)}en una esquina, y (K,K){\displaystyle (K,K')}como la esquina diagonalmente opuesta. Como en el diagrama, las cuatro esquinas del rectángulo auxiliar se denominans{\displaystyle \mathrm {s} },do{\displaystyle \mathrm {c} },d{\displaystyle \mathrm {d} }, ynorte{\displaystyle \mathrm {n} }, yendo en sentido contrario a las agujas del reloj desde el origen. La funciónpq(,metro){\displaystyle \operatorname {pq} (u,m)}tendrá un cero en elpag{\displaystyle \mathrm {p} }esquina y un poste en elq{\displaystyle \mathrm {q} }esquina. Las doce funciones corresponden a las doce maneras de disponer estos polos y ceros en las esquinas del rectángulo.

Cuando el argumento{\displaystyle u}y parámetrometro{\displaystyle m}son reales, con0<metro<1{\displaystyle 0<m<1},K{\displaystyle K}yK{\displaystyle K'}será real y el paralelogramo auxiliar será de hecho un rectángulo, y las funciones elípticas de Jacobi serán todas de valor real en la recta real.

Dado que las funciones elípticas de Jacobi son doblemente periódicas en{\displaystyle u}, se factorizan a través de un toro ; en efecto, su dominio puede considerarse un toro, al igual que el coseno y el seno se definen en un círculo. En lugar de tener un solo círculo, ahora tenemos el producto de dos círculos, uno real y otro imaginario. El plano complejo puede ser reemplazado por un toro complejo . La circunferencia del primer círculo es4K{\displaystyle 4K}y el segundo4K{\displaystyle 4K'}, dóndeK{\displaystyle K}yK{\displaystyle K'}son los cuartos de período . Cada función tiene dos ceros y dos polos en posiciones opuestas en el toro. Entre los puntos0{\displaystyle 0},K{\displaystyle K},K+iK{\displaystyle K+iK'},iK{\displaystyle iK'}Hay un cero y un polo.

Las funciones elípticas de Jacobi son entonces funciones meromorfas doblemente periódicas que satisfacen las siguientes propiedades:

  • Hay un simple cero en la esquinapag{\displaystyle \mathrm {p} }y un simple poste en la esquina q{\displaystyle \mathrm {q} }.
  • El número complejopagq{\displaystyle \mathrm {p} -\mathrm {q} }es igual a la mitad del período de la funciónpq{\displaystyle \operatorname {pq} u}; es decir, la funciónpq{\displaystyle \operatorname {pq} u}es periódica en la direcciónpq{\displaystyle \operatorname {pq} }, siendo el período el2(pagq){\displaystyle 2(\mathrm {p} -\mathrm {q} )}. La funciónpq{\displaystyle \operatorname {pq} u}También es periódica en las otras dos direcciones.pagpag{\displaystyle \mathrm {pp} '}ypagq{\displaystyle \mathrm {pq} '}, con períodos tales quepagpag{\displaystyle \mathrm {p} -\mathrm {p} '}ypagq{\displaystyle \mathrm {p} -\mathrm {q} '}son períodos trimestrales.
alt=Función elíptica de Jacobi sn {\displaystyle \operatorname {sn} } , k = 0.8 {\displaystyle k=0.8}
función elíptica de Jacobisn{\displaystyle \operatorname {sn} }
alt=Función elíptica de Jacobi cn {\displaystyle \operatorname {cn} } , k = 0.8 {\displaystyle k=0.8}
función elíptica de Jacobicn{\displaystyle \operatorname {cn} }
alt=Función elíptica de Jacobi dn {\displaystyle \operatorname {dn} } , k = 0.8 {\displaystyle k=0.8}
función elíptica de Jacobidn{\displaystyle \operatorname {dn} }
alt=Función elíptica de Jacobi sc {\displaystyle \operatorname {sc} } , k = 0.8 {\displaystyle k=0.8}
función elíptica de JacobiCarolina del Sur{\displaystyle \operatorname {sc} }
Gráficas de cuatro funciones elípticas de Jacobi en el plano complejo de{\displaystyle u}, ilustrando su comportamiento doblemente periódico. Imágenes generadas utilizando una versión del método de coloración de dominio . [ 4 ] Todos tienen valores dek=metro{\displaystyle k={\sqrt {m}}}igual a0,8{\displaystyle 0.8}.

Notación

Las funciones elípticas se pueden expresar en diversas notaciones, lo que puede hacer que el tema resulte innecesariamente confuso. Las funciones elípticas son funciones de dos variables. La primera variable se puede expresar en términos de la amplitud.φ{\displaystyle \varphi }, o más comúnmente, en términos de{\displaystyle u}dado a continuación. La segunda variable podría darse en términos del parámetrometro{\displaystyle m}, o como el módulo elípticok{\displaystyle k}, dóndek2=metro{\displaystyle k^{2}=m}, o en términos del ángulo modularα{\displaystyle \alpha }, dóndemetro=pecado2α{\displaystyle m=\sin ^{2}\alpha }. Los complementos dek{\displaystyle k}ymetro{\displaystyle m}se definen comometro=1metro{\displaystyle m'=1-m}yk=metro{\textstyle k'={\sqrt {m'}}}Estos cuatro términos se utilizan a continuación sin comentarios para simplificar diversas expresiones.

Las doce funciones elípticas de Jacobi se escriben generalmente comopq(,metro){\displaystyle \operatorname {pq} (u,m)}dóndepag{\displaystyle \mathrm {p} }yq{\displaystyle \mathrm {q} }¿Alguna de las letras?do{\displaystyle \mathrm {c} },s{\displaystyle \mathrm {s} },norte{\displaystyle \mathrm {n} }, yd{\displaystyle \mathrm {d} }. Funciones de la formapáginas(,metro){\displaystyle \operatorname {pp} (u,m)}se establecen trivialmente a la unidad para una completitud notacional. Las funciones “principales” generalmente se toman comocn(,metro){\displaystyle \operatorname {cn} (u,m)},sn(,metro){\displaystyle \operatorname {sn} (u,m)}ydn(,metro){\displaystyle \operatorname {dn} (u,m)}A partir de estas tres funciones se pueden derivar todas las demás, y las expresiones suelen escribirse únicamente en términos de ellas; sin embargo, diversas simetrías y generalizaciones se expresan con mayor facilidad utilizando el conjunto completo. (Esta notación se debe a Gudermann y Glaisher y no es la notación original de Jacobi).

A lo largo de este artículo,pq(,t2)=pq(;t){\displaystyle \operatorname {pq} (u,t^{2})=\operatorname {pq} (u;t)}.

Las funciones están relacionadas entre sí mediante la regla de multiplicación: (argumentos suprimidos)

pqpagq=pagqpagq{\displaystyle \operatorname {pq} \cdot \operatorname {p'q'} =\operatorname {pq'} \cdot \operatorname {p'q} }

de las cuales se pueden derivar otras relaciones de uso común:

prcódigo QR=pq{\displaystyle {\frac {\operatorname {pr} }{\operatorname {qr} }}=\operatorname {pq} }prrq=pq{\displaystyle \operatorname {pr} \cdot \operatorname {rq} =\operatorname {pq} }1qp=pq{\displaystyle {\frac {1}{\operatorname {qp} }}=\operatorname {pq} }

La regla de multiplicación se deduce inmediatamente de la identificación de las funciones elípticas con las funciones theta de Neville [ 5 ].

pq(,metro)=θpag(,metro)θq(,metro){\displaystyle \operatorname {pq} (u,m)={\frac {\theta _{\operatorname {p} }(u,m)}{\theta _{\operatorname {q} }(u,m)}}} Tenga en cuenta también que:

K(metro)=K(k2)=01dt(1t2)(1metrot2)=01dt(1t2)(1k2t2).{\displaystyle K(m)=K(k^{2})=\int _{0}^{1}{\frac {dt}{\sqrt {(1-t^{2})(1-mt^{2})}}}=\int _{0}^{1}{\frac {dt}{\sqrt {(1-t^{2})(1-k^{2}t^{2})}}}.}

Definición en términos de inversas de integrales elípticas

Modelo de la amplitud de Jacobi (medida a lo largo del eje vertical) en función de las variables independientes u y el módulo k.

Existe una definición que relaciona las funciones elípticas con la inversa de la integral elíptica incompleta de primera especie.F{\displaystyle F}Estas funciones toman los parámetros{\displaystyle u}ymetro{\displaystyle m}como entradas. Elφ{\displaystyle \varphi }que satisface

=F(φ|metro)=0φdθ1metropecado2θ{\displaystyle u=F(\varphi |m)=\int _{0}^{\varphi }{\frac {\mathrm {d} \theta }{\sqrt {1-m\sin ^{2}\theta }}}}

se denomina amplitud de Jacobi :

soy(|metro)=φ.{\displaystyle \operatorname {am} (u|m)=\varphi .}

En este marco, el seno elíptico sn u (en latín: sinus amplitudinis ) viene dado por 

sn(|metro)=pecadosoy(|metro){\displaystyle \operatorname {sn} (u|m)=\sin \operatorname {am} (u|m)}

y el coseno elíptico cn u (en latín: cosinus amplitudinis ) viene dado por 

cn(|metro)=porquesoy(|metro){\displaystyle \operatorname {cn} (u|m)=\cos \operatorname {am} (u|m)}

y la amplitud delta dn u (latín: delta amplitudinis ) [ nota 1 ] 

dn(|metro)=ddsoy(|metro).{\displaystyle \operatorname {dn} (u|m)={\frac {\mathrm {d} }{\mathrm {d} u}}\operatorname {am} (u|m).} En lo anterior, el valormetro{\displaystyle m}es un parámetro libre , generalmente tomado como real tal que0metro1{\displaystyle 0\leq m\leq 1}(pero pueden ser complejas en general), y por lo tanto las funciones elípticas pueden pensarse como dadas por dos variables,{\displaystyle u}y el parámetro metro{\displaystyle m}. Las nueve funciones elípticas restantes se construyen fácilmente a partir de las tres anteriores (sn{\displaystyle \operatorname {sn} },cn{\displaystyle \operatorname {cn} },dn{\displaystyle \operatorname {dn} }), y se dan en una sección más abajo. Tenga en cuenta que cuandoφ=π/2{\displaystyle \varphi =\pi /2}, eso{\displaystyle u}entonces es igual al cuarto período K{\displaystyle K}.

En el contexto más general,soy(|metro){\displaystyle \operatorname {am} (u|m)}es una función multivaluada (en{\displaystyle u}) con infinitos puntos de ramificación logarítmicos (las ramas difieren en múltiplos enteros de2π{\displaystyle 2\pi }), es decir, los puntos2sK(metro)+(4t+1)K(1metro)i{\displaystyle 2sK(m)+(4t+1)K(1-m)i}y2sK(metro)+(4t+3)K(1metro)i{\displaystyle 2sK(m)+(4t+3)K(1-m)i}dóndes,tZ{\displaystyle s,t\in \mathbb {Z} }. [ 6 ] Esta función multivaluada puede convertirse en univaluada cortando el plano complejo a lo largo de los segmentos de línea que unen estos puntos de ramificación (el corte puede hacerse de maneras no equivalentes, dando funciones univaluadas no equivalentes), haciendo asísoy(|metro){\displaystyle \operatorname {am} (u|m)}analítico en todas partes excepto en los cortes de rama . En contraste,pecadosoy(|metro){\displaystyle \sin \operatorname {am} (u|m)}y otras funciones elípticas no tienen puntos de ramificación, dan valores consistentes para cada rama desoy{\displaystyle \operatorname {am} }y son meromorfas en todo el plano complejo. Dado que toda función elíptica es meromorfa en todo el plano complejo (por definición),soy(|metro){\displaystyle \operatorname {am} (u|m)}(considerada como una función unívoca) no es una función elíptica.

Sin embargo, un recorte en particular parasoy(|metro){\displaystyle \operatorname {am} (u|m)}se puede hacer en el{\displaystyle u}-plano por segmentos de línea desde2sK(metro)+(4t+1)K(1metro)i{\displaystyle 2sK(m)+(4t+1)K(1-m)i}a2sK(metro)+(4t+3)K(1metro)i{\displaystyle 2sK(m)+(4t+3)K(1-m)i}cons,tZ{\displaystyle s,t\in \mathbb {Z} }; entonces solo queda definirsoy(|metro){\displaystyle \operatorname {am} (u|m)}en la rama corta por continuidad desde alguna dirección. Entoncessoy(|metro){\displaystyle \operatorname {am} (u|m)}se vuelve unívoca y singularmente periódica en{\displaystyle u}con el período mínimo4iK(1metro){\displaystyle 4iK(1-m)}y tiene singularidades en los puntos de ramificación logarítmica mencionados anteriormente. SimetroR{\displaystyle m\in \mathbb {R} }ymetro1{\displaystyle m\leq 1},soy(|metro){\displaystyle \operatorname {am} (u|m)}es continuo en{\displaystyle u}en la línea real. Cuandometro>1{\displaystyle m>1}, los cortes de ramas desoy(|metro){\displaystyle \operatorname {am} (u|m)}en el{\displaystyle u}-el avión cruza la línea real en2(2s+1)K(1/metro)/metro{\displaystyle 2(2s+1)K(1/m)/{\sqrt {m}}}parasZ{\displaystyle s\in \mathbb {Z} }; por lo tanto parametro>1{\displaystyle m>1},soy(|metro){\displaystyle \operatorname {am} (u|m)}no es continuo en{\displaystyle u}en la línea real y salta por2π{\displaystyle 2\pi }en las discontinuidades.

Pero definirsoy(|metro){\displaystyle \operatorname {am} (u|m)}De esta forma se producen cortes de ramas muy complicados en elmetro{\displaystyle m}-avión ( no el{\displaystyle u}-avión); aún no han sido descritos completamente.

Dejar mi(φ|metro)=0φ1metropecado2θdθ{\displaystyle E(\varphi |m)=\int _{0}^{\varphi }{\sqrt {1-m\sin ^{2}\theta }}\,\mathrm {d} \theta } sea ​​la integral elíptica incompleta de segundo tipo con parámetrometro{\displaystyle m}.

Entonces, la función épsilon de Jacobi se puede definir como mi(|metro)=mi(soy(|metro)|metro){\displaystyle {\mathcal {E}}(u|m)=E(\operatorname {am} (u|m)|m)} paraR{\displaystyle u\in \mathbb {R} }y0<metro<1{\displaystyle 0<m<1}y por continuación analítica en cada una de las variables en caso contrario: la función épsilon de Jacobi es meromorfa en todo el plano complejo (en ambos{\displaystyle u}ymetro{\displaystyle m}). Alternativamente, a lo largo de ambos{\displaystyle u}-avión ymetro{\displaystyle m}-plano, [ 7 ]mi(|metro)=0dn2(t|metro)dt;{\displaystyle {\mathcal {E}}(u|m)=\int _{0}^{u}\operatorname {dn} ^{2}(t|m)\,\mathrm {d} t;}mi{\displaystyle {\mathcal {E}}}está bien definido de esta manera porque todos los residuos detdn(t|metro)2{\displaystyle t\mapsto \operatorname {dn} (t|m)^{2}}son cero, por lo que la integral es independiente de la trayectoria. Así, el épsilon de Jacobi relaciona la integral elíptica incompleta de primera especie con la integral elíptica incompleta de segunda especie: mi(φ|metro)=mi(F(φ|metro)|metro).{\displaystyle E(\varphi |m)={\mathcal {E}}(F(\varphi |m)|m).} La función épsilon de Jacobi no es una función elíptica, pero aparece al diferenciar las funciones elípticas de Jacobi con respecto al parámetro.

La función zn de Jacobi se define por zinc(|metro)=mi(|metro)mi(metro)K(metro).{\displaystyle \operatorname {zn} (u|m)={\mathcal {E}}(u|m)-{\frac {E(m)}{K(m)}}u.} Es una función periódica simple que es meromorfa en{\displaystyle u}, pero no enmetro{\displaystyle m}(debido a los cortes de ramas demi{\displaystyle E}yK{\displaystyle K}). Su período mínimo en{\displaystyle u}es2K(metro){\displaystyle 2K(m)}Está relacionada con la función zeta de Jacobi porZ(φ|metro)=zinc(F(φ|metro)|metro).{\displaystyle Z(\varphi |m)=\operatorname {zn} (F(\varphi |m)|m).}

Históricamente, las funciones elípticas de Jacobi se definieron inicialmente utilizando la amplitud. En textos más modernos sobre funciones elípticas, estas se definen por otros medios, por ejemplo, mediante cocientes de funciones theta (véase más adelante), y se ignora la amplitud.

En términos modernos, la relación con las integrales elípticas se expresaría mediantesn(F(φ|metro)|metro)=pecadoφ{\displaystyle \operatorname {sn} (F(\varphi |m)|m)=\sin \varphi }(ocn(F(φ|metro)|metro)=porqueφ{\displaystyle \operatorname {cn} (F(\varphi |m)|m)=\cos \varphi }) en lugar desoy(F(φ|metro)|metro)=φ{\displaystyle \operatorname {am} (F(\varphi |m)|m)=\varphi }.

Definición trigonométrica: la elipse de Jacobi

Gráfico de la elipse de Jacobi ( x 2  + y 2 / b 2 = 1, b real) y las doce funciones elípticas de Jacobi pq ( u , m ) para valores particulares del ángulo φ y el parámetro b . La curva continua es la elipse, con m = 1 − 1/ b 2 y u = F ( φ , m ) donde F ( , ) es la integral elíptica de primera especie (con parámetro           metro=k2{\displaystyle m=k^{2}}La curva punteada es el círculo unitario. Las líneas tangentes del círculo y la elipse en x  =  cd que cruzan el eje x en dc se muestran en gris claro.

porqueφ,pecadoφ{\displaystyle \cos \varphi ,\,\sin \varphi }se definen en el círculo unitario con radior=1{\displaystyle r=1}y ánguloφ={\displaystyle \varphi =}longitud de arco del círculo unitario medida desde el eje x positivo . De manera similar, las funciones elípticas de Jacobi se definen en la elipse unitaria cona=1{\displaystyle a=1}yb1{\displaystyle b\geq 1}. Dejar

incógnita2+y2b2=1,metro=11b2,incógnita=rporqueφ,y=rpecadoφ.{\displaystyle {\begin{aligned}&x^{2}+{\frac {y^{2}}{b^{2}}}=1,\quad m=1-{\frac {1}{b^{2}}},\\&x=r\cos \varphi ,\quad y=r\sin \varphi .\end{aligned}}}

Entonces0metro<1{\displaystyle 0\leq m<1}y

r(φ,metro)=11metropecado2φ.{\displaystyle r(\varphi ,m)={\frac {1}{\sqrt {1-m\sin ^{2}\varphi }}}\,.}

Para cada ánguloφ{\displaystyle \varphi }el parámetro =(φ,metro)=0φr(θ,metro)dθ{\displaystyle u=u(\varphi ,m)=\int _{0}^{\varphi }r(\theta ,m)\,d\theta } (la integral elíptica incompleta de primera especie) se calcula. En el círculo unitario (a=b=1{\displaystyle a=b=1}),{\displaystyle u}sería una longitud de arco. Sin embargo, la relación de{\displaystyle u}Calcular la longitud del arco de una elipse es más complicado. [ 8 ]

DejarPAG=(incógnita,y)=(rporqueφ,rpecadoφ){\displaystyle P=(x,y)=(r\cos \varphi ,r\sin \varphi )}Sea un punto en la elipse y dejemos quePAG=(incógnita,y)=(porqueφ,pecadoφ){\displaystyle P'=(x',y')=(\cos \varphi ,\sin \varphi )}sea ​​el punto donde el círculo unitario interseca la línea entrePAG{\displaystyle P}y el origenO{\displaystyle O}Luego, las relaciones familiares del círculo unitario

incógnita=porqueφ,y=pecadoφ{\displaystyle x'=\cos \varphi ,\quad y'=\sin \varphi }

leer para la elipsis

incógnita=cn(,metro),y=sn(,metro).{\displaystyle x'=\operatorname {cn} (u,m),\quad y'=\operatorname {sn} (u,m).}

Entonces, las proyecciones del punto de intersecciónPAG{\displaystyle P'}de la líneaOPAG{\displaystyle OP}con el círculo unitario en los ejes x e y son simplementecn(,metro){\displaystyle \operatorname {cn} (u,m)}ysn(,metro){\displaystyle \operatorname {sn} (u,m)}Estas proyecciones pueden interpretarse como una "definición como trigonometría". En resumen,

cn(,metro)=incógnitar(φ,metro),sn(,metro)=yr(φ,metro),dn(,metro)=1r(φ,metro).{\displaystyle \operatorname {cn} (u,m)={\frac {x}{r(\varphi ,m)}},\quad \operatorname {sn} (u,m)={\frac {y}{r(\varphi ,m)}},\quad \operatorname {dn} (u,m)={\frac {1}{r(\varphi ,m)}}.}

Para elincógnita{\displaystyle x}yy{\displaystyle y}valores del puntoPAG{\displaystyle P}con {\displaystyle u}y parámetrometro{\displaystyle m}obtenemos, insertando la relación

r(φ,metro)=1dn(,metro){\displaystyle r(\varphi ,m)={\frac {1}{\operatorname {dn} (u,m)}}}

enincógnita=r(φ,metro)porqueφ{\displaystyle x=r(\varphi ,m)\cos \varphi }yy=r(φ,metro)pecadoφ{\displaystyle y=r(\varphi ,m)\sin \varphi }, encontramos

incógnita=cn(,metro)dn(,metro),y=sn(,metro)dn(,metro).{\displaystyle x={\frac {\operatorname {cn} (u,m)}{\operatorname {dn} (u,m)}},\quad y={\frac {\operatorname {sn} (u,m)}{\operatorname {dn} (u,m)}}.}

Las últimas relaciones para las coordenadas x e y de los puntos en la elipse unitaria pueden considerarse como una generalización de las relacionesincógnita=porqueφ{\displaystyle x=\cos \varphi }yy=pecadoφ{\displaystyle y=\sin \varphi }para las coordenadas de los puntos en el círculo unitario.

La siguiente tabla resume las expresiones para todas las funciones elípticas de Jacobi pq(u,m) en las variables ( x , y , r ) y ( φ ,dn) conr=incógnita2+y2{\textstyle r={\sqrt {x^{2}+y^{2}}}}.

Definición en términos de las funciones theta de Jacobi

Utilizando integrales elípticas

De forma equivalente, las funciones elípticas de Jacobi pueden definirse en términos de las funciones theta . [ 9 ] Conz,τdo{\displaystyle z,\tau \in \mathbb {C} }de tal manera queSoyτ>0{\displaystyle \operatorname {Im} \tau >0}, dejar

θ1(z|τ)=norte=(1)norte12mi(2norte+1)iz+πiτ(norte+12)2,{\displaystyle \theta _{1}(z|\tau )=\displaystyle \sum _{n=-\infty }^{\infty }(-1)^{n-{\frac {1}{2}}}e^{(2n+1)iz+\pi i\tau \left(n+{\frac {1}{2}}\right)^{2}},}θ2(z|τ)=norte=mi(2norte+1)iz+πiτ(norte+12)2,{\displaystyle \theta _{2}(z|\tau )=\displaystyle \sum _{n=-\infty }^{\infty }e^{(2n+1)iz+\pi i\tau \left(n+{\frac {1}{2}}\right)^{2}},}θ3(z|τ)=norte=mi2norteiz+πiτnorte2,{\displaystyle \theta _{3}(z|\tau )=\displaystyle \sum _{n=-\infty }^{\infty }e^{2niz+\pi i\tau n^{2}},}θ4(z|τ)=norte=(1)nortemi2norteiz+πiτnorte2{\displaystyle \theta _{4}(z|\tau )=\displaystyle \sum _{n=-\infty }^{\infty }(-1)^{n}e^{2niz+\pi i\tau n^{2}}} y dejarθ2(τ)=θ2(0|τ){\displaystyle \theta _{2}(\tau )=\theta _{2}(0|\tau )},θ3(τ)=θ3(0|τ){\displaystyle \theta _{3}(\tau )=\theta _{3}(0|\tau )},θ4(τ)=θ4(0|τ){\displaystyle \theta _{4}(\tau )=\theta _{4}(0|\tau )}. Luego conK=K(metro){\displaystyle K=K(m)},K=K(1metro){\displaystyle K'=K(1-m)},ζ=π/(2K){\displaystyle \zeta =\pi u/(2K)}yτ=iK/K{\displaystyle \tau =iK'/K},

sn(,metro)=θ3(τ)θ1(ζ|τ)θ2(τ)θ4(ζ|τ),cn(,metro)=θ4(τ)θ2(ζ|τ)θ2(τ)θ4(ζ|τ),dn(,metro)=θ4(τ)θ3(ζ|τ)θ3(τ)θ4(ζ|τ).{\displaystyle {\begin{aligned}\operatorname {sn} (u,m)&={\frac {\theta _{3}(\tau )\theta _{1}(\zeta |\tau )}{\theta _{2}(\tau )\theta _{4}(\zeta |\tau )}},\\\operatorname {cn} (u,m)&={\frac {\theta _{4}(\tau )\theta _{2}(\zeta |\tau )}{\theta _{2}(\tau )\theta _{4}(\zeta |\tau )}},\\\operatorname {dn} (u,m)&={\frac {\theta _{4}(\tau )\theta _{3}(\zeta |\tau )}{\theta _{3}(\tau )\theta _{4}(\zeta |\tau )}}.\end{aligned}}}

La función zn de Jacobi también puede expresarse mediante funciones theta: zinc(,metro)=π2Kθ4(ζ|τ)θ4(ζ|τ)=π2Kθ3(ζ|τ)θ3(ζ|τ)+metrosn(,metro)cn(,metro)dn(,metro)=π2Kθ2(ζ|τ)θ2(ζ|τ)+dn(,metro)sn(,metro)cn(,metro)=π2Kθ1(ζ|τ)θ1(ζ|τ)cn(,metro)dn(,metro)sn(,metro){\displaystyle {\begin{aligned}\operatorname {zn} (u,m)&={\frac {\pi }{2K}}{\frac {\theta _{4}'(\zeta |\tau )}{\theta _{4}(\zeta |\tau )}}\\&={\frac {\pi }{2K}}{\frac {\theta _{3}'(\zeta |\tau )}{\theta _{3}(\zeta |\tau )}}+m{\frac {\operatorname {sn} (u,m)\operatorname {cn} (u,m)}{\operatorname {dn} (u,m)}}\\&={\frac {\pi }{2K}}{\frac {\theta _{2}'(\zeta |\tau )}{\theta _{2}(\zeta |\tau )}}+{\frac {\operatorname {dn} (u,m)\operatorname {sn} (u,m)}{\operatorname {cn} (u,m)}}\\&={\frac {\pi }{2K}}{\frac {\theta _{1}'(\zeta |\tau )}{\theta _{1}(\zeta |\tau )}}-{\frac {\operatorname {cn} (u,m)\operatorname {dn} (u,m)}{\operatorname {sn} (u,m)}}\end{aligned}}} dónde{\displaystyle '}denota la derivada parcial con respecto a la primera variable.

Utilizando la inversión modular

De hecho, la definición de las funciones elípticas de Jacobi en Whittaker & Watson se enuncia de forma ligeramente diferente a la dada anteriormente (pero es equivalente a ella) y se basa en la inversión modular: La funciónλ{\displaystyle \lambda }, definido por

La regiónF1{\displaystyle F_{1}}en el plano complejo. Está delimitado por dos semicírculos desde abajo, por un rayo desde la izquierda y por un rayo desde la derecha.

λ(τ)=θ2(τ)4θ3(τ)4,{\displaystyle \lambda (\tau )={\frac {\theta _{2}(\tau )^{4}}{\theta _{3}(\tau )^{4}}},} asume cada valor endo{0,1}{\displaystyle \mathbb {C} -\{0,1\}}una y solo una vez [ 10 ] en F1(F1{τH:Reτ<0}){\displaystyle F_{1}-(\partial F_{1}\cap \{\tau \in \mathbb {H} :\operatorname {Re} \tau <0\})} dondeH{\displaystyle \mathbb {H} }es el semiplano superior en el plano complejo,F1{\displaystyle \partial F_{1}}es el límite deF1{\displaystyle F_{1}}y F1={τH:|Reτ|1,|Re(1/τ)|1}.{\displaystyle F_{1}=\{\tau \in \mathbb {H} :\left|\operatorname {Re} \tau \right|\leq 1,\left|\operatorname {Re} (1/\tau )\right|\leq 1\}.} De esta manera, cadametro=definiciónλ(τ)do{0,1}{\displaystyle m\,{\overset {\text{def}}{=}}\,\lambda (\tau )\in \mathbb {C} -\{0,1\}}puede estar asociado con uno y solo unoτ{\displaystyle \tau }. Luego, Whittaker y Watson definen las funciones elípticas de Jacobi mediante sn(,metro)=θ3(τ)θ1(ζ|τ)θ2(τ)θ4(ζ|τ),cn(,metro)=θ4(τ)θ2(ζ|τ)θ2(τ)θ4(ζ|τ),dn(,metro)=θ4(τ)θ3(ζ|τ)θ3(τ)θ4(ζ|τ){\displaystyle {\begin{aligned}\operatorname {sn} (u,m)&={\frac {\theta _{3}(\tau )\theta _{1}(\zeta |\tau )}{\theta _{2}(\tau )\theta _{4}(\zeta |\tau )}},\\\operatorname {cn} (u,m)&={\frac {\theta _{4}(\tau )\theta _{2}(\zeta |\tau )}{\theta _{2}(\tau )\theta _{4}(\zeta |\tau )}},\\\operatorname {dn} (u,m)&={\frac {\theta _{4}(\tau )\theta _{3}(\zeta |\tau )}{\theta _{3}(\tau )\theta _{4}(\zeta |\tau )}}\end{aligned}}} dóndeζ=/θ3(τ)2{\displaystyle \zeta =u/\theta _{3}(\tau )^{2}}En el libro, imponen una restricción adicional ametro{\displaystyle m}(esometro(,0)(1,){\displaystyle m\notin (-\infty ,0)\cup (1,\infty )}), pero de hecho no es una restricción necesaria (véase la referencia de Cox). Además, simetro=0{\displaystyle m=0}ometro=1{\displaystyle m=1}Las funciones elípticas de Jacobi degeneran en funciones no elípticas, lo cual se describe a continuación.

Definición en términos de funciones theta de Neville

Las funciones elípticas de Jacobi se pueden definir de forma muy sencilla utilizando las funciones theta de Neville : [ 11 ]

pq(,metro)=θpag(,metro)θq(,metro){\displaystyle \operatorname {pq} (u,m)={\frac {\theta _{\operatorname {p} }(u,m)}{\theta _{\operatorname {q} }(u,m)}}}

A menudo, la simplificación de productos complejos de las funciones elípticas de Jacobi se facilita mediante el uso de estas identidades.

transformaciones de Jacobi

Las transformaciones imaginarias de Jacobi

Gráfico de la curva de Jacobi degenerada ( x 2  + y 2 / b 2 = 1, b = ∞) y las doce funciones elípticas de Jacobi pq( u ,1) para un valor particular del ángulo φ . La curva continua es la elipse degenerada ( x 2 = 1) con m = 1 y u = F ( φ ,1) donde F ( , ) es la integral elíptica de primera especie. La curva punteada es el círculo unitario. Dado que estas son las funciones de Jacobi para m = 0 (funciones trigonométricas circulares) pero con argumentos imaginarios, corresponden a las seis funciones trigonométricas hiperbólicas.              

Las transformaciones imaginarias de Jacobi relacionan diversas funciones de la variable imaginaria iu o, equivalentemente, relaciones entre diversos valores del parámetro m . En términos de las funciones principales: [ 12 ] : 506

cn(,metro)=Carolina del Norte(i,1metro){\displaystyle \operatorname {cn} (u,m)=\operatorname {nc} (i\,u,1\!-\!m)}sn(,metro)=iCarolina del Sur(i,1metro){\displaystyle \operatorname {sn} (u,m)=-i\operatorname {sc} (i\,u,1\!-\!m)}dn(,metro)=corriente continua(i,1metro){\displaystyle \operatorname {dn} (u,m)=\operatorname {dc} (i\,u,1\!-\!m)}

Utilizando la regla de la multiplicación, todas las demás funciones pueden expresarse en términos de las tres anteriores. Las transformaciones pueden escribirse generalmente comopq(,metro)=γpqpq(i,1metro){\displaystyle \operatorname {pq} (u,m)=\gamma _{\operatorname {pq} }\operatorname {pq} '(i\,u,1\!-\!m)}La siguiente tabla muestra laγpqpq(i,1metro){\displaystyle \gamma _{\operatorname {pq} }\operatorname {pq} '(i\,u,1\!-\!m)}para el pq( u,m ) especificado. [ 11 ] (Los argumentos(i,1metro){\displaystyle (i\,u,1\!-\!m)}se suprimen)

Dado que las funciones trigonométricas hiperbólicas son proporcionales a las funciones trigonométricas circulares con argumentos imaginarios, se deduce que las funciones de Jacobi producirán las funciones hiperbólicas para m=1. [ 5 ] : 249 En la figura, la curva de Jacobi se ha degenerado en dos líneas verticales en x  =  1 y x  =  −1.

Las transformaciones reales de Jacobi

Las transformaciones reales de Jacobi [ 5 ] : 308 proporcionan expresiones para las funciones elípticas en términos de valores alternos de m . Las transformaciones pueden escribirse generalmente comopq(,metro)=γpqpq(k,1/metro){\displaystyle \operatorname {pq} (u,m)=\gamma _{\operatorname {pq} }\operatorname {pq} '(k\,u,1/m)}La siguiente tabla muestra laγpqpq(k,1/metro){\displaystyle \gamma _{\operatorname {pq} }\operatorname {pq} '(k\,u,1/m)}para el pq( u,m ) especificado. [ 11 ] (Los argumentos(k,1/metro){\displaystyle (k\,u,1/m)}se suprimen)

Otras transformaciones de Jacobi

Las transformaciones reales e imaginarias de Jacobi se pueden combinar de varias maneras para obtener tres transformaciones más simples. [ 5 ] : 214 Las transformaciones reales e imaginarias son dos transformaciones en un grupo ( D 3 o grupo anarmónico ) de seis transformaciones. Si

μR(metro)=1/metro{\displaystyle \mu _{R}(m)=1/m}

es la transformación para el parámetro m en la transformación real, y

μI(metro)=1metro=metro{\displaystyle \mu _{I}(m)=1-m=m'}

es la transformación de m en la transformación imaginaria, entonces las otras transformaciones se pueden construir mediante la aplicación sucesiva de estas dos transformaciones básicas, lo que produce solo tres posibilidades más:

μIR(metro)=μI(μR(metro))=metro/metroμRI(metro)=μR(μI(metro))=1/metroμRIR(metro)=μR(μI(μR(metro)))=metro/metro{\displaystyle {\begin{aligned}\mu _{IR}(m)&=&\mu _{I}(\mu _{R}(m))&=&-m'/m\\\mu _{RI}(m)&=&\mu _{R}(\mu _{I}(m))&=&1/m'\\\mu _{RIR}(m)&=&\mu _{R}(\mu _{I}(\mu _{R}(m)))&=&-m/m'\end{aligned}}}

Estas cinco transformaciones, junto con la transformación identidad ( μ U ( m )  = m ), dan como resultado el grupo de seis elementos. Con respecto a las funciones elípticas de Jacobi, la transformación general se puede expresar utilizando solo tres funciones: 

cs(,metro)=γidos(γi,μi(metro)){\displaystyle \operatorname {cs} (u,m)=\gamma _{i}\operatorname {cs'} (\gamma _{i}u,\mu _{i}(m))}ns(,metro)=γinortes(γi,μi(metro)){\displaystyle \operatorname {ns} (u,m)=\gamma _{i}\operatorname {ns'} (\gamma _{i}u,\mu _{i}(m))}ds(,metro)=γids(γi,μi(metro)){\displaystyle \operatorname {ds} (u,m)=\gamma _{i}\operatorname {ds'} (\gamma _{i}u,\mu _{i}(m))}

donde i = U, I, IR, R, RI o RIR, identificando la transformación, γ i es un factor de multiplicación común a estas tres funciones, y la prima indica la función transformada. Las otras nueve funciones transformadas se pueden construir a partir de las tres anteriores. La razón por la que se eligieron las funciones cs, ns y ds para representar la transformación es que las demás funciones serán cocientes de estas tres (excepto sus inversas) y los factores de multiplicación se cancelarán.

La siguiente tabla enumera los factores de multiplicación para las tres funciones ps, los m transformados y los nombres de las funciones transformadas para cada una de las seis transformaciones. [ 5 ] : 214 (Como es habitual, k 2  = m , 1 − k 2 = k 1 2 = m y los argumentos (       γi,μi(metro){\displaystyle \gamma _{i}u,\mu _{i}(m)}) se suprimen)

Así, por ejemplo, podemos construir la siguiente tabla para la transformación RIR. [ 11 ] La transformación se escribe generalmentepq(,metro)=γpqpagq(k,metro/metro){\displaystyle \operatorname {pq} (u,m)=\gamma _{\operatorname {pq} }\,\operatorname {pq'} (k'\,u,-m/m')}(Los argumentos(k,metro/metro){\displaystyle (k'\,u,-m/m')}se suprimen)

El valor de las transformaciones de Jacobi es que cualquier conjunto de funciones elípticas de Jacobi con cualquier parámetro de valor real m puede convertirse en otro conjunto para el cual0<metro1/2{\displaystyle 0<m\leq 1/2}y, para valores reales de u , los valores de la función serán reales. [ 5 ] : pág. 215

Transformaciones de amplitud

A continuación, la segunda variable se suprime y es igual ametro{\displaystyle m}:

pecado(soy(+v)+soy(v))=2sncndnv1metrosn2sn2v,{\displaystyle \sin(\operatorname {am} (u+v)+\operatorname {am} (u-v))={\frac {2\operatorname {sn} u\operatorname {cn} u\operatorname {dn} v}{1-m\operatorname {sn} ^{2}u\operatorname {sn} ^{2}v}},}porque(soy(+v)soy(v))=cn2vsn2vdn21metrosn2sn2v{\displaystyle \cos(\operatorname {am} (u+v)-\operatorname {am} (u-v))={\dfrac {\operatorname {cn} ^{2}v-\operatorname {sn} ^{2}v\operatorname {dn} ^{2}u}{1-m\operatorname {sn} ^{2}u\operatorname {sn} ^{2}v}}} donde ambas identidades son válidas para todos,v,metrodo{\displaystyle u,v,m\in \mathbb {C} }de manera que ambas partes estén bien definidas.

Con

metro1=(1metro1+metro)2,{\displaystyle m_{1}=\left({\frac {1-{\sqrt {m'}}}{1+{\sqrt {m'}}}}\right)^{2},}

tenemos

porque(soy(,metro)+soy(K,metro))=sn((1metro),1/metro1),{\displaystyle \cos(\operatorname {am} (u,m)+\operatorname {am} (K-u,m))=-\operatorname {sn} ((1-{\sqrt {m'}})u,1/m_{1}),}pecado(soy(metro,metro/metro)+soy((1metro),1/metro1))=sn(,metro),{\displaystyle \sin(\operatorname {am} ({\sqrt {m'}}u,-m/m')+\operatorname {am} ((1-{\sqrt {m'}})u,1/m_{1}))=\operatorname {sn} (u,m),}pecado(soy((1+metro),metro1)+soy((1metro),1/metro1))=pecado(2soy(,metro)){\displaystyle \sin(\operatorname {am} ((1+{\sqrt {m'}})u,m_{1})+\operatorname {am} ((1-{\sqrt {m'}})u,1/m_{1}))=\sin(2\operatorname {am} (u,m))}

donde todas las identidades son válidas para todos,metrodo{\displaystyle u,m\in \mathbb {C} }de manera que ambas partes estén bien definidas.

La hipérbola de Jacobi

Gráfico de la hipérbola de Jacobi ( x 2  + y 2 / b 2 = 1, b imaginario) y las doce funciones elípticas de Jacobi pq( u , m ) para valores particulares del ángulo φ y el parámetro b . La curva continua es la hipérbola, con m = 1 − 1/ b 2 y u = F ( φ , m ) donde F ( , ) es la integral elíptica de primera especie. La curva punteada es el círculo unitario. Para el triángulo ds-dc, σ = sin( φ )cos( φ ).           

Al introducir los números complejos, nuestra elipse tiene asociada una hipérbola:

incógnita2y2b2=1{\displaystyle x^{2}-{\frac {y^{2}}{b^{2}}}=1} aplicando la transformación imaginaria de Jacobi [ 11 ] a las funciones elípticas en la ecuación anterior para x e y . 

incógnita=1dn(,1metro),y=sn(,1metro)dn(,1metro){\displaystyle x={\frac {1}{\operatorname {dn} (u,1-m)}},\quad y={\frac {\operatorname {sn} (u,1-m)}{\operatorname {dn} (u,1-m)}}}

De ello se deduce que podemos ponerincógnita=dn(,1metro),y=sn(,1metro){\displaystyle x=\operatorname {dn} (u,1-m),y=\operatorname {sn} (u,1-m)}. Así que nuestra elipse tiene una elipse dual con m reemplazado por 1-m. Esto lleva al toro complejo mencionado en la Introducción. [ 13 ] Generalmente, m puede ser un número complejo , pero cuando m es real y m<0, la curva es una elipse con eje mayor en la dirección x. En m=0 la curva es un círculo, y para 0<m<1, la curva es una elipse con eje mayor en la dirección y. En m  =  1, la curva degenera en dos líneas verticales en x  =  ±1. Para m  >  1, la curva es una hipérbola. Cuando m es complejo pero no real, x o y o ambos son complejos y la curva no puede describirse en un diagrama x - y real .

Funciones menores

Al invertir el orden de las dos letras del nombre de la función, se obtienen las funciones recíprocas de las tres funciones anteriores:

ns()=1sn(),Carolina del Norte()=1cn(),Dakota del Norte()=1dn().{\displaystyle \operatorname {ns} (u)={\frac {1}{\operatorname {sn} (u)}},\qquad \operatorname {nc} (u)={\frac {1}{\operatorname {cn} (u)}},\qquad \operatorname {nd} (u)={\frac {1}{\operatorname {dn} (u)}}.}

De manera similar, las razones de las tres funciones primarias corresponden a la primera letra del numerador seguida de la primera letra del denominador:

Carolina del Sur()=sn()cn(),Dakota del Sur()=sn()dn(),corriente continua()=dn()cn(),ds()=dn()sn(),cs()=cn()sn(),cd()=cn()dn().{\displaystyle {\begin{aligned}\operatorname {sc} (u)&={\frac {\operatorname {sn} (u)}{\operatorname {cn} (u)}},&\operatorname {sd} (u)&={\frac {\operatorname {sn} (u)}{\operatorname {dn} (u)}},&\operatorname {dc} (u)&={\frac {\operatorname {dn} (u)}{\operatorname {cn} (u)}},\\[1ex]\operatorname {ds} (u)&={\frac {\operatorname {dn} (u)}{\operatorname {sn} (u)}},&\operatorname {cs} (u)&={\frac {\operatorname {cn} (u)}{\operatorname {sn} (u)}},&\operatorname {cd} (u)&={\frac {\operatorname {cn} (u)}{\operatorname {dn} (u)}}.\end{aligned}}} De forma más compacta, tenemos

pq()=pn()qn(){\displaystyle \operatorname {pq} (u)={\frac {\operatorname {pn} (u)}{\operatorname {qn} (u)}}}

donde p y q son cualquiera de las letras s, c, d.

Periodicidad, polos y residuos

Gráficas de la fase para las doce funciones elípticas de Jacobi pq(u,m) como función de argumento complejo u, con polos y ceros indicados. Las gráficas abarcan un ciclo completo en las direcciones real e imaginaria, con la porción coloreada indicando la fase según la rueda de colores en la parte inferior derecha (que reemplaza la función dd trivial). Las regiones con valor absoluto inferior a 1/3 están coloreadas de negro, indicando aproximadamente la ubicación de un cero, mientras que las regiones con valor absoluto superior a 3 están coloreadas de blanco, indicando aproximadamente la posición de un polo. Todas las gráficas usan m  =  2/3 con K  = K ( m ), K = K (1 − m ), K ( ) siendo la integral elíptica completa de primera especie. Las flechas en los polos apuntan en la dirección de fase cero. Las flechas derecha e izquierda implican residuos reales positivos y negativos respectivamente. Las flechas arriba y abajo implican residuos imaginarios positivos y negativos respectivamente.     

En el plano complejo del argumento u , las funciones elípticas de Jacobi forman un patrón repetitivo de polos (y ceros). Los residuos de los polos tienen el mismo valor absoluto, diferenciándose únicamente en el signo. Cada función pq( u , m ) tiene una "función inversa" (en sentido multiplicativo) qp( u , m ) en la que se intercambian las posiciones de los polos y los ceros. Los periodos de repetición suelen ser diferentes en las direcciones real e imaginaria, de ahí el uso del término "doblemente periódicas" para describirlas.

Para la amplitud de Jacobi y la función épsilon de Jacobi: soy(+2K,metro)=soy(,metro)+π,{\displaystyle \operatorname {am} (u+2K,m)=\operatorname {am} (u,m)+\pi ,}soy(+4iK,metro)=soy(,metro),{\displaystyle \operatorname {am} (u+4iK',m)=\operatorname {am} (u,m),}mi(+2K,metro)=mi(,metro)+2mi,{\displaystyle {\mathcal {E}}(u+2K,m)={\mathcal {E}}(u,m)+2E,}mi(+2iK,metro)=mi(,metro)+2imiKKπiK{\displaystyle {\mathcal {E}}(u+2iK',m)={\mathcal {E}}(u,m)+2iE{\frac {K'}{K}}-{\frac {\pi i}{K}}} dóndemi(metro){\displaystyle E(m)}es la integral elíptica completa de segundo tipo con parámetrometro{\displaystyle m}.

La doble periodicidad de las funciones elípticas de Jacobi puede expresarse como:

pq(+2αK(metro)+2iβK(1metro),metro)=(1)γpq(,metro){\displaystyle \operatorname {pq} (u+2\alpha K(m)+2i\beta K(1-m)\,,\,m)=(-1)^{\gamma }\operatorname {pq} (u,m)}

donde α y β son pares cualesquiera de enteros. K ( ) es la integral elíptica completa de primera especie, también conocida como cuarto de período . La potencia de la unidad negativa ( γ ) se da en la siguiente tabla:

Cuando el factor (−1) γ es igual a −1, la ecuación expresa cuasiperiodicidad. Cuando es igual a la unidad, expresa periodicidad completa. Se puede ver, por ejemplo, que para las entradas que contienen solo α cuando α es par, la periodicidad completa se expresa mediante la ecuación anterior, y la función tiene periodos completos de 4 K ( m ) y 2 iK (1  m ). Del mismo modo, las funciones con entradas que contienen solo β tienen periodos completos de 2K(m) y 4 iK (1 − m ), mientras que aquellas con α + β tienen periodos completos de 4 K ( m ) y 4 iK (1 − m ).     

En el diagrama de la derecha, que representa una unidad repetitiva para cada función, indicando la fase junto con la ubicación de los polos y ceros, se pueden observar varias regularidades: La inversa de cada función se encuentra frente a la diagonal y tiene la misma celda unitaria, con los polos y ceros intercambiados. La disposición de polos y ceros en el rectángulo auxiliar formado por (0,0), ( K ,0), (0, K ) y ( K , K ) coincide con la descripción de la ubicación de polos y ceros descrita en la introducción. Además, el tamaño de los óvalos blancos que indican los polos es una medida aproximada del valor absoluto del residuo para ese polo. Los residuos de los polos más cercanos al origen en la figura (es decir, en el rectángulo auxiliar) se enumeran en la siguiente tabla:

Cuando corresponda, los polos desplazados hacia arriba en 2 K o hacia la derecha en 2 K tienen el mismo valor pero con signos invertidos, mientras que los opuestos en diagonal tienen el mismo valor. Cabe destacar que los polos y ceros en los bordes izquierdo e inferior se consideran parte de la celda unitaria, mientras que los de los bordes superior y derecho no.

La información sobre los polos puede utilizarse, de hecho, para caracterizar las funciones elípticas de Jacobi: [ 14 ]

La funciónsn(,metro){\displaystyle u\mapsto \operatorname {sn} (u,m)}es la única función elíptica que tiene polos simples en2rK+(2s+1)iK{\displaystyle 2rK+(2s+1)iK'}(conr,sZ{\displaystyle r,s\in \mathbb {Z} }) con residuos(1)r/metro{\displaystyle (-1)^{r}/{\sqrt {m}}}tomando el valor0{\displaystyle 0}en0{\displaystyle 0}.

La funcióncn(,metro){\displaystyle u\mapsto \operatorname {cn} (u,m)}es la única función elíptica que tiene polos simples en2rK+(2s+1)iK{\displaystyle 2rK+(2s+1)iK'}(conr,sZ{\displaystyle r,s\in \mathbb {Z} }) con residuos(1)r+s1i/metro{\displaystyle (-1)^{r+s-1}i/{\sqrt {m}}}tomando el valor1{\displaystyle 1}en0{\displaystyle 0}.

La funcióndn(,metro){\displaystyle u\mapsto \operatorname {dn} (u,m)}es la única función elíptica que tiene polos simples en2rK+(2s+1)iK{\displaystyle 2rK+(2s+1)iK'}(conr,sZ{\displaystyle r,s\in \mathbb {Z} }) con residuos(1)s1i{\displaystyle (-1)^{s-1}i}tomando el valor1{\displaystyle 1}en0{\displaystyle 0}.

Valores especiales

Configuraciónmetro=1{\displaystyle m=-1}proporciona las funciones elípticas de la lemniscatasl{\displaystyle \operatorname {sl} }ycl{\displaystyle \operatorname {cl} }:

sl=sn(,1),cl=cd(,1)=cn(,1)dn(,1).{\displaystyle \operatorname {sl} u=\operatorname {sn} (u,-1),\quad \operatorname {cl} u=\operatorname {cd} (u,-1)={\frac {\operatorname {cn} (u,-1)}{\operatorname {dn} (u,-1)}}.}

Cuandometro=0{\displaystyle m=0}ometro=1{\displaystyle m=1}Las funciones elípticas de Jacobi se reducen a funciones no elípticas:

Para la amplitud de Jacobi,soy(,0)={\displaystyle \operatorname {am} (u,0)=u}ysoy(,1)=buen día{\displaystyle \operatorname {am} (u,1)=\operatorname {gd} u}dóndebuen día{\displaystyle \operatorname {gd} }es la función de Gudermann .

En general, si ninguno de p,q es d, entoncespq(,1)=pq(buen día(),0){\displaystyle \operatorname {pq} (u,1)=\operatorname {pq} (\operatorname {gd} (u),0)}.

Identidades

Fórmula del semiángulo

sn(2,metro)=±1cn(,metro)1+dn(,metro){\displaystyle \operatorname {sn} \left({\frac {u}{2}},m\right)=\pm {\sqrt {\frac {1-\operatorname {cn} (u,m)}{1+\operatorname {dn} (u,m)}}}}cn(2,metro)=±cn(,metro)+dn(,metro)1+dn(,metro){\displaystyle \operatorname {cn} \left({\frac {u}{2}},m\right)=\pm {\sqrt {\frac {\operatorname {cn} (u,m)+\operatorname {dn} (u,m)}{1+\operatorname {dn} (u,m)}}}}dn(2,metro)=±metro+dn(,metro)+metrocn(,metro)1+dn(,metro){\displaystyle \operatorname {dn} \left({\frac {u}{2}},m\right)=\pm {\sqrt {\frac {m'+\operatorname {dn} (u,m)+m\operatorname {cn} (u,m)}{1+\operatorname {dn} (u,m)}}}}

Fórmulas K

Fórmula de media K

sn[12K(k);k]=21+k+1k{\displaystyle \operatorname {sn} \left[{\tfrac {1}{2}}K(k);k\right]={\frac {\sqrt {2}}{{\sqrt {1+k}}+{\sqrt {1-k}}}}}

cn[12K(k);k]=21k241+k+1k{\displaystyle \operatorname {cn} \left[{\tfrac {1}{2}}K(k);k\right]={\frac {{\sqrt {2}}\,{\sqrt[{4}]{1-k^{2}}}}{{\sqrt {1+k}}+{\sqrt {1-k}}}}}

dn[12K(k);k]=1k24{\displaystyle \operatorname {dn} \left[{\tfrac {1}{2}}K(k);k\right]={\sqrt[{4}]{1-k^{2}}}}

Tercera fórmula K

sn[13K(incógnita3incógnita6+1+1);incógnita3incógnita6+1+1]=2incógnita4incógnita2+1incógnita2+2+incógnita2+112incógnita4incógnita2+1incógnita2+2+incógnita2+1+1{\displaystyle \operatorname {sn} \left[{\frac {1}{3}}K\left({\frac {x^{3}}{{\sqrt {x^{6}+1}}+1}}\right);{\frac {x^{3}}{{\sqrt {x^{6}+1}}+1}}\right]={\frac {{\sqrt {2{\sqrt {x^{4}-x^{2}+1}}-x^{2}+2}}+{\sqrt {x^{2}+1}}-1}{{\sqrt {2{\sqrt {x^{4}-x^{2}+1}}-x^{2}+2}}+{\sqrt {x^{2}+1}}+1}}}

Para obtener x 3 , tomamos la tangente del doble de la arcotangente del módulo.

Además, esta ecuación conduce al valor sn del tercio de K :

k2s42k2s3+2s1=0{\displaystyle k^{2}s^{4}-2k^{2}s^{3}+2s-1=0}

s=sn[13K(k);k]{\displaystyle s=\operatorname {sn} \left[{\tfrac {1}{3}}K(k);k\right]}

Estas ecuaciones conducen a los demás valores de las funciones de Jacobi:

cn[23K(k);k]=1sn[13K(k);k]{\displaystyle \operatorname {cn} \left[{\tfrac {2}{3}}K(k);k\right]=1-\operatorname {sn} \left[{\tfrac {1}{3}}K(k);k\right]}

dn[23K(k);k]=1/sn[13K(k);k]1{\displaystyle \operatorname {dn} \left[{\tfrac {2}{3}}K(k);k\right]=1/\operatorname {sn} \left[{\tfrac {1}{3}}K(k);k\right]-1}

Quinta fórmula K

La siguiente ecuación tiene la siguiente solución:

4k2incógnita6+8k2incógnita5+2(1k2)2incógnita(1k2)2=0{\displaystyle 4k^{2}x^{6}+8k^{2}x^{5}+2(1-k^{2})^{2}x-(1-k^{2})^{2}=0}

incógnita=1212k2sn[25K(k);k]2sn[45K(k);k]2=sn[45K(k);k]2sn[25K(k);k]22sn[25K(k);k]sn[45K(k);k]{\displaystyle x={\frac {1}{2}}-{\frac {1}{2}}k^{2}\operatorname {sn} \left[{\tfrac {2}{5}}K(k);k\right]^{2}\operatorname {sn} \left[{\tfrac {4}{5}}K(k);k\right]^{2}={\frac {\operatorname {sn} \left[{\frac {4}{5}}K(k);k\right]^{2}-\operatorname {sn} \left[{\frac {2}{5}}K(k);k\right]^{2}}{2\operatorname {sn} \left[{\frac {2}{5}}K(k);k\right]\operatorname {sn} \left[{\frac {4}{5}}K(k);k\right]}}}

Para obtener los valores sn, introducimos la solución x en las siguientes expresiones:

sn[25K(k);k]=(1+k2)1/22(1incógnitaincógnita2)(incógnita2+1incógnitaincógnita2+1){\displaystyle \operatorname {sn} \left[{\tfrac {2}{5}}K(k);k\right]=(1+k^{2})^{-1/2}{\sqrt {2(1-x-x^{2})(x^{2}+1-x{\sqrt {x^{2}+1}})}}}

sn[45K(k);k]=(1+k2)1/22(1incógnitaincógnita2)(incógnita2+1+incógnitaincógnita2+1){\displaystyle \operatorname {sn} \left[{\tfrac {4}{5}}K(k);k\right]=(1+k^{2})^{-1/2}{\sqrt {2(1-x-x^{2})(x^{2}+1+x{\sqrt {x^{2}+1}})}}}

Relaciones entre cuadrados de las funciones

Las relaciones entre los cuadrados de las funciones se pueden derivar de dos relaciones básicas (Argumentos ( u , m ) omitidos): cn2+sn2=1{\displaystyle \operatorname {cn} ^{2}+\operatorname {sn} ^{2}=1}cn2+metrosn2=dn2{\displaystyle \operatorname {cn} ^{2}+m'\operatorname {sn} ^{2}=\operatorname {dn} ^{2}} donde m  +  m'  =  1. Multiplicando por cualquier función de la forma nq se obtienen ecuaciones más generales:

cq2+cuadrado2=nq2{\displaystyle \operatorname {cq} ^{2}+\operatorname {sq} ^{2}=\operatorname {nq} ^{2}}cq2+metrocuadrado2=dq2{\displaystyle \operatorname {cq} ^{2}{}+m'\operatorname {sq} ^{2}=\operatorname {dq} ^{2}}

Con q  = d , estas corresponden trigonométricamente a las ecuaciones del círculo unitario ( incógnita2+y2=r2{\displaystyle x^{2}+y^{2}=r^{2}}) y la elipse unitaria (incógnita2+metroy2=1{\displaystyle x^{2}{}+m'y^{2}=1}), con x  = cd , y = sd y r = nd . Usando la regla de la multiplicación, se pueden derivar otras relaciones. Por ejemplo:     

dn2+metro=metrocn2=metrosn2metro{\displaystyle -\operatorname {dn} ^{2}{}+m'=-m\operatorname {cn} ^{2}=m\operatorname {sn} ^{2}-m}

metroDakota del Norte2+metro=metrometroDakota del Sur2=metrocd2metro{\displaystyle -m'\operatorname {nd} ^{2}{}+m'=-mm'\operatorname {sd} ^{2}=m\operatorname {cd} ^{2}-m}

metroCarolina del Sur2+metro=metroCarolina del Norte2=corriente continua2metro{\displaystyle m'\operatorname {sc} ^{2}{}+m'=m'\operatorname {nc} ^{2}=\operatorname {dc} ^{2}-m}

cs2+metro=ds2=ns2metro{\displaystyle \operatorname {cs} ^{2}{}+m'=\operatorname {ds} ^{2}=\operatorname {ns} ^{2}-m}

Teoremas de la suma

Las funciones satisfacen las dos relaciones cuadráticas (dependencia de m suprimida). cn2()+sn2()=1,{\displaystyle \operatorname {cn} ^{2}(u)+\operatorname {sn} ^{2}(u)=1,\,}

dn2()+metrosn2()=1.{\displaystyle \operatorname {dn} ^{2}(u)+m\operatorname {sn} ^{2}(u)=1.\,}

De esto vemos que (cn, sn, dn) parametriza una curva elíptica que es la intersección de las dos cuádricas definidas por las dos ecuaciones anteriores. Ahora podemos definir una ley de grupo para puntos en esta curva mediante las fórmulas de adición para las funciones de Jacobi [ 3 ].

cn(incógnita+y)=cn(incógnita)cn(y)sn(incógnita)sn(y)dn(incógnita)dn(y)1metrosn2(incógnita)sn2(y),sn(incógnita+y)=sn(incógnita)cn(y)dn(y)+sn(y)cn(incógnita)dn(incógnita)1metrosn2(incógnita)sn2(y),dn(incógnita+y)=dn(incógnita)dn(y)metrosn(incógnita)sn(y)cn(incógnita)cn(y)1metrosn2(incógnita)sn2(y).{\displaystyle {\begin{aligned}\operatorname {cn} (x+y)&={\operatorname {cn} (x)\operatorname {cn} (y)-\operatorname {sn} (x)\operatorname {sn} (y)\operatorname {dn} (x)\operatorname {dn} (y) \over {1-m\operatorname {sn} ^{2}(x)\operatorname {sn} ^{2}(y)}},\\[8pt]\operatorname {sn} (x+y)&={\operatorname {sn} (x)\operatorname {cn} (y)\operatorname {dn} (y)+\operatorname {sn} (y)\operatorname {cn} (x)\operatorname {dn} (x) \over {1-m\operatorname {sn} ^{2}(x)\operatorname {sn} ^{2}(y)}},\\[8pt]\operatorname {dn} (x+y)&={\operatorname {dn} (x)\operatorname {dn} (y)-m\operatorname {sn} (x)\operatorname {sn} (y)\operatorname {cn} (x)\operatorname {cn} (y) \over {1-m\operatorname {sn} ^{2}(x)\operatorname {sn} ^{2}(y)}}.\end{aligned}}}

Las funciones épsilon y zn de Jacobi satisfacen un teorema de cuasi-adición: mi(incógnita+y,metro)=mi(incógnita,metro)+mi(y,metro)metrosn(incógnita,metro)sn(y,metro)sn(incógnita+y,metro),zinc(incógnita+y,metro)=zinc(incógnita,metro)+zinc(y,metro)metrosn(incógnita,metro)sn(y,metro)sn(incógnita+y,metro).{\displaystyle {\begin{aligned}{\mathcal {E}}(x+y,m)&={\mathcal {E}}(x,m)+{\mathcal {E}}(y,m)-m\operatorname {sn} (x,m)\operatorname {sn} (y,m)\operatorname {sn} (x+y,m),\\\operatorname {zn} (x+y,m)&=\operatorname {zn} (x,m)+\operatorname {zn} (y,m)-m\operatorname {sn} (x,m)\operatorname {sn} (y,m)\operatorname {sn} (x+y,m).\end{aligned}}}

Las fórmulas del ángulo doble se pueden derivar fácilmente de las ecuaciones anteriores haciendo x  = y . [ 3 ] Las fórmulas del ángulo medio [ 11 ] [ 3 ] son ​​todas de la forma: 

pq(12,metro)2=Fpag/Fq{\displaystyle \operatorname {pq} ({\tfrac {1}{2}}u,m)^{2}=f_{\mathrm {p} }/f_{\mathrm {q} }}

dónde: Fdo=cn(,metro)+dn(,metro){\displaystyle f_{\mathrm {c} }=\operatorname {cn} (u,m)+\operatorname {dn} (u,m)}Fs=1cn(,metro){\displaystyle f_{\mathrm {s} }=1-\operatorname {cn} (u,m)}Fnorte=1+dn(,metro){\displaystyle f_{\mathrm {n} }=1+\operatorname {dn} (u,m)}Fd=(1+dn(,metro))metro(1cn(,metro)){\displaystyle f_{\mathrm {d} }=(1+\operatorname {dn} (u,m))-m(1-\operatorname {cn} (u,m))}

Funciones elípticas de Jacobi como soluciones de ecuaciones diferenciales ordinarias no lineales

Derivadas con respecto a la primera variable

Las derivadas de las tres funciones elípticas básicas de Jacobi (con respecto a la primera variable, conmetro{\displaystyle m}fijos) son: ddzsn(z)=cn(z)dn(z),{\displaystyle {\frac {\mathrm {d} }{\mathrm {d} z}}\operatorname {sn} (z)=\operatorname {cn} (z)\operatorname {dn} (z),}ddzcn(z)=sn(z)dn(z),{\displaystyle {\frac {\mathrm {d} }{\mathrm {d} z}}\operatorname {cn} (z)=-\operatorname {sn} (z)\operatorname {dn} (z),}ddzdn(z)=metrosn(z)cn(z).{\displaystyle {\frac {\mathrm {d} }{\mathrm {d} z}}\operatorname {dn} (z)=-m\operatorname {sn} (z)\operatorname {cn} (z).}

Estos se pueden utilizar para derivar las derivadas de todas las demás funciones como se muestra en la tabla siguiente (argumentos (u,m) omitidos):

También ddzmi(z)=dn(z)2.{\displaystyle {\frac {\mathrm {d} }{\mathrm {d} z}}{\mathcal {E}}(z)=\operatorname {dn} (z)^{2}.}

Con los teoremas de adición anteriores y para un m dado con 0  < m < 1, las funciones principales son, por lo tanto, soluciones de las siguientes ecuaciones diferenciales ordinarias no lineales :   

  • soy(incógnita){\displaystyle \operatorname {am} (x)}Resuelve las ecuaciones diferencialesd2ydincógnita2+metropecado(y)porque(y)=0{\displaystyle {\frac {\mathrm {d} ^{2}y}{\mathrm {d} x^{2}}}+m\sin(y)\cos(y)=0}y

(dydincógnita)2=1metropecado(y)2{\displaystyle \left({\frac {\mathrm {d} y}{\mathrm {d} x}}\right)^{2}=1-m\sin(y)^{2}}(paraincógnita{\displaystyle x}(no en una rama cortada)

  • sn(incógnita){\displaystyle \operatorname {sn} (x)}Resuelve las ecuaciones diferencialesd2ydincógnita2+(1+metro)y2metroy3=0{\displaystyle {\frac {\mathrm {d} ^{2}y}{\mathrm {d} x^{2}}}+(1+m)y-2my^{3}=0}y(dydincógnita)2=(1y2)(1metroy2){\displaystyle \left({\frac {\mathrm {d} y}{\mathrm {d} x}}\right)^{2}=(1-y^{2})(1-my^{2})}
  • cn(incógnita){\displaystyle \operatorname {cn} (x)}Resuelve las ecuaciones diferencialesd2ydincógnita2+(12metro)y+2metroy3=0{\displaystyle {\frac {\mathrm {d} ^{2}y}{\mathrm {d} x^{2}}}+(1-2m)y+2my^{3}=0}y(dydincógnita)2=(1y2)(1metro+metroy2){\displaystyle \left({\frac {\mathrm {d} y}{\mathrm {d} x}}\right)^{2}=(1-y^{2})(1-m+my^{2})}
  • dn(incógnita){\displaystyle \operatorname {dn} (x)}Resuelve las ecuaciones diferencialesd2ydincógnita2(2metro)y+2y3=0{\displaystyle {\frac {\mathrm {d} ^{2}y}{\mathrm {d} x^{2}}}-(2-m)y+2y^{3}=0}y(dydincógnita)2=(y21)(1metroy2){\displaystyle \left({\frac {\mathrm {d} y}{\mathrm {d} x}}\right)^{2}=(y^{2}-1)(1-m-y^{2})}

La función que resuelve exactamente la ecuación diferencial del péndulo , d2θdt2+dopecadoθ=0,{\displaystyle {\frac {\mathrm {d} ^{2}\theta }{\mathrm {d} t^{2}}}+c\sin \theta =0,} con ángulo inicialθ0{\displaystyle \theta _{0}}y velocidad angular inicial cero es

θ=2arcoseno(metrocd(dot,metro))=2soy(1+metro2(dot+K),4metro(1+metro)2)2soy(1+metro2(dotK),4metro(1+metro)2)π{\displaystyle {\begin{aligned}\theta &=2\arcsin({\sqrt {m}}\operatorname {cd} ({\sqrt {c}}t,m))\\&=2\operatorname {am} \left({\frac {1+{\sqrt {m}}}{2}}({\sqrt {c}}t+K),{\frac {4{\sqrt {m}}}{(1+{\sqrt {m}})^{2}}}\right)-2\operatorname {am} \left({\frac {1+{\sqrt {m}}}{2}}({\sqrt {c}}t-K),{\frac {4{\sqrt {m}}}{(1+{\sqrt {m}})^{2}}}\right)-\pi \end{aligned}}} dóndemetro=pecado(θ0/2)2{\displaystyle m=\sin(\theta _{0}/2)^{2}},do>0{\displaystyle c>0}ytR{\displaystyle t\in \mathbb {R} }.

Derivadas con respecto a la segunda variable

Con el primer argumentoz{\displaystyle z}fijas, las derivadas con respecto a la segunda variablemetro{\displaystyle m}son los siguientes: [ 15 ]

ddmetrosn(z)=dn(z)cn(z)((1metro)zmi(z)+metrocd(z)sn(z))2metro(1metro),ddmetrocn(z)=sn(z)dn(z)((metro1)z+mi(z)metrosn(z)cd(z))2metro(1metro),ddmetrodn(z)=sn(z)cn(z)((metro1)z+mi(z)dn(z)Carolina del Sur(z))2(1metro),ddmetromi(z)=cn(z)(sn(z)dn(z)cn(z)mi(z))2(1metro)z2sn(z)2.{\displaystyle {\begin{aligned}{\frac {\mathrm {d} }{\mathrm {d} m}}\operatorname {sn} (z)&={\frac {\operatorname {dn} (z)\operatorname {cn} (z)((1-m)z-{\mathcal {E}}(z)+m\operatorname {cd} (z)\operatorname {sn} (z))}{2m(1-m)}},\\{\frac {\mathrm {d} }{\mathrm {d} m}}\operatorname {cn} (z)&={\frac {\operatorname {sn} (z)\operatorname {dn} (z)((m-1)z+{\mathcal {E}}(z)-m\operatorname {sn} (z)\operatorname {cd} (z))}{2m(1-m)}},\\{\frac {\mathrm {d} }{\mathrm {d} m}}\operatorname {dn} (z)&={\frac {\operatorname {sn} (z)\operatorname {cn} (z)((m-1)z+{\mathcal {E}}(z)-\operatorname {dn} (z)\operatorname {sc} (z))}{2(1-m)}},\\{\frac {\mathrm {d} }{\mathrm {d} m}}{\mathcal {E}}(z)&={\frac {\operatorname {cn} (z)(\operatorname {sn} (z)\operatorname {dn} (z)-\operatorname {cn} (z){\mathcal {E}}(z))}{2(1-m)}}-{\frac {z}{2}}\operatorname {sn} (z)^{2}.\end{aligned}}}

Expansión en términos del nomo

Que el nombre seaq=exp(πK(metro)/K(metro))=miiπτ{\displaystyle q=\exp(-\pi K'(m)/K(m))=e^{i\pi \tau }},Soy(τ)>0{\displaystyle \operatorname {Im} (\tau )>0},metro=k2{\displaystyle m=k^{2}}y dejarv=π/(2K(metro)){\displaystyle v=\pi u/(2K(m))}Entonces, las funciones tienen expansiones como series de Lambert .

soy(,metro)=π2K(metro)+2norte=1qnortenorte(1+q2norte)pecado(2nortev),{\displaystyle \operatorname {am} (u,m)={\frac {\pi u}{2K(m)}}+2\sum _{n=1}^{\infty }{\frac {q^{n}}{n(1+q^{2n})}}\sin(2nv),}

sn(,metro)=2πkK(metro)norte=0qnorte+1/21q2norte+1pecado((2norte+1)v),{\displaystyle \operatorname {sn} (u,m)={\frac {2\pi }{kK(m)}}\sum _{n=0}^{\infty }{\frac {q^{n+1/2}}{1-q^{2n+1}}}\sin((2n+1)v),}

cn(,metro)=2πkK(metro)norte=0qnorte+1/21+q2norte+1porque((2norte+1)v),{\displaystyle \operatorname {cn} (u,m)={\frac {2\pi }{kK(m)}}\sum _{n=0}^{\infty }{\frac {q^{n+1/2}}{1+q^{2n+1}}}\cos((2n+1)v),}

dn(,metro)=π2K(metro)+2πK(metro)norte=1qnorte1+q2norteporque(2nortev),{\displaystyle \operatorname {dn} (u,m)={\frac {\pi }{2K(m)}}+{\frac {2\pi }{K(m)}}\sum _{n=1}^{\infty }{\frac {q^{n}}{1+q^{2n}}}\cos(2nv),}

zinc(,metro)=2πK(metro)norte=1qnorte1q2nortepecado(2nortev){\displaystyle \operatorname {zn} (u,m)={\frac {2\pi }{K(m)}}\sum _{n=1}^{\infty }{\frac {q^{n}}{1-q^{2n}}}\sin(2nv)}

cuando|Soy(/K)|<Soy(iK/K).{\displaystyle \left|\operatorname {Im} (u/K)\right|<\operatorname {Im} (iK'/K).}

Schett ha publicado expansiones en series de potencias bivariadas . [ 16 ]

Computación rápida

Las razones de la función theta proporcionan una forma eficiente de calcular las funciones elípticas de Jacobi. Existe un método alternativo, basado en la media aritmético-geométrica y las transformaciones de Landen : [ 6 ]

Inicializar a0=1,b0=1metro{\displaystyle a_{0}=1,\,b_{0}={\sqrt {1-m}}} dónde0<metro<1{\displaystyle 0<m<1}. Definir anorte=anorte1+bnorte12,bnorte=anorte1bnorte1,donorte=anorte1bnorte12{\displaystyle a_{n}={\frac {a_{n-1}+b_{n-1}}{2}},\,b_{n}={\sqrt {a_{n-1}b_{n-1}}},\,c_{n}={\frac {a_{n-1}-b_{n-1}}{2}}} dóndenorte1{\displaystyle n\geq 1}. Luego define φnorte=2norteanorte{\displaystyle \varphi _{N}=2^{N}a_{N}u} paraR{\displaystyle u\in \mathbb {R} }y un fijonortenorte{\displaystyle N\in \mathbb {N} }. Si φnorte1=12(φnorte+arcoseno(donorteanortepecadoφnorte)){\displaystyle \varphi _{n-1}={\frac {1}{2}}\left(\varphi _{n}+\arcsin \left({\frac {c_{n}}{a_{n}}}\sin \varphi _{n}\right)\right)} paranorte1{\displaystyle n\geq 1}, entonces soy(,metro)=φ0,zinc(,metro)=norte=1nortedonortepecadoφnorte{\displaystyle \operatorname {am} (u,m)=\varphi _{0},\quad \operatorname {zn} (u,m)=\sum _{n=1}^{N}c_{n}\sin \varphi _{n}} comonorte{\displaystyle N\to \infty }Esto es notable por su rápida convergencia. Entonces es trivial calcular todas las funciones elípticas de Jacobi a partir de la amplitud de Jacobi.soy{\displaystyle \operatorname {am} }en la línea real. [ nota 2 ]

En combinación con los teoremas de adición para funciones elípticas (que son válidos para números complejos en general) y las transformaciones de Jacobi, el método de cálculo descrito anteriormente puede utilizarse para calcular todas las funciones elípticas de Jacobi en todo el plano complejo.

Otro método de cálculo rápido de las funciones elípticas de Jacobi mediante la media aritmético-geométrica, evitando el cálculo de la amplitud de Jacobi, se debe a Herbert E. Salzer: [ 17 ]

Dejar 0metro1,0K(metro),a0=1,b0=1metro,{\displaystyle 0\leq m\leq 1,\,0\leq u\leq K(m),\,a_{0}=1,\,b_{0}={\sqrt {1-m}},}anorte+1=anorte+bnorte2,bnorte+1=anortebnorte,donorte+1=anortebnorte2.{\displaystyle a_{n+1}={\frac {a_{n}+b_{n}}{2}},\,b_{n+1}={\sqrt {a_{n}b_{n}}},\,c_{n+1}={\frac {a_{n}-b_{n}}{2}}.} Colocar ynorte=anortepecado(anorte)ynorte1=ynorte+anortedonorteynorteynorte2=ynorte1+anorte1donorte1ynorte1=y0=y1+metro4y1.{\displaystyle {\begin{aligned}y_{N}&={\frac {a_{N}}{\sin(a_{N}u)}}\\y_{N-1}&=y_{N}+{\frac {a_{N}c_{N}}{y_{N}}}\\y_{N-2}&=y_{N-1}+{\frac {a_{N-1}c_{N-1}}{y_{N-1}}}\\\vdots &=\vdots \\y_{0}&=y_{1}+{\frac {m}{4y_{1}}}.\end{aligned}}} Entonces sn(,metro)=1y0cn(,metro)=11y02dn(,metro)=1metroy02{\displaystyle {\begin{aligned}\operatorname {sn} (u,m)&={\frac {1}{y_{0}}}\\\operatorname {cn} (u,m)&={\sqrt {1-{\frac {1}{y_{0}^{2}}}}}\\\operatorname {dn} (u,m)&={\sqrt {1-{\frac {m}{y_{0}^{2}}}}}\end{aligned}}} comonorte{\displaystyle N\to \infty }.

Sin embargo, a continuación se muestra otro método para un cálculo rápido y de convergencia rápida de la función seno elíptico de Jacobi que se encuentra en la literatura. [ 18 ]

Dejar:

a0=b0=11metro1+1metroa1=a01+b0b1=11b021+1b02==anorte=anorte11+bnorte1bnorte=11bnorte121+1bnorte12{\displaystyle {\begin{aligned}&a_{0}=u&b_{0}={\frac {1-{\sqrt {1-m}}}{1+{\sqrt {1-m}}}}\\&a_{1}={\frac {a_{0}}{1+b_{0}}}&b_{1}={\frac {1-{\sqrt {1-b_{0}^{2}}}}{1+{\sqrt {1-b_{0}^{2}}}}}\\&\vdots =\vdots &\vdots =\vdots \\&a_{n}={\frac {a_{n-1}}{1+b_{n-1}}}&b_{n}={\frac {1-{\sqrt {1-b_{n-1}^{2}}}}{1+{\sqrt {1-b_{n-1}^{2}}}}}\\\end{aligned}}}

Luego, configure:

ynorte+1=pecado(anorte)ynorte=ynorte+1(1+bnorte)1+ynorte+12bnorte=y0=y1(1+b0)1+y12b0{\displaystyle {\begin{aligned}y_{n+1}&=\sin(a_{n})\\y_{n}&={\frac {y_{n+1}(1+b_{n})}{1+y_{n+1}^{2}b_{n}}}\\\vdots &=\vdots \\y_{0}&={\frac {y_{1}(1+b_{0})}{1+y_{1}^{2}b_{0}}}\\\end{aligned}}}

Entonces:

sn(,metro)=y0 como norte{\displaystyle \operatorname {sn} (u,m)=y_{0}{\text{ as }}n\rightarrow \infty }.

Aproximación en términos de funciones hiperbólicas

Las funciones elípticas de Jacobi pueden expandirse en términos de funciones hiperbólicas. Cuandometro{\displaystyle m}está cerca de la unidad, de tal manera quemetro2{\displaystyle m'^{2}}y poderes superiores de metro{\displaystyle m'}puede ser ignorado, tenemos: [ 19 ] [ 20 ]

  • sn( u ):sn(,metro)tanh()+14metro(sinh()aporrear())sech2().{\displaystyle \operatorname {sn} (u,m)\approx \tanh(u)+{\frac {1}{4}}m'(\sinh(u)\cosh(u)-u)\operatorname {sech} ^{2}(u).}
  • cn( u ):cn(,metro)sech()14metro(sinh()aporrear())tanh()sech().{\displaystyle \operatorname {cn} (u,m)\approx \operatorname {sech} (u)-{\frac {1}{4}}m'(\sinh(u)\cosh(u)-u)\tanh(u)\operatorname {sech} (u).}
  • dn( u ):dn(,metro)sech()+14metro(sinh()aporrear()+)tanh()sech().{\displaystyle \operatorname {dn} (u,m)\approx \operatorname {sech} (u)+{\frac {1}{4}}m'(\sinh(u)\cosh(u)+u)\tanh(u)\operatorname {sech} (u).}

Para la amplitud de Jacobi, soy(,metro)buen día()+14metro(sinh()aporrear())sech().{\displaystyle \operatorname {am} (u,m)\approx \operatorname {gd} (u)+{\frac {1}{4}}m'(\sinh(u)\cosh(u)-u)\operatorname {sech} (u).}

fracciones continuas

Suponiendo que los números son realesa,pag{\displaystyle a,p}con0<a<pag{\displaystyle 0<a<p}y el nombreq=miπiτ{\displaystyle q=e^{\pi i\tau }},Soy(τ)>0{\displaystyle \operatorname {Im} (\tau )>0}con módulo elípticok(τ)=1k(τ)2=(ϑ10(0;τ)/ϑ00(0;τ))2{\textstyle k(\tau )={\sqrt {1-k'(\tau )^{2}}}=(\vartheta _{10}(0;\tau )/\vartheta _{00}(0;\tau ))^{2}}. SiK[τ]=K(k(τ)){\displaystyle K[\tau ]=K(k(\tau ))}, dóndeK(incógnita)=π/22F1(1/2,1/2;1;incógnita2){\displaystyle K(x)=\pi /2\cdot {}_{2}F_{1}(1/2,1/2;1;x^{2})}es la integral elíptica completa de primera especie , entonces se cumple la siguiente expansión en fracción continua [ 21 ].dn((pag/2a)τK[pagτ2];k(pagτ2))k(pagτ2)=norte=qpag/2norte2+(pag/2a)nortenorte=(1)norteqpag/2norte2+(pag/2a)norte=1+21qa+qpaga1qpag+(qa+q2paga)(qa+pag+qpaga)1q3pag+qpag(qa+q3paga)(qa+2pag+qpaga)1q5pag+q2pag(qa+q4paga)(qa+3pag+qpaga)1q7pag+{\displaystyle {\begin{aligned}&{\frac {{\textrm {dn}}\left((p/2-a)\tau K\left[{\frac {p\tau }{2}}\right];k\left({\frac {p\tau }{2}}\right)\right)}{\sqrt {k'\left({\frac {p\tau }{2}}\right)}}}={\frac {\sum \limits _{n=-\infty }^{\infty }q^{p/2n^{2}+(p/2-a)n}}{\sum \limits _{n=-\infty }^{\infty }(-1)^{n}q^{p/2n^{2}+(p/2-a)n}}}\\[4pt]&=-1+{\cfrac {2}{1-{\cfrac {q^{a}+q^{p-a}}{1-q^{p}+{\cfrac {(q^{a}+q^{2p-a})(q^{a+p}+q^{p-a})}{1-q^{3p}+{\cfrac {q^{p}(q^{a}+q^{3p-a})(q^{a+2p}+q^{p-a})}{1-q^{5p}+{\cfrac {q^{2p}(q^{a}+q^{4p-a})(q^{a+3p}+q^{p-a})}{1-q^{7p}+\cdots }}}}}}}}}}\end{aligned}}} Fracciones continuas conocidas que involucransn(t),cn(t){\displaystyle {\textrm {sn}}(t),{\textrm {cn}}(t)}ydn(t){\displaystyle {\textrm {dn}}(t)}con módulo elípticok{\displaystyle k}son

Parazdo{\displaystyle z\in \mathbb {C} },|k|<1{\displaystyle |k|<1}: [ 22 ] pág. 374 0sn(t)mitzdt=112(1+k2)+z21223k232(1+k2)+z23425k252(1+k2)+z2{\displaystyle \int _{0}^{\infty }{\textrm {sn}}(t)e^{-tz}\,\mathrm {d} t={\frac {1}{1^{2}(1+k^{2})+z^{2}-{}}}\,{\frac {1\cdot 2^{2}\cdot 3k^{2}}{3^{2}(1+k^{2})+z^{2}-{}}}\,{\frac {3\cdot 4^{2}\cdot 5k^{2}}{5^{2}(1+k^{2})+z^{2}-{}}}\cdots }

Parazdo{0}{\displaystyle z\in \mathbb {C} \setminus \{0\}},|k|<1{\displaystyle |k|<1}: [ 22 ] pág. 375 0sn2(t)mitzdt=2z122(1+k2)+z22324k242(1+k2)+z24526k262(1+k2)+z2{\displaystyle \int _{0}^{\infty }{\textrm {sn}}^{2}(t)e^{-tz}\,\mathrm {d} t={\frac {2z^{-1}}{2^{2}(1+k^{2})+z^{2}-{}}}\,{\frac {2\cdot 3^{2}\cdot 4k^{2}}{4^{2}(1+k^{2})+z^{2}-{}}}\,{\frac {4\cdot 5^{2}\cdot 6k^{2}}{6^{2}(1+k^{2})+z^{2}-{}}}\cdots }

Parazdo{0}{\displaystyle z\in \mathbb {C} \setminus \{0\}},|k|<1{\displaystyle |k|<1}: [ 23 ] pág. 220 0cn(t)mitzdt=1z+12z+22k2z+32z+42k2z+52z+{\displaystyle \int _{0}^{\infty }{\textrm {cn}}(t)e^{-tz}\,\mathrm {d} t={\frac {1}{z+{}}}\,{\frac {1^{2}}{z+{}}}\,{\frac {2^{2}k^{2}}{z+{}}}\,{\frac {3^{2}}{z+{}}}\,{\frac {4^{2}k^{2}}{z+{}}}\,{\frac {5^{2}}{z+{}}}\cdots }

Parazdo{0}{\displaystyle z\in \mathbb {C} \setminus \{0\}},|k|<1{\displaystyle |k|<1}: [ 22 ] pág. 374 0dn(t)mitzdt=1z+12k2z+22z+32k2z+42z+52k2z+{\displaystyle \int _{0}^{\infty }{\textrm {dn}}(t)e^{-tz}\,\mathrm {d} t={\frac {1}{z+{}}}\,{\frac {1^{2}k^{2}}{z+{}}}\,{\frac {2^{2}}{z+{}}}\,{\frac {3^{2}k^{2}}{z+{}}}\,{\frac {4^{2}}{z+{}}}\,{\frac {5^{2}k^{2}}{z+{}}}\cdots }

Parazdo{\displaystyle z\in \mathbb {C} },|k|<1{\displaystyle |k|<1}: [ 22 ] pág. 375 0sn(t)cn(t)dn(t)mitzdt=1212(2k2)+z21223k4232(2k2)+z23425k4252(2k2)+z2{\displaystyle \int _{0}^{\infty }{\frac {{\textrm {sn}}(t){\textrm {cn}}(t)}{{\textrm {dn}}(t)}}e^{-tz}\,\mathrm {d} t={\frac {1}{2\cdot 1^{2}(2-k^{2})+z^{2}-{}}}\,{\frac {1\cdot 2^{2}\cdot 3k^{4}}{2\cdot 3^{2}(2-k^{2})+z^{2}-{}}}\,{\frac {3\cdot 4^{2}\cdot 5k^{4}}{2\cdot 5^{2}(2-k^{2})+z^{2}-{}}}\cdots }

Funciones inversas

Las inversas de las funciones elípticas de Jacobi se pueden definir de forma similar a las funciones trigonométricas inversas ; siincógnita=sn(ξ,metro){\displaystyle x=\operatorname {sn} (\xi ,m)},ξ=arcosn(incógnita,metro){\displaystyle \xi =\operatorname {arcsn} (x,m)}Se pueden representar como integrales elípticas, [ 24 ] [ 25 ] [ 26 ] y se han encontrado representaciones en series de potencias. [ 27 ] [ 3 ]

  • arcosn(incógnita,metro)=0incógnitadt(1t2)(1metrot2){\displaystyle \operatorname {arcsn} (x,m)=\int _{0}^{x}{\frac {\mathrm {d} t}{\sqrt {(1-t^{2})(1-mt^{2})}}}}
  • arccn(incógnita,metro)=incógnita1dt(1t2)(1metro+metrot2){\displaystyle \operatorname {arccn} (x,m)=\int _{x}^{1}{\frac {\mathrm {d} t}{\sqrt {(1-t^{2})(1-m+mt^{2})}}}}
  • arco(incógnita,metro)=incógnita1dt(1t2)(t2+metro1){\displaystyle \operatorname {arcdn} (x,m)=\int _{x}^{1}{\frac {\mathrm {d} t}{\sqrt {(1-t^{2})(t^{2}+m-1)}}}}

Proyección cartográfica

La proyección quincuncial de Peirce es una proyección cartográfica basada en funciones elípticas jacobianas.

Véase también

Notas

  1. SiR{\displaystyle u\in \mathbb {R} }ymetro{\displaystyle m}está restringido a[0,1]{\displaystyle [0,1]}, entoncesdn(|metro){\displaystyle \operatorname {dn} (u|m)}también se puede escribir como1metropecado2soy(|metro).{\displaystyle {\sqrt {1-m\sin ^{2}\operatorname {am} (u|m)}}.}
  2. Para eldn{\displaystyle \operatorname {dn} }función, dn(,metro)=cn(,metro)sn(K(metro),metro){\displaystyle \operatorname {dn} (u,m)={\frac {\operatorname {cn} (u,m)}{\operatorname {sn} (K(m)-u,m)}}}puede utilizarse.

Citas

  1. Armitage, JV; Eberlein, WF (2006). Funciones elípticas (Primera  ed.). Cambridge University Press. ISBN 978-0-521-78078-0.pág. 48
  2. 1 2 Walker, Peter (2003). "La analiticidad de las funciones jacobianas con respecto al parámetro k" . Actas de la Royal Society . 459 (2038): 2569–2574 . doi : 10.1098/rspa.2003.1157 . JSTOR 3560143 . 
  3. 1 2 3 4 5 Olver, FWJ; et al., eds. (2017-12-22). "Biblioteca digital de funciones matemáticas del NIST (versión 1.0.17)" . Instituto Nacional de Estándares y Tecnología . Recuperado el 26 de febrero de 2018 . 
  4. "cplot, paquete de Python para graficar funciones de valores complejos" . GitHub .
  5. 1 2 3 4 5 6 Neville, Eric Harold (1944). Funciones elípticas jacobianas . Oxford: Oxford University Press.
  6. 1 2 Sala, Kenneth L. (noviembre de 1989). "Transformaciones de la función de amplitud jacobiana y su cálculo mediante la media aritmético-geométrica". SIAM Journal on Mathematical Analysis . 20 (6): 1514– 1528. doi : 10.1137/0520100 .
  7. Reinhardt, WP; Walker, PL (2010), "Funciones elípticas jacobianas" , en Olver, Frank WJ ; Lozier, Daniel M.; Boisvert, Ronald F.; Clark, Charles W. (eds.), NIST Handbook of Mathematical Functions , Cambridge University Press, ISBN 978-0-521-19225-5, MR 2723248 .
  8. Carlson, BC (2010), "Integrales elípticas" , en Olver, Frank WJ ; Lozier, Daniel M.; Boisvert, Ronald F.; Clark, Charles W. (eds.), NIST Handbook of Mathematical Functions , Cambridge University Press, ISBN 978-0-521-19225-5, MR 2723248 .
  9. Whittaker, Edmund Taylor ; Watson, George Neville (1927). Un curso de análisis moderno (4.ª ed.). Cambridge University Press. pág. 492.  
  10. Cox, David Archibald (enero de 1984). "La media aritmético-geométrica de Gauss" . L'Enseignement Mathématique . 30 (2): 290.
  11. 1 2 3 4 5 6 "Introducción a las funciones elípticas de Jacobi" . El sitio web de funciones de Wolfram . Wolfram Research, Inc. 2018. Consultado el 7 de enero de 2018 .
  12. Whittaker, ET ; Watson, GN (1940). Un curso de análisis moderno . Nueva York, EE. UU.: The MacMillan Co.
  13. "Funciones elípticas: variables complejas" .
  14. Whittaker, Edmund Taylor ; Watson, George Neville (1927). Un curso de análisis moderno (4.ª ed.). Cambridge University Press. págs. 504–505 .  
  15. Dixon, Alfred Cardew (1894). Las propiedades elementales de las funciones elípticas, con ejemplos . Macmillan. págs. 73–75 . 
  16. Schett, Alois (1976). "Propiedades de los coeficientes de la expansión en serie de Taylor de las funciones elípticas jacobianas" . Math . Comp . 30 (133): 143–147 . doi : 10.1090/S0025-5718-1976-0391477-3.MR 0391477 . 
  17. Salzer, Herbert E. (julio de 1962). "Cálculo rápido de funciones elípticas jacobianas" . Communications of the ACM . 5 (7): 399. doi : 10.1145/368273.368573 .
  18. Smith, John I. (5 de mayo de 1971). "Los parámetros de capacitancia de modo par e impar para líneas acopladas en sustrato suspendido". IEEE Transactions on Microwave Theory and Techniques . MTT-19 (5): 430. doi : 10.1109/TMTT.1971.1127543 .
  19. Reinhardt, WP; Walker, PL (2010), "Funciones elípticas jacobianas" , en Olver, Frank WJ ; Lozier, Daniel M.; Boisvert, Ronald F.; Clark, Charles W. (eds.), NIST Handbook of Mathematical Functions , Cambridge University Press, ISBN 978-0-521-19225-5, MR 2723248 .
  20. Reinhardt, WP; Walker, PL (2010), "Funciones elípticas jacobianas" , en Olver, Frank WJ ; Lozier, Daniel M.; Boisvert, Ronald F.; Clark, Charles W. (eds.), NIST Handbook of Mathematical Functions , Cambridge University Press, ISBN 978-0-521-19225-5, MR 2723248 .
  21. Bagis, N. (2018). Evaluaciones de series relacionadas con la función elíptica de Jacobi (Preimpresión). arXiv : 1803.09445 .
  22. 1 2 3 4 H.S. Wall. (1948). "Teoría analítica de las fracciones continuas", Van Nostrand, Nueva York.
  23. ^ Perron, O. (1957). "Die Lehre von den Kettenbruchen", Banda II, BG Teubner, Stuttgart.
  24. Reinhardt, WP; Walker, PL (2010), "§22.15 Funciones inversas" , en Olver, Frank WJ ; Lozier, Daniel M.; Boisvert, Ronald F.; Clark, Charles W. (eds.), NIST Handbook of Mathematical Functions , Cambridge University Press, ISBN 978-0-521-19225-5, MR 2723248 .
  25. Ehrhardt, Wolfgang. "Funciones especiales de AMath y DAMath: Manual de referencia y notas de implementación" (PDF) . pág. 42. Archivado del original (PDF) el 31 de julio de 2016. Recuperado el 17 de julio de 2013 . 
  26. Byrd, PF; Friedman, MD (1971). Manual de integrales elípticas para ingenieros y científicos (2.ª ed.). Berlín: Springer-Verlag. 
  27. Carlson, BC (2008). "Series de potencias para funciones elípticas jacobianas inversas" (PDF) . Mathematics of Computation . 77 (263): 1615– 1621. doi : 10.1090/s0025-5718-07-02049-2 . Recuperado el 17 de julio de 2013 .

Referencias

  • Abramowitz, Milton ; Stegun, Irene Ann , eds. (1983) [junio de 1964]. «Capítulo 16» . Manual de funciones matemáticas con fórmulas, gráficas y tablas matemáticas . Serie de Matemáticas Aplicadas. Vol.  55 (novena reimpresión con correcciones adicionales de la décima edición original con correcciones (diciembre de 1972); primera  ed.). Washington D. C.; Nueva York: Departamento de Comercio de los Estados Unidos, Oficina Nacional de Normas; Dover Publications. pág.  569. ISBN 978-0-486-61272-0. LCCN 64-60036 . MR 0167642 . LCCN 65-12253 .   
  • NI Akhiezer , Elementos de la teoría de las funciones elípticas (1970) Moscú, traducido al inglés como AMS Translations of Mathematical Monographs Volumen 79 (1990) AMS, Rhode Island ISBN 0-8218-4532-2
  • AC Dixon, Las propiedades elementales de las funciones elípticas, con ejemplos (Macmillan, 1894)
  • Alfred George Greenhill, Las aplicaciones de las funciones elípticas (Londres, Nueva York, Macmillan, 1892)
  • Edmund T. Whittaker, George Neville Watson: Un curso de análisis moderno . 4.ª ed. Cambridge, Inglaterra: Cambridge University Press, 1990. Págs. 469-470.
  • Conferencias de H. Hancock sobre la teoría de las funciones elípticas (Nueva York, J. Wiley & Sons, 1910)
  • Jacobi, CGJ (2012) [1829], Fundamenta nova theoriae functionum ellipticarum (en latín), Königsberg, ISBN 978-1-108-05200-9Reimpreso por Cambridge University Press 2012
  • Reinhardt, William P.; Walker, Peter L. (2010), "Funciones elípticas jacobianas" , en Olver, Frank WJ ; Lozier, Daniel M.; Boisvert, Ronald F.; Clark, Charles W. (eds.), NIST Handbook of Mathematical Functions , Cambridge University Press, ISBN 978-0-521-19225-5, MR 2723248 .
  • (en francés) P. Appell y E. Lacour Principes de la théorie des fonctions elliptiques et apps (París, Gauthier Villars, 1897)
  • (en francés) GH Halphen Traité des fonctions elliptiques et de leurs apps (vol. 1) (París, Gauthier-Villars, 1886–1891)
  • (en francés) GH Halphen Traité des fonctions elliptiques et de leurs apps (vol. 2) (París, Gauthier-Villars, 1886–1891)
  • (en francés) GH Halphen Traité des fonctions elliptiques et de leurs apps (vol. 3) (París, Gauthier-Villars, 1886–1891)
  • (en francés) J. Tannery y J. Molk Eléments de la théorie des fonctions elliptiques. Tomo I, Introducción. Calculo diferencial. Ire partie (París  : Gauthier-Villars et fils, 1893)
  • (en francés) J. Tannery y J. Molk Eléments de la théorie des fonctions elliptiques. Tomo II, Calcul différentiel. IIe partie (París  : Gauthier-Villars et fils, 1893)
  • (en francés) J. Tannery y J. Molk Eléments de la théorie des fonctions elliptiques. Tomo III, Cálculo integral. Ire partie, Théorèmes généraux. Inversión (París  : Gauthier-Villars et fils, 1893)
  • (en francés) J. Tannery y J. Molk Eléments de la théorie des fonctions elliptiques. Tomo IV, Cálculo integral. IIe partie, Aplicaciones (París  : Gauthier-Villars et fils, 1893)
  • (en francés) C. Briot y JC Bouquet Théorie des fonctions elliptiques (París  : Gauthier-Villars, 1875)
  • Toshio Fukushima: Cálculo rápido de integrales elípticas completas y funciones elípticas jacobianas . 2012, Observatorio Astronómico Nacional de Japón (国立天文台)
  • Lowan, Blanch y Horenstein: Sobre la inversión de la serie q asociada a las funciones elípticas jacobianas . Bull. Amer. Math. Soc. 48, 1942
  • H. Ferguson, DE Nielsen, G. Cook: Una fórmula de partición para los coeficientes enteros de la función theta . Matemáticas de la computación, Volumen 29, Número 131, Julio de 1975.
  • JD Fenton y RS Gardiner-Garden: Métodos de convergencia rápida para evaluar integrales elípticas y funciones elípticas y theta . J. Austral. Math. Soc. (Serie B) 24, 1982, pág. 57
  • Adolf Kneser: Neue Untersuchung einer Reihe aus der Theorie der elliptischen Funktionen . J. reina u. ángulo. Matemáticas. 157, 1927. páginas 209 – 218