Fonction elliptique de Jacobi

Un article de Wikipédia, l'encyclopédie libre.

Cet article a une forme trop académique.

La forme ressemble trop à un extrait de cours et nécessite d'être réécrit afin de correspondre aux standards de l'encyclopédie. N'hésitez pas à l'améliorer.

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, cn et dn, qui rappellent cos et sin. Si les fonctions elliptiques thêta de Weierstrass semblent mieux adaptées aux considérations théoriques, les problèmes pratiques (physiques) font plus appel aux fonctions de Jacobi, ce qui en justifie l'étude.

Rappelons qu'une fonction elliptique est une fonction méromorphe définie sur le plan complexe, uniforme, non-constante et doublement périodique (une période est réelle, l'autre imaginaire).

[modifier] Introduction

Construction auxiliaire : le rectangle de Jacobi.

Il existe 12 fonctions elliptiques de Jacobi. Chacune correspond à une flèche reliant un sommet du rectangle à un autre. Ces sommets sont, par convention, notés s, c, d et n.

Le rectangle peut être considéré dans le plan complexe, en placant s à l'origine, c à un point K sur l'axe des réels, d à un point K + iK' et n à un point iK' sur l'axe imaginaire. Les nombres K et K' sont appelés « quarts de période ».

Les 12 fonctions elliptiques de Jacobi sont alors sc, sd, sn, cs, cd, cn, ds, dc, dn, ns, nc, et nd.

Plus formellement, les fonctions elliptiques de Jacobi sont les seules fonctions doublement périodiques et méromorphes vérifiant les propriétés suivantes :

Il existe un zéro simple au coin p, et un pôle simple au coin q.
La distance de p à q est égale à la moitié de la période de la fonction pq u; c'est-à-dire que la fonction pq u est périodique dans la direction pq, de période double de la distance de p à q. Il en est de même dans les deux autres directions.
Si la fonction pq u est étendue à l'un des sommets, le terme principal de l'expansion est multiplié par un coefficient 1.

Dans le cas général, on peut travailler à partir d'un parallélogramme. Cependant, si K et iK' restent sur l'axe réel et imaginaire respectivement, alors la fonction elliptique de Jacobi pq u sera réelle si u est réel.

[modifier] K(k) : intégrale elliptique complète de première espèce

Soit k positif inférieur à l'unité, appelé « module ». Soit $k' = \sqrt{1-k^2}$ , appelé « comodule ».

Soit la fonction impaire croissante u = F(A, k) :

$F(A,k) = \int_0^{A} \frac{dx}{\sqrt{1-k^2\sin^2x}}$

c'est l' « intégrale elliptique de première espèce ».

Soit A(u,k) la fonction réciproque de cette fonction croissante u : c'est l'« amplitude » de u. Il s'agit d'une fonction impaire et croissante, qui augmente de π lorsque u augmente de 2K(k), avec :

$K(k) = u(\pi/2,k) = \int_0^1\frac{dx}{\sqrt{(1-x^2)(1-k^2x^2)}}$

K(k) est l'« intégrale elliptique complète de première espèce de Legendre ». Elle intervient dans nombre de problèmes de physique mathématique.

pour les petites valeurs de k, K(k) = $\frac{\pi}{2} (1 + k^2(\frac{1}{2})^2(1 + k^2(\frac{ 3}{ 4})^2 (1 + )))$ .
quand k tend vers 1, K tend vers l'infini, comme 4 ln (4/k')
$K(k) = \frac{\pi}{2}\frac{1}{AGM(1,k')}$ , avec AGM(a,b) la moyenne arithmético-géométrique de (a,b) : on peut faire efficacement des calculs numériques.

Une table des valeurs de K(k) pourra être trouvée dans l'article « pendule pesant ».

[modifier] Signification de K

Il est pratique , mnémotechniquement, de retenir que le module k est $\sin\frac{\theta_0}{2}$ dans le cas du pendule simple, et que K joue le rôle de $π$ /2, le « K-Quart de tour », si bien que « l'amplitude du temps » A(t,k), fonction croissante de t , joue le rôle d'« échelle de temps » (voir mesure en physique) adaptée au problème : à chaque période de temps 4K, le « temps en radians » aura augmenté de 2 $π$ , c’est-à-dire d'un aller-retour du pendule.

L'anisochronicité du mouvement est patente, puisque 4K = 4K(k) dépend du module k, donc de $θ 0$ et devient infinie comme 4ln(4/k'); c'est-à-dire en amplitude : T = To . $\frac{2}{\pi}(\ln 8-\ln (\pi-\theta_0))$ .

[modifier] Remarques d'ordre général

Les tables de K(k) ne sont établies que pour k<0,5 car on donne simultanément la valeur de K'(k) = K(k'), et le « nome » q(k) = exp(- $π$ K'/K). Il existe des cas où le nome se calcule exactement à l'aide de la fonction gamma.

La variation de K(k) est tout à fait inhabituelle : très plate au début (K' = K = 1,18 pour k = 0,707, soit une corde BB' = 21) elle n'augmente vers l'infini logarithmiquement que pour des angles très proches de 180°, si bien que l'on manque expérimentalement le phénomène si l'on n'y est pas attentif.

La donnée d'un point de cette courbe est une indication utile : pour K = 5,43, dont le rapport à K(0) est 3,45, theta = 178° ( k'~ $π$ /180) ! On peut aussi retenir des rapports particuliers : K'/K = √3, pour une corde BB'= 1,l; K'/K = √7 pour une corde BB' = 1/8.

[modifier] Les trois fonctions basiques de Jacobi (1827)

la fonction sinus de Jacobi : sn (u,k) = sin A(u,k), de période 4K
la fonction cosinus de Jacobi : cn (u,k) = cos A(u,k), de période 4K
la fonction dn de Jacobi : dn (u,k) = √(1- k^2.sn^2), de période 2K

Si k = 0, on retrouve la trigonométrie ordinaire. Si k = 1,

sn(u,1) = th(u) ;
cn(u,1) = 1/ch(u) ;
dn(u,1) = 1/ch(u) ;
tan(A/2) = th (u/2) (formule de Gudermann).

Gudermann(1838), puis Glaisher(1882) introduiront les 9 autres fonctions avec les notations « modernes ».

[modifier] Fonctions réciproques

Arcsn(x,k) = $\int_0^x \frac{dx}{\sqrt(1-x^2)(1-k^2x^2)}$
Arccn(x,k) = $\int_x^1 \frac{dx}{\sqrt(1-x^2)(1-k^2+k^2x^2)}$

[modifier] Dérivées

$\frac{d sn(u)}{du} = cn(u)dn(u)$

$\frac{d cn(u)}{du} = -sn(u)dn(u)$

$\frac{d dn(u)}{du} = -k^2sn(u)cn(u)$

[modifier] Relations trigonométriques

[modifier] Addition

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

[modifier] Carrés

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

[modifier] Valeurs remarquables

pour u = 0 , on a $s n = 0$ , $c n = 1$ , $d n = 1$
pour u = K , on a $s n = 1$ , $c n = 0$ , $d n = k'$
pour u = K/2, on a $sn =1/\sqrt{1+k'}$ , $cn = \sqrt{k'/(1+k')}$ , $dn = \sqrt{k'}$

[modifier] Translation de K

$sn(u + K) = cn(u)/dn(u)\,$
$cn(u + K) = -k\,' sn(u)/dn(u)$
$dn(u + K) = k'/dn(u)\,$

[modifier] Transformées des carrés en arc double

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

[modifier] Équations différentielles

Le formulaire précédent - assez succinct - permet alors de montrer que cn, dn et sn sont liées aux équations différentielles suivantes :

$\ddot{x} + (1+k^2)x -2k^2 x^3 = 0$ (ressort mou)
$\ddot{x} + (1-2k^2)x +2k^2 x^3 = 0$ (ressort dur)
$\ddot{x} - (2-k^2)x +2 x^3 = 0$ (origine instable)

[modifier] Équation différentielle du pendule

On a :

s i n (x / 2) = k s n (x, k)

; donc

c o s (x / 2) = d n (x, k)

En dérivant cette relation, on trouve :

$cos(x/2) \dot{x} = 2k cn(x,k)dn(x,k)$ : c'est bien $\dot{x}$ .

En dérivant à nouveau :

$\ddot{x}= -2k sn(x,k)dn(x,k)= -2sin(x/2)cos(x/2) = -sin(x)$ .

[modifier] Démonstration des formules du pendule simple

Partant de $v 2 + 2 g h = 2 g H = 2 g l . k 2$ , il existe deux cas:

k<1 : le pendule oscille ;
k>1 : le pendule tournoie.

[modifier] Oscillations

Soit A le point le plus bas, et B le point le plus haut. Il se trouve que la « bonne fonction inconnue » est la corde AM, rapportée à AB. On peut la trouver à partir de u(t), sa racine carrée, qui vérifie l'équation différentielle :

$\dot{u^2} = \omega_0^2 (1-u^2)(1-k^2u^2)$

la solution est :

u(t) = sn(

ω 0

t,k).

La vitesse v(t) s'obtient par dérivation de sin( $θ$ /2) :

v 2 (t) = 2 g (H - h) = 2 g H . c n 2 (ω 0 t, k)

Une autre manière de retrouver ce résultat est de remarquer que h.2l = AM².

[modifier] Tournoiements

La corde est rapportée à 2l, et sa racine carrée est posée égale à y(t). On trouve :

$\dot{y^2} = k^2\omega_0^2 (1-y^2)(1-y^2/k^2)$

La solution est :

y(t) = sn (k

ω 0

t,1/k).

La vitesse v(t) s'obtient par dérivation de sin( $θ$ /2) : $v 2 (t) = 2 g (H - h) = 2 g H . d n 2 (k ω 0 t,1 / k)$ ,

De même, une autre manière de le retrouver est de noter que h.2l = AM².

La fonction dn(t,k) ne s'annulant jamais et variant de 1 à k', v² varie bien de 2gH à 2g(H-l).

[modifier] Fascination expérimentale

En prenant deux pendules de Mach identiques, et en lançant les deux boules de telle sorte que k1 > 1 > k2, avec K(1/k1)/k1 = K(k2), si les deux plateaux sont maintenus horizontaux et éclairés par une lumière rasante de sorte à projeter sur un écran l'ombre des deux boules, on observe la différence entre cn²(t) et dn²(t).

Ces « chutes ralenties » auraient vraisemblablement ravi Galilée^{[réf. nécessaire]}.

[modifier] Autres utilisations

La fonction cn(B(x-Ct))² 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.

[modifier] Références

Greenhill, ed. Dover, 1959.

Portail des mathématiques

See also ebooksgratis.com: no banners, no cookies, totally FREE.