ラプラス変換

関数変換

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

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

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

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

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

定義

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

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

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

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

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

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

ファイル:ブロムウィッチ積分.png

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

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

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

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

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

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

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

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

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

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

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

両側ラプラス変換

"「両側ラプラス変換」"

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

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

母関数との関係

数列 $a n$ の（通常型）母関数

[math]G(a_n;x)=\sum_n a_nx^n\,[/math]

において $x = e - s$ とすると、

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

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

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

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

[math]G(a_n;x)G(b_n;x)=G(a_n*b_n;x)\,[/math]

はラプラス変換の性質

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

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

性質

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

[math]\mathcal{L}\mathcal{L}^{-1} = \mathcal{L}^{-1}\mathcal{L} = I [/math]

ここで、 $I$ は恒等変換を表わす。

線型性

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

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

が成り立つということである。ただし、 $a, b$ は $t$ に関係しない定数。逆ラプラス変換も同様に線形性を持ち、

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

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

相似性

$a > 0$ のとき、

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

が成立する。

微分式

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

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

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

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

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

[math] \mathcal{L}\left[ \frac{\mathrm d^nf(t)}{\mathrm dt^n} \right] = s^nF(s) - \sum_{k=0}^{n-1} s^{n-k-1}f^{(k)}(0) [/math]
[math]= 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) [/math]

積分式

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

畳み込み

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

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

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

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

ならば

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

と書くこともできる。

初期値の定理・最終値の定理

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

初期値の定理: $t$ の関数 $f (t)$ が $t = 0$ で連続ならば; [math] f(0)=\lim_{t \rightarrow 0} f(t) = \lim_{s\rightarrow \infty} sF(s) [/math]; が成り立つ。特に、 $f$ が微分可能なときは部分積分により容易に証明できる。
最終値の定理: $t$ の関数 $f (t)$ が $t \to \infty$ で収束するなら; [math] f(\infty)=\lim_{t \rightarrow\infty} f(t) = \lim_{s\rightarrow 0} sF(s)\qquad s \in \Delta_0 [/math]; が成り立つ。ただし、 $Δ 0$ は $s > 0$ を含む角領域である。

性質一覧表

表中の凡例

[math]u(t)[/math] : ヘビサイド関数

[math](f * g)(t)[/math] :

f

と

g

の畳み込み

[math]f'(t)[/math] :

f (t)

の

1

階微分

[math]f^{(n)}(t)[/math] :

f (t)

の

n

階微分

**片側ラプラス変換の性質（その 1）**
性質	原関数 [math]f(t) = \mathcal{L}^{-1} \left\{ F(s) \right\}[/math] （' $t$ ' 領域 / 時間領域）	像関数 [math]F(s) = \mathcal{L}\left\{ f(t) \right\}[/math] （' $s$ ' 領域 / 周波数領域）	備考
線形性	[math]a f(t) + b g(t) ~[/math]	[math]a F(s) + b G(s) ~[/math]
相似性	[math] f(at) ~[/math]	[math] {1 \over a} F \left ( {s \over a} \right )[/math]	ただし、 $a > 0$
移動	[math] f(t - a) u(t - a) ~[/math]	[math] \mathrm{e}^{-as} F(s) ~[/math]
移動	[math] f(t + \lambda) ~[/math]	[math] \mathrm{e}^{\lambda s}\left\{F(s) - \int_0^{\lambda}\mathrm{e}^{-st}f(t)\,\mathrm dt\right\}~[/math]	移動第 2 則ただし、 $λ > 0$
$1$ 階微分	[math] f'(t) ~[/math]	[math] s F(s) - f(0) ~[/math]	ただし、 $ƒ$ は $1$ 階微分可能とする。
$2$ 階微分	[math]f''(t) \,[/math]	[math]s^{2}F(s)-sf(0)-f'(0) \,[/math]	ただし、 $ƒ$ は $2$ 階微分可能とする。
$n$ 階微分	[math] f^{(n)}(t) ~[/math]	[math] s^n F(s) - s^{n - 1} f(0) - \cdots - f^{(n - 1)}(0) ~[/math]	ただし、 $ƒ$ は $n$ 階微分可能とする。
積分	[math] \int_0^t f(\tau)\, \mathrm d\tau = (u * f)(t) ~[/math]	[math] {1 \over s} F(s) ~[/math]
	[math] \frac{1}{(n-1)!}\int_0^t (t-q)^{n-1}f(q)\mathrm{d}q ~[/math]	[math] \frac{1}{s^n} F(s) ~[/math]	ただし、 $n \geq 1$
	[math] \int_0^t \int_0^{\tau_{n-1}}\cdots \int_0^{\tau_1} f(\tau)\,\mathrm d\tau \mathrm d\tau_1 \cdots \mathrm d\tau_{n-1} ~[/math]	[math] \frac{1}{s^n}F(s) ~[/math]
畳み込み	[math] (f * g)(t) = \int_0^t f(u)\,g(t-u)\,du~[/math]	[math] F(s) \cdot G(s) ~[/math]
周期関数	[math] f(t) = f(t + T) ~[/math]	[math]{1 \over 1 - \mathrm{e}^{-Ts}} \int_0^T \mathrm{e}^{-st} f(t)\,\mathrm dt [/math]	$f (t)$ は周期 $T$ の周期関数。

**片側ラプラス変換の性質（その 2）**
性質	像関数 [math]F(s) = \mathcal{L}\left\{ f(t) \right\}[/math] （' $s$ ' 領域 / 周波数領域）	原関数 [math]f(t) = \mathcal{L}^{-1} \left\{ F(s) \right\}[/math] （' $t$ ' 領域 / 時間領域）	備考
移動	[math] F(s - a) ~[/math]	[math] \mathrm{e}^{at} f(t) ~[/math]
$1$ 階微分	[math] F'(s) ~[/math]	[math] -t f(t) ~[/math]	ただし、 $F$ は $1$ 階微分可能とする。
$2$ 階微分	[math]F^{\prime\prime}(s)\,[/math]	[math]t^{2}f(t) \,[/math]	ただし、 $F$ は $2$ 階微分可能とする。
$n$ 階微分	[math] F^{(n)}(s) ~[/math]	[math] (-t)^{n} f(t) ~[/math]	ただし、 $F$ は $n$ 階微分可能とする。
積分	[math] \int_s^\infty F(\sigma)\, \mathrm d\sigma ~[/math]	[math] \frac{f(t)}{t} ~[/math]
積分	[math] \int_s^{\infty}\int_{\sigma_{n-1}}^{\infty}\cdots\int_{\sigma_1}^{\infty}F(\sigma)\,\mathrm d\sigma \mathrm d\sigma_1 \cdots \mathrm d\sigma_{n-1} ~ [/math]	[math] \frac{f(t)}{t^n} ~[/math]
畳み込み	[math] (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 [/math]	[math] f(t) \,g(t) [/math]

変換表

原関数
[math]f(t) = \mathcal{L}^{-1} \left\{ F(s) \right\}[/math]
'

t

' 領域 / 時間領域

像関数
[math]F(s) = \mathcal{L}\left\{ f(t) \right\}[/math]
'

s

' 領域 / 周波数領域

収束域

単位インパルス

[math] \delta(t) ~[/math]

[math] 1 [/math]

[math] \mathrm{all}~ s \,[/math]

単位ステップ関数

[math] u(t) ~[/math]

[math] { 1 \over s } [/math]

[math] \operatorname{Re} \{ s \} \gt 0 \, [/math]

ランプ関数

[math] t \cdot u(t)~[/math]

[math]\frac{1}{s^2}[/math]

[math] \operatorname{Re} \{ s \} \gt 0 \, [/math]

n

乗
（

n

は整数）

[math]{ t^n \over n! } \cdot u(t) [/math]

[math] { 1 \over s^{n+1} } [/math]

[math] \operatorname{Re} \{ s \} \gt 0 \, [/math]
[math] (n \gt -1) \, [/math]

q

乗
（

q

は複素数）

[math]{ t^q \over \Gamma(q+1) } \cdot u(t) [/math]

[math] { 1 \over s^{q+1} } [/math]

[math] \operatorname{Re} \{ s \} \gt 0 \, [/math]
[math] (\operatorname{Re}\{q\} \gt -1) \, [/math]

n

乗根

[math] \sqrt[n]{t} \cdot u(t) = t^{1/n} \cdot u(t) [/math]

[math] \frac{1}{s^{1+1/n}} \cdot \Gamma\left(1+\frac{1}{n}\right)[/math]

[math] \operatorname{Re} \{ s \} \gt 0 \, [/math]

指数減衰

[math] \mathrm{e}^{-\alpha t} \cdot u(t) ~[/math]

[math] { 1 \over s+\alpha } [/math]

[math] \operatorname{Re} \{ s \} \gt - \alpha ~[/math]

n

乗の指数減衰

[math]\frac{t^{n}}{n!}\mathrm{e}^{-\alpha t} \cdot u(t) [/math]

[math]\frac{1}{(s+\alpha)^{n+1}}[/math]

[math] \operatorname{Re} \{ s \} \gt - \alpha \, [/math]

理想遅延

[math] \delta(t-\tau) ~[/math]

[math] \mathrm{e}^{-\tau s} ~[/math]

遅延付き単位ステップ関数

[math] u(t-\tau) ~[/math]

[math] { \frac{1}{s} \cdot \mathrm{e}^{-\tau s} } [/math]

[math] \operatorname{Re} \{ s \} \gt 0 \, [/math]

遅延付き

n

乗の指数減衰

[math]\frac{(t-\tau)^n}{n!} \mathrm{e}^{-\alpha (t-\tau)} \cdot u(t-\tau) [/math]

[math] \frac{1}{(s+\alpha)^{n+1}} \cdot \mathrm{e}^{-\tau s} [/math]

[math] \operatorname{Re} \{ s \} \gt - \alpha \, [/math]

指数関数的接近

[math]( 1-\mathrm{e}^{-\alpha t}) \cdot u(t) ~[/math]

[math]\frac{\alpha}{s(s+\alpha)} [/math]

[math] \operatorname{Re} \{ s \} \gt 0~[/math]

正弦関数

[math] \sin(\omega t) \cdot u(t) ~[/math]

[math] { \omega \over s^2 + \omega^2 } [/math]

[math] \operatorname{Re} \{ s \} \gt 0 ~[/math]

余弦関数

[math] \cos(\omega t) \cdot u(t) ~[/math]

[math] { s \over s^2 + \omega^2 } [/math]

[math] \operatorname{Re} \{ s \} \gt 0 ~[/math]

双曲線正弦関数
（ハイパボリックサイン）

[math] \sinh(\alpha t) \cdot u(t) ~[/math]

[math] { \alpha \over s^2 - \alpha^2 } [/math]

[math] \operatorname{Re} \{ s \} \gt | \alpha | ~[/math]

双曲線余弦関数
（ハイパボリックコサイン）

[math] \cosh(\alpha t) \cdot u(t) ~[/math]

[math] { s \over s^2 - \alpha^2 } [/math]

[math] \operatorname{Re} \{ s \} \gt | \alpha | ~[/math]

正弦波の指数減衰

[math]\mathrm{e}^{-\alpha t} \sin(\omega t) \cdot u(t) ~[/math]

[math] { \omega \over (s+\alpha )^2 + \omega^2 } [/math]

[math] \operatorname{Re} \{ s \} \gt -\alpha ~[/math]

余弦波の指数減衰

[math]\mathrm{e}^{-\alpha t} \cos(\omega t) \cdot u(t) ~[/math]

[math] { s+\alpha \over (s+\alpha )^2 + \omega^2 } [/math]

[math] \operatorname{Re} \{ s \} \gt -\alpha ~[/math]

自然対数

[math] \ln \left ( { t \over t_0 } \right ) \cdot u(t) [/math]

[math] - {1 \over s} \left[\ln(t_0 s)+\gamma \right] [/math]

[math] \operatorname{Re} \{ s \} \gt 0 \, [/math]

第 1 種ベッセル関数

[math] J_n( \omega t) \cdot u(t)[/math]

[math]\frac{ \omega^n \left(s+\sqrt{s^2+ \omega^2}\right)^{-n}}{\sqrt{s^2 + \omega^2}}[/math]

[math] \operatorname{Re} \{ s \} \gt 0 \, [/math]
[math] (n \gt -1) \, [/math]

第 1 種変形ベッセル関数

[math]I_n(\omega t) \cdot u(t)[/math]

[math] \frac{ \omega^n \left(s+\sqrt{s^2-\omega^2}\right)^{-n}}{\sqrt{s^2-\omega^2}} [/math]

[math] \operatorname{Re} \{ s \} \gt | \omega | \, [/math]

第 2 種ベッセル関数
（次数が

0

の場合）

[math] Y_0(\alpha t) \cdot u(t)[/math]

[math]-{2 \sinh^{-1}(s/\alpha) \over \pi \sqrt{s^2+\alpha^2}}[/math]

[math]\operatorname{Re} \{ s \} \gt 0 \, [/math]

第 2 種変形ベッセル関数
（次数が

0

の場合）

[math] K_0(\alpha t) \cdot u(t)[/math]

誤差関数

[math] \mathrm{erf}(t) \cdot u(t) [/math]

[math] {\mathrm{e}^{s^2/4} \left(1 - \operatorname{erf} \left(s/2\right)\right) \over s}[/math]

[math] \operatorname{Re} \{ s \} \gt 0 \, [/math]

凡例

[math] u(t) \, [/math]：ヘビサイド関数。
[math] \delta(t) \, [/math]：ディラックのデルタ関数。
[math] \Gamma (z) \, [/math]：ガンマ関数。 $n$ が自然数の場合、 $Γ(n + 1) = n!$ 。
[math] \gamma \, [/math]：オイラー・マスケローニ定数.

$t$ は時間に対応する実数。
$s$ は複素数で複素角周波数と呼ばれる。 $Re{s}$ はその実部。
$α, β, τ, ω$ は実数。
$n$ は整数。

外部リンク

Weisstein, Eric W. “Laplace Transform”. MathWorld（英語）. Template:Cite webの呼び出しエラー：引数 accessdate は必須です。

ラプラス変換

Contents

定義

両側ラプラス変換

母関数との関係

性質

線型性

相似性

微分式

積分式

畳み込み

初期値の定理・最終値の定理

性質一覧表

変換表

関連文献

関連項目

外部リンク

案内メニュー

表示

個人用ツール

navigation

search

tool