ラプラス変換

出典: フリー百科事典『ウィキペディア(Wikipedia)』
移動: 案内検索
関数変換

関数解析学において、ラプラス変換(ラプラスへんかん、: Laplace transform)とは、積分で定義される関数空間の間の写像(線型作用素)の一種。関数変換。

ラプラス変換の名はピエール=シモン・ラプラスにちなむ。

ラプラス変換によりある種の微分・積分は積などの代数的な演算に置き換わるため、制御工学などにおいて時間領域の(とくに超越的な)関数を別の領域の(おもに代数的な)関数に変換することにより、計算方法の見通しを良くするための数学的な道具として用いられる。

フーリエ変換を発展させて、より実用本位で作られた計算手法である。1899年に電気技師であったオリヴァー・ヘヴィサイドが回路方程式を解くための実用的な演算子を経験則として考案して発表し、後に数学者がその演算子に対し厳密に理論的な裏付けを行った経緯がある。理論的な根拠が曖昧なままで発表されたため、この計算手法に対する懐疑的な声も多かった。この「ヘヴィサイドの演算子」の発表の後に、多くの数学者達により数学的な基盤は1780年の数学者ピエール=シモン・ラプラスの著作にある事が指摘された(この著作においてラプラス変換の公式が頻繁に現れていた)。

従って、数学の中ではかなり応用寄りの分野である。ラプラス変換の理論は微分積分、線形代数、ベクトル解析、フーリエ解析、複素解析を基盤としているため、理解するためにはそれらの分野を習得するべきである。

これと類似の解法として、より数学的な側面から作られた演算子法がある。こちらは演算子の記号を多項式に見立て、代数的に変形し、公式に基づいて特解を求める方法である。

定義[編集]

実数 t ≥ 0 について定義された関数 f (t)ラプラス変換とは

F(s) = \int_{0}^\infty f(t)\mathrm{e}^{-st}dt

で定義される s の関数 F(s) のことである。ここで s複素数であり、2 つの実数 σ, ω を用いて s = σ + と表すことができる(i虚数単位)。右辺の積分はラプラス積分 (Laplace integral) と呼ばれる。これは時間領域から複素平面への写像である。

また、c > 0 として、関数 F(s) から元の関数 f (t) を計算することを逆ラプラス変換 (inverse Laplace transform) といい、

f(t) = \lim_{p\to\infty}
 \frac{1}{2\pi i} \int_{c-ip}^{c+ip} F(s)\mathrm{e}^{st}\,ds

のように定義されている。右辺の積分はブロムウィッチ積分 (Bromwich integral) と呼ばれる。これは複素平面から時間領域への写像である。

複素平面を用いたブロムウィッチ積分の解説

これは複素積分となっている。定義通りの積分経路では計算が難しくなるが、閉曲線となるように積分経路を変更して留数を計算することにより簡単に逆ラプラス変換を求める事が可能となる。 結果を言えば複素平面上の全ての特異点の留数の総和となる。 ここで、f (t)原関数 (original function)F(s)像関数 (image function) という。

ラプラス変換の他の記述の仕方として、次のようなものもある。

F(s) = \mathcal{L}[f(t)].

同様に逆ラプラス変換は、次のようにも記述される。

f(t) = \mathcal{L}^{-1}[F(s)].

また、これらの記号を用いた写像

\begin{align}
\mathcal{L} & \colon ~ f \mapsto F,\\
\mathcal{L}^{-1} & \colon F \mapsto f
\end{align}

のことも、それぞれラプラス変換逆ラプラス変換と呼ぶ。

普通、ラプラス変換および逆ラプラス変換を行う際には変換表を参照して計算する場合が多いので、前述した定義式にしたがって計算することは少ない。だが場合によっては定義式から計算したほうが簡単なときもある。たとえば逆ラプラス変換をする際に部分分数分解をしなければならない場合、むしろブロムウィッチ積分を計算したほうが早いことも多い。

注:
ラプラス変換は、関数 f (t) にいったん eσtθ(t) を乗じてからフーリエ変換する操作であると考えることができる(ここで θ(t)ステップ関数である)。
F(s):=F(\sigma,\omega) = \int_{-\infty}^{\infty} \theta(t)f(t)\mathrm{e}^{-\sigma t} \mathrm{e}^{-i\omega t}\,dt
\overset{s=\sigma+i\omega}{=\!=}\int_0^{\infty} f(t)\mathrm{e}^{-st}\,dt.

両側ラプラス変換[編集]

両側ラプラス変換は積分区間を全実数域へと拡張したもので、以下のように定義される。

F(s)  = \mathcal{L}\left\{f(t)\right\}  =\int_{-\infty}^{\infty} \mathrm{e}^{-st} f(t)\,dt

母関数との関係[編集]

数列 an の(通常型)母関数

G(a_n;x)=\sum_n a_nx^n\,

において x = es とすると、

G(a_n;\mathrm{e}^{-s})=\sum_n a_n \mathrm{e}^{-sn}\,

となる。 ここで和を積分に変えれば

G(a_t;\mathrm{e}^{-s})=\int a_t \mathrm{e}^{-st}dt\,

となり、関数 at のラプラス変換と一致する。 この意味においてラプラス変換は母関数の「連続版」とみなすことができる。 こうした理由により、母関数とラプラス変換は同種の性質を満たすことがある。 たとえば母関数の性質

G(a_n;x)G(b_n;x)=G(a_n*b_n;x)\,

はラプラス変換の性質

\mathcal{L}[f](s)\mathcal{L}[g](s)=\mathcal{L}[f*g](s)\,

に対応する。ここで *畳み込み積。

性質[編集]

ラプラス変換と逆ラプラス変換は互いに他の逆変換である。

\mathcal{L}\mathcal{L}^{-1} =
 \mathcal{L}^{-1}\mathcal{L} = \mbox{identical}.

線型性[編集]

ラプラス変換は線型性を持ち、したがって特に重ね合わせの原理 を用いて計算することが可能である。ラプラス変換が線型性を持つとは、任意の関数 f(t), g(t) に対して

 \mathcal{L}[af(t) + bg(t)] = 
  aF(s) + bG(s)

が成り立つということである。ただし、a, bt に関係しない定数。 逆ラプラス変換も同様で、

 \mathcal{L}^{-1}[aF(s) + bG(s)] = 
  af(t) + bg(t)

が成り立つ。したがって、与えられた関数を部分分数分解できるとき、各因子がラプラス変換の表にあるものに合致すれば、その変換が求められる。

相似性[編集]

a > 0 のとき、

\mathcal{L}\left[f(at)\right]=\frac{1}{a}F\left(\frac{s}{a}\right)

が成立する。

微分式[編集]

時間 t に関する導関数のラプラス変換は多項式の差となって現れる。実際に、一階の導関数をラプラス変換すると以下のように f (0)(元の式に 0 を代入した値)が現れる。

 \mathcal{L}\left[ \frac{df(t)}{dt} \right] 
  = sF(s) - f(0)

また、二階導関数の場合は f (0) に加え、t = 0 における微分係数 f'(0) が現れる。

 \mathcal{L}\left[ \frac{d^2f(t)}{dt^2} \right] 
  = s^2F(s) - sf(0) - f'(0)

これを繰り返すと、一般の n 階の導関数のラプラス変換は以下のようになる。

 \mathcal{L}\left[ \frac{d^nf(t)}{dt^n} \right] =
  s^nF(s) - \sum_{k=0}^{n-1} s^{n-k-1}f^{(k)}(0)
= s^nF(s) - s^{n-1}f(0) - s^{n-2}f^{(1)}(0) - 
 s^{n-3}f^{(2)}(0) - \cdots - f^{(n-1)}(0).

積分式[編集]

 \mathcal{L}\left[\int_0^t f(u)du \right] = \frac{1}{s} F(s).

畳み込み[編集]

関数の畳み込みはラプラス変換で積(値ごとの積)に写される。

\mathcal{L}[f(t)*g(t)] = 
 F(s)\,G(s).

これは、H(s) = F(s)G(s) かつ

F(s) = \mathcal{L}[f(t)],\, G(s) = \mathcal{L}[g(t)]

ならば

\mathcal{L}^{-1}[H(s)] = f(t)*g(t)

と書くこともできる。

初期値の定理・最終値の定理[編集]

ラプラス変換の原関数の初期値(t = 0 での値)や最終値(t → ∞ における極限値)を表す初期値の定理 (initial value theorem) および最終値の定理 (final value theorem) と呼ばれる公式が以下のような式によって与えられる。

初期値の定理
t の関数 f (t)t = 0 で連続ならば

f(0)=\lim_{t \rightarrow 0} f(t) = \lim_{s\rightarrow \infty} sF(s)
が成り立つ。特に、f が微分可能なときは部分積分により容易に証明できる。
最終値の定理
t の関数 f (t)t → ∞ で収束するなら

f(\infty)=\lim_{t \rightarrow\infty} f(t) = \lim_{s\rightarrow 0} sF(s)\qquad s \in \Delta_0
が成り立つ。ただし、Δ0s > 0 を含む角領域である。

性質一覧表[編集]

  • 表中の凡例
u(t) : ヘビサイド関数
(f * g)(t) : fg畳み込み
f'(t) : f (t)1 階微分
f^{(n)}(t) : f (t)n 階微分
片側ラプラス変換の性質(その 1)
性質 原関数
f(t) = \mathcal{L}^{-1} \left\{ F(s) \right\}
('t' 領域 / 時間領域)
像関数
F(s) = \mathcal{L}\left\{ f(t) \right\}
('s' 領域 / 周波数領域)
備考
線形性 a f(t) + b g(t) ~ a F(s)  + b G(s) ~
相似性  f(at) ~   {1 \over a} F \left ( {s \over a} \right ) ただし、a > 0
移動  f(t - a) u(t - a) ~   \mathrm{e}^{-as} F(s) ~
 f(t + \lambda) ~  \mathrm{e}^{\lambda s}\left\{F(s) - \int_0^{\lambda}\mathrm{e}^{-st}f(t)\,dt\right\}~ 移動第 2 則
ただし、λ > 0
1微分  f'(t) ~  s F(s) - f(0) ~ ただし、ƒ1微分可能とする。
2 階微分 f''(t) \, s^{2}F(s)-sf(0)-f'(0) \, ただし、ƒ2 階微分可能とする。
n 階微分  f^{(n)}(t)  ~   s^n F(s) - s^{n - 1} f(0) - \cdots - f^{(n - 1)}(0) ~ ただし、ƒn 階微分可能とする。
積分  \int_0^t f(\tau)\, d\tau  =  (u * f)(t) ~   {1 \over s} F(s) ~
 \frac{1}{(n-1)!}\int_0^t (t-q)^{n-1}f(q)\mathrm{d}q ~  \frac{1}{s^n} F(s) ~ ただし、n ≥ 1
 \int_0^t \int_0^{\tau_{n-1}}\cdots \int_0^{\tau_1} f(\tau)\,d\tau d\tau_1 \cdots d\tau_{n-1} ~  \frac{1}{s^n}F(s) ~
畳み込み  (f * g)(t) = \int_0^t f(u)\,g(t-u)\,du~  F(s) \cdot G(s) ~
周期関数  f(t) = f(t + T) ~ {1 \over 1 - \mathrm{e}^{-Ts}} \int_0^T \mathrm{e}^{-st} f(t)\,dt f (t) は周期 T の周期関数。
片側ラプラス変換の性質(その 2)
性質 像関数
F(s) = \mathcal{L}\left\{ f(t) \right\}
s' 領域 / 周波数領域)
原関数
f(t) = \mathcal{L}^{-1} \left\{ F(s) \right\}
('t' 領域 / 時間領域)
備考
移動  F(s - a) ~  \mathrm{e}^{at} f(t) ~
1 階微分  F'(s) ~  -t f(t) ~ ただし、F1 階微分可能とする。
2 階微分 F^{\prime\prime}(s)\, t^{2}f(t) \, ただし、F2 階微分可能とする。
n 階微分  F^{(n)}(s) ~  (-t)^{n} f(t) ~ ただし、Fn 階微分可能とする。
積分  \int_s^\infty F(\sigma)\, d\sigma ~  \frac{f(t)}{t} ~
 \int_s^{\infty}\int_{\sigma_{n-1}}^{\infty}\cdots\int_{\sigma_1}^{\infty}F(\sigma)\,d\sigma d\sigma_1 \cdots d\sigma_{n-1} ~  \frac{f(t)}{t^n} ~
畳み込み  (F*G)(s) = \frac{1}{2 \pi \mathrm{i}}\int_{c-\mathrm{i}\infty}^{c+\mathrm{i}\infty} F(\sigma) \,G(s-\sigma)\mathrm{d}\sigma  f(t) \,g(t)

変換表[編集]

変換表 原関数
f(t) = \mathcal{L}^{-1} \left\{ F(s) \right\}
't' 領域 / 時間領域
像関数
F(s) = \mathcal{L}\left\{ f(t) \right\}
's' 領域 / 周波数領域
収束域
単位インパルス  \delta(t) ~  1   \mathrm{all}~ s \,
単位ステップ関数  u(t) ~  { 1 \over s }  \operatorname{Re} \{ s \} > 0 \,
ランプ関数  t \cdot u(t)~ \frac{1}{s^2}  \operatorname{Re} \{ s \} > 0 \,
n
n は整数)
{  t^n \over n! } \cdot u(t)  { 1 \over s^{n+1} }  \operatorname{Re} \{ s \} > 0 \,
 (n > -1) \,
q
q は複素数)
{  t^q \over \Gamma(q+1) } \cdot u(t)  { 1 \over s^{q+1} }  \operatorname{Re} \{ s \} > 0 \,
 (\operatorname{Re}\{q\} > -1) \,
n 乗根  \sqrt[n]{t} \cdot u(t) = t^{1/n} \cdot u(t)  \frac{1}{s^{1+1/n}} \cdot \Gamma\left(1+\frac{1}{n}\right)  \operatorname{Re} \{ s \} > 0 \,
指数減衰  \mathrm{e}^{-\alpha t} \cdot u(t)  ~  { 1 \over s+\alpha }   \operatorname{Re} \{ s \} > - \alpha ~
n 乗の指数減衰 \frac{t^{n}}{n!}\mathrm{e}^{-\alpha t} \cdot u(t) \frac{1}{(s+\alpha)^{n+1}}  \operatorname{Re} \{ s \} > - \alpha \,
理想遅延  \delta(t-\tau) ~  \mathrm{e}^{-\tau s} ~
遅延付き単位ステップ関数  u(t-\tau) ~  { \frac{1}{s} \cdot \mathrm{e}^{-\tau s} }  \operatorname{Re} \{ s \} > 0 \,
遅延付き n 乗の指数減衰 \frac{(t-\tau)^n}{n!} \mathrm{e}^{-\alpha (t-\tau)} \cdot u(t-\tau)  \frac{1}{(s+\alpha)^{n+1}} \cdot \mathrm{e}^{-\tau s}  \operatorname{Re} \{ s \} > - \alpha \,
指数関数的接近 ( 1-\mathrm{e}^{-\alpha t})  \cdot u(t) ~ \frac{\alpha}{s(s+\alpha)}   \operatorname{Re} \{ s \} > 0~
正弦関数  \sin(\omega t) \cdot u(t) ~  { \omega \over s^2 + \omega^2  }  \operatorname{Re} \{ s \} > 0  ~
余弦関数  \cos(\omega t) \cdot u(t) ~  { s \over s^2 + \omega^2  }  \operatorname{Re} \{ s \} > 0 ~
双曲線正弦関数
(ハイパボリックサイン)
 \sinh(\alpha t) \cdot u(t) ~  { \alpha \over s^2 - \alpha^2 }  \operatorname{Re} \{ s \} > | \alpha | ~
双曲線余弦関数
(ハイパボリックコサイン)
 \cosh(\alpha t) \cdot u(t) ~  { s \over s^2 - \alpha^2  }  \operatorname{Re} \{ s \} > | \alpha | ~
正弦波の指数減衰 \mathrm{e}^{-\alpha t}  \sin(\omega t) \cdot u(t) ~  { \omega \over (s+\alpha )^2 + \omega^2  }  \operatorname{Re} \{ s \} > -\alpha ~
余弦波の指数減衰 \mathrm{e}^{-\alpha t}  \cos(\omega t) \cdot u(t) ~  { s+\alpha \over (s+\alpha )^2 + \omega^2  }  \operatorname{Re} \{ s \} > -\alpha ~
自然対数  \ln \left (  { t \over t_0 } \right ) \cdot u(t)  - {1 \over s} \left[\ln(t_0 s)+\gamma \right]  \operatorname{Re} \{ s \} > 0 \,
第 1 種ベッセル関数  J_n( \omega t) \cdot u(t) \frac{ \omega^n \left(s+\sqrt{s^2+ \omega^2}\right)^{-n}}{\sqrt{s^2 + \omega^2}}  \operatorname{Re} \{ s \} > 0 \,
 (n > -1) \,
第 1 種変形ベッセル関数 I_n(\omega t) \cdot u(t)  \frac{ \omega^n \left(s+\sqrt{s^2-\omega^2}\right)^{-n}}{\sqrt{s^2-\omega^2}}  \operatorname{Re} \{ s \} > | \omega | \,
第 2 種ベッセル関数
(次数が 0 の場合)
 Y_0(\alpha t) \cdot u(t) -{2 \sinh^{-1}(s/\alpha) \over \pi \sqrt{s^2+\alpha^2}} \operatorname{Re} \{ s \} > 0 \,
第 2 種変形ベッセル関数
(次数が 0 の場合)
 K_0(\alpha t) \cdot u(t)    
誤差関数  \mathrm{erf}(t) \cdot u(t)  {\mathrm{e}^{s^2/4} \left(1 - \operatorname{erf} \left(s/2\right)\right) \over s}  \operatorname{Re} \{ s \} > 0 \,
凡例

関連項目[編集]