logo

周期関数のラプラス変換 📂微分方程式

周期関数のラプラス変換

公式

ffを周期がTT周期関数としよう。するとf(t+T)=f(t)f(t+T)=f(t)であり、f(t)f(t)のラプラス変換は以下のようになる。

L{f(t)}=0estf(t)dt=0Testf(t)dt1est \mathcal{L} \left\{ f(t) \right\} = \int_{0}^\infty e^{-st}f(t)dt = \frac{\displaystyle \int_{0}^T e^{-st}f(t)dt}{1-e^{-st}}

導出

ラプラス変換の定義から、次のように積分を分けよう。

0estf(t)dt=0Testf(t)dt+T2Testf(t)dt+2T3Testf(t)dt+ \int_{0}^\infty e^{-st}f(t)dt = \int_{0}^T e^{-st}f(t)dt + \int_{T}^{2T} e^{-st}f(t)dt + \int_{2T}^{3T}e^{-st}f(t)dt + \cdots

この時、2番目の項の積分範囲を1番目の項と同じにするために、t=t+Tt=t^{\prime}+Tで置換しよう。すると

T2Testf(t)dt=0Tes(t+T)f(t+T)dt \int_{T}^{2T} e^{-st}f(t)dt=\int_{0}^T e^{-s(t^{\prime}+T)}f(t^{\prime}+T)dt^{\prime}

ffは周期がTTの周期関数であるからf(t+T)=f(t)f(t^{\prime}+T)=f(t^{\prime})であり、定数を前に出すと以下のように整理される。

esT0Testf(t)dt=esT0Testf(t)dt e^{-sT}\int_{0}^T e^{-st^{\prime}}f(t^{\prime})dt^{\prime}=e^{-sT}\int_{0}^T e^{-st}f(t)dt

3番目の項もその後もt=t+2T,t=t+3T,t=t^{\prime}+2T, t=t^{\prime}+3T, \cdotsのように置換すると同じ方法で整理できる。そうすると、f(t)f(t)のラプラス変換は以下の形で表せる。

0estf(t)dt=0Testf(t)dt+esT0Testf(t)dt+e2sT0Testf(t)dt+=(1+esT+e2sT+)0Testf(t)dt \begin{align*} & \int_{0}^\infty e^{-st}f(t)dt \\ &= \int_{0}^T e^{-st}f(t)dt + e^{-sT}\int_{0}^T e^{-st}f(t)dt + e^{-2sT}\int_{0}^T e^{-st}f(t)dt + \cdots \\ &= \left( 1+e^{-sT} + e^{-2sT} + \cdots \right) \int_{0}^T e^{-st}f(t)dt \end{align*}

無限等比級数

r<1|r|<1の時、

n=1arn1=a1r \sum_{n=1}^{\infty} a r^{n-1} = { a \over {1-r}}

前にまとめた項を無限等比級数の公式で表すと以下のようになる。

11esT \frac{1}{1-e^{-sT}}

したがって、以下を得る。

0estf(t)dt=0Testf(t)dt1esT \int_{0}^\infty e^{-st}f(t)dt = \dfrac{\displaystyle \int_{0}^T e^{-st}f(t)dt}{1-e^{-sT}}

参考