ランベルトのW関数
![](http://upload.wikimedia.org/wikipedia/commons/thumb/4/4f/Lambert-w.svg/288px-Lambert-w.svg.png)
ランベルトのW函数(ランベルトのWかんすう、英: Lambert W function)あるいはオメガ函数 (ω function)、対数積(product logarithm; 乗積対数)は、函数 f(z) = zez の逆関係の分枝として得られる函数 W の総称である。ここで、ez は指数函数、z は任意の複素数とする。すなわち、W は z = f−1(zez) = W(zez) を満たす。
上記の方程式で、z' = zez と置きかえれば、任意の複素数 z' に対する W 函数(一般には W 関係)の定義方程式
を得る。
函数 ƒ は単射ではないから、関係 W は(0 を除いて)多価である。仮に実数値の W に注意を制限するとすれば、複素変数 z は実変数 x に取り換えられ、関係の定義域は区間 x ≥ −1/e に限られ、また開区間 (−1/e, 0) 上で二価の函数になる。さらに制約条件として W ≥ −1 を追加すれば一価函数 W0(x) が定義されて、W0(0) = 0 および W0(−1/e) = −1 を得る。それと同時に、下側の枝は W ≤ −1 であって、W−1(x) と書かれる。これは W−1(−1/e) = −1 から W−1(−0) = −∞ まで単調減少する。
ランベルト W 関係は初等函数では表すことができない[1]。ランベルト W は組合せ論において有用で、例えば木の数え上げに用いられる。指数函数を含む様々な方程式(例えばプランク分布、ボーズ–アインシュタイン分布、フェルミ–ディラック分布などの最大値)を解くのに用いられ、またy'(t) = ay(t − 1) のような遅延微分方程式 の解としても生じる。生化学において、また特に酵素動力学において、ミカエリス–メンテン動力学の経時動力学解析に対する閉じた形の解はランベルト W 函数によって記述される。
![](http://upload.wikimedia.org/wikipedia/commons/thumb/1/19/Product_Log.jpg/288px-Product_Log.jpg)
用語について
[編集]![](http://upload.wikimedia.org/wikipedia/commons/thumb/5/5d/Diagram_of_the_real_branches_of_the_Lambert_W_function.png/220px-Diagram_of_the_real_branches_of_the_Lambert_W_function.png)
ランベルト W-函数はヨハン・ハインリヒ・ランベルトに因んで名づけられた。Digital Library of Mathematical Functions では主枝 W0 を Wp, 分枝 W−1 は Wm と書いている。ここでの表記の規約(つまり W0, W−1)はランベルト W に関する標準的な参考文献Corless et al. (1996)[2]に従った。
歴史
[編集]ランベルトは初め「ランベルトの超越方程式」に関連して1758年に考察した[3]。これはレオンハルト・オイラーの1783年の wew の特別な場合を論じた論文[4]に繋がる。
ランベルト W-函数は、特殊化された応用において、十年程度毎に「再発見」されてきた[要出典]。1993年には、等電荷に対する量子力学的二重井戸型ディラックデルタ函数モデル(物理学における基本問題)の厳密解をランベルト W-函数が与えることが報告されたとき、コーレスら計算機代数システムMapleの開発者たちはライブラリを精査して、この函数が自然界に遍く存在することを発見した[2][5]。
微分積分学
[編集]導函数
[編集]を満たすことが示せる(z = −1/e では W は微分できない)。従って、W の導函数は
を満たす。ここで恒等式 eW(z) = z/W(z) を用いるならば、
と書きなおすこともできる。
原始函数
[編集]函数 W(x)(およびそれを含む多くの式)は、w = W(x), (x = wew) と置いた置換積分によって
と積分できる。
したがって、(W(e) = 1 であることも考慮して)等式
が得られる。
漸近展開
[編集]ダランベールの収束判定法によると、収束半径は 1⁄e である。この級数の定める函数は、区間 (−∞, −1/e]に沿って分岐切断を入れれば、ガウス平面の全域で定義される正則函数に延長することができる。この正則函数をランベルト W 函数の主値と定める。
x が十分大きければ、W0 は漸近的に
と展開される。ただし、L1 = ln(x), L2 = ln(ln(x)) であり、[k
n ] は非負の第一種スターリング数である[6]。
もう一つの、区間 (−∞, −1/e] 上で定義される実函数な枝 W−1 は、L1 = ln(−x), L2 = ln(−ln(−x)) と書けば、x が十分 0 に近いとき同じ形の漸近展開を持つ。
x ≥ e なるとき、
という上下の評価が成り立つ[7]。また もう一つの枝 W−1 の評価は u > 0 に対して
となる[8]。
整数冪・複素数冪の展開
[編集]W0 の整数乗もまた 0 において単純なテイラー級数(あるいはローラン級数)展開を持つ。例えば
より一般に、ラグランジュの反転公式を用いれば、r ∈ Z に対して
となることが示せる(これは一般に、位数 r のローラン級数になっている)。あるいは同じことだが、この式を W0(x)/x の冪に関するテイラー級数として
と書くことができる。これは任意の r ∈ C と |x| < e−1 に対して成立する。
特殊値
[編集]任意の非零代数的数 x に対して W(x) は超越数になる。実際、W(x) が零ならば x も零でなければならず、また W(x) が非零代数的数ならばリンデマン–ワイエルシュトラスの定理により eW(x) は超越的でなければならず、従ってx = W(x)eW(x) もまた超越的でなければならない。
- (オメガ定数)
等式
[編集]いくつかの等式は定義から直ちに得られる:
ここで、f(x) = x⋅ex は単射でないから、W(f(x)) = x は常には成り立たないことに注意すべきである。x < 0 かつ x ≠ -1 なる x を固定して、方程式 x⋅ex = y⋅ey は y に関して二つの解を持ち、その一方はもちろん y = x である。もう一方の解は、W0 の場合 x < -1 に、W−1 の場合 x ∈ (-1, 0) にある。これらを踏まえて、次の式を導くことができる。
- [9]
- (これは正しく枝を選べば他の n, x に対しても拡張できる)
f(ln(x)) を反転すれば
を得る。
オイラーの反復指数函数 h(x) を用いれば
を得る。
W を含む有用な積分公式がいくつか存在し、例えば以下のようなものが挙げられる:
分岐切断 (−∞, 1/e] に沿う z を除けば(そのような z では以下の積分が確定しない)、ランベルト W 函数の主枝は、以下の積分
によって計算できる[10]。この二つの積分の値が等しいことは被積分函数の対称性による。
応用
[編集]指数関数を含む方程式の多くは、W関数を用いることで解くことができる。主な方針は、未知数を含む項を方程式の左辺(あるいは右辺)に寄せ、W関数で解を表現できる の形にすることである。例えば、方程式 を解くには、両辺に を掛け、を得る。
ここで、両辺のW関数をとれば、、即ち となる。
同様の方法で、xx = z の解は、
あるいは
となる。
複素数の無限回の累乗
が収束するとき、ランベルトのW関数を用いれば、その極限値を次のように表現できる。
一般化
[編集]通常のランベルト W は x に関する
の形(ただし、a0, c, r は実定数)の「超越代数」方程式の厳密解 x = r + 1/cW(ce−cr/a0) を記述することができる。
ランベルト W 函数の一般化[11][12][13][14]として以下のようなものを挙げることができる:
- を考える。ここで、この二次多項式の根 r1, r2 は相異なる実定数とする。この方程式の解は一つの引数 x を持つ函数だが、ri や a0 のような項は解函数のパラメータとして働く。そのような側面で見れば、この一般化は超幾何函数やメイヤーG函数を作るのと似たような方式とも思えるが、これらの函数とは異なる「クラス」に属する。r1 = r2 のときは、式 2 の両辺は因数分解できて、1 に帰着されるから、解函数も通常の W 函数に還元される。式 2 はディラトン場を支配する方程式を表すから、それにより不等静止質量の場合に対する 1+1-次元(空間一次元・時間一次元)における R=Tあるいは「直列」(lineal) 二体重力問題の計量や、一次元の不等電荷に対する量子力学的二重井戸型デルタ函数モデルの固有エネルギーが導かれる。
(1)式で表現される標準的な場合においても、原子・分子・光物理学[17] などの分野>からリーマン仮説[18]に対するKeiper-Li基準に至るまで、ランベルトのW函数の応用分野についての議論は十分尽くされたとは言えていない。
複素平面上のグラフ
[編集]- Plots of the Lambert W function on the complex plane
-
z = Re(W0(x + iy))
-
z = Im(W0(x + iy))
-
z = |W0(x + iy)|
-
Superimposition of the previous three plots
数値的評価
[編集]W 函数はニュートン法を用いて近似することができて、w = W(z)(したがって z = wew)に対する逐次近似は
として与えられる。また、ハレー法 を用いた近似
を Corless et al. (1996) は W の計算において与えている。
ソフトウェア実装
[編集]W 函数の実装は:
- LambertW in Maple,
lambertw
in GP (glambertW
in PARI),lambertw
in MATLAB,[19]lambertw
in octave with the 'specfun' package,lambert_w
in Maxima,[20]ProductLog
(with a silent aliasLambertW
) in Mathematica,[21]lambertw
in Python scipy's special function package[22]gsl_sf_lambert_W0
andgsl_sf_lambert_Wm1
functions in special functions section of the GNU Scientific Library - GSL.
などがある。
関連項目
[編集]脚注
[編集]- ^ Chow, Timothy Y. (1999), “What is a closed-form number?”, American Mathematical Monthly 106 (5): 440–448, doi:10.2307/2589148, MR1699262.
- ^ a b Corless, R. M.; Gonnet, G. H.; Hare, D. E. G.; Jeffrey, D. J.; Knuth, D. E. (1996). “On the Lambert W function” (PostScript). Advances in Computational Mathematics 5: 329–359. doi:10.1007/BF02124750 .
- ^ Lambert JH, "Observationes variae in mathesin puram", Acta Helveticae physico-mathematico-anatomico-botanico-medica, Band III, 128–168, 1758 (facsimile)
- ^ 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. (facsimile)
- ^ Corless, R. M.; Gonnet, G. H.; Hare, D. E. G.; Jeffrey, D. J. (1993). “Lambert's W function in Maple”. The Maple Technical Newsletter (MapleTech) 9: 12–22. CiteSeerx: 10.1.1.33.2556.
- ^ Approximation of the Lambert W function and the hyperpower function, Hoorfar, Abdolhossein; Hassani, Mehdi.
- ^ http://www.emis.de/journals/JIPAM/images/107_07_JIPAM/107_07_www.pdf
- ^ Chatzigeorgiou, I. (2013). “Bounds on the Lambert function and their Application to the Outage Analysis of User Cooperation” (PDF). IEEE Communications Letters 17 (8): 1505–1508. doi:10.1109/LCOMM.2013.070113.130972 .
- ^ Weisstein, Eric W. "Lambert W-Function". mathworld.wolfram.com (英語).
- ^ The Lambert W Function, Ontario Research Centre
- ^ Scott, T. C.; Mann, R. B.; Martinez Ii, Roberto E. (2006). “General Relativity and Quantum Mechanics: Towards a Generalization of the Lambert W Function”. AAECC (Applicable Algebra in Engineering, Communication and Computing) 17 (1): 41–47. arXiv:math-ph/0607011. doi:10.1007/s00200-006-0196-1.
- ^ Scott, T. C.; Fee, G.; Grotendorst, J. (2013). “Asymptotic series of Generalized Lambert W Function”. SIGSAM (ACM Special Interest Group in Symbolic and Algebraic Manipulation) 47 (185): 75–83. doi:10.1145/2576802.2576804 .
- ^ Scott, T. C.; Fee, G.; Grotendorst, J.; Zhang, W.Z. (2014). “Numerics of the Generalized Lambert W Function”. SIGSAM 48 (188): 42–56. doi:10.1145/2644288.2644298 .
- ^ Maignan, Aude; Scott, T. C. (2016). “Fleshing out the Generalized Lambert W Function”. SIGSAM 50 (2): 45–60. doi:10.1145/2992274.2992275.
- ^ Farrugia, P. S.; Mann, R. B.; Scott, T. C. (2007). “N-body Gravity and the Schrödinger Equation”. Class. Quantum Grav. 24 (18): 4647–4659. arXiv:gr-qc/0611144. doi:10.1088/0264-9381/24/18/006.
- ^ Scott, T. C.; Aubert-Frécon, M.; Grotendorst, J. (2006). “New Approach for the Electronic Energies of the Hydrogen Molecular Ion”. Chem. Phys. 324 (2–3): 323–338. arXiv:physics/0607081. doi:10.1016/j.chemphys.2005.10.031.
- ^ Scott, T. C.; Lüchow, A.; Bressanini, D.; Morgan, J. D. III (2007). “The Nodal Surfaces of Helium Atom Eigenfunctions”. Phys. Rev. A 75 (6): 060101. doi:10.1103/PhysRevA.75.060101.
- ^ McPhedran, R. C.; Scott, T.C.; Maignan, Aude (2023). “The Keiper-Li Criterion for the Riemann Hypothesis and Generalized Lambert Functions”. ACM Commun. Comput. Algebra 57 (3): 85-110. doi:10.1145/3637529.3637530.
- ^ lambertw - MATLAB
- ^ Maxima, a Computer Algebra System
- ^ ProductLog at WolframAlpha
- ^ [1]
参考文献
[編集]- Corless, R.; Gonnet, G.; Hare, D.; Jeffrey, D.; Knuth, Donald (1996). “On the Lambert W function”. Advances in Computational Mathematics (Berlin, New York: Springer-Verlag) 5: 329–359. doi:10.1007/BF02124750. ISSN 1019-7168
- Chapeau-Blondeau, F.; Monir, A. (2002). “Evaluation of the Lambert W Function and Application to Generation of Generalized Gaussian Noise With Exponent 1/2”. IEEE Trans. Signal Processing 50 (9) .
- Francis (2000). “Quantitative General Theory for Periodic Breathing”. Circulation 102 (18): 2214–21. doi:10.1161/01.cir.102.18.2214. PMID 11056095 . (Lambert function is used to solve delay-differential dynamics in human disease.)
- Hayes, B. (2005). “Why W?”. American Scientist 93: 104–108. doi:10.1511/2005.2.104 .
- Roy, R.; Olver, F. W. J. (2010), “Lambert W function”, in Olver, Frank W. J.; Lozier, Daniel M.; Boisvert, Ronald F. et al., NIST Handbook of Mathematical Functions, Cambridge University Press, ISBN 978-0521192255
- Stewart, Seán M. (2005). “A New Elementary Function for Our Curricula?” (PDF). Australian Senior Mathematics Journal (Australian Association of Mathematics Teachers) 19 (2): 8–26. ISSN 0819-4564. ERIC EJ720055 .
- Veberic, D., "Having Fun with Lambert W(x) Function" arXiv:1003.1628 (2010); Veberic, D. (2012). “Lambert W function for applications in physics”. Computer Physics Communications 183: 2622–2628. doi:10.1016/j.cpc.2012.07.008.
- Chatzigeorgiou, I. (2013). “Bounds on the Lambert function and their Application to the Outage Analysis of User Cooperation”. IEEE Communications Letters 17 (8): 1505–1508. doi:10.1109/LCOMM.2013.070113.130972 .
外部リンク
[編集]- National Institute of Science and Technology Digital Library - Lambert W
- Weisstein, Eric W. "Lambert W-Function". mathworld.wolfram.com (英語).
- Computing the Lambert W function
- Corless et al. Notes about Lambert W research
- GPL C++ implementation with Halley's and Fritsch's iteration.
- Special Functions of the GNU Scientific Library - GSL
- An implementation of the Lambert W function for C99