主题
Search

Lambert W 函数


LambertWFunction

Lambert W-函数,也称为 omega 函数,是 反函数,其原函数为:

 f(W)=We^W.
(1)

上图显示了该函数在 实轴 上的图像。Lambert W-函数的主值在 Wolfram 语言 中实现为ProductLog[z]。该函数的不同分支在 Wolfram 语言 中可用,表示为ProductLog[k, z],其中 k 是任意整数,k=0 对应于主值。虽然未记录在案,LambertW[k, z] 自动求值为ProductLog[k, z] 在 Wolfram 语言 中。

Lambert (1758) 考虑了以下方程的解

 x^a-x^b=(a-b)vx^(a+b),
(2)

现在被称为 兰伯特超越方程。1764 年,当兰伯特从苏黎世搬到柏林时,欧拉听说了兰伯特的论文。在 1770/1771 年关于一些相关级数展开优先权的私人争论之后,欧拉 (1783) 写了一篇关于 兰伯特超越方程 的论文,其中他介绍了一个特殊情况,简化为 wa^w=lx,这几乎是 W(x) 的定义,尽管欧拉建议定义的函数更像 -W(-x)。欧拉在这篇论文中考虑了级数解,并在第一段明确引用兰伯特是第一个考虑这个方程的人。

Eisenstein (1844) 考虑了无穷 幂塔 的级数

 h(z)=z^(z^(·^(·^·))),
(3)

它可以表示为闭合形式:

 h(z)=-(W(-lnz))/(lnz).
(4)

Pólya 和 Szegö (1925) 是第一个使用符号 W 表示兰伯特函数的人。

Banwell 和 Jayakumar (2000) 表明,W-函数描述了二极管中电压、电流和电阻之间的关系,Packel 和 Yuen (2004) 将 W-函数应用于存在空气阻力的弹道射弹。在统计力学、量子化学、组合数学、酶动力学、视觉生理学、薄膜工程、水文学和算法分析中也发现了其他应用 (Hayes 2005)。

LambertWReImAbs
最小值 最大值
实部
虚部 Powered by webMathematica

上图显示了复平面中的 Lambert W-函数。

LambertWRiemannSurfaces

上图显示了 W(z) 在复平面上的解析延拓的实部(左)和虚部(右)(M. Trott,私人通信)。

W(x)x>=-1/e 时,W(x) 为实数。它具有以下特殊值:

W(-1/2pi)=1/2ipi
(5)
W(-e^(-1))=-1
(6)
W(0)=0
(7)
W(1)=0.567143....
(8)

W(1)=0.567143... (OEIS A030178) 称为 omega 常数,可以被认为是指数的某种“黄金比例”,因为

 exp[-W(1)]=W(1),
(9)

给出

 ln[1/(W(1))]=W(1).
(10)

Lambert W-函数服从以下恒等式:

 W(x)+W(y)=W(xy(1/(W(x))+1/(W(y))))
(11)

(R. Corless 给 O. Marichev 的私人通信,2015 年 9 月 29 日)。

LambertWUnityReImAbs
最小值 最大值
实部
虚部 Powered by webMathematica

函数 W(ze^z)/z 在复平面中具有非常复杂的结构,但对于 R[z]>=1 以及实轴上方和下方略微扩展的区域,它简单地等于 1。

Lambert W-函数具有以下级数展开式:

W(x)=sum_(n=1)^(infty)((-1)^(n-1)n^(n-2))/((n-1)!)x^n
(12)
=x-x^2+3/2x^3-8/3x^4+(125)/(24)x^5-(54)/5x^6+(16807)/(720)x^7+....
(13)

拉格朗日反演定理 给出了等价的级数展开式:

 W(z)=sum_(n=1)^infty((-n)^(n-1))/(n!)z^n,
(14)

其中 n! 是一个 阶乘。然而,对于 实数 z>~0.4,此级数在越来越大的 值和 值之间振荡,因此不能用于实际的数值计算。

对于 z>~3,一个渐近 公式 可以产生相当精确的结果:

W(z)=lnz-lnlnz+sum_(k=0)^(infty)sum_(m=0)^(infty)c_(km)(lnlnz)^(m+1)(lnz)^(-k-m-1)
(15)
=L_1-L_2+(L_2)/(L_1)+(L_2(-2+L_2))/(2L_1^2)+(L_2(6-9L_2+2L_2^2))/(6L_1^3)+(L_2(-12+36L_2-22L_2^2+3L_2^3))/(12L_1^4)+(L_2(60-300L_2+350L_2^2-125L_2^3+12L_2^4))/(60L_1^5)+O[((L_2)/(L_1))^6],
(16)

其中

L_1=lnz
(17)
L_2=lnlnz
(18)

(Corless et al. 1996),纠正了 de Bruijn (1981) 中的一个印刷错误。Gosper (私人通信,1996 年 7 月 22 日) 的另一个展开式是 双重级数

 W(x)=a+sum_(n=0)^infty{sum_(k=0)^n(S_1(n,k))/([ln(x/a)-a]^(k-1)(n-k+1)!)}[1-(ln(x/a))/a]^n,
(19)

其中 S_1非负 第一类斯特林数a 是可以用于在分支之间进行选择的第一个近似值。对于 -1/e<=x<0,Lambert W-函数是双值的。对于 W(x)>=-1,该函数表示为 W_0(x) 或简称为 W(x),这称为 主分支。对于 W(x)<=-1,该函数表示为 W_(-1)(x)W导数

W^'(x)=1/([1+W(x)]exp[W(x)])
(20)
=(W(x))/(x[1+W(x)])
(21)

对于 x!=0。对于 z>0 时的 主分支

 ln[W(z)]=lnz-W(z).
(22)

Lambert W-函数的第 n 阶导数由下式给出:

 W^((n))(z)=(W^(n-1)(z))/(z^n[1+W(z)]^(2n-1))sum_(k=1)^na_(kn)W^k(z),
(23)

其中 a_(kn) 是数三角形

 1
-2  -1
9    8    2
-64  -79  -36  -6
625    974    622    192    24
(24)

(OEIS A042977)。这具有 指数生成函数

f(x)=(W(e^x(x+y(1+x)^2))-x)/(1+x)
(25)
=y-1/(2!)(x+2)y^2+1/(3!)(2x^2+8x+9)y^3-1/(4!)(6x^3+36x^2+79x+64)y^4+....
(26)

另请参阅

Abel 多项式, 位移常数, 兰伯特超越方程, Omega 常数, 幂塔

相关 Wolfram 站点

http://functions.wolfram.com/ElementaryFunctions/ProductLog/

使用 Wolfram|Alpha 探索

参考文献

Banwell, T. C. and Jayakumar, A. "Exact Analytical Solution for Current Flow Through Diode with Series Resistance." Electronics Lett. 36, 291-292, 2000.Barry, D. J., Culligen-Hensley, P. J.; and Barry, S. J. "Real Values of the W Function." ACM Trans. Math. Software 21, 161-171, 1995.Borwein, J. M. and Corless, R. M. "Emerging Tools for Experimental Mathematics." Amer. Math. Monthly 106, 899-909, 1999.Briggs, K. "W-ology, or, Some Exactly Solvable Growth Models." http://keithbriggs.info/W-ology.html.Briggs, K. "Graph Theory and Lambert's W function." http://keithbriggs.info/graph_theory_and_W.html.Corless, R. M. "Material on the Lambert W Function." (Ed. W. W. Küchlin). New York: ACM, pp. 197-204, 1997. http://www.apmaths.uwo.ca/~rcorless/frames/PAPERS/LambertW/.Corless, R. M.; Gonnet, G. H.; Hare, D. E. G.; Jeffrey, D. J.; and Knuth, D. E. "On the Lambert W Function." Adv. Comput. Math. 5, 329-359, 1996.Corless, R. M.; Gonnet, G. H.; Hare, D. E. G.; and Jeffrey, D. J. "Lambert's W Function in Maple." Maple Technical Newsletter 9, 12-22, Spring 1993.Corless, R. M. and Jeffrey, D. J. "The Wright omega Function." In Artificial Intelligence, Automated Reasoning, and Symbolic Computation (Ed. J. Calmet, B. Benhamou, O. Caprotti, L. Henocque and V. Sorge). Berlin: Springer-Verlag, pp. 76-89, 2002.Corless, R. M.; Jeffrey, D. J.; and Knuth, D. E. "A Sequence of Series for the Lambert W Function." In Proceedings of the 1997 International Symposium on Symbolic and Algebraic Computation, Maui, Hawaii. New York: ACM Press, pp. 197-204, 1997.de Bruijn, N. G. Asymptotic Methods in Analysis. New York: Dover, pp. 27-28, 1981.Eisenstein, G. "Entwicklung von alpha^(alpha^(alpha^...))." J. reine angew. Math. 28, 49-52, 1844.Euler, L. "De serie Lambertina Plurimisque eius insignibus proprietatibus." Acta Acad. Scient. Petropol. 2, 29-51, 1783. Reprinted in Euler, L. Opera Omnia, Series Prima, Vol. 6: Commentationes Algebraicae. Leipzig, Germany: Teubner, pp. 350-369, 1921.Fritsch, F. N.; Shafer, R. E.; and Crowley, W. P. "Algorithm 443: Solution of the Transcendental Equation we^w=x." Comm. ACM 16, 123-124, 1973.Gosper, R. W. Jr. "The Solutions of ye^(y^2)=x and ye^y=x." ACM SIGSAM Bull. 32, 8-10, 1998.Gosper, R. W. "Re: 2nd Order Eulerians." [email protected] posting, July 22, 1996.Gray, J. J. and Tilling, L. "Johann Heinrich Lambert, Mathematician and Scientist 1728-1777." Historia Math. 5, 13-41, 1978.Hayes, B. "Why W?" Amer. Sci. 93, 104-108, 2005.Jeffrey, D. J.; Hare, D. E. G.; and Corless, R. M. "Unwinding the Branches of the Lambert W Function." Math. Scientist 21, 1-7, 1996.Jeffrey, D. J.; Hare, D. E. G.; and Corless, R. M. "Exact Rational Solutions of a Transcendental Equation." C. R. Math. Acad. Sci. Canada 20, 71-76, 1998.Jeffrey, D. J.; Corless, R. M.; Hare, D. E. G.; and Knuth, D. E. "Sur l'inversion de y^ae^y au moyen des nombres de Stirling associes." Comptes Rendus Acad. Sci. Paris 320, 1449-1452, 1995.Kalman, D. "A Generalized Logarithm for Exponential-Linear Equations." College Math. J. 32, 2-14, 2001.Lambert, J. H. "Observations variae in Mathesin Puram." Acta Helvitica, physico-mathematico-anatomico-botanico-medica 3, 128-168, 1758.ORCCA. "The Lambert W Function." http://www.orcca.on.ca/LambertW.Packel, E. and Yuen, D. "Projectile Motion with Resistance and the Lambert W Function." College Math. J. 35, 337-350, 2004.Pólya, G. and Szegö, G. Aufgaben und Lehrsätze der Analysis. Berlin, 1925. Reprinted as Problems and Theorems in Analysis I. Berlin: Springer-Verlag, 1998.Sloane, N. J. A. Sequences A030178 and A042977 in "The On-Line Encyclopedia of Integer Sequences.""Time for a New Elementary Function?" FOCUS: Newsletter Math. Assoc. Amer. 20, 2, Feb. 2000.Valluri, S. R.; Jeffrey, D. J.; and Corless, R. M. "Some Applications of the Lambert W Function to Physics." Canad. J. Phys. 78, 823-831, 2000.Wright, E. M. "Solution of the Equation ze^z=a." Bull. Amer. Math. Soc. 65, 89-93, 1959.

在 Wolfram|Alpha 上引用

Lambert W 函数

引用为

Weisstein, Eric W. "Lambert W 函数。" 来自 MathWorld--Wolfram Web 资源。 https://mathworld.net.cn/LambertW-Function.html

主题分类