공업수학 시리즈 4편: 2차 선형 미분방정식

1차 미분방정식이 "현재 상태의 변화율"을 다뤘다면, 2차 미분방정식은 가속도나 곡률까지 포함된 시스템을 다룹니다. 그래서 진동, 회로, 기계 시스템처럼 관성이 있는 현상을 설명할 때 자주 등장합니다.

표준형

가장 기본적인 2차 선형 미분방정식은 다음과 같습니다.

$ay'' + by' + cy = g(x)$

여기서 $a$, $b$, $c$가 상수이고 $g(x)$가 외부 입력입니다. 먼저 외부 입력이 없는 동차식부터 보면

$ay'' + by' + cy = 0$

이 됩니다.

왜 지수함수를 먼저 떠올릴까

상수계수 선형 방정식에서는 지수함수 $y = e^{rx}$를 넣어 보면 미분할 때 모양이 유지됩니다. 이 때문에

$ar^2 + br + c = 0$

이라는 특성방정식이 나오고, 해의 구조는 이 2차 방정식의 근에 따라 결정됩니다.

세 가지 기본 경우

서로 다른 두 실근

근이 $r_1$, $r_2$이면

$y = C_1 e^{r_1x} + C_2 e^{r_2x}$

입니다.

중근

근이 $r$ 하나로 겹치면

$y = (C_1 + C_2 x)e^{rx}$

입니다.

켤레복소근

근이 $\alpha \pm i\beta$이면

$y = e^{\alpha x}\left(C_1 \cos \beta x + C_2 \sin \beta x\right)$

가 됩니다.

이 마지막 형태가 바로 감쇠 진동과 연결됩니다.

손으로 풀어보는 예제

다음 문제를 보겠습니다.

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

특성방정식은

$r^2 - 3r + 2 = 0$

이고 인수분해하면

$(r-1)(r-2) = 0$

이므로 근은 $r=1$, $r=2$입니다. 따라서 일반해는

$y = C_1 e^x + C_2 e^{2x}$

입니다.

초기조건이

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

이라면

$C_1 + C_2 = 1$

$C_1 + 2C_2 = 0$

을 얻고, 풀면

$C_1 = 2, \quad C_2 = -1$

이므로

$y = 2e^x - e^{2x}$

가 됩니다.

공학 응용

질량-스프링-감쇠기 시스템의 자연응답은 보통

$m x'' + c x' + kx = 0$

으로 씁니다. 여기서

• $m$은 질량
• $c$는 감쇠 계수
• $k$는 스프링 상수

입니다.

회로에서는 같은 구조가

$Lq'' + Rq' + \frac{1}{C}q = 0$

처럼 RLC 회로에 나타납니다. 즉 기계 시스템과 전기 시스템이 같은 수학 구조를 공유합니다.

해를 해석하는 법

공업수학에서는 해를 구한 뒤 반드시 거동을 읽어야 합니다.

• 지수함수 양의 실근이면 발산
• 음의 실근이면 감쇠
• 복소근이면 진동
• 실수부가 음수인 복소근이면 감쇠 진동

즉 수식은 결국 시스템의 안정성과 응답 속도를 알려주는 도구입니다.

자주 하는 실수

특성방정식까지만 풀고 해를 잘못 쓴다

중근과 복소근의 해 형태는 서로 다릅니다. 근만 구하고 일반해 형태를 기계적으로 쓰면 틀리기 쉽습니다.

초기조건 적용에서 미분을 대충 한다

$y'$를 정확히 구하지 않고 상수만 맞추려 하면 오차가 납니다. 일반해를 먼저 정리한 뒤 차분하게 미분하는 편이 안전합니다.

해의 의미를 읽지 않는다

공학 문제에서는 발산, 감쇠, 진동 여부가 핵심입니다. 계산만 끝내고 시스템 해석을 빼먹지 않는 습관이 중요합니다.

한 줄 요약

2차 선형 미분방정식은 특성방정식의 근을 통해 자연응답의 형태를 읽는 도구입니다.

다음 편 예고

다음 글에서는 2차 방정식을 실제 물리 시스템에 연결해 감쇠, 강제진동, 공진이 어떻게 나타나는지 살펴보겠습니다.

참고자료

• Erwin Kreyszig, Advanced Engineering Mathematics, 10th Edition
• Steven H. Strogatz, Nonlinear Dynamics and Chaos
• MIT OpenCourseWare, Vibrations and Waves 관련 자료