主題
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/

使用 探索

參考文獻

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." math-fun@cs.arizona.edu 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.

在 上引用

Lambert W 函式

引用為

Weisstein, Eric W. "Lambert W 函式。" 來自 Web 資源。 https://mathworld.tw/LambertW-Function.html

主題分類