LambertのW関数 (第k解）をカシオの高精度計算サイト keisan.casio.jpにUP!
z=w*exp(w)の逆関数、w=W(z)であるLambert（ランベルト）のW関数については、何回か自作式を作ったけど、結構間違っていたりした、、、
今回、MapleやMatlabでの実装を参考に、WolframAlphaでも検証したものに作り替えたので再UP!
ちゃんと第k分岐の解まで計算できるようになった。
リンクはこちら。
説明文：
ランベルトのW関数は z=w*exp(w)の逆関数 w=W(z)です。主要解の他に第k解まで計算できます。
|
最終的に画面イメージはこんな感じに。
参考文献は、大本はこれ。
”On the Lambert W function”
Matlabのも参考にしました。
そしてWolframAlphaについては
ProductLog[k,z]
で計算できる。
例えばこんな感じ。
« キャプテンサンダーボルト（上下、文庫版）（阿部和重＆伊坂幸太郎）を読んだ。期待通りめちゃくちゃ面白かった。 | トップページ | 橿原神宮へ行ってきた。神武天皇の絵巻物が公開されていた。 »
「パソコン・インターネット」カテゴリの記事
- Wordleを100回やってみて正解率の分布を見た。1回目で当たる確率をpとして、2回目からはk倍ずつ確率が上がっていく(kp,k²p,k³p…）とすると、私の結果はp=0.0089, k=2.4になった。(2022.05.27)
- 高周波回路シミュレータQucsStudioを使ってみる（その５）マルチポート（４ポートだとs4pみたいな）のSパラメータデータをTouchstoneファイルに出力しようとするとExport to SnPでは2ポートになってしまう。Project→Import Dataからポート数を手で入れないと。(2022.05.24)
- LTspiceで高周波解析（２）伝送線路のマイクロストリップライン素子を作ったのでそれをスタブとして使ってスタブBPF（バンドパスフィルタ）のSパラメータを計算してみる。QucsStudioのツールを使ってフィルタ合成したもの。(2022.05.09)
- LTspiceで高周波解析（１）伝送線路のマイクロストリップライン素子（Microstrip LineのWheelerやHammaerstadモデル）をサブサーキットとして作ってみて、かつSパラメータ解析まで.net文をつけてやってみる。線幅、基板厚、比誘電率などをパラメータで計算できる。(2022.05.06)
- ExcelでLAMBDA関数が使えるようになった(17) Sパラメータの標準フォーマットTouchstoneが3ポート以上で並びが変なので、LAMBDA,LET,SEQUENCE,IFS,INDEXなどを組み合わせてフラットにするTouchstoneFlatten(s,n)を作った。ただし今はs3p, s4p(3ポート,4ポート）のみ。(2022.04.22)
「学問・資格」カテゴリの記事
- Wordleを100回やってみて正解率の分布を見た。1回目で当たる確率をpとして、2回目からはk倍ずつ確率が上がっていく(kp,k²p,k³p…）とすると、私の結果はp=0.0089, k=2.4になった。(2022.05.27)
- 高周波回路シミュレータQucsStudioを使ってみる（その５）マルチポート（４ポートだとs4pみたいな）のSパラメータデータをTouchstoneファイルに出力しようとするとExport to SnPでは2ポートになってしまう。Project→Import Dataからポート数を手で入れないと。(2022.05.24)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2022年5月19日 IEEE Microwave Magazineで磁性体を使わない非可逆デバイス、Microwave Journalは5Gスモールセル、THz、Movanoがミリ波血圧＆血糖値計発表、三菱の3Dプリンタ衛星アンテナ、など。(2022.05.19)
- シン・ウルトラマンをIMAXで観てきた。面白かった！物理学者の有岡君がラグランジアン,AdS,プランクブレーンとかしゃべってる！Randall–Sundrumも。マグカップは標準理論＋重力で、物理監修は予想通り橋本幸士さんでした！ブラックホールの構造もタイムリー！(2022.05.14)
- 個人の正解確率pとして多数決をとったときの正解確率が、pが0.5に近づくとどのくらいの人数で0.9を超えるか計算してみた。p=0.50005だと日本の人口くらいの多数決。（神とさざなみの密室読んで正当性確率が気になったので）(2022.05.10)
「日記・コラム・つぶやき」カテゴリの記事
- 京都府では免許の更新時、優良運転者講習がオンラインでできるようになった！早速やってみたが、PCでやろうとするとマイナンバーカードをNFCリーダで読まないとだめで断念、iPhoneはChromeはだめでSafariじゃないと駄目とかいろいろ罠が…でも最終的にはよかった。(2022.05.28)
- Wordleを100回やってみて正解率の分布を見た。1回目で当たる確率をpとして、2回目からはk倍ずつ確率が上がっていく(kp,k²p,k³p…）とすると、私の結果はp=0.0089, k=2.4になった。(2022.05.27)
- 高周波回路シミュレータQucsStudioを使ってみる（その５）マルチポート（４ポートだとs4pみたいな）のSパラメータデータをTouchstoneファイルに出力しようとするとExport to SnPでは2ポートになってしまう。Project→Import Dataからポート数を手で入れないと。(2022.05.24)
- 新型コロナウイルス、日本の陽性者数＆ワクチン接種者数総計をプロット＆中国、韓国、アメリカ、ドイツ、フランス、イギリスの陽性者数もプロット（5/22更新）どの国もじわじわと増加は続いている。マスクをしなくなったり油断してるからかな…(2022.05.22)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2022年5月19日 IEEE Microwave Magazineで磁性体を使わない非可逆デバイス、Microwave Journalは5Gスモールセル、THz、Movanoがミリ波血圧＆血糖値計発表、三菱の3Dプリンタ衛星アンテナ、など。(2022.05.19)
コメント
« キャプテンサンダーボルト（上下、文庫版）（阿部和重＆伊坂幸太郎）を読んだ。期待通りめちゃくちゃ面白かった。 | トップページ | 橿原神宮へ行ってきた。神武天皇の絵巻物が公開されていた。 »
民視眼鏡 - 媒體認證台北市士林最便宜眼鏡配到好 http://www.people-eye.com.tw/
投稿: Arnoldkasia | 2017年11月29日 (水) 23時10分