Fonction W de Lambert
En mathématiques, et plus précisément en analyse, la fonction W de Lambert, nommée ainsi d'après Jean-Henri Lambert, et parfois aussi appelée la fonction Oméga, est la réciproque de la fonction de variable complexe Modèle:Mvar définie par Modèle:Math, c'est-à-dire que pour tous nombres complexes z et w, nous avons :
Puisque la fonction Modèle:Mvar n'est pas injective, Modèle:Mvar est une fonction multivaluée ou « multiforme » qui comprend deux branches pour les valeurs réelles <math>x\geqslant -\frac{1}{\mathrm{e}}</math>. Une des branches, la branche principale, Modèle:Math peut être prolongée analytiquement en dehors de Modèle:Math. Pour tout nombre complexe Modèle:Math, on a :
La fonction Modèle:Mvar de Lambert ne peut pas être exprimée à l'aide de fonctions élémentaires.
Historique
Lambert s'est intéressé à l'équation connue sous le nom d'équation transcendante de Lambert en 1758<ref>Modèle:Article.</ref>, ce qui conduisit à une note de Leonhard Euler en 1783<ref>Modèle:Article, réimprimée dans Modèle:Ouvrage.</ref> qui discutait le cas particulier de Modèle:Math. La première description de la fonction Modèle:Mvar semble due à George Pólya et Gábor Szegő en 1925<ref> Modèle:Ouvrage.</ref>. La fonction de Lambert fut « redécouverte » tous les dix ans environ dans des applications spécialisées, mais son importance ne fut pas vraiment appréciée avant les années 1990. Lorsqu'il fut annoncé que la fonction de Lambert donnait une solution exacte aux valeurs propres de l'énergie du système quantique correspondant au modèle décrit par l'opérateur de Dirac à puits double pour le cas de charges égales — un problème physique fondamental —, Corless et d'autres développeurs du système Maple firent une recherche bibliographique et découvrirent que cette fonction apparait un peu partout dans des applications pratiques<ref>Modèle:Article.</ref>.
Branches de la « fonction » de Lambert
Si nous nous limitons aux arguments réels Modèle:Math, il existe une fonction et une seule Modèle:Math à valeurs réelles <math>\ge -1</math> telle que
c'est la branche principale de Modèle:Mvar dans ce domaine. La représentation graphique de Modèle:Math figure à droite.
On note généralement Modèle:Math l'autre branche à valeurs réelles, c'est-à-dire la branche correspondant aux arguments x tels que <math>-1/\mathrm{e}\le x<0</math>, et à valeurs <math>\le -1</math>.
Propriétés élémentaires
Expression de Modèle:Math
On a Modèle:Math, donc, si Modèle:Mvar désigne une des deux branches Modèle:Math ou Modèle:Math :
<math>\forall y \neq 0, \ \mathrm{e}^{W(y)}=\frac{y}{W(y)}.</math>
Conséquences de la définition
De l'égalité de la définition, on peut déduire :
- <math>W(-\mathrm{e}^{-1})=-1</math> (où Modèle:Mvar désigne l'une quelconque des deux branches)
- <math>W_0(x \cdot \mathrm{e}^x) = x, \quad</math> si x ≥ - 1 .
- <math>W_{-1}(x \cdot \mathrm{e}^x) = x, \quad</math> si x ≤ - 1 .
- <math>W(x) = \ln\left(\frac{x}{W(x)}\right)</math> (où Modèle:Mvar désigne l'une quelconque des deux branches et Modèle:Mvar est non nul)
- <math>\ln W_0(x) = \ln(x) - W_0(x)</math> si Modèle:Math<ref>Modèle:Mathworld</ref>
- <math>W_0\left(-\frac{\ln a}{a}\right)=-\ln a\quad \left(0 < a\leqslant \mathrm{e}\right)</math>
- <math>W_{-1}\left(-\frac{\ln a}{a}\right)=-\ln a\quad \left(a\geqslant \mathrm{e}\right)</math>
- <math>W_0\left(A{\ln A}\right)=\ln A\quad \textrm{ et } \quad \mathrm{e}^{W_0(A\ln A)}=A \quad \left(A\geqslant \frac{1}{\mathrm{e}}\right)</math>
- <math>W_{-1}\left(x{\ln x}\right)=\ln x\quad \textrm{ et } \quad \mathrm{e}^{W_{-1}(x\ln x)}=x \quad \left(0<x\leqslant \frac{1}{\mathrm{e}}\right)</math>
Valeurs particulières
Voici quelques valeurs remarquables de Modèle:Mvar, obtenues simplement en remarquant que Modèle:Math, Modèle:Math, Modèle:Math, etc. :
- <math>W_0(0)=0\,</math>
- <math>W_0(\mathrm{e})=1\,</math>
- <math>W_0\left(-\frac{1}{\mathrm{e}}\right)=W_{-1}\left(-\frac{1}{\mathrm{e}}\right) = -1\,</math>
- <math>W_0(1)=\Omega\simeq 0,56714329\dots\,\,</math> où Modèle:Math est la constante oméga
- <math>W_0(1) = \mathrm{e}^{-W_0(1)} = \ln\left(\frac{1}{W_0(1)}\right) = -\ln W_0(1)</math>
On peut obtenir de même des valeurs complexes de Modèle:Math pour certains Modèle:Math ; ainsi <math>W\left(-\frac{\pi}{2}\right)=\pm\frac{\mathrm{i}\pi}{2}</math>
Dérivée
Si Modèle:Mvar désigne une des deux branches Modèle:Math ou Modèle:Math, la formule de dérivation des bijections réciproques montre que sa dérivée est :
ce qui a pour conséquence que chacune des deux branches de Modèle:Mvar satisfait l'équation différentielle :
Cette équation est d'ailleurs à variables séparables, et ses solutions sont toutes de la forme <math>x\mapsto W_{-1}(kx)</math> (avec k ≠ 0) ou <math>x\mapsto W_{0}(kx)</math>.
Primitives
La fonction Modèle:Mvar désignant une des deux branches Modèle:Math ou Modèle:Math, beaucoup de fonctions impliquant Modèle:Mvar, peuvent être intégrées en utilisant le changement de variable Modèle:Math, i.e. Modèle:Math :
\int W(x)\,\mathrm dx = xW(x) -x +\mathrm{e}^{W(x)} +C = x \left( W(x) - 1 + \frac{1}{W(x)} \right) + C
</math>Méthodes de calcul de Modèle:Math
Par la série de Taylor
La série de Taylor de Modèle:Math au voisinage de 0 peut être obtenue par l'utilisation du théorème d'inversion de Lagrange<ref>Modèle:Article.</ref> et est donnée par
W_0 (x) = \sum_{n=1}^\infty
\frac{(-n)^{n-1}}{n!}\ x^n= x - x^2 + \frac 32 x^3 - \frac 83 x^4 + \frac{125}{24}x^5 - \dotsb</math>.Le rayon de convergence est égal à Modèle:Sfrac. Cette série peut être prolongée en une fonction holomorphe définie en tout nombre complexe n'appartenant pas à l'intervalle réel Modèle:Math ; cette fonction holomorphe est aussi appelée la branche principale de la fonction Modèle:Mvar de Lambert.
Nous déduisons de la série de Taylor l'équivalent suivant de Modèle:Math en 0 : <math>W_0(x)\,\underset0\sim\,x</math>
Comme limite d'une suite
On peut calculer Modèle:Math de manière itérative, en commençant avec une valeur initiale Modèle:Math égale à 1 et en calculant les termes de la suite
- <math>w_{n+1}=x\mathrm{e}^{-w_n}</math>.
Si cette suite converge, on voit aisément que sa limite est Modèle:Math. On démontre que c'est en effet le cas si <math>x \in [-\mathrm{e}^{-1};\mathrm{e}]</math> :
- <math>\lim_{n\to+\infty}w_n =W_0(x).</math>
Il est moins simple, mais beaucoup plus efficace, d'utiliser la méthode de Newton, partant de Modèle:Math, et posant
- <math>w_{n+1}=w_n-\frac{w_n\mathrm{e}^{w_n} - x}{(1+w_n)\mathrm{e}^{w_n}}\ ;</math>
cette suite converge (très rapidement) vers Modèle:Math pour tout Modèle:Math.
Développements asymptotiques de Modèle:Math
On a, pour Modèle:Mvar tendant vers <math>+\infty</math>, le développement asymptotique à trois termes suivant<ref>On trouvera beaucoup plus de termes de ce développement dans Modèle:MathWorld</ref> :
On a pour x tendant vers Modèle:Math, le développement asymptotique de Modèle:Math :
Développement asymptotique de Modèle:Math
On peut également obtenir un développement asymptotique pour Modèle:Math avec Modèle:Mvar tendant vers Modèle:Math :
Paramétrisation des deux branches réelles de la fonction W de Lambert
Les deux branches réelles Modèle:Math et Modèle:Math de la fonction Modèle:Mvar de Lambert peuvent s'écrire de façon paramétrée.
En effet il existe <math>\alpha\in[0,1]</math> qui permet d'écrire <math>\forall x\in \left[- \frac{1}{\mathrm{e}},0 \right]</math>:<math>\left\{ \begin{array}{rl} W_{0}(x) & = \ \frac{\alpha\ln(\alpha)}{1-\alpha} \\ W_{-1}(x) & = \ \frac{\ln(\alpha)}{1-\alpha} \\ x & = \ \alpha^{\tfrac1{1-\alpha}}\ln\left(\alpha^{\tfrac1{1-\alpha}}\right) \end{array} \right.</math><ref name=":0">Modèle:Lien web</ref>
Utilisation
Beaucoup d'équations impliquant des exponentielles peuvent être résolues par l'utilisation de la fonction Modèle:Mvar. La stratégie générale est de déplacer toutes les instances de l'inconnue d'un côté de l'équation et de le faire ressembler à Modèle:Math. À ce point, la Modèle:Citation Modèle:Mvar nous fournit les solutions :
(chaque branche différente de la Modèle:Citation Modèle:Mvar donne une solution différente).
Exemples d'applications
Équation Modèle:Math
Par exemple, pour résoudre l'équation Modèle:Math, nous divisons par <math>-\frac5{\ln 2}2^t</math> pour obtenir <math> \frac{-\ln 2}{5}= -\ln(2)t \mathrm{e}^{-\ln(2)t}.</math> La définition de la Modèle:Citation Modèle:Mvar donne alors <math>-\ln(2){t}= W\left(\frac{-\ln(2)}{5}\right)</math>, soit <math>t = -\frac{W(-\ln(2)/5)}{\ln(2)}.</math>
Comme <math>-\frac{\ln 2}5<0,</math> cette formule donne deux solutions réelles : <math>t_0= -\frac{W_0(-\ln(2)/5)}{\ln(2)}</math> et <math>t_1= -\frac{W_{-1}(-\ln(2)/5)}{\ln(2)}.</math>
Équations Modèle:Mvar et Modèle:Math
Avec la « fonction » Modèle:Mvar de Lambert, on peut résoudre des équations du type Modèle:Math (avec <math>z \geq \mathrm{e}^{-1/\mathrm{e}} </math> et <math>x\geq 1/\mathrm{e} </math>) par :
donc
et, si <math>z\not =1</math>, <math>\;x=\frac{\ln(z)}{W(\ln z)}</math>.
Les solutions de l'équation :
(avec <math> b>0 </math> et <math> a\ln b \geq -1/\mathrm{e} </math>), équivalente à <math>x^x=b^a\,</math>, sont données avec la Modèle:Citation Modèle:Mvar de Lambert :
et, si <math>a\ln b\not =0</math>, <math>\quad x = \frac{a \ln(b)}{W(a \ln(b))}</math>
.
La tétration infinie
En général, la tour de puissances infinie <math>x^{x^{x^{..}}}</math> converge si et seulement si <math>\mathrm e^{-\mathrm e}\le x\le\mathrm \mathrm{e}^{1/\mathrm e}</math>.
Si Modèle:Mvar est un nombre réel avec <math>\mathrm e^{-1}\le r\le\mathrm e</math> et Modèle:Mvar le nombre <math>x=r^{1/r}</math>, alors la limite de la suite définie par <math>u_{n+1}=x^{u_n}</math> et <math>u_0=x</math> est Modèle:Mvar :
Quand une tétration infinie <math>x^{x^{x^{\cdot^{\cdot^{\cdot}}}}}</math> converge, la fonction Modèle:Math de Lambert fournit la valeur de la limite réelle comme : <math>r=\mathrm e^{-W_0(-\ln x)}</math>
(<math>r=\frac{W_0(-\ln x)}{-\ln x}</math> si Modèle:Math).
Cela peut être étendu aux nombres complexes Modèle:Mvar avec la définition :
où Modèle:Math représente la branche principale de la fonction logarithme complexe.
Équation Modèle:Math
La bijection réciproque de <math>g : x\mapsto x+\mathrm{e}^x</math> peut être obtenue explicitement : résolvant l'équation <math>g(x)=x+\mathrm{e}^x=y,</math> on remarque d'abord qu'elle équivaut, en posant <math>X=\mathrm{e}^x,</math> à <math>X\mathrm{e}^X=\mathrm{e}^y,</math> et donc <math>X=W_0(\mathrm{e}^y),</math> soit :
Équations Modèle:Math et Modèle:Math
Résolution des équations de forme : <math>a e^x + bx +c = 0</math> avec <math>(a,b,c) \in \mathbb{R}^{*2} \times \mathbb{R}</math> et Modèle:Mvar dans <math>\mathbb{R}</math>.
On pose <math>\Delta = \frac{a}{b} \mathrm{e}^{-\frac{c}{b}}</math>, ce nombre est appelé le discriminant. Il intervient dans la détermination du nombre de solutions de l'équation.
Modèle:Théorème</math> alors l'équation admet une solution dans <math>\mathbb{R}</math>.
- Si <math>\Delta \in ] -\frac{1}{\mathrm{e}} ; 0[</math> alors l'équation admet deux solutions dans <math>\mathbb{R}</math>.
- Si <math>\Delta < -\frac{1}{\mathrm{e}}</math> alors <math>E_2</math> n'admet pas de solution dans <math>\mathbb{R}</math>.}}
Modèle:Démonstration = \frac{a}{b} \mathrm{e}^{-\frac{c}{b}} </math> (produit par <math>\frac{a}{b} \mathrm{e}^{-\frac{c}{b}}</math>)
Posons <math>X= -x - \frac{c}{b} </math> et <math>\Delta = \frac{a}{b} \mathrm{e}^{-\frac{c}{b}}</math>. On a alors : <math>X \mathrm{e}^X = \Delta</math>. On trouve alors : <math>X= W(\Delta)</math>. On obtient finalement :
La fonction Modèle:Mvar de Lambert étant multivaluée sur <math>\left]\frac{-1}{\mathrm{e}};0\right[</math>, on en déduit la quantité de solutions réelles de l'équation.
}}
À l'aide du changement de variable Modèle:Math, on en déduit la résolution des équations de la forme : <math>a \ln(x) + bx + c = 0</math> avec <math>(a,b,c) \in \mathbb{R}^{*2} \times \mathbb{R}</math> et Modèle:Mvar dans <math>\mathbb{R}^{*+}</math>. Les solutions sont alors (en n'oubliant pas que Modèle:Mvar est multivaluée) de la forme :
x = \exp\left( -{W\left( \Delta \right)} - \frac{c}{a} \right)
</math> ,où <math>{ \Delta = \frac{b}{a} \mathrm{e}^{-\frac{c}{a}}}</math>.
Équations Modèle:Math et Modèle:Math
Plus généralement, la fonction Modèle:Mvar de Lambert permet de résoudre les équations de la forme : <math>a \lambda^x + bx +c = 0</math> et <math>a \log_{\lambda}(x) + bx + c = 0</math> avec <math>\lambda > 0</math> et <math>\lambda \neq 1</math>, Modèle:Mvar dans <math>\mathbb{R}</math> et <math>(a,b,c) \in \mathbb{R}^{*2} \times \mathbb{R}</math>.
Il suffit pour cela de considérer une fonction <math>W_\lambda</math> tel que <math>X \lambda^X = Y \Longleftrightarrow X = W_\lambda(Y)</math> de répéter la démonstration ci-dessus et de considérer la formule de changement de base : <math>W_\lambda(x) = \frac{W(x \ln(\lambda))}{\ln(\lambda)}</math>. On obtient alors, avec <math>\Delta = \frac{a}{b} \lambda^{-\frac{c}{b}}</math> :
a \lambda^x +bx +c = 0 \Leftrightarrow x = - W_\lambda(\Delta) - \frac{c}{b} = - \frac{W(\Delta \ln(\lambda))}{\ln(\lambda)} - \frac{c}{b}
</math>et avec <math>\Delta^\prime = \frac{b}{a} \lambda^{-\frac{c}{a}}</math>:
a \log_\lambda(x) +bx +c = 0 \Leftrightarrow x = \lambda^{\left( -\frac{W\left( \Delta^\prime \ln(\lambda) \right)}{\ln(\lambda)} - \frac{c}{a} \right)}
</math>Il faut alors considérer le nombre <math> \Delta \ln(\lambda) </math> pour déterminer la quantité de solutions des équations.
Applications en physique
Constante de Wien
Dans la loi du déplacement de Wien : <math> \sigma_w =\mathrm{ T \cdot \lambda_{max} } = 2,898 \cdot 10^{-3}\; \mathrm{m \cdot K} </math>. La constante de Wien, noté <math> \sigma_w </math> peut être déterminée explicitement à l'aide de la fonction W de Lambert.
Elle vaut : <math> \sigma_w = \frac{1}{5 + W_0 (-5e^{-5})} \frac{\mathrm{hc}}{\mathrm{k_B}} </math>, avec Modèle:Math la constante de Planck, Modèle:Math la vitesse de la lumière dans le vide et Modèle:Math la constante de Boltzmann.
Courant dans un circuit diode-résistance
La solution pour connaître la valeur du courant dans un circuit en série de diode/résistance peut être donnée par la fonction Modèle:Mvar de Lambert. Voir la Modèle:Lien.
Diverses formules intégrales
- <math>\int_{0}^{\pi} W_0\left( 2\cot^2(x) \right)\sec^2(x)\;\mathrm dx = 4\sqrt{\pi}</math> (intégrale de Gauss en coordonnées polaires)
On obtient alors par changements de variable les égalités remarquables :
- <math>\int_{0}^{+\infty} W_0\left(\frac{1}{x^2}\right)\;\mathrm dx = \sqrt{2\pi}</math>
- <math>
\int_{0}^{+\infty} \frac{W_0(x)}{x\sqrt{x}}\mathrm dx = 2\sqrt{2\pi}. </math>
Représentations graphiques
- Représentation de la partie réelle, de la partie imaginaire et du module de la fonction W de Lambert dans le plan complexe
-
<math>\scriptstyle z=\Re(W_0(x+{\rm i}y))</math>
-
<math>\scriptstyle z=\Im(W_0(x+{\rm i}y))</math>
-
</math>
-
Les trois fonctions réunies
Généralisations
La fonction W de Lambert fournit des solutions exactes aux équations « algébriques-transcendantes » (en x) de la forme :
- <math> \mathrm{e}^{-c x} = a_0 (x-r) ~~\quad\qquad\qquad\qquad\qquad(1)</math>
ou a0, c et r sont des constantes réelles. La solution est <math> x = r + W( c \mathrm{e}^{-c r}/a_o )/c. </math>
Les généralisations de la fonction W de Lambert<ref>Modèle:Article </ref>,<ref>Modèle:Article</ref>,<ref>Modèle:Article</ref> incluent :
- un lien entre la relativité générale et la mécanique quantique (gravité quantique) en dimensions réduites<ref>Modèle:Article
</ref> où la partie de droite de l'équation (1) est maintenant un polynôme quadratique en x :
- <math> \mathrm{e}^{-c x} = a_0 (x-r_1 ) (x-r_2 ) ~~\qquad\qquad(2)</math>
- où r1 et r2 sont des constantes réelles, les racines du polynôme quadratique. Dans ce cas, la solution est une fonction avec un seul argument x mais les termes comme ri et a0 sont des paramètres de la fonction. De ce point de vue, la généralisation ressemble à la série hypergéométrique et la fonction de Meijer G mais appartient pourtant à une « classe » différente de fonctions. Quand r1 = r2, chaque côté de (2) peut être factorisé et réduit à (1) et donc la solution se réduit à celle de la fonction standard de W.
L'équation (2) est celle gouvernant le champ d'un dilaton parvenant du modèle R=T- par lequel est dérivée la métrique du système gravitationnel de deux corps dans les dimensions 1+1 (c’est-à-dire une dimension spatiale et une dimension temporelle) pour le cas des masses (au repos) inégales - ainsi que les valeurs propres de l'énergie du système quantique qui est constitué du modèle décrit par l'opérateur de Dirac à puits double pour le cas de charges inégales en une dimension.
- les solutions analytiques pour les valeurs propres de l'énergie d'un cas spécial de la version quantique du problème des trois corps, c’est-à-dire l’ion hydrogène moléculaire (en trois dimensions)<ref>Modèle:Article
</ref>.
La partie de droite de (1) (ou (2)) est maintenant un quotient de « polynômes » d'ordre infini en x :
- <math> \mathrm{e}^{-c x} = a_o \frac{\prod_{i=1}^{\infty} (x-r_i )}{ \prod_{i=1}^{\infty} (x-s_i)} \qquad \qquad\qquad(3)</math>
- où ri et si sont des constantes réelles distinctes et x est une fonction de la valeur propre de l'énergie et la distance internucléaire R. L'équation (3) avec ces cas spécialisés et exprimés dans (1) et (2) correspond à une classe considérable d'équations à délai différentiel. La « fausse dérivée » de Hardy fournit des racines exactes pour des cas spéciales de (3)<ref>Modèle:Article</ref>.
Les applications de la fonction W de Lambert dans les problèmes de la physique fondamentale ne sont pas épuisées même pour le cas standard exprimé dans (1), comme on vient de le voir dans les domaines de la physique atomique et moléculaire, ainsi qu'en optique<ref>Modèle:Article </ref>.
Notes
Références
- Modèle:Article ou là
