두 함수 합성곱(convolution)은 하나의 함수 $f$와 또 다른 함수 $g$를 반전 이동한 값을 곱한 다음, 구간에 대해 적분하여 새로운 함수를 구하는 수학 연산자이다.
일반적으로 라플라스 변환은 분배법칙이 성립하지 않는다.
$$ \mathcal{L}(fg) \neq \mathcal{L}(f) \mathcal{L}(g) $$
간단한 예로, 두 함수 $f = e^t$와 $g = 1$를 보면,
$$ \mathcal{L}(fg) = \frac{1}{s-1} \neq \frac{1}{s} \frac{1}{s-1} = \mathcal{L}(f) \mathcal{L}(g) $$
임을 보일 수 있다. 하지만 두 함수 $f$와 $g$의 합성곱 $f \ast g$에 대해서 다음이 성립한다.
증명. 첫째로, 다음을 확인하자.
$$ \begin{align*} G(s) & = \int_{0}^{\infty} e^{-s\sigma} g(\sigma) \, d\sigma \\[5px] & = \int_{\tau}^{\infty} e^{-s(t-\tau)} g(t-\tau) \, d\sigma \qquad (\sigma = t-\tau) \\[5px] & = e^{s\tau}\int_{\tau}^{\infty} e^{-st} g(t-\tau) \, d\sigma \end{align*} $$
따라서
$$ \begin{align*} F(s)G(s) & = \left( \int_{0}^{\infty} e^{-s\tau} f(\tau) \, d\tau \right) G(s) \\[5px] & = \int_{0}^{\infty} \int_{\tau}^{\infty} f(\tau) e^{-st} g(t-\tau) \, dt d\tau \\[5px] & = \int_{0}^{\infty} \int_{0}^{t} f(\tau) e^{-st} g(t-\tau) \, d\tau dt \\[5px] & = \int_{0}^{\infty} e^{-st} (f \ast g)(t) \, dt \\[5px] & = \mathcal{L}(f \ast g) \end{align*} $$
예제 19. 함수 $H(s) = (s^2 + w^2)^{-2}$가 있다고 하자. $\mathcal{L}^{-1}(H)$을 구하기 위하여, $F(s) = (s^2 + w^2)^{-1}$를 정의하자. 그러면, $H(s) = F(s)F(s)$임을 알수있고, 따라서
$$ \begin{align*} \mathcal{L}^{-1}(H) & = f(t) \ast f(t) \\[5px] & = \frac{1}{w} \sin wt \ast \frac{1}{w} \sin wt \\[5px] & = \frac{1}{w^2} \int_{0}^{t} \sin w\tau \sin w(t-\tau) \, d\tau \\[5px] & = \frac{1}{2w^2} \int_{0}^{t} -\cos wt + \cos w(2\tau-t) \, d\tau \\[5px] & = \frac{1}{2w^2} \left( \left. -\tau \cos wt + \frac{\sin w(2\tau-t)}{2w} \right|_{0}^{t} \right) \\[5px] & = \frac{1}{2w^2} \left( -t \cos wt + \frac{\sin wt}{w} \right) \end{align*} $$
합성곱의 기본적인 성질들은 다음과 같다.
이제 합성곱을 이용한 미분방정식(ordinary differential equation)에의 응용에 대해 알아보자.
다음의 미분방정식
$$ \begin{cases} y'' + ay' + by = r(t), \\[5px] y(0) = y'(0) = 0. \end{cases} $$
에 라플라스 변환을 적용하면,
$$ Y = \frac{1}{s^2 + as + b} R =: QR $$
따라서, 합성곱의 성질을 이용하면, 미분방정식의 해 $y$는 다음과 같이 구할 수 있다.
$$ y = q \ast r = \int_{0}^{t} q(t-\tau) r(t) \, d\tau $$
다시 말해 위의 미분방정식의 해를 구하는데 있어서 $R$을 구할 필요가 없다는 사실을 알 수 있다.
예제 20. 다음의 미분방정식
$$ y'' + 3y' + 2y = r(t) = \begin{cases} 1, & \qquad 1 < t < 2 \\[5px] 0, & \qquad \text{otherwise} \end{cases} $$
에 대해 생각해 보자. 이때 초기값(initial condition)은 $0$이라 하자. 그러면,
$$ \begin{align*} Y & = \frac{1}{s^2 + 3s + 2}R \\[5px] & = \left( \frac{1}{s+1} - \frac{1}{s+2} \right) R =: QR \end{align*} $$
이 때, $q = \mathcal{L}^{-1}(Q) = e^{-t} - e^{-2t}$ 이므로,
$$ \begin{align*} y & = q \ast r \\[5px] & = \int_{0}^{t} q(t-\tau) r(\tau) \, d\tau \\[5px] & = \int_{0}^{t} e^{-(t-\tau)} - e^{-2(t-\tau)} r(\tau) \, d\tau \end{align*} $$
이제, $t$의 범위에 따라 세 가지 경우를 고려해야 한다.
- $t<1$인 경우, $y=0$.
- $1<t<2$인 경우, $$ \begin{align*} y & = \int_{1}^{t} e^{-(t-\tau)} - e^{-2(t-\tau)} \, d\tau \\[5px] & = e^{-(t-\tau)} - \frac{1}{2}e^{-2(t-\tau)} \bigg|_{1}^{t} \\[5px] & = \frac{1}{2} + e^{-(t-1)} - \frac{1}{2}e^{-2(t-1)} \end{align*} $$
- $t>2$인 경우, $$ \begin{align*} y & = \int_{1}^{2} e^{-(t-\tau)} - e^{-2(t-\tau)} \, d\tau \\[5px] & = e^{-(t-2)} - e^{-(t-1)} - \frac{1}{2}e^{-2(t-2)} + \frac{1}{2}e^{-2(t-1)} \end{align*} $$