3

現在、割り当ての一部として、MATLABで指数方程式を解こうとしています。方程式が

exp(-t)+t*exp(-t)-n=0

1つはゼロより大きく、もう1つは小さい2つのソリューションがあります。

ただし、solve関数だけを使用すると、MATLABはlambertw関数と呼ばれるものを返し、ゼロ未満の解に対してのみeval（）を返すことができます。これは、たまたま私が答えたいものではありません。誰かがここで私を助けてもらえますか？

すべての回答とコメントを事前に感謝します！

ps別の方法として、ニュートンラプソン法を使用してそれを解決することを考えていますが、solve（）と比較して速度はどうですか？

Ziyao Wei

4

3 に答える 3

2

以下のコードでは、(定数) の方程式を数値的に解いてn=0.5いますが、選択した他の値についても同様である必要があります。

SOLVE 関数が最初に見つかった解のみを返したことに注意してください。そのため、MuPAD エンジンを直接呼び出し、解を検索する間隔を毎回指定しています。

%# lets plot the function: f(x) = exp(-x)+x*exp(-x)
h(1) = ezplot('0.5', [-1.5 10]); hold on
h(2) = ezplot('exp(-x)+x.*exp(-x)', [-1.5 10]);
set(h(1), 'LineStyle',':', 'Color','r')
legend(h, 'y = 0.5', 'y = exp(-x)+x.*exp(-x)')

%# The numeric solver only returns the first solution that it finds
x = solve('exp(-x)+x*exp(-x)=0.5')
x = vpa(x)

%# we can call the MuPAD solver and give the interval where solution can be found
x1 = evalin(symengine, 'numeric::solve(exp(-x)+x*exp(-x)=0.5, x = -1..0)')
x2 = evalin(symengine, 'numeric::solve(exp(-x)+x*exp(-x)=0.5, x = 0..3)')

%# show the solutions on the plot
plot([x1 x2], 0.5, 'ro')

SOLVE によって返される解:

x =
- 1.0*lambertw(0, -1/(2*exp(1))) - 1.0
x =
-0.76803904701346556525568352607755

MuPAD 数値解:

x1 =
-0.76803904701346556525568352607755
x2 =
1.6783469900166606534128845120945

ここに画像の説明を入力

于 2011-06-29T16:55:21.787 に答える