![](http://oshiete.xgoo.jp/images/v2/pc/qa/question_title.png?e8efa67)
No.3ベストアンサー
- 回答日時:
まだ締切られていないようなので補足します。
ANo.2の式(1)と式(2)が同じかどうかですが、計算したら同じでした。
【A = x*ln(x)*B*x の解(2通り)】
x = e^[W(A*e^B) - B] --- (1)
x = A/W(A*e^B) --- (2)
ただし-1/e≦A*e^B → -e^{-(B+1)}≦A
W(x)はLambertのW関数[1]で、x=W(x)*e^W(x) の解。
【(1)=(2)の証明】
xを使うと混乱するので、sを使ってW関数の定義を書くと
s = W(s)*e^W(s)
s=A*e^B とおくと、A*e^B = W(A*e^B)*e^W(A*e^B)
両辺をe^B (≠0)で割ると、A = W(A*e^B)*e^{W(A*e^B)-B}
A=0のときA*e^B = 0となるので、A≠0のとき両辺をW(A*e^B) で割ると、A/W(A*e^B) = e^{W(A*e^B)-B}。したがって式(1)=式(2)が成り立つ。A=0のときは、式(1)からx = 1が解となるが、式(2)では0/0となってしまうので、式(1)のほうがAがゼロかどうか気にしなくてもいいという点で使いやすいと思います。
[1] LambertのW関数 http://ja.wikipedia.org/wiki/%E3%83%A9%E3%83%B3% …
【LambertのW関数の数値計算法】
LambertのW関数の参考URLに出ている漸化式を使えば、-1/e≦xに対する関数値 W(x)を計算することができる。参考まで、ExcelVBAを使って計算する方法を示します。
Function LambertW(x As Double) As Variant
If x < -Exp(-1) Then
LambertW = ""
Exit Function
Else
If x = 0 Then
LambertW = 0
Exit Function
End If
End If
Dim W As Double, W1 As Double, eps As Double, ew As Double
eps = 10 ^ (-15)
W = 1
While Abs((W - W1) / W) > eps
W1 = W
ew = Exp(W)
W = W - (W * ew - x) / (ew * (W + 1) - (W + 2) * (W * ew - x) / (W + 1) / 2)
Wend
LambertW = W1
End Function
Excelの「ツール→マクロ→Visual Basic Editor→挿入→標準モジュール」で出たコード画面に上のプログラムを貼り付け、Excelシート上で「=LambertW()」で呼び出せば計算できます。この関数を使えば、元の方程式の解は、= exp(LambertW(A*exp(B))-B)で計算できます。
No.2
- 回答日時:
inaraです。
(解法1)
xがlnの中にあるのでx>0を考えて、x = exp(t-B)とおきます(-Bをつけているのがミソ)。
すると、与式は
A = x*ln(x) - B*x = exp(t-B)*(t-B) + B*exp(t-B) = t*exp(t-B)
となります。始めに x=exp(t-B) としたのは、 B*exp(t-B) の項を消すためです。計算を続けると、
= t*exp(t)*exp(-B) → A*exp(B) = t*exp(t)
となりますから、W関数を使うと、t = W{A*exp(B)}
したがって、x = exp(t-B) = exp[W{A*exp(B)} - B] --- (1)
(解法2)
これはストレートな方法です。
A = x*ln(x) + B*x の両辺を x で割ると
A/x = ln(x) + B = ln{x*exp(B)} → x*exp(B) = exp(A/x) → x = exp(A/x)*exp(-B)
ここまではkitakazetotaiyouさんと同じですが、その後が計算違いです。
両辺にA*exp(B)/xをかけると、A*exp(B) = A/x* exp(A/x)
W関数を使うと、A/x = W{A*exp(B)} → x = A/[W{A*exp(B)}] --- (2)
(1)と(2)は形が違いますが、私はW関数の専門家でないので、(1)=(2)が成り立つのかよく分かりません。
No.1
- 回答日時:
初等関数では解けませんが、解は x = exp(LambertW(A*exp(B))-B) です。
LambertW()は "The Lambert W function" といって、次の方程式の解です。
LambertW(x) * exp(LambertW(x)) = x
数値解を求めるには、Excelのソルバーを使ったほうが良いと思います。
ありがとうございます。
教えていただいたとおりランベルトのW関数で与式を変形したところ、
inara様の解と違う式が出てきました。
どのような変形をしたらなるのでしょうか。
途中計算を教えていただけませんか。
また、W関数の中はどう計算すればいいのですか?
wikiペディアをにある練習問題を見ると
2^t=5t
をW関数で示すと
t=W{-ln(2)/5}/{-ln(2)}
とあり、W{-ln(2)/5}の計算値が分かりません。
計算方法を教えていただけませんか。
すいませんが上の2点を教えていただけないでしょうか。
よろしくお願いします。
______________________________
以下に自分の式変形を示します。
A=xln(x)+xB
両辺をxで割って
A/x=ln(x)+B
移項して
ln(x)=A/x-B
lnをexpであらわすと
x=exp(A/x-B)
expを分けると
x=exp(A/x)exp(-B)
両辺をexp(A/x)で割ると
xexp(x/A)=exp(-B)
両辺をAで割ると
(x/A)exp(x/A)=(1/A)exp(-B)
W関数を用いると
x/A=W{(1/A)exp(-B)}
両辺にAをかけて
x=(A)W{(1/A)exp(-B)}
となりました。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 物理学 積分について(2) 2 2022/11/22 13:20
- 数学 積分計算を使った漸化式とその極限 4 2023/07/04 15:40
- 数学 微分方程式の非線形2階微分方程式が解けないので教えてください!特殊解とその見つけ方だけでもお願いしま 4 2022/11/21 23:35
- 数学 3次方程式の解で実部が正のものが存在する条件の調べ方 0 2023/03/23 15:07
- 数学 xの2次方程式x2+5x-2m+1=0が異なる二つの実数解をもつような定数mの範囲を求めたいです。 2 2022/05/27 22:05
- 数学 xの2次方程式4x2+(k-1)x+1=0がただ1つの実数解をもつような、定数kの値を求めたいです。 5 2022/05/27 22:14
- 物理学 至急お願いします。高1力学です。 添付写真の問題で、(d)まで解きすすめたのですが最後方程式を解くだ 1 2022/08/01 23:07
- 数学 不定方程式について 不定方程式を解く上で合同式を利用した解法があると思いますが、合同式を使って解けな 6 2022/10/25 16:23
- 宇宙科学・天文学・天気 AIが答えた方程式 1 2023/02/20 00:12
- 数学 連立方程式についての疑問 7 2022/06/19 19:48
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
高校数学の計算(指数)
-
よく解けない漸化式の問題で絶...
-
次の微分方程式の特殊解を求め...
-
指数関数の展開方法を教えて下さい
-
常微分方程式の問題
-
外イキはなぜ1回しか出来ないの...
-
X2乗+Xの解き方について
-
内田伏一著 集合と位相 例19.2
-
√32の解き方を教えてください。...
-
再質問 5x+7y(x,yは自然数)の...
-
数学2の問題で、x2024乗をx2乗+...
-
「区分的に連続」と「区分的に...
-
数学
-
因数分解
-
1kgの10%は?
-
数学の問題について 1個200円の...
-
3次元座標から体積を求める参考書
-
X二乗の解き方。
-
写真の赤線のところがわかりま...
-
三次関数、四次関数の概形について
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
おすすめ情報