トップQs
タイムライン
チャット
視点

ランベルトのW関数

ウィキペディアから

ランベルトのW関数
Remove ads

ランベルトのW函数(ランベルトのWかんすう、: Lambert W function)あるいはオメガ函数 (ω function)、対数積(product logarithm; 乗積対数)は、複素函数 f(z) = zez逆関係分枝として得られる函数 W の総称である。ここで、ez複素指数函数Wz= f(W(z)) を満たすことに注意すれば、任意の複素数 z に対する W 函数(一般には W 関係)の定義方程式

Thumb
W(x) のグラフの W > 4 および x < 6 の部分。W −1 なる上の枝を主枝 W0 と言い、W −1 なる下側の分枝を W1 という。

を得る。

函数 ƒ単射ではないから、関係 W は(0 を除いて)多価である。仮に実数値の W に注意を制限するとすれば、複素変数 z は実変数 x に取り換えられ、関係の定義域は区間 x1/e に限られ、また開区間 (1/e, 0) 上で二価の函数になる。さらに制約条件として W1 を追加すれば一価函数 W0(x) が定義されて、W0(0) = 0 および W0(1/e) = 1 を得る。それと同時に、下側の枝は W1 であって、W1(x) と書かれる。これは W1(1/e) = 1 から W1(0) = まで単調減少する。

ランベルト W 関係は初等函数では表すことができない[1]。ランベルト W組合せ論において有用で、例えばの数え上げに用いられる。指数函数を含む様々な方程式(例えばプランク分布ボーズアインシュタイン分布フェルミディラック分布などの最大値)を解くのに用いられ、またy'(t) = ay(t 1) のような遅延微分方程式英語版 の解としても生じる。生化学において、また特に酵素動力学において、ミカエリスメンテン動力学の経時動力学解析に対する閉じた形の解はランベルト W 函数によって記述される。

Thumb
複素数平面におけるランベルト W 函数の主枝。負の実軸に沿った分岐切断1/e を端点に持つ。この図では、点 z における色相を W偏角で、輝度を W絶対値で決定している。
Remove ads

用語について

Thumb
二つの主枝 W0, W1

ランベルト W-函数はヨハン・ハインリヒ・ランベルトに因んで名づけられた。Digital Library of Mathematical Functions では主枝 W0Wp, 分枝 W1Wm と書いている。ここでの表記の規約(つまり W0, W1)はランベルト W に関する標準的な参考文献Corless et al. (1996)[2]に従った。

歴史

ランベルトは初め「ランベルトの超越方程式」に関連して1758年に考察した[3]。これはレオンハルト・オイラーの1783年の wew の特別な場合を論じた論文[4]に繋がる。

ランベルト W-函数は、特殊化された応用において、十年程度毎に「再発見」されてきた[要出典]。1993年には、等電荷に対する量子力学的二重井戸型ディラックデルタ函数モデル英語版(物理学における基本問題)の厳密解をランベルト W-函数が与えることが報告されたとき、コーレスら計算機代数システムMapleの開発者たちはライブラリを精査して、この函数が自然界に遍く存在することを発見した[2][5]

微分積分学

要約
視点

導函数

陰函数微分法により、W の任意の枝が常微分方程式

を満たすことが示せる(z = 1/e では W は微分できない)。従って、W の導函数は

を満たす。ここで恒等式 eW(z) = z/W(z) を用いるならば、

と書きなおすこともできる。

原始函数

函数 W(x)(およびそれを含む多くの式)は、w = W(x), (x = wew) と置いた置換積分によって

と積分できる。

したがって、(W(e) = 1 であることも考慮して)等式

が得られる。

Remove ads

漸近展開

要約
視点

W00 を中心とするテイラー級数は、逆に解いて英語版

ダランベールの収束判定法によると、収束半径1e である。この級数の定める函数は、区間 (∞, 1/e]に沿って分岐切断英語版を入れれば、ガウス平面の全域で定義される正則函数に延長することができる。この正則函数をランベルト W 函数の主値と定める。

x が十分大きければ、W0 は漸近的に

と展開される。ただし、L1 = ln(x), L2 = ln(ln(x)) であり、[k
n
 
] は非負の第一種スターリング数である[6]

もう一つの、区間 (∞, 1/e] 上で定義される実函数な枝 W1 は、L1 = ln(x), L2 = ln(ln(x)) と書けば、x が十分 0 に近いとき同じ形の漸近展開を持つ。

x e なるとき、

という上下の評価が成り立つ[7]。また もう一つの枝 W1 の評価は u > 0 に対して

となる[8]

整数冪・複素数冪の展開

W0 の整数乗もまた 0 において単純なテイラー級数(あるいはローラン級数)展開を持つ。例えば

より一般に、ラグランジュの反転公式を用いれば、r Z に対して

となることが示せる(これは一般に、位数 r のローラン級数になっている)。あるいは同じことだが、この式を W0(x)/x の冪に関するテイラー級数として

と書くことができる。これは任意の r C|x| < e1 に対して成立する。

Remove ads

特殊値

要約
視点

任意の非零代数的数 x に対して W(x)超越数になる。実際、W(x) が零ならば x も零でなければならず、また W(x) が非零代数的数ならばリンデマンワイエルシュトラスの定理により eW(x) は超越的でなければならず、従ってx = W(x)eW(x) もまた超越的でなければならない。

(オメガ定数)
Remove ads

等式

要約
視点

いくつかの等式は定義から直ちに得られる:

ここで、f(x) = x⋅ex は単射でないから、W(f(x)) = x は常には成り立たないことに注意すべきである。x < 0 かつ x ≠ -1 なる x を固定して、方程式 x⋅ex = y⋅eyy に関して二つの解を持ち、その一方はもちろん y = x である。もう一方の解は、W0 の場合 x < -1 に、W1 の場合 x ∈ (-1, 0) にある。これらを踏まえて、次の式を導くことができる。

[9]
(これは正しく枝を選べば他の n, x に対しても拡張できる)

f(ln(x)) を反転すれば

を得る。

オイラーの反復指数函数 h(x) を用いれば

を得る。

W を含む有用な積分公式がいくつか存在し、例えば以下のようなものが挙げられる:

  1. 一つ目の等式
    ガウス積分極座標で書き表すときに現れる。
  2. 三つ目の式
    は、二つ目の式で u = 1/x2 と置き換えることによって得られる。また一つ目の式はこの三つ目の式で z = tan(x)2 と置くことでも得られる。

分岐切断 (∞, 1/e] に沿う z を除けば(そのような z では以下の積分が確定しない)、ランベルト W 函数の主枝は、以下の積分

によって計算できる[10]。この二つの積分の値が等しいことは被積分函数の対称性による。

Remove ads

応用

要約
視点

指数関数を含む方程式の多くは、W関数を用いることで解くことができる。主な方針は、未知数を含む項を方程式の左辺(あるいは右辺)に寄せ、W関数で解を表現できる の形にすることである。例えば、方程式 を解くには、両辺に を掛け、を得る。

ここで、両辺のW関数をとれば、、即ち となる。

同様の方法で、xx = z の解は、

あるいは

となる。

複素数の無限回の累乗

が収束するとき、ランベルトのW関数を用いれば、その極限値を次のように表現できる。

ただし、log(z) は複素対数函数主値とする。

Remove ads

一般化

要約
視点

通常のランベルト Wx に関する

(1)

の形(ただし、a0, c, r は実定数)の「超越代数」方程式の厳密解 x = r + 1/cW(cecr/a0) を記述することができる。

ランベルト W 函数の一般化[11][12][13][14]として以下のようなものを挙げることができる:

  • 低次元における一般相対論および量子力学量子重力)(実は、両分野の繋がりは、以前には知られていなかった[15])への応用では、式 1 の右辺を今度は x の二次多項式とした
(2)
を考える。ここで、この二次多項式の根 r1, r2 は相異なる実定数とする。この方程式の解は一つの引数 x を持つ函数だが、ria0 のような項は解函数のパラメータとして働く。そのような側面で見れば、この一般化は超幾何函数メイヤーG函数英語版を作るのと似たような方式とも思えるが、これらの函数とは異なる「クラス」に属する。r1 = r2 のときは、式 2 の両辺は因数分解できて、1 に帰着されるから、解函数も通常の W 函数に還元される。式 2ディラトン場を支配する方程式を表すから、それにより不等静止質量の場合に対する 1+1-次元(空間一次元・時間一次元)における R=T英語版あるいは「直列」(lineal) 二体重力問題の計量や、一次元の不等電荷に対する量子力学的二重井戸型デルタ函数モデル英語版の固有エネルギーが導かれる。
  • 量子力学的三体問題英語版の特別の場合、具体的には(三次元)水素分子イオン[16]の固有エネルギーの解析解の場合、今度は式 1 の(あるいは式 2 の)右辺を x に関する無限次多項式の比とした
(3)
を考える。各 ri, si は実定数、x は固有エネルギーと核間距離 R の函数である。式 3 は、その特別の場合として式 1 および 2 も含めて、遅延微分方程式英語版の成す大きなクラスに関係する。

(1)式で表現される標準的な場合においても、原子・分子・光物理学[17] などの分野>からリーマン仮説[18]に対するKeiper-Li基準に至るまで、ランベルトのW函数の応用分野についての議論は十分尽くされたとは言えていない。

Remove ads

複素平面上のグラフ

数値的評価

要約
視点

W 函数はニュートン法を用いて近似することができて、w = W(z)(したがって z = wew)に対する逐次近似は

として与えられる。また、ハレー法英語版 を用いた近似

Corless et al. (1996)W の計算において与えている。

Remove ads

ソフトウェア実装

W 函数の実装は:

などがある。

関連項目

脚注

参考文献

外部リンク

Loading related searches...

Wikiwand - on

Seamless Wikipedia browsing. On steroids.

Remove ads