7
現代数学解説
文献あり

ランベルトのW関数の積分を解いたり, テトレーションをW関数で表したりした

339
0

ランベルトのW関数

定義

ランベルトのW関数は関数f(x)=xexの逆関数のことです. 次のように定義します.

ランベルトのW関数

関数f:R(1/e,),xxexの逆関数を
f1(x)=W0(x)(1x)
f1(x)=W1(x)(x1)
と定める. このとき, W0を主枝, W1を分枝と言う.

関数のグラフは以下のようになります.

ランベルトのW関数 ランベルトのW関数

W0の級数展開

W0の級数展開

1ex1eのとき, 次の式が成り立つ.
W0(x)=n=1(n)n1n!xn

英語版WikipediaのLagrange inversion theorem(ラグランジュの反転公式)を参考に級数展開を求めていきます.

ラグランジュの反転公式 (0中心の展開)

f(x)が点0で解析的でf(0)=0,f(0)0を満たすとき
f1(x)=n=1limt0(dn1dtn1(tf(t))n)xnn!

f(x)=xex(1x)とすると, f1(x)=W0(x)である.
W0(x)=n=1limt0(dn1dtn1ent)xnn!=n=1(n)n1n!xn
を得る. ダランベールの収束判定法を用いると
|(n1)n(n+1)!xn+1(n)n1n!xn|=|(n+1)n1nn1x|=|(1+1n)n1x|<e|x|<1
よって, 収束半径は1/eです.
また, Wolfram Mathematicaでx=±1/eのときを調べると級数は収束したので, 1/ex1/eで級数は収束します.

ランベルトのW関数の積分

1W0(W11(x))dx+10W1(W01(x))dx=π23

W1の積分をW0の積分で表す

W0(W11(x))(x1)W1(W01(x))(1x)y=xで線対称であることを示す.
x1のとき1W0(W11(x))より, W0(W11(x))W1(W01(x))に代入して
W1(W01(W0(W11(x))))=W1(W11(x))=x
W0(W11(x))(x1)W1(W01(x))(1x)の逆関数となっているから, y=xで線対称である.
よって, 次の等式を得られる.
10W1(W01(x))dx=1W0(W11(x))dx1
これを求める積分に代入して
1W0(W11(x))dx+10W1(W01(x))dx=21W0(W11(x))dx1=20W0(W11(x))dx=20W0(xex)dx
となる.

級数展開で積分を級数にする

W0(x)x=0中心のテイラー展開は
W0(x)=n=1(n)n1n!xn(|x|1e)
であり, x1のとき, 1exex<0なので
20W0(xex)dx=20n=1(n)n1n!(xex)ndx

積分と無限和の入れ替えのために, 有界であることを示したい.

優関数f
f(x)={x2ex(x2)x+2+4e2(2<x0)
と定めると, 任意のx0,NNに対して
0<|n=1N(n)n1n!(xex)n|<f(x)となる.
0f(x)dx=2x2exdx+20(x+2+4e2)dx=2+18e2<


これより, ルベーグの収束定理が使えて, 総和と積分を入れ換えると
20n=1(n)n1n!(xex)ndx=2n=10(n)n1n!(xex)ndx=2n=1(n)n1n!0(xex)ndx
となる.

残りの積分部分を計算する.

0(xex)ndx
y=nxとおくと, dx=dynより
=1n0(yn)neydy=(1)nn1+n0y(n+1)1eydy=(1)nn1+nΓ(n+1)=(1)nn!n1+n
を得る.

よって, 求めた積分値から無限和を整理すると
2n=1(n)n1n!0(xex)ndx=2n=1(n)n1n!(1)nn!n1+n
=2n=11n2=π23
となる.

また, この積分より
0W0(xex)dx=π26,10W1(xex)dx=π2612
を得る.

xの無限テトレーションxxxW関数で考える

y=xy,y=xxyの概形の比較

y=xyの概形

y=xyはランベルトのW関数を用いて, yについて解ける.
y=xylogy=ylogx(1y)log(1y)=logx
ランベルトのW関数を取ると
log(1y)=logy=W(logx)y=eW(logx)

W関数は2価の関数なので, それぞれの場合でグラフを考えると, 次のようになります.

!FORMULA[72][1114525241][0] y=xy

y=xyを構成する関数
  • y=eW0(logx):(0,e1e](0,e]
  • y=eW1(logx):(1,e1e][e,)

y=xxyの概形

y=xxyyについては解けないが, xについては解けるので, 逆関数を求めることでどんな関数で構成されるのか考えます.
y=xxylogy=xylogxylogy=xylogxy=(logxy)elogxy
ランベルトのW関数を取ると
W((logxy)elogxy)=logxy=W(ylogy)x=eW(ylogy)y

!FORMULA[83][-1309638227][0] y=xxy

y=xxyを構成する関数
  • x=eW0(ylogy)y:(0,e1e](0,)
  • x=eW1(ylogy)y:(0,1e](0,1]

赤と青の曲線が交わる共有点は(1ee,1e)です.

テトレーション関連の定理

xxxが収束する実数xの範囲とその極限値

xRのとき, 0<xe1eで収束し, その極限値は
xxx=eW0(logx)=W0(logx)logx
である.

証明は前セクションからW関数を使った表示を得ているので, そこを参照してください. Wikipediaには代入するxが複素数の場合でも収束するものを選べば, この値に収束することが書かれていますが, 収束する条件が書いてませんでした. 0<|x|e1eらへんが妥当かな?

逆関数から導かれる無限テトレーションの値

f(x)=xxx(0<xe1e)に対して
f(x1x)=x

定理1のきれいな値になるバージョンです.
f(2)=222=2
とかが有名なんじゃないですかね.

y=f(x)とすると, y=xylogy=ylogxf1(y)=x=y1y
前述より, x1xfの逆関数となっている.
f(x1x)=f(f1(x))
y=f1(x)とおくと, y=f(x)が収束する範囲にx=f(y)となるxを一つ取れる.
f(x1x)=f(f1(x))=f(y)=x

おわりに

W関数の積分はxexと置換したり, 級数展開したりすれば, 上手くいくことが多いです. テトレーションをW関数で表したものを積分とかしたいのですが, 飽きてしまったので少ないですがこれで終わりにします.

書くことがなくなって題材を探すようになれば追加するかも......

参考文献

投稿日:20241220
更新日:20241221
OptHub AI Competition

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。
バッチを贈って投稿者を応援しよう

バッチを贈ると投稿者に現金やAmazonのギフトカードが還元されます。

投稿者

▽X/Twitter▽ | TyLite(トワイライト)です。認知欲求のためにSNSとMathlogしてる。フォロバはほぼ返すと思う。

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中
  1. ランベルトのW関数
  2. 定義
  3. W0の級数展開
  4. ランベルトのW関数の積分
  5. xの無限テトレーションxxxW関数で考える
  6. y=xy,y=xxyの概形の比較
  7. テトレーション関連の定理
  8. おわりに
  9. 参考文献