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 Modèle:Math et Modèle:Math, qui rappellent Modèle:Math et Modèle:Math. 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
Il existe 12 fonctions elliptiques de Jacobi.
Ce sont des fonctions d'une variable complexe mais qui dépendent d'un paramètre Modèle:Math élément de ]0,1[, sous-entendu dans les notations. Modèle:Math s'appelle le module des fonctions de Jacobi. À ce paramètre Modèle:Math, on associe les deux nombres Modèle:Math et Modèle:Math, définis par les intégrales elliptiques <math>K = \int_0^{\pi/2} \frac{\mathrm{d}\theta}{\sqrt{1-k^2 \sin^2\theta}}</math> et <math>K' = \int_0^{\pi/2} \frac{\mathrm{d}\theta}{\sqrt{1-(1-k^2)\sin^2\theta}}</math>, ainsi que le nombre <math>k'= \sqrt{1-k^2}</math>, appelé comodule.
Dans le plan complexe, on dispose un rectangle dont les quatre sommets sont conventionnellement notés Modèle:Math, Modèle:Math, Modèle:Math et Modèle:Math, de façon que Modèle:Math soit à l'origine, Modèle:Math au point d'abscisse Modèle:Math sur l'axe des réels, Modèle:Math au point d'affixe complexe Modèle:Math, et Modèle:Math au point d'affixe Modèle:Math 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 : Modèle:Math, Modèle:Math, Modèle:Math, Modèle:Math, Modèle:Math, Modèle:Math, Modèle:Math, Modèle:Math, Modèle:Math, Modèle:Math, Modèle:Math et Modèle:Math.
Pour tout sommet Modèle:Math parmi les quatre sommets Modèle:MathModèle:MathModèle:MathModèle:Math, et pour tout sommet Modèle:Math pris parmi les trois sommets restants, la fonction de Jacobi Modèle:Math 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<ref name="ref01">Modèle:Harvsp</ref> :
- Elle admet un zéro simple au sommet Modèle:Math, et un pôle simple au sommet Modèle:Math.
- Elle est périodique de période Modèle:Math selon l'axe réel, et périodique de période Modèle:Math selon l'axe imaginaire. Les nombres Modèle:Math et Modèle:Math sont appelés « quarts de période ».
- Elle est périodique dans la direction Modèle:Math, de période double de la distance de Modèle:Math à Modèle:Math.
- Le coefficient du premier terme de son développement en série au voisinage de Modèle:Math = 0 vaut 1. Autrement dit, ce premier terme vaut Modèle:Math, 1/Modèle:Math ou 1 selon que le sommet correspondant à Modèle:Math = 0 est un zéro, un pôle ou un point ordinaire de la fonction.
Dans un cadre plus général, Modèle:Math est complexe, ainsi que Modèle:Math et Modèle:Math, et on travaille à partir d'un parallélogramme. Cependant, si Modèle:Math et Modèle:Math 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 Modèle:Math, Modèle:Math et Modèle:Math. 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 Modèle:Math est la fonction impaire croissante sur les réels définie par :
- <math>F(a,k) = \int_0^{a} \frac{\mathrm dx}{\sqrt{1-k^2\sin^2x}}</math>
On remarque que la constante Modèle:Math définie précédemment n'est autre que <math display=inline>F(\frac{\pi}{2},k)</math>. Elle est appelée intégrale elliptique complète de première espèce.
On appelle fonction amplitude de Jacobi<ref name="ref01" /> la fonction réciproque de Modèle:Math, notée Modèle:Math :
- <math>u = F(a,k) \iff a = A(u,k)</math>
Elle est elle-même impaire et croissante sur les réels, et augmente de π lorsque Modèle:Math augmente de Modèle:Math.
Les trois fonctions de base de Jacobi (1827)
On les définit comme suit<ref name="ref01" /> :
- la fonction sinus de Jacobi : <math>{\rm sn}(u,k) = \sin(A(u,k))</math>. Sur les réels, elle est périodique de période Modèle:Math.
- la fonction cosinus de Jacobi : <math>{\rm cn}(u,k) = \cos(A(u,k))</math>. Sur les réels, elle est périodique de période Modèle:Math.
- la fonction dn de Jacobi : <math>{\rm dn}(u,k) = \sqrt{1 - k^2 {\rm sn}(u,k)^2}</math>. Sur les réels, elle est périodique de période Modèle:Math.
Modèle:Math est une fonction impaire, alors que Modèle:Math et Modèle:Math sont paires.
Cas limites
On retrouve les fonctions trigonométriques circulaires et hyperboliques<ref name="ref02">Modèle:Harvsp</ref> pour les valeurs limites 0 et 1 de Modèle:Math :
- Si Modèle:Math = 0, on retrouve la trigonométrie ordinaire. En effet, <math>F(a,0) = a</math>, <math>A(u) = u</math>, <math>K=\frac{\pi}{2}</math> et Modèle:Math est envoyé à l'infini. Modèle:Math est le sinus, Modèle:Math le cosinus et Modèle:Math la fonction constante 1.
- Si Modèle:Math = 1, on voit apparaître les fonctions de la trigonométrie hyperbolique. En effet, <math>u = F(a,1) = {\rm argth}(\sin(a))</math> de sorte que <math>\sin(a)=\tanh(u)</math> (formule de Gudermann), Modèle:Math est envoyé à l'infini et <math>K'= \frac{\pi}{2}</math>. Modèle:Math est la fonction Modèle:Math, Modèle:Math et Modèle:Math la fonction Modèle:Math.
Les autres fonctions
Gudermann (1838), puis Glaisher (1882) introduiront les neuf autres fonctions<ref>Modèle:Harvsp</ref> :
- <math>{\rm ns}(u,k) = \frac{1}{{\rm sn}(u,k)}</math>, <math>{\rm nc}(u,k) = \frac{1}{{\rm cn}(u,k)}</math>, <math>{\rm nd}(u,k) = \frac{1}{{\rm dn}(u,k)}</math>
- <math>{\rm sc}(u,k) = \frac{{\rm sn}(u,k)}{{\rm cn}(u,k)}</math>, <math>{\rm cs}(u,k) = \frac{{\rm cn}(u,k)}{{\rm sn}(u,k)}</math>
- <math>{\rm sd}(u,k) = \frac{{\rm sn}(u,k)}{{\rm dn}(u,k)}</math>, <math>{\rm ds}(u,k) = \frac{{\rm dn}(u,k)}{{\rm sn}(u,k)}</math>
- <math>{\rm cd}(u,k) = \frac{{\rm cn}(u,k)}{{\rm dn}(u,k)}</math>, <math>{\rm dc}(u,k) = \frac{{\rm dn}(u,k)}{{\rm cn}(u,k)}</math>
Fonctions réciproques
On peut définir les fonctions réciproques des fonctions elliptiques de Jacobi, pour x entre -1 et 1<ref>Modèle:Lien web, §22.15, Inverse Functions</ref> :
- <math>\mathrm{Arcsn}\,(x,k) = \int_0^x \frac{\mathrm{d}t}{\sqrt{(1-t^2)(1-k^2t^2)}}</math>
- <math>\mathrm{Arccn}\,(x,k) =\int_x^1 \frac{\mathrm{d}t}{\sqrt{(1-t^2)(1-k^2+k^2t^2)}}</math>
Formulaire
Valeurs remarquables
Pour les valeurs réelles de la variable<ref name="ref02" /> :
- pour u = 0, on a <math>\text{sn}(0) = 0,\,\text{cn}(0) = 1,\,\text{dn}(0) = 1.</math>
- pour u = K/2, on a <math>\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'}.</math>
- pour u = K, on a <math>\text{sn}(K) = 1,\,\text{cn}(K) = 0,\,\text{dn}(K) = k'</math>
Dérivées
Les dérivées des fonctions de base sont<ref name="ref03">Modèle:Harvsp</ref> :
- <math> A'(u) = \text{dn}(u)</math>
- <math> \text{sn}'(u) = \text{cn}(u)\text{dn}(u)</math>
- <math> \text{cn}'(u) = -\text{sn}(u)\text{dn}(u)</math>
- <math> \text{dn}'(u) = -k^2\text{sn}(u)\text{cn}(u)</math>
Translation
On dispose des relations suivantes<ref>Modèle:Harvsp</ref> :
- <math>\text{sn}(u + K) = \frac{\text{cn}(u)}{\text{dn}(u)}\,</math>
- <math>\text{cn}(u + K) = -k' \frac{\text{sn}(u)}{\text{dn}(u)}</math>
- <math>\text{dn}(u + K) = \frac{k'}{\text{dn}(u)}\,</math>
- <math>\text{sn}(u + 2K) = -\text{sn}(u)</math>
- <math>\text{cn}(u + 2K) = -\text{cn}(u)</math>
- <math>\text{dn}(u + 2K) = \text{dn}(u)</math>
Relations trigonométriques
Addition
On dispose des formules d'addition suivantes, généralisant les formules d'addition trigonométriques<ref name="ref03" /> :
- <math>\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)}</math>
- <math>\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)}</math>
- <math>\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)}</math>
Carrés
- <math>\text{sn}^2(u) + \text{cn}^2(u) = 1\,</math>
- <math>k^2 \text{sn}^2(u) + \text{dn}^2(u) = 1\,</math>
- <math>k^2 \text{cn}^2(u) + k'^2 = \text{dn}^2(u)\,</math> avec <math>k' = \sqrt{1-k^2}</math> le complément du module Modèle:Math.
- <math>\text{cn}^2(u) + k'^2 \text{sn}^2(u) = \text{dn}^2(u)\,</math>
Transformées des carrés en arc double
- <math>\text{sn}^2(u) = \frac{1-\text{cn}(2u)}{1+\text{dn}(2u)}</math>
- <math>\text{cn}^2(u) = \frac{\text{dn}(2u)+\text{cn}(2u)}{1+\text{dn}(2u)}</math>
- <math>\text{dn}^2(u) = \frac{\text{dn}(2u)+\text{cn}(2u)}{1+\text{cn}(2u)}</math>
Équations différentielles
Les règles de dérivation des fonctions de Jacobi permettent de montrer que Modèle:Math, Modèle:Math et Modèle:Math sont respectivement solutions des équations différentielles suivantes<ref>Modèle:Lien web, §22.13, Derivatives and Differential Equations</ref> :
- Modèle:Math : <math>\ddot{x} + (1+k^2)x -2k^2 x^3 = 0</math>
- Modèle:Math : <math>\ddot{x} + (1-2k^2)x +2k^2 x^3 = 0</math>
- Modèle:Math : <math>\ddot{x} - (2-k^2)x +2 x^3 = 0</math>
Applications
Le pendule simple oscillant
On considère un pendule simple, de longueur Modèle:Mvar, oscillant dans un champ de pesanteur Modèle:Mvar. Soit Modèle:Mvar l'angle qu'il forme avec la verticale descendante, et Modèle:Math son amplitude maximale. Modèle:Mvar vérifie l'équation du mouvement suivante (provenant de la conservation de l'énergie mécanique du pendule) :
- <math>{\dot\theta^2} = 2\frac{g}{l} (\cos \theta - \cos \theta_0)</math>
La solution de cette équation qui s'annule au temps Modèle:Mvar = 0 vérifie :
- <math>\sin\left(\frac{\theta}{2}\right) = \sin\left(\frac{\theta_0}{2}\right)\,\mathrm{sn}(\omega_0t,k)</math>
où l'on a donné au module de la fonction de Jacobi la valeur <math>k = \sin\left(\frac{\theta_0}{2}\right)</math>, et où <math>\omega_0 = \sqrt{\frac{g}{l}}</math> est la pulsation du pendule simple pour les petites amplitudes.
Modèle:Démonstration</math>, on obtient bien <math>{\dot\theta^2} = 2\frac{g}{l} (\cos \theta - \cos \theta_0)</math>.
Le calcul ci-dessus montre également que :
- <math>\dot\theta= 2k \, \omega_0 \,\mathrm{cn}(\omega_0t)</math>
}}
La période du pendule est <math>T = \frac{4K(k)}{\omega_0}</math>. La fonction amplitude <math>A(\omega_0 t,k)</math> croît avec Modèle:Mvar, 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 Modèle:Math. L'anisochronicité du mouvement est patente, puisque la période du pendule dépend du module Modèle:Mvar, donc de Modèle:Math.
Pour les petites oscillations, Modèle:Mvar est très petit, de sorte que la fonction Modèle:Math est assimilable au sinus. En approximant le sinus de Modèle:Mvar par Modèle:Mvar et en faisant de même pour Modèle:Math, on retrouve la formule classique <math>\theta = \theta_0 \sin(\omega_0 t)</math>.
Quand Modèle:Math tend vers Modèle:Math, Modèle:Mvar tend vers 1 et Modèle:Math tend vers l'infini comme <math>\ln\left(\frac{4}{\sqrt{1-k^2}}\right)</math>. Si Modèle:Math est la période <math>\frac{2\pi}{\omega_0}</math> du pendule simple pour les petites oscillations, alors, la période du pendule devient :
- <math>T = T_0 \, \frac{2}{\pi} \ln \left( \frac{8}{\pi-\theta_0} \right)</math>.
Quand la limite est atteinte, Modèle:Math est égale à la fonction Modèle:Math. On a alors :
- <math>\theta = 2\arcsin(\tanh(\omega_0 t)) = 4\arctan(\exp(\omega_0 t)) - \pi</math>
qui tend vers Modèle:Mvar quand Modèle:Mvar tend vers l'infini.
Le pendule simple tournoyant
Dans le cas d'un pendule animé d'une vitesse suffisamment grande pour le faire tournoyer, l'équation du mouvement s'écrit :
- <math>\dot\theta^2 = \frac{2g}{l^2}[H - l + l\cos(\theta)]</math>
où Modèle:Mvar est une constante homogène à une longueur et strictement supérieure à Modèle:Math. La solution Modèle:Mvar s'exprime alors à l'aide de la fonction amplitude de Jacobi sous la forme :
- <math>\theta = 2A\left(\sqrt{\frac{Hg}{2l^2}}\,t,k\right)</math>
où l'on donne au module de la fonction de Jacobi la valeur <math>k=\sqrt{\frac{2l}{H}}</math>.
Modèle:Démonstration \mathrm{dn}\left(\sqrt{\frac{Hg}{2l^2}}\,t\right)</math> donc, en élevant au carré, et en tenant compte du fait que <math>\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)</math> :
- <math>\begin{align}
\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{align} </math> comme attendu. }}
Le mouvement à la Poinsot d'un solide
Modèle:Voir Ce mouvement est celui d'un solide en rotation, pris relativement à son centre d'inertie Modèle:Mvar, quand le moment par rapport à Modèle:Mvar 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 à Modèle:Math, Modèle:Math, Modèle:Math<ref> Modèle:Ouvrage</ref>.
Propagation d'onde
La fonction <math>\xi : (x,t) \to \mathrm{cn}^2(B(x-ct))</math> 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, Modèle:Mvar est la hauteur de la vague, Modèle:Mvar est l'abscisse où l'on mesure cette hauteur, Modèle:Mvar est le temps, et Modèle:Mvar un paramètre prenant en compte la profondeur du milieu. C'est en effet une des solutions de l'équation de Korteweg-de Vries. L'onde ainsi modélisée s'appelle onde cnoïdale.
Pompage optique
La fonction Modèle:Math intervient pour modéliser la déplétion de la pompe dans le mélange à trois ondes optiques<ref>Modèle:Lien web, ch.5 (« Parametric amplification and generation »)</ref>, qui est utilisé dans les Oscillateurs paramétriques optiques.
Voir aussi
Bibliographie
- Modèle:Ouvrage, chapitre 16, par L. M. Milne-Thomson.
- Modèle:Ouvrage
- Modèle:Ouvrage
- Modèle:Ouvrage
Liens externes
- Modèle:Lien web. 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.
