ラプラス変換の定義と具体例・性質

定義

関数 f(x)f(x) に対して L[f](s)=0esxf(x)dx \mathscr{L} [f] (s)=\int_0^{\infty} e^{-sx}f(x)\, dx という「ss についての関数を返す変換」をラプラス変換という。

この記事では,ラプラス変換の基本的な性質と公式について紹介します。

ラプラス変換表

主要な関数のラプラス変換を紹介します。

1. 直接計算できるもの

この記事の後半で証明します。

f(t)L[f](s)f(t)L[f](s)aaseat1s+aatas2sinat1s2+a2tnn!sn+1cosatss2+a2tαΓ(α+1)sα+1δ(s)1 \begin{array}{c|c||c|c} f(t) & \mathscr{L} [f] (s) & f(t) & \mathscr{L} [f] (s)\\ \hline\hline a & \dfrac{a}{s} & e^{-at} & \dfrac{1}{s+a}\\ \hline at & \dfrac{a}{s^2} & \sin at & \dfrac{1}{s^2+a^2} \\ \hline t^n & \dfrac{n!}{s^{n+1}} & \cos at & \dfrac{s}{s^2+a^2} \\ \hline t^{\alpha} & \dfrac{\Gamma (\alpha + 1)}{s^{\alpha + 1}} & \delta(s) & 1 \end{array}

aa は複素数,nn は正整数,α\alpha1-1 より大きい実数の定数とする。

2. ラプラス変換の性質を組み合わせるもの

前の表の関数と,この後説明する基本的な性質の組み合わせで計算できるものです。

f(t)L[f](s)f(t)L[f](s)teat1(s+a)2cos(at+θ)scosθasinθs2+a2tneatn!(s+a)n+1sin(at+θ)acosθ+ssinθs2+a2sinhatas2a2eatsinbta(s+a)2+b2coshatss2a2eatcosats+a(s+a)2+b2 \begin{array}{c|c||c|c} f(t) & \mathscr{L} [f] (s) & f(t) & \mathscr{L} [f] (s)\\ \hline\hline te^{-at} & \dfrac{1}{(s+a)^2} & \cos (at + \theta) & \dfrac{s \cos \theta - a \sin \theta}{s^2 + a^2} \\ \hline t^n e^{-at} & \dfrac{n!}{(s+a)^{n+1}} & \sin (at + \theta) & \dfrac{a \cos \theta + s \sin \theta}{s^2 + a^2} \\ \hline \sinh a t & \dfrac{a}{s^2-a^2} &e^{-at} \sin b t & \dfrac{a}{(s+a)^2+b^2}\\ \hline \cosh at& \dfrac{s}{s^2-a^2} & e^{-at} \cos at & \dfrac{s+a}{(s+a)^2+b^2} \\ \end{array}

基本的な性質

線型性

定理

L[af+bg]=aL[f]+bL[g] \mathscr{L}[af+bg] = a \mathscr{L}[f] + b \mathscr{L}[g]

証明

積分の線型性より明らか。

スケール変換

定理

L[f(ax)](s)=1aL[f(x)](sa) \mathscr{L}[f(ax)] (s)= \frac{1}{a} \mathscr{L} [f(x)] \left( \frac{s}{a} \right)

証明

定義より L[f(ax)](s)=0esxf(ax)dx \mathscr{L}[f(ax)] (s) = \int_0^{\infty} e^{-sx}f(ax) dx である。y=axy = ax と置換すると L[f(ax)](s)=0esaxf(y)dya=1aL[f(x)](sa)\begin{aligned} \mathscr{L}[f(ax)] (s) &= \int_0^{\infty} e^{-\frac{s}{a} x}f(y) \dfrac{dy}{a}\\ &= \frac{1}{a} \mathscr{L} [f(x)] \left( \frac{s}{a} \right) \end{aligned} である。

微分・積分との関係

導関数・原始関数のラプラス変換は,元の関数のラプラス変換を用いて記述できます。

定理

L[f(x)](s)=sL[f](s)f(0)\begin{aligned} \mathscr{L} [f'(x)] (s) &= s \mathscr{L}[f] (s) - f(0)\\ \end{aligned}

証明

L[f(x)](s)=0ddtf(t)estdt=[f(t)est]00(s)f(t)estdt=sL[f](s)f(0)\begin{aligned} &\mathscr{L} [f'(x)] (s)\\ &= \int_{0}^{\infty} \dfrac{d}{dt} f(t) e^{-st} dt\\ &= \Big[ f(t) e^{-st} \Big]_0^{\infty} - \int_0^{\infty} (-s) f(t)e^{-st} dt\\ &= s \mathscr{L}[f] (s) - f(0) \end{aligned}

定理

L[F(x)](s)=1sL[f(x)](s) \mathscr{L} [ F(x) ] (s) = \frac{1}{s} \mathscr{L}[f(x)] (s)

ただし,FFff の原始関数である。

証明

L[0xf(t)dt](s)=0(0tf(u)du)estdt=[(0tf(u)du)(1sest)]0+01sf(t)estdt=1sL[f(x)](s)\begin{aligned} &\mathscr{L} \Big[ \int_0^x f(t) dt \Big] (s)\\ &= \int_{0}^{\infty} \left( \int_0^t f(u) du \right) e^{-st} dt\\ &= \Big[\left( \int_0^t f(u) du \right) \left(-\frac{1}{s} e^{-st} \right) \Big]_0^{\infty} + \int_0^{\infty} \frac{1}{s} f(t) e^{-st} dt\\ &= \frac{1}{s} \mathscr{L}[f(x)] (s) \end{aligned}

「平行移動」

ラプラス変換の平行移動は,指数関数倍によって記述できます。

定理

L[eaxf(x)](s)=L[f](sa) \mathscr{L} [e^{ax} f(x)] (s) = \mathscr{L} [f] (s-a)

証明は難しくないため省きます。

計算例

様々な例を計算してみましょう。

定数・多項式

定理

L[1](s)=1sL[x](s)=1s2L[xa](s)=Γ(a+1)sa+1\begin{aligned} \mathscr{L}[1] (s) &= \frac{1}{s}\\ \mathscr{L}[x] (s) &= \frac{1}{s^2}\\ \mathscr{L}[x^{a}] (s) &= \frac{\Gamma (a+1)}{s^{a+1}} \end{aligned}

証明
  1. 定数関数 L[1]=0estdt=1s\begin{aligned} \mathscr{L}[1] &= \int_{0}^{\infty} e^{-st} dt\\ &=\frac{1}{s} \end{aligned}

  2. 変数 L[x]=0testdt=[1stest1s2est]0=1s2\begin{aligned} \mathscr{L}[x] &= \int_{0}^{\infty} te^{-st} dt\\ &= \Big[ -\dfrac{1}{s} te^{-st} - \dfrac{1}{s^2} e^{-st} \Big]_0^{\infty}\\ &= \frac{1}{s^2}\\ \end{aligned}

  3. 冪関数 L[xa]=0taestdt=01sa(st)aestdt=01sa+1(st)(a+1)1estdt=Γ(a+1)sa+1\begin{aligned} \mathscr{L}[x^{a}] &= \int_{0}^{\infty} t^{a} e^{-st} dt\\ &= \int_{0}^{\infty} \frac{1}{s^{a}} (st)^{a} e^{-st} dt\\ &= \int_{0}^{\infty} \frac{1}{s^{a+1}} (st)^{(a+1)-1} e^{-st} dt\\ &= \frac{\Gamma (a+1)}{s^{a+1}}\\ \end{aligned} → ガンマ関数

三角関数・指数関数

定理

L[eax]=1s+aL[sinx]=1s2+1L[cosx]=ss2+1L[sinhx](s)=1s21L[coshx](s)=ss21\begin{aligned} \mathscr{L}[e^{-ax}] &= \frac{1}{s+a}\\ \mathscr{L}[\sin x] &= \frac{1}{s^2+1} \\ \mathscr{L}[\cos x] &= \frac{s}{s^2+1}\\ \mathscr{L} [\sinh x] (s) &= \dfrac{1}{s^2-1}\\ \mathscr{L} [\cosh x] (s) &= \dfrac{s}{s^2-1} \end{aligned}

証明
  • 指数関数 L[eax](s)=0eatestdt=1s+a\begin{aligned} \mathscr{L}[e^{-ax}] (s) &= \int_{0}^{\infty} e^{-at}e^{-st} dt\\ &= \frac{1}{s+a}\\ \end{aligned}

  • 三角関数 L[sinax](s)=0estsintdt=as2+1L[cosax](s)=0estcostdt=ss2+1\begin{aligned} \mathscr{L}[\sin ax] (s) &= \int_{0}^{\infty} e^{-st} \sin t dt\\ &= \frac{a}{s^2+1} \\ \mathscr{L}[\cos ax] (s) &= \int_{0}^{\infty} e^{-st} \cos t dt\\ &= \frac{s}{s^2+1} \end{aligned} → 三角関数と指数関数の積の積分

  • 双曲関数

定義 sinht=etet2\sinh t = \dfrac{e^{t} - e^{-t}}{2}cosht=et+et2\cosh t = \dfrac{e^t + e^{-t}}{2} と指数関数の結果から従う。

デルタ関数

定理

L[δ(x)](s)=1 \mathscr{L}[\delta (x)] (s) = 1

証明

L[δ(x)](s)=0δ(t)estdt=e0=1\begin{aligned} \mathscr{L}[\delta (x)] (s) &= \int_{0}^{\infty} \delta (t) e^{-st} dt \\ &= e^0 = 1 \end{aligned}

その他の性質

畳み込み

ラプラス変換は畳み込みと相性が良いです。

定理

L[fg](s)=L[f](s)L[g](s) \mathscr{L} [f \ast g] (s) = \mathscr{L} [f] (s) \mathscr{L} [g] (s)

ただし 畳み込み (fg)(x)=f(t)g(xt)dt(f \ast g)(x) = \int_{-\infty}^{\infty} f(t)g(x-t) dt である。

証明

L[fg](s)=L[f(t)g(xt)dt](s)=(f(t)g(xt)dt)esxdx\begin{aligned} &\mathscr{L} [f*g] (s)\\ &=\mathscr{L}\left[ \int_{-\infty}^{\infty} f(t) g(x-t) dt \right] (s)\\ &= \int_{-\infty}^{\infty}\left(\int_{-\infty}^{\infty}f(t)g(x-t)dt\right)e^{-sx}dx \end{aligned}

フビニの定理より積分の順序を交換すると,上式は

(g(xt)esxdx)f(t)dt \int_{-\infty}^{\infty}\left(\displaystyle\int_{-\infty}^{\infty}g(x-t)e^{-sx}dx\right)f(t)dt

となる。この xx についての積分は xt=yx-t=y と置換すると

g(y)esydy×est=estL[g](s) \int_{-\infty}^{\infty}g(y)e^{-sy}dy\times e^{-st}=e^{-st}\mathscr{L} [g] (s)

となる。よって L[fg]=estL[g](s)f(t)dt=L[f](s)L[g](s)\begin{aligned} &\mathscr{L} [f \ast g] \\ &= \int_{-\infty}^{\infty} e^{-st} \mathscr{L} [g] (s) f(t) dt\\ &= \mathscr{L} [f] (s) \mathscr{L} [g] (s) \end{aligned} である。

ラプラス逆変換

定義

ラプラス逆変換 L1\mathscr{L}^{-1}L1[f](s)=limp12πicpic+pif(t)estdt \mathscr{L}^{-1} [f] (s) = \lim_{p \to \infty} \dfrac{1}{2\pi i} \int_{c-pi}^{c+pi} f(t) e^{st} dt と定める。ただし,cc は十分大きい実数とする。(※ 参照)

このとき L1L[f]=f, LL1[F]=F \mathscr{L}^{-1} \mathscr{L} [f] = f, \ \mathscr{L} \mathscr{L}^{-1} [F] = F となる。

cc は,Re(z)>c\mathrm{Re} (z) > cF(z)F(z) が正則となるように取る。

証明や計算方法については次回紹介します。

余談

関数の似たような変換として,フーリエ変換があります。

フーリエ変換

関数 f(x)f(x) のフーリエ変換 f^(ξ)\hat{f} (\xi)f^(ξ)=f(x)eixξdx \hat{f} (\xi) = \int_{-\infty}^{\infty} f(x) e^{-ix \xi} dx と定める。

ラプラス変換と同様に線型性が成り立ちます。また,微分・積分についても良いふるまいをします。

詳しくは フーリエ変換の意味と応用例 をご覧ください。

微分方程式の計算のときなどに用いられます。