私は(f(t) * g(t))^a
のような簡単な方程式を持っています。ここで、a
はパラメータであり、f
とg
はtの関数です。私は複製しようとしている方法はf(t), g(t), f'(t)
と表現、および `G '(t)となるべき、t
に対する式を差別シンピはなぜ値を正しく置き換えないのですか?
です。アップ上記の単純な例では、結果は今
a * (f(t) * g(t))**(a - 1) * (f'(t) * g(t) + f(t) * g'(t))
する必要があり、我々は、この特定の問題のいくつかの知識一つの特定の定常値のみで、私たちが知っている(経済学の問題を)、使用します
f(t)
およびg(t)
の値である。f(tss) = 1
とg(tss) = 100
としましょう。ここでtss
は定常状態の値です。これはtss = 7
と任意に設定されます。これらは、であり、fおよびgの一般的な機能形式ではありません。です。これらの値を代入すると、
f'(tss)
とg'(tss)
という2つの未知数の式があります。この時点で、それが派生かどうかは関係ありません。彼らは単なる未知数であり、これと組み合わせると、scipy.optimize.fsolve
またはsympyのソルバーの1つを使って解くことができる非線形システムを私に与えるという別の式があります。
質問は、私はステップ1と2に固執しています。下記のコードは正しく値を置き換えていないようです。
from sympy import *
t = symbols('t')
a = symbols('a')
f, g = symbols('f g', cls=Function)
eq = (f(t) * g(t))**a
eq_diff = eq.diff(t)
output = eq_diff.evalf(subs={f:1, g:100, a:0.5})
output
を出力します。私は間違って何をしていますか?
この場合も、これはほんの数学的な例ですが、問題をうまく示しています。
「d」はどこから来たのですか? – user2357112
@ user2357112 Typo。それを私が直した。値はまだ正しく置き換えられません。 –