応用数学 - ラプラス逆変換

概要

ラプラス逆変換について理解する。
部分分数展開について理解する。

ラプラス逆変換の定義

定義:
$F(s)$ に対し、 ${\mathcal {L}}[f(t)]=F(s)$となる $f(t)$ が存在する時、
$f(t)={\mathcal {L}}^{-1}[F(s)]$
と書き、これを $F(s)$ のラプラス逆変換(inverse Laplace transform)という。

説明:
ラプラス変換は全単射ではないので、通常の関数の逆関数と同じようにラプラス逆変換を定義できるわけではない。

ラプラス逆変換の定義の考え方

$F(s)$ に対し、${\mathcal {L}}[f(t)]=F(s)$ となる $f(t)$ が存在する場合、それらは無数に存在して1つに決まらない。
そのような、${\mathcal {L}}[f(t)]=F(s)$ となる $f(t)$ を全部まとめて1グループにして、その中から代表を1つ選んで $f(t)={\mathcal {L}}^{-1}[F(s)]$ とする。
もし、グループの中に連続な関数が存在すれば、それを代表として選び、$f(t)={\mathcal {L}}^{-1}[F(s)]$ と書くことにする。

$F(s)$ に対し、${\mathcal {L}}[f(t)]=F(s)$ となる $f(t)$ が無数に存在することは、以下の定理によって示される。
すなわち、高々有限個の不連続な点で異なる原関数たちは、同一の像関数(ラプラス変換結果)をもつ。

定理:
$f(t)$ と $g(t)$は、$t\geqq 0$ において区分的に連続な関数であるとする。
この時、${\mathcal {L}}[f(t)]={\mathcal {L}}[g(t)]$ ならば、$f(t)$ と $g(t)$ は不連続な点を除いて一致する。

このような不連続な点を除けば一致するような原関数の中に連続な関数があれば、それを代表させてラプラス逆変換の値として採用する。

$t\geqq 0$ において、下図の関数 $f(t)$、$g(t)$、$h(t)$ のラプラス変換の値は一致する。
${\mathcal {L}}[f(t)]={\mathcal {L}}[g(t)]={\mathcal {L}}[h(t)]=F(s)$

この中で $f(t)=1\quad (t\geqq 0)$ をラプラス逆変換の値として採用する。
${\mathcal {L}}^{-1}[F(s)]=f(t)$

ブロムウィチ積分

ラプラス逆変換を、ラプラス変換の定義式のような積分式として与えることも可能である。
そのような定義式は、ブロムウィチ積分(Bromwich integral)あるいは反転公式と呼ばれる。

ブロムウィチ積分は、フーリエの積分定理を使用することにより導くことができる。

ラプラス逆変換を行う場合、通常は変換公式を使用する。
そのため、ブロムウィチ積分の定義式を使用する必要はあまり生じない。

定理:
以下のブロムウィチ積分により、像関数 $F(s)$ から原関数 $f(t)$ を求めることができる。
$f(t)={\frac {1}{2\pi i}}\int _{p-i\infty }^{p+i\infty }{e^{st}F(s)ds}$

ここで、$s=p+i\alpha \quad (p,\,\alpha :{\mbox{ 実 数 }},\,\,i:{\mbox{ 虚 数 単 位 }})$ とする。

説明:
ここまでの話では変数sを実数として扱ってきた。
しかし、$F(s)$ から $f(t)$ を求めるブロムウィチ積分ではsは複素数である。

ブロムウィチ積分は、複素数平面上の虚軸に平行な直線に沿った区間 $p-i\infty$ から $p+i\infty$ の無限積分となる。

ブロムウィチ積分:
$f(t)={\frac {1}{2\pi i}}\int _{p-i\infty }^{p+i\infty }{e^{st}F(s)\,ds}$

$\qquad s=p+i\alpha \quad (p,\,\alpha :{\mbox{ 実 数 }},\,\,i:{\mbox{ 虚 数 単 位 }})$

下図に、ブロムウィチ積分の積分経路を示す。

ブロムウィチ積分は、複素数平面から時間領域への写像である。

複素数平面を用いたブロムウィチ積分は、定義通りの積分経路では計算が難しくなる。
そのため、閉曲線となるように積分経路を変更して留数(residue)を計算する。

ブロムウィチ積分の値は、複素数平面上の全ての特異点の留数の総和となる。

ラプラス逆変換公式の活用

ラプラス逆変換は、ラプラス逆変換公式を使用して計算する。

ラプラス逆変換公式は、先のラプラス逆変換の定義やブロムウィチ積分を使用して基本的な関数に対してあらかじめ求められた値を公式として表したものである。

ラプラス変換とラプラス逆変換を使用した微分方程式の解法では、そのようなラプラス逆変換公式を活用する。
ラプラス逆変換公式とラプラス変換公式は本質的に同じものである。

ラプラス逆変換公式

下表のような原関数 $f(t)$ と像関数 $F(s)$ の対応関係を変換表として使用すれば、ラプラス逆変換を機械的に実行できる。

$f(t)$	$F(s)$	条件
$1$	${\frac {1}{s}}$	$s>0$
$t$	${\frac {1}{s^{2}}}$	$s>0$
${\frac {t^{n-1}}{(n-1)!}}$	${\frac {1}{s^{n}}}$	$n$ は正整数
${\frac {t^{p-1}}{\Gamma (p)}}$	${\frac {1}{s^{p}}}$	$p>0$
$e^{at}$	${\frac {1}{s-a}}$	$s>\alpha$
$\sin {at}$	${\frac {a}{s^{2}+a^{2}}}$	$s>0,\quad a:{\mbox{ 実 定 数　}}$
$\cos {at}$	${\frac {s}{s^{2}+a^{2}}}$	$s>0,\quad a:{\mbox{ 実 定 数　}}$
$\sinh {at}$	${\frac {a}{s^{2}-a^{2}}}$	$s>0,\quad a:{\mbox{ 実 定 数　}}$
$\cosh {at}$	${\frac {s}{s^{2}-a^{2}}}$	$s>0,\quad a:{\mbox{ 実 定 数　}}$
単位階段関数 $U(t-a)$	${\frac {e^{-as}}{s}}$	$s>0,\quad a>0$
デルタ関数 $\delta (t)$	$1$	$-\infty <s<\infty$

法則

$f(t)$	$F(s)$	条件
線形法則 $af(t)+bg(t)$	$aF(s)+bF(s)$
相似法則 ${\frac {1}{a}}f\left({\frac {t}{a}}\right)$	$F(as)$
第1移動法則 $e^{at}f(t)$	$F(s-a)$
第2移動法則 $U(t-a)\,f(t-a)$	$e^{-as}\,F(s)$	$a>0$
微分法則 ${\frac {df(t)}{dt}}$ ${\frac {d^{2}f(t)}{dt^{2}}}$	$sF(s)-f(0)$ $s^{2}F(s)-sf(0)-{\frac {df(0)}{ds}}$
積分法則 $\int _{0}^{t}{f(u)\,du}$	${\frac {F(s)}{s}}$
$t^{n}$ 積法則 $t^{n}f(t)$	$(-1)^{n}\,{\frac {dF^{n}(s)}{ds}}$	$n:{\mbox{ 正 整 数 }}$
合成積 $(f*g)(t)$	$F(s)\,G(s)$
$-tf(t)$	構文解析に失敗 (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle {\frac {dF(s)}{ds}}}
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle (-1)^{n} \, t^{n} \, f(t)}	構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{dF^{n} (s)}{ds^{n}}}	構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle n = 1, 2, 3, \cdots}
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle erf(\sqrt{t})}	構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{1}{s \sqrt{s + 1}}}

ラプラス逆変換公式は、ラプラス変換公式を基にして作成されている。
ラプラス逆変換公式とラプラス変換公式は本質的に同じものであり、ラプラス逆変換公式で示されているように、ラプラス逆変換もラプラス変換と同様の性質(線形法則、移動法則、微分・積分法則、等)をもつ。

ラプラス逆変換公式における逆変換の値は、なるべく 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle t > 0} において連続な関数を選んである。 (構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle t < 0} ではどんな値でもよい)

なお、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{1}{s}} のラプラス逆変換は 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle f(t) = 1} でも 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle f(t) = U(t)} でもよい。

部分分数展開

複雑な関数のラプラス逆変換を行う際、部分分数展開を使用するとうまくいく場合がある。

以下の公式における係数A、B、C、Dは定数とする。
A、B、C、Dを計算して求める。

• 公式 1

  構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{1}{(s - a)(s - b)} = \frac{A}{s - a} + \frac{B}{s - b}}

  
  

• 公式 2

  構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{1}{(s - a)^{2} \, (s - b)} = \frac{A}{s - a} + \frac{B}{(s - a)^{2}} + \frac{C}{s - b}}

  
  

• 公式 3

  構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{1}{(s - a)(s^{2} + b)} = \frac{A}{s - a} + \frac{Bs + C}{s^{2} + b} \quad (b > 0)}

  
  

• 公式 4
  分母が全てsの1次式の積の形に因数分解される場合
  構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle P(s)} は 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle s} の2次以下の式とする。

  構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{P(s)}{(s - a)(s - b)(s - c)} = \frac{A}{s - a} + \frac{B}{s - b} + \frac{C}{s - c}}

  
  

• 公式 5
  分母 = 0が多重根をもつ場合
  構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle P(s)} は 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle s} の3次以下の式とする。

  構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{P(s)}{(s - a)^{3} \, (s - b)} = \frac{A}{s - a} + \frac{B}{(s - a)^{2}} + \frac{C}{(s - a)^{3}} + \frac{D}{s - b}}

  
  

• 公式 6
  分母 = 0が多重根をもち、さらにsの2次式の項も存在する場合
  構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle P(s)} は 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle s} の4次以下の式とする。

  構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle P(s)/(s - a)3(s2 + bs + c) = A/(s - a) + B/(s - a)2 + C/(s - a)3 + (Ds + E)/(s2 + bs + c)}