全記事一覧 > 数学 > フーリエ変換／ラプラス変換

【例題で学ぶ】ラプラス逆変換（線形微分方程式の初期値問題）

　例題を使って微分方程式の初期値問題をラプラス変換で解く。 初期値が与えられている微分方程式はラプラス変換によってをかんたんに解くことができる。 ここでは下の例題のような単純な斉次線形微分方程式を解いていく（非斉次微分方程式のラプラス変換による解法は こちら ）。

例題

以下のx(t)に関する微分方程式をラプラス変換を用いて解け。ただし初期値はx(0)=x_0,\,\dot{x}(0)=v_0とする。

    \begin{eqnarray*} &&(1)\quad\ddot{x}(t)+3x(t)=0\\ &&(2)\quad\ddot{x}(t)+5\dot{x}(t)+4x(t)=0\\ &&(3)\quad\ddot{x}(t)+4\dot{x}(t)+4x(t)=0\\ &&(4)\quad\ddot{x}(t)+4\dot{x}(t)+13x(t)=0 \end{eqnarray*}



1. ラプラス変換・逆変換の基本事項

　基本事項をざっと復習する。

微分のラプラス変換

　微分方程式をラプラス変換で解く上で大事な x'(t),x''(t) のラプラス変換を簡単にまとめる。

１回微分のラプラス変換：部分積分する

    \begin{eqnarray*} {\mathcal L}[x'(t)](s)&=&\int_0^{\infty}x'(t)e^{-st}\,dt\\\\ &=&\left[x(t)e^{-st}\right]_0^{\infty}-(-s)\int_0^{\infty}x(t)e^{-st}\,dt\\\\ &=&-x(0)+sX(s)\quad(X(s)\equiv{\mathcal L}[x(t)](s))\quad\blacksquare \end{eqnarray*}


２回微分のラプラス変換：上の結果を使う（x\rightarrowx'と置き換える）

    \begin{eqnarray*} {\mathcal L}[x''(t)](s) &=&-x'(0)+s\textcolor{red}{\tilde{X}(s)} \quad(\textcolor{red}{\tilde{X}(s)}\equiv{\mathcal L}[x'(t)](s))\\\\ &=&-x'(0)+s\left[\textcolor{red}{-x(0)+sX(s)}\right] \quad(X(s)\equiv{\mathcalL}[x(t)](s))\\\\ &=&-x'(0)+sx(0)+s^2X(s)\quad\blacksquare \end{eqnarray*}


s移動を用いたラプラス逆変換

　s 移動について簡単にまとめる。以下では「⇒」をラプラス逆変換とする。


    \begin{eqnarray*} X(s)\quad\Rightarrow\quad x(t) \end{eqnarray*}

に対して

    \begin{eqnarray*} X(\textcolor{red}{s-a})\quad\Rightarrow\quad\textcolor{red}{e^{at}}x(t) \end{eqnarray*}


例えば、

    \begin{eqnarray*} \frac{1}{(s-a)^2}\quad\Rightarrow\quad te^{at} \end{eqnarray*}

である。


2. 解答

　以下の解答ではx(t)のラプラス変換を X(s)={\mathcalL}[x(t)](s) とする。 また、与えられた初期値は x(0)=x_0,\,\dot{x}(0)=v_0 である。


例題(1)の解答

　与えられた微分方程式の両辺をラプラス変換すると

    \begin{eqnarray*} &&s^2X(s)-sx_0-v_0+3X(s)=0\\\\ \Leftrightarrow\quad&&(s^2+3)X(s)=sx_0+v_0\\\\ \Leftrightarrow\quad&&X(s)=\frac{s}{s^2+\sqrt{3}^2}\,x_0 +\frac{\textcolor{blue}{\sqrt{3}}}{s^2+\sqrt{3}^2} \cdot\frac{v_0}{\textcolor{blue}{\sqrt{3}}} \end{eqnarray*}

\sin,\cos のラプラス逆変換が使えるように 3=\sqrt{3}^2 とした。 最後の行を「s 移動」を用いてラプラス逆変換すると

    \begin{eqnarray*} x(t)=x_0\cos\sqrt{3}t+\frac{v_0}{\sqrt{3}}\sin\sqrt{3}t \quad\blacksquare \end{eqnarray*}


※ここで 3=(-\sqrt{3})^2 と考えても同じ結果になる。つまり、

    \begin{eqnarray*} x(t)&=&x_0\cos(-\sqrt{3})t+\frac{v_0}{-\sqrt{3}}\sin(-\sqrt{3})t\\\\ &=&x_0\cos\sqrt{3}t+\frac{v_0}{\sqrt{3}}\sin\sqrt{3}t\quad\blacksquare \end{eqnarray*}


例題(2)の解答

　与えられた微分方程式の両辺をラプラス変換すると

    \begin{eqnarray*} &&s^2X(s)-sx_0-v_0+5\left\{sX(s)-x_0\right\}+4X(s)=0\\\\ \Leftrightarrow\quad&& (s+1)(s+4)X(s)=sx_0+v_0+5x_0\\\\ \Leftrightarrow\quad&& X(s)=\textcolor{blue}{\frac{s}{(s+1)(s+4)}}x_0 +\textcolor{red}{\frac{1}{(s+1)(s+4)}}\cdot\left\{v_0+5x_0\right\} \end{eqnarray*}


ここで、部分分数分解

    \begin{eqnarray*} &&\textcolor{blue}{\frac{s}{(s+1)(s+4)}} =\left(\frac{4}{s+4}-\frac{1}{s+1}\right)\cdot\frac{1}{3}\\\\ &&\textcolor{red}{\frac{1}{(s+1)(s+4)}} =\left(\frac{1}{s+4}-\frac{1}{s+1}\right)\cdot\left(-\frac{1}{3}\right) \end{eqnarray*}

を用いて、

    \begin{eqnarray*} X(s)&=&\left(\frac{4}{s+4}-\frac{1}{s+1}\right)\cdot\frac{x_0}{3}\\\\ &+&\left(\frac{1}{s+4}-\frac{1}{s+1}\right)\cdot\left(-\frac{v_0+5x_0}{3}\right)\cdot\\\\ &=&\frac{1}{s+4}\cdot\left(-\frac{v_0+x_0}{3}\right) +\frac{1}{s+1}\cdot\left(\frac{v_0+4x_0}{3}\right) \end{eqnarray*}

となる。\frac{1}{s+4} と \frac{1}{s+1} に対して「s 移動」を用いたラプラス逆変換をおこなうと、

    \begin{eqnarray*} x(t)=-\frac{v_0+x_0}{3}\,e^{-4t}+\frac{v_0+4x_0}{3}\,e^{-t}\quad\blacksquare \end{eqnarray*}


例題(3)の解答

　 与えられた微分方程式の両辺をラプラス変換すると

    \begin{eqnarray*} &&s^2X(s)-sx_0-v_0+4\left\{sX(s)-x_0\right\}+4X(s)=0\\\\ \Leftrightarrow\quad&& (s+2)^2X(s)=sx_0+v_0+4x_0\\\\ \Leftrightarrow\quad&& X(s)=\textcolor{blue}{\frac{1}{(s+2)^2}}x_0 +\frac{1}{(s+2)^2}\cdot\left\{v_0+4x_0\right\}\\\\ \Leftrightarrow\quad&& \textcolor{blue}{\left[\frac{1}{s+2}-\frac{2}{(s+2)^2}\right]}x_0 +\frac{1}{(s+2)^2}\left\{v_0+4x_0\right\}\\\\ \Leftrightarrow\quad&& \frac{x_0}{s+2}+\frac{v_0+2x_0}{(s+2)^2} \end{eqnarray*}


「s 移動」を用いてラプラス逆変換すると、

    \begin{eqnarray*} x(t) &=&x_0e^{-2t}+\left(v_0+2x_0\right)\,te^{-2t}\\\\ &=&\Bigl[x_0+(v_0+2x_0)t\Bigr]e^{-2t}\quad\blacksquare \end{eqnarray*}


例題(4)の解答

　与えられた微分方程式の両辺をラプラス変換すると

    \begin{eqnarray*} &&s^2X(s)-sx_0-v_0+4\left\{sX(s)-x_0\right\}+13X(s)=0\\\\ \Leftrightarrow\quad&& (s^2+4s+13)X(s)=sx_0+v_0+4x_0\\\\ \Leftrightarrow\quad&& \Bigl[(s+2)^2+3^2\Bigr]X(s)=sx_0+v_0+4x_0\\\\ \Leftrightarrow\quad&& X(s)=\frac{s}{(s+2)^2+3^2}x_0+\frac{v_0+4x_0}{(s+2)^2+3^2}\\\\ \Leftrightarrow\quad&& X(s)=\frac{\textcolor{red}{(s+2)}-2}{(s+2)^2+3^2}x_0 +\frac{v_0+4x_0}{(s+2)^2+3^2}\\\\ \Leftrightarrow\quad&& X(s)=\frac{\textcolor{red}{(s+2)}}{(s+2)^2+3^2} x_0+\frac{v_0+2x_0}{(s+2)^2+3^2}\\\\ \Leftrightarrow\quad&& X(s)=\frac{\textcolor{red}{(s+2)}}{(s+2)^2+3^2}x_0 +\frac{\textcolor{blue}{3}}{(s+2)^2+3^2}\cdot \frac{v_0+2x_0}{\textcolor{blue}{3}} \end{eqnarray*}


「s移動」と\sin3t,\cos3tのラプラス変換（逆変換）より、

    \begin{eqnarray*} x(t)=x_0 e^{-2t}\cos3t+\frac{v_0+2x_0}{3}e^{-2t}\sin3t\quad\blacksquare \end{eqnarray*}



3. これらの解の物理的意味

　得られた解について簡単に触れておく。 (1)から(4)の微分方程式は振動の運動方程式と関係する。

(1)は１次元の単振動（自由振動）を表す運動方程式である。 したがって得られた解は

    \begin{eqnarray*} x(t)=A\cos(\sqrt{3}t+\delta) \end{eqnarray*}

のように初期位相 \delta と振幅 A で表すことができる。


　 (2)(3)(4)の微分方程式に含まれる \dot{x}(t) の項は 係数が0より大きいため摩擦項に対応する（下図）。

それぞれ、\dot{x}(t) の係数と x(t) の係数の関係より解の意味が異なる。 (2)の解はもはや振動的な解ではない。これは過減衰(over damping)に対応する運動である。 (3)の解は臨界減衰(critical damping)に対応する。 (4)の解は三角関数の振動的性質の残った減衰振動（under damping）である。 これは解を

    \begin{eqnarray*} x(t)=Ae^{-2t}\cos(3t+\delta) \end{eqnarray*}

と書いてみればわかる。



4. まとめ

一般に、ここで挙げた４タイプは以下のように分類できる。

    \begin{eqnarray*} &&(1)\quad\ddot{x}(t)+\omega^2 x(t)=0\\ &&(2)\quad\ddot{x}(t)-(\alpha+\beta)\dot{x}(t)+\alpha\beta x=0 \quad(\beta\neq \alpha)\\&&(3)\quad\ddot{x}(t)-2\alpha\dot{x}(t)+\alpha^2 x=0\\ &&(4)\quad\ddot{x}(t)-2\alpha\dot{x}(t)+(\alpha^2+\beta^2)x=0 \quad(\beta\neq 0)\end{eqnarray*}

一般的には、\alpha と \beta の組み合わせによって、上のように過減衰、減衰振動、臨界振動の解にならないことがある。

また、本例題の微分方程式はラプラス変換を使わなくても解くことができる。x(t)=e^{\lambda t} とおいて特性方程式を作って解くやり方などがある。この方法は力学の教科書にも載っている よくある解き方になる。



2件のコメント

  1. mark

    ボケ防止目的で頭の運動のため利用させていただいておりありがとうございます。
    例題４の回答のCos3tの係数にX0を掛ける必要がありませんでしょうか？
    私の計算ミスであれば申し訳ありません。

    返信
    • batapara

      係数x0は必要です。途中式で消えていたようです。
      ご指摘ありがとうございました。

      返信

コメントを残す

メールアドレスが公開されることはありません。