物理を解説 ♪
ツイッター用のシェアボタン フェイスブック用のシェアボタン はてなブックマーク用のシェアボタン ライン用のシェアボタン
ツイッター用のシェアボタン フェイスブック用のシェアボタン はてなブックマーク用のシェアボタン ライン用のシェアボタン

ラプラス変換

フィギュアスケートの羽生くんもやってるって?!
作成:2013/7/20

定義など

ラプラス変換は,フーリエ変換に似ている.他に説明する機会がないのでここでやっておこう.フーリエ変換ほど物理的なイメージは明確ではないが,式の形は似ている.

数式 関数f(t)をラプラス変換するとF(s)になる.この関数F(s)のことを「関数f(t)のラプラス変換」と呼び,L\big(f(t)\big)のように表記することもある.このLというのはラプラス(Laplace)の頭文字だ.

参考までに書いておくと,フーリエ変換は次のようなものだった. 数式 (フーリエ変換には色々な流儀があって,これはその一例である.

ラプラス変換はフーリエ変換を拡張したものだと言う人もいるし,いや,ラプラス変換を拡張したものがフーリエ変換だと言う人もいるが,そんなことはあまり気にしなくても良いのではなかろうか.良く見ると色々な点で違っている.

ちなみに,ラプラスとフーリエは同じ頃に活躍した人物であるが,ラプラス変換が 1780 年頃,フーリエ変換が 1807 年頃の発表であり,両者の発想にそれほど深い関わりはないように思われる.

なぜラプラス変換の積分範囲が正の範囲に限られているかというと,e^{-x}という関数はx\to-∞で急速に増大する形になっているからである.指数関数よりも急速に 0 に近付くような関数を使わない限り,積分範囲を負の方向に広げると発散してしまうことになる.それでは限られた関数しかラプラス変換できないので,積分範囲を制限するのである.

sは複素数の範囲を考えて構わないが,面倒なので実数の範囲だけで議論している教科書もある.曖昧な点は残したくないので,この記事では複素数だとして話をしよう.

ラプラス逆変換は次のような形で与えられ,ブロムウィッチ積分と呼ばれている. 数式 要するに複素平面を縦に一直線に積分するのである.それ以外の点では「フーリエ逆変換」に似ていないこともない.しかしこのような複素積分を使う必要はほとんどないだろう.ラプラス変換には決まったパターンがあるので,わざわざブロムウィッチ積分を使わなくても逆変換を推測することは難しくないからである.


変換の具体例

では,ラプラス変換の幾つかの例を見てみよう. 数式 これよりもっと単純な,例えばf(t)=1/tという関数のラプラス変換はここに載ってないがどうなるのだろう,と思うかも知れないが,これはt\to0のところで発散してしまうためにラプラス変換が存在していないのである.このように,あらゆる関数にラプラス変換が存在しているわけではなく,割りと厳しいのである.簡単に求められるのはここに挙げられている関数くらいであるし,とりあえずはこれくらい知っていれば問題ない.

これらを求める計算は定義通りにやりさえすればいいのだが,補足しておきたいことがあるので少しだけ具体的な計算をやってみせよう.

f(t)=1については,とりあえずsが実数だとして計算すると次のようになる. 数式 3 行目の極限のところではs<0だったら発散してしまうし,s=0だったら分母が 0 になってしまうという事情でs>0に制限しなくてはならないのが分かるだろう.

ではsが複素数だとしたらどうだろうかs=u+viだとして,上と同じ計算をしてみよう. 数式 4 行目の極限の部分では,もし実数部分が負(u<0)だったなら発散してしまうし,実数部分が 0 だったならe^{ivt}が振動してしまって一定値を持たず,これも良くないだろう.それで結局,実数部分が正だという条件付きで,同じ結果が得られるのである.

このように,ラプラス変換の結果として得られる関数F(s)にはsの変域に制限がある.sが実数である場合について言えば,同様にL(e^{at})にはs>aという制限があるし,L(sinat)L(cosat)にはs>0という制限があるし,L(sinhat)L(coshat)にはs>|a|という制限がある.これらは見た目では判断できるようなものでなく,計算をしてみることで分かる.

しかしsの範囲についてはそれほど心配するようなものではない.大事なのは,ラプラス変換によって,どんな関数とどんな関数とが結びつけられているのかという関係である.

自分はf(t)の変域とF(s)の変域にまるで関係がないように見えることを奇妙に思ったものだった.実際tsには直接の関係がない.まるでf(t)F(s)とはそれぞれ別の世界であって,ラプラス変換はその二つの世界をつなぐトンネルのようなイメージだ.

F(s) の変域に制限があるのは自分にとってはあたかも関数 F(s) に欠陥があるようなイメージである.それなのに,そのような「欠陥世界」から逆変換で f(t) に正しく戻ってこれるというのが不思議な気がしていたのだった.しかしブロムウィッチ積分というのは虚軸方向の無限積分なので,s の変域の制限には邪魔されることなくラプラス逆変換の計算ができるというのである.s を実数の範囲だけで考えているとこのようなことを考えることが出来ず,よくもまぁ,うまく出来ているものだと感心したのだった.


ラプラス変換の性質

ラプラス変換が二つの関数の世界を繋ぐトンネルだというイメージを話したが,このトンネルの働きは無秩序ではない.変換前後の世界を結びつける法則がある.

一言で言えば「線形性」である.ある関数を定数倍した関数をラプラス変換すると,その関数をラプラス変換したものを定数倍した関数の世界へとたどり着く. 数式 また,ある関数と別の関数の和をラプラス変換すると,それぞれの関数をラプラス変換したものの和を取って作った関数の世界へとたどり着く. 数式 さらに幾つもの法則が見い出せる.

その一つは変数を定数倍した時にどうなるか,というもので「相似法則」と呼ばれている.言葉で言うより式で表した方が早いだろう. 数式 これは逆変換をする時にも全く似た形の法則が成り立っていて, 数式 となっている.

さらに「像の移動法則」というものがあり,元の関数にe^{at}を掛けておけば,変換後の変数をaだけずらすことができる. 数式 この証明はひときわ簡単なので書いておこう. 数式 この法則を使えば上で紹介したL(e^{at})をわざわざ計算で求めることもなくて,L(1)の結果から簡単に導くことが出来ただろう.

このような法則を使うことで,わざわざ毎回複雑な計算をせずとも,上に紹介したラプラス変換の例以外にもちょっとしたバリエーションを作ることが可能である.


微分法則

ラプラス変換の性質は他にもあって,次の「微分法則」は今回の記事で特に強調したいものである. 数式 なんと,微分された関数f'(t)をラプラス変換したものは,変換後の世界では,微分せずに変換した関数F(s)に変数sを掛けることで表現されてしまうのである.-f(0)が余分に付いている辺りがどうもスッキリしないが,それはまぁ,そういうものだから仕方ない.

実はこの法則には少し条件があって,t\to∞においてe^{-st}f(t)\to0となることが必要なのである.この法則を求める過程でその条件を使う必要があるからだ.この法則の導出には部分積分を使えば良くて,それほど難しくはないので省略しよう.

余談だが,この法則さえ知っていればL\big(e^{at}\big)=1/(s-a)という関係を導くことは簡単だった.なぜなら, 数式 であり,これと同じものは次のようにも計算できて, 数式 となり,両方の結果を等式で結べば, 数式 となるからである.先ほど話した「微分法則が使える条件」を満たすためにはe^{-st}e^{at}\to0 (t\to∞)でなければならず,そこからs-a>0でなければならないこともちゃんと出てくる.

このように,バラバラに見えていたものが色々なところで繋がっていることが見えてくる.

この微分法則を繰り返し用いることで,n階微分した時の法則も導かれる.例えば 2 階微分を考えたければ, 数式 のようになる.繰り返せばどうなるかは想像が付くだろう. 数式 微分法則があるなら,積分法則も作れそうだ.ラプラス変換後に変数sを掛ければ微分をしたものをラプラス変換したものに等しいのだから,おそらくsで割ってやったものは,積分したものをラプラス変換したものに相当するはずだ.それは次のようにして示すことができる.まず次のように置く. 数式 すると,g'(t)=f(t)という関係になっている.そこで微分法則を次のように使えば良い. 数式 ここから次の法則が言える. 数式 微分法則とは違ってこちらには余計な項が出てこなくて気持ちいいが,それは積分範囲を 0 からtまでに設定することでうまくスッキリ見せているだけである.


微分方程式への応用

関数f(t)を微分したり積分したりしたものは,ラプラス変換後の世界では変数sを掛けたり,変数sで割ったりすることで表せるのである.これは前回やった,演算子法を思い出させる.演算子法では微分演算子Dを掛けることで微分を表現し,微分演算子Dによる割り算で積分を表現していたのだった.

どれくらい似ているものなのかを確認するために,前回の演算子法の実例でやってみたのと同じ問題をラプラス変換を使って解いてみることを試してみよう.次のような微分方程式を解く問題だった. 数式 この両辺をラプラス変換してやる. 数式 この左辺はラプラス変換の線形性を利用している.ここまでの説明で使ってきた関数f(t)の代わりにy(x)を使っており,変数tだったものがxになっているだけだから混乱しないようにして欲しい.続けよう. 数式 確かにsDに置き換えれば演算子法に似た風にも見えてきた.しかし右辺は演算子法では見かけなかった形である.敢えてL(x)のままにしておいた方が良かったのだろうかそれにy(0)y'(0)が邪魔である.演算子法ではこういうものにわずらわされることもなかったのだった.できるだけ計算が楽になるように,これらは 0 だということにして続けてみよう. 数式 この右辺を部分分数に分解するのは面倒だが,高校の数学でもできることなので途中は省略しよう.次のようになる. 数式 ここまで来たら両辺を逆変換してやればいい. 数式 ところで前回求めた一般解は, 数式 というものだったから,この任意定数のところに具体的な値を代入した解が求められたことが分かる.今回得た解は途中で使った条件y(0)=0かつy'(0)=0を満たすものなのだろう. 数式 確かにそうなっている.


演算子法と等価だと言えるのか

残る疑問は,ラプラス変換を使ったこのような解法と演算子法とが本当に等価だと言えるのかどうかだ.微分を掛け算で,積分を割り算で表現して,代数的に処理しようとする辺りは非常に似通っているのだが,全く同じかと言われるとよく分からない.

1920 年頃になってからブロムウィッチやミクシンスキーなどの数学者によって厳密に証明されたというのだから,そう簡単な話ではないのだろう.ヘヴィサイドによって演算子法が発表されたのが 1890 年代なので,誰にでも簡単に説明が付くような話ならもっと早くに解決していたはずだと思われる.

言われてみれば,演算子法の公式とラプラス変換の公式には似通っているものがある.前回,衒学的だと言いながら紹介した演算子法の公式と,今回紹介した三角関数のラプラス変換はそっくりだ. 数式 それを言うなら,前回の最後に書き直してみた演算子法の基本的な公式も全く同じ規則になっており,e^{at}のラプラス変換の公式と対応が付くのである. 数式 中途半端なところで申し訳ないが,現時点で自分にできるのはこのような類似点を指摘することくらいである.



趣味の物理学書店

趣味で量子力学2の広告バナー