忍者ブログ

研究の掃溜ノオト

since 2011/2/13 知能ロボ研究の合間に思ったこととか書いてます。

[PR]

×

[PR]上記の広告は3ヶ月以上新規記事投稿のないブログに表示されています。新しい記事を書く事で広告が消えます。

ランベルトのW関数

 指数関数と多項式が混ざった方程式を解くときによく出てくる”ランベルトのW関数”。
これは
f(x)=x*e^x
の逆関数として定義されています。
反復計算による関数値の計算法はよく知られているようです。
ちなみに
x=c*e^x
の解は
W(-c)
となります。
(両辺をexpで割って-1をかければ出るはず!)

参考文献:ランベルトのW関数(wikipedia)

W(c)を求める。10^1オーダーの計算で大体収束します。速いです!
Scilabソース
**********************************************************
//calculate W(z)
z=10;
//initial value
w1=0;
 
for i=1:20
    w2=w1-(w1*exp(w1)-z)/((exp(w1)*(w1+1)-(w1+2)*(w1*exp(w1)-z)/(2*w1+2)))
    w1=w2;
end
 
w1
PR

この記事へのコメント

Vodafone絵文字 i-mode絵文字 Ezweb絵文字
管理人のみ閲覧できます
 

この記事へのトラックバック

トラックバックURL

プロフィール

HN: 相馬 豊
所属:KU
連絡先(Twitter): @i-horse
インタビューはこちら

カレンダー

04 2017/05 06
S M T W T F S
1 2 3 4 5 6
7 8 9 10 11 12 13
14 15 16 17 18 19 20
21 22 23 24 25 26 27
28 29 30 31

Twitter

アンケート

マクロミルへ登録

Google Adsence

アクセス解析

リンク

Copyright ©  -- 研究の掃溜ノオト --  All Rights Reserved

Design by CriCri / Photo by momo111 / powered by NINJA TOOLS / 忍者ブログ / [PR]