タイトルの通りですが、Pythonで微分方程式を定義してラプラス変換をしたいと考えています。
sympyを使って、微分方程式を定義して解くことまではできていますが、その式をラプラス変換する方法が調べてもわかりませんでした。
python
1import sympy as sp 2sp.init_printing() 3sp.var('a, t, s') 4y = sp.Function('y')(t) 5eq = sp.Eq( sp.diff(y, t), -a*y ) 6ans = sp.dsolve(eq) 7print(sym.latex(eq)) 8print(eq) 9ans 10
出力:
Eq(Derivative(y(x), x), -a*y(x))
これをラプラス変換するにはsympyのlaplase_transformを使えば良いそうですが、実際使ってみると、
python
1sym.laplace_transform(eq, t, s) 2
出力:
AttributeError: 'Equality' object has no attribute 'as_coeff_mul'
とエラーが出力されます。
ラプラス変換するにはどのようにすれば良いでしょうか?
もしsympyではできないのであれば他のやり方でも構いませんので、ご教示いただけると幸いです。
よろしくお願い申し上げます。
バッドをするには、ログインかつ
こちらの条件を満たす必要があります。
2021/07/02 14:54
2021/07/04 23:07