# 2차 미분방정식의 라플라스 해법

#### 2차 선형 미분방정식의 일반 형태

초등적이면서도 매우 중요한 미분방정식의 한 예로서 2차 선형 미분방정식을 들 수 있다. 미분방정식이 선형이고 계수들이 상수라고 할 때, 일반적으로 다음과 같은 형태를 생각할 수 있다

$$
y'' + a\_1 y' + a\_2 y = f(t)
$$

이때 초기조건으로서 $y(0) = y\_0$, $y'(0) = v\_0$ 등이 주어지면, 이 미분방정식의 해를 완전히 특정할 수 있다. 라플라스 변환은 이러한 초기조건이 주어진 2차 상수계수 선형 미분방정식을 효과적으로 해결하는 한 가지 중요한 방법이다.

#### 초기조건 문제와 라플라스 변환의 기본 아이디어

라플라스 변환을 통해 미분방정식을 해석하는 관건은, 시간영역에서의 미분 연산을 주파수 영역(또는 복소 영역)에서 곱셈으로 단순화하는 것이다. 특히 라플라스 변환에서

$$
\mathcal{L}{y'(t)}(s) = sY(s) - y(0),
\\
\mathcal{L}{y''(t)}(s) = s^2 Y(s) - s y(0) - y'(0),
$$

와 같은 관계가 성립하므로, 미분항을 가진 방정식이 $Y(s)$(즉, 미지함수 $y(t)$의 라플라스 변환)만을 포함하는 대수적 방정식으로 바뀌게 된다. 이 대수 방정식을 풀고, 다시 역라플라스 변환을 취하면 시간영역에서의 해 $y(t)$를 구할 수 있다.

#### 2차 방정식의 라플라스 변환 적용 과정

상술한 일반식

$$
y''(t) + a\_1 y'(t) + a\_2 y(t) = f(t)
$$

에 라플라스 변환을 취해보자. $y(0)=y\_0$, $y'(0)=v\_0$라고 할 때, 선형성과 초기조건을 반영하여 변환하면

$$
s^2 Y(s) - s y\_0 - v\_0 + a\_1 \bigl(sY(s) - y\_0 \bigr) + a\_2 Y(s) = F(s),
$$

여기서 $Y(s) = \mathcal{L}{y(t)}(s)$, $F(s) = \mathcal{L}{f(t)}(s)$이다. 위 식을 $Y(s)$에 대하여 정리하면

$$
\bigl(s^2 + a\_1 s + a\_2\bigr) Y(s) - s y\_0 - v\_0 - a\_1 y\_0 = F(s).
$$

곧

$$
Y(s) = \frac{F(s) + s y\_0 + v\_0 + a\_1 y\_0}{s^2 + a\_1 s + a\_2}.
$$

이는 $Y(s)$에 대한 해이며, 결국 다음 단계는 $Y(s)$를 적절히 분해하고 역라플라스 변환을 취하여 $y(t)$를 구하는 과정으로 귀결된다.

#### 특성방정식과 해의 구조

상수계수 선형 미분방정식을 해석할 때에는, 대응하는 특성방정식을

$$
r^2 + a\_1 r + a\_2 = 0
$$

라 쓰고, 그 근의 분포에 따라 다음과 같은 세 가지 형식적 해를 얻는 사실을 미리 알고 있다.

실제로 라플라스 해법을 적용해도 결과적으로 같은 형태가 도출된다. 즉 특성근 $r\_1, r\_2$가 실수이거나 복소수일 때, 해는 지수함수 혹은 진동함수(사인, 코사인) 형태로 나타난다.

라플라스 변환을 통한 해석에서도, 분모가 $s^2 + a\_1 s + a\_2$인 표현을 부분분수로 분해하거나, $s$ 평면에서의 극점 구조를 해석함으로써 결과가 $e^{\alpha t}$나 $e^{\alpha t}\sin(\beta t)$, $e^{\alpha t}\cos(\beta t)$ 형태로 나타난다. 따라서 라플라스 해법은 특성방정식 해법과 동일한 분석적 구조를 가지면서도, 초기값 문제가 있는 실제 공학적 상황에서 계산하기가 더 편리할 수 있다.

#### 예시: 강제진동계 문제

물리학이나 공학에서 빈번히 등장하는 예로, 질량-스프링-댐퍼 시스템의 운동방정식을 꼽을 수 있다. 마찰계수(또는 감쇠계수)를 $c$, 스프링상수를 $k$, 질량을 $m$이라 하면, 뉴턴의 운동법칙에 의해 주어지는 미분방정식은

$$
m y''(t) + c y'(t) + k y(t) = F(t)
$$

와 같은 형태가 된다. 이를 표준화하기 위해 양변을 $m$으로 나누면

$$
y''(t) + \frac{c}{m}y'(t) + \frac{k}{m}y(t) = \frac{F(t)}{m}
$$

가 된다. 이때 $\displaystyle a\_1 = \frac{c}{m}$, $\displaystyle a\_2 = \frac{k}{m}$, 그리고 $f(t) = \frac{F(t)}{m}$라 할 수 있다. 위 식을 라플라스 변환 기법으로 풀면, 초기속도나 초기변위를 가정했을 때 진동계의 움직임이 정해진다.

#### 해석 절차의 시각화

아래는 2차 상수계수 선형 미분방정식을 라플라스 변환으로 해석하는 과정을 단순화한 흐름도이다.

{% @mermaid/diagram content="flowchart LR
A(("2차 미분방정식")) --> B\[라플라스 변환하기]
B --> C\[변형된 대수 방정식 풀기]
C --> D\[역라플라스 변환]
D --> E(("최종 해 y(t)"))" %}

시간영역에서 다루기 까다로운 미분 연산은, 라플라스 영역에서 변환된 대수적 조작으로 간단히 나타난다. 그 후 역변환만 잘 취하면 시간영역에서의 완전한 해를 구할 수 있다. 이러한 과정을 체계적으로 익혀두면, 점점 복잡한 형태의 2차 미분방정식도 손쉽게 해석할 수 있다.

#### 부분분수 분해의 중요성

분수 형태로 나타난 $Y(s)$를 역라플라스 변환하기 위해서는 부분분수 분해(partial fraction decomposition)가 중요하게 사용된다. 예를 들어,

$$
Y(s) = \frac{F(s)}{s^2 + a\_1 s + a\_2}
$$

형태가 주어졌다고 하자. 특성다항식 $s^2 + a\_1 s + a\_2$가 서로 다른 두 실근을 갖는 경우, 다음과 같은 방식으로 분해할 수 있다

$$
\frac{F(s)}{(s-r\_1)(s-r\_2)} = \frac{A}{s-r\_1} + \frac{B}{s-r\_2}.
$$

복소근이 존재할 때는 묶음 형태의 분해를 고려하여 사인과 코사인 항을 얻는다. 구체적인 분해 과정은 $F(s)$의 형태에 따라 크게 달라지므로, 실제 계산 상황에 맞게 적절한 기법을 사용한다.

#### 초기조건 항의 반영

라플라스 변환에서 미분항이

$$
\mathcal{L}{y''}(s) = s^2Y(s) - sy(0) - y'(0)
$$

와 같이 전개되면, 분수식에 $s y(0)$, $y'(0)$, $a\_1 y(0)$ 등이 추가로 나타난다. 이들은 역라플라스 변환을 통해 시간영역에서 초기치에 해당하는 지수함수나 델타함수 등의 조합으로 다시 표현된다. 따라서 해 $y(t)$는 항상 균일해(즉, 제로 초기조건일 때의 해)와 초기치에 의한 자유응답의 합으로 구분해서 볼 수 있다. 자세히 쓰면

$$
Y(s) = \frac{F(s)}{s^2 + a\_1 s + a\_2} + \text{(초기조건 항에 대응되는 부분)}
$$

이 되고, 역라플라스 변환 과정에서 이 항들을 분리해서 해석하면, 초기치에 의해 형성되는 자연응답과 강제응답을 좀 더 명확히 파악할 수 있다.

#### 라플라스 변환 해법의 장점

라플라스 변환을 통한 2차 미분방정식 풀이가 갖는 장점은 크게 두 가지이다.

첫째, 미분연산이 라플라스 영역에서 단순히 $s$와의 곱셈 및 초기조건 보정항으로 표현되기 때문에, 복잡한 적분 및 미분 계산을 대수적 조작으로 바꿀 수 있다.

둘째, 초기조건이 자연스럽게 식에 반영되어, 별도의 결정을 해야 하는 번거로운 단계 없이도 최종해에서 초기조건이 만족되도록 자동 조절된다.

#### 일반해의 구성

2차 미분방정식이기 때문에, 시간영역에서의 해 $y(t)$는 언제나 두 개의 선형독립해들의 조합으로 표현된다. 동일한 표현을 라플라스 공간에서도 선형성의 원리를 이용해 증명할 수 있다. 특성근이 단순근 두 개이면 지수함수 두 개, 중근이면 지수함수와 $t$가 곱해진 항, 복소근이면 감쇠 혹은 증가하는 사인과 코사인 형태가 얻어진다. 결국 최종적으로는

$$
y(t) = y\_h(t) + y\_p(t)
$$

꼴로서, 균일해 $y\_h(t)$와 특정해 $y\_p(t)$를 합친 형태임을 알 수 있다. 이는 라플라스 변환에서 분모인 특성다항식과 분자 $F(s)$를 분해하는 과정에서 자연스럽게 분리된다.

#### 비균일 항의 역라플라스 변환 기법

일반적으로 2차 미분방정식에서 비균일 항 $f(t)$가 존재할 때, 그 라플라스 변환 $F(s)$의 형태가 단순하지 않을 수 있다. 예컨대 $f(t) = e^{\alpha t}$, $\sin(\omega t)$, $\delta(t - T)$ 등 다양한 함수가 올 수 있으며, 이들에 대한 라플라스 변환은 다음과 같은 표준 식을 활용할 수 있다

$$
\mathcal{L}{ e^{\alpha t} }(s) = \frac{1}{s-\alpha}
\quad
\text{(단, } s>\alpha \text{)}
\\
\mathcal{L}{ \sin(\omega t) }(s) = \frac{\omega}{s^2 + \omega^2}
\\
\mathcal{L}{ \delta(t - T) }(s) = e^{-Ts}
$$

이들을 문제에 맞게 대입하여 부분분수를 시도한 뒤 역라플라스 변환을 취하면, 비균일 항에 대응되는 특정해( particular solution )를 얻는다.

$f(t)$가 특히 복잡한 조합으로 주어졌을 때에는, 합성 정리를 이용하여 적분형태의 해를 구성하기도 한다. 예를 들어 $f(t)$의 라플라스 변환이 단순 분수 형태로 나타나지 않는다면, 부분분수 분해가 곤란할 수 있는데, 이때 컨볼루션(합성곱) 정리를 활용하면 간결하게 표현될 수 있다. 1차, 2차를 막론하고 선형 미분방정식을 해결할 때, 이 컨볼루션 정리는 매우 중요한 역할을 한다.

#### 컨볼루션 정리와 그 적용

컨볼루션 정리는

L{f∗g}(s)=F(s)G(s)\mathcal{L}{f\*g}(s) = F(s) G(s)

와 같은 형태로 잘 알려져 있다. 여기서 $f\*g$는

$$
(f\*g)(t) = \int\_{0}^{t} f(\tau) g(t-\tau), d\tau
$$

를 의미한다. 만약 비균일 항이 라플라스 공간에서 곱셈으로 표현된다고 하면, 시간영역에서는 이 합성곱 형태로 나타난다.

2차 미분방정식 해를 $y(t)$라고 할 때, 균일해 $y\_h(t)$에 더하여 특정해 $y\_p(t)$가 컨볼루션으로 표현되는 경우가 종종 발생한다. 예를 들어

$$
y\_p(t) = (h\*g)(t)
$$

같은 식이 얻어질 수 있는데, 여기서 $h(t)$는 대응하는 계의 임펄스응답(즉, 그린 함수에 대응)이고, $g(t)$는 비균일 항 $f(t)$에 관련된 표현이다. 이는

$$
Y(s) = H(s) \cdot F(s)
$$

에서 기인하는데, $H(s)$를 시스템(계수 $a\_1, a\_2$로 특징지어지는 2차 시스템)의 전달함수로 간주하면, $F(s)$는 외부 구동 항에 대한 라플라스 변환으로 해석될 수 있다. 물리학이나 공학의 시스템 해석 관점에서 매우 중요한 통찰점이며, 임펄스응답과 컨볼루션 연산의 일반성이 라플라스 변환의 유용성을 뒷받침한다.

#### 델타함수 구동항과 충격응답

비균일 항이 델타함수(Dirac delta) 형태일 때, 예를 들어 $f(t) = \delta(t - T)$와 같은 형태가 주어지면, 라플라스 공간에서는 $F(s) = e^{-Ts}$가 된다. 이 경우

$$
Y(s) = \frac{e^{-Ts}}{s^2 + a\_1 s + a\_2}
$$

형태를 얻게 되며, 이를 부분분수로 분해하거나 고유모드 해석을 적용하면 시간지연 $T$가 반영된 충격응답(impulse response)을 구할 수 있다. 시간영역에서 살펴보면, $t=T$ 시점에 계에 순간적인 충격이 가해지고 그 이후로 시스템이 (자유응답 형태로) 진동하거나 지수적으로 감쇠하는 양상을 보인다.

덧붙여, $\delta(t)$ 자체가 초기조건과는 무관하게 특정 시점에서 계를 강제하는 효과이므로, 이것이 라플라스 변환에서 $e^{-Ts}$와 연결된다는 사실은 인과성(causality)의 관점에서도 해석이 가능하다.

#### 중근을 갖는 경우의 해

2차 미분방정식의 특성방정식이 $r^2 + a\_1 r + a\_2 = 0$인데, 판별식이 0이 되어 한 근 $r\_0$만 중복되면, 균일해는

$$
y\_h(t) = C\_1 e^{r\_0 t} + C\_2,t,e^{r\_0 t}
$$

의 형태가 된다. 이 사실은 지수함수 두 개가 동일하므로, 두 번째 독립해를 얻기 위해 $t$를 곱해주어야 함을 미분방정식 이론에서 배운다. 라플라스 해법을 적용하면,

$$
s^2 + a\_1 s + a\_2 = (s - r\_0)^2
$$

이 되므로, $Y(s)$가 $(s-r\_0)^2$ 형태를 분모에 포함하는 양상의 부분분수식을 가지게 된다. 즉

$$
\frac{A}{s-r\_0} + \frac{B}{(s-r\_0)^2}
$$

와 같은 식으로 분해된다. 이때 역라플라스 변환을 취하면 $A e^{r\_0 t} + B, t, e^{r\_0 t}$ 형태가 자연스럽게 도출된다. 초기조건이나 비균일 항이 존재할 때에도, 이 중복 근의 경우에는 $t$가 곱해진 항이 반드시 등장한다는 점에 유의해야 한다.

#### 복소근을 갖는 경우의 해

판별식이 음수가 되어 복소근 $r\_{1,2} = \alpha \pm i\beta$를 가지면, 미분방정식의 균일해는

$$
e^{\alpha t}\bigl(C\_1 \cos(\beta t) + C\_2 \sin(\beta t)\bigr)
$$

의 꼴이 된다. 라플라스 공간에서는 분모가 $s^2 + a\_1 s + a\_2$를 완전히 제곱꼴로 묶거나, 혹은 $s - \alpha$와 $(s - \alpha)^2 + \beta^2$ 같은 형태로 재배열하는 방식으로 처리한다. 예를 들어

$$
s^2 + 2\alpha s + (\alpha^2 + \beta^2) = (s + \alpha)^2 + \beta^2
$$

와 같은 식으로 표현되어, 부분분수 분해에서 $\frac{1}{(s+\alpha)^2 + \beta^2}$, $\frac{s+\alpha}{(s+\alpha)^2 + \beta^2}$ 등으로 나뉘게 된다. 이 항들의 역라플라스 변환이 지수함수와 사인·코사인 함수를 복합적으로 결합하는 형태로 연결되어, 위와 같은 시간영역 해를 준다.

비균일 항이 존재하면, 해당 항도 마찬가지 방식으로 분해 및 역변환되어 최종 해에 사인이나 코사인 항이 포함될 수 있다.

#### 라플라스 변환에서의 초기조건 해석

초기조건이 $y(0) = y\_0$, $y'(0) = v\_0$로 주어졌을 때, 라플라스 식에 $-s y\_0 - v\_0$ 등으로 나타나는 것을 종종 ‘초기값 편차’라고 부른다. 이는

$$
\mathcal{L}{y''}(s) = s^2Y(s) - s y(0) - y'(0)
$$

에서 기원하는 것으로서, 시간영역에서 미분 연산이 0부터 시작되는 적분 형태와 연동되어 있음을 시사한다.

또한, 실제 물리계 해석에서는 초기조건을 별도로 설정하기보다는, $t=0$ 이전에 시스템이 받았던 충격이나 운동이 $t=0$에 이르는 동안의 효과로서 반영된다고 볼 수 있다. 라플라스 도메인에서는 이 모든 것을 단순 대수항으로 옮겨 놓을 수 있으므로, 초기상태 문제를 해결하기가 훨씬 수월하다.

#### 고차 미분방정식으로의 확장

2차 미분방정식에 대한 라플라스 변환 접근은 고차(3차, 4차 등) 미분방정식으로 쉽게 확장된다. $n$차 선형 상수계수 미분방정식

$$
y^{(n)} + a\_{1}y^{(n-1)} + \cdots + a\_{n}y = f(t)
$$

에서도, 라플라스 변환은

$$
s^n Y(s) - s^{n-1}y(0) - s^{n-2}y'(0) - \cdots - y^{(n-1)}(0)
$$

과 같은 구조를 형성한다. 따라서 2차인 경우와 마찬가지로 $Y(s)$에 관한 대수 방정식을 풀고, 역라플라스를 구하면 해를 얻을 수 있다. 이러한 일관성 덕분에, 라플라스 기법은 높은 차수의 미분방정식 문제에서도 주로 활용된다.

#### 단계별 요약의 부재와 직관

2차 미분방정식을 풀 때, 흔히 전형적인 풀이 절차(균일해-특정해-상수 결정)를 단계적으로 정리한다. 라플라스 변환은 이러한 단계에 일일이 이름을 붙이지 않고, 변환-대수식 해법-역변환으로 단순화한다. 그러나 본질적으로는 균일해와 특정해를 모두 포함하며, 각각의 해가 최종적으로 어떻게 결정되는지 추적하면, 전통적 해법과 동일한 논리임을 알 수 있다.

다만, 복잡한 비균일 항이 주어졌을 때 라플라스 변환 쪽이 훨씬 빠르고 편리한 계산 과정을 제공한다. 특히, 초기조건을 설정하거나 시간 지연, 충격, 주기적 신호가 섞인 문제의 경우에도, 라플라스 도메인에서 단순 대수 조작만 잘 수행하면 되므로 일반적인 풀이가 용이하다.

#### 라플라스 변환에서의 통합인자 비교

전통적인 미분방정식 풀이에서는 1차 방정식에 통합인자를 쓰고, 2차 이상의 경우에는 그린 함수나 미분연산자의 역연산 형태 등으로 확장한다. 라플라스 변환 기법은 이 통합인자(또는 그린 함수)를 자동적으로 만들어내는 셈이다. 실제로 2차 계에 대해 시스템함수 혹은 임펄스응답을 구해놓으면, 외부 입력을 이 임펄스응답과 컨볼루션하여 $y\_p(t)$를 얻을 수 있다.

이것은 전통적 방법에서의 특정해 구하기와 본질적으로 동일하며, 오히려 라플라스 방법이 보다 투명하고 기계적인 계산 절차를 제시한다.

#### 예시 해법: 구체적인 문제 풀이 과정

이번에는 실제 예시를 들어 2차 미분방정식을 라플라스 변환으로 푸는 전 과정을 좀 더 자세히 살펴보자. 초기조건이 있는 비균일 미분방정식을 다뤄보면, 앞서 언급된 여러 개념들이 어떻게 맞물려 작용하는지 확인할 수 있다.

예를 들어 다음과 같은 미분방정식을 생각한다

$$
y''(t) + 3y'(t) + 2y(t) = e^{-t},
$$

초기조건으로서

$$
y(0) = 1,\quad y'(0) = 0
$$

이 주어졌다고 하자. 이 문제를 라플라스 변환 기법으로 풀어본다.

먼저, 양변에 라플라스 변환을 취한다. $Y(s) = \mathcal{L}{y(t)}(s)$를 쓰면, 선형성과 초기조건 반영으로 다음 식을 얻는다

$$
\mathcal{L}{y''(t)}(s) = s^2 Y(s) - s y(0) - y'(0),
\\
\mathcal{L}{y'(t)}(s) = s Y(s) - y(0),
\\
\mathcal{L}{y(t)}(s) = Y(s),
\\
\mathcal{L}{e^{-t}}(s) = \frac{1}{s+1}.
$$

따라서 미분방정식

$$
y''(t) + 3y'(t) + 2y(t) = e^{-t}
$$

는 라플라스 영역에서

$$
\bigl(s^2 Y(s) - s y(0) - y'(0)\bigr) + 3 \bigl(s Y(s) - y(0)\bigr) + 2 Y(s) = \frac{1}{s+1}
$$

로 표현된다.

이를 초기조건 $y(0)=1$, $y'(0)=0$에 대입하면,

$$
s^2 Y(s) - s \cdot 1 - 0 + 3\bigl(sY(s) - 1\bigr) + 2Y(s) = \frac{1}{s+1}.
$$

즉

$$
s^2 Y(s) - s + 3sY(s) - 3 + 2Y(s) = \frac{1}{s+1}.
$$

이제 $Y(s)$에 관한 항들을 모으고 정리한다

$$
\bigl(s^2 + 3s + 2\bigr) Y(s) - s - 3 = \frac{1}{s+1}.
$$

좌변에 등장한 다항식 $s^2 + 3s + 2$는 $(s+1)(s+2)$로 인수분해되므로,

$$
(s+1)(s+2),Y(s) - s - 3 = \frac{1}{s+1}.
$$

곧

$$
(s+1)(s+2),Y(s) = \frac{1}{s+1} + s + 3.
$$

오른쪽을 통분하거나 여러 방식으로 정리할 수 있다. 우선 간단히 분리한다

$$
(s+1)(s+2),Y(s) = s + 3 + \frac{1}{s+1}.
$$

이를 $(s+1)(s+2)$로 나누어 $Y(s)$를 구하면

$$
Y(s) = \frac{s + 3}{(s+1)(s+2)} + \frac{1}{(s+1)(s+1)(s+2)} \quad\text{(분자/분모 분할)}
$$

이와 같이 표현된다.

첫 항

$$
\frac{s + 3}{(s+1)(s+2)}
$$

는 이미 분자 차수가 1이므로, 이를 부분분수 분해하면 된다. 예를 들어

$$
\frac{s + 3}{(s+1)(s+2)} = \frac{A}{s+1} + \frac{B}{s+2}.
$$

이때 $A(s+2) + B(s+1) = s+3$ 식을 비교하거나 $s$의 계수 및 상수항을 맞추면 $A$와 $B$를 찾을 수 있다. 실제로

$$
s+3 = (A+B),s + (2A + B).
$$

계수를 비교하면

$$
A + B = 1, \quad 2A + B = 3.
$$

연립하여 풀면 $A=1$, $B=0$을 얻게 된다. 따라서 첫 항은

$$
\frac{s + 3}{(s+1)(s+2)} = \frac{1}{s+1}.
$$

즉, 첫 항이 단순히 $\displaystyle \frac{1}{s+1}$임이 드러난다.

두 번째 항

$$
\frac{1}{(s+1)(s+1)(s+2)}
$$

는 $(s+1)^2(s+2)$를 분모로 갖는다. 이를

$$
\frac{1}{(s+1)^2(s+2)} = \frac{C}{s+1} + \frac{D}{(s+1)^2} + \frac{E}{s+2}
$$

와 같이 부분분수 분해한다. 곧

$$
1 = C(s+1)(s+2) + D(s+2) + E(s+1)^2.
$$

이를 전개하고 계수를 비교하여 $C, D, E$를 구하면 된다. 전개하면

$$
C(s^2 + 3s + 2) + D(s+2) + E(s^2 + 2s + 1).
$$

보다 구체적으로 쓰면

$$
\= C s^2 + 3C s + 2C + D s + 2D + E s^2 + 2E s + E
$$

이고, 이를 $s^2$, $s$, 상수항 순으로 정리하면

$$
(s^2),(C + E) + s,(3C + D + 2E) + (2C + 2D + E).
$$

이것이 항등적으로 1과 같아야 하므로,

$$
C + E = 0,
\\
3C + D + 2E = 0,
\\
2C + 2D + E = 1.
$$

연립하여 해를 구해본다.

첫 식에서 $E = -C$.

이를 두 번째 식에 대입하면

$$
3C + D + 2(-C) = 0 \quad\Rightarrow\quad C + D = 0 \quad\Rightarrow\quad D = -C.
$$

세 번째 식에 대입하면

$$
2C + 2(-C) + (-C) = 1 \quad\Rightarrow\quad 2C - 2C - C = 1 \quad\Rightarrow\quad -C = 1,
$$

즉 $C = -1$이 된다. 그러면 $D = -(-1)=1$, $E = -(-1)=1$을 얻는다.

결국

$$
\frac{1}{(s+1)^2(s+2)} = \frac{-1}{s+1} + \frac{1}{(s+1)^2} + \frac{1}{s+2}.
$$

위 결과를 종합하면,

$$
Y(s) = \frac{1}{s+1} + \Bigl(\frac{-1}{s+1} + \frac{1}{(s+1)^2} + \frac{1}{s+2}\Bigr).
$$

간단히 묶어서 정리하면

$$
Y(s) = \frac{1}{(s+1)^2} + \frac{1}{s+2}.
$$

왜냐하면 $\displaystyle \frac{1}{s+1}$와 $\displaystyle \frac{-1}{s+1}$가 서로 상쇄되고, 남은 항들이 $\displaystyle \frac{1}{(s+1)^2} + \frac{1}{s+2}$가 된다.

따라서

$$
Y(s) = \frac{1}{(s+1)^2} + \frac{1}{s+2}.
$$

이제 이를 역라플라스 변환한다.

첫 항 $\displaystyle \frac{1}{(s+1)^2}$의 역변환은 $t e^{-t}$이고, 두 번째 항 $\displaystyle \frac{1}{s+2}$의 역변환은 $e^{-2t}$이다. 즉

$$
y(t) = \mathcal{L}^{-1}\Bigl{\frac{1}{(s+1)^2}\Bigr} + \mathcal{L}^{-1}\Bigl{\frac{1}{s+2}\Bigr} = t e^{-t} + e^{-2t}.
$$

이것이 초기조건 $y(0)=1,, y'(0)=0$를 만족하면서, $y'' + 3y' + 2y = e^{-t}$를 만족하는 미분방정식의 해다. 실제로 $t=0$에서 $y(0) = 0 + 1 = 1$, $y'(t) = (1\cdot e^{-t} + t(-e^{-t})) + (-2e^{-2t})$ 등을 대입하면 $y'(0)=0$임을 확인할 수 있다.

이 과정을 통해, 라플라스 변환을 쓰면 비균일 항과 초기조건을 자동으로 처리하면서도, 부분분수 분해 및 역변환 한 번으로 해를 구할 수 있다는 사실이 명확히 드러난다.

#### 계단함수(스텝함수) 구동의 예시

비균일 항이 헤비사이드 계단함수(Heaviside step function) $u(t - T)$ 형태를 가질 수도 있다. 예컨대

$$
f(t) = u(t - T) =  \begin{cases} 0, & t < T, \ 1, & t \ge T \end{cases}
$$

라 하자. 그 라플라스 변환은 $F(s)=\frac{e^{-Ts}}{s}$가 된다. 즉 $t=T$ 이후부터 1인 상수함수이므로, 라플라스 공간에서는 $e^{-Ts}$ 형태가 붙는다. 이를 2차 미분방정식의 우변에 대입하고, $Y(s)$를 구한 뒤 역변환을 취하면, $t\<T$ 구간에서는 0이었던 구동력이 $t=T$ 이후에 작용하기 시작하는 해석적 구조가 도출된다.

이를 시간영역에서 직접 풀면 구간별로 해를 달리 구하고 적절한 매칭 조건을 적용해야 하지만, 라플라스 변환에서는 단순히 $F(s)=\frac{e^{-Ts}}{s}$를 이용하여 부분분수 분해하고, $e^{-Ts}$가 역변환에서 시간지연 연산으로 나타난다는 점만 고려하면 손쉽게 해를 얻을 수 있다.

#### 시뮬레이션 관점에서의 라플라스 해석

공학적 시뮬레이션(예: 전기회로 해석, 기계계 해석)에서는, 2차 선형 미분방정식은 시스템의 상태공간 표현을 직접 시뮬레이션하거나, 주파수 응답 함수를 사용하여 각종 입력에 대한 시스템 출력을 분석한다. 라플라스 변환은 이 두 관점 모두에 다리를 놓아주는 도구다.

시스템 전달함수 $G(s)=\frac{Y(s)}{U(s)}$를 정의한다면, $U(s)$가 주어진 입력(비균일 항의 라플라스 변환)일 때 $Y(s)=G(s),U(s)$로 표현되며, 이후 역라플라스 변환을 통해 시간영역 해를 얻는다. 이는 임펄스응답 $g(t)$와의 컨볼루션 $y(t) = (g\*u)(t)$로 재해석되며, $g(t) = \mathcal{L}^{-1}{G(s)}$가 그린 함수(그 계의 임펄스응답)로서 활용된다.

#### 미분연산자의 다항식 해석

2차 미분방정식을 간결하게 쓸 때, 다음과 같은 ‘미분연산자’ 표기도 자주 쓴다

$$
(D^2 + a\_1 D + a\_2),y(t) = f(t).
$$

여기에서 $D=\frac{d}{dt}$로 표기한다. 이를 라플라스 변환하면 $D$가 $s$로, $D^2$가 $s^2$로 대응된다고 볼 수 있다. 즉,

$$
\mathcal{L}{(D^2 + a\_1 D + a\_2),y(t)}(s) = (s^2 + a\_1 s + a\_2), Y(s) - \text{(초기조건 항)}.
$$

이처럼 미분연산자 다항식을 $s$에 대한 다항식으로 치환하는 방식은, 라플라스 변환이 고차 미분방정식을 풀 때 자동으로 적용되는 ‘연산자 대수’의 예라 할 수 있다. 결과적으로, 시간 영역에서의 어려운 미분·적분 조합을 $s$ 영역에서 곱셈·나눗셈 문제로 치환해주는 점이, 라플라스 변환 해법의 핵심이다.

#### 물리적 적용과 자유응답·강제응답 분리

2차 시스템에 외력이 작용하는 상황에서, 해 $y(t)$는 초기조건에 의한 자유응답(free response)과 비균일 항에 의한 강제응답(forced response)이 합쳐진 형태로 나타난다. 라플라스 변환으로 구한 $Y(s)$에서 초기조건 항과 $F(s)$ 항을 각각 분리해 보면, 자유응답과 강제응답을 구분할 수 있다.

예컨대, 동일한 시스템

$$
y'' + 3y' + 2y = 0
$$

에 대한 해(즉 $f(t)=0$)를 구하면, 이것이 자유응답에 해당한다. 특정 시점 이후에 외력이 주어진다면(예: $u(t-T)$ 또는 $\delta(t-T)$), 그에 따른 강제응답 항이 라플라스 영역에서 자동 분리되어 나타난다.

이처럼 라플라스 변환을 통해 자유응답은 특성방정식의 근(계의 고유모드)로부터 만들어지며, 강제응답은 비균일 항의 변환으로부터 생성된다는 구조적 사실이 더욱 명료해진다.

#### 주파수 영역 해석과 관계

라플라스 변환의 변수 $s = \sigma + i\omega$를 복소 주파수라 할 때, 이 시스템이 가진 특성근(즉, 분모를 0으로 만드는 극점)이 곧 $s$-평면 상에서의 극점이 된다. 2차 시스템의 극점은 실축에 있는 경우 감쇠형 지수, 좌반평면 복소영역에 있는 경우 감쇠진동, 우반평면 복소영역에 있는 경우 발산진동 등을 예측케 한다.

실제로, 선형 2차 계의 전달함수 $G(s)=\frac{Y(s)}{U(s)}$의 극점이 곧 특성근이므로, 시스템의 안정성이나 응답특성을 조사할 때 라플라스 해법이 직접적인 이론적 토대가 된다. 시간영역 해석과 주파수영역 해석이 이처럼 긴밀히 연결되어 있다는 점이, 2차 미분방정식(또는 일반적인 $n$차 미분방정식) 해석을 단순한 상미분방정식 풀이 이상의 의미로 끌어올린다.

#### 전기회로 해석과 RLC 회로 모델

라플라스 변환은 전기회로 해석에서도 매우 유용한 도구로 활용된다. 대표적인 예가 RLC 회로다. 저항(resistor) $R$, 인덕터(inductor) $L$, 커패시터(capacitor) $C$로 이루어진 직렬 회로를 생각하면, 커패시터 양단 전압이나 인덕터 전류, 혹은 특정 지점의 전압과 전류 사이에 2차 상미분방정식이 성립한다. 예를 들어, 직렬 RLC 회로에서 커패시터 전압을 $v\_C(t)$라 하고, 외부에서 공급되는 전압원을 $v\_{\text{in}}(t)$라고 할 때,

$$
\frac{d^2 i(t)}{dt^2} ;+; R,\frac{d i(t)}{dt} ;+; \frac{1}{C},i(t) ;=; \frac{d v\_{\text{in}}(t)}{dt}
$$

등의 형태의 방정식을 얻을 수 있다. 인덕터와 커패시터의 상호작용으로 인해 $i(t)$나 $v\_C(t)$가 2차 미분항을 포함하게 되고, 이때 초기조건은 커패시터에 축적된 전하나 인덕터에 흐르던 전류 등으로 설정된다.

라플라스 변환을 쓰면, 인덕터의 임피던스가 $sL$, 커패시터의 임피던스가 $\frac{1}{sC}$로 대응되며, 회로방정식이 대수적으로 변환된다. 전기회로 과목에서 빈번히 사용하는 고전적 임피던스 해석이 사실상 라플라스 해법과 일맥상통한다. 이를 통해 회로 응답을 간단히 구하거나, 공진주파수와 댐핑 특성을 추적하기가 수월해진다.

특히 2차 시스템인 RLC 회로는 전압과 전류가 지수함수적(또는 지수함수와 사인·코사인 결합의) 응답을 보이며, 판별식에 따라 저댐핑(underdamped), 임계댐핑(critically damped), 과댐핑(overdamped) 등의 양상이 결정된다. 전기회로에서 저댐핑은 고속의 링잉(ringing)과 진동현상을 야기하고, 과댐핑은 반응이 느리지만 진동 없이 점진적으로 수렴한다. 임계댐핑은 진동 없이 가장 빠른 속도로 최종값에 도달하는 경우다.

라플라스 해석에서는 $s$-영역에서 극점(pole)의 위치로 이 현상을 정확히 파악할 수 있다. 극점이 복소공액쌍으로 좌반평면에 있으면 저댐핑, 극점이 실수 중복근이면 임계댐핑, 서로 다른 실근 두 개가 음수이면 과댐핑 등으로 분류한다.

이에 따라 초기조건 또는 스텝 입력 등에 대해 커패시터 전압 $v\_C(t)$나 인덕터 전류 $i\_L(t)$를 간단히 표현할 수 있으며, 특정 시점 이후에 가해지는 펄스 입력, 스위치 동작에 따른 전압 전류 변화 등도 $e^{-Ts}$ 형태의 지연항으로 모델링하여 해를 구할 수 있다.

#### 기계계 해석과 해석적 유사성

질량-스프링-댐퍼 시스템이 보이는 2차 운동방정식과 RLC 직렬회로가 보이는 2차 전기회로방정식은 구조가 매우 유사하다. 실제로 전기회로에서는 $L, R, C$가 각각 관성, 마찰(감쇠), 탄성에 대응된다고 볼 수 있다.

이런 대응관계 하에서, 라플라스 공간에서 해를 구하는 과정은 전기·기계 시스템 할 것 없이 같은 원리로 작동한다. $m y'' + c y' + k y = F(t)$를 풀 때와, $L i'' + R i' + \frac{1}{C} i = \dots$ 식을 풀 때가 형식적으로 유사하므로, 한 쪽 해석 기법을 다른 쪽에도 그대로 적용할 수 있다. 이는 물리학에서 말하는 ‘모델링의 유사성(analogy)’이 라플라스 변환을 매개로 간단히 표현되는 한 예라 할 수 있다.

#### 감쇠와 공진

2차 선형계에서 매우 중요한 현상 중 하나가 공진(resonance)이다. 외부로부터의 주기적 구동(예: $\sin(\omega t)$)이 계의 고유진동수(또는 고유진동수 부근)에 가까울 때, 응답이 크게 증폭될 수 있다. 라플라스 해법에서도, 분모의 극점과 구동 주파수의 근접 여부에 따라 역라플라스 변환에서 사인·코사인 항이 증폭되는 상황이 드러난다.

만약 감쇠항 $a\_1 y'(t)$이 작다면, 공진 때 응답 진폭이 상당히 커지면서 긴 과도진동이나 큰 진폭의 진동이 발생한다. 전자공학적으로는 LC 공진, 기계공학적으로는 spring-mass 계의 공진이 대표 사례다. 공진을 방지하거나 조절하기 위해 감쇠를 적절히 설정하는 것이 매우 중요하며, 라플라스 해석을 통해 공진주파수 주변의 주파수응답 곡선을 분석하여 최적의 파라미터를 찾기도 한다.

#### 안정성과 극점의 위치

2차 선형계의 안정성은 특성방정식 근(즉, 라플라스 영역에서의 극점)의 실부가 음수인지 여부로 결정된다. 모든 근의 실부가 음수이면, 계는 시간에 따라 무한히 진동하거나 발산하지 않고, 결국 평형점으로 수렴한다. 만약 실부가 양수인 근이 하나라도 존재하면, 계는 지수적으로 발산하여 불안정해진다.

2차 계에서는 판별식에 따라 근의 위치가 다양하게 달라지며, 그 중에서도 실근이 0에 정확히 위치하거나 양수가 되면 한계안정, 불안정을 각각 의미한다. 즉, 감쇠나 저항이 전혀 없다면 에너지를 계속 보존하거나 증폭해 진동이 꺼지지 않는다. 라플라스 해석에서는 $s$-평면에서의 극점(근) 배치를 한눈에 볼 수 있어, 안정성 판단 및 응답특성 이해가 수월하다.

#### 미분방정식 해와 전달함수의 연결

2차 상수계수 미분방정식을

$$
y'' + a\_1 y' + a\_2 y = f(t)
$$

로 쓸 때, $F(s) = \mathcal{L}{f(t)}(s)$, $Y(s) = \mathcal{L}{y(t)}(s)$라 하면,

$$
(s^2 + a\_1 s + a\_2)Y(s) = F(s) + \dots
$$

형태가 되어, 시스템 전달함수(transfer function)를

$$
G(s) = \frac{Y(s)}{F(s)}
$$

로 정의할 수 있다. 계단함수, 지수함수, 사인 등 다양한 입력에 대한 해를 얻으려면, 이 전달함수에 해당 입력의 라플라스 변환을 곱한 뒤 역변환을 구하면 된다.

이 방식은 제어공학에서 주로 사용하는 블록선도 기법과 동일한 사고방식을 제공한다. 2차 미분방정식으로 표현되는 물리계(또는 회로)는 $G(s)$로 특징지어지고, 그 입력이 $F(s)$로 주어지면, 출력이 $Y(s)=G(s),F(s)$가 된다. 라플라스 변환은 이 과정의 기저이며, 시간영역 응답을 구하고자 할 때는 역변환만 수행하면 되므로, 입력이 달라져도 $G(s)$만 한 번 구해두면 응답 해석이 손쉽게 이루어진다.

#### 선형성과 중첩의 원리

2차 미분방정식이 선형이므로, 입력함수 여러 개가 동시에 작용할 때 전체 해는 각 입력이 따로 작용했을 때의 해를 모두 더한 것과 동일하다. 라플라스 변환은 이 선형성을 매우 간단히 보장해준다. 실제로

$$
\mathcal{L}{f\_1(t) + f\_2(t)} = F\_1(s) + F\_2(s),
\\
G(s),(F\_1(s) + F\_2(s)) = G(s),F\_1(s) + G(s),F\_2(s).
$$

역변환 과정에서도 동일한 선형성으로 인해, 각각의 역변환을 합치면 최종 해가 얻어진다. 물리나 공학에서 복잡한 입력신호를 단순 신호들의 합으로 분해하는 슈퍼포지션(superposition) 개념이, 라플라스 해석에 자연스레 녹아 있다.

#### 불연속 입력과 구간별 해석

구동신호가 시각별로 달라지는 piecewise 함수이거나, 특정 구간에서만 존재하는 펄스 신호, 주기적으로 반복되는 신호 등이 주어질 수도 있다. 이때 시간영역에서 직접 미분방정식을 구간별로 풀고 경계조건을 맞추는 방식은 복잡하다. 하지만 라플라스 변환에서는 스위치가 켜지는 순간이나 함수가 바뀌는 지점을 ‘시프트 연산’과 ‘계단함수’ 등으로 모델링할 수 있어서, 한 번에 대수 방정식으로 해석 가능해진다.

대표적 예가 스텝 입력 $u(t - T)$나 사각파(square wave)와 같은 신호이며, 라플라스 도메인에서 지연항 $e^{-Ts}$나 주기적 확장에 따른 적절한 형식($\frac{1}{1 - e^{-Ts}}$ 류)이 잘 정리되어 있으면, 간단히 곱셈으로 표현해낼 수 있다. 역변환도 부분분수 분해와 표준 라플라스 역변환 공식들을 활용하면 어렵지 않게 얻을 수 있다.

#### 초기값 문제와 경계값 문제

2차 미분방정식 해석에서, 경계값 문제(boundary value problem)와 초기값 문제(initial value problem)는 본질이 다르다. 라플라스 변환은 대개 초기값 문제를 다루는 데 특화되어 있다. $t=0$에서 주어지는 초기상태로부터 $t>0$ 범위의 시간응답을 구하는 문제가 주된 대상이다.

경계값 문제, 예컨대 $y(0)=0$과 $y(T)=1$ 같은 조건이 주어지는 경우에는 라플라스 변환이 조금 덜 직접적으로 쓰인다. 그럼에도, 구간 $\[0,T]$에서의 해를 구하고, $t=T$를 기준으로 적분 불연속 등을 고려할 수도 있긴 하지만, 통상 푸리에 급수나 그린 함수 전개 등을 더 많이 활용한다.

따라서 2차 상수계수 미분방정식에서 라플라스 해법은 크게 초기값 문제를 해석하는 수단으로 집중적으로 사용된다고 볼 수 있다.

#### 그린 함수 관점에서의 라플라스 해법

2차 선형 연산자

$$
L\[y] = y'' + a\_1 y' + a\_2 y
$$

에 대해 임펄스응답(그린 함수) $g(t)$를 정의하면, 모든 외부입력 $f(t)$에 대해

$$
y(t) = (g \* f)(t)
$$

와 같은 컨볼루션 구조가 성립한다. 라플라스 변환에서 $G(s)=\frac{1}{s^2 + a\_1 s + a\_2}$에 해당하므로, $f(t)$의 라플라스 변환이 $F(s)$일 때

$$
Y(s) = G(s),F(s).
$$

이를 역변환하면 컨볼루션 형태가 얻어진다. 이러한 관점에서, 라플라스 해법이란 결국 그린 함수를 간단한 분수 형태로 얻고, 이를 필요한 입력과 곱한 뒤 역변환하여 컨볼루션 적분을 계산하는 과정과 동일하다. 미분방정식을 직접 풀 때에도 특정해 구하기와 균일해 구하기를 조합하면 본질적으로 같은 결과에 도달한다.

#### 더 나아간 활용과 복잡계 연장

실제로는 비선형 항이 추가되거나, 계수들이 시간에 따라 달라지는 시변계(time-varying system)일 수 있다. 그 경우 라플라스 변환이 직접적으로 적용되지 않을 수 있으며, 미분방정식 풀이도 훨씬 까다로워진다. 하지만 수많은 물리·공학적 문제들이 선형성과 상수계수를 전제하거나, 적어도 국소적으로 선형근사를 적용할 수 있는 형태이므로, 2차 미분방정식에 대한 라플라스 해법은 여전히 중요한 위치를 차지한다.

또한, 전형적인 2차 계를 벗어나서 고차 계나 편미분방정식 등으로 확장해도, 부분적으로 라플라스 변환이 적용 가능한 장면이 많다. 예를 들어, 라플라스-푸리에 이중 변환으로 2차원 파동방정식이나 열방정식을 다루는 기법들이 대표적이다. 이렇게 라플라스 변환은 1차원 상미분방정식 수준을 넘어서는 폭넓은 응용으로도 이어진다.

\---과 심화 참고

지금까지 2차 상수계수 선형 미분방정식에 대해 라플라스 변환을 어떻게 적용하고, 어떤 특징적 결과가 나타나는지 여러 관점에서 논의해 왔다. 초기조건이 주어진 문제를 풀거나, 물리·공학적 시스템(기계계, 전기회로 등)을 해석할 때 라플라스 방법이 매우 강력한 도구임을 확인했다. 일반 해법 과정을 다시 간단히 배열해 보자.

(1) 미분방정식을 시간영역에서 설정한다:

$$
y''(t) + a\_1 y'(t) + a\_2 y(t) = f(t),
$$

초기조건 $y(0)=y\_0$, $y'(0)=v\_0$.

(2) 라플라스 변환을 취한다:

$$
s^2 Y(s) - s y\_0 - v\_0 + a\_1 \bigl(sY(s)-y\_0\bigr) + a\_2 Y(s) = F(s).
$$

(3) $Y(s)$에 대한 대수 방정식을 푼다:

$$
Y(s) = \frac{F(s) + \dots(\text{초기조건항})}{s^2 + a\_1 s + a\_2}.
$$

(4) 부분분수 분해 등으로 적절히 $Y(s)$를 단순화한다.

(5) 역라플라스 변환 $\mathcal{L}^{-1}{Y(s)}$을 통해 시간영역 해 $y(t)$를 구한다.

이 방식은 외부 강제항의 종류나 초기조건이 얼마든지 달라져도 동일하게 적용된다. 실제로 문제를 풀 때는 주어진 $f(t)$가 $\delta(t)$, $u(t-T)$, $\sin(\omega t)$ 등일 수 있으며, 때때로 여러 항이 조합된 형태를 띨 수도 있다. 하지만 선형성과 부분분수 분해, 그리고 라플라스 변환 표나 컨볼루션 정리를 활용하면, 최종 해를 구하는 과정은 늘 일정한 절차로 귀결된다.

2차 방정식 해법 자체가 워낙 널리 알려져 있는 만큼, 라플라스 변환 기법을 정교하게 익혀두면 고차나 더욱 복잡한 시스템에도 자연스럽게 확장 적용할 수 있게 된다. 또, 해석 과정에서 특성방정식 해석, 그린 함수, 임펄스응답, 전달함수 등 유용한 개념들과 서로 교차 참조하면서 공부하면, 실제 응용에서 보다 정확하고 빠르게 문제를 해결할 수 있다.