Fonction elliptique de Jacobi

En mathématiques, les fonctions elliptiques de Jacobi sont des fonctions elliptiques d'une grande importance historique.

Introduites par Carl Gustav Jakob Jacobi vers 1830, elles ont des applications directes, par exemple dans l'équation du pendule. Elles présentent aussi des analogies avec les fonctions trigonométriques, qui sont mises en valeur par le choix des notations $sn$ et $cn$ , qui rappellent $sin$ et $cos$ . Si les fonctions elliptiques thêta de Weierstrass semblent mieux adaptées aux considérations théoriques, les problèmes physiques pratiques font plus appel aux fonctions de Jacobi.

Introduction

Construction auxiliaire : le rectangle de Jacobi.

Il existe 12 fonctions elliptiques de Jacobi.

Ce sont des fonctions d'une variable complexe mais qui dépendent d'un paramètre $k$ élément de ]0,1[, sous-entendu dans les notations. $k$ s'appelle le module des fonctions de Jacobi. À ce paramètre $k$ , on associe les deux nombres $K$ et $K'$ , définis par les intégrales elliptiques $K=\int _{0}^{\pi /2}{\frac {\mathrm {d} \theta }{\sqrt {1-k^{2}\sin ^{2}\theta }}}$ et $K'=\int _{0}^{\pi /2}{\frac {\mathrm {d} \theta }{\sqrt {1-(1-k^{2})\sin ^{2}\theta }}}$ , ainsi que le nombre $k'={\sqrt {1-k^{2}}}$ , appelé comodule.

Dans le plan complexe, on dispose un rectangle dont les quatre sommets sont conventionnellement notés $s$ , $c$ , $d$ et $n$ , de façon que $s$ soit à l'origine, $c$ au point d'abscisse $K$ sur l'axe des réels, $d$ au point d'affixe complexe $K + i K'$ , et $n$ au point d'affixe $i K'$ sur l'axe imaginaire.

Le nom de chacune des fonctions de Jacobi est alors associé à un couple formé de deux sommets du rectangle. Ainsi, les noms des 12 fonctions elliptiques de Jacobi sont : $sc$ , $sd$ , $sn$ , $cd$ , $cn$ , $cs$ , $dn$ , $ds$ , $dc$ , $ns$ , $nc$ et $nd$ .

Pour tout sommet $p$ parmi les quatre sommets $s$ $c$ $d$ $n$ , et pour tout sommet $q$ pris parmi les trois sommets restants, la fonction de Jacobi $pq$ est la seule fonction de la variable complexe qui soit doublement périodique et méromorphe, et qui vérifie les propriétés suivantes[1] :

Elle admet un zéro simple au sommet $p$ , et un pôle simple au sommet $q$ .
Elle est périodique de période $4 K$ selon l'axe réel, et périodique de période $4 K'$ selon l'axe imaginaire. Les nombres $K$ et $K'$ sont appelés « quarts de période ».
Elle est périodique dans la direction $pq$ , de période double de la distance de $p$ à $q$ .
Le coefficient du premier terme de son développement en série au voisinage de $u$ = 0 vaut 1. Autrement dit, ce premier terme vaut $u$ , 1/ $u$ ou 1 selon que le sommet correspondant à $u$ = 0 est un zéro, un pôle ou un point ordinaire de la fonction.

Fonction

sn

Fonction

cn

Fonction

dn

Représentation des trois fonctions elliptiques de Jacobi

sn

,

dn

et

dn

dans le plan complexe, en utilisant la méthode de coloration de régions. Ces représentations mettent en évidence leur double comportement périodique. La valeur du module

k

est 0.8.

Dans un cadre plus général, $k$ est complexe, ainsi que $K$ et $K'$ , et on travaille à partir d'un parallélogramme. Cependant, si $K$ et $K'$ sont réels , alors les fonctions elliptiques de Jacobi prennent des valeurs réelles lorsqu'elles sont appliquées à une variable réelle.

Définition

Parmi les douze fonctions de Jacobi, on en distingue trois, appelées fonctions de base de Jacobi. Ce sont $sn$ , $cn$ et $dn$ . On définit les autres fonctions de Jacobi à partir de celles-ci. Pour définir les trois fonctions de base, on introduit une fonction intermédiaire, la fonction amplitude de Jacobi.

Intégrale elliptique incomplète de première espèce et fonction amplitude

On rappelle que l'intégrale elliptique incomplète de première espèce associée au module $k$ est la fonction impaire croissante sur les réels définie par :

F(a,k)=\int _{0}^{a}{\frac {\mathrm {d} x}{\sqrt {1-k^{2}\sin ^{2}x}}}

On remarque que la constante $K$ définie précédemment n'est autre que $F({\frac {\pi }{2}},k)$ . Elle est appelée intégrale elliptique complète de première espèce.

On appelle fonction amplitude de Jacobi[1] la fonction réciproque de $F$ , notée $A$ :

u=F(a,k)\iff a=A(u,k)

Elle est elle-même impaire et croissante sur les réels, et augmente de π lorsque $u$ augmente de $2 K$ .

Les trois fonctions de base de Jacobi (1827)

On les définit comme suit[1] :

la fonction sinus de Jacobi : ${\rm {sn}}(u,k)=\sin(A(u,k))$ . Sur les réels, elle est périodique de période $4 K$ .
la fonction cosinus de Jacobi : ${\rm {cn}}(u,k)=\cos(A(u,k))$ . Sur les réels, elle est périodique de période $4 K$ .
la fonction dn de Jacobi : ${\rm {dn}}(u,k)={\sqrt {1-k^{2}{\rm {sn}}(u,k)^{2}}}$ . Sur les réels, elle est périodique de période $2 K$ .

$sn$ est une fonction impaire, alors que $cn$ et $dn$ sont paires.

Cas limites

On retrouve les fonctions trigonométriques circulaires et hyperboliques[2] pour les valeurs limites 0 et 1 de $k$ :

Si $k$ = 0, on retrouve la trigonométrie ordinaire. En effet, $F(a,0)=a$ , $A(u)=u$ , $K={\frac {\pi }{2}}$ et $K'$ est envoyé à l'infini. $sn$ est le sinus, $cn$ le cosinus et $dn$ la fonction constante 1.
Si $k$ = 1, on voit apparaître les fonctions de la trigonométrie hyperbolique. En effet, $u=F(a,1)={\rm {argth}}(\sin(a))$ de sorte que $\sin(a)=\tanh(u)$ (formule de Gudermann), $K$ est envoyé à l'infini et $K'={\frac {\pi }{2}}$ . $sn$ est la fonction $tanh$ , $cn$ et $dn$ la fonction $1/cosh$ .

Fonction

sn

Fonction

cn

Fonction

dn

Fonction

amplitude A

Représentation des trois fonctions elliptiques de base de Jacobi et de la fonction amplitude, entre -5 et 5, pour diverses valeurs du module. Les cas limites correspondant au module égal à 0 ou 1 sont en vert.

Les autres fonctions

Gudermann (1838), puis Glaisher (1882) introduiront les neuf autres fonctions[3] :

{\rm {ns}}(u,k)={\frac {1}{{\rm {sn}}(u,k)}}

,

{\rm {nc}}(u,k)={\frac {1}{{\rm {cn}}(u,k)}}

,

{\rm {nd}}(u,k)={\frac {1}{{\rm {dn}}(u,k)}}

{\rm {sc}}(u,k)={\frac {{\rm {sn}}(u,k)}{{\rm {cn}}(u,k)}}

,

{\rm {cs}}(u,k)={\frac {{\rm {cn}}(u,k)}{{\rm {sn}}(u,k)}}

{\rm {sd}}(u,k)={\frac {{\rm {sn}}(u,k)}{{\rm {dn}}(u,k)}}

,

{\rm {ds}}(u,k)={\frac {{\rm {dn}}(u,k)}{{\rm {sn}}(u,k)}}

{\rm {cd}}(u,k)={\frac {{\rm {cn}}(u,k)}{{\rm {dn}}(u,k)}}

,

{\rm {dc}}(u,k)={\frac {{\rm {dn}}(u,k)}{{\rm {cn}}(u,k)}}

Fonctions réciproques

On peut définir les fonctions réciproques des fonctions elliptiques de Jacobi, pour x entre -1 et 1[4] :

$\mathrm {Arcsn} \,(x,k)=\int _{0}^{x}{\frac {\mathrm {d} t}{\sqrt {(1-t^{2})(1-k^{2}t^{2})}}}$
$\mathrm {Arccn} \,(x,k)=\int _{x}^{1}{\frac {\mathrm {d} t}{\sqrt {(1-t^{2})(1-k^{2}+k^{2}t^{2})}}}$

Formulaire

Valeurs remarquables

Pour les valeurs réelles de la variable[2] :

pour u = 0, on a ${\text{sn}}(0)=0,\,{\text{cn}}(0)=1,\,{\text{dn}}(0)=1.$
pour u = K/2, on a ${\text{sn}}\left({\frac {K}{2}}\right)={\frac {1}{\sqrt {1+k'}}},\,{\text{cn}}\left({\frac {K}{2}}\right)={\sqrt {\frac {k'}{1+k'}}},\,{\text{dn}}\left({\frac {K}{2}}\right)={\sqrt {k'}}.$
pour u = K, on a ${\text{sn}}(K)=1,\,{\text{cn}}(K)=0,\,{\text{dn}}(K)=k'$

Dérivées

Les dérivées des fonctions de base sont[5] :

A'(u)={\text{dn}}(u)

{\text{sn}}'(u)={\text{cn}}(u){\text{dn}}(u)

{\text{cn}}'(u)=-{\text{sn}}(u){\text{dn}}(u)

{\text{dn}}'(u)=-k^{2}{\text{sn}}(u){\text{cn}}(u)

Translation

On dispose des relations suivantes[6] :

${\text{sn}}(u+K)={\frac {{\text{cn}}(u)}{{\text{dn}}(u)}}\,$
${\text{cn}}(u+K)=-k'{\frac {{\text{sn}}(u)}{{\text{dn}}(u)}}$
${\text{dn}}(u+K)={\frac {k'}{{\text{dn}}(u)}}\,$
${\text{sn}}(u+2K)=-{\text{sn}}(u)$
${\text{cn}}(u+2K)=-{\text{cn}}(u)$
${\text{dn}}(u+2K)={\text{dn}}(u)$

Addition

On dispose des formules d'addition suivantes, généralisant les formules d'addition trigonométriques[5] :

${\text{sn}}(u+v)={\frac {{\text{sn}}(u){\text{cn}}(v){\text{dn}}(v)+{\text{sn}}(v){\text{cn}}(u){\text{dn}}(u)}{1-k^{2}{\text{sn}}^{2}(u){\text{sn}}^{2}(v)}}$
${\text{cn}}(u+v)={\frac {{\text{cn}}(u){\text{cn}}(v)-{\text{sn}}(u){\text{dn}}(u){\text{sn}}(v){\text{dn}}(v)}{1-k^{2}{\text{sn}}^{2}(u){\text{sn}}^{2}(v)}}$
${\text{dn}}(u+v)={\frac {{\text{dn}}(u){\text{dn}}(v)-k^{2}{\text{sn}}(u){\text{cn}}(u){\text{sn}}(v){\text{cn}}(v)}{1-k^{2}{\text{sn}}^{2}(u){\text{sn}}^{2}(v)}}$

Carrés

${\text{sn}}^{2}(u)+{\text{cn}}^{2}(u)=1\,$
$k^{2}{\text{sn}}^{2}(u)+{\text{dn}}^{2}(u)=1\,$
$k^{2}{\text{cn}}^{2}(u)+k'^{2}={\text{dn}}^{2}(u)\,$ avec $k'={\sqrt {1-k^{2}}}$ le complément du module $k$ .
${\text{cn}}^{2}(u)+k'^{2}{\text{sn}}^{2}(u)={\text{dn}}^{2}(u)\,$

Transformées des carrés en arc double

${\text{sn}}^{2}(u)={\frac {1-{\text{cn}}(2u)}{1+{\text{dn}}(2u)}}$
${\text{cn}}^{2}(u)={\frac {{\text{dn}}(2u)+{\text{cn}}(2u)}{1+{\text{dn}}(2u)}}$
${\text{dn}}^{2}(u)={\frac {{\text{dn}}(2u)+{\text{cn}}(2u)}{1+{\text{cn}}(2u)}}$

Équations différentielles

Les règles de dérivation des fonctions de Jacobi permettent de montrer que $sn$ , $cn$ et $dn$ sont respectivement solutions des équations différentielles suivantes[7] :

$sn$ : ${\ddot {x}}+(1+k^{2})x-2k^{2}x^{3}=0$
$cn$ : ${\ddot {x}}+(1-2k^{2})x+2k^{2}x^{3}=0$
$dn$ : ${\ddot {x}}-(2-k^{2})x+2x^{3}=0$

Applications

Le pendule simple oscillant

Schéma du pendule simple.

On considère un pendule simple, de longueur $l$ , oscillant dans un champ de pesanteur $g$ . Soit $θ$ l'angle qu'il forme avec la verticale descendante, et $θ 0$ son amplitude maximale. $θ$ vérifie l'équation du mouvement suivante (provenant de la conservation de l'énergie mécanique du pendule) :

{{\dot {\theta }}^{2}}=2{\frac {g}{l}}(\cos \theta -\cos \theta _{0})

La solution de cette équation qui s'annule au temps $t$ = 0 vérifie :

\sin \left({\frac {\theta }{2}}\right)=\sin \left({\frac {\theta _{0}}{2}}\right)\,\mathrm {sn} (\omega _{0}t,k)

où l'on a donné au module de la fonction de Jacobi la valeur $k=\sin \left({\frac {\theta _{0}}{2}}\right)$ , et où $\omega _{0}={\sqrt {\frac {g}{l}}}$ est la pulsation du pendule simple pour les petites amplitudes.

Démonstration

Vérifions que la solution $θ$ proposée satisfait l'équation du mouvement. On dérive l'égalité $\sin \left({\frac {\theta }{2}}\right)=\sin \left({\frac {\theta _{0}}{2}}\right)\,\mathrm {sn} (\omega _{0}t)=k\,\mathrm {sn} (\omega _{0}t)$ par rapport au temps, ce qui donne :

{\frac {1}{2}}\cos \left({\frac {\theta }{2}}\right)\,{\dot {\theta }}=k\,\omega _{0}\,\mathrm {cn} (\omega _{0}t)\,\mathrm {dn} (\omega _{0}t)

On élève au carré, et on exprime $cn$ et $dn$ en fonction de $sn$ , puis on utilise des identités trigonométriques :

{\begin{aligned}{\frac {1}{4}}\cos ^{2}\left({\frac {\theta }{2}}\right){\dot {\theta }}^{2}&=k^{2}\omega _{0}^{2}\left(1-\mathrm {sn} ^{2}(\omega _{0}t)\right)\left(1-k^{2}\mathrm {sn} ^{2}(\omega _{0}t)\right)\\&=\omega _{0}^{2}\left(k^{2}-k^{2}\,\mathrm {sn} ^{2}(\omega _{0}t)\right)\left(1-\sin ^{2}\left({\frac {\theta }{2}}\right)\right)\\&=\omega _{0}^{2}\left(\sin ^{2}\left({\frac {\theta _{0}}{2}}\right)-\sin ^{2}\left({\frac {\theta }{2}}\right)\right)\cos ^{2}\left({\frac {\theta }{2}}\right)\\&=\omega _{0}^{2}{\frac {\cos \theta -\cos \theta _{0}}{2}}\cos ^{2}\left({\frac {\theta }{2}}\right)\\\end{aligned}}

Sachant que $\omega _{0}={\sqrt {\frac {g}{l}}}$ , on obtient bien ${{\dot {\theta }}^{2}}=2{\frac {g}{l}}(\cos \theta -\cos \theta _{0})$ .

Le calcul ci-dessus montre également que :

{\dot {\theta }}=2k\,\omega _{0}\,\mathrm {cn} (\omega _{0}t)

La période du pendule est $T={\frac {4K(k)}{\omega _{0}}}$ . La fonction amplitude $A(\omega _{0}t,k)$ croît avec $t$ , et joue le rôle d'« échelle de temps » adaptée au problème : à chaque période de temps réel du pendule, l'amplitude aura augmenté de $2π$ . L'anisochronicité du mouvement est patente, puisque la période du pendule dépend du module $k$ , donc de $θ 0$ .

Pour les petites oscillations, $k$ est très petit, de sorte que la fonction $sn$ est assimilable au sinus. En approximant le sinus de $θ$ par $θ$ et en faisant de même pour $θ 0$ , on retrouve la formule classique $\theta =\theta _{0}\sin(\omega _{0}t)$ .

Quand $θ 0$ tend vers $π$ , $k$ tend vers 1 et $K (k)$ tend vers l'infini comme $\ln \left({\frac {4}{\sqrt {1-k^{2}}}}\right)$ . Si $T 0$ est la période ${\frac {2\pi }{\omega _{0}}}$ du pendule simple pour les petites oscillations, alors, la période du pendule devient :

T=T_{0}\,{\frac {2}{\pi }}\ln \left({\frac {8}{\pi -\theta _{0}}}\right)

.

Quand la limite est atteinte, $sn$ est égale à la fonction $tanh$ . On a alors :

\theta =2\arcsin(\tanh(\omega _{0}t))=4\arctan(\exp(\omega _{0}t))-\pi

qui tend vers $π$ quand $t$ tend vers l'infini.

Le pendule simple tournoyant

Pendule tournoyant (dans le sens inverse au sens trigonométrique). La courbe paramétrée en haut de l'image représente le mouvement du pendule dans l'espace des phases :

\theta

en abscisse (proportionnelle à la fonction amplitude

A

de Jacobi), l'axe étant orienté vers la gauche, et

{\dot {\theta }}

en ordonnée (proportionnelle à la fonction de Jacobi périodique

dn

).

Dans le cas d'un pendule animé d'une vitesse suffisamment grande pour le faire tournoyer, l'équation du mouvement s'écrit :

{\dot {\theta }}^{2}={\frac {2g}{l^{2}}}[H-l+l\cos(\theta )]

où $H$ est une constante homogène à une longueur et strictement supérieure à $2 l$ . La solution $θ$ s'exprime alors à l'aide de la fonction amplitude de Jacobi sous la forme :

\theta =2A\left({\sqrt {\frac {Hg}{2l^{2}}}}\,t,k\right)

où l'on donne au module de la fonction de Jacobi la valeur $k={\sqrt {\frac {2l}{H}}}$ .

Démonstration

Si on dérive l'égalité par rapport au temps, on obtient :

{\dot {\theta }}=2{\sqrt {\frac {Hg}{2l^{2}}}}\mathrm {dn} \left({\sqrt {\frac {Hg}{2l^{2}}}}\,t\right)

donc, en élevant au carré, et en tenant compte du fait que $\sin \left({\frac {\theta }{2}}\right)=\sin \left[A\left({\sqrt {\frac {Hg}{2l^{2}}}}\,t\right)\right]=\mathrm {sn} \left({\sqrt {\frac {Hg}{2l^{2}}}}\,t\right)$ :

{\begin{aligned}{\dot {\theta }}^{2}&={\frac {2Hg}{l^{2}}}\left[1-k^{2}\mathrm {sn} ^{2}\left({\sqrt {\frac {Hg}{2l^{2}}}}\,t\right)\right]\\&={\frac {2Hg}{l^{2}}}\left[1-{\frac {2l}{H}}\sin ^{2}\left({\frac {\theta }{2}}\right)\right]\\&={\frac {2g}{l^{2}}}[H-l(1-\cos(\theta ))]\end{aligned}}

comme attendu.

Le mouvement à la Poinsot d'un solide

Ce mouvement est celui d'un solide en rotation, pris relativement à son centre d'inertie $G$ , quand le moment par rapport à $G$ des forces extérieures est nul. Pour un solide quelconque sans symétrie particulière, les équations du mouvement se résolvent à l'aide des fonctions elliptiques de Jacobi. En particulier, les trois composantes du vecteur de rotation instantanée dans le référentiel lié au solide constitué des axes principaux d'inertie sont proportionnelles respectivement à $cn$ , $sn$ , $dn$ [8].

Propagation d'onde

La fonction ${\displaystyle \xi$ permet de modéliser la surélévation de la surface de l'eau au passage d'un soliton, tel qu'un tsunami par exemple, où, à changement d'unité près, $ξ$ est la hauteur de la vague, $x$ est l'abscisse où l'on mesure cette hauteur, $t$ est le temps, et $B$ un paramètre prenant en compte la profondeur du milieu. C'est en effet une des solutions de l'équation de Korteweg et de Vries. L'onde ainsi modélisée s'appelle onde cnoïdale.

Pompage optique

La fonction $sn$ intervient pour modéliser la déplétion de la pompe dans le mélange à trois ondes optiques[9], qui est utilisé dans les Oscillateurs paramétriques optiques.

Voir aussi

Bibliographie

M. Abramowitz, I. A. Stegun, Handbook of Mathematical Functions, National Bureau of Standards, 1972 (lire en ligne), chapitre 16, par L. M. Milne-Thomson.
Hermann Laurent, Théorie élémentaire des fonctions elliptiques, Gauthier-Villars, Paris, 1880 (lire en ligne)
Alfred George Greenhill, Les fonctions elliptiques et leurs applications, G. Carré, Paris, 1895 (lire en ligne)
Paul Appell, Émile Lacour, Principes de la théorie des fonctions elliptiques et applications, Gauthier-Villars et fils, Paris, 1897 (lire en ligne)

Liens externes

W. P. Reinhardt, P. L. Walker, « Jacobian Elliptic Functions ». Parmi les nombreuses propriétés des fonctions elliptiques de Jacobi que ce site donne, on trouvera en particulier au chapitre 22.20 des méthodes de calcul numérique rapide de ces fonctions.

Notes et références

Abramowitz-Stegun 1972, p. 569
Abramowitz-Stegun 1972, p. 571
Abramowitz-Stegun 1972, p. 570
W. P. Reinhardt, P. L. Walker, « Jacobian Elliptic Functions », sur dlmf.nist.gov, §22.15, Inverse Functions
Abramowitz-Stegun 1972, p. 574
Abramowitz-Stegun 1972, p. 572
W. P. Reinhardt, P. L. Walker, « Jacobian Elliptic Functions », sur dlmf.nist.gov, §22.13, Derivatives and Differential Equations
L. Landau, E. Lifchitz, Physique théorique, mécanique, Ellipses, 1994, p. 176
Paul Elwyn Britton, « Fibre laser pumped periodically poled lithium niobate based nonlinear devices », sur University of Southampton, 2000, p. 101, ch.5 (« Parametric amplification and generation »)

Portail de l'analyse

Cet article est issu de Wikipedia. Le texte est sous licence Creative Commons - Attribution - Partage dans les Mêmes. Des conditions supplémentaires peuvent s'appliquer aux fichiers multimédias.

[ref01-1] Abramowitz-Stegun 1972, p. 569

[ref02-2] Abramowitz-Stegun 1972, p. 571

[3] Abramowitz-Stegun 1972, p. 570

[4] W. P. Reinhardt, P. L. Walker, « Jacobian Elliptic Functions », sur dlmf.nist.gov, §22.15, Inverse Functions

[ref03-5] Abramowitz-Stegun 1972, p. 574

[6] Abramowitz-Stegun 1972, p. 572

[7] W. P. Reinhardt, P. L. Walker, « Jacobian Elliptic Functions », sur dlmf.nist.gov, §22.13, Derivatives and Differential Equations

[8] L. Landau, E. Lifchitz, Physique théorique, mécanique, Ellipses, 1994, p. 176

[9] Paul Elwyn Britton, « Fibre laser pumped periodically poled lithium niobate based nonlinear devices », sur University of Southampton, 2000, p. 101, ch.5 (« Parametric amplification and generation »)