質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

新規登録して質問してみよう
ただいま回答率
85.46%
Matplotlib

MatplotlibはPythonのおよび、NumPy用のグラフ描画ライブラリです。多くの場合、IPythonと連携して使われます。

NumPy

NumPyはPythonのプログラミング言語の科学的と数学的なコンピューティングに関する拡張モジュールです。

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

Q&A

1回答

3009閲覧

Pythonで微分方程式をラプラス変換したい

ToyGames

総合スコア1

Matplotlib

MatplotlibはPythonのおよび、NumPy用のグラフ描画ライブラリです。多くの場合、IPythonと連携して使われます。

NumPy

NumPyはPythonのプログラミング言語の科学的と数学的なコンピューティングに関する拡張モジュールです。

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

0グッド

0クリップ

投稿2021/06/27 12:54

タイトルの通りですが、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ではできないのであれば他のやり方でも構いませんので、ご教示いただけると幸いです。
よろしくお願い申し上げます。

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

guest

回答1

0

質問欄のコードをコピペして、spとsymが混在していたのでsymに統一して実行したところ、エラーもなく実行されました。

python3

1import sympy as sym 2sym.init_printing() 3sym.var('a, t, s') 4y = sym.Function('y')(t) 5eq = sym.Eq( sym.diff(y, t), -a*y ) 6ans = sym.dsolve(eq) 7print(sym.laplace_transform(eq, t, s)) 8# LaplaceTransform(Eq(Derivative(y(t), t), -a*y(t)), t, s)

投稿2021/06/27 23:49

jeanbiego

総合スコア3966

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

ToyGames

2021/07/02 14:54

ちなみにsympyのバージョンを教えていただいてもよろしいですか?
jeanbiego

2021/07/04 23:07

返答遅くなりましてすみません。 sympyは1.4でした。(python 3.7)
guest

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

まだベストアンサーが選ばれていません

会員登録して回答してみよう

アカウントをお持ちの方は

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.46%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問