Intégrale de Borwein

De testwiki
Aller à la navigation Aller à la recherche

Modèle:Voir homonymes

En mathématiques, une intégrale de Borwein est une intégrale mettant en jeu des produits de sinc(ax), où sinc est la fonction sinus cardinal, définie par sinc(x) = sin(x)/x. Les intégrales de Borwein, découvertes par David Borwein et Jonathan Borwein en 2001, présentent des régularités apparentes qui finissent par cesser. Ainsi,

0+sinxxdx=π20+sinxxsin(x/3)x/3dx=π20+sinxxsin(x/3)x/3sin(x/5)x/5dx=π2.

Ce schéma continue jusqu'à

0+sinxxsin(x/3)x/3sin(x/13)x/13dx=π2.

Cependant, à l'étape suivante, on a l'étrange résultat[1]

0+sinxxsin(x/3)x/3sin(x/15)x/15dx=467807924713440738696537864469935615849440640907310521750000π=π26879714958723010531935615849440640907310521750000ππ22,31×1011.

Plus généralement, des intégrales similaires ont pour valeur Modèle:Math chaque fois que les nombres 3, 5, ... sont remplacés par des réels positifs dont la somme des inverses est inférieure à 1. Dans l'exemple précédent, Modèle:Math, mais Modèle:Math.

En ajoutant un facteur supplémentaire, Modèle:Math, dans le produit, le schéma peut être prolongé : on a

0+cos(x)sin(x)xsin(x/3)x/3sin(x/111)x/111dx=π4,

mais

0+cos(x)sin(x)xsin(x/3)x/3sin(x/111)x/111sin(x/113)x/113dx<π4,

Dans ce cas, on a Modèle:Math, mais Modèle:Math

Les démonstrations de ces schémas ont été établies par des démonstrations intuitives[2]. En particulier, une reformulation en termes de marche aléatoire, couplée à un argument de causalité, éclaire le changement de comportement des intégrales de Borwein, et permet des généralisations à des familles reliées[3].

Explication du phénomène

La théorie de Fourier donne des éléments de réponse : la fonction sinus cardinal est la transformée de Fourier d'une fonction porte et la fonction intégrée est la transformée de Fourier du produit de convolution de plusieurs fonctions porte. Cette convolution multiple va entraîner un lissage de la fonction porte initiale, mais surtout une « érosion » du plateau de celle-ci au fil des convolutions, jusqu'au point où celui-ci n'est plus visible, provoquant le décrochage de la suite d'intégrales[4].

Formule générale

Pour une suite de nombres réels non nuls Modèle:Math, on peut associer une intégrale de la forme[5]

0+k=0nsin(akx)akxdx

Pour établir la formule, on devra considérer des sommes à partir des Modèle:Mvar. En particulier, si Modèle:Math est un n-uplet où chaque terme vaut Modèle:Math, alors on écrit Modèle:Math, qui est une variation de la somme alternée, et le produit Modèle:Math. Avec ces notations, l'intégrale se réécrit :

0+k=0nsin(akx)akxdx=π2a0Cn

avec Modèle:Math si a0>|a1|+|a2|++|an|, et

Cn=112nn!k=1nakγ{±1}nεγbγnsgn(bγ) sinon.

Le cas de Borwein correspond à la suite Modèle:Math.

Pour Modèle:Math on a Modèle:Math, et Modèle:Math mais Modèle:Math. Ainsi, puisque Modèle:Math, on trouve bien

0sin(x)xsin(x/3)x/3sin(x/13)x/13dx=π2

(pour n = 6, et de même pour toutes les intégrales avec n < 7), mais

0+sin(x)xsin(x/3)x/3sin(x/15)x/15dx=π2(1(31+51+71+91+111+131+1511)7267!(1/31/51/71/91/111/131/15)).

Extensions

Afin de prolonger la suite d'intégrales constante, on peut ajouter un facteur Modèle:Math dans l'intégrande[4]:

0+cos(x)sinxxdx=π40+cos(x)sinxxsin(x/3)x/3dx=π40+cos(x)sinxxsin(x/3)x/3sin(x/5)x/5dx=π40+cos(x)sinxxsin(x/3)x/3sin(x/111)x/111dx=π40+cos(x)sinxxsin(x/3)x/3sin(x/111)x/111sin(x/113)x/113dx<π4.

Notes

Modèle:Références

Liens externes

  • Modèle:Youtube, présentation du problème et résolution à un niveau abordable.

Références

Modèle:Traduction/Référence

Modèle:Portail

  1. Le site de MathOverflow mentionne ces intégrales Modèle:En comme ayant fait craindre une erreur dans un système de calcul formel, jusqu'à ce que les programmeurs réalisent que le résultat était correct.
  2. Modèle:Article
  3. Modèle:Article
  4. 4,0 et 4,1 Modèle:Article
  5. Modèle:Article