ラプラス変換・逆ラプラス変換をパソコンで自動的に計算する方法 - Wolframalpha編
mathはじめに
今回は理系の人たちにおなじみのラプラス変換・逆ラプラス変換の面倒くさい計算を、ブラウザで無料で計算できてしまう魔法のツール、Wolframalphaを紹介します。
数学の教科書の演習問題をやっていると、奇数番号しか解答がなかったり、そもそも全く解答が用意されていなかったりすることがありますよね。まあ考え方が重要であって答えが合ってるかどうかはどうでもいいともいえるけど、どうしても答えが合ってるのか気になる!という時によくお世話になっています。
以前にも記事で紹介しましたが、このWolframalphaというのは、Mathematicaをつくっている会社がやっている、ある事実についての質問に対して直接答えを返す質問応答システムが使えるWEBサービスです。 検索ボックスに数式を入れることで答えと考えられるものを提示し、グラフが描けそうなときは描いてくれたりします。
では早速計算してみましょう。
計算する
まずはラプラス変換のリファレンスを見てみましょう。
リファレンスによると、
- ラプラス変換は・・・
LaplaceTransform[expr,t,s]
- 逆ラプラス変換は・・・
InverseLaplaceTransform[expr,s,t]
とあります。最初の引数”expr”に変換したい式、”t”には変換したい式の変数(時間なのでtにしときましょう)、”s”にはとりあえずsを入れときましょう。
ラプラス変換
まずはラプラス変換の定番の以下を試してみましょう。
$$\mathcal{L}[exp(at)]= \frac{1}{s-a}$$
こんな風に入力します。
LaplaceTransform[exp(at),t,s]
すると以下のように確認でき、グラフが描かれたり勝手にいろいろやってくれます。
どんなわけわからん数式にも頑張って対応してくれます。
逆ラプラス変換
次はみんな大好き逆ラプラス変換です。
手で計算するの厳しいですよね・・・。
$$ \mathcal{L}^{-1} \left [ \frac{10s^2+8}{s^5+41s^3+13s} \right ]= ?$$
吐き気がしそうな式ですが、このコードで一発です。
InverseLaplaceTransform[ (10s^2+8)/(s^5+41s^3+13s) ,s,t]
これを手計算する機会はほぼないと思いますが大体のパターンは計算できますし、いくつかグラフも出力されるのでとても便利です。
おわりに
ほかにも色々な計算が出来るので、もしお手元に解答のない問題集とかで困っていて答えが知りたくなったら、是非WolframのDocumentationを是非見てみてください。だいたいはありますよ。